Stable carbon and nitrogen isotopic composition of fine mode aerosols (PM2.5) over the Bay of Bengal: impact of continental sources

Stable carbon and nitrogen isotopic composition of fine mode aerosols (PM2.5) over the Bay of Bengal: impact of continental sources

By SRINIVAS BIKKINA1,2, KIMITAKA KAWAMURA1* and MANMOHAN SARIN2, 1Institute of Low Temperature Science, Hokkaido University, Japan; 2Geosciences Division, Physical Research Laboratory, Ahmedabad, India


This study reports on stable carbon (δ13CTC) and nitrogen (δ15NTN) isotopic composition of total carbon and nitrogen (TC and TN) in the fine mode aerosols (PM2.5; N=31) collected over the Bay of Bengal (BoB). The samples represent two distinct wind regimes during the cruise (27 December 2008–28 January 2009); one from the Indo-Gangetic Plain (referred as IGP-outflow) and another from Southeast Asia (SEA-outflow). The PM2.5 samples from the IGP-outflow show higher δ13CTC (−25.0 to −22.8 ‰; −23.8±0.6 ‰) than those from the SEA-outflow (−27.4 to −24.7 ‰; −25.3±0.9 ‰). Similarly, δ15NTN varied from +11.8 to +30.6 ‰ (+20.4±5.4 ‰) and +10.4 to +31.7 ‰ (+19.4±6.1 ‰) for IGP- and SEA-outflows, respectively. Based on the literature data, MODIS-derived fire hotspots and back trajectories, we infer that higher δ13CTC in the IGP-outflow is predominantly associated with fossil fuel and biofuel combustion. In contrast, contribution of primary organic aerosols from the combustion of C3 plants or secondary organic aerosol (SOA) formation from biomass/biofuel-burning emissions (BBEs) can explain the lower δ13CTC values in the SEA-outflow. This inference is based on the significant linear correlations among δ13CTC, water-soluble organic carbon and non-sea-salt potassium (nss-K+, a proxy for BBEs) in the SEA-outflow. A significant linear relationship of δ15N with and equivalent mass ratio of / is evident in both the continental outflows. Since abundance dominates the TN over the BoB (>90 %), atmospheric processes affecting its concentration in fine mode aerosols can explain the observed large variability of δ15NTN.

Keywords: stable C- and N-isotopes, Bay of Bengal, aerosols, biomass burning, fossil-fuel combustion, South Asia, South-east Asia, Indo-Gangetic Plain

Citation: Tellus B 2016, 68, 31518,

Responsible Editor: Kaarle Hämeri, Helsinki University, Finland.

Copyright: © 2016 S. Bikkina et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution 4.0 International License, allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material for any purpose, even commercially, provided the original work is properly cited and states its license.

Received: 5 March 2016; Accepted: 3 June 2016; Published: 12 July 2016

*Correspondence to: email:

To access the supplementary material to this article, please see Supplementary files under ‘Article Tools’.

1. Introduction

The rapidly growing human activities and significant increase in emissions of airborne pollutants from South Asia are of considerable interest owing to their impact on regional climate (Ramanathan et al., 2001a; Lawrence and Lelieveld, 2010). In this context, organic aerosols are considered to be one of the important components of atmospheric particulate matter for their crucial role in regional radiative forcing through direct and indirect effects (Kanakidou et al., 2005; Jimenez et al., 2009). Several studies have attempted to characterize the sources, transport and transformation pathways of organic aerosols over South Asia (Jayaraman et al., 1998; Novakov et al., 2000; Lelieveld et al., 2001; Ramanathan et al., 2001b; Mayol-Bracero et al., 2002a; Rengarajan et al., 2007; Lawrence and Lelieveld, 2010; Ram et al., 2010; Rajput et al., 2011; Srinivas and Sarin, 2013a, 2014). Jimenez et al. (2009) demonstrated that the identified fraction of atmospheric organic aerosols is no more than 20–30 % of total organic matter. These observations emphasize our ability to assess the impact of organic aerosols on climate forcing, mainly due to limited observations. Therefore, the chemical characterization of various organics and their sources from different geographical locations could provide a better means to reassess their realistic aerosol radiative forcing (Jacobson, 2012).

Recent studies on stable isotopic composition of bulk carbon and nitrogen (δ13C and δ15N) of airborne particulatematter have demonstrated their usefulness to apportion the sources and formation pathways of organic aerosols (Kawamura et al., 2004; Aggarwal and Kawamura, 2008; Narukawa et al., 2008; Pavuluri et al., 2011; Aggarwal et al., 2013; Kirillova et al., 2013). Previous studies have used δ13C and δ15N of bulk particulate matter to assess the contribution of biomass/biofuel-burning sources (Martinelli et al., 2002; Kundu et al., 2006; Ometto et al., 2006; Widory, 2007; Aggarwal and Kawamura, 2008; Agnihotri et al., 2011; Aggarwal et al., 2013; Mkoma et al., 2014). Although the stable isotopic composition of OA produced from various sources are significantly different, several processes can affect their δ13C and δ15N, viz., ageing of air masses, condensation/evaporation of volatile organic compounds (VOCs) on/from a pre-existing particles, oxidation or chemical reactions of primary organic aerosols (POA) by O3 or NOx (Kirillova et al., 2013 and references therein).

However, studies from South Asia are limited in assessing the stable C- and N-isotopic composition of ambient aerosols (Agnihotri et al., 2011; Pavuluri et al., 2011; Hegde and Kawamura, 2012; Aggarwal et al., 2013). Based on the physicochemical properties of aerosols, earlier studies have demonstrated the continental impact from South and Southeast Asia over the northern Indian Ocean (NIO) (Lelieveld et al., 2001; Ramanathan et al., 2001b; Mayol-Bracero et al., 2002a). Furthermore, the predominant continental influence on the eastern part of the NIO [the Bay of Bengal (BoB)] has been emphasized by subsequent studies compared to that over the Arabian Sea, western side of the NIO (Kumar et al., 2008; Kedia et al., 2010; Sarin et al., 2010; Srinivas et al., 2011a; Srinivas and Sarin, 2012, 2013a). These studies have further suggested that the strength of the continental sources decreases from winter (December–February) to spring–intermonsoon (March–April) (Srinivas et al., 2011b; Srinivas and Sarin, 2013b, 2014).

A recent study by Agnihotri et al. (2011) has assessed stable C- and N-isotopic composition of marine aerosols from the BoB and Arabian Sea during spring to intermonsoon (March–May 2006). However, no single study deals with the assessment of δ13C of bulk carbon and δ15N of total nitrogen in marine aerosols from the BoB during the wintertime when marine atmospheric boundary layer (MABL) is influenced by continental air masses from South and Southeast Asia. Here, we have studied the stable isotopic composition of total carbon and nitrogen (TC and TN) of fine mode aerosols (PM2.5) collected from the BoB during winter season. We have also made a comparison of our results from the Bay region with those from the nearby continental sites as well as other oceanic regions. The overall objective of this study is to understand the measured stable C- and N-isotopic composition of aerosol-TC and TN, respectively, in terms of their source contribution as well as changes in composition (if any) during atmospheric transport to the oceanic region (BoB) located downwind of pollution sources.

2. Methodology

2.1. Cruise track and prevailing meteorology

Aerosol samples (PM2.5) were collected during a cruise (SK-254) conducted in the BoB as part of Indian national programme on integrated campaign of aerosols, trace gases and radiation budget (ICARB) during 27 December 2008–28 January 2009. During the study period, prevailing winds were mostly north-easterly/westerly. The shallow boundary layer along with the relatively weak winds brings continental air masses to the MABL during the wintertime (December–February). We obtained the meteorological parameters on board (viz., relative humidity, ambient temperature, wind speed and wind direction) every hour, which were averaged for daily means. The daily average of relative humidity, ambient temperature and wind speed during the cruise varied from 49.1 to 79.8 %, 22.3 to 26.8 °C and 1.2 to 6.3 m s−1, respectively (Kumar et al., 2010). Further details on wind regimes and meteorological conditions are described in our earlier publications (Kumar et al., 2010; Srinivas et al., 2011b).

2.2. Aerosol sampling and measurement of chemical constituents

In this study, we used a total number of 31 PM2.5 samples that were collected using a high volume air sampler (HVS, Thermo-Anderson Tech., flow rate: 1.13 m3 minute−1) on board ORV Sagar Kanya. The sampler was placed on the upper deck in front of the ship’s navigation room in order to avoid the contamination from the ship’s smoke stack and was operated when the wind direction is from the bow and ship cruising at a speed of more than 10 knots h−1. The sampler was calibrated before and after the cruise, and the variation in the flow rate was within 3 %. On average, the samples were collected for a period of ~20 h, and the volume of filtered air ranged between 1300 and 1400 m3. All PM2.5 samples were collected on pre-combusted (at 450 °C, for 3–4 h) quartz filters (PALLFLEX®™, 2500 QAT-UP). After collection, aerosol filters were wrapped in aluminium foils, sealed in zip-lock bags and stored in deep freezer at −19 °C until the chemical analyses.

In PM2.5, water-soluble inorganic ions (Na+, , K+, Mg2+, Ca2+, Cl, and ) were determined withDionex-500 ion chromatograph equipped with suppressed conductivity detector (Srinivas et al., 2011b). Similarly, carbonaceous components (EC/OC) were measured on Sunset EC/OC analyser using NIOSH-5040 method (Srinivas and Sarin, 2013a), whereas water-soluble organic carbon (WSOC) was measured using Schimadzu TOC-5000 analyser (Srinivas and Sarin, 2013a). A part of the chemical composition data (K+, , , , EC, OC and WSOC) from earlier publications have been used here to support the inferences related to stable C- and N-isotopic composition of PM2.5 over the BoB. For more details, see our earlier publications (Rengarajan et al., 2007; Kumar et al., 2010; Srinivas and Sarin, 2014).

To determine the blank levels, one quarter of the filter aliquot (ca. ~100 cm2) was extracted with 50 ml of Milli-Q water (specific resistivity ~18.2 MΩ-cm) like samples and the extracts were measured for water-soluble inorganic constituents. Except for Na+, K+ and Cl, no signals were detected on ion chromatogram. The detection limits, as defined by three times of standard deviation of the concentrations of procedural filter blanks (N=8) normalized to average filtered volume of air for the PM2.5 (~1400 m3), were 30 ng m−3 for Na+, 18 ng m−3 for K+ and 24 ng m−3 for Cl. We have corrected the K+ and concentrations in PM2.5 for the sea-salt contribution to derive the non-sea-salt (nss) K+ and , respectively, using measured Na+ in aerosols (Keene et al., 1986). All PM2.5 samples were analysed to assess the mass concentrations of total carbon and nitrogen (TC and TN) along with their stable isotopic composition (δ13CTC and δ15NTN) using an elemental analyser (EA) and EA/isotope ratio mass spectrometer as reported in Kawamura et al. (2004) and Kundu et al. (2010).

2.3. Stable C- and N-isotopic composition

For this study, stable isotopic composition of δ13CTC and δ15NTN were determined using an EA interfaced online to an isotope ratio mass spectrometer (EA/irMS, model: Carlo Erba NA 1500 EA+Finnigan MAT Delta Plus), and the protocol described in Kawamura et al. (2004). Briefly, the filter aliquot (ca. 2.0 cm in diameter) of aerosol sample was sealed in a tin cup, which is pre-cleaned with acetone to remove organic contaminants and dried prior to use, and introduced into EA where it is oxidized in the combustion column packed with copper oxide (CuO) filings at a temperature of 1020 °C. The evolved gases (CO2 and NOx) were subsequently introduced into a reduction column, where NOx is converted to molecular N2 and thereby separated by a gas chromatograph equipped within EA. To measure stable isotopic composition of evolved CO2 and N2, an interface (ConFlo II) is used to transfer the gases from EA to irMS. The measured isotopic composition of aerosol-TC and TN was expressed as δ13CTC and δ15NTN relative to Pee Dee Beleminite carbon standard and atmospheric N2, respectively, using the following equations:



2.4. Quality assurance

Prior to the measurements of δ13CTC and δ15NTN, instrumental blanks are checked with pre-cleaned tin cups (N=3) and a five-point calibration, having standard amount between 0.2 and 0.6 mg of acetanilide (Thermo Scientific) with known stable carbon and nitrogen isotopic composition as −27.26 ‰ (Cao et al., 2016) and +11.88 ‰ (Kundu et al., 2010), respectively. Filter (lab+field) blanks were also analysed and the analytical signal of samples was suitably corrected for blanks using isotopic mass balance equations described in Turekian et al. (2003). The overall analytical uncertainties in the measurement of δ13CTC and δ15NTN were 0.15 ‰ and 0.25 ‰, respectively, as ascertained by repeat analyses of acetanilide and samples, being consistent with those reported in Kawamura et al. (2004). The uncertainties associated with each measurement were propagated and discussed when comparing the significant differences between the Indo-Gangetic Plain (IGP)- and Southeast Asia (SEA)-outflows (Supplementary Table 1).

In this study, TC concentrations (OC+EC) were measured in PM2.5 soon after their collection on Sunset EC–OC analyser using NIOSH protocol (Rengarajan et al., 2007). These measurements show consistency with total carbon content obtained from the EA-irMS. This comparison clearly demonstrates that artefacts, if any, related to shelf storage of aerosol filters have not caused any error in the isotopic measurements. Similarly, the closure of mass concentration of water-soluble inorganic nitrogen (mainly as -N; >90 %) in the TN (Fig. 2) further establishing the consistency between previous measurements (Srinivas et al., 2011b) and this study. This observation clearly demonstrates the lack of artefacts related to the storage of aerosol filters.

3. Results and discussion

3.1. Temporal variability

During a winter cruise, the PM2.5 collected over the N-BoB is influenced by the continental outflow from the IGP-outflow (27 December 2008 to 10 January 2009, Fig. 1), whereas those sampled over the S-BoB has influence from SEA-outflow (11–28 January 2009; Fig. 1). This is inferred based on the origin of 7-d isentropic air mass back trajectories (AMBTs) for the sampling days, computed using hybrid single particle Lagrangian-integrated trajectory model, HYSPLIT-4 (Draxler et al., 1999). The impact of IGP- and SEA-outflows on the chemical composition of aerosols over the BoB is a conspicuous seasonal meteorological phenomenon, occurring in winter and spring, documented through large-scale field campaigns (INDOEX and ICARB-I & II) using various optical and chemical properties of aerosols (Kumar et al., 2010; Lawrence and Lelieveld, 2010; Srinivas et al., 2011b).

Fig 1
Fig. 1.   Seven-day isentropic air mass back trajectories (AMBTs) computed using hybrid single particle Lagrangian-integrated trajectory model, HYSPLIT-4 (Draxler et al., 1999) for the sampling days (27 December 2008 to 26 January 2009) along the cruise track in the Bay of Bengal. Blue pins and red dots over the Indo-Gangetic Plain represent the locations of coal-fired power plants and MODIS-derived fire spots, respectively.

Organic aerosols in winter over the IGP are mostly emitted from the fossil fuel combustion sources (FF-comb) and biomass/biofuel-burning emissions (BBEs) (wood/crop-residue combustion) (Sudheer and Sarin, 2008; Rajput et al., 2011; Srinivas and Sarin, 2014; Rastogi et al., 2015). However, with the exception of a few studies emphasizing the contribution of organic aerosols from the forest fires over South-east Asia (Narukawa et al., 1999; Streets et al., 2003; Srinivas and Sarin, 2013a), no such detailed information exists in the literature about other pollution sources. To understand the influence of different pollution sources to atmospheric organic aerosols sampled over the BoB (this study), we have superimposed the locations of coal-fired thermal power plants in the IGP (see blue pins in Fig. 1) and MODIS-derived fire count data (red dots in Fig. 1) from South and Southeast Asia on the back trajectory cluster information (Fig. 1).

From Fig. 1, it is implicit that the SEA-outflow has a contribution from intense forest fires occurring over Indonesia, Thailand, Vietnam and Myanmar. However, no such influence from BBEs was observed in the IGP-outflow during winter. On the other hand, the IGP-outflow has a contribution of organic aerosols from coal-fired thermal power plants in the IGP [also referred as ‘coal belt of India’ (Prasad et al., 2006; Nair et al., 2007)]. Besides this, the residential biofuel (wood) combustion emissions (particularly significant in winter) also contribute to organic aerosols in the IGP-outflow (Gustafsson et al., 2009; Ram and Sarin, 2012). Therefore, we try to infer the relative significance of various sources and compositional changes in the continental outflows (consisting of emissions from coal-fired power plants, vehicular traffic and biomass burning) by comparing the observed δ13CTC over the BoB with source-specific signatures and regression analysis between isotopic composition and other chemical tracers measured in the MABL.

TC and TN in PM2.5 exhibit strong latitudinal variability with relatively high concentrations in the IGP-outflow (close to the continent) compared to those in the SEA-outflow (open ocean samples; Table 1). Similar north (high in the IGP-outflow) to south (low in the SEA-outflow) concentration gradient is also noteworthy for other chemical species in the PM2.5 (e.g. OC, EC, WSOC and ; Supplementary Fig. 1). The δ13CTC has shown pronounced spatial variability over the BoB (Fig. 3) with significantly higher δ13CTC (t-score: 5.5, df=29; p<0.001) in the IGP-outflow (−25.0 to −22.8 ‰; average −23.8±0.6 ‰) than SEA-outflow (−26.9 to −24.2 ‰; average −25.3±0.9 ‰). This discrepancy is due to varying contributions of aerosol-TC from FF-comb and BBEs over the North and South BoB, respectively. However, no significant differences were observed for δ15NTN (t-score: 0.5, df=29, p>0.05) in both the outflows (IGP: +20.4±5.4 ‰ and SEA: +19.4±6.1 ‰; Fig. 3).

Table 1. Statistical description of mass concentrations (in µg m−3) of chemical constituents and stable carbon and nitrogen isotopic compositions (δ15N and δ13C) in PM2.5 collected over the Bay of Bengal during January 2009
Parameter IGP-outflow (N=15) SEA-outflow (N=16) t-Score, df, p-value
aPM2.5 13.2–76.7 (38.0±20.2) 2.0–35.3 (22.3±9.9) 2.8, 29, <0.05
TN 0.6–8.8 (3.7±2.5) 0.2–3.0 (1.4±0.8) 3.7, 29, <0.05
TC 3.6–15.4 (8.1±4.3) 1.1–6.5 (3.8±1.8) 3.4, 29, <0.05
TC/TN 1.3–9.4 (3.0±2.1) 1.9–5.4 (2.9±0.8) 0.1, 29, >0.05
δ13CTC (‰) −25.0 to −22.8 (−23.8±0.6) −26.9 to −24.2 (−25.3±0.9) 7.3, 29, >0.05
δ15NTN (‰) +11.8 to +30.6 (20.4±5.4) +10.4 to +31.7 (19.4±6.1) 0.5, 29, <0.05
aOC 1.6–11.6 (5.9±3.7) 0.2–5.3 (2.6±1.9) 3.0, 27, <0.05
aEC 0.8–5.0 (2.0±1.3) 0.2–1.8 (1.2±0.4) 2.4, 29, <0.05
WSOC 1.5–8.3 (4.4±2.4) 0.5–4.1 (2.3±1.4) 2.9, 29, <0.05
WIOC 0.0–4.5 (1.6±1.4) 0.0–1.8 (0.6±0.5) 2.6, 29, <0.05
a 0.4–10.1 (4.4±3.0) 0.3–4.1 (2.1±1.2) 2.7, 29, <0.05
anss- 1.6–28.5 (11.8±7.8) 1.3–12.1 (7.0 ±3.1) 2.3, 29, <0.05
(/nss-)equivalent 0.62–1.12 (0.92±0.15) 0.60–0.93 (0.78±0.11) 3.2, 29, <0.05
nss-K+/EC 0.19–0.51 (0.29±0.09) 0.23–0.52 (0.39±0.09) 3.0, 29, <0.05
nss-/EC 1.3–14.5 (6.2±3.6) 2.7–10.3 (6.2±2.2) 0.01, 29, >0.05
aSupporting parameters were obtained from Kumar et al. (2010).

In comparison to δ13CTC, δ15NTN show large variability in both outflows (IGP: +11.8 to +31.6 ‰; SEA: +10.4 to +31.7 ‰; Table 1). The wide spread combustion of biofuels and agricultural crop residue from South and Southeast Asia is a major source of aerosol-TN over the BoB during winter and spring seasons (Srinivas et al., 2011b). In this regard, Turekian et al. (1998) documented that unlike δ13CTC, δ15NTN in particles emitted from the combustion of C3 and C4 vegetation overlap within spread of the data. However, their study suggested that changes in combustion temperature during BBEs lead to large variability in the stable N-isotopic composition of aerosols. This is partly due to accessibility of various nitrogenous pools in the source plant at different temperatures during biomass combustion (Turekian et al., 1998). Therefore, observed large variability in the δ15NTN over the BoB is attributed to variations in the combustion temperature during BBEs in the continental outflows (i.e. wood burning in the IGP or forest fires in the SEA).

Although wood-burning emissions usually involved with low-temperature combustion processes than forest fires (i.e. high-temperature combustion), however, we observe no significant differences (t-score: 0.5, df=29, p>0.05) in the mean δ15NTN between both the outflows (IGP: +20.4±5.4 ‰ and SEA: +19.4±6.1 ‰; Fig. 3). This feature indicates that δ15NTN over the BoB in winter season is mostly governed by the atmospheric processes affecting aerosol-TN concentrations during long-range transport rather than contribution from different sources (i.e., δ15NTN is ‘process dependent’ over the BoB). From Fig. 2, it is obvious that the concentrations of water-soluble inorganic nitrogen species ( and ) in the MABL account for more than 90 % of the TN mass. Therefore, processes affecting the relative abundances of and in aerosol-TN might control the δ15NTN (Section 3.3 for a more detailed discussion).

Fig 2
Fig. 2.   (a) Comparison of mass concentrations of total carbon (TC) measured by elemental analyser coupled to isotope ratio mass spectrometer (EA-irMS) versus estimated TC concentration based on the atmospheric abundances of elemental and organic carbon (as EC+OC, measured on Sunset Carbon Analyzer), and (b) correlation plot between total nitrogen (TN) measured by EA-irMS and water-soluble inorganic nitrogen (i.e. and ) in PM2.5 samples collected over the Bay of Bengal during January 2009 cruise.

The sudden change (t-score: 5.5, df=29; p<0.001) in wind regimes from the IGP- to SEA-outflow was characterized by distinct δ13CTC over the BoB (−23.8±0.6 ‰ and −25.3±0.9 ‰, respectively). This observation indicates that stable C-isotopic composition of aerosols over the BoB is ‘source dependent’. A significant linear relationship is noteworthy between TN and δ15NTN in both the IGP- and SEA-outflows (Fig. 4a). However, such a significant linear relationship between TC and δ13CTC is observed only in the SEA-outflow but not in the IGP-outflow (Fig. 4b). This could be due to variability in the sources of aerosol-TN and TC in the IGP-outflow (wood burning and FF-comb, respectively) than those in the SEA-outflow (forest fires). Thus, a significant depletion of 13C in TC for the SEA-outflow compared with that for the IGP-outflow, perhaps, indicates a source-specific contribution from BBEs in Southeast Asia during the cruise.

Overall, 13C enrichment and depletion of TC in the IGP- and SEA-outflows, respectively, can be explained based on the differences in the source contribution of carbonaceous species and their transport history to the BoB. A noteworthy feature of the data is evident from a change in the TC concentration in the IGP-outflow with no variability in δ13CTC Fig. 4b). This feature suggests that dominant contribution of aerosol-TC in the IGP-outflow is derived from point sources in the IGP. Based on the analyses of various aerosol properties (e.g. AOD, Angstrom exponent), derived from Multiangle Imaging SpectroRadiometer (MISR) level 3 remote sensing data, Prasad et al. (2006) have documented that emissions from coal-fired power plant in the IGP is a dominant source of absorbing organic aerosols and other airborne particulate matters in winter season. Furthermore, there exists a notable declining trend in the emission of organic aerosols from the biofuel (wood) combustion and an ongoing rise in coal combustion based on the comparison of energy consumption between 1990 and 2000 (Prasad et al., 2006). Therefore, we attribute coal-fired power plant emissions in the IGP is a major source of organic aerosols sampled over the N-BoB (i.e. in the IGP-outflow).

The forest fires, occurring over Southeast Asia during the cruise in the BoB (January 2009), contribute to large-scale emissions of POA as well as VOCs. Thus, the observed depletion in 13C of TC in the SEA-outflow (Fig. 3b) can be explained in terms of contribution of POA from C3 vegetation burning and/or from the secondary organic aerosol (SOA) formation through the photochemical oxidation of VOCs. Similar to our study, lower values of δ13CTC (−26.6 to −23.2 ‰, average −25.0 ‰) were observed for the wheat residue-burning aerosols sampled over Mt. Tai in East Asia (Fu et al., 2012). Furthermore, Fu et al. (2012) documented a non-linear (exponential) relationship between δ13CTC and levoglucosan (a tracer of BBEs), suggesting that an increased contribution of BBEs causes a significant depletion of 13C in aerosols (δ13C: −26 to −27 ‰) (Fu et al., 2012). Narukawa et al. (1999) also reported a conspicuous decrease in δ13CTC from −25.5 to −27.5 ‰ for the fine mode aerosols sampled from the Indonesian forest fires, in response to an increased contribution of organics from C3 vegetation. Also, the fractionation of stable carbon isotopes during the combustion of C3 plants is not significant (~0.5 ‰ enrichment between unburned plant and aerosols (Turekian et al., 1998).

Fig 3
Fig. 3.   Spatial distributions of (a) stable carbon (δ13CTC) and (b) nitrogen (δ15NTN) isotopic composition of fine mode (PM2.5) aerosols collected from the marine atmospheric boundary layer of the Bay of Bengal during a winter cruise (27 December 2008 to 28 January 2009).

Based on the above arguments, the lower δ13CTC in the SEA-outflow can be attributed to contribution of organic aerosols from the combustion of C3 plants. However, SOAs produced from the photooxidation of biogenic VOCs could also contribute to lower the δ13CTC (Kirillova et al., 2013). A recent laboratory photooxidation/smog chamber experiment of toluene (a VOC emitted from motor vehicles) has demonstrated formation of SOA, whose δ13C values were ~6 ‰ lower than that of the parent toluene (Irei et al., 2006, 2011). Furthermore, δ13C fractionation between gas and aerosol phases is found to be small during photooxidation (Irei et al., 2006, 2011). Similarly, the laboratory photooxidation on gaseous isoprene and ozone system results in the formation of less volatile methacrolein and methyl vinyl ketone (i.e. known gaseous isoprene oxidation products), whose δ13C values are depleted by more than 3.6 and 4.5 ‰, respectively (Iannone et al., 2010). These observations, thus, support the argument that photooxidation of anthropogenic/biogenic VOCs results in not only more oxidized semi-volatile organics for which isotopic fractionation effects are found to be minimal (Irei, 2008) but also lowers the δ13C values than those of precursor fuel/vegetation.

These semi-VOCs subsequently partition to SOAs and, hence, can explain their lower δ13C signatures. In this regard, lower δ13C values in the SEA-outflow can, therefore, be attributed to the contribution of SOAs derived from atmospheric oxidation of VOCs emitted from the BBEs in Southeast Asia. We have also observed a significant linear relationship of δ13CTC with nss-K+ (R2=0.58, p<0.05) and WSOC (R2=0.60, p<0.05) in the SEA-outflow (Fig. 4c and d) but not in the IGP-outflow (δ13CTC vs. WSOC: R2=0.01; p>0.05 and δ13CTC vs. nss-K+: R2=0.03; p>0.05). In the literature, the occurrence of nss-K+ in the fine mode aerosols and linear relationship with WSOC are used to assess the contribution of carbonaceous aerosols from the BBEs (Andreae, 1983; Andreae and Merlet, 2001; Paris et al., 2010; Pachon et al., 2014; Srinivas and Sarin, 2014). Therefore, a lack of significant linear relationship between nss-K+ and δ13CTC in the IGP-outflow sampled over the BoB (Fig. 4c), which is consistent with that of TC and δ13CTC (Fig. 4b), indicate the role of other sources of organic aerosols over the IGP (e.g. coal-fired power plants or vehicular emissions) than BBEs. Based on the literature, SOAs contribute significantly to atmospheric WSOC (Graham et al., 2002; Mayol-Bracero et al., 2002b; Sciare et al., 2008; Ram and Sarin, 2010; Hegde and Kawamura, 2012; Rajput et al., 2014; Srinivas and Sarin, 2014). In this regard, the BBE-derived POA and/or SOA causes 13C depletion of TC in the SEA-outflow, which can also be supported by the strong correlations of δ13C with nss-K+ and WSOC (Fig. 4c and d).

Fig 4
Fig. 4.   Correlation plots for (a) TN versus δ15NTN, (b) TC versus δ13CTC, (c) nss-K+ versus δ13CTC and (d) WSOC versus δ13CTC for fine mode aerosols collected over the Bay of Bengal during a winter cruise (27 December 2008 to 28 January 2009).

The impact of BBEs over FF-comb sources over the BoB can also be inferred based on the relative increase in mass ratio of nss-K+ to EC (Table 1) as the latter contains negligible potassium content (Novakov et al., 2000; Mayol-Bracero et al., 2002b). Furthermore, higher nss-K+/EC ratios in the SEA-outflow (0.39±0.09) compared with those in the IGP-outflow (0.29±0.09) suggest the dominance of BBEs in the former samples. A significant linear relationship is noteworthy between δ13CTC and δ15NTN in the SEA-outflow (slope=− 5.5; R2=0.63; p<0.05; Supplementary Fig. 3); however, no such correlation is observed for the IGP-outflow (slope=− 3.5; R2=0.15; p=0.15; Supplementary Fig. 3). The observed weak (or no) and strong correlation of δ13CTC with δ15NTN for the IGP- and SEA-outflows, respectively, could be due to variable contribution of TC and TN from FF-comb and BBEs over the BoB. Furthermore, moderate linear correlation is found between nss- and EC for the IGP-outflow samples (R2=0.57; p<0.05), which indicates an impact of FF-combustion sources over the north BoB. However, rather weak linear correlation was observed for SEA-outflow between nss- and EC (R2=0.26; p=0.03). This observation emphasizes a significant contribution of TC from FF-comb sources in the IGP over the BoB during the study period.

3.2. Relative contribution of fossil fuel versus biomass/biofuel-burning source

The δ13CTC in PM2.5 over the BoB during a winter cruise can be explained based on two end-member mixing between FF-comb and BBEs. Since the IGP is a densely populated and industrialized region of India with intense ongoing anthropogenic activities (viz., coal-fired power plants, vehicular emissions, wood/post-harvest crop-residue burning), it is not logical to assume that biogenic emissions mainly contribute to aerosol-TC in the IGP-outflow. The dominant contribution of aerosol-TC in the SEA-outflow samples is from forest fires occurring in Southeast Asia (as inferred from MODIS satellite-derived fire count data; Supplemantary Fig. 3).

Organic compounds, derived from phytoplankton debris and zooplankton excreta in the sea surface microlayer, are coated on sea-salt particles, contributing to marine aerosols collected over the BoB. Therefore, marine phytoplankton-derived organic aerosols have an average δ13C value of ca. −22 ‰ (Miyazaki et al., 2011). However, the BoB is not a productive oceanic basin and is perennially influenced by fresh water influx from the Ganga–Brahmaputra and other peninsular rivers. Moreover, the persistence of weak NE-monsoonal winds (<3 m second−1) are not likely to trigger sea-to-air transfer of organics from the surface waters of the BoB during the continental outflow. Therefore, contribution from marine organic matter is unlikely due to prevailing winds during the study period. This is further supported by the relative contribution of sea-salts and hence particulate organic matter from the ocean surface by bubble bursting to ambient aerosols over the BoB is no more than 5 % (Srinivas and Sarin, 2014). Therefore, we believe that a two end-member mixing model is a reasonable first-order pseudo-approximation to assess the relative contribution of aerosol-TC from BBEs and FF-comb over the BoB.

The rationale for the contribution from these two sources is further supported by the observed nss-K+/EC and nss-/EC ratios in PM2.5 over the BoB (for additional information, see Section 3.1). Higher nss-K+/EC ratios in the SEA-outflow than those in the IGP-outflow samples indicate a dominance of BBEs over the South BoB (Table 1) (Srinivas and Sarin, 2014). The nss-/EC ratio, a proxy for fossil fuel combustion with higher values documented for coal-fired power plants (~6) (Ramana et al., 2010), showed somewhat consistency between the IGP- and SEA-outflows (Srinivas and Sarin, 2014) (Table 1), indicating a contribution from point sources in the IGP (e.g. coal-fired power plants). Also, the 7-d AMBTs for the PM2.5 samples collected over the north BoB have the source origin from the IGP (Fig. 1), which intercept the receptor sites dominated by carbonaceous aerosols (and hence TC) derived from coal-fired power plants (see blue pins in Fig. 1).

Widory (2006) and Mori et al. (1999) have documented that δ13C of aerosols from coal combustion as −23.6±0.7 ‰ (−22.9 to −24.9 ‰) and −23.4±1.2 ‰, respectively, with almost zero fractionation (or slightly negative). On the other hand, combustion-derived particles from liquid fossil fuels (e.g. unleaded/leaded gasoline; δ13C: −26.4 to −28.0 ‰) are depleted in 13C relative to those originated from solid fuel (coal: −24 to −21 ‰) (Widory, 2006) but are enriched with 13C relative to gaseous fuels (−40 to −28 ‰) (Huang et al., 2006; Widory, 2006; López-Veneroni, 2009; Agnihotri et al., 2011; Kawashima and Haneishi, 2012). A comparison of δ13CTC (−23.8±0.6 ‰) over the north BoB with those representative of various fossil fuels as mentioned above (coal, diesel, gasoline, etc.) suggest a dominant contribution from coal-fired power plant emissions in the IGP (Fig. 5).

Fig 5
Fig. 5.   Comparison of δ13C of aerosol-TC over the BoB during a winter cruise (January 2009) with literature values of various sources such as combustion/emissions from C3, C4 and marine phytoplankton (Turekian et al., 1998; Cao et al., 2013), coal combustion (Mori et al., 1999; Widory, 2006) and vehicular emissions (Widory, 2006; Cao et al., 2013).

Although we have observed a decrease in the number (intensity) of fire counts (derived from MODIS data) over the IGP than those over SEA, the contribution of wood-burning emissions to aerosol-TC (Gustafsson et al., 2009; Ram and Sarin, 2012) cannot be ruled out in the IGP-outflow sampled over the BoB. In this context, laboratory and natural field burning chars from C3 grasses and trees (wood and leaves) are enriched and/or depleted in 13C (e.g. δ13C for Eucalyptus wood: −27.2 ‰; Sphagnum: −26.0 ‰; C3 grasses: −28.2 to −32.4 ‰; Eucalyptus leaf: −28.2 ‰) compared with plant materials (Krull et al., 2003). On a similar note, based on the laboratory-based combustion, Agnihotri et al. (2011) documented that aerosols from typical BBEs in India exhibit a very narrow range of δ13C (−27.9±1.1 ‰). However, the observed δ13CTC in the IGP-outflow (−23.8±0.6 ‰) is consistent with that of coal-derived emissions compared with other sources (i.e. gasoline combustion and wood-burning emissions; Fig. 5. Therefore, we have used coal-burning emissions as one end-member (δ13CTC: −23 ‰; Widory, 2006) for explaining the observed stable carbon isotopic signatures of TC in the continental outflow sampled over the BoB during winter. We can further support the significance of coal-fired power plants as a potential source of organic aerosols in the IGP-outflow sampled over the BoB as follows.

As mentioned earlier, FF-comb (coal combustion and vehicular emissions: VEs) and BBEs (wood/crop-residue) in the IGP contribute to organic aerosols over the BoB during the winter period (Rengarajan et al., 2007; Sudheer and Sarin, 2008; Kumar et al., 2010; Srinivas and Sarin, 2013, 2014). In such cases, the δ13CTC of PM2.5 over the BoB can be expressed as a product of individual fractions of TC from each of these three respective sources times their stable C-isotopic composition as below

Based on the 14C-analysis of EC, Gustafsson et al. (2009) documented near equal contribution of wood-burning emissions and fossil fuel combustion sources to Atmospheric Brown Clouds over South Asia. Therefore, we have assumed 50 % contribution from biofuels (Gustafsson et al. 2009), 30 % from coal-fired thermal power plants and 20 % from vehicular emissions (Prasad et al., 2006) to atmospheric organic aerosols in the IGP-outflow. Agnihotri et al. (2011) documented that source-specific δ13CTC of various biofuels and crop-residue combustion-derived aerosols from South Asia centred on −27.9±1.1 ‰. Similarly, a source-specific δ13CTC of −23.6±0.7 ‰ from coal combustion-derived particles has been documented by Widory (2006). Also, the δ13CTC of particles generated from diesel combustion emissions span from −28 to −26 ‰ (Widory, 2006).

Using the respective source contributions and their corresponding δ13CTC in eq. (1), we have estimated δ13C of organic aerosols collected in the IGP-outflow to be −26 to −27 ‰. Based on the narrow spread of δ13CTC in the IGP-outflow (−23.8±0.6 ‰), we estimated δ13CTC from the contribution of above-mentioned three major sources in the IGP. However, it may be argued that atmospheric oxidation, and hence, enrichment in 13C could be compensated by SOA formation (with overall 13C depletion) during transport. Therefore, the observed overlap of δ13CTC over the N-BoB with coal-fired power plant emissions in the IGP emphasizes the source significance in the IGP-outflow. However, further studies on seasonal variability of the IGP-outflow are needed along with the δ13CTC.

The δ13CTC in the SEA-outflow (−26.9 to −24.2 ‰) overlap with those derived from both vehicular emissions and forest fires (Fig. 5). However, the SEA-outflow samples are collected over the pelagic (deep) waters (i.e. away from the coast and unlike samples from the IGP-outflow, where vehicular emissions could contribute) and, therefore, not likely to have contribution from vehicular emissions. Furthermore, AMBTs originating from Southeast Asia have a significant contribution from forest fires, inferred from the MODIS fire count data (Supplemantary Fig. 2). Therefore, we have considered BBEs as another end-member to account for the variability in δ13CTC in the SEA-outflow.

Aerosol-TC in the SEA-outflow (−25.3±0.9 ‰) is significantly depleted in 13C than that derived from C4 plants (−17 to −9 ‰; average −13 ‰), being somewhat consistent with those originated from C3 vegetation burning (−32 to −20 ‰; average −28 ‰) (Smith and Epstein, 1971). In a previous study, Turekian et al. (1998) suggested that particles from C3 plant burning emissions showed a slight enrichment in 13C in aerosol-TC (or higher in δ13C by only 0.5 ‰ or less) than host plant, while C4 plants resulted in a depletion of 13C (or lower in δ13C by 3.5 ‰) (Turekian et al., 1998). According to Turekian et al. (1998), the increased relative proportion of refractory components (e.g. lipids; 13C depleted than bulk plant) over labile carbon pools (13C enrichment) in aerosols during C4 plant combustion causes 13C depletion of TC; however, no sucheffects were observed for C3 plant combustion-derived aerosols. A similar result has been reported for the smoke and ash that are derived from the combustion of C3 plants which show negligible (or no) fractionation (Krull et al., 2003; Das et al., 2010).

Taken together of all these findings, we have assumed −23 and −28 ‰ to be the representative end-member δ13C signatures for FF-comb source and BBEs, respectively, based on the two end-member mixing model, which is described as follows:

Here, fFF-comb and fBBEs refer to the fractional contribution of TC from fossil fuel combustion and biomass-burning emission sources, respectively. The end-member δ13C for FF-comb and BBEs are −23.0 and −28.0 ‰, respectively. Using eq. (4), we estimated the contribution of fossil fuel combustion to both the IGP- and SEA-outflows. In this study, the estimated contribution of FF-comb source to δ13CTC varied from 61 to 100 % (average 84±12 %) for the IGP-outflow samples and 12 to 66 % (average 44±18 %) for the SEA-outflow samples.

3.3. Elevated δ15N of TN: transformation of gaseous NH3 to particulate

In both the IGP- and SEA-outflows, the water-soluble inorganic nitrogen (i.e. -N accounting for 81–100 %) dominates the total soluble nitrogen (Srinivas et al., 2011b). Further, a moderate-to-strong linear correlation of δ15NTN with and equivalent mass ratio of / is observed in both continental outflows (Fig. 6). Since dominates aerosol-TN over the BoB (i.e. >90 %; R2=0.95; Fig. 2), atmospheric processes affecting the concentration of particulate in fine mode aerosols might explain the observed δ15NTN. The particulate concentration in the MABL, however, is influenced by the atmospheric reactions of gaseous NH3 with acidic species (e.g. H2SO4, HNO3) to form NH4HSO4, (NH4)2SO4 and NH4NO3. Furthermore, NH4NO3 is formed mostly under ammonia-rich conditions (i.e. for which to equivalent mass ratio is greater than unity) because the reaction between ammonia and sulphate is preferable. However, NH4NO3 is thermodynamically not stable and, therefore, often dissociated into precursor chemical constituents during atmospheric transport, which also affects the ambient concentration. As a result, mainly exists in coarse mode due to the reaction with mineral dust and sea-salt (Jickells, 2006). All these processes/reactions could influence δ15NTN in ambient aerosols. It is implicit from Fig. 6 that linear relationship of concentration and the equivalent mass ratios of to with δ15NTN can explain the observed enrichment of 15N in the PM2.5 sampled over the BoB.

Fig 6
Fig. 6.   Linear regression analysis between (a) concentration of and δ15N, and (b) equivalent mass ratio of to and δ15N in PM2.5 collected over the Bay of Bengal during a winter cruise (January 2009).

The mean equivalent mass ratio of / in the samples from IGP- and SEA-outflow corresponds to 0.92±0.15 and 0.77±0.12, respectively, suggesting an ammonia-deficit (poor) condition. Moreover, the observed ambient temperature during the cruise averages around 26.7 °C (IGP-outflow: 24.4–27.4 °C; SEA-outflow: 25.9–29.0 °C) and high nss- concentrations facilitate the dissociation of NH4NO3. Interestingly, Srinivas et al. (2011b) have documented the predominant occurrence of in coarse mode and in fine mode aerosols during the January 2009 cruise. This implies that particulate in the MABL mainly occur due to the formation of more stable NH4HSO4 or (NH4)2SO4. In this regard, we observed a very strong correlation between and (Srinivas et al., 2011b), suggesting the formation (NH4)2SO4 and NH4HSO4 over the BoB during the study period. Therefore, we believe that atmospheric transformation of gaseous NH3 to particulate might be responsible for the observed large variability in δ15N of TN over the BoB.

Similar to our study, Pavuluri et al. (2010) reported an elevated δ15N of TN and its significant positive linear correlation with sulphate to ammonium molar ratio (or negative correlation with that of / molar ratio) for the urban coastal aerosols from Chennai, India, during the late winter (January). Further, they observed a positive relationship between ambient temperature and δ15N; in particular, higher δ15N are noteworthy for those sampling days with temperatures >20 °C (Pavuluri et al., 2010). The ambient temperature range over the BoB in this study is higher than that documented in Pavuluri et al. (2010). During a laboratory experiment, an exposure of gaseous NH3 to aerosol filters loaded with sulphuric acid results in the formation of (NH4)2SO4 (i.e. ascertained through stoichiometric equilibrium between and ), which is enriched in 15N with isotope enrichment factor of +33 ‰ (Heaton et al., 1997). Combining all the above arguments, we believe that gas to particle formation of NH3 to particulate in the MABL can explain the enrichment of 15N in the PM2.5 sampled over the BoB during winter.

The mean value of δ15NTN in the IGP-outflow sampled over the BoB during a winter cruise (this study) is higher than those obtained during a spring–intermonsoon cruise [+10.6±2.7 ‰ (Agnihotri et al., 2011)]. Although the continental outflow from the IGP persists during spring intermonsoon, the observed depletion of 15N could be explained based on the decrease in source strength/variability. The weakening of the IGP-outflow during spring–intermonsoon, compared with winter months (Table 1), is inferred by a decrease in the atmospheric concentrations of (0.01–2.3 µg m−3; 0.9±0.6 µg m−3) and nss- (2.5–10.3 µg m−3; 10.3±2.5 µg m−3), and equivalent mass ratios of to nss- (0.01–0.66; average 0.41±0.20) over the BoB. Therefore, we believe that the observed decrease in particulate in the marine aerosols during spring intermonsoon might indicate the non-availability of gaseous NH3 to form NH4HSO4 and (NH4)2SO4 aerosols and, thus, explain the observed lower δ15NTN.

Similar to our study, higher δ15N values (+13.4 to +22.1 ‰) have been reported for fine mode aerosols sampled over Tanzania during dry season, which is attributed to a significant 15N enrichment by isotopic exchange reaction of N between gaseous NH3 and aerosol during atmospheric transport (Mkoma et al., 2014). Similarly, Pavuluri et al. (2010) have also reported similar δ15N values in PM10 aerosols, ranging from +18.0 to +27.8 ‰ and +18.6 to +25.7 ‰ in day and night samples, respectively, collected during winter from a coastal site in India (Chennai). Therefore, observed variability in δ15NTN for the IGP-outflow samples (+11.8 to +30.6 ‰) can be attributed to pronounced continental impact of anthropogenic sources (reflected from their higher ambient abundances of and over the north BoB) to TN and also the exchange reactions between gaseous NH3 and particulate during transport. It is important to emphasize that δ13CTC and δ15NTN showed a good correlation for the SEA-outflow (influenced by BBEs; see Fig. 1), but not in the IGP-outflow. Similar to this study, Agnihotri et al. (2011) have documented a strong correlation between δ13C and δ15N over the BoB during spring–intermonsoon cruise, which is also influenced by the BBEs. In a previous study, Turekian et al. (1998) documented that various stages of biomass combustion resulted in a large spread in stable nitrogen isotopic signatures of aerosol-TN. Therefore, variability in δ15NTN in the SEA-outflow (+10.4to +30.7 ‰; average +19.4±6.1 ‰) can be explained by the highly variable temperature conditions that prevail during the forest fires in Southeast Asia.

3.4. Comparison with studies from the Northern Indian Ocean

The stable isotopic composition of TC and TN from this study (January 2009) was compared with those from earlier cruises conducted in the NIO (Agnihotri et al., 2011) and other sampling sites from South Asia. A significant difference in the δ13CTC is noteworthy during the cruises conducted in the BoB in January 2009 (Table 2) and March–April 2006 [range: −26.6 to −24.1 ‰; average −25.6±0.6 ‰ for TSP samples (Agnihotri et al., 2011)]. Although δ13C of TC for both cruises refers to different size fractions of aerosols, it is relevant to state that during long-range transport the mass size distribution of atmospheric particulates changes from coarse to fine mode due to the settling of coarser particles (Duce et al., 1991). This is most relevant for our sampling from the MABL of the BoB and the chemical composition mostly comprised dominant fine particles (Srinivas et al., 2011b). Although source strength of carbonaceous aerosols in the IGP-outflow is varied in between the cruises (i.e. higher abundances in winter than spring–intermonsoon) over the BoB, their sources remain the same, as inferred from the diagnostic mass ratios of chemical composition of aerosols (Sarin et al., 2010; Srinivas and Sarin, 2012).

Table 2. Comparison of stable carbon and nitrogen isotopic compositions of aerosols from this study with those documented in the literature over South Asia
Region Season Sample type δ13C (‰) δ 15N (‰) Reference
Bay of Bengal Winter (IGP-outflow) PM2.5 −23.8±0.6 20.4±5.4 This study
Bay of Bengal Winter (SEA-outflow) PM2.5 −25.8±0.9 19.4±6.1 This study
Bay of Bengal Spring–intermonsoon TSP −25.6±0.6 10.6±2.7 Agnihotri et al. (2011)
Arabian Sea Spring–intermonsoon TSP −26.5±0.8 1.4±3.2 Agnihotri et al. (2011)
Mumbai Summer TSP −26.5±0.3 20.2±1.2 Aggarwal et al. (2013)
Mumbai Winter TSP −25.9±0.3 22.8±1.4 Aggarwal et al. (2013)
Chennai Summer and winter PM10 −25.0±0.6 Pavuluri et al. (2011)
Chennai Winter (day) PM10 25.5±2.4 Pavuluri et al. (2010)
Chennai Winter (night) PM10 23.0±2.3 Pavuluri et al. (2010)
Chennai Summer (day) PM10 22.2±4.3 Pavuluri et al. (2010)
Chennai Summer (night) PM10 24.3±3.6 Pavuluri et al. (2010)
Africa Summer PM10 −25.9±0.3 13.7±2.2 Mkoma et al. (2014)

The δ13CTC in the IGP-outflow from a winter cruise (January 2009) in the BoB (this study) is significantly higher than that obtained during spring–intermonsoon (Agnihotri et al., 2011) (Table 2). However, δ13CTC in the SEA-outflow sampled over the BoB in winter (11–28 January, 2009) is somewhat similar with that from spring–intermonsoon cruise (Agnihotri et al., 2011). The observed differences in δ13CTC during the two cruises could be related to temporal shift in the source contributions (anthropogenic aerosols accounting for more than 80 % in winter vis-à-vis near equal contribution of dust and anthropogenic species during spring–intermonsoon, March–April 2006) (Srinivas and Sarin, 2012). Analysis of MODIS fire images (Supplementary Fig. 3) clearly showed hotspots over the IGP during spring–intermonsoon (March–April 2006) and Southeast Asia during winter (January 2009; Fig. 1). Therefore, 13C depleted values of TC in the IGP-outflow samples in spring–intermonsoon and the SEA-outflow samples in winter can be explained by a significant contribution of carbonaceous species from BBEs. The impact of BBEs in the IGP-outflow sampled over the BoB during spring–intermonsoon and the SEA-outflow during a winter cruise was also reflected through consistent nss-K+/EC ratios (0.35±0.11and 0.39±0.09, respectively).

We have also observed consistency in the δ13CTC from this study with those documented for wintertime aerosols (−25.9±0.3 ‰) (Aggarwal et al., 2013) from Mumbai (a continental site in India). The biofuel/biomass burning is a dominant source of TC over Mumbai (Venkataraman et al., 2002, 2005) and also over the IGP during winter (Rengarajan et al., 2007; Ram and Sarin, 2012; Rajput et al., 2014; Srinivas and Sarin, 2014). Therefore, the overlapping δ13CTC in the SEA-outflow samples in winter with those of the IGP-outflow samples from the BoB during spring–intermonsoon and of continental aerosols from Mumbai might correspond to biofuel/biomass-burning emissions.

4. Conclusions

This study reports the stable C- and N-isotopic composition of wintertime aerosols (PM2.5) over the BoB. During the study period (December 2008–January 2009), the MABL is influenced by the continental outflow from the IGP- and SEA-outflows, respectively, as inferred from the 7-d AMBTs over the BoB. The δ13CTC in the IGP-outflow (−23.8±0.6 ‰) is distinctly different than that in the SEA-outflow (−25.3±0.9 ‰). A significant linear relationship of δ13CTC with TC, nss-K+ and WSOC is noteworthy in the SEA-outflow but not in the IGP-outflow, which indicates the impact of BBEs in the former samples (which is also supported by the MODIS fire count data). However, no such influence from BBEs was observed for the IGP-outflow. Besides, the temporal variability of δ13CTC in the IGP-outflow is rather small unlike TC concentration and overlaps with those derived from coal combustion sources (−23.6±0.7 ‰) in the IGP. The observed 13C-depletion of aerosol-TC in the SEA-outflow is attributed to the contribution of POA directly emitted from the BBEs or due to SOAs formed by the oxidation of BVOCs from the forest fires. We observe no significant differences in the δ15NTN between IGP- and SEA-outflows. Moreover, mass concentration of dominates the aerosol-TN over the BoB, showing significant linear relationship with δ15NTN. Thus, we infer that δ15NTN in aerosols over the BoB during winter is governed by the atmospheric processes affecting the concentrations of aerosol water-soluble inorganic nitrogen species (e.g. , ) during transport.

5. Acknowledgements

The authors thank all the crewmembers of ORV Sagar Kanya for their help during aerosol collection. The authors acknowledge the partial financial support provided from the ISRO-Geosphere Biosphere Programme. The logistic support provided by Drs. K. Krishnamoorthy and C.B.S. Dutt is thankfully acknowledged. The financial support provided by the Japanese Society for the Promotion of Science (JSPS) through grant-in-aid no. 24221001 is acknowledged for carrying out C and N isotope analyses at Hokkaido University.


Aggarwal, S. G. and Kawamura, K. 2008. 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. Publisher Full Text

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. Publisher Full Text

Agnihotri, R., Mandal, T. K., Karapurkar, S. G., Naja, M., Gadi, R. and co-authors. 2011. Stable carbon and nitrogen isotopic composition of bulk aerosols over India and northern Indian Ocean. Atmos. Environ. 45, 2828–2835. Publisher Full Text

Andreae, M. O. 1983. Soot carbon and excess fine potassium: long-range transport of combustion-derived aerosols. Science. 220, 1148–1151. Publisher Full Text

Andreae, M. O. and Merlet, P. 2001. Emission of trace gases and aerosols from biomass burning. Glob. Biogeochem. Cycles. 15, 955–966. Publisher Full Text

Cao, F., Zhang, S.-C., Kawamura, K. and Zhang, Y.-L. 2016. Inorganic markers, carbonaceous components and stable carbon isotope from biomass burning aerosols in Northeast China. Sci. Total Environ. DOI:

Cao, J.-J., Zhu, C.-S., Tie, X.-X., Geng, F.-H., Xu, H.-M. and co-authors. 2013. Characteristics and sources of carbonaceous aerosols from Shanghai, China. Atmos. Chem. Phys. 13, 803–817. Publisher Full Text

Das, O., Wang, Y. and Hsieh, Y.-P. 2010. Chemical and carbon isotopic characteristics of ash and smoke derived from burning of C3 and C4 grasses. Org. Geochem. 41, 263–269. Publisher Full Text

Draxler, R. R., Stunder, B., Rolph, G. and Taylor, A. 1999. HYSPLIT4 users guide. NOAA Tech. Memo. ERL ARL. 230, 35.

Duce, R., Liss, P., Merrill, J., Atlas, E., Buat-Menard, P. and co-authors, A. 1991. The atmospheric input of trace species to the world ocean. Glob. Biogeochem. Cycles. 5, 193–259. Publisher Full Text

Fu, P., Kawamura, K., Chen, J., Li, J., Sun, Y. and co-authors. 2012. Diurnal variations of organic molecular tracers and stable carbon isotopic composition in atmospheric aerosols over Mt. Tai in the North China Plain: an influence of biomass burning. Atmos. Chem. Phys. 12, 8359–8375. Publisher Full Text

Graham, B., Mayol-Bracero, O. L., Guyon, P., Roberts, G. C., Decesari, S. and co-authors. 2002. Water-soluble organic compounds in biomass burning aerosols over Amazonia 1. Characterization by NMR and GC-MS. J. Geophys. Res. Atmos. 107, LBA 14-11–LBA 14-16. Publisher Full Text

Gustafsson, O., Krusa, M., Zencak, Z., Sheesley, R. J., Granat, L. and co-authors. 2009. Brown clouds over South Asia: biomass or fossil fuel combustion? Science. 323, 495–498. Publisher Full Text

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. Publisher Full Text

Hegde, P. and Kawamura, K. 2012. Seasonal variations of water-soluble organic carbon, dicarboxylic acids, ketocarboxylic acids, and α-dicarbonyls in Central Himalayan aerosols. Atmos. Chem. Phys. 12, 6645–6665. Publisher Full Text

Huang, L., Brook, J., Zhang, W., Li, S., Graham, L. and co-authors, A. 2006. Stable isotope measurements of carbon fractions (OC/EC) in airborne particulate: a new dimension for source characterization and apportionment. Atmos. Environ. 40, 2690–2705. Publisher Full Text

Iannone, R., Koppmann, R. and Rudolph, J. 2010. Stable carbon kinetic isotope effects for the production of methacrolein and methyl vinyl ketone from the gas-phase reactions of isoprene with ozone and hydroxyl radicals. Atmos. Environ. 44, 4135–4141. Publisher Full Text

Irei, S., Huang, L., Collin, F., Zhang, W., Hastie, D. and co-authors. 2006. Flow reactor studies of the stable carbon isotope composition of secondary particulate organic matter generated by OH-radical-induced reactions of toluene. Atmos. Environ. 40, 5858–5867. Publisher Full Text

Irei, S. 2008. Laboratory Studies of Stable Carbon Isotope Ratio of Secondary Particulate Organic Matter in the Gas-Phase. PhD Dissertation, Department of Chemistry, York University, Toronto, Canada.

Irei, S., Rudolph, J., Huang, L., Auld, J. and Hastie, D. 2011. Stable carbon isotope ratio of secondary particulate organic matter formed by photooxidation of toluene in indoor smog chamber. Atmos. Environ. 45, 856–862. Publisher Full Text

Jacobson, M. Z. 2012. Investigating cloud absorption effects: global absorption properties of black carbon, tar balls, and soil dust in clouds and aerosols. J. Geophys. Res. 117, D06205. Publisher Full Text

Jayaraman, A., Lubin, D., Ramachandran, S., Ramanathan, V., Woodbridge, E. and co-authors. 1998. Direct observations of aerosol radiative forcing over the tropical Indian Ocean during the January–February 1996 pre-INDOEX cruise. J. Geophys. Res. 103, 13827–13836. Publisher Full Text

Jickells, T. 2006. The role of air-sea exchange in the marine nitrogen cycle. Biogeosci. Discuss. 3, 183–210. Publisher Full Text

Jimenez, J. L., Canagaratna, M. R., Donahue, N. M., Prevot, A. S. H., Zhang, Q. and co-authors. 2009. Evolution of organic aerosols in the atmosphere. Science. 326, 1525–1529. Publisher Full Text

Kanakidou, M., Seinfeld, J. H., Pandis, S. N., Barnes, I., Dentener, F. J. and co-authors. 2005. Organic aerosol and global climate modelling: a review. Atmos. Chem. Phys. 5, 1053–1123. Publisher Full Text

Kawamura, K., Kobayashi, M., Tsubonuma, N., Mochida, M., Watanabe, T. and co-authors. 2004. Organic and inorganic compositions of marine aerosols from East Asia: seasonal variations of water-soluble dicarboxylic acids, major ions, total carbon and nitrogen, and stable C and N isotopic composition. In: The Geochemical Society Special Publications (eds. R. J. Hill et al.). Elsevier Science, New York, pp. 243–265.

Kawashima, H. and Haneishi, Y. 2012. Effects of combustion emissions from the Eurasian continent in winter on seasonal δ13C of elemental carbon in aerosols in Japan. Atmos. Environ. 46, 568–579. Publisher Full Text

Kedia, S., Ramachandran, S., Kumar, A. and Sarin, M. 2010. Spatiotemporal gradients in aerosol radiative forcing and heating rate over Bay of Bengal and Arabian Sea derived on the basis of optical, physical, and chemical properties. J. Geophys. Res. 115, D07205. Publisher Full Text

Keene, W. C., Pszenny, A. A. P., Galloway, J. N. and Hawley, M. E. 1986. Sea-salt corrections and interpretation of constituent ratios in marine precipitation. J. Geophys. Res. 91, 6647–6658. Publisher Full Text

Kirillova, E. N., Andersson, A., Sheesley, R. J., Kruså, M., Praveen, P. and co-authors. 2013. 13C- and 14C-based study of sources and atmospheric processing of water-soluble organic carbon (WSOC) in South Asian aerosols. J. Geophys. Res. 118, 614–626.

Krull, E. S., Skjemstad, J. O., Graetz, D., Grice, K., Dunning, W. and co-authors. 2003. 13C-depleted charcoal from C4 grasses and the role of occluded carbon in phytoliths. Org. Geochem. 34, 1337–1352. Publisher Full Text

Kumar, A., Sudheer, A. and Sarin, M. 2008. Chemical characteristics of aerosols in MABL of Bay of Bengal and Arabian Sea during spring inter-monsoon: a comparative study. J. Earth Syst. Sci. 117, 325–332. Publisher Full Text

Kumar, A., Sarin, M. M. and Srinivas, B. 2010. Aerosol iron solubility over Bay of Bengal: role of anthropogenic sources and chemical processing. Mar. Chem. 121, 167–175. Publisher Full Text

Kundu, S., Kawamura, K., Andreae, T. W., Hoffer, A. and Andreae, M. O. 2006. Diurnal variation in the water-soluble inorganic ions, organic carbon and isotopic compositions of total carbon and nitrogen in biomass burning aerosols from the LBA-SMOCC campaign in Rondônia, Brazil. J. Aerosol Sci. 41, 118–133. Publisher Full Text

Kundu, S., Kawamura, K. and Lee, M. 2010. Seasonal variation of the concentrations of nitrogenous species and their nitrogen isotopic ratios in aerosols at Gosan, Jeju Island: implications for atmospheric processing and source changes of aerosols. J. Geophys. Res. 115, D20305. Publisher Full Text

Lawrence, M. G. and Lelieveld, J. 2010. Atmospheric pollutant outflow from southern Asia: a review. Atmos. Chem. Phys. 10, 11017–11096. Publisher Full Text

Lelieveld, J., Crutzen, P., Ramanathan, V., Andreae, M., Brenninkmeijer, C. and co-authors. 2001. The Indian Ocean experiment: widespread air pollution from South and Southeast Asia. Science. 291, 1031–1036. Publisher Full Text

López-Veneroni, D. 2009. The stable carbon isotope composition of PM2.5 and PM10 in Mexico City Metropolitan Area air. Atmos. Environ. 43, 4491–4502. Publisher Full Text

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

Mayol-Bracero, O., Gabriel, R., Andreae, M., Kirchstetter, T., Novakov, T. and co-authors. 2002a. Carbonaceous aerosols over the Indian Ocean during the Indian Ocean Experiment (INDOEX): chemical characterization, optical properties, and probable sources. J. Geophys. Res. 107, INX2 29-21–INX22 29-21.

Mayol-Bracero, O., Guyon, P., Graham, B., Roberts, G., Andreae, M. and co-authors. 2002b. 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. 107, 8091. Publisher Full Text

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

Mkoma, S. L., Kawamura, K., Tachibana, E. and Fu, P. 2014. Stable carbon and nitrogen isotopic compositions of tropical atmospheric aerosols: sources and contribution from burning of C 3 and C 4 plants to organic aerosols. Tellus B. 66, 20176. DOI:

Mori, I., Nishikawa, M., Quan, H. and Iwasaka, Y. 1999. Regional characteristics of Chinese atmospheric aerosols in terms of the stable carbon isotope ratios. In: Proceedings of the Abstract of Annual Meeting of Japan Society for Atmospheric Environment, Tokyo, 40, p. 419. (in Japanese).

Nair, V. S., Moorthy, K. K., Alappattu, D. P., Kunhikrishnan, P., George, S. and co-authors. 2007. Wintertime aerosol characteristics over the Indo-Gangetic Plain (IGP): impacts of local boundary layer processes and long-range transport. J. Geophys. Res. 112, D13205.

Narukawa, M., Kawamura, K., Takeuchi, N. and Nakajima, T. 1999. Distribution of dicarboxylic acids and carbon isotopic ratios in aerosols from 1997 Indonesian forest fires. Geophys. Res. Lett. 26, 3101–3104. Publisher Full Text

Narukawa, M., Kawamura, K., Li, S. M. and Bottenheim, J. 2008. Stable carbon isotopic ratios and ionic composition of the high-Arctic aerosols: an increase in (13C values from winter to spring. J. Geophys. Res. 113, D18307. Publisher Full Text

Novakov, T., Andreae, M., Gabriel, R., Kirchstetter, T., Mayol-Bracero, O. and co-authors. 2000. Origin of carbonaceous aerosols over the tropical Indian Ocean: biomass burning or fossil fuels? Geophys. Res. Lett. 27, 4061–4064. Publisher Full Text

Ometto, J. P., Ehleringer, J. R., Domingues, T. F., Berry, J. A., Ishida, F. Y. and co-authors. 2006. The stable carbon and nitrogen isotopic composition of vegetation in tropical forests of the Amazon Basin, Brazil. Biogeochemistry. 79, 251–274. Publisher Full Text

Pachon, J. E., Weber, R. J., Zhang, X., Mulholland, J. A. and Russell, A. G. 2014. Revising the use of potassium (K) in the source apportionment of PM2.5. Atmos. Pollut. Res. 4, 14–21. Publisher Full Text

Paris, R., Desboeufs, K., Formenti, P., Nava, S. and Chou, C. 2010. Chemical characterisation of iron in dust and biomass burning aerosols during AMMA-SOP0/DABEX: implication for iron solubility. Atmos. Chem. Phys. 10, 4273–4282. Publisher Full Text

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. Publisher Full Text

Pavuluri, C. M., Kawamura, K., Swaminathan, T. and Tachibana, E. 2011. Stable carbon isotopic compositions of total carbon, dicarboxylic acids and glyoxylic acid in the tropical Indian aerosols: implications for sources and photochemical processing of organic aerosols. J. Geophys. Res. 116, D18307. Publisher Full Text

Prasad, A. K., Singh, R. P. and Kafatos, M. 2006. Influence of coal based thermal power plants on aerosol optical properties in the Indo-Gangetic basin. Geophys. Res. Lett. 33, L05805. Publisher Full Text

Rajput, P., Sarin, M. M., Rengarajan, R. and Singh, D. 2011. Atmospheric polycyclic aromatic hydrocarbons (PAHs) from post-harvest biomass burning emissions in the Indo-Gangetic Plain: isomer ratios and temporal trends. Atmos. Environ. 45, 6732–6740. Publisher Full Text

Rajput, P., Sarin, M., Sharma, D. and Singh, D. 2014. Characteristics and emission budget of carbonaceous species from post-harvest agricultural-waste burning in source region of the Indo-Gangetic Plain. Tellus B. 66, 21026. DOI: Publisher Full Text

Ram, K. and Sarin, M. 2010. Spatio-temporal variability in atmospheric abundances of EC, OC and WSOC over Northern India. J. Aerosol Sci. 41, 88–98. Publisher Full Text

Ram, K., Sarin, M. and Tripathi, S. 2010. A 1 year record of carbonaceous aerosols from an urban site in the Indo-Gangetic Plain: characterization, sources, and temporal variability. J. Geophys. Res. 115, D24313. Publisher Full Text

Ram, K. and Sarin, M. 2012. Carbonaceous aerosols over Northern India: sources and spatio-temporal variability. Proc. Indian Natl. Sci. Acad. 78, 523–533.

Ramana, M. V., Ramanathan, V., Feng, Y., Yoon, S. C., Kim, S. W. and co-authors. 2010. Warming influenced by the ratio of black carbon to sulphate and the black-carbon source. Nat. Geosci. 3, 542–545. Publisher Full Text

Ramanathan, V., Crutzen, P., Kiehl, J. and Rosenfeld, D. 2001a. Aerosols, climate, and the hydrological cycle. Science. 294, 2119–2124. Publisher Full Text

Ramanathan, V., Crutzen, P. J., Lelieveld, J., Mitra, A., Althausen, D. and co-authors. 2001b. Indian Ocean experiment: an integrated analysis of the climate forcing and effects of the great Indo-Asian haze. J. Geophys. Res. 106, 28371–28398. Publisher Full Text

Rastogi, N., Singh, A., Sarin, M. M. and Singh, D. 2015. Temporal variability of primary and secondary aerosols over northern India: impact of biomass burning emissions. Atmos. Environ. 125, 396–403. Publisher Full Text

Rengarajan, R., Sarin, M. M. and Sudheer, A. K. 2007. Carbonaceous and inorganic species in atmospheric aerosols during wintertime over urban and high-altitude sites in North India. J. Geophys. Res. 112, D21307. Publisher Full Text

Sarin, M., Kumar, A., Srinivas, B., Sudheer, A. and Rastogi, N. 2010. Anthropogenic sulphate aerosols and large Cl-deficit in marine atmospheric boundary layer of tropical Bay of Bengal. J. Atmos. Chem. 66, 1–10. Publisher Full Text

Sciare, J., Oikonomou, K., Favez, O., Liakakou, E., Markaki, Z. and co-authors. 2008. Long-term measurements of carbonaceous aerosols in the Eastern Mediterranean: evidence of long-range transport of biomass burning. Atmos. Chem. Phys. 8, 5551–5563. Publisher Full Text

Smith, B. N. and Epstein, S. 1971. Two categories of 13C/12C ratios for higher plants. Plant Physiol. 47, 380–384. Publisher Full Text

Srinivas, B., Sarin, M. and Kumar, A. 2011a. Impact of anthropogenic sources on aerosol iron solubility over the Bay of Bengal and the Arabian Sea. Biogeochemistry. 110, 257–268. Publisher Full Text

Srinivas, B., Sarin, M. M. and Sarma, V. V. S. S. 2011b. Atmospheric dry deposition of inorganic and organic nitrogen to the Bay of Bengal: impact of continental outflow. Mar. Chem. 127, 170–179. Publisher Full Text

Srinivas, B. and Sarin, M. 2012. Atmospheric pathways of phosphorous to the Bay of Bengal: contribution from anthropogenic sources and mineral dust. Tellus B. 64, 17174. DOI:

Srinivas, B. and Sarin, M. 2013a. Light absorbing organic aerosols (brown carbon) over the tropical Indian Ocean: impact of biomass burning emissions. Environ. Res. Lett. 8, 044042. Publisher Full Text

Srinivas, B. and Sarin, M. M. 2013b. Atmospheric dry-deposition of mineral dust and anthropogenic trace metals to the Bay of Bengal. J. Mar. Syst. 126, 56–68. Publisher Full Text

Srinivas, B. and Sarin, M. M. 2014. PM2.5, EC and OC in atmospheric outflow from the Indo-Gangetic Plain: temporal variability and aerosol organic carbon-to-organic mass conversion factor. Sci. Total Environ. 487, 196–205. Publisher Full Text

Streets, D. G., Bond, T. C., Carmichael, G. R., Fernandes, S. D., Fu, Q. and co-authors. 2003. An inventory of gaseous and primary aerosol emissions in Asia in the year 2000. J. Geophys. Res. 108, 8809. Publisher Full Text

Sudheer, A. K. and Sarin, M. M. 2008. Carbonaceous aerosols in MABL of Bay of Bengal: influence of continental outflow. Atmos. Environ. 42, 4089–4100. Publisher Full Text

Turekian, V. C., Macko, S., Ballentine, D., Swap, R. J. and Garstang, M. 1998. Causes of bulk carbon and nitrogen isotopic fractionations in the products of vegetation burns: laboratory studies. Chem. Geol. 152, 181–192. Publisher Full Text

Turekian, V. C., Macko, S. A. and Keene, W. C. 2003. Concentrations, isotopic compositions, and sources of size-resolved, particulate organic carbon and oxalate in near-surface marine air at Bermuda during spring. J. Geophys. Res. 108, 4157. Publisher Full Text

Venkataraman, C., Reddy, C. K., Josson, S. and Reddy, M. S. 2002. Aerosol size and chemical characteristics at Mumbai, India, during the INDOEX-IFP (1999). Atmos. Environ. 36, 1979–1991. Publisher Full Text

Venkataraman, C., Habib, G., Eiguren-Fernandez, A., Miguel, A. and Friedlander, S. 2005. Residential biofuels in South Asia: carbonaceous aerosol emissions and climate impacts. Science. 307, 1454–1456. Publisher Full Text

Widory, D. 2006. Combustibles, fuels and their combustion products: a view through carbon isotopes. Combust. Theory Model. 10, 831–841. Publisher Full Text

Widory, D. 2007. Nitrogen isotopes: tracers of origin and processes affecting PM10 in the atmosphere of Paris. Atmos. Environ. 41, 2382–2390. Publisher Full Text


Now at: Chubu Institute for Advanced Studies, Chubu University, Kasugai, Japan

About The Authors

Srinivas Bikkina
Post-doctoral research fellow Dept of Applied Environmental Sciences Bolin Center for Climate Research Stokcholm University, Sweden E-mail:

Post-doctoral research fellow
Dept of Applied Environmental Sciences
Bolin Center for Climate Research
Stokcholm University, Sweden



Kimitaka Kawamura

Manmohan Sarin

Article Metrics

Metrics Loading ...

Metrics powered by PLOS ALM

Related Content