Start Submission Become a Reviewer

Reading: Source appointment of nitrogen in PM2.5 based on bulk δ15N signatures and a Bayesian isotope...

Download

A- A+
Alt. Display

Original Research Papers

Source appointment of nitrogen in PM2.5 based on bulk δ15N signatures and a Bayesian isotope mixing model

Authors:

Yan-Li Wang,

State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences; Chinese Academy for Environmental Planning, CN
X close

Xue-Yan Liu ,

Institute of Surface-Earth System Science, Tianjin University; State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences, CN
X close

Wei Song,

Institute of Surface-Earth System Science, Tianjin University, CN
X close

Wen Yang,

State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, CN
X close

Bin Han,

State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, CN
X close

Xiao-Yan Dou,

Qinghai Environmental Monitoring Center, CN
X close

Xu-Dong Zhao,

Qinghai Environmental Monitoring Center, CN
X close

Zhao-Liang Song,

Institute of Surface-Earth System Science, Tianjin University; State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences, CN
X close

Cong-Qiang Liu,

Institute of Surface-Earth System Science, Tianjin University; State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences, CN
X close

Zhi-Peng Bai

State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, CN
X close

Abstract

Nitrogen isotope (δ15N) has been employed to differentiate major sources of atmospheric N. However, it remains a challenge to quantify contributions of multiple sources based on δ15N values of the N mixture in atmospheric samples. This study measured δ15N of bulk N in PM2.5 at an urban site of Beijing during a severe haze episode of 22–30 January 2013 and a background site of Qinghai, north-western China from 6 September to 15 October 2013, then applied a Bayesian isotope mixing model (SIAR, Stable Isotope Analysis in R) to analyse the N sources. At Beijing site, δ15N values of PM2.5 (−4.1‰ to +13.5‰, +2.8 ± 6.4‰) were distributed within the range of major anthropogenic sources (including NH3 and NO2 from coal combustion, vehicle exhausts and domestic wastes/sewage). At Menyuan site, δ15N values of PM2.5 (+8.0‰ to +27.9‰, +18.5 ± 5.8‰) were significantly higher than that of potential sources (including NH3 and NO2 from biomass burning, animal wastes, soil N cycle, fertilizer application and dust N). High molar ratios of NH4+ to NO3- and/or SO42- in PM2.5 at the background site suggested that the equilibrium of NH3 ↔ NH4+ caused apparent 15N enrichments in ammonium. Results of the SIAR model showed that 39 and 32% of bulk N in PM2.5 of Beijing site were contributed from N emissions of coal combustion and vehicle exhausts, respectively, whereas N in PM2.5 at Menyuan site was derived mainly from N emissions of biomass burning (46%) and NH3 volatilization (34%). These results revealed that the stoichiometry between NH3 and acidic gases plays an important role in controlling the bulk δ15N signatures of PM2.5 and emissions of reactive N from coal combustion and urban transportation should be strictly controlled to advert the risk of haze episodes in Beijing.

How to Cite: Wang, Y.-L., Liu, X.-Y., Song, W., Yang, W., Han, B., Dou, X.-Y., Zhao, X.-D., Song, Z.-L., Liu, C.-Q. and Bai, Z.-P., 2017. Source appointment of nitrogen in PM2.5 based on bulk δ15N signatures and a Bayesian isotope mixing model. Tellus B: Chemical and Physical Meteorology, 69(1), p.1299672. DOI: http://doi.org/10.1080/16000889.2017.1299672
6
Views
3
Downloads
  Published on 01 Jan 2017
 Accepted on 20 Feb 2017            Submitted on 20 Jan 2017
1.  

Introduction

Urban air pollution is a globally challenging issue. Nitrogen (N) emissions play a key role in the formation of atmospheric particulates, especially secondary N-containing aerosols (Huang et al., 2014). Ammonia is the precursor of ammonium (NH4+) and readily reacts with available acids formed by SO2 and NO2, and also it can be transformed to organic N or amines (Ge et al., 2011). Nitrogen oxides are major precursors of both inorganic (as nitrate ions (NO3-)) and organic (as organic NO3-) N aerosols (Berkemeier et al., 2016). Therefore, the source apportionment of N in PM2.5 is always of significance for better understanding origins of particulates and haze pollution (Guo et al., 2014).

Stable isotopes of N (i.e. δ15N values) have been used to trace major sources and processes of atmospheric N (Heaton, 1986; Michalski et al., 2004; Kendall et al., 2007; Pavuluri et al., 2010; Savarino et al., 2013). The analysis of bulk δ15N in PM2.5 is a quick method compared to δ15N measurements of inorganic and organic N components (Widory, 2007; Hegde et al., 2015; Bikkina et al., 2016), and it does also provide valuable information on δ15N of dry N deposition (Yeatman et al., 2001a; Heaton et al., 2004; Elliott et al., 2007, 2009). At a background site, the δ15N of PM2.5 allows us to examine the impacts of emissions from non-point sources’ agricultural N emissions on the N chemistry of regional atmosphere, whereas at an urban site, it can imprint the anthropogenic N emissions.

The bulk δ15N in atmospheric particulates is mainly determined by the δ15N of N precursors (Aggarwal et al., 2013; Hegde et al., 2015). Often, reported δ15N values of typical inorganic N sources (Table 1) can be used in studies on sources and fates of atmospheric N (Elliott et al., 2007, 2009; Kendall et al., 2007; Kawashima and Kurahashi, 2011; Michalski et al., 2014). For PM2.5, dust is a primary N source (Zhang, 2010; Huang et al., 2014). At background site, NO2 and/or NH3 from microbial N cycle, fertilization application and animal wastes are strongly 15N-depleted (Elliott et al., 2007; Li and Wang, 2008; Felix et al., 2014), while N emissions from biomass burning (Kawashima and Kurahashi, 2011; Divers et al., 2014) are typically 15N-enriched. At urban site, most N sources of PM2.5 are anthropogenic. The NH3 from animal wastes (including sewages; Heaton, 1986), coal combustion and vehicle exhausts (Felix et al., 2013), as well as NO2 from vehicle exhausts (Walters et al., 2015) showed negative δ15N values, but NO2 from coal combustion had exclusively positive δ15N values (Felix et al., 2012).

Besides, the bulk δ15N in atmospheric particulates is influenced by the isotopic fractionations during gas (g) ↔ particle (p) exchange processes. However, isotope effects between N precursors and the aerosol N remain unclear, especially in field conditions. Isotope effects are assumed to be more significant between NH3 and NH4+, and much smaller in the case of NO2 and aerosol N (Yeatman et al., 2001a; Kawashima and Kurahashi, 2011). This assumption was supported by small differences in mean δ15N values between roadside NO2 (5.7‰) and local aerosol N (6.8‰) (Ammann et al., 1999; Pearson et al., 2000). Although the kinetic isotope effect of NH3 → NH4+ reaction is small at the beginning, it becomes significant when NH3 ↔ NH4+ equilibrium attains and causes a preferential enrichment of 14N in NH3 and 15N in NH4+ of aerosols (Heaton et al., 1997; Fukuzaki and Hayasaka, 2009; Li et al., 2012). This explained generally higher δ15N values of NH4+ in aerosols than that in rain NH4+ and gaseous NH3 (Yeatman et al., 2001a, 2001b; Jia and Chen 2010; Felix et al., 2013). In a hypothetical model by Heaton et al. (1997), the δ15N of particulate NH4+ stabilized at values of 33‰ (an enrichment coefficient) higher than that of NH3 when NH3 ↔ NH4+ equilibrium was achieved at 25 °C. However, mechanisms for atmospheric NH3 ↔ NH4+ equilibrium in the field circumstances are poorly understood, which is particularly important for interpreting the δ15N variations of PM2.5 at locations dominated by NH4+-N.

This study measured bulk δ15N of PM2.5 at an urban site (Chinese Research Academy of Environmental Sciences (CRAES), Beijing, northern China) and a national atmospheric background monitoring station (Menyuan, Qinghai province, northwestern China). Based on bulk δ15N of PM2.5 and major N sources, a Bayesian isotope mixing model (SIAR, Stable Isotope Analysis in R) (Parnell and Jackson, 2008) was used to estimate the proportions of different source contributions to N in PM2.5 and to evaluate anthropogenic N emissions during the haze events in Beijing. As inorganic N in the atmosphere of both sites was dominated by NH4-N, we hypothesized that 15N enrichments in PM2.5 relative to dominant sources were mainly derived from the NH3 ↔ NH4+ equilibrium (assumed as 33‰) (Heaton et al., 1997; Li et al., 2012).

2.  

Materials and methods

2.1.  

Study sites

The sampling site in Beijing is located in the courtyard of CRAES (40°04′ N, 116°42′ E), at Lishuiqiao South of Beiyuan Road. Due to rapid urbanization and economic development, the vehicle exhausts and energy consumption are large in Beijing, resulting in deterioration of air quality. Atmospheric PM2.5 in Beijing was characterized by high contributions of secondary components from anthropogenic origins (Sun et al., 2006). Secondary inorganic ions (such as SO42-, NH4+ and NO3-) were the dominant contributors in PM2.5 of Beijing (Zhang et al., 2013). During the sampling period (January 2013), Beijing suffered from the worst PM2.5 pollutions in history (http://cleanairinitiative.org/portal/node/11599), registering the highest PM2.5 hourly concentration of 886 μg/m3 (http://www.nasa.gov/multimedia/imagegallery/imagefeature2425.html).

The background site is located on the Daban Mountain (37°36′ N, 101°15′ E) in Menyuan county, north-eastern of Qinghai province, which is one of the 14 National Background Stations established by the Chinese Ministry of Environmental Protection in 2012. It has a typical Plateau continental climate, with an altitude of 3295 m above sea level, lower than the average of the Tibetan Plateau (about 4000 m). The mean annual temperature and precipitation amount are 0.8 °C and 520 mm, respectively. Agricultural activity is not intensive locally, except in low-altitude areas far away from the Daban Mountain in Menyuan. The sampling period (6 September–15 October 2013) is within the harvesting period after an intensive fertilization and pronounced biomass burning. The mean hourly temperature was 6.5 °C (3–11 °C) during the study period. There is no fossil fuel emission locally, with limited road traffic on the national highway of G227.

2.2.  

Sample collection and chemical analyses

PM2.5 was collected using a pre-baked quartz filter (47 mm in diameter) and aerosol sampler (Leckel, MVS6, Germany) equipped with a size-segregating impactor. The operating air flow rate was 38.3 L/min. To collect sufficient PM2.5 sample for bulk δ15N analyses, sampling was conducted for every 47–71 h at Menyuan (n = 14) and for 23 h at Beijing (n = 14). Filter blanks were also collected following the same procedure. The PM2.5 mass on each filter was gravimetrically measured using microbalance (AWS-1, COMDE DERENDA, Germany, approved by European Standard) after being desiccated for at least 24 h under controlled temperature (20 ± 1 °C) and humidity (50 ± 5%). All filter samples were immediately stored at −20 °C prior to chemical analyses.

Concentrations of bulk N in PM2.5 (mainly including NH4+, NO3- and organic N) were measured using three punches (ca. 0.53 cm2 for each) of the filter in a vario MACRO cube (Elementar Analysensysteme GmbH, Germany) with an analytical precision of 0.02%. Based on N contents, bulk δ15N values of about 50 μg N in each PM2.5 sample were determined by a Thermo MAT 253 isotope ratio mass spectrometer (Thermo Scientific, Bremen, Germany) coupled with an elemental analyzer (Flash EA 2000). IAEA-N-1 (Ammonium sulphate; δ15N = 0.4‰), USGS25 (Ammonium sulphate, δ15N = −30.4‰) and IAEA-NO-3 (Potassium nitrate; δ15N = +4.7‰) were used as standards for the calibration of δ15N values. The average standard deviation for replicate analyses of an individual sample was ±0.1‰. The δ15N in PM2.5 was expressed in parts per thousand (per mille) by multiplying them by 1000:

δ15N=(Rsample/Rstandard)-1,

where R = 15N/14N for samples and standard (atmospheric N2).

The concentrations of NO3-, NH4+ and SO42- in PM2.5 were measured during the sampling period at both sites by an ambient ion monitor (AIM-IC system: Model URG 9000B, URG Corporation, USA). It draws air in through a PM2.5 sharp-cut cyclone at a volumetric flow controlled rate of 3 L/min to remove the larger particles from the air stream. The real-time instruments installed at both the stations have a detection limit of 0.05 μg/m3. Gases such as SO2, NH3 and HNO3 are stripped from the air stream by passing through a liquid parallel plate denuder with continuously replenished solvent flowing across the surface. Then, the PM2.5 air stream is constrained into a supersaturated steam condensation coil and cyclone assembly and grown hygroscopically for collection. Enlarged particles are dissolved in water solutions for anion chromatographic analysis every hour following 60 min of ambient sampling. Concentrations of NO2 were measured using a NO–NO2–NOx chemiluminescence analyzer (Model 42i, Thermo-Fisher Scientific). The instruments were operated and maintained properly to ensure data integrity. Scheduled quality control procedures included daily zero and span checks, weekly precision checks and data validations.

3.  

Results

The PM2.5 levels at Beijing varied from 43.0 to 433.6 μg/m3 (mean = 264.3 ± 118.0 μg/m3) (Tables 1 and S2). Volumetric concentrations of elements and ions in PM2.5 differed distinctly between the two study sites, thus N contents were presented in the unit of N mass in PM2.5 mass for comparison. The bulk N and δ15N values of PM2.5 at Beijing averaged 16.7 ± 4.6% (8.2% to 29.3%) and +2.8 ± 6.4‰ (−4.1‰ to +13.5‰), respectively (Tables 1 and S2; Fig. 1). The NH4+-N, NO3--N and SO42--S in PM2.5 at Beijing averaged 7.4 ± 3.4%, 5.0 ± 3.0% and 5.5 ± 2.4%, respectively. The mean molar ratio of NH4+ to (NO3-+2SO42-) was 0.8 (Table 2). Ambient concentrations of NO2 (this study), NH3 (during April of 2013) and SO2 (during January of 2013) averaged 89.2 ± 21.2 μg/m3, 14.1 μg/m3 and 22.9 μg/m3 (He et al., 2014; Wei et al., 2015), respectively, showing a mean molar ratio of ambient NH3 to (NO2 + 2 * SO2) of 0.3 (Table 2).

Fig. 1.  

δ15N values of bulk N in PM2.5 and dominant N sources assigned for PM2.5 at the Beijing CRAES site (in red) and the Menyuan site (in blue). The box encompasses the 25th–75th percentiles, whiskers are SD values. The line and square in each box mark the median and arithmetic mean values, respectively. The number of jittered replicate δ15N data (dots around the boxes) is 1–34. Mean and SD values of source δ15N data were used in the SIAR model. δ15N values of N from dust were assumed as those of surface soils (Wang et al., 2014) according to the air mass backward trajectories (Fig. 2).

The filter-based average concentrations of PM2.5 at the background site (Menyuan, Qinghai province) varied from 7.0 to 17.8 μg/m3 (mean = 13.0 ± 3.2 μg/m3) (Tables 1 and S2), and (13.0 ± 4.8 μg/m3; 4.6–22.7 μg/m3) measured using an ambient monitor (AIM-IC system: Model URG 9000B, URG Corporation, USA). The bulk N concentrations and δ15N values of PM2.5 at the background site were 8.6 ± 5.6% and +18.5 ± 5.8‰ (+8.0‰ to +27.9‰), respectively (Tables 1 and S1; Fig. 1). Concentrations of NH4+-N, NO3--N and SO42--S in PM2.5 at Menyuan averaged 5.9 ± 1.8%, 1.9 ± 0.4% and 0.2 ± 0.0%, respectively (Table 2), showing a mean molar ratio of NH4+ to (NO3-+2SO42-) as 2.9 ± 1.0 (Table 2). Ambient concentrations of NO2 averaged 4.3 ± 1.3 μg/m3 at the background site (Table 2). Ambient NH3 and SO2 concentrations were not available at the Menyuan site (37°36′ N, 101°15′ E; 3295 m); however, these concentrations were reported as 4.8 μg/m3 and 0.31 μg/m3, respectively, at Waliguan (a global baseline station, 36°30′ N, 100°10′E, 3816 m), another background site in Qinghai (Carmichael et al., 2003). The estimated molar ratio of ambient NH3 to (NO2 + 2 * SO2) averaged 2.7 at the background site (Table 2).

4.  

Discussion

4.1.  

Major sources of N in PM2.5 of Beijing

According to the source appointment of PM2.5 at Beijing during the severe haze episode of January 2013 (Huang et al., 2014; Zhang et al., 2015), the following six dominant sources can be assigned for bulk N of PM2.5.

  • S1: N from dust,
  • S2: NO2 from coal combustion,
  • S3: NH3 from coal combustion,
  • S4: NO2 from vehicle exhausts,
  • S5: NH3 from vehicle exhausts,
  • S6: NH3 from animal wastes (mainly domestic wastes and sewages).

It should be explained that NO is the initial precursor for NOx emission sources, but NO is quite reactive and readily oxidized to NO2 which is more often taken as the precursor of NO3- in the atmosphere. Thus, NO2 was used in this work uniformly and its δ15N values were assumed as those of corresponding NOx emissions.

In this study, agricultural and biogenic N emissions were not considered as the major sources of bulk N in PM2.5 of Beijing for two main reasons. First, the urban site is located in the centre of Beijing city cluster in CRAES. During the severe haze events occurring in Beijing, several studies have shown that aerosols have been mainly influenced by anthropogenic sources. Second, as the sampling of Beijing PM2.5 was conducted in the winter time, contributions of NO2 from microbial N cycle, NH3 emission from seawater (δ15N = −8‰ to −5‰ in Jickells et al., 2003) and lightening NOx15N = −0.5‰ to +1.4‰; Hoering, 1957) were quite small, with relatively lower contribution than anthropogenic N sources to the formation of near-surface PM2.5, especially in urban circumstances.

To date, δ15N values of various NO2 and NH3 emissions are unavailable in China. However, according to source δ15N data compiled from previous studies (Table 1, Fig. 1), δ15N values were distinctive between most typical sources, which have been broadly used in isotopic tracing or partitioning of atmospheric N deposition (e.g. Elliott et al., 2007, 2009; Kawashima and Kurahashi, 2011). In this study, we did not use δ15N data of emissions influenced by post-emission processes and measured through controlled tests or simulation, e.g. the δ15N values of NH3 near highway (−5.0‰ to +0.4‰ in Smirnoff et al., 2012), NO2 near highway (+2‰ to +10‰ in Moore, 1977; Ammann et al., 1999; Pearson et al., 2000; −13.3‰ to +0.4‰ in Smirnoff et al., 2012), NO2 in tunnels (+15.0 ± 1.6‰ for NO2; +5.7 ± 2.8‰ for HNO3; Felix et al., 2014), NO2 from vehicle engine (−13.0‰ to +3.7‰; Moore, 1977; Freyer, 1978, 1991; Heaton, 1990), NO2 from controlled experiments of diesel combustion (+3.9‰ to +5.4‰; Widory, 2007) and coal combustion (−5.3‰; Widory, 2007). According to the air mass backward trajectories (Fig. 2), the δ15N values of surface soils in northern China (+4.3 ± 1.8‰; Wang et al., 2014) were used as the value of N from dust in this study.

Fig. 2.  

Seventy-two-h air mass backward trajectories for all sampling dates at the Beijing CRAES site and the Menyuan site, based on NOAA HYSPLIT model back trajectories.

As bulk δ15N values of PM2.5 at Beijing were distributed within those of major sources (Fig. 1), no substantial isotopic effect between N sources and bulk N of PM2.5 at Beijing was assumed. In particular, as inorganic N of PM2.5 was dominated by NH4+ (with a mean molar ratio of NH4+ to NO3- of 2.5; Table 2), the isotope effect of NH3 ↔ NH4+ equilibrium is considered quite low in the PM2.5 of Beijing. First, the low molar ratios of ambient NH3 to (NO2 + 2 * SO2) as 0.3 (Table 2) reflected a relatively thorough neutralization of NH3 by acidic gases, producing relatively more stable ammonium salts of NH4NO3, NH4HSO4 and (NH4)2SO4. Second, the molar ratios of NH4+ to (NO3-+2SO42-) were calculated as 0.8 (Table 2), indicating a full fixation of NH4+ by existing NO3- and SO42- for PM2.5 of Beijing. In the calculation, NH4+ is the actual molar concentrations of NH4+ in PM2.5 while the (NO3-+2SO42-) in PM2.5 represents the concentrations of NH4+ that can be fixed by NO3- and SO42-. More often, due to the high emissions of anthropogenic SO2 and NO2 in urban environments, NH3, after converting to NH4+, reacts mainly with acids formed by SO2 and NO2, with little opportunity of NH3 losses from PM2.5; thus, no substantial 15N enrichment in NH4+ of PM2.5 (Yeatman et al., 2001a; Pavuluri et al., 2010; Kawashima and Kurahashi, 2011) is observed. Consequently, bulk δ15N values of PM2.5 at Beijing were mainly controlled by the mixing of N sources with inappreciable isotopic effects.

4.2.  

Major sources of N in PM2.5 of Menyuan

According to the molar ratios of ambient NH3 to NO2 (ca. 3.0) or NH4+ to NO3- (ca. 3.3) in PM2.5 at Menyuan (Table 2), inorganic N in both ambient atmosphere and PM2.5 were dominated by NH3 and NH4+, respectively. Moreover, δ15N values of PM2.5 did not assemble those of dust N and/or natural N (mainly NO2 from N cycle) emissions; instead, they were much higher than those of potential sources (Table 1, Fig. 1). More likely, agricultural and biogenic NH3 sources should be important to bulk N of the background PM2.5. Hence, we assigned major N sources of PM2.5 at the background site as follows:

  • S7: N from dust,
  • S8: NO2 from biomass burning,
  • S9: NH3 from biomass burning,
  • S10: NH3 from animal wastes,
  • S11: NH3 from fertilizer application,
  • S12: NO2 from microbial N cycle.

The stoichiometry between ambient NH3 and acidic gases (NO2 + 2 * SO2), NH4+ and (NO3-+2SO42-) in PM2.5 allowed us to further interpret different patterns of bulk δ15N values of PM2.5 against those major sources (Figs. 1 and 2). At Menyuan, the molar ratios of ambient NH3 to (NO2 + 2 * SO2) averaged 2.7 (Table 2), illustrating an incomplete neutralization of ambient NH3 by NO2 and SO2. Molar ratios of NH4+ to (NO3-+2SO42-) in PM2.5 (ca. 2.9; Table 2) also suggested that part of NH4+ existed as relatively less stable ammonium salts (e.g. NH4Cl). The diffusion of NH3 back to the atmosphere during the reversible reaction and strong equilibrium between NH3 and NH4+ caused significant 15N enrichment in NH4+ of PM2.5, thus exhibiting much higher δ15N values of PM2.5 than potential sources (Fig. 1). The regulation of acidic gases-to-NH3 stoichiometry on the reaction and isotopic effect between NH3 and NH4+ was supported by a positive correlation between δ15N values and NH4+/(NO3-+2SO42-) ratios in PM2.5 (Fig. 3). Accordingly, a net isotopic effect of NH3 (g) ↔ NH4+ (p) at equilibrium (εeq) (33‰; Heaton et al., 1997) was considered in the SIAR model for the background PM2.5 (details down in Section 4.2). However, it should be noted that isotope effects for the atmospheric NH3(g) ↔ NH4+(p) equilibrium in the field circumstances remain unclear. The value of 33‰ is the only empirical one for 15N enrichment in particulate NH4+ (Heaton et al., 1997). Experimental studies have been conducted on the isotope fractionations of NH3 volatilization (e.g. Li et al., 2012), but it is uncertain what factors can be used to make corrections of the isotope effects for background PM2.5. Further studies are strongly needed to verify the relationships between the isotope effects and the ratio of NH3(g) to NH4+(p), which may be a feasible factor to make a correction of the isotope effects.

Fig. 3.  

Correlations between δ15N values of PM2.5 and molecular ratios of NH4+ to (NO3-+2SO42-) (expressed as n-NH4+/n-(NO3-+2NO3-+2SO42-)) in PM2.5 at Beijing CRAES site and Menyuan site. The regression line was drawn on data of both sites.

4.3.  

Using the SIAR model to partition bulk N in PM2.5

The proportional contributions (F, %) of major sources to N in PM2.5 are estimated using the SIAR model. This model uses a Bayesian framework to establish a logical prior distribution based on Dirichlet distribution (Evans et al., 2000), and then to determine the probability distribution for the contribution of each source to the mixture (Parnell and Jackson, 2008). It can substantially incorporate the uncertainties associated with multiple sources, fractionations and isotope signatures (Moore and Semmens, 2008; Davis et al., 2015). In our estimations, uncertainties should be evaluated for the δ15N variabilities of bulk N in PM2.5 and N sources, isotopic effect of the NH3 (g) ↔ NH4+(p) equilibrium.

By defining a set of N mixture measurements on J isotopes with K source contributors, the mixing model can be expressed as follows (Parnell et al., 2010):

Xij=k=1KFk(Sjk+cjk)+εijSijN(μjk,ωjk2)cijN(λjk,τjk2)εijN(0,σj2)

where all F values sum to 1 (unity), Xij is the isotope value j of the mixture i, in which i = 1, 2, 3, …, N and j = 1, 2, 3, …, J; Sjk is the source value k on isotope j (k = 1, 2, 3, …, K) and is normally distributed with mean μjk and standard deviation ωjk; Fk is the proportion of source k estimated by the SIAR model; cjk is the fractionation factor for isotope j on source k and is normally distributed with mean λjk and standard deviation τjk; and εij is the residual error representing the additional unquantified variation between individual mixtures and is normally distributed with mean 0 and standard deviation σj. A detailed description of this model can be found in Moore and Semmens (2008), Jackson et al. (2009) and Parnell et al. (2010). To estimate the contributions of N sources in the PM2.5 samples at two study sites (n = 14 for each), one isotope (j = 1) (δ15N of bulk N) and six potential N sources (as discussed in Sections 4.1 and 4.2: S1–S6 for Beijing and S7–S12 for Menyuan) (Fig. 1) are utilized. δ15N values of replicate PM2.5 samples at each study site were analysed in the SIAR model as one group.

Our estimation showed that the contribution of NO2 (FNO2) reached 41 ± 11% in bulk N of PM2.5 in Beijing, which was much higher than FNO2 at the background site (22 ± 10%) (Table 3). The mean ratios of FNH3 to FNO2 were about 1.6 and 4.4 for PM2.5 at Beijing and at the background site, respectively (Table 3), which generally followed the molar ratios of NH4+ to NO3- in PM2.5 (Table 2). Aqueous phase reaction experiments have shown that atmospheric NO2 and NH3 potentially react with organic compounds to form organic N (Ge et al., 2011; Pavuluri et al., 2015), which might contribute to the high secondary organic aerosols during the study haze event in Beijing (Huang et al., 2014).

In Beijing, anthropogenic N in PM2.5 averaged 81% of its bulk N and was mainly derived from N emissions of fossil fuel combustions, with the highest contribution (ca. 25%) from NO2 of coal combustion (Table 3; Fig. 4). The N emissions from coal combustion showed higher contributions (ca. 39%) than traffic emissions (ca. 32%), fossil-derived NO2 contributed more N (ca. 39%) than fossil-derived NH3 (ca. 30%) (Table 3; Fig. 4). Comparable contributions (ca. 14–16%) were observed between NH3 from coal combustion and NH3 from traffic emissions, between NH3 and NO2 from vehicle exhausts (Fig. 4). Accordingly, fossil-derived NH3 emissions substantially contributed to urban PM2.5 pollution; regulatory controls of N emissions from coal combustion and urban transportation are important to advert the risk of severe haze episodes in Beijing.

Fig. 4.  

Fractional contributions of dominant N sources to bulk N in PM2.5 at the Beijing CRAES site and the Menyuan site. Dots around the boxes (n = 104) show the percentages estimated by the SIAR model. The box encompasses the 25th–75th percentiles, whiskers are the 5th and 95th percentiles. The line and cross in each box mark the median and arithmetic mean values, respectively.

The N in PM2.5 at the background site was mainly contributed by N emissions from biomass burning (46 ± 10%) and NH3 volatilization (34 ± 12%) (Table 3). The contribution of NH3 from biomass burning (29 ± 6%) was comparable with the total contributions of NH3 from animal wastes and fertilizer application (ca. 35%) (Table 3). Biomass burning contributed less N as NO2 (17 ± 10%) than as NH3 (29 ± 6%) to N of PM2.5 at the background site (Table 3). Higher production of NH3 than NO2 from biomass burning has been documented previously (Crutzen and Andreae, 1990). A burning experiment by Lobert et al. (1990) showed higher emission ratios of NH3 (ca. 3.8%) than that of SO2 (ca. 0.3%) during biomass burning. The emission factors of NH3 were ca. 2–5 times higher than that of SO2 from various types of biomass burning (Andreae and Merlet, 2001).

5.  

Remarks

This study attempted to quantify major sources of N in PM2.5 based on bulk δ15N analysis using a Bayesian isotope mixing model. The isotopic effect of NH3 ↔ NH4+ equilibrium was recognized under the condition of lower acid gases relative to ambient NH3, which was a main reason for higher bulk δ15N of PM2.5 than potential sources at the background site. Based on the estimations of SIAR model, PM2.5 of Beijing derived N mainly from coal combustion and vehicle exhausts, while background PM2.5 derived N mainly from biomass burning and NH3 volatilization. Regulatory controls of N emissions from coal burning and urban transportation are important and effective steps to reduce the risk of severe haze episodes in Beijing. However, emissions of N from non-fossil emissions (particularly biomass burning) in broad rural areas should be stressed to meet a rigorous reduction of reactive N emissions in China.

Although δ15N interpretation using the SIAR model provided proportional contributions of major sources to bulk N in PM2.5, further investigations are needed to validate the assumptions and boundary conditions in this work. Particularly, δ15N analyses of gaseous N emissions should be implemented for reducing the uncertainties of source δ15N values. So far, isotopic studies on gaseous N emissions from typical anthropogenic and natural emissions are still sparse globally, especially in China. Isotope effects revealed in conversions between NO and NO2, NOx and NO3- (Monse et al., 1969; Walters and Michalski, 2015), NH3 and NH4+ (Heaton et al., 1997) and the regulatory mechanisms behind the kinetic and equilibrium isotope effects should be explored and properly considered in future studies.

Disclosure statement

No potential conflict of interest was reported by the authors.

Funding

This work was supported by the State Environmental Protection Commonweal Trade Scientific Research, Ministry of Environmental Protection of China [grant number 2013467010] and the National Natural Science Foundation of China [grant number 41273026], [grant number 41473081], [grant number 41522301], [grant number 41603007]. Xue-Yan Liu was also supported by the 11st Recruitment Program of Global Experts (the Thousand Talents Plan) for Young Professionals granted by the central budget of China, and Youth Innovation Promotion Association of Chinese Academy of Sciences [grant number 2015327].

Acknowledgements

All the financial support from fund and research support from the staff of CRAES are gratefully acknowledged.

References

  1. Aggarwal , S. G. , Kawamura , K. , Umarji , G. S. , Tachibana , E. , Patil , R. S. and co-authors . 2013 . Organic and inorganic markers and stable C-, N-isotopic compositions of tropical coastal aerosols from megacity Mumbai: sources of organic aerosols and atmospheric processing . Atmos. Chem. Phys. 13 , 4667 – 4680 . https://doi.org/10.5194/acp-13-4667-2013  

  2. Ammann , M. , Siegwolf , R. , Pichlmayer , F. , Suter , M. , Saurer , M. and co-authors . 1999 . Estimating the uptake of traffic-derived NO2 from 15N abundance in Norway spruce needles . Oecologia 118 , 124 – 131 . https://doi.org/10.1007/s004420050710  

  3. Andreae , M. O. and Merlet , P. 2001 . Emission of trace gases and aerosols from biomass burning . Glob. Biogeochem. Cy. 15 , 955 – 966 . https://doi.org/10.1029/2000GB001382  

  4. Berkemeier , T. , Ammann , M. , Mentel , T. F. , Pöschl , U. and Shiraiwa , M. 2016 . Organic nitrate contribution to new particle formation and growth in secondary organic aerosols from α-pinene ozonolysis . Environ. Sci. Technol. 50 , 6334 – 6342 .  

  5. Bikkina , S. , Kawamura , K. and Sarin , M. 2016 . Stable carbon and nitrogen isotopic composition of fine mode aerosols (PM2.5) over the Bay of Bengal: impact of continental sources . Tellus B. 68 , 31518 . https://doi.org/10.3402/tellusb.v68.31518  

  6. Carmichael , G. R. , Ferm , M. , Thongboonchoo , N. , Woo , J. H. , Chan , L. Y. and co-authors . 2003 . Measurements of sulfur dioxide, ozone and ammonia concentrations in Asia, Africa, and South America using passive samplers . Atmos. Environ. 37 , 1293 – 1308 . https://doi.org/10.1016/S1352-2310(02)01009-9  

  7. Crutzen , P. J. and Andreae , M. O. 1990 . Biomass burning in the tropics: impact on atmospheric chemistry and biogeochemical cycles . Science. 250 , 1669 – 1678 . https://doi.org/10.1126/science.250.4988.1669  

  8. Davis , P. , Syme , J. , Heikoop , J. , Fessenden-Rahn , J. , Perkins , G. and co-authors . 2015 . Quantifying uncertainty in stable isotope mixing models . J. Geophys. Res.: Biogeosci. 120 , 903 – 923 . https://doi.org/10.1002/2014JG002839  

  9. Divers , M. T. , Elliott , E. M. and Bain , D. J. 2014 . Quantification of nitrate sources to an urban stream using dual nitrate isotopes . Environ. Sci. Technol. 48 , 10580 – 10587 . https://doi.org/10.1021/es404880j  

  10. Elliott , E. M. , Kendall , C. , Boyer , E. W. , Burns , D. A. , Lear , G. and co-authors . 2009 . Dual nitrate isotopes in actively and passively collected dry deposition: utility for partitioning NOx sources contributing to landscape nitrogen deposition . J. Geophys. Res. Biogeosci. 114 , G04020 .  

  11. Elliott , E. M. , Kendall , C. , Wankel , S. D. , Burns , D. A. , Boyer , E. W. and co-authors . 2007 . Nitrogen isotopes as indicators of NOx source contributions to atmospheric nitrate deposition across the Midwestern and Northeastern United States . Environ. Sci. Technol. 41 , 7661 – 7667 .  

  12. Evans , J. S. B. T. , Handley , S. J. , Perham , N. , Over , D. E. and Thompson , V. A. 2000 . Frequency versus probability formats in statistical word problems . Cognition . 77 , 197 – 213 . https://doi.org/10.1016/S0010-0277(00)00098-6  

  13. Felix , J. D. and Elliott , E. M. 2014 . The isotopic composition of passively collected nitrogen dioxide emissions: vehicle, soil and livestock source signatures . Atmos. Environ. 92 , 359 – 366 . https://doi.org/10.1016/j.atmosenv.2014.04.005  

  14. Felix , J. D. , Elliott , E. M. , Gish , T. , Maghirang , R. , Cambal , L. and co-authors . 2014 . Examining the transport of ammonia emissions across landscapes using nitrogen isotope ratios . Atmos. Environ. 95 , 563 – 570 . https://doi.org/10.1016/j.atmosenv.2014.06.061  

  15. Felix , J. D. , Elliott , E. M. , Gish , T. , McConnell , L. and Shaw , S. 2013 . Characterizing the isotopic composition of atmospheric ammonia emission sources using passive samplers and a combined oxidation-bacterial denitrifier isotope ratio mass spectrometer method . Rapid. Commun. Mass. Sp. 27 , 2239 – 2246 . https://doi.org/10.1002/rcm.6679  

  16. Felix , J. D. , Elliott , E. M. and Shaw , S. L. 2012 . The isotopic composition of coal-fired power plant NOx: the influence of emission controls and implications for global emission inventories . Environ. Sci. Technol. 46 , 3528 – 3535 . https://doi.org/10.1021/es203355v  

  17. Freyer , H. D. 1978 . Preliminary 15N studies on atmospheric nitrogenous trace gases . Pure Appl. Geophys. 116 , 393 – 404 . https://doi.org/10.1007/BF01636894  

  18. Freyer , H. 1991 . Seasonal variation of 15N/14N ratios in atmospheric nitrate species . Tellus B. 43 , 30 – 44 . https://doi.org/10.3402/tellusb.v43i1.15244  

  19. Fukuzaki , N. and Hayasaka , H. 2009 . Seasonal variations of nitrogen isotopic ratios of ammonium and nitrate in precipitations collected in the Yahiko-Kakuda Mountains Area in Niigata Prefecture, Japan . Water Air. Soil Pollut. 203 , 391 – 397 . https://doi.org/10.1007/s11270-009-0026-8  

  20. Ge , X. L. , Wexler , A. S. and Clegg , S. L. 2011 . Atmospheric amines Part I. A review . Atmos. Environ. 45 , 524 – 546 . https://doi.org/10.1016/j.atmosenv.2010.10.012  

  21. Guo , S. , Hu , M. , Zamora , M. L. , Peng , J. , Shang , D. and co-authors . 2014 . Elucidating severe urban haze formation in China . Proc. Natl. Acad. Sci. USA . 111 , 17373 – 17378 . https://doi.org/10.1073/pnas.1419604111  

  22. Hastings , M. G. , Jarvis , J. C. and Steig , E. J. 2009 . Anthropogenic impacts on nitrogen isotopes of ice-core nitrate . Science . 324 , 1288 . https://doi.org/10.1126/science.1170510  

  23. Heaton , T. H. E. 1986 . Isotopic studies of nitrogen pollution in the hydrosphere and atmosphere: a review . Chem. Geol. 59 , 87 – 102 . https://doi.org/10.1016/0168-9622(86)90059-X  

  24. Heaton , T. H. E. 1987 . 15N/14N ratio of nitrate and ammonium in rain of Pretoria, south Africa . Atmos. Environ. 21 , 843 – 852 . https://doi.org/10.1016/0004-6981(87)90080-1  

  25. Heaton , T. H. E. 1990 . 15N/14N ratios of NOx from vehicle engines and coal-fired power stations . Tellus . 42 , 304 – 307 .  

  26. Heaton , T. H. E. , Spiro , B. and Robertson , S. M. C. 1997 . Potential canopy influences on the isotopic composition of nitrogen and sulphur in atmospheric deposition . Oecologia . 109 , 600 – 607 . https://doi.org/10.1007/s004420050122  

  27. Heaton , T. H. E. , Wynn , P. and Tye , A. M. 2004 . Low 15N/14N ratios for nitrate in snow in the High Arctic (79°N) . Atmos. Environ. 38 , 5611 – 5621 . https://doi.org/10.1016/j.atmosenv.2004.06.028  

  28. Hegde , P. , Kawamura , K. , Joshi , H. and Naja , M. 2015 . Organic and inorganic components of aerosols over the central Himalayas: winter and summer variations in stable carbon and nitrogen isotopic composition . Environ. Sci. Pollut. Res. 23 , 5997 – 6001 .  

  29. He , H. , Wang , Y. , Ma , Q. , Ma , J. , Chu , B. and co-authors., 2014 . Mineral dust and NOx promote the conversion of SO2 to sulfate in heavy pollution days . Sci. Rep. 4 , 4172 – 4176 .  

  30. Hoering , T. 1957 . The isotopic composition of ammonia and the nitrate ion in rain . Geochim. Cosmochim. Acta. 12 , 97 – 102 . https://doi.org/10.1016/0016-7037(57)90021-2  

  31. Huang , R. J. , Zhang , Y. L. , Bozzetti , C. , Ho , K. F. , Cao , J. J. and co-authors . 2014 . High secondary aerosol contribution to particulate pollution during haze events in China . Nature . 514 , 218 – 222 .  

  32. Jackson , A. L. , Inger , R. , Bearhop , S. and Parnell , A. 2009 . Erroneous behaviour of MixSIR, a recently published Bayesian isotope mixing model: a discussion of Moore & Semmens (2008) . Ecol. Lett. 12 , E1 – E5 . https://doi.org/10.1111/ele.2009.12.issue-3  

  33. Jia , G. and Chen , F. 2010 . Monthly variations in nitrogen isotopes of ammonium and nitrate in wet deposition at Guangzhou, south China . Atmos. Environ. 44 , 2309 – 2315 . https://doi.org/10.1016/j.atmosenv.2010.03.041  

  34. Jickells , T. D. , Kelly , S. D. , Baker , A. R. , Biswas , K. , Dennis , P. F. and co-authors . 2003 . Isotopic evidence for a marine ammonia source . Geo. Res. Lett. 30 , 359 – 376 .  

  35. Kawashima , H. and Kurahashi , T. 2011 . Inorganic ion and nitrogen isotopic compositions of atmospheric aerosols at Yurihonjo, Japan: implications for nitrogen sources . Atmos. Environ. 45 , 6309 – 6316 . https://doi.org/10.1016/j.atmosenv.2011.08.057  

  36. Kendall , C. , Elliott , E. M. and Wankel , S. D. 2007 . Tracing anthropogenic inputs of nitrogen to ecosystems . In: Stable Isotopes in Ecology and Environmental Science (eds. R.M. Michener and K.E. Lajtha ). Blackwell , Oxford , pp. 375 – 449 . https://doi.org/10.1002/9780470691854  

  37. Li , L. , Lollar , B. S. , Li , H. , Wortmann , U. G. and Lacrampe-Couloume , G. 2012 . Ammonium stability and nitrogen isotope fractionations for NH4+-NH3(aq)-NH3 (gas) systems at 20–70 °C and pH of 2–13: applications to habitability and nitrogen cycling in low-temperature hydrothermal systems . Geochim. Cosmochi. Acta . 84 , 280 – 296 . https://doi.org/10.1016/j.gca.2012.01.040  

  38. Li , D. J. and Wang , X. M. 2008 . Nitrogen isotopic signature of soil-released nitric oxide (NO) after fertilizer application . Atmos. Environ. 42 , 4747 – 4754 . https://doi.org/10.1016/j.atmosenv.2008.01.042  

  39. Lobert , J. M. , Scharffe , D. H. , Hao , W. M. and Crutzen , P. J. 1990 . Importance of biomass burning in the atmospheric budgets of nitrogen-containing gases . Nature . 346 , 552 – 554 . https://doi.org/10.1038/346552a0  

  40. Michalski , G. , Bhattacharya , S. and Girsch , G. 2014 . NOx cycle and the tropospheric ozone isotope anomaly: an experimental investigation . Atm. Chem. Phys. 14 , 4935 – 4953 . https://doi.org/10.5194/acp-14-4935-2014  

  41. Michalski , G. , Meixner , T. , Fenn , M. , Hernandez , L. , Sirulnik , A. , Allen , E. and Thiemens , M. 2004 . Tracing atmospheric nitrate deposition in a complex semiarid ecosystem using Δ17O . Environ. Sci. Technol. 38 , 2175 – 2181 . https://doi.org/10.1021/es034980+  

  42. Monse , E. U. , Spindel , W. and Stern , M. J. 1969 . Analysis of isotope-effect calculations illustrated with exchange equilibria among oxynitrogen compounds, ACS Adv . Chem. Ser. 89 , 148 – 184 .  

  43. Moore , H. 1977 . The isotopic composition of ammonia, nitrogen dioxide, and nitrate in the atmosphere . Atmos. Environ. 11 , 1239 – 1243 . https://doi.org/10.1016/0004-6981(77)90102-0  

  44. Moore , J. W. and Semmens , B. X. 2008 . Incorporating uncertainty and prior information into stable isotope mixing models . Ecol. Lett. 11 , 470 – 480 . https://doi.org/10.1111/j.1461-0248.2008.01163.x  

  45. Parnell , A. C. , Inger , R. , Bearhop , S. and Jackson , A. L. 2010 . Source partioning using stable isotopes: coping with too much variation . PLoS ONE . 5 , e9672 . https://doi.org/10.1371/journal.pone.0009672  

  46. Parnell , A. C. and Jackson , A. 2008 . SIAR: stable isotope analysis in R . Online at: http://cran.r-project.org/web/packages/siar/index.html (accessed 10 December 2012).  

  47. Pavuluri , C. M. , Kawamura , K. and Fu , P. Q. 2015 . Atmospheric chemistry of nitrogenous aerosols in northeastern Asia: biological sources and secondary formation . Atmos. Chem. Phys. 15 , 9883 – 9896 . https://doi.org/10.5194/acp-15-9883-2015  

  48. Pavuluri , C. M. , Kawamura , K. , Tachibana , E. and Swaminathan , T. 2010 . Elevated nitrogen isotope ratios of tropical Indian aerosols from Chennai: implication for the origins of aerosol nitrogen in South and Southeast Asia . Atmos. Environ. 44 , 3597 – 3604 . https://doi.org/10.1016/j.atmosenv.2010.05.039  

  49. Pearson , J. , Wells , D. , Seller , K. J. , Bennett , A. , Soares , A. and co-authors . 2000 . Traffic exposure increases natural 15N and heavy metal concentrations in mosses . New Phytol. 147 , 317 – 326 . https://doi.org/10.1046/j.1469-8137.2000.00702.x  

  50. Savarino , J. , Morin , S. , Erbland , J. , Grannec , F. , Patey , M. D. and co-authors . 2013 . Isotopic composition of atmospheric nitrate in a tropical marine boundary layer . Proc. Natl. Acad. Sci. USA . 110 , 17668 – 17673 . https://doi.org/10.1073/pnas.1216639110  

  51. Smirnoff , A. , Savard , M. M. , Vet , R. and Simard , M. C. 2012 . Nitrogen and triple oxygen isotopes in near-road air samples using chemical conversion and thermal decomposition . Rapid. Commun. Mass. SP. 26 , 2791 – 2804 . https://doi.org/10.1002/rcm.6406  

  52. Sun , Y. L. , Zhuang , G. S. , Tang , A. H. , Wang , Y. and An , Z. S. 2006 . Chemical characteristics of PM2.5 and PM10 in haze-fog episodes in Beijing . Environ. Sci. Technol. 40 , 3148 – 3155 . https://doi.org/10.1021/es051533g  

  53. Walters , W. W. , Goodwin , S. R. and Michalski , G. 2015 . Nitrogen stable isotope composition (15N) of vehicle emitted NOx . Environ. Sci. Technol. 49 , 2278 – 2285 . https://doi.org/10.1021/es505580v  

  54. Walters , W. W. and Michalski , G. 2015 . Theoretical calculation of nitrogen isotope equilibrium exchange fractionation factors for various NOy molecules . Geochim. Cosmochim. Acta 164 , 284 – 297 . https://doi.org/10.1016/j.gca.2015.05.029  

  55. Wang , C. , Wang , X. B. , Liu , D. W. , Wu , H. H. , Lü , X. T. and co-authors . 2014 . Aridity threshold in controlling ecosystem nitrogen cycling in arid and semi-arid grasslands . Nat. Commun. 5 , 4799 . https://doi.org/10.1038/ncomms5799  

  56. Wei , L. F. , Duan , J. C. , Tan , J. H. , Yongliang , M. A. , Kebin , H. E. and co-authors . 2015 . Gas-to-particle conversion of atmospheric ammonia and sampling artifacts of ammonium in spring of Beijing . Sci. China. Earth. Sci. 63 , 186 – 187 .  

  57. Widory , D. 2007 . Nitrogen isotopes: tracers of origin and processes affecting PM10 in the atmosphere of Paris . Atmos. Environ. 41 , 2382 – 2390 . https://doi.org/10.1016/j.atmosenv.2006.11.009  

  58. Yeatman , S. G. , Spokes , L. J. , Dennis , P. F. and Jickells , T. D. 2001a . Comparisons of aerosol nitrogen isotopic composition at two polluted coastal sites . Atmos. Environ. 35 , 1307 – 1320 . https://doi.org/10.1016/S1352-2310(00)00408-8  

  59. Yeatman , S. G. , Spokes , L. J. and Jickells , T. D. 2001b . Comparisons of coarse-mode aerosol nitrate and ammonium at two polluted coastal sites . Atmos. Environ. 35 , 1321 – 1335 . https://doi.org/10.1016/S1352-2310(00)00452-0  

  60. Zhang , R. Y. 2010 . Getting to the critical nucleus of aerosol formation . Science . 328 , 1366 – 1367 . https://doi.org/10.1126/science.1189732  

  61. Zhang , R. J. , Jing , J. , Tao , J. , Hsu , S.-C. , Wang , G. and co-authors . 2013 . Chemical characterization and source apportionment of PM2.5 in Beijing: seasonal perspective . Atmos. Chem. Phys. 13 , 7053 – 7074 . https://doi.org/10.5194/acp-13-7053-2013  

  62. Zhang , L. , Wang , T. , Lv , M. Y. and Zhang , Q. 2015 . On the severe haze in Beijing during January 2013: unraveling the effects of meteorological anomalies with WRF-Chem . Atmos. Environ. 104 , 11 – 21 . https://doi.org/10.1016/j.atmosenv.2015.01.001  

comments powered by Disqus