Large contribution of fossil

Apr 29, 2023

Nature Communications volume 13, Article number: 5115 (2022) Cite this article

4655 Accesses

8 Citations

11 Altmetric

Metrics details

Incomplete understanding of the sources of secondary organic aerosol (SOA) leads to large uncertainty in both air quality management and in climate change assessment. Chemical reactions occurring in the atmospheric aqueous phase represent an important source of SOA mass, yet, the effects of anthropogenic emissions on the aqueous SOA (aqSOA) are not well constrained. Here we use compound-specific dual-carbon isotopic fingerprints (δ13C and Δ14C) of dominant aqSOA molecules, such as oxalic acid, to track the precursor sources and formation mechanisms of aqSOA. Substantial stable carbon isotope fractionation of aqSOA molecules provides robust evidence for extensive aqueous-phase processing. Contrary to the paradigm that these aqSOA compounds are largely biogenic, radiocarbon-based source apportionments show that fossil precursors produced over one-half of the aqSOA molecules. Large fractions of fossil-derived aqSOA contribute substantially to the total water-soluble organic aerosol load and hence impact projections of both air quality and anthropogenic radiative forcing. Our findings reveal the importance of fossil emissions for aqSOA with effects on climate and air quality.

Organic aerosol (OA) is the dominant component of atmospheric fine particulate mass1. The majority of OA is secondary (SOA), formed in the atmosphere from the oxidation of precursor gases2. Most SOA are composed of oxygenated and hygroscopic compounds, thereby, have harmful impacts on respiratory health and large, and uncertain effects on atmospheric radiative forcing3. However, current global models typically underpredict SOA magnitude, distribution, and dynamics, suggesting a limited understanding of their sources and formation processes4. An increasing number of modeling and experimental studies point toward aqueous chemical reactions occurring in cloud droplets and wet aerosols as an important missing pathway for SOA formation2,3,4,5,6,7,8. Aqueous SOA (aqSOA) formation often involves mixed effects of anthropogenic and biogenic emissions2, 3. The extent to which aqSOA is controllable or "natural" is an important subject of ongoing research for both air quality and climate issues2, 3. However, the poorly understood aqueous chemical processes and the complexity of anthropogenic-biogenic interactions renders any firm conclusion hard to come by2.

In model simulations, aqSOA are put to largely derive from the natural biogenic precursors, such as the oxidation products of isoprene3, 4, 8. Fossil precursors are noticed to be less polar and less hydrophilic than biogenic precursors, thus diminishing the possibility of partitioning to the aqueous phase, and subsequent aqSOA formation9. However, a few recent studies in East Asia suggest that fossil precursors may substantially contribute to the formation of aqSOA. For example, fast aqueous-phase oxidation of polycyclic aromatic hydrocarbons emitted from fossil fuel combustion was observed in Beijing winter haze and may explain the observed SOA10. Radiocarbon-based estimations attributed ~50% of water-soluble OA to fossil sources in the East Asia outflow11, 12, compared to <30% in Europe, U.S., and South Asia13,14,15,16,17,18. Nevertheless, on aqSOA formation, a quantitative estimation of the fossil contribution remains challenging, largely due to the lack of specific and reliable technical tools in tracking the sources and formation mechanism of aqSOA in ambient aerosols.

Here we leverage the recent advent of compound-specific dual-carbon isotope fingerprint (δ13C−Δ14C) of aqSOA molecules19 to quantify and characterize aqSOA sources and atmospheric chemical processes. The molecular-level Δ14C contents of aqSOA tracers will provide direct constrain on their origins19, whereas δ13C fingerprints of molecular tracers can differentiate various atmospheric processes/reactions20. The technique was applied to purified oxalic acid and other similarly abundant organic acids (e.g., glyoxylic acid). Oxalic acid, likely represents the highest oxidation state of OA, has frequently been used as a proxy for aqueous processing in both laboratory and field studies21,22,23. Organic acids that have high O/C ratios (around 1−2), are among the most abundant SOA components and are key end products in the aqueous-phase photochemical oxidation of various volatile organic compounds (VOCs) and intermediates such as glyoxal (Gly) and methylglyoxal (MeGly) in clouds or wet aerosols4, 5, 24, 25. Therefore, these organic acids can serve as signature compounds to trace the aqSOA formation pathways21,22,23, 26.

We used a set of year-round aerosol samples from a regional receptor site (Heshan Atmospheric Environmental Monitoring Superstation, Supplementary Fig. 1) in South China. We measured the aqSOA compounds and their δ13C-Δ14C isotope fingerprints during two contrasting phases of the East Asia monsoon system, i.e., continental outflow vis-à-vis South China Sea (SCS) air masses (Fig. 1). By combining the dual-carbon isotopic compositions and other chemical indicators together with meteorological parameters, we elucidate the relative contribution of fossil and biomass aqSOA precursors and their subsequent evolution processes, and provide ambient isotopic observational evidence for a large aqSOA formation from anthropogenic fossil emissions. Finally, the dual δ13C-Δ14C characterization of aqSOA molecules are extended to wider spatial coverage of China, which further confirm a ubiquitous contribution of aqueous-phase transformation of fossil precursors to organic aerosols.

The average aerosol optical depth (AOD) at 550 nm during June 2017 to May 2018 over the East Asia region are shown. The mean air mass back trajectories (BT) based on HYSPLIT cluster analysis show two dominant air mass transport pathways (dashed line and arrows): continental outflow (62% of the total clusters) and oceanic air masses (38% of total clusters). See Supplementary Fig. 2 for the 72-h BT with a 6-h interval and Supplementary Fig. 3 for relative cluster contribution. The bar charts depict the mass concentrations of total carbonaceous aerosol [TCA = organic matter (OM) + elemental carbon (EC)], anthropogenic WSIC [WSICanth = non-sea-salt (nss) K+ + nss-SO42- + NH4+ + NO3− + nss-Cl−], and aerosol liquid water (ALW). The AOD data is obtained from Moderate Resolution Imaging Spectroradiometer (MODIS) observations (https://giovanni.gsfc.nasa.gov/giovanni/). The coastline boundaries in the map are originated from Natural Earth free vector map data (https://www.naturalearthdata.com/). The administration boundaries in the map are originated from map products of National Geomatics Center of China (https://www.webmap.cn/).

The Heshan receptor site features an East Asia Monsoon Climate. The air is controlled by continental monsoon in winter and oceanic monsoon in summer, respectively (Fig. 1). Backward trajectories revealed that clean SCS air mass dominated the airshed during June-August 2017, while continental outflow rendered by enhanced anthropogenic emissions from the Pearl River Delta (PRD) city clusters, dominated during September 2017 to March 2018 (Fig. 1). Therefore, we categorized the 2017-2018 samples into two groups based on their distinctly different air mass transport regimes (Supplementary Table 1), i.e., "coastal background (n = 11)" and "continental outflow (n = 21)" (see Supplementary Fig. 2 for air mass clusters and Supplementary Fig. 3 for relative cluster contribution).

The chemical composition of aerosols during the campaign varied greatly with the air mass transport regimes. Concentrations of gas-phase pollutants (SO2, NO2, O3, CO) and major components of PM2.5 during the continental outflow were 2-5 times higher than during the coastal background (Supplementary Table 2 and Supplementary Fig. 4). The observed organic matter (OM) and anthropogenic water-soluble inorganic constituents (WSICanth) (Supplementary Method 1) increased in continental outflow samples, along with the aerosol liquid water (ALW) content associated with hygroscopic aerosol components (Supplementary Method 2). The concentrations of OM, WSICanth, and ALW show a clear increase in continental outflow samples with 3, 5, and 3 times higher values than the coastal background, respectively (Fig. 1).

The molecular composition of diacids and related aqueous processing tracers are compiled in Supplementary Table 3. The malonic to succinic acid ratios (C3/C4) are significantly higher in the coastal background samples (1.5 ± 0.4) than the continental outflow (0.9 ± 0.3) (Supplementary Table 4). Since malonic acid is produced by photochemical oxidation of succinic acid in the atmosphere, this indicates that the coastal background aerosols have undergone more extensive photochemical aging27, 28. The ratio of water-soluble organic carbon (WSOC) to organic carbon (WSOC/OC) is another indicator of atmospheric photochemical aging29. For example, the WSOC/OC ratios increased from ~40% in South Asian source region to ~70% in the Indian Ocean receptor site due to the photochemical aging of OC during long-range atmospheric transport17, 29. However, in this study, the WSOC/OC ratios of the continental outflow samples (44 ± 12%) cannot be distinguished statistically from the coastal background (39 ± 10%) (Supplementary Table 4), implying that additional mechanisms other than atmospheric aging may exist.

The uptake of water-soluble gaseous organic precursors to wet aerosol particles and fog/clouds, with subsequent aqueous-phase transformation to lower volatile compounds could be a significant pathway for WSOC aerosol formation23, 26, 30, 31. ALW mass was driven primarily by inorganic salts, especially by nitrate and sulfate31. As shown in Supplementary Fig. 5, the ALW mass well correlated with the concentrations of WSICanth and fractions of WSICanth in PM2.5 during the campaign. The nitrate in the continental outflow samples display concentrations 20 times higher than the coastal background samples, and account for nearly half of the anthropogenic WSIC (Supplementary Fig. 4). Because nitrate is more hygroscopic than sulfate30, such changes in aerosol compositions would be a dominant factor driving the ALW content. Given the high relative humidity (RH%) (74 ± 9%; Supplementary Table 2) and the abundant hygroscopic particles (particularly for particulate nitrate) in continental outflow samples, we expect that the changes in the gas-particle partitioning driven by the increased ALW could have enhanced the formation of aqSOA species, such as oxalic acid23, 26, 32. This can be supported by the sixfold higher aqSOA precursors (Gly + MeGly) mass fractions in OC observed in the continental outflow samples than the coastal background (Supplementary Table 4). Observation in continental outflow of higher aqSOA precursor concentrations and higher WSOC/OC ratios is consistent with aqSOA formation processes. Yet more precise and robust evidence/information is needed to give insights into the aqSOA formation.

The δ13C signatures reflect both sources and atmospheric processing. Emissions from different sources often have different end-member δ13C values. Once formed, the δ13C signature is then influenced by kinetic isotope effects (KIE) during atmospheric reactions/processing15. Secondary aerosols stemming from the oxidation of precursor gases due to the KIE typically lead to lower δ13C values in particulate products compared to the gaseous reactants, whereas aerosol aging (i.e., oxidation of organic matter in the aerosol phase) can lead to a release of isotopically lighter carbon such as in CO2 and CO, leaving residual (aged) aerosols enriched in 13C12, 15, 17, 29. Figure 2 shows the temporal variation of the mass concentrations and dual-carbon isotopic signatures of oxalic acid during the year-round observation. The δ13C values for oxalic acid are on average −24.6 ± 2.7‰ in the continental outflow samples and substantially more enriched in 13C in the coastal background samples (−19.9 ± 2.3‰). This clear contrast can likely be attributable to different emission sources and/or atmospheric processing.

The unshaded periods represent coastal background air mass regime, while the shadowed periods represent continental outflow air mass regime. Red and blue bars represent the concentration of oxalic acid during coastal background air mass regime and continental outflow air mass regime, respectively. The orange and green symbols show the δ13C composition and 14C-based fraction of non-fossil for oxalic acid.

Source mixing is not expected to have a major influence on the δ13C composition of oxalic acid in this study. First, the δ13C compositions of azelaic acid (C9 diacids), an oxidation product of biogenic unsaturated fatty acids33, in the coastal background samples (−34.1 ± 2.8‰) and continental outflow (−34.7 ± 2.6‰) are very close (Fig. 3a). These δ13C values point to unsaturated fatty acids of C3 plant origin (−38.5‰ to −32.5‰)34, while excluding marine phytoplankton and C4 plant sources that both instead are more enriched in δ13C35, 36. This is consistent with the limited marine-biogenic contributions as indicated by the mass ratios of nss-SO42- to SO42- (99 ± 1%) and nss-K+ to K+ (94 ± 4%) in all samples29. Given contributions from marine and C4 plants sources are mostly excluded at the Heshan site, the potential precursors of oxalic acid would then be anthropogenic VOCs (AVOCs; e.g., aromatic hydrocarbons) and C3 plant-derived biogenic VOCs (BVOCs; e.g., isoprene), of which the δ13C values are similar (Supplementary Table 5). Therefore, the isotopic range of various sources of oxalic acid at the supersite is minimal, leaving kinetic isotopic fractionation the most probable processes controlling the δ13C composition of oxalic acid.

a Average δ13C values of saturated aliphatic dicarboxylic acids (C2−C9), aromatic dicarboxylic acids [phthalic acid (Ph) and terephthalic acid (tPh)], and oxocarboxylic acids [glyoxylic acid (ωC2) and pyruvic acid (Pyr)] for the two air mass regimes. The error bars indicate the 95% confidence interval of δ13C values in each air mass regime. Pink shadows represent the compounds for which the δ13C values are statistically significantly different (p < 0.05) between the two air mass origins. b Relationships between oxalic acid mass concentration and aerosol liquid water (ALW). c Relationships between δ13C-oxalic acid and mass concentration of water-soluble organic carbon (WSOC) as a proxy for aqSOA abundance.

Oxalic acid is an SOA end product and can be formed by two pathways. The first pathway is a stepwise photochemical breakdown of longer-chain aliphatic diacids27, 28, which leads to relatively enriched δ13C values in the residual shorter-chain diacids37. Recently, significant amounts of oxalic acid followed by succinic acid were also found to be produced via ozone oxidation of isoprene under dry condition38. The second pathway is the oxidation of BVOCs and AVOCs to semi-volatile organic compounds (SVOCs, e.g., Gly and MeGly), followed by partitioning of SVOCs to wet aerosols or cloud/fog droplets, with subsequent aqueous-phase transformation to pyruvic acid (Pyr) and glyoxylic acid (ωC2) and finally oxidized to oxalic acid24,25,26. In this case, oxalic acid-C carries a depleted δ13C signal as the oxalic acid is an oxidation product, as opposed to a residual C pool. Our dataset displays a significant correlation of oxalic acid with ALW (r2 = 0.78, P < 0.001, Fig. 3b), suggesting that oxalic acid was dominantly formed via gas-to-liquid transfer of precursors and subsequent aqueous-phase reactions. This is consistent with a previous model study revealing that aqSOA formation is limited by ALW rather than by the availability of precursors39. As shown in Fig. 3b, the five points with excessive ALW are off the trend, signifying that ALW was no longer the limiting parameter in the specific samples. Oxalic acid and its aqueous-phase intermediate product, glyoxylic acid, also display a significant correlation (r2 = 0.65, P < 0.001), which stands for the aqueous-phase pathway.

The compound-specific δ13C values of oxalic acids and relevant SOA species can help to explore the formation pathway. First, the higher homologs aliphatic diacids from C3 to C9, in contrast to C2, show no significant statistical difference (t-test, P > 0.05, Supplementary Table 6) between the two air mass source regimes (Fig. 3a). Instead, ωC2 and Pyr show significantly lower (t-test, P < 0.05, Supplementary Table 6) δ13C values of (by 4.4‰ and 3.2‰, respectively, Fig. 3a) during continental outflow than coastal background. As aqSOA formed from VOCs or SVOCs typically leads to lower δ13C values relative to the precursors since lighter-isotope-containing precursors are preferentially oxidized to form reaction products. The continental outflow regime was characterized by high RH% and substantial anthropogenic hygroscopic particles, which could facilitate aqSOA formation. Therefore, the vast difference in the δ13C composition of oxalic acid between the two air mass source regimes could be mainly resulted from the aqueous-phase reaction pathway, i.e., AVOCs/BVOCs →SVOCs/WSOC → Pyr → ωC2 → C240, rather than gas-phase photochemical oxidation (or breakdown) processes. Note that the source mixing contribution is minimal as mentioned above.

This aqueous reaction effect on isotope fractionation may also apply to water-soluble OA, which is a cocktail of primary OA (e.g., sugars) and SOA species. Figure 3c reveals a clear decrease in δ13C signatures of oxalic acid against an increase in WSOC concentrations (r2 = 0.70, P < 0.001). Atmospheric aqueous-phase processing would simultaneously cause a decrease in the δ13C signature of reaction products and an increase in the WSOC mass concentrations. This is precisely what was observed during this campaign. Overall, the results indicate that aqSOA formation processes may be a significant route to form water-soluble OA in the atmosphere of South China and similar settings.

Besides δ13C signature, which undergoes strong isotope fractionation during atmospheric processing, radiocarbon isotopic fingerprint (Δ14C) is conservative and highly precise in distinguishing fossil from biogenic/biomass sources19. Individual water-soluble diacids and related compounds (C2, C3, C4, ωC2, and MeGly) were isolated by preparative capillary gas chromatography techniques (Methods). It is possible to harvest a sufficient mass of C2 oxalic acid from aerosol samples, while the other isolated compounds had to each be pooled per month to obtain sufficient carbon mass for the 14C measurement. The determined Δ14C compositions reported as the "fraction of modern carbon" (Fm) are shown in Supplementary Table 7 and Supplementary Table 8. It is further possible to translate the Fm signatures into fractional contribution from fossil (ffossil) versus non-fossil (including biogenic emissions: bio and biomass burning: bb; hereafter referred as fbio/bb) sources (Supplementary Method 3).

The excellent sensitivity and precision of molecular-level 14C analysis reveal a clear diversity in 14C signatures of oxalic acid for different samples and regimes (Fig. 2 and Supplementary Table 7). In the coastal background air mass regime, 67 ± 9% of oxalic acid is from biogenic emission and biomass burning. Although summer monsoon originating from SCS (Fig. 1) results in efficient precipitation, there is still a substantial fossil contribution (ffossil = 19–48%; 33 ± 9%) from anthropogenic activities. For the continental outflow samples, the ffossil−C2 (ffossil = 42–78%; 55 ± 10%) are even higher than biogenic/biomass source contributions. Back trajectory analysis reveals that during the latter period the receptor site was downwind from highly polluted continental source regions, which thus apparently is a source of oxalic acid in the regional aerosol (Fig. 1). The high ratios of NO3- in atmospheric particles (~10%; Supplementary Fig. 4) indicate a significant contribution of traffic emissions41 to atmospheric particles. Moreover, terephthalic acid (tPh) was one of the most abundant anthropogenic diacids accounting for 5% of total diacids (Supplementary Table 3). As tPh is a combustion tracer of plastic polymers42, this highlights the involvement of precursors from plastic burning to the formation of diacids. Thus, both traffic and plastic burning emissions were readily among the contributors to the fossil source-derived oxalic acid.

The mass ratio of adipic acid to azelaic acid (C6/C9) is often used to discriminate the relative significance of biogenic versus anthropogenic source contribution to diacids. Here, C6 adipic acid is the photochemical oxidation product of cyclic hexene emitted from the fossil fuel combustion sources27, 43, whereas C9 azelaic acid is a photochemical degradation product of biogenic unsaturated fatty acids33. Lower C6/C9 ratios (0.37 ± 0.13) were observed in the coastal background samples than continental outflow (0.76 ± 0.21). Furthermore, C6/C9 ratios correlates well with the parameter of ffossil−C2 (r2 = 0.63, P < 0.001; Supplementary Fig. 6). Hence, C6/C9 ratios gives strong and consistent support to the result of diacids 14C-based source apportionment.

The fossil contents of oxalic acid in the continental outflow samples are compared with its presumed precursors of longer-chain diacids (C3 and C4 diacids) and its aqueous processing precursors (ωC2 and MeGly) (Fig. 4c and Supplementary Table 8). The result indicates that fossil sources were the overwhelming contributor to both oxalic acid (55 ± 10%) and to its aqueous processing precursors, i.e., glyoxylic acid (ωC2; 69 ± 4%) and methylglyoxal (MeGly; 67 ± 5%). Therefore, a dominance from fossil fuel sources to their aqSOA formation in the continental outflow is obvious.

a Relationships between biogenic/biomass burning (bio/bb) source-derived oxalic acid vs. levoglucosan. b Relationships between bio/bb source-derived oxalic acid vs. biogenic secondary organic aerosol (SOA) tracer (include isoprene SOA tracers and α/β−pinene SOA tracers; Supplementary Method 1). c Monthly variability of non-fossil contribution to oxalic acid, malonic acid, succinic acid, glyoxylic acid, and methylglyoxal. The error bars represent ±1 standard deviation for total propagated uncertainty (Supplementary Method 3). d Relative contribution of biogenic emission, biomass burning, and fossil-fuel emission to oxalic acid in coastal background and continental outflow.

Although malonic C3 and succinic C4 acids are believed to be partly produced from aqueous processing22, 26, there is for these C3 and C4 a somewhat lower but still large contribution from fossil sources (C3: 31 ± 4%; C4: 36 ± 2%). This demonstrates that C3 and C4 diacids instead were mainly originated from the stepwise photochemical breakdown of biogenic unsaturated fatty acids even in highly polluted continental outflow air mass. The large contrast of fossil fraction among C2 diacids and its higher homologs (i.e., C3 and C4 diacids) further indicate that substantial fossil precursors contribute to the formation of C2 through aqueous-phase processing, which results in a much higher fossil content of C2 than of its higher homologs.

The 14C-derived concentrations of oxalic acid from non-fossil sources (bio/bb−C2) span from 40 to 417 ng m−3 (207 ± 88 ng m−3), whereas from fossil source (fossil−C2) are 21 to 650 ng m-3 (230 ± 168 ng m−3) (Supplementary Fig. 7). We observed a close similarity in the temporal variability of bio/bb−C2 with nss-K+ and levoglucosan (Supplementary Fig. 8), especially during fall to winter, when open burning of post-harvest agricultural crop residues is widespread in rural China44. For all data from the whole campaign, bio/bb−C2 concentration are positively correlated with levoglucosan (r2 = 0.67, P < 0.001; Fig. 4a). The ratios of bio/bb−C2 and nss-K+ also show significant correlation (r2 = 0.58, P < 0.001; Supplementary Fig. 9). In contrast, the coastal background samples record low concentrations of levoglucosan and nss-K+ (Fig. 4a and Supplementary Fig. 9), signifying an insignificant contribution from biomass burning. Further, a strong linear correlation (r2 = 0.81, P < 0.001, Fig. 4b) was found between bio/bb−C2 concentrations and biogenic SOA tracers (i.e., isoprene tracers and α/β−pinene tracers) for the coastal background samples, suggesting that the bio/bb−C2 in low biomass burning scenario was formed mainly from biogenic precursors. These independent lines of evidence are suggestive of and consistent with that background C2 were biogenic, while bb−C2 likely represent additional non-fossil C2 concentrations. We therefore assume that, when there is little biomass burning contribution (levoglucosan ≈ 0), the concentration of bio/bb−C2 (128 ng m-3; Fig. 4a) would approximately represent the background C2 level, that is, the bio−C2. With this estimation, the total concentration of oxalic acid can be quantitatively distributed into three sources, i.e., fossil−C2, bb−C2 and bio−C2.

The anthropogenic C2 (including fossil−C2 and bb−C2) accounted for 78% of total C2 in the continental outflow samples, whereas its contribution is 39% in the coastal background (Fig. 4d). Note that natural biomass burning (e.g., wildfires) was ignored since it contributed <4% of total biomass-burning-derived organic pollutants in China45. A tracer-based estimation suggests that anthropogenic gasSOA, which are formed by photochemical reaction in gas phase followed by thermodynamic partition to the condensed phase4, account for 79% of the total gasSOA in the PRD during fall-winter46. This anthropogenic content of gasSOA is consistent with our 14C-constrained source apportionment of aqSOA (fanthtopogenic−C2 = 78%). By extension, the contribution of anthropogenic emissions in aqSOA and gasSOA are likely to be consistent under ambient conditions with high ALW content. Our observation-based isotope constraints question the paradigm, largely based on model simulations, which instead have held that aqSOA, which may contribute as much mass as gasSOA to SOA budget, are largely formed from biogenic emissions4. Previous study considered that fossil precursors may be more hydrophobic than biogenic precursors, and thus are likely to form SOA preferentially through gas-phase processing9. However, the oxidation products of fossil precursors are largely water-soluble and thus their chemical reaction in the aqueous phase constitutes an unneglected pathway, which may be competitive to SOA formation via the gas-phase chemistry.

Figure 5 shows the dual effects of precursor sources and atmospheric processing of aqSOA by combining the 13C and 14C signatures of oxalic acid. Lower oxalic acid concentrations were found in the coastal background samples, with higher non-fossil contributions, and more enriched in 13C. Higher oxalic acid concentrations in the continental outflow samples are characterized by higher fossil contributions and more depleted in 13C. The hygroscopic particles and precursors in continental outflow air masses facilitate the aqueous-phase chemical processes, which results in lighter δ13C values and higher fossil contribution to oxalic acid.

The color of the symbols represents the oxalic acid concentration. The δ13C source signatures of oxalic acid precursor gases from biogenic volatile organic compounds (BVOCs; i.e., isoprene) and from anthropogenic volatile organic compounds (AVOCs; i.e., anthropogenic emitted non-methane hydrocarbons) are obtained from reported literature values (Supplementary Table 5). The two ovals include the coastal background samples and the majority of the continental outflow samples, respectively. The shadowed pink circle remarks the outliers of continental samples subjected to strong atmospheric aging.

Conversely, such an aqSOA formation process became less important during the coastal background period, when there are less significant anthropogenic activities, and hence less precursors and hygroscopic particles in the atmosphere. In this case, the photochemical degradation (aging) of oxalic acid would surpass its aqueous-phase formation, which is evident in the increasing δ13C values of oxalic acid (Fig. 5). The interpretation may also apply to several outliers of the continental outflow samples (pink shadow in Fig. 5), of which the OC/EC ratios are much lower than the others, referring to exhaustive atmospheric aging.

Many studies have investigated the formation of aqSOA from biogenic/biomass precursors (e.g., isoprene)3, 7, 8, 22, 47, 48, while much less attention has been paid to fossil/anthropogenic precursors (e.g., aromatic hydrocarbons). Fossil precursors have been thought to consist mostly of hydrophobic components, and thus fossil-derived SOA have been expected to be less soluble in water9. Smog chamber experiments show that water-soluble organic aerosols exhibited more similarities to biogenic SOA than fossil SOA49, 50. In fact, low fossil contributions to WSOC (0−29%) have been reported in Europe, North America, and South Asia across different seasons13,14,15,16,17,18, and 14C analyses of aerosol OC and its sub-fractions have been interpreted such that the majority of fossil-derived OA was water-insoluble51, 52. Hence, these earlier field and laboratory studies show that water-soluble OA is only to very limited extent stemmed from fossil sources. In contrast, the results from the current study, using an extensive set of molecular 13C/14C data in China, strongly point to anthropogenic-fossil emissions as an important source of water-soluble OA. This is also supported by a large contribution of fossil sources to WSOC (30−60%) reported for both Chinese urban areas52,53,54,55 and a receptor station in the southeast Yellow Sea intercepting a large regional footprint11, 12. Taken together, these two sets of findings imply a significant missing pathway responsible for the enhancement in fossil fuel-derived water-soluble OA in China.

The enhanced atmospheric aqueous-phase processing (δ13C of oxalic acid is used as the proxy for aqSOA formation) simultaneously cause an increasing in the WSOC mass concentration, suggesting that aqSOA constitute a substantial water-soluble OA source (Fig. 3c). A recent study also showed that the oxidation of gas-phase components and subsequent uptake of water-soluble oxidation products triggered by the increase of ALW can explain the increase in SOA mass in Chinese megacities32. This is further corroborated by radiocarbon signatures of individual aqSOA tracers (i.e., C2, ωC2, and MeGly), which reveals a dominant fossil source to aqSOA (Fig. 4c). The aqSOA even showed a fossil content (55−69%) comparable to the water-insoluble OC in other Chinese urban areas such as Beijing (62−68%) and Guangzhou (53−56%) in winter haze56. Our dual carbon isotope data of aqSOA tracers indicate that various fossil precursors may be oxidized to form aqSOA in the atmosphere and contribute substantially to the WSOC aerosols. We suggest that this atmospheric mechanism is responsible for the enhanced formation of fossil-derived WSOC in China.

To test the spatial extent of the findings, we retrieved aerosol samples from five emission hotspot megacities of China (Fig. 6a). The difference of δ13C and 14C-derived fbio/bb values between oxalic acid and bulk WSOC were compared in winter and in summer, respectively (Fig. 6b). Although oxalic acid accounts for a large fraction of WSOC (~5.2%), probably the single most abundant compound, we observed significant but opposite difference in the δ13C and Δ14C compositions between oxalic acid and WSOC during the two seasons. In winter, oxalic acid was more depleted in 13C and 14C than the WSOC pool in each of the cities, red-above-blue in Fig. 6b (refer to Supplementary Table 9). This agrees with the dominant role of aqueous-phase formation processes triggered by the increased ALW in winter32. As mentioned above, overwhelming fossil contributions to aqSOA compounds were observed, that were consistent with or even higher than the fossil content in water-insoluble OA reported in Chinese urban areas during winter haze52, 56. Substantial fossil-derived precursors are likely oxidized to WSOC aerosol through secondary aqueous processing in winter, resulting in products (such as oxalic acid) with more negative δ13C values and more fossil contributions than the bulk WSOC pool.

a Locations of the five megacities (Beijing, Guangzhou, Wuhan, Chengdu and Shanghai) and average fossil fuel CO2 (FF CO2) emissions during the year 2018. FF CO2 emissions indicating high levels of anthropogenic activities over the five megacities. The FF CO2 emission data is obtained from Open-source Data Inventory for Anthropogenic CO2 (https://db.cger.nies.go.jp/dataset/ODIAC/DL_odiac2020b.html). The coastline boundaries in the map are originated from Natural Earth free vector map data (https://www.naturalearthdata.com/). The administration boundaries in the map are originated from map products of National Geomatics Center of China (https://www.webmap.cn/). b The difference of δ13C and 14C-based fraction of non-fossil sources between oxalic acid and WSOC in winter (January 2018) and in summer (July 2018), respectively. The error bars for δ13C composition represent the relative standard deviation for repeat analyses (n = 3). The error bars for 14C-based fraction of non-fossil sources represent ±1 standard deviation for total propagated uncertainty (Supplementary Method 3).

On the contrary, in summer, oxalic acid was more enriched in 13C and 14C than the WSOC pool in Guangzhou, Wuhan and Shanghai. As shown in Fig. 6b, a general trend of blue-above-red is displayed, despite the very close values in Beijing and Chengdu which are within the analytical uncertainties. A probable interpretation is that biomass/biogenic component from the WSOC pool is preferentially subjected to oxidative aging to small molecules12, 57, which results in oxalic acid enriched both in δ13C and Δ14C. The differing isotopic signatures of different emission sources would also contribute to the enrichment of carbon isotopes of oxalic acid. Note here that biogenic contributions to oxalic acid can be expected to increase during summer due to it being the growing season.

The large contribution of fossil-derived carbon to aqSOA highlighted in this study has important implications for aerosol climate forcing and regional air quality. Aqueous-phase processing of fossil precursors, competitive to gas-phase pathway, could amplify the impact of anthropogenic emissions on water-soluble OA, alter the chemical property and hygroscopicity of aerosols. The derived hygroscopic SOA affect not only direct radiative forcing, but also cloud condensation nuclei (CCN) activity, and the subsequent cloud-mediated effect of anthropogenic aerosols on the climate system58. As for oxalic acid, its interactions with inorganic salts in the aqueous phase may modify the CCN activity of aerosol particle59,60,61.

A mutual promotion effect may take place between ALW and inorganic particles (e.g., sulfate and nitrate)31, 62. Our year-round observation at the Heshan receptor site witnesses a positive feedback loop among ALW, inorganics, as well as aqSOA compounds. The dataset supports that the inorganic particles induce ALW growth (Supplementary Fig. 5), facilitating the partitioning of gas-phase oxidation products into aqueous-phase medium, and promoting the formation of aqSOA (Fig. 3b)6, 30. The formation of aqSOA in the organic aerosol fraction would decrease aerosol viscosity and increase ALW63, as in turn promotes further gas-to-liquid transfer of water-soluble organic precursors23, 32. Oxalic acid is identified as the upper bound of the hygroscopic properties of OA, with a k-Kohler of 0.487. Using the Zdanovskii−Stokes−Robinson (ZSR) mixing model64, the effect of anthropogenic oxalic acid on aerosol water accounts for 3−38% (average: 10%) of the ALW contributed by organic compounds (Supplementary Method 2). The inorganic particles can be considered anthropogenic in most continental atmosphere3, 30, 39. We provide 14C-based evidence that, besides inorganic species, substantial aqSOA compounds are also derived from anthropogenic precursors. Therefore, a broad control of various anthropogenic emissions such as SO2, NOx, and VOC precursors is vital for reducing the organic particulate pollution.

Nitrate is the dominant atmospheric hygroscopic particles in China30, 31, 41. A model study predicted that 92% of increased ALW in eastern China by the year 2100 can be attributed to the increase in anthropogenic nitrate aerosol30. Here, the average nitrate concentration in the continental outflow samples was nearly 20 times higher than coastal background (Supplementary Table 2). Thus, this nitrate enhancement in the particulate phase could be the dominant factor that facilitated the increase of ALW in the continental outflow regime. Furthermore, there has been a rapid translation from coal combustion to natural gas in China. However, natural gas combustion can produce more than three times as much water vapor as coal burning65. The combustion-derived water was interpreted to constitute 6.2% of the atmospheric moisture and added 5.1% of the anthropogenic PM2.5 in northwestern China65. Therefore, the enhancing effects of particulate nitrate and combustion-derived water vapor on aqSOA formation should be addressed when changing energy structure in future climate and air quality scenarios.

The present findings highlight the relevance of aqueous-phase chemistry for processing of fossil-fuel emissions over China. These scenarios are not unique to East Asia and could be relevant in other high fossil-fuel consumption regions with humid weather. In a study in northern Georgia of the US, Weber et al.13 observed strong correlations between WSOC and anthropogenic precursors, while in contrast the radiocarbon analysis revealed dominant biogenic contributions to WSOC. Our results suggest that aqSOA do not necessarily share the same source patterns with bulk WSOC. WSOC have a large primary contribution from biomass burning66, which would mask the fossil content in aqSOA. Therefore, the fossil contributions to aqSOA would be much higher than expected in locations such as North America, Europe, and South Asia where the aqueous-phase chemistry pathway has been demonstrated significant7, 39, 67. Even in pristine environments where biogenic emissions are dominant, the importance of anthropogenic emissions to aqSOA formation would also be significant since the ALW are largely anthropogenic2, 39.

In this study, the δ13C differences of oxalic acid was dominantly attributed to kinetic isotope effects in atmospheric processing. It should be noted that the differences in the isotopic signatures of emission sources may also contribute to the isotopic variations of oxalic acid. However, these factors are not quantitively constrained. The extent to which the atmospheric processing is representative for the differences in isotopic signatures remains to be investigated in more details. We also stress that apart from aqueous-phase processes, both biomass burning and gas-phase photochemical aging could also be sources for oxalic acid. To differentiating these pathways in field measurements is still challenging. Therefore, cautious interpretations of oxalic acid isotopic signals are needed where there are intensive biomass burning, extensive atmospheric aging, and seasonality of C3/C4 vegetation changes. Future dual-carbon isotope studies on other aqSOA constituents (e.g., glyoxal) are strongly warranted to better constrain the chemical mechanism responsible for fossil-derived aqSOA formation.

Overall, compound-specific dual carbon isotopic evidences from field observations show that fossil anthropogenic precursors contribute to the formation of aqSOA to a much higher extent than expected. Atmospheric aqueous-phase processes are expected to increase due to the enhanced evapotranspiration in a warmer world64. An absence of accounting for such aqueous processing incurred by fossil precursors could lead to an underestimation of the anthropogenic contribution to organic aerosols. Understanding the role of anthropogenic emissions in aqSOA formation and the chemical mechanisms involved are of critical concern for future climate and air quality projections, and within the context of energy needs and choices.

The regional receptor site was positioned in Heshan Atmospheric Environmental Monitoring Superstation (22.711° N, 112.927° E, 60 m asl), a rural site located at 50 km southwest to the megacity of Guangzhou in the central PRD. The surrounding area of the sampling site is dominated by farmlands and forests and is far from local anthropogenic emissions68. During the winter northeast-monsoon season, this specific site well intercepts high anthropogenically dominated outflow airmass from Chinese continental. The strength of anthropogenic impact decreases in the summer southeast/southwest monsoon, with the arrival of clean air mass originating from the Western Pacific/Indian Ocean. Field experiments started in June 2017 and ended in May 2018. PM2.5 aerosol samples (n = 32) were collected onto pre-combusted quartz-fiber filters by using a high-volume sampler at an air-flow rate of 1 m3 /min for 48 h.

PM2.5 samples were additionally studied in five emission hot spot megacities of China (Beijing, Shanghai, Guangzhou, Chengdu, Wuhan), which were collected continuously every 24 h for 1 week in winter (January 2018) and in summer (July 2018). Continuous sets of 7 filters collected in one week were combined to represent the seasonal characteristic of each urban site. The sample filters were stored in darkness at −20 °C prior to analysis. Field blank filters were collected for each campaign.

Data of atmospheric concentration of meteorological parameters and gaseous pollutants were obtained from the hourly monitor by Heshan Atmospheric Environmental Monitoring Superstation. To characterize different geographical sources of air masses encountered at Heshan, 3-day isentropic backward air mass trajectories were computed at an elevation of 100 m for every 12 h using the HYSPLIT model (version-4) and meteorological datasets from the NOAA air resources laboratory (http://ready.arl.noaa.gov/HYSPLIT.php). The back trajectories are grouped into two major transport pathways with vastly different source regions during the sampling campaigns.

The concentration of typical aqSOA species (diacids, oxoacids, and α-dicarbonyls) was measured following earlier protocols19, 27. An aliquot of filter samples was extracted with Milli-Q water. The extracts were concentrated to dryness and derivatized with 10% BF3 in 1-butanol (Sigma-Aldrich, St Louis, MO) in a hot water bath at 100 °C for about 1 h. The derivatives were extracted with n-hexane and finally quantified using GC-MS (Agilent 7890 A gas chromatography and Agilent 5975C mass spectrometry). The concentrations reported here are corrected for field blank, but are not corrected for recovery efficiencies.

Relevant details of the analysis method for aerosol carbonaceous components (organic carbon: OC, elemental carbon: EC; water-soluble organic carbon: WSOC), water-soluble inorganic constituents (WSIC, i.e., Na+, NH4+, K+, Mg2+, Ca2+, Cl−, NO3−, and SO42−), biogenic SOA tracers; and sugar compounds are provided in the Supplementary Method 1.

The total aerosol liquid water (ALW) is taken as the sum of water associated with individual aerosol chemical components, with the assumption that the particles are internally mixed69. ALW from inorganic species was estimated using a thermodynamic equilibrium model, ISORROPIA-II; http://nenes.eas.gatech.edu/ISORROPIA/ (Supplementary Method 2) using particle concentration of inorganic ions and meteorological parameters (temperature and RH%). The organic contribution for ALW was calculated by Zdanovskii−Stokes−Robinson (ZSR) mixing rule as discussed in Supplementary Method 2.

The δ13C values of diacids derivatives were determined using GC ISOLink2/IRMS (Thermo Fisher Scientific). The δ13C values of free diacids and related compounds were then calculated by an isotopic mass balance approach, based on the measured δ13C values of derivatives and the derivatizing agent (1-butanol, δ13C = −30.21‰). The isotopic fractionation during the derivatization step has been reported to be less than 0.73‰ in our earlier publications19. Each sample was measured in duplicate and the mean δ13C values were reported. The difference in δ13C of diacids (i.e., C2, C3, C4, C5, C6, Ph, tPh, C9) in replicate analyses was generally <1‰, though the analytical accuracy was < 2‰ for other compounds (i.e., Pyr, ωC2).

The analytical method for the isolation of individual diacids to determine Δ14C compositions builds on our previously publications19. In brief, the filter was extracted with the Milli-Q water and derivatized with BF3/1-butanol. Microgram quantities of C2, C3, C4, ωC2, and MeGly were isolated and collected by preparative capillary gas chromatography (pcGC) through around 50 consecutive runs (5 μl per injection), which was sufficient for off-line natural abundance radiocarbon analysis. The pcGC isolates were combusted at 920 °C and the resulting CO2 was finally converted to graphite via hydrogen reduction method. Measurements of 14C were performed using Accelerator Mass Spectrometry (1.5SDH-1, 0.5MV, NEC, USA)19. In all cases, 14C results are reported as "fraction of modern carbon" (Fm) normalized to a common δ13C value of −25‰. The Fm values were further converted into the fraction of non-fossil source (fbio/bb) with a conversion factor of 1.06 to compensate for the excess 14C produced by nuclear bomb testing in the 1950s–1960s. The concentration of diacids in blank filters was <0.06% of the real samples, therefore, no filter blank subtraction was performed for the isotope analysis in this study. Sample results were corrected for procedural blanks using a pair of processing standards with modern or fossil radiocarbon composition. An isotopic mass-balance approach was adopted to correct the carbon contribution of butanol groups (1-butanol, Fm = 0.0029 ± 0.001) added in the derivatives of diacids, where appropriate. Calculations for blank correction, Fm results and error propagation are shown in Supplementary Method 3.

The analytical method for the measurements of δ13C and Δ14C composition of WSOC builds on ultrasonication extraction protocols were reported in Supplementary Method 4. The validity of the ultrasonication method during extraction of water-soluble organic materials was demonstrated by comparing ultrasonication with another extraction protocol, i.e., soaking (see Supplementary Table 10 and Supplementary Fig. 10). Measurements of the δ13C composition of WSOC were performed on a Flash 2000 elemental analyzer connected to a Thermo Scientific Delta V isotope ratio mass spectrometer. Measurements of the Δ14C composition were performed using accelerator mass spectrometry facility (1.5SDH-1, 0.5MV, NEC, USA) of the Guangzhou Institute of Geochemistry of the Chinese Academy of Sciences (GIGCAS).

The data supporting the findings of this study are available at the Figshare digital repository (https://doi.org/10.6084/m9.figshare.20469540). The aerosol optical depth data are available through the Moderate Resolution Imaging Spectroradiometer (MODIS) observations (https://giovanni.gsfc.nasa.gov/giovanni/). The isentropic backward air mass trajectories are available through the NOAA air resources laboratory (http://ready.arl.noaa.gov/HYSPLIT.php). The fossil fuel CO2 emission data were accessed through the Open-source Data Inventory for Anthropogenic CO2 (https://db.cger.nies.go.jp/dataset/ODIAC/DL_odiac2020b.html). Source data are provided with this paper.

Jimenez, J. L. et al. Evolution of organic aerosols in the atmosphere. Science 326, 1525–1529 (2009).

Article ADS CAS PubMed Google Scholar

Carlton, A. G., Christiansen, A. E., Flesch, M. M., Hennigan, C. J. & Sareen, N. Mulitphase atmospheric chemistry in liquid water: impacts and controllability of organic aerosol. Acc. Chem. Res. 53, 1715–1723 (2020).

Article CAS PubMed Google Scholar

McNeill, V. F. Aqueous organic chemistry in the atmosphere: sources and chemical processing of organic aerosols. Environ. Sci. Technol. 49, 1237–1244 (2015).

Article ADS CAS PubMed Google Scholar

Ervens, B., Turpin, B. J. & Weber, R. J. Secondary organic aerosol formation in cloud droplets and aqueous particles (aqSOA): a review of laboratory, field and model studies. Atmos. Chem. Phys. 11, 11069–11102 (2011).

Article ADS CAS Google Scholar

Lim, Y. B., Tan, Y., Perri, M. J., Seitzinger, S. P. & Turpin, B. J. Aqueous chemistry and its role in secondary organic aerosol (SOA) formation. Atmos. Chem. Phys. 10, 10521–10539 (2010).

Article ADS CAS Google Scholar

Hennigan C. J., Bergin M. H., Dibb J. E., Weber R. J. Enhanced secondary organic aerosol formation due to water uptake by fine particles. Geophys. Res. Lett. 35, L18801 (2008).

Gilardoni, S. et al. Direct observation of aqueous secondary organic aerosol from biomass-burning emissions. Proc. Natl Acad. Sci. USA 113, 10013–10018 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Lamkaddam, H. et al. Large contribution to secondary organic aerosol from isoprene cloud chemistry. Sci. Adv. 7, eabe2952 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Favez O., Sciare J., Cachier H., Alfaro S. C., Abdelwahab M. M. Significant formation of water-insoluble secondary organic aerosols in semi-arid urban environment. Geophys. Res. Lett. 35, L15801 (2008).

Wang, J. et al. Aqueous production of secondary organic aerosol from fossil-fuel emissions in winter Beijing haze. Proc. Natl Acad. Sci. USA 118, e2022179118 (2021).

Article CAS PubMed PubMed Central Google Scholar

Fang, W. et al. Divergent evolution of carbonaceous aerosols during dispersal of East Asian Haze. Sci. Rep. 7, 10422 (2017).

Article ADS PubMed PubMed Central CAS Google Scholar

Kirillova, E. N., Andersson, A., Han, J., Lee, M. & Gustafsson, Ö. Sources and light absorption of water-soluble organic carbon aerosols in the outflow from northern China. Atmos. Chem. Phys. 14, 1413–1422 (2014).

Article ADS CAS Google Scholar

Weber, R. J. et al. A study of secondary organic aerosol formation in the anthropogenic-influenced southeastern United States. J. Geophys. Res. Atmos. 112, D13302 (2007).

Zhang, Y. L. et al. Fossil and non-fossil sources of different carbonaceous fractions in fine and coarse particles by radiocarbon measurement. Radiocarbon 55, 1510–1520 (2013).

Article CAS Google Scholar

Kirillova, E. N. et al. 13C- and14C-based study of sources and atmospheric processing of water-soluble organic carbon (WSOC) in South Asian aerosols. J. Geophys. Res.: Atmos. 118, 614–626 (2013).

Article ADS CAS Google Scholar

Szidat, S. et al. Fossil and non-fossil sources of organic carbon (OC) and elemental carbon (EC) in Goteborg, Sweden. Atmos. Chem. Phys. 9, 1521–1535 (2009).

Article ADS CAS Google Scholar

Bosch, C. et al. Source-diagnostic dual-isotope composition and optical properties of water-soluble organic carbon and elemental carbon in the South Asian outflow intercepted over the Indian Ocean. J. Geophys. Res.: Atmos. 119, 11,743–711,759 (2014).

Article CAS Google Scholar

Wozniak, A. S., Bauer, J. E. & Dickhut, R. M. Characteristics of water-soluble organic carbon associated with aerosol particles in the eastern United States. Atmos. Environ. 46, 181–188 (2012).

Article ADS CAS Google Scholar

Xu, B. et al. Compound-specific radiocarbon analysis of low molecular weight dicarboxylic acids in ambient aerosols using preparative gas chromatography: method development. Environ. Sci. Technol. Lett. 8, 135–141 (2021).

Article CAS Google Scholar

Cao, F. et al. New insights into the sources and formation of carbonaceous aerosols in China: potential applications of dual-carbon isotopes. Natl Sci. Rev. 4, 804–+ (2017).

Article CAS Google Scholar

Sorooshian, A. et al. Particulate organic acids and overall water-soluble aerosol composition measurements from the 2006 Gulf of Mexico Atmospheric Composition and Climate Study (GoMACCS). J. Geophys. Res.: Atmos. 112, D13201 (2007).

Chen, Y. et al. Low-molecular-weight carboxylic acids in the Southeastern U.S.: Formation, partitioning, and implications for organic aerosol aging. Environ. Sci. Technol. 55, 6688–6699 (2021).

Article ADS CAS PubMed Google Scholar

Wang, J. et al. Enhanced aqueous-phase formation of secondary organic aerosols due to the regional biomass burning over North China Plain. Environ. Pollut. 256, 113401 (2020).

Article CAS PubMed Google Scholar

Carlton, A. G. et al. Atmospheric oxalic acid and SOA production from glyoxal: Results of aqueous photooxidation experiments. Atmos. Environ. 41, 7588–7602 (2007).

Article ADS CAS Google Scholar

Lim, Y. B., Tan, Y. & Turpin, B. J. Chemical insights, explicit chemistry, and yields of secondary organic aerosol from OH radical oxidation of methylglyoxal and glyoxal in the aqueous phase. Atmos. Chem. Phys. 13, 8651–8667 (2013).

Article ADS CAS Google Scholar

Bikkina, S., Kawamura, K. & Sarin, M. Secondary organic aerosol formation over coastal ocean: inferences from atmospheric water-soluble low molecular weight organic compounds. Environ. Sci. Technol. 51, 4347–4357 (2017).

Article ADS CAS PubMed Google Scholar

Kawamura, K. & Ikushima, K. Seasonal changes in the distribution of dicarboxylic acids in the urban atmosphere. Environ. Sci. Technol. 27, 2227–2235 (1993).

Article ADS CAS Google Scholar

Kawamura, K. & Sakaguchi, F. Molecular distributions of water-soluble dicarboxylic acids in marine aerosols over the Pacific Ocean including tropics. J. Geophys. Res.: Atmos. 104, 3501–3509 (1999).

Article ADS CAS Google Scholar

Dasari, S. et al. Photochemical degradation affects the light absorption of water-soluble brown carbon in the South Asian outflow. Sci. Adv. 5, eaau8066 (2019).

Article ADS PubMed PubMed Central CAS Google Scholar

Hodas, N. et al. Aerosol liquid water driven by anthropogenic nitrate: implications for lifetimes of water-soluble organic gases and potential for secondary organic aerosol formation. Environ. Sci. Technol. 48, 11127–11136 (2014).

Article ADS CAS PubMed Google Scholar

Wu, Z. et al. Aerosol liquid water driven by anthropogenic inorganic salts: implying its key role in haze formation over the North China Plain. Environ. Sci. Technol. Lett. 5, 160–166 (2018).

Article CAS Google Scholar

Gkatzelis, G. I. et al. Uptake of water-soluble gas-phase oxidation products drives organic particulate pollution in Beijing. Geophys. Res. Lett. 48, e2020GL091351 (2021).

Article ADS Google Scholar

Kawamura, K. & Gagosian, R. B. Implications of ω-oxocarboxylic acids in the remote marine atmosphere for photo-oxidation of unsaturated fatty acids. Nature 325, 330–332 (1987).

Article ADS CAS Google Scholar

Ballentine, D. C., Macko, S. A. & Turekian, V. C. Variability of stable carbon isotopic compositions in individual fatty acids from combustion of C4 and C3 plants: implications for biomass burning. Chem. Geol. 152, 151–161 (1998).

Article ADS CAS Google Scholar

Miyazaki, Y., Kawamura, K., Jung, J., Furutani, H. & Uematsu, M. Latitudinal distributions of organic nitrogen and organic carbon in marine aerosols over the western North Pacific. Atmos. Chem. Phys. 11, 3037–3049 (2011).

Article ADS CAS Google Scholar

Martinelli, L. A., Camargo, P. B., Lara, L., Victoria, R. L. & Artaxo, P. Stable carbon and nitrogen isotopic composition of bulk aerosol particles in a C4 plant landscape of southeast Brazil. Atmos. Environ. 36, 2427–2432 (2002).

Article ADS CAS Google Scholar

Aggarwal, S. G. & Kawamura K. Molecular distributions and stable carbon isotopic compositions of dicarboxylic acids and related compounds in aerosols from Sapporo, Japan: Implications for photochemical aging during long-range atmospheric transport. J. Geophys. Res.: Atmos. 113, D14301 (2008).

Bikkina, S., Kawamura, K., Sakamoto, Y. & Hirokawa, J. Low molecular weight dicarboxylic acids, oxocarboxylic acids and α-dicarbonyls as ozonolysis products of isoprene: Implication for the gaseous-phase formation of secondary organic aerosols. Sci. Total Environ. 769, 144472 (2021).

Article ADS CAS PubMed Google Scholar

Carlton, A. G. & Turpin, B. J. Particle partitioning potential of organic compounds is highest in the Eastern US and driven by anthropogenic water. Atmos. Chem. Phys. 13, 10203–10214 (2013).

Article ADS CAS Google Scholar

Zhang, Y.-L., Kawamura, K., Cao, F. & Lee, M. Stable carbon isotopic compositions of low-molecular-weight dicarboxylic acids, oxocarboxylic acids,α-dicarbonyls, and fatty acids: Implications for atmospheric processing of organic aerosols. J. Geophys. Res.: Atmos. 121, 3707–3717 (2016).

Article ADS CAS Google Scholar

Zhang, Q. et al. Drivers of improved PM2.5 air quality in China from 2013 to 2017. Proc. Natl Acad. Sci. USA 116, 24463–24469 (2019).

Article ADS CAS PubMed PubMed Central Google Scholar

Pavuluri C. M., Kawamura K., Swaminathan T. Water-soluble organic carbon, dicarboxylic acids, ketoacids, andα-dicarbonyls in the tropical Indian aerosols. J. Geophys. Res. Atmos 115, D11302 (2010).

Hatakeyama, S. et al. Ozone-cyclohexene reaction in air: quantitative analysis of particulate products and the reaction mechanism. Environ. Sci. Technol. 19, 935–942 (1985).

Article ADS CAS Google Scholar

Yuan, B., Liu, Y., Shao, M., Lu, S. & Streets, D. G. Biomass burning contributions to ambient VOCs species at a receptor site in the Pearl River Delta (PRD), China. Environ. Sci. Technol. 44, 4577–4582 (2010).

Article ADS CAS PubMed Google Scholar

Zhou, Y. et al. A comprehensive biomass burning emission inventory with high spatial and temporal resolution in China. Atmos. Chem. Phys. 17, 2839–2864 (2017).

Article ADS CAS Google Scholar

Ding, X. et al. Tracer-based estimation of secondary organic carbon in the Pearl River Delta, south China. J. Geophys. Res.: Atmos. 117, D05313 (2012).

Lim, H.-J., Carlton, A. G. & Turpin, B. J. Isoprene forms secondary organic aerosol through cloud processing: Model simulations. Environ. Sci. Technol. 39, 4441–4446 (2005).

Article ADS CAS PubMed Google Scholar

Carlton, A. G. et al. CMAQ model performance enhanced when in-cloud secondary organic aerosol is included: comparisons of organic carbon predictions with measurements. Environ. Sci. Technol. 42, 8798–8802 (2008).

Article ADS CAS PubMed Google Scholar

Baltensperger, U. et al. Secondary organic aerosols from anthropogenic and biogenic precursors. Faraday Discuss. 130, 265–278 (2005).

Article ADS CAS PubMed Google Scholar

Updyke, K. M., Nguyen, T. B. & Nizkorodov, S. A. Formation of brown carbon via reactions of ammonia with secondary organic aerosols from biogenic and anthropogenic precursors. Atmos. Environ. 63, 22–31 (2012).

Article ADS CAS Google Scholar

Wozniak A. S., Bauer J. E., Dickhut R. M., Xu L., McNichol A. P. Isotopic characterization of aerosol organic carbon components over the eastern United States. J. Geophys. Res. Atmos. 117, D13303 (2012).

Zhang, Y.-L. et al. Large contribution of fossil fuel-derived secondary organic carbon to water soluble organic aerosols in winter haze in China. Atmos. Chem. Phys. 18, 4005–4017 (2018).

Article ADS CAS Google Scholar

Mo, Y. et al. Dual carbon isotope-based source apportionment and light absorption properties of water-soluble organic carbon in PM2.5 over China. J. Geophys. Res.: Atmos., 126, e2020JD033920 (2021).

Liu, J. et al. Source apportionment and dynamic changes of carbonaceous aerosols during the haze bloom-decay process in China based on radiocarbon and organic molecular tracers. Atmos. Chem. Phys. 16, 2985–2996 (2016).

Article ADS CAS Google Scholar

Liu, J. et al. Radiocarbon-derived source apportionment of fine carbonaceous aerosols before, during, and after the 2014 Asia-Pacific Economic Cooperation (APEC) summit in Beijing, China. J. Geophys. Res.: Atmos. 121, 4177–4187 (2016).

Article ADS CAS Google Scholar

Zhang, Y.-l et al. Micro-scale (μg) radiocarbon analysis of water-soluble organic carbon in aerosol samples. Atmos. Environ. 97, 1–5 (2014).

Article ADS CAS Google Scholar

Kirillova, E. N. et al. Water-soluble organic carbon aerosols during a full New Delhi winter: Isotope-based source apportionment and optical properties. J. Geophys. Res.: Atmos. 119, 3476–3485 (2014).

Article ADS CAS Google Scholar

Glassmeier, F. et al. Aerosol-cloud-climate cooling overestimated by ship-track data. Science 371, 485–489 (2021).

Article ADS CAS PubMed Google Scholar

Prenni, A. J. et al. The Effects of low molecular weight dicarboxylic acids on cloud formation. J. Phys. Chem. A 105, 11240–11248 (2001).

Article CAS Google Scholar

Jing, B. et al. Hygroscopic behavior of multicomponent organic aerosols and their internal mixtures with ammonium sulfate. Atmos. Chem. Phys. 16, 4101–4118 (2016).

Article ADS CAS Google Scholar

Boreddy, S. K. R. & Kawamura, K. Investigation on the hygroscopicity of oxalic acid and atmospherically relevant oxalate salts under sub- and supersaturated conditions. Environ. Sci.: Process. Impacts 20, 1069–1080 (2018).

CAS Google Scholar

Wang, Y. et al. Mutual promotion between aerosol particle liquid water and particulate nitrate enhancement leads to severe nitrate-dominated particulate matter pollution and low visibility. Atmos. Chem. Phys. 20, 2161–2175 (2020).

Article ADS CAS Google Scholar

Shiraiwa, M., Ammann, M., Koop, T. & Pöschl, U. Gas uptake and chemical aging of semisolid organic aerosol particles. Proc. Natl Acad. Sci. USA 108, 11003–11008 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Nguyen, T. K. V., Zhang, Q., Jimenez, J. L., Pike, M. & Carlton, A. G. Liquid water: ubiquitous contributor to aerosol mass. Environ. Sci. Technol. Lett. 3, 257–263 (2016).

Article CAS Google Scholar

Xing, M. et al. Vapor isotopic evidence for the worsening of winter air quality by anthropogenic combustion-derived water. Proc. Natl Acad. Sci. USA 117, 33005–33010 (2020).

Article ADS CAS PubMed Central Google Scholar

Mayol-Bracero, O. L. et al. Water-soluble organic compounds in biomass burning aerosols over Amazonia 2. Apportionment of the chemical composition and importance of the polyacidic fraction. J. Geophys. Res. Atmos. 107, LBA 59-1–LBA 59-15 (2002).

Article CAS Google Scholar

Chen, Y. et al. Ammonium chloride associated aerosol liquid water enhances daze in Delhi, India. Environ. Sci. Technol. 56, 7163–7173 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Kong, L. et al. Aerosol optical properties under different pollution levels in the Pearl River Delta (PRD) region of China. J. Environ. Sci. 87, 49–59 (2020).

Article Google Scholar

Guo, H. et al. Fine-particle water and pH in the southeastern United States. Atmos. Chem. Phys. 15, 5211–5228 (2015).

Article ADS CAS Google Scholar

Download references

This study was funded by the Natural Science Foundation of China (42030715), the Alliance of International Science Organizations research and cooperative projects (ANSO-CR-KP−2021-05), the Guangdong Foundation for Program of Science and Technology Research (2017BT01Z134), and the National Key R&D Program of China (2017YFC0212000). We would like to thank Dr Duohong Chen for providing the meteorological data, Dr Weiwei Hu for helpful discussions, and Dr Lili Ming for support with instruments.

Srinivas Bikkina

Present address: CSIR-National Institute of Oceanography, Dona Paula, 403004, Goa, India

Ambarish Pokhrel

Present address: Institute of Science and Technology, Tribhuvan University, Kathmandu, 44600, Nepal

State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou, 510640, China

Buqing Xu, Gan Zhang, Jun Li, Guangcai Zhong, Shizhen Zhao, Jing Li, Chen Huang, Zhineng Cheng, Sanyuan Zhu, Pingan Peng & Guoying Sheng

CAS Center for Excellence in Deep Earth Science, Guangzhou, 510640, China

Buqing Xu, Gan Zhang, Jun Li, Guangcai Zhong, Shizhen Zhao, Jing Li, Chen Huang, Zhineng Cheng, Sanyuan Zhu, Pingan Peng & Guoying Sheng

Department of Environment Science and the Bolin Centre for Climate Research, Stockholm University, Stockholm, 10691, Sweden

Örjan Gustafsson & August Andersson

Chubu Institute for Advanced Studies, Chubu University, Kasugai, 487-8501, Japan

Kimitaka Kawamura, Srinivas Bikkina, Bhagawati Kunwar & Ambarish Pokhrel

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

B.X. and G.Z. designed the experiment. JUN.L., Z.C., and S.-Z.Z. provided the samples. B.X., A.A., B.K., A.P., JING.L., C.H., Z.C. and S.-Y.Z. processed data and performed analyses. B.X. wrote the paper. G.Z., Ö.G., K.K., A.A., S.B., G.-C.Z., P.P., and G.S. commented on the manuscript.

Correspondence to Gan Zhang or Örjan Gustafsson.

The authors declare no competing interests.

Nature Communications thanks Andrew Wozniak and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Xu, B., Zhang, G., Gustafsson, Ö. et al. Large contribution of fossil-derived components to aqueous secondary organic aerosols in China. Nat Commun 13, 5115 (2022). https://doi.org/10.1038/s41467-022-32863-3

Download citation

Received: 09 December 2021

Accepted: 22 August 2022

Published: 31 August 2022

DOI: https://doi.org/10.1038/s41467-022-32863-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Chemical Research in Chinese Universities (2023)

Advances in Atmospheric Sciences (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.