3.1. Overview of Oxalic-Acid-Vapor* and Particulate Oxalate in the Coastal Atmosphere
Figure 2a,b showed the hourly averaged molecular concentrations of oxalic-acid-vapor* and mass concentrations of particulate oxalate in PM
2.5 during Period 1 and Period 2. The concentrations of oxalic-acid-vapor* were (4.4 ± 1.0) × 10
8 molecules cm
-3 in February and (3.5 ± 0.7) × 10
8 molecules cm
-3 in April during Period 1. Although the difference between two average values was small, a significant difference in oxalic-acid-vapor* concentrations between February and April was observed through a t-test (P < 0.05). The concentrations decreased to (1.5 ± 0.5) × 10
8 molecules cm
-3 during Period 2, with a significant difference compared to those in February and April during Period 1 (P < 0.01). In contrast, the coefficient of variation (CV) for oxalic-acid-vapor* concentrations slightly increased during Period 2 relative to Period 1. All these molecule concentrations of oxalic-acid-vapor* were 1-2 orders of magnitude higher than those of sulfuric acid vapor reported in ambient air [
30,
31], suggesting that either oxalic acid vapor is substantially more volatile than sulfuric acid regardless of atmospheric NH
3, or the true concentration of oxalic acid vapor is lower than the detected value. For particulate oxalate in PM
2.5, the mass concentrations were 0.05 ± 0.03 μg m
-3 in February and 0.07 ± 0.03 μg m
-3 in April during Period 1, with a significant difference between the two (P < 0.05). During Period 2, the mass concentrations of oxalate in PM
2.5 were 0.03 ± 0.02 μg m
-3, showing a significant difference from Period 1 (P < 0.01). When the values measured simultaneously in both phases were compared, the average mass concentration of oxalic-acid-vapor* was almost equal to that of oxalate in PM
2.5 during Period 1 (0.059 μg m
-3 vs. 0.057 μg m
-3). During Period 2, the two averages were also comparable, with 0.021 μg m
-3 for oxalic-acid-vapor* and 0.029 μg m
-3 for oxalate in PM
2.5.
Section 3.3 will examine whether the ratios of oxalic-acid-vapor* to oxalate in PM
2.5 are theoretically supported by thermodynamic equilibrium considerations.
The observed average concentration of oxalic-acid-vapor* during Period 1 was slightly lower than the 0.074 μg m
-3 measured at a semi-urban site at Sydney, FL, from 27 April to 31 May 2002 [
21], but nearly double the value observed in Yorkville, Geogia, from mid-August to mid-October 2006 [
24,
25]. In contrast, the observed average concentration of oxalate in PM
2.5 during Period 1 was only approximately 20% of the average 0.29 μg m
-3 measured in Sydney. The average concentration of oxalic-acid-vapor* during Period 2 was close to the value of 22 ± 8 ng m
−3 measured at a rural site in Shanghai, China during the cold season from 8 November to 2 December 2020, using an identical AIM-IC system [
22]. However, Yao et al. [
22]. reported oxalate concentrations in PM
2.5 as high as 475 ± 266 ng m
−3. Yao et al.. measured oxalic-acid-vapor* at an urban site in Hong Kong using Na
2CO
3-coated dry denuders in December 2000 [
20]. They reported only four values in the range of 0.02-0.05 μg m
-3, with the rest of the measurements being below twice the field blank. The Na
2CO
3-coated dry denuder likely absorbed fewer gaseous precursors of oxalic acid and oxidants in comparison with the wet denuder used in the studies above, although no direct comparison test has been conducted.
Periodic analysis of the two observed variables in
Figure 2a,b revealed that the mass concentrations of oxalate in PM
2.5 exhibited short-period variations (3-7 days), superimposed by higher-frequency diurnal variations in different extents during both Period 1 and Period 2. The (3-7)-day variations likely correspondeded to planetary waves in the northern hemisphere, reflecting long-range transport of oxalate in PM
2.5. For example, the strongest diurnal variations occurred during 18-22 April 2023, coinciding a 5-day event (
Figure 2a,b and
Figure S3a). In contrast, much weaker diurnal variations were observed during 11-13 February and 18-23 February 2023, with other periods falling between these extremes. The concentrations of oxalic-acid-vapor*, however, exhibited distinct diurnal fluctuations (
Figure 2a,b). During Period 1, the concentrations of oxalic-acid-vapor* peaked in the afternoon, in line with the O
3 concentrations shown in
Figure S3c. In Period 2, however, the highest concentrations of oxalic-acid-vapor* occurred later in the evening, except an increase during the morning rush hours, corresponding to a trough in O
3 levels and a slight rise in PM
2.5 (
Figure S3b,d). The lower mixing boundary layer in the morning favored the accumulation of air pollutants, enhancing the rush-hour effect. However, O
3 concentrations still reached their peak later in the afternoon during Period 2 (
Figure S3d). Further discussion of the diurnal fluctuations of oxalic-acid-vapor* and oxalate in PM
2.5 will be provided in
Section 3.2.
When analyzing the correlation between oxalic-acid-vapor* and particulate oxalate in PM
2.5 mass concentration, a significant negative correlation was obtained during Period 1 (
Figure 2c, black markers, and the black regression curve). This implied a source-sink relationship between them, which is further supported by the time delay in the peak concentration of oxalate in PM
2.5 relative to that of oxalic-acid-vapor* in the diurnal variation curve (
Figure S3a). In a closed system, the uptake of oxalic-acid-vapor* as the primary source of oxalate in PM
2.5 would theoretically lead to a perfect negative correlation. However, the correlation observed was relatively weak, which could be attributed to several factors, e.g., 1) the presence of oxidizable precursors in oxalic-acid-vapor*; 2) fluctuating external inputs of particulate oxalate over the 3-7 day periods; and 3) secondary formation of oxalate in PM
2.5 from its precursors, etc. Notably, a significant positive correlation between them was observed during certain short periods, such as 20-23 April 2023 when secondary formation of oxalate in PM
2.5 may have occured simultaneously with that of oxalic-acid-vapor* (
Figure 2c, red markers).
During Period 2, a significant positive correlation was obtained between oxalic-acid-vapor* and particulate oxalate in PM
2.5, with a low R
2 value (
Figure 2d). This low R
2 values occurred throughout Period 2, even during the period with the strongest diurnal variations in oxalate in PM
2.5 (
Figure 2d, red markers). In these instances, freshly formed oxalic-acid-vapor* likely overwhelmed its condensation sink on atmospheric particles, eliminating their inherent negative correlation. Again, this implied that oxalic acid has a high vapor pressure and does not exist in a quasi steady state. Theoretically, an increase in the the source term relative to the sink term would strengthen the positive correlation. This positive correlation could also be explained by: 1) the formation rates of oxalic-acid-vapor* being the determining step for both species; 2) the inherent limitations of Euler observation at the fixed site,where oxalic-acid-vapor* and oxalate in PM
2.5 might be affected by confounding factors, regardless of their origins or formation pathways. Unfortunately, obtaining observations of moving atmospheric flows for either variable is practically impossible.
3.2. Daytime and Nighttime Formation of Oxalic-Acid-Vapor* and Particulate Oxalate
Based on the distinguishable diurnal variations of oxalic-acid-vapor* concentrations during Period 1 and its positive correlation with O
3 (
Figure S3a,c), it can be reasonably inferred that oxalic-acid-vapor* is primarily produced through daytime photochemical reactions [
11]. The lack of a positive response of oxalic-acid-vapor* to most (3-7)-day variations in oxalate in PM
2.5 suggested that long-range transport of oxalic-acid-vapor* and its precursors made only a negligible contribution to the observed oxalic-acid-vapor*. For example, the lower concentrations of oxalate in PM
2.5 were observed with higher concentrations of oxalic-acid-vapor* during 10-13 February than other moments during Period 1. In addition, the increase in oxalic-acid-vapor* during morning rush hours implied a signifiant primary source for it
11 (
Figure S3a). In fact, the primary source of oxalic-acid-vapor* in morning rush hours was also identifiable during Period 2 (
Figure S3b). However, the higher concentrations of oxalic-acid-vapor* occurred at nighttime during Period 2 (
Figure 2d and
Figure S3b), which contrasted with the even larger maximum concentrations of O
3 occurring in the afternoon than Period 1 (
Figure S3c,d). The strong northeast monsoon winds blowing from the sea during the daytime of Period 2 swept local precursors (
Figure S4–S6), while at night, the mixing layer height decreased significantly, allowing the accumulation of local precursors and facilitating their chemical conversion to oxalic-acid-vapor* through dark chemistry [
11]. The role of dark chemistry in this process warrants further investigation.
During Period 1, the (3-7)-day variations in oxalate concentrations in PM
2.5 were likely driven by the long-range transport, as the air masses could have traveled over thousands of kilometers (
Figure S4–S6). The superimposed diurnal variations in oxalate in PM
2.5 were more likely attributed to the condensation of oxalic-acid-vapor* from local precursors through photochemical reactions, which is supported by the negative correlation and time delay between oxalic-acid-vapor* and oxalate in PM
2.5, as mentioned eralier. When the diurnal variations of oxalate in PM
2.5 were negligible, the concentrations of oxalate in PM
2.5 were found to be highly correlated with SO
42- concentrations in PM
2.5, e.g., during 18-22 February 2023 (
Figure 3a, blue markers). Following the logic from previous studies [
11,
35,
43,
44], this strong correlation highly indicated in-cloud processing of particulate oxalate, although it is unclear whether this process predominantly occurred during the daytime or nighttime. However, during periods of strong photochemical formation of oxalic acid vapor, followed by condensation on PM
2.5 (e.g., 20-23 April), this significant but weak correlation between oxalate and SO
42- in PM
2.5 existed for the whole Period 1 was even disrupted (
Figure 3a).
During Period 2, a moderate correlation between oxalate and SO
42- in PM
2.5 also suggested that in-cloud processing of particulate oxalate was the dominant formation pathway. The percentage increase in oxalate concentrations at nighttime was notably larger than that of the total mass of PM
2.5 (
Figure S3b,d), which can be attributed to dark chemistry conversion, compounded by the low mixing layer height that facilitated the accumulation (
Figure S3b). In addition, the strong northeast monsoon normally dominated over the weakened sea-land breeze in this season (
Figures S4–S6a). As with Period 1, the (3-7)-day variations in oxalate concentrations in PM
2.5 during Period 2 were likely influenced by either daytime or nighttime in-cloud processing of particulate oxalate during long-range transport. Moreover, when we correlated the concentrations of oxalate in PM
2.5 with ambient RH (
Figure S7), no positive correlation was observed as Yang et al. [
37] suggested, e.g., lower concentrations of particulate oxalate corresponding to higher RH (11-14 February), in contrasting with higher concentrations corresponding to lower RH (21 February). This lack of correlation implied that aerosol droplet chemistry was not a major factor in the formation of particulate oxalate, at least in this study. Further analysis correlating the concentrations of oxalate in PM
2.5 + oxalic acid vapor* with those of SO
42- in PM
2.5 (
Figure S8) showed that the correlations were expectedly weaker during Period 1, as oxalic-acid-vapor* was likely generated via photochemical reactions in the gas phase. However, during Period 2, the determination coefficient remained nearly the same, with a doubled slope, indicating that oxalic-acid-vapor* might have been evaporated from the particulate phase and generated through a similar process as particulate oxalate during Period 2.
3.3. Thermodynamic Examination for True Oxalic Acid Vapor and Its Volatility
The E-AIM model was further employed to examine the gas-aerosol partitioning of oxalic acid species using the full ambient data during the study period (
Table S1). For approximately 60% of the observational data (883 hours), the model predicted that over 95% of (oxalic-acid-vapor + particulate oxalate) would be partitioned into the gas phase. The reason will be analyzed later. For the remaining ~ 40% of the observational data (495 hours), the model predicated (oxalic-acid-vapor + particulate oxalate) were more evenly distributed between the two phases. The parts of the equilibrated concentrations of oxalate in PM
2.5 predicted by the model in different periods were plotted against the observed concentrations, as shown in
Figure 4a-d, and the differences between predicted oxalic-acid-vapor concentrations and the observed values were also plotted against the observed concentrations, as shown in
Figure 4e-h. When considering a 10% analytical error for all chemical inputs, only a small fraction (~9%) of the equilibrated concentrations of oxalate in PM
2.5 were consistent with the observed values within a 20% margin along the 1:1 line during different periods (
Figure 4a-d). The consistency indicated that 1) gas-aerosol equilibrium has been achieved; 2) for certain cases in the coastal atmosphere, the vapor pressure of oxalic acid reported by Booth et al. [
27] is valid. Considering that the estimated timescale for partitioning inorganic and organic volatility compounds into liquid particles in hot ambient air ranges from seconds to minutes [
50,
51] this short equilibration timescale implies that gas-aerosol equilibrium should be reached. However, a significant discrepancy was observed in 85% of cases, where the equilibrated concentrations of oxalate in PM
2.5 were over 20% larger than the observed values. In these cases, the equilibrated concentrations of oxalic acid vapor were lower than the observations. The deviation indicated that the observed oxalic-acid-vapor* might include a fraction of precursors to oxalic acid, rather than only the oxalic acid itself. The maximum deviation was observed in February 2023, when episodically high concentrations of oxalic-acid-vapor* were detected. By assuming that the overestimation of equilibrated oxalate concentrations in PM
2.5 (more than 20%) was entirely due to the presence of oxalic-acid-vapor precursors, the contributions of the precursors to oxalic acid vapor was estimated to be approximately 69%, 57%, 73% and 70% for February, April, and before and after 15:00 on 13 December during Period 2, respectively. These results were partially supported by additional observations after 15:00 on 13 December 2023, where the observed oxalic-acid-vapor* concentrations, measured with pure H
2O in the wet denuder, were (1.0 ± 0.3) × 10
8 molecules cm
-3 and significantly lower than (1.6 ± 0.4) × 10
8 molecules cm
-3 measured during other times in Period 2 (P < 0.01). Assuming that the difference in the average values was completely due to the presence of non-oxalic-acid-vapor (its precursors), it was estimated that these precursors accounted for approximately 1/3 of the observed oxalic-acid-vapor* before 15:00 on 13 December during Period 2. The estimated precursors of oxalic acid vapor might be considered the lower values. The pure-H
2O-added wet denuder can still absorb oxidants and precursors from the sampled ambient air, leading to the aqueous conversion of these precursors into oxalic acid. In the future, a reductant will be tested during specific periods to further investigate potential artifact.
The equilibrated concentrations of oxalate in PM
2.5 were more than 20% lower than observed in 6% of cases, where the modeled liquid water content (LWC) was generally less than 1 μg m
-3. This underestimation strongly indicated the vapor pressure of ammnium oxalate, or a combination of ammnium oxalate and oxalic acid vapor pressures, should be used in place of oxalic acid alone for predicting the solid phase. A similar conclusion may apply to approximately 60% of cases, where the model predicted that over 95% of the (oxalic-acid-vapor + particulate oxalate) were partitioned in the gas phase. We analyzed the inputs and outputs for ~60% versus ~40% observational cases in
Table S1, with the major difference being lower RH in the ~60% cases. However, the potential limitations of the E-AIM model in simulating dry aerosols cannot be ruled out.
When the modeled ratios of oxalic acid to particulate oxalate were plotted against aerosol pH in
Figure 5a-d, the data with modeled LWC < 1 μg m
-3 showed considerable scattered. Excluding these data, the ratios increased as pH decreased, ranging from < 0.01 to 1.0, with a median value of 0.3 in February 2023. The corresponding values for April were <0.01, 4.1 and 0.4; for November and December, they were <0.01, 2.1 and 0.2; and for December with pure H
2O used in wet denuder, the values were 0.01. 1.0 and 0.2. When comparing the modeled ratios to the corresponding observational values, the modeled ratios were consistently larger with the modeled LWC < 1 μg m
-3. Again, the vapor pressure of ammonium oxalate, rather than oxalic acid, should be used in such cases [
27,
28]. For the cases with modeled LWC ≧ 1 μg m
-3, smaller ratios were generally obtained, possibly due to artifacts related to the collected oxalic acid vapor precursors. Interestingly, the observed ratios in April 2023 were more narrowly distributed in the low-value range of 0.4-1.4, with a median value of 0.8, compared to other months, regardless of black markers or color markers in
Figure 5f. It is noteworthy that the photochemical formation of oxalic acid vapor was significantly stronger in April than that in the November and December. This suggested a stronger source of particulate oxalate during April, which will be further explored in the next section.
3.4. Size Distributions of Oxalate in Atmospheric Particles over the SCS
To investigate the notably lower ratios of oxalic acid vapor to particulate oxalate observed in April than other months, we further analyzed the formation of particulate oxalate in the upwind marine atmosphere using historic data from April of a different year.
Figure 6 showed the mass size contributions of water-extracted oxalate in atmospheric particles over the SCS in April 2017. Although the mass concentrations of oxalate generally exhibited a bimodal distribution across particle sizes, the contribution of supermicron oxalate to the total particulate mass concentration varied significantly between samples. Specifically, supermicron oxalate exceeded submicron oxalate on the dates 2-6, 7-10, 15-17, 18 and 23-25 April 2017, with this distribution primarily occuring south of 15°N latitude (
Figure 6a,h). Conversely, the reverse was observed on other dates (
Figure 6b,h). To investigate the origins and formation pathways of both submicron and supermicron oxalates, whether anthropogenic or natural, we performed a comparative analysis involving particulate non-sea-salt-SO
42- (nss-SO
42-), dimethylaminium (DMA
+), Na
+, NO
3- and nss-K
+ (
Figure 6c-g).
The possible primary anthropogenic sources of particulate oxalate might include emissions from vessel combustion. To examine the impact of the vessel combustion emissions, samples were collected while the vessel was moored with its engines running on 15-17 and 18 April. Extremely high concentrations of nss-SO
42-, apparently resulting from vessel combustion, were observed, predominantly in the size range below 0.32 µm (marked with blue inverted triangles and diamonds in
Figure 6c). The contribution of nss-SO
42- in the size ranges above 0.32 µm in the two moored samples were negligible, as these concentrations were lower than those observed in some size bins from other samples unaffected by vessel combustion emissions (
Figure 6c). In contrast, an opposite pattern was observed for oxalate, with negligible amounts in the size range below 0.32 µm in the two moored samples, indicating that the contribution of the vessel combustion emissions to oxalate levels was minimal relative to the oxalate measured in this study. Similar to the high concentrations nss-SO
42-, DMA
+ was also predominantly distributed in the size range below 0.32 µm (
Figure 6d). DMA
+ likely neutralized particulate inorganic and organic acids through gas-particle condensation [
7,
9]. In addition, heterogeneous conversions of glycolaldehyde and glyoxal, etc., on sea-salt or crustal aerosols, potentially enhanced by the photochemical degradation of NO
3- to generate OH free radical in deliquesced sea-salt aerosol droplets under extremely strong UV radiation [
52] (
Figure S9b), could theoretically produce particulate oxalate in a supermicron mode, similar to NO
3-, with a mode ranging between 2-6 µm [
11]. However, this was not observed in this study. The observed particulate oxalate over the SCS was dominantly distributed in the 0.5-2.0 µm range, suggesting it was primarily formed through in-cloud processing of oxalate, i.e., 0.5-1.0 µm via mid-high cloud and 1.0-2.0 µm via low cloud (or fog) [
32,
34,
35,
43,
54].
Rinaldi et al. [
34] reported a similar finding of a supermicron mode of oxalate at 1-3 µm in coastal atmospheres, attributing it to the oxidation of natural marine precursors. A similar bi-modal pattern of particulate oxalate was also reported in marine atmospheres over the northwest Pacific Ocean (NWPO) [
43,
55]. As shown in
Figure S9a, most oceanic zones over the SCS are part of the ocean desert zone, characterized by low chl-a concentrations and limited marine biogenic emissions. This restriction in oceanic biogenic VOCs emissions to specific hotspots during that month (
Figures S4–S6), making locally emitted VOCs the primary precursors of oxalic-acid-vapor* and causing its strong diurnal variation in relation to O
3 levels. Moreover, the absence of a supermicron mode of nss-SO
42- with oxalate at 1-3 µm might explain the lack of a significant correlation between oxalate and SO
42- in PM
2.5 observed at the coastal site in April 2023. Whether the supermicron oxalate mode is related to fog processing of sea-spray organic aerosols requires further investigation. However, the strong supermicron mode of oxalate likely contributed to the low ratios of oxalic acid vapor to particulate oxalate observed in that month. In terms of natural sources, biomass burning has been considered a significant contributor to oxalate in marine atmospheric particles [
55,
56]. However, in this study, biomass burning appeared to be a negligible source of particulate oxalate over the SCS, as indicated by the low concentrations of nss-K
+. The analysis over the SCS strongly supported the hypothesis that the observed (3-7)-day variations in particulate oxalate at the coastal island site were related to in-cloud formation of oxalate in the upwind marine atmosphere, despite the fact that the former and latter observations were conducted at different times.