India Is Overtaking China as the World’s Largest Emitter of Anthropogenic Sulfur Dioxide

India Is Overtaking China as the World’s Largest Emitter of Anthropogenic Sulfur Dioxide www.nature.com/scientificreports Correction: Author Correction OPEN India Is Overtaking China as the World’s Largest Emitter of Anthropogenic Sulfur Dioxide Received: 7 July 2017 1,2 3 3 2 4 Can Li , Chris McLinden , Vitali Fioletov , Nickolay Krotkov , Simon Carn , Joanna 2 5 6 6,7 1,6,8 1,6 Accepted: 12 October 2017 Joiner , David Streets , Hao He , Xinrong Ren , Zhanqing Li & Russell R. Dickerson Published: xx xx xxxx Severe haze is a major public health concern in China and India. Both countries rely heavily on coal for energy, and sulfur dioxide (SO ) emitted from coal-fired power plants and industry is a major pollutant contributing to their air quality problems. Timely, accurate information on SO sources is a required input to air quality models for pollution prediction and mitigation. However, such information has been difficult to obtain for these two countries, as fast-paced changes in economy and environmental regulations have often led to unforeseen emission changes. Here we use satellite observations to show that China and India are on opposite trajectories for sulfurous pollution. Since 2007, emissions in China have declined by 75% while those in India have increased by 50%. With these changes, India is now surpassing China as the world’s largest emitter of anthropogenic SO . This finding, not predicted by emission scenarios, suggests effective SO control in China and lack thereof in India. Despite this, haze remains severe in China, indicating the importance of reducing emissions of other pollutants. In India, ~33 million people now live in areas with substantial SO pollution. Continued growth in emissions will adversely affect more people and further exacerbate morbidity and mortality. China and India are the top two consumers of coal in the world . Coal typically contains a few percent of sulfur by weight, and its combustion emits large amounts of SO , a toxic air pollutant. SO forms sulfate aerosols, the 2 2 principal component of the historic “London Smog” and a major contributor to the two countries’ current haze 2,3 4,5 problem that causes over one million premature deaths each year . Sulfate commonly makes up >10% of the 2 3 6 fine particles in China and India , oen m ft uch more during heavy pollution episodes . To predict and mitigate air pollution, air quality models require accurate information on the emissions of SO and other pollutants. In the conventional approach, one compiles bottom-up emission inventories based on activity rates and emission factors. These inventories are normally updated every 3–5 years and oen h ft ave to be projected for very recent years. Substantial uncertainties can therefore exist in the estimated or projected emissions, especially for regions 8 9 experiencing rapid changes in economy and environmental regulations such as China and India . Advances in satellite measurements have yielded new data and techniques that help to evaluate and improve 10–13 bottom-up inventories . For SO , the Ozone Monitoring Instrument (OMI) has been particularly useful owing to its superior ground resolution . OMI SO measurements uncovered the first evidence that China had started to reduce emissions through the installation of flue gas desulfurization devices , and also observed large changes 16,17 9 in SO emissions from power plants in the U.S. and India . More recently, a new technique that combines wind and improved SO data was employed to develop an OMI-based emission catalogue for nearly 500 sources around 18–21 the globe . This technique enabled the detection of ~40 sources missing from the conventional bottom-up 18 22 inventories and provided the first emission estimates for a number of degassing volcanoes in remote locations . 1 2 Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD, 20742, USA. Atmospheric Chemistry and Dynamics Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD, 20771, USA. Air Quality Research Division, Environment and Climate Change Canada, Toronto, M3H 5T4, Canada. Department of Geological and Mining Engineering and Sciences, Michigan Technological University, Houghton, MI, 49931, USA. Energy Systems Division, Argonne National Laboratory, Argonne, IL, 60439, USA. Department of Atmospheric and Oceanic Science, University of Maryland, College Park, MD, 20742, USA. Air Resources Laboratory, National Oceanic and Atmospheric Administration, College Park, MD, 20740, USA. State Key Laboratory of Earth Surface Processes and Resource Ecology and College of Global Change and Earth System Science, Beijing Normal University, Beijing, 100875, China. Correspondence and requests for materials should be addressed to C.L. (email: can.li@nasa.gov) SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 1 www.nature.com/scientificreports/ Figure 1. Changes in SO loading over India and China between 2005 and 2016. (a) Average SO vertical column 2 2 amounts over India and China for 2005 from the OMI instrument on the Aura satellite, expressed in Dobson Units 16 −2 (1 DU = 2.69 × 10 molecules cm ). (b) Same as (a) but for 2016, showing significant increase and decrease of SO over India and China, respectively, during the 12-year span. The maps in the figure were generated by Chris McLinden using Matlab (version 2016a; https://www.mathworks.com/products/matlab.html). Here we analyze OMI SO data to study the changes in SO pollution in China and India from 2005 to 2016. 2 2 We examine several recent emission projections to determine whether our observed changes were predicted in any emission scenarios. To investigate the underlying causes for the die ff rent trends between China and India, we compare emissions to coal consumption. Finally, we investigate the implications of these changes in SO pollution in terms of their health impacts. Results Changes in SO loading. For both China and India, OMI data show large differences in SO loading 2 2 between 2005 and 2016, and in Fig. 1a, one can identify isolated hot spots with SO column amount >0.5 Dobson 16 −2 Units (DU, 1 DU = 2.69 × 10 molecules cm ) over India in 2005. Several are associated with large coal-fired power plants in the northeastern Indian states of Odisha, Jharkhand, and Chhattisgarh, the southeastern state of Tamil Nadu (which includes Chennai), and the western state of Maharashtra (which includes Mumbai). By 2016 (Fig. 1b), these hotspots in northeastern India have grown into a cluster covering a large area, likely due to 9,23 emissions from new power plants constructed over the past decade . SO columns in other parts of the country have also increased, particularly near Jamnagar on the west coast, where expansion of a large oil refinery and construction of the largest power plant in India took place in 2008–2012. As for China, SO exceeded 0.5 DU over almost the entire North China Plain in 2005 (Fig. 1a). SO columns of 2 2 >2 DU are found over Hebei and Henan, two provinces just east of the Taihang Mountains and home to numer- ous power plants as well as coking and cement industries. Additional SO hotspots can be seen over the populous Sichuan Basin, the Pearl River Delta in southern China (which includes Guangzhou and Hong Kong), the Yangtze River Delta in eastern China (which includes Shanghai), as well as Inner Mongolia. By 2016, almost no hotspots with SO >0.5 DU can be found outside of the North China Plain. Even for areas near the Taihang Mountains, the SO loading has generally decreased to below 1 DU (Fig. 1b). Changes in SO emissions. We estimate the annual SO emissions from China and India during 2005–2016 2 2 (Tables S1 and S2, Supplementary Material), by first summing up the sources in the OMI catalogue for the two countries. e Th catalogue includes emissions estimated based on OMI data (see Methods) for 47 sources in India and 82 sources in China. One caveat is that OMI can only detect large point sources and globally, the catalogue represents approximately 50% of all emissions in bottom-up inventories . Here we compare the OMI-derived catalogue emissions with those from several recent regional inventories (Table 1). For China, the ratio between OMI catalogue emissions and inventories ranges within 40–62%; for India, this ratio is 36–48%. We then use the average of these ratios to adjust or normalize the OMI-derived emissions to reflect the national total. As shown in Fig. 2a, the normalized estimates reveal opposite trends in China and India. The 6 −1 emissions from China peaked at 36.6 Mt (10 tonnes) yr in 2007 and have since been on a generally decreasing −1 −1 trajectory. At 8.4 Mt yr , the level in 2016 is 26% of that in 2005 (31.8 Mt yr ). The decrease reflects stricter pollution control measures, coupled with a gradual shift to other, non-coal-based energy sources, and the recent slowdown of the Chinese economy. Since the early 2000s, the Chinese government has introduced, for example, 24 8 policies to reduce SO emissions and a new national air quality standard for fine particles . Electricity generation in China grew by more than 100% during 2005–2015, but coal consumption increased by ~50% . The brief period of emission growth in 2009–2011 can probably be attributed to government stimulus in response to the global financial crisis of 2007–2008. e es Th timated emissions for India, on the other hand, indicate relatively steady growth throughout the entire −1 −1 period. For 2016, the emissions from India (11.1 Mt yr , 9.5–12.6 Mt yr considering the range of OMI/ −1 bottom-up ratios) are at approximately the same level as China (7.5–11.6 Mt yr ). If the current trends continue, SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 2 www.nature.com/scientificreports/ * −1 ** −1 Emission Estimates (Mt yr ) Emission Projections (Mt yr ) 2005 2006 2007 2008 2009 2010 2011 2015 2020 2025 2030 Source 7.8 (48%) 10.4 12.4–12.9 2.9–15.3 3.3–18.8 37 6.4 (39%) 16.0 31.5 38 India 5.8 (43%) 7.5 (45%) 7.9 (47%) 9.1–9.7 9.3–11.9 8.6–13.5 8.7–15.6 34 6.8 (36%) 7.2 (40%) 7.7 (41%) 8.4 (38%) 9.1 (38%) 9.5 (39%) 10.1 (36%) 25 8.0 (40%) 8.8 (42%) 39 34.4 (51%) 33.3 32.9 38 30.4 (58%) 33.2 (40%) 33.8 (45%) 33.6–34.6 28.8–33.3 22.4–30.2 17.7–27.7 34 32.4 (54%) 33.3 (58%) 32.6 (62%) 31.2 (55%) 31.0 (43%) 29.8 (51%) 29.1 (60%) 25 32.1 (53%) 30.8 (49%) 39 China 28.6 (61%) 22.9–33.0 40 28.7 (61%) 24.4 (62%) 15.7–29.1 8.3–30.7 28 27.7 (55%) 19.6–33.8 19.6–36.3 16.6–37.8 15.5–38.1 41 32.3 (54%) 33.2 (58%) 32.3 (62%) 31.3 (55%) 24 31.0 (62%) 42 Table 1. Recent bottom-up estimates and projections of SO emissions for China and India. Percentages in parentheses are the fraction of bottom-up emissions observed by OMI. The fraction for China is 40–62%, ** with a mean of 55%. The fraction for India is 36–48%, with a mean of 41%. Lower end of the ranges typically represents emission scenarios with stricter emission control and energy policies that facilitate the shift to alternate energy sources rather than coal. Higher end typically represents emission scenarios with current environmental and energy policies (business as usual). Figure 2. Emissions, loading, and potential impact of SO in India and China. (a) Total annual SO emissions 2 2 for India and China during 2005–2016 based on OMI measurements. To account for the sources that are undetectable by OMI, the top-down estimates from the OMI emission catalogue are normalized using the average ratio (55% for China, 41% for India) between the catalogue and various bottom-up inventories in Table 1. The lower and upper bounds of the error bars are the same OMI estimates normalized with the maximum and minimum ratios between OMI and bottom-up inventories, respectively. Black lines represent annual emissions from one of the bottom-up inventories . Vertical bars show the range of projected emissions for 2015 (Table 1). (b) The ratio of the unnormalized OMI emission estimates to coal consumption during 2005–2015 (unit: tonne SO /tonne oil equivalent). (c) Population-weighted SO loading in Dobson Units. 2 2 (d) Population living in areas with annual mean SO of at least 0.5 DU during 2013–2016. SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 3 www.nature.com/scientificreports/ India will emit significantly more SO than China in the coming years. For both countries, the ratio between the OMI catalogue emissions and total emissions may change over time. We also estimate SO emissions based on the OMI-observed total SO mass, calculated from observations over the entire country and therefore less ae ff cted by the detection limit. We arrive at the same conclusion that India is becoming, if it is not already, the world’s top SO emitting country (Tables S3 and S4). It is enlightening to compare the OMI catalogue emissions with coal consumption (Fig. 2b). The ratio between the two is similar for China and India for 2005–2007 at ~0.012–0.013 tonne SO /tonne oil equivalent. Considering that OMI observes ~50% of all SO sources, that ~70–90% of SO emissions are from coal , and upon accounting 2 2 for standard conversions (1 tonne oil equivalent = 1.4 tonnes coal), we arrive at an effective emission factor of 24,26 12–16 g SO /kg coal. This is largely consistent with previously used values in bottom-up inventories and sug- gests little control on SO in either country before 2007. Since then, the ratio has remained virtually unchanged for India, implying continued absence of SO emission control . The ratio for China, however, has dropped to ~0.002–0.003 (corresponding emission factor: 2–3 g SO /kg coal), suggesting ee ff ctive control measures that have eliminated ~80% of potential SO emissions. Assuming that carbon makes up ~60–80% of coal by weight, our estimated emission factor for 2015 implies −4 a molar ratio of ~5–9 × 10 between SO and CO emitted from coal combustion in China. This is compara- 2 2 −4 ble with the recently measured ΔSO /ΔCO ratio of ~3–10 × 10 (ΔSO and ΔCO represent the observed 2 2 2 2 enhancements within plumes over background levels) in the boundary layer over Hebei in spring 2016 (Fig. S1), confirming the efficient SO removal in China. Satellite observations also point to a ~25% increase in NH over 2 3 China during our study period. This relatively modest growth (as compared with the decrease in SO ) is partially attributable to reductions in SO as a sink for NH . It also suggests that there must be excess NH and other alka- 2 3 3 line cations that neutralize sulfate; otherwise the growth rate in NH would have been much greater. Indeed recent measurements in northern China seem to indicate complete neutralization of sulfate and nitrate in aerosols. In Table 1, we examine projections of SO emissions from several studies published between 2009 and 2015. −1 −1 For India, the projected emissions for 2015 are 9.1–10.4 Mt yr , close to our estimate of 8.5–11.3 Mt yr −1 (Table S1). For China, the projected emissions for 2015 (19.6–33.8 Mt yr ) are a factor of 1.5–4 greater than our −1 estimate (8.7–13.5 Mt yr ). In fact, all but one study predicted that SO emissions from China would still exceed −1 28 −1 15 Mt yr even in 2030. In the only exception (8 Mt yr in 2030), it is assumed that lifestyle-changing energy policies and the most efficient emission control technology would be fully implemented in China. e Th difference between our observation and projections suggests that there are currently much more efficient SO controls in China than assumed in the various emission scenarios. Population exposure to SO pollution. Population-weighted SO loading (Fig. 2c, Table S5) closely fol- lows OMI-estimated emissions. Over the past 10 years, the SO loading over China decreased by a factor of v fi e, from 0.89 DU in 2007 to 0.17 DU in 2016. At the same time, the loading over India climbed by nearly 50%, reaching 0.13 DU in 2016. There is no simple relationship between the OMI-observed column amount and the concentration at ground level. If we assume that all the SO is within the lowest 1000 m of the atmosphere and well mixed at 1:30 pm local time (OMI overpass time), an SO column of 0.5 DU corresponds to a mass concentration −3 −3 of ~14 μg m . Given that the World Health Organization’s guideline for SO is 20 μg m (for a 24-hour mean), column amounts of 0.5–1 DU represent sufficiently high SO loading to adversely ae ff ct human health both as a toxic gas and as a precursor to sulfate aerosols. In China, over 450 million people were exposed to >0.5 DU of SO in 2013, but this number decreased to 99 million in 2016 (Fig. 2d). Similarly, the population exposed to >1.0 DU of SO decreased from ~190 million in 2013 to 13 million in 2016, a remarkable drop of over 90% (Table S6). As for India, 13 (0.7) million people were exposed to >0.5 (1.0) DU of SO in 2013. In just three years, this has grown to 33 (3.8) million people (Table S7). Discussion Our findings have important implications for future environmental policies in both countries. Despite the large reductions in SO , haze in China remains a severe environmental issue . This may be partly due to the shift in the thermodynamic equilibrium of the sulfate-nitrate-ammonium system . It will be critical to better understand the benefits of SO reductions before viable and balanced policies can be devised to further improve air quality in China. To a certain extent, the impact of SO emissions is presently limited in India, as SO loading is relatively 2 2 low over the densely populated Indo-Gangetic Plain. But this may change as the demand for electricity continues to grow. In the various Representative Concentration Pathways for the latest Assessment Report (AR5) by the Intergovernmental Panel on Climate Change , SO emissions from Asia were projected to increase until the 2020s before starting to decrease. The sooner-than-expected reductions in SO could also accelerate regional warming, as they would reduce the loading of sulfate aerosols that scatter sunlight and partially offset the warm- ing effects of greenhouse gases. Methods OMI SO data. SO data used in this study were retrieved from earthshine radiances in the wavelength range 2 2 of 310.5–340 nm measured by the Ozone Monitoring Instrument (OMI) aboard the NASA Aura spacecraft. 16 −2 e r Th esults are in Dobson Units (1 DU = 2.69 × 10 molecules cm ), and represent the estimated total number of SO molecules in the entire atmospheric column above a unit area (or simply, column amount). The current retrieval algorithm applies a principal component analysis technique to OMI radiances to minimize spectral interferences and maximize the quality of SO data. A detailed description of the retrieval technique can be found 20,32 elsewhere . Because the OMI SO sensitivity varies with altitude, the retrieved total column amount depends on the assumed vertical distribution of SO . Several different fixed SO profiles are assumed in operational OMI 2 2 SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 4 www.nature.com/scientificreports/ retrievals. The present study uses version 1.3 level 2 (orbital level) OMI retrievals assuming that all SO is in the planetary boundary layer (PBL, or the lowest 1 km of the atmosphere). For the present study, OMI pixels with a radiative cloud fraction >0.2 or solar zenith angle >70° were excluded from data analysis. Data from OMI cross-track positions (or rows) ae ff cted by the row anomaly (http:// projects.knmi.nl/omi/research/product/rowanomaly-background.php) or near the edge of the swath (rows 1–10 and 51–60) were also excluded. Additionally, data from days potentially influenced by large transient volcanic plumes were excluded. Details on the data filtering can be found elsewhere . The SO columns from the remain- ing OMI pixels were then averaged to a spatial resolution of 0.1 × 0.1° for the maps in Fig. 1. OMI-based SO source detection and estimate. e Th methods for source detection and emission esti- 18,19,21 mate are based on a previously described algorithm that combines satellite measurements with reanalysis wind data (ECMWF interim reanalysis was used here). Wind information is matched with each OMI pixel. Emissions from about 500 continuously emitting point sources (or clusters of sources in close proximity), includ- ing 47 in India and 82 in China, are derived from OMI and wind data by tracking the downwind decay of the −1 −1 plumes. These sources have estimated SO emissions ranging from about 30 kt yr to more than 4000 kt yr . Due to the coarse spatial resolution of OMI (relative to a point source) and the limited precision of individual SO column observations, data spanning a year are analyzed together using a wind rotation scheme to align the wind 19 19 vectors of all overpasses considered . The emissions were estimated by fitting OMI columns to a plume function 18,21 consisting of coordinates, wind speed and direction, with a single parameter representing the total mass . Other fitting parameters, including an effective lifetime (5.9 hours), are specified . Emissions are then calculated as the ratio of mass to lifetime, effectively assuming a steady-state. The operational OMI retrievals use an effective air mass factor of 0.36 for all locations. In the emission estimate algorithm, OMI data for each emission source were adjusted using an air mass factor calculated for the location based on its elevation, surface albedo and sun/ 18,21 viewing geometry to better represent OMI sensitivities to the local source . Wet removal of pollutants by the summer monsoon rainfall causes a strong seasonality in air pollution in India, especially for aerosols . OMI generally also observes smaller SO columns over India during summer months (see monthly maps at https://avdc.gsfc.nasa.gov/pub/data/satellite/Aura/OMI/V03/L3/OMSO2m/ Monthly_mean_jpeg/). This seasonality may be partially attributed to the washout effects and the shorter lifetime of SO in summer, but may also reflect reduced coverage by OMI due to increased cloud cover during the mon - soon. For OMI-based emission estimates, since pixels with a small cloud fraction from an entire year are analyzed, the collection of data used to derive emissions for a given source may be more representative of non-summer conditions. The impact of this seasonal change in sampling on estimated emissions is currently unclear, but it is unlikely to significantly ae ff ct their long-term trend. Aircraft measurements of ΔSO /ΔCO ratio. Between 8 May and 11 June 2016, a twin engine Y-12 2 2 research aircraft was flown on 11 missions over the heavily industrialized Hebei Province of China. A modified, commercial pulsed-fluorescence detector (TEI Model 43 C) was used to measure ambient SO . A Picarro cavity ring-down spectrometer (Model G2401-m) was used to measure CO . Profiles were flown from near the surface to the top of the planetary boundary layer, at ~1500 m above ground. The Δ SO /ΔCO ratio was determined 2 2 from the deviation from background in plumes. Only data with significant correlation between ∆SO and ∆CO 2 2 (R > 0.6) are included in Fig. S1. Population data and population exposure to SO pollution. Population data for 2005, 2010, and 2015 from the Gridded Population of the World, Version 4 (GPWv4) were used in this study. For each of the three years, the (30 arc seconds) GPWv4 population count and nation identifier data were used to calculate the counts of Chinese and Indian population for each grid cell in Fig. 1. An annual growth rate was then estimated for each grid cell between 2005 and 2010, and between 2010 and 2015, to interpolate population data to other years. With OMI SO (Ω) and population count (P) data now on the same grid, the population-weighted SO column 2 2 amount (Ω ) for the entire domain with n grid cells can be calculated as: in = PΩ ii i=1 Ω= , in = ∑ (1) i=1 th where P and Ω are population count and OMI SO column for the i grid cell, respectively. i i 2 Data Availability. Level 2 Principal Component Analysis SO data from OMI are available from the Goddard Earth Science Data and Information Service Center (http://disc.sci.gsfc.nasa.gov/). Wind reanalysis data are avail- able from ECMWF (http://apps.ecmwf.int/datasets/data/interim-full-daily). Derived SO emissions are available from the global SO monitoring website at NASA Goddard Space Flight Center (https://so2.gsfc.nasa.gov). The GPWv4 population data are available from the Socioeconomic Data and Applications Center (SEDAC) in NASA’s Earth Observing System Data and Information System and hosted by Center for International Earth Science Information Network at Columbia University (http://sedac.ciesin.columbia.edu/data/collection/gpw-v4). Aircra ft measurements acquired during the ARIAS campaign are available upon request from X. Ren (ren@umd.edu). References 1. BP P.L.C. BP Statistical Review of WorldEnergy 2016. London, UK: BP p.l.c. https://www.bp.com/content/dam/bp/pdf/energy- economics/statistical-review-2016/bp-statistical-review-of-world-energy-2016-full-report.pdf. Accessed 27 Mar. 2017 (2016). 2. Huang, R.-J. et al. High secondary aerosol contribution to particulate pollution during haze events in China. Nature 514, 218–222 (2014). SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 5 www.nature.com/scientificreports/ 3. Chakraborty, A., Gupta, T. & Tripathi, S. N. Chemical composition and characteristics of ambient aerosols and rainwater residues during Indian summer monsoon: Insight from aerosol mass spectrometry. Atmos. Environ. 136, 144–155 (2016). 4. Lelieveld, J., Evans., J. S., Fnais, M., Giannadaki, D. & Pozzer, A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature 525, 367–371 (2015). 5. Zhang, Q. et al. Transboundary health impacts of transported global air pollution and international trade. Nature 543, 705–709 (2017). 6. Wang, G. et al. Persistent sulfate formation from London Fog to Chinese haze. Proc. Natl. Acad. Sci. 113, 13630–13635 (2016). 7. Environmental Protection Agency. National Emissions Inventory. Research Triangle Park, NC: Environmental Protection Agency, https://www.epa.gov/air-emissions-inventories/national-emissions-inventory-nei. Accessed 15 Mar. 2017 (2014). 8. Zhang, Q., He, K. B. & Huo, H. Cleaning China’s air. Nature 484, 161–162 (2012). 9. Lu, Z., Streets, D. G., de Foy, B. & Krotkov, N. A. Ozone monitoring instrument observations of interannual increases in SO emissions from Indian coal-fired power plants during 2005-2012. Environ. Sci. Technol. 47, 13993–14000 (2013). 10. Streets, D. G. et al. Emissions estimation from satellite retrievals: A review of current capability. Atmos. Environ. 77, 1011–1042 (2013). 11. Martin, R. V. Satellite remote sensing of surface air quality. Atmos. Environ. 42, 7823–7843 (2008). 12. Beirle, S. et al. Megacity emissions and lifetimes of nitrogen oxides probed from space. Science 333, 1737–1739 (2011). 13. Lamsal, L. N. et al. Application of satellite observations for timely updates to global anthropogenic NOx emission inventories. Geophys. Res. Lett. 38, L05810 (2011). 14. Fioletov, V. E. et al. Application of OMI, SCIAMACHY, and GOME-2 satellite SO retrievals for detection of large emission sources. J. Geophys. Res. Atmos. 118, 11399–11418 (2013). 15. Li, C. et al. Recent large reduction in sulfur dioxide emissions from Chinese power plants observed by the Ozone Monitoring Instrument. Geophys. Res. Lett. 37, L08807 (2010). 16. Fioletov, V. E., McLinden, C. A., Krotkov, N., Moran, M. D. & Yang, K. Estimation of SO emissions using OMI retrievals. Geophys. Res. Lett. 38, L21811 (2011). 17. He, H. et al. Response of SO and particulate air pollution to local and regional emission controls: A case study in Maryland. Earth’s Future 4, 94–109 (2016). 18. McLinden, C. A. et al. Space-based detection of missing sulfur dioxide sources of global air pollution. Nature Geoscience 9, 496–500 (2016). 19. Fioletov, V. E., McLinden, C. A., Krotkov, N. & Li, C. Lifetimes and emissions of SO from point sources estimated from OMI. Geophys. Res. Lett. 42, 1969–1976 (2015). 20. Li, C., Joiner, J., Krotkov, N. A. & Bhartia, P. K. A fast and sensitive new satellite SO retrieval algorithm based on principal component analysis: application to the Ozone Monitoring Instrument. Geophys. Res. Lett. 40, 6314–6318 (2013). 21. Fioletov, V. E. et al. A global catalogue of large SO sources and emissions derived from the Ozone Monitoring Instrument. Atmos. Chem. Phys. 16, 11497–11519 (2016). 22. Carn, S. A., Fioletov, V. E., McLinden, C. A., Li, C. & Krotkov, N. A. A decade of global volcanic SO emissions measured from space. Scientific Reports 7, 44095 (2017). 23. Krotkov, N. A. et al. Aura OMI observations of regional SO and NO pollution changes from 2005 to 2015. Atmos. Chem. Phys. 16, 2 2 4605–4629 (2016). 24. Lu, Z. et al. Sulfur dioxide emissions in China and sulfur trends in East Asia since 2000. Atmos. Chem. Phys. 10, 6311–6331 (2010). 25. Klimont, Z., Smith, S. J. & Cofala, J. The last decade of global anthropogenic sulfur dioxide: 2000–2011 emissions. Environ. Res. Lett. 8, 014003 (2013). 26. Shekar Reddy, M. & Venkataraman, C. Inventory of aerosol and sulphur dioxide emissions from India: I—Fossil fuel combustion. Atmos. Environ. 36, 677–697 (2002). 27. Warner, J. X. et al. Increased atmospheric ammonia over the world’s major agricultural areas detected from space. Geophys. Res. Lett. 44, https://doi.org/10.1002/2016GL072305 (2017). 28. Wang, S. X. et al. Emission trends and mitigation options for air pollutants in East Asia. Atmos. Chem. Phys. 14, 6571–6603 (2014). 29. Song, C. et al. Health burden attributable to ambient PM in China. Environ. Pollution 223, 575–586 (2017). 2.5 30. Stocker, T. F et al. (eds)Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge Univ. Press, 2013). 31. Levelt, P. F. et al. The Ozone Monitoring Instrument. IEEE Trans. Geosci. Remote 44, 1093–1101 (2006). 32. Li, C. et al. New-generation NASA aura ozone monitoring instrument (OMI) volcanic SO dataset: Algorithm description, initial results, and continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS). Atmos. Meas. Tech. 10, 445–458 (2017). 33. Dee, D. P. et al. e Th ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 137, 553–597 (2011). 34. International Institute for Applied Systems Analysis. Emissions of Air Pollutants for the World Energy Outlook 2011 Energy Scenarios, Laxenburg, Austria: International Institute for Applied Systems Analysis, http://pure.iiasa.ac.at/9766/1/XO-11-028.pdf. Accessed 27 Mar. 2017 (2011). 35. Hyvärinen, A.-P. et al. Effect of the summer monsoon on aerosols at two measurement stations in Northern India – Part 1: PM and BC concentrations. Atmos. Chem. Phys. 11, 8271–8282 (2011). 36. Center for International Earth Science Information Network. Gridded Population of the World, Version 4 (GPWv4): Population Count. Palisades, NY: NASA Socioeconomic Data and Applications Center. http://dx.doi.org/10.7927/H4X63JVC. Accessed 16 Feb. 2017 (2016). 37. International Institute for Applied Systems Analysis. Implications of energy trajectories from the World Energy Outlook 2015 for India’s air pollution, Laxenburg, Austria: International Institute for Applied Systems Analysis, http://www.worldenergyoutlook.org/media/ weowebsite/2015/Air_pollution_emissions_impacts_India_WEO2015_IIASA.pdf. Accessed 27 Mar. 2017 (2015). 38. Klimont, Z. et al. Projections of SO , NO and carbonaceous aerosols emissions in Asia. Tellus B 61, 602–617 (2009). 2 x 39. Lu, Z., Zhang, Q. & Streets, D. G. Sulfur dioxide and primary carbonaceous aerosol emissions in China and India, 1996–2010. Atmos. Chem. Phys. 11, 9839–9864 (2011). 40. Xing, J. et al. Projections of air pollutant emissions and its impacts on regional air quality in China in 2020. Atmos. Chem. Phys. 11, 3119–3136 (2011). 41. Zhao, Y., Zhang, J. & Nielsen, C. P. The effects of energy paths and emission controls and standards on future trends in China’s emissions of primary air pollutants. Atmos. Chem. Phys. 14, 8849–8868 (2014). 42. Zhang, Q. et al. Asian emissions in 2006 for the NASA INTEX-B mission. Atmos. Chem. Phys. 9, 5131–5153 (2009). Acknowledgements We thank the NASA Earth Science Division (ESD) Aura Science Team program for funding of OMI SO2 product development and analysis (Grant # 80NSSC17K0240). The Dutch and Finnish built OMI instrument is part of the NASA’s Earth Observing System (EOS) Aura satellite payload. The OMI project is managed by the Royal Meteorological Institute of the Netherlands (KNMI) and the Netherlands Space Agency (NSO). We thank Dr. SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 6 www.nature.com/scientificreports/ Steve Fetter (University of Maryland) for helpful comments. We thank the National Science Foundations of China (91544217) and US (Grant # 1558259) for supporting aircraft measurements over China presented in this study. C.L. acknowledges partial support from NASA’s Earth Science New Investigator Program (Grant # NNX14AI02G). Author Contributions C.L., C.M., V.F. and S.C. planned the research. C.L., N.K. and J.J. provided satellite retrievals. C.L., C.M., V.F. and D.S. analyzed the satellite data. H.H., X.R., Z.L. and R.D. provided and analyzed aircraft data. C.L. wrote the paper. All authors commented on the manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-14639-8. Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 7 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

India Is Overtaking China as the World’s Largest Emitter of Anthropogenic Sulfur Dioxide

Free
7 pages

Loading next page...
 
/lp/springer_journal/india-is-overtaking-china-as-the-world-s-largest-emitter-of-JHOnTua832
Publisher
Springer Journals
Copyright
Copyright © 2017 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
D.O.I.
10.1038/s41598-017-14639-8
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports Correction: Author Correction OPEN India Is Overtaking China as the World’s Largest Emitter of Anthropogenic Sulfur Dioxide Received: 7 July 2017 1,2 3 3 2 4 Can Li , Chris McLinden , Vitali Fioletov , Nickolay Krotkov , Simon Carn , Joanna 2 5 6 6,7 1,6,8 1,6 Accepted: 12 October 2017 Joiner , David Streets , Hao He , Xinrong Ren , Zhanqing Li & Russell R. Dickerson Published: xx xx xxxx Severe haze is a major public health concern in China and India. Both countries rely heavily on coal for energy, and sulfur dioxide (SO ) emitted from coal-fired power plants and industry is a major pollutant contributing to their air quality problems. Timely, accurate information on SO sources is a required input to air quality models for pollution prediction and mitigation. However, such information has been difficult to obtain for these two countries, as fast-paced changes in economy and environmental regulations have often led to unforeseen emission changes. Here we use satellite observations to show that China and India are on opposite trajectories for sulfurous pollution. Since 2007, emissions in China have declined by 75% while those in India have increased by 50%. With these changes, India is now surpassing China as the world’s largest emitter of anthropogenic SO . This finding, not predicted by emission scenarios, suggests effective SO control in China and lack thereof in India. Despite this, haze remains severe in China, indicating the importance of reducing emissions of other pollutants. In India, ~33 million people now live in areas with substantial SO pollution. Continued growth in emissions will adversely affect more people and further exacerbate morbidity and mortality. China and India are the top two consumers of coal in the world . Coal typically contains a few percent of sulfur by weight, and its combustion emits large amounts of SO , a toxic air pollutant. SO forms sulfate aerosols, the 2 2 principal component of the historic “London Smog” and a major contributor to the two countries’ current haze 2,3 4,5 problem that causes over one million premature deaths each year . Sulfate commonly makes up >10% of the 2 3 6 fine particles in China and India , oen m ft uch more during heavy pollution episodes . To predict and mitigate air pollution, air quality models require accurate information on the emissions of SO and other pollutants. In the conventional approach, one compiles bottom-up emission inventories based on activity rates and emission factors. These inventories are normally updated every 3–5 years and oen h ft ave to be projected for very recent years. Substantial uncertainties can therefore exist in the estimated or projected emissions, especially for regions 8 9 experiencing rapid changes in economy and environmental regulations such as China and India . Advances in satellite measurements have yielded new data and techniques that help to evaluate and improve 10–13 bottom-up inventories . For SO , the Ozone Monitoring Instrument (OMI) has been particularly useful owing to its superior ground resolution . OMI SO measurements uncovered the first evidence that China had started to reduce emissions through the installation of flue gas desulfurization devices , and also observed large changes 16,17 9 in SO emissions from power plants in the U.S. and India . More recently, a new technique that combines wind and improved SO data was employed to develop an OMI-based emission catalogue for nearly 500 sources around 18–21 the globe . This technique enabled the detection of ~40 sources missing from the conventional bottom-up 18 22 inventories and provided the first emission estimates for a number of degassing volcanoes in remote locations . 1 2 Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD, 20742, USA. Atmospheric Chemistry and Dynamics Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD, 20771, USA. Air Quality Research Division, Environment and Climate Change Canada, Toronto, M3H 5T4, Canada. Department of Geological and Mining Engineering and Sciences, Michigan Technological University, Houghton, MI, 49931, USA. Energy Systems Division, Argonne National Laboratory, Argonne, IL, 60439, USA. Department of Atmospheric and Oceanic Science, University of Maryland, College Park, MD, 20742, USA. Air Resources Laboratory, National Oceanic and Atmospheric Administration, College Park, MD, 20740, USA. State Key Laboratory of Earth Surface Processes and Resource Ecology and College of Global Change and Earth System Science, Beijing Normal University, Beijing, 100875, China. Correspondence and requests for materials should be addressed to C.L. (email: can.li@nasa.gov) SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 1 www.nature.com/scientificreports/ Figure 1. Changes in SO loading over India and China between 2005 and 2016. (a) Average SO vertical column 2 2 amounts over India and China for 2005 from the OMI instrument on the Aura satellite, expressed in Dobson Units 16 −2 (1 DU = 2.69 × 10 molecules cm ). (b) Same as (a) but for 2016, showing significant increase and decrease of SO over India and China, respectively, during the 12-year span. The maps in the figure were generated by Chris McLinden using Matlab (version 2016a; https://www.mathworks.com/products/matlab.html). Here we analyze OMI SO data to study the changes in SO pollution in China and India from 2005 to 2016. 2 2 We examine several recent emission projections to determine whether our observed changes were predicted in any emission scenarios. To investigate the underlying causes for the die ff rent trends between China and India, we compare emissions to coal consumption. Finally, we investigate the implications of these changes in SO pollution in terms of their health impacts. Results Changes in SO loading. For both China and India, OMI data show large differences in SO loading 2 2 between 2005 and 2016, and in Fig. 1a, one can identify isolated hot spots with SO column amount >0.5 Dobson 16 −2 Units (DU, 1 DU = 2.69 × 10 molecules cm ) over India in 2005. Several are associated with large coal-fired power plants in the northeastern Indian states of Odisha, Jharkhand, and Chhattisgarh, the southeastern state of Tamil Nadu (which includes Chennai), and the western state of Maharashtra (which includes Mumbai). By 2016 (Fig. 1b), these hotspots in northeastern India have grown into a cluster covering a large area, likely due to 9,23 emissions from new power plants constructed over the past decade . SO columns in other parts of the country have also increased, particularly near Jamnagar on the west coast, where expansion of a large oil refinery and construction of the largest power plant in India took place in 2008–2012. As for China, SO exceeded 0.5 DU over almost the entire North China Plain in 2005 (Fig. 1a). SO columns of 2 2 >2 DU are found over Hebei and Henan, two provinces just east of the Taihang Mountains and home to numer- ous power plants as well as coking and cement industries. Additional SO hotspots can be seen over the populous Sichuan Basin, the Pearl River Delta in southern China (which includes Guangzhou and Hong Kong), the Yangtze River Delta in eastern China (which includes Shanghai), as well as Inner Mongolia. By 2016, almost no hotspots with SO >0.5 DU can be found outside of the North China Plain. Even for areas near the Taihang Mountains, the SO loading has generally decreased to below 1 DU (Fig. 1b). Changes in SO emissions. We estimate the annual SO emissions from China and India during 2005–2016 2 2 (Tables S1 and S2, Supplementary Material), by first summing up the sources in the OMI catalogue for the two countries. e Th catalogue includes emissions estimated based on OMI data (see Methods) for 47 sources in India and 82 sources in China. One caveat is that OMI can only detect large point sources and globally, the catalogue represents approximately 50% of all emissions in bottom-up inventories . Here we compare the OMI-derived catalogue emissions with those from several recent regional inventories (Table 1). For China, the ratio between OMI catalogue emissions and inventories ranges within 40–62%; for India, this ratio is 36–48%. We then use the average of these ratios to adjust or normalize the OMI-derived emissions to reflect the national total. As shown in Fig. 2a, the normalized estimates reveal opposite trends in China and India. The 6 −1 emissions from China peaked at 36.6 Mt (10 tonnes) yr in 2007 and have since been on a generally decreasing −1 −1 trajectory. At 8.4 Mt yr , the level in 2016 is 26% of that in 2005 (31.8 Mt yr ). The decrease reflects stricter pollution control measures, coupled with a gradual shift to other, non-coal-based energy sources, and the recent slowdown of the Chinese economy. Since the early 2000s, the Chinese government has introduced, for example, 24 8 policies to reduce SO emissions and a new national air quality standard for fine particles . Electricity generation in China grew by more than 100% during 2005–2015, but coal consumption increased by ~50% . The brief period of emission growth in 2009–2011 can probably be attributed to government stimulus in response to the global financial crisis of 2007–2008. e es Th timated emissions for India, on the other hand, indicate relatively steady growth throughout the entire −1 −1 period. For 2016, the emissions from India (11.1 Mt yr , 9.5–12.6 Mt yr considering the range of OMI/ −1 bottom-up ratios) are at approximately the same level as China (7.5–11.6 Mt yr ). If the current trends continue, SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 2 www.nature.com/scientificreports/ * −1 ** −1 Emission Estimates (Mt yr ) Emission Projections (Mt yr ) 2005 2006 2007 2008 2009 2010 2011 2015 2020 2025 2030 Source 7.8 (48%) 10.4 12.4–12.9 2.9–15.3 3.3–18.8 37 6.4 (39%) 16.0 31.5 38 India 5.8 (43%) 7.5 (45%) 7.9 (47%) 9.1–9.7 9.3–11.9 8.6–13.5 8.7–15.6 34 6.8 (36%) 7.2 (40%) 7.7 (41%) 8.4 (38%) 9.1 (38%) 9.5 (39%) 10.1 (36%) 25 8.0 (40%) 8.8 (42%) 39 34.4 (51%) 33.3 32.9 38 30.4 (58%) 33.2 (40%) 33.8 (45%) 33.6–34.6 28.8–33.3 22.4–30.2 17.7–27.7 34 32.4 (54%) 33.3 (58%) 32.6 (62%) 31.2 (55%) 31.0 (43%) 29.8 (51%) 29.1 (60%) 25 32.1 (53%) 30.8 (49%) 39 China 28.6 (61%) 22.9–33.0 40 28.7 (61%) 24.4 (62%) 15.7–29.1 8.3–30.7 28 27.7 (55%) 19.6–33.8 19.6–36.3 16.6–37.8 15.5–38.1 41 32.3 (54%) 33.2 (58%) 32.3 (62%) 31.3 (55%) 24 31.0 (62%) 42 Table 1. Recent bottom-up estimates and projections of SO emissions for China and India. Percentages in parentheses are the fraction of bottom-up emissions observed by OMI. The fraction for China is 40–62%, ** with a mean of 55%. The fraction for India is 36–48%, with a mean of 41%. Lower end of the ranges typically represents emission scenarios with stricter emission control and energy policies that facilitate the shift to alternate energy sources rather than coal. Higher end typically represents emission scenarios with current environmental and energy policies (business as usual). Figure 2. Emissions, loading, and potential impact of SO in India and China. (a) Total annual SO emissions 2 2 for India and China during 2005–2016 based on OMI measurements. To account for the sources that are undetectable by OMI, the top-down estimates from the OMI emission catalogue are normalized using the average ratio (55% for China, 41% for India) between the catalogue and various bottom-up inventories in Table 1. The lower and upper bounds of the error bars are the same OMI estimates normalized with the maximum and minimum ratios between OMI and bottom-up inventories, respectively. Black lines represent annual emissions from one of the bottom-up inventories . Vertical bars show the range of projected emissions for 2015 (Table 1). (b) The ratio of the unnormalized OMI emission estimates to coal consumption during 2005–2015 (unit: tonne SO /tonne oil equivalent). (c) Population-weighted SO loading in Dobson Units. 2 2 (d) Population living in areas with annual mean SO of at least 0.5 DU during 2013–2016. SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 3 www.nature.com/scientificreports/ India will emit significantly more SO than China in the coming years. For both countries, the ratio between the OMI catalogue emissions and total emissions may change over time. We also estimate SO emissions based on the OMI-observed total SO mass, calculated from observations over the entire country and therefore less ae ff cted by the detection limit. We arrive at the same conclusion that India is becoming, if it is not already, the world’s top SO emitting country (Tables S3 and S4). It is enlightening to compare the OMI catalogue emissions with coal consumption (Fig. 2b). The ratio between the two is similar for China and India for 2005–2007 at ~0.012–0.013 tonne SO /tonne oil equivalent. Considering that OMI observes ~50% of all SO sources, that ~70–90% of SO emissions are from coal , and upon accounting 2 2 for standard conversions (1 tonne oil equivalent = 1.4 tonnes coal), we arrive at an effective emission factor of 24,26 12–16 g SO /kg coal. This is largely consistent with previously used values in bottom-up inventories and sug- gests little control on SO in either country before 2007. Since then, the ratio has remained virtually unchanged for India, implying continued absence of SO emission control . The ratio for China, however, has dropped to ~0.002–0.003 (corresponding emission factor: 2–3 g SO /kg coal), suggesting ee ff ctive control measures that have eliminated ~80% of potential SO emissions. Assuming that carbon makes up ~60–80% of coal by weight, our estimated emission factor for 2015 implies −4 a molar ratio of ~5–9 × 10 between SO and CO emitted from coal combustion in China. This is compara- 2 2 −4 ble with the recently measured ΔSO /ΔCO ratio of ~3–10 × 10 (ΔSO and ΔCO represent the observed 2 2 2 2 enhancements within plumes over background levels) in the boundary layer over Hebei in spring 2016 (Fig. S1), confirming the efficient SO removal in China. Satellite observations also point to a ~25% increase in NH over 2 3 China during our study period. This relatively modest growth (as compared with the decrease in SO ) is partially attributable to reductions in SO as a sink for NH . It also suggests that there must be excess NH and other alka- 2 3 3 line cations that neutralize sulfate; otherwise the growth rate in NH would have been much greater. Indeed recent measurements in northern China seem to indicate complete neutralization of sulfate and nitrate in aerosols. In Table 1, we examine projections of SO emissions from several studies published between 2009 and 2015. −1 −1 For India, the projected emissions for 2015 are 9.1–10.4 Mt yr , close to our estimate of 8.5–11.3 Mt yr −1 (Table S1). For China, the projected emissions for 2015 (19.6–33.8 Mt yr ) are a factor of 1.5–4 greater than our −1 estimate (8.7–13.5 Mt yr ). In fact, all but one study predicted that SO emissions from China would still exceed −1 28 −1 15 Mt yr even in 2030. In the only exception (8 Mt yr in 2030), it is assumed that lifestyle-changing energy policies and the most efficient emission control technology would be fully implemented in China. e Th difference between our observation and projections suggests that there are currently much more efficient SO controls in China than assumed in the various emission scenarios. Population exposure to SO pollution. Population-weighted SO loading (Fig. 2c, Table S5) closely fol- lows OMI-estimated emissions. Over the past 10 years, the SO loading over China decreased by a factor of v fi e, from 0.89 DU in 2007 to 0.17 DU in 2016. At the same time, the loading over India climbed by nearly 50%, reaching 0.13 DU in 2016. There is no simple relationship between the OMI-observed column amount and the concentration at ground level. If we assume that all the SO is within the lowest 1000 m of the atmosphere and well mixed at 1:30 pm local time (OMI overpass time), an SO column of 0.5 DU corresponds to a mass concentration −3 −3 of ~14 μg m . Given that the World Health Organization’s guideline for SO is 20 μg m (for a 24-hour mean), column amounts of 0.5–1 DU represent sufficiently high SO loading to adversely ae ff ct human health both as a toxic gas and as a precursor to sulfate aerosols. In China, over 450 million people were exposed to >0.5 DU of SO in 2013, but this number decreased to 99 million in 2016 (Fig. 2d). Similarly, the population exposed to >1.0 DU of SO decreased from ~190 million in 2013 to 13 million in 2016, a remarkable drop of over 90% (Table S6). As for India, 13 (0.7) million people were exposed to >0.5 (1.0) DU of SO in 2013. In just three years, this has grown to 33 (3.8) million people (Table S7). Discussion Our findings have important implications for future environmental policies in both countries. Despite the large reductions in SO , haze in China remains a severe environmental issue . This may be partly due to the shift in the thermodynamic equilibrium of the sulfate-nitrate-ammonium system . It will be critical to better understand the benefits of SO reductions before viable and balanced policies can be devised to further improve air quality in China. To a certain extent, the impact of SO emissions is presently limited in India, as SO loading is relatively 2 2 low over the densely populated Indo-Gangetic Plain. But this may change as the demand for electricity continues to grow. In the various Representative Concentration Pathways for the latest Assessment Report (AR5) by the Intergovernmental Panel on Climate Change , SO emissions from Asia were projected to increase until the 2020s before starting to decrease. The sooner-than-expected reductions in SO could also accelerate regional warming, as they would reduce the loading of sulfate aerosols that scatter sunlight and partially offset the warm- ing effects of greenhouse gases. Methods OMI SO data. SO data used in this study were retrieved from earthshine radiances in the wavelength range 2 2 of 310.5–340 nm measured by the Ozone Monitoring Instrument (OMI) aboard the NASA Aura spacecraft. 16 −2 e r Th esults are in Dobson Units (1 DU = 2.69 × 10 molecules cm ), and represent the estimated total number of SO molecules in the entire atmospheric column above a unit area (or simply, column amount). The current retrieval algorithm applies a principal component analysis technique to OMI radiances to minimize spectral interferences and maximize the quality of SO data. A detailed description of the retrieval technique can be found 20,32 elsewhere . Because the OMI SO sensitivity varies with altitude, the retrieved total column amount depends on the assumed vertical distribution of SO . Several different fixed SO profiles are assumed in operational OMI 2 2 SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 4 www.nature.com/scientificreports/ retrievals. The present study uses version 1.3 level 2 (orbital level) OMI retrievals assuming that all SO is in the planetary boundary layer (PBL, or the lowest 1 km of the atmosphere). For the present study, OMI pixels with a radiative cloud fraction >0.2 or solar zenith angle >70° were excluded from data analysis. Data from OMI cross-track positions (or rows) ae ff cted by the row anomaly (http:// projects.knmi.nl/omi/research/product/rowanomaly-background.php) or near the edge of the swath (rows 1–10 and 51–60) were also excluded. Additionally, data from days potentially influenced by large transient volcanic plumes were excluded. Details on the data filtering can be found elsewhere . The SO columns from the remain- ing OMI pixels were then averaged to a spatial resolution of 0.1 × 0.1° for the maps in Fig. 1. OMI-based SO source detection and estimate. e Th methods for source detection and emission esti- 18,19,21 mate are based on a previously described algorithm that combines satellite measurements with reanalysis wind data (ECMWF interim reanalysis was used here). Wind information is matched with each OMI pixel. Emissions from about 500 continuously emitting point sources (or clusters of sources in close proximity), includ- ing 47 in India and 82 in China, are derived from OMI and wind data by tracking the downwind decay of the −1 −1 plumes. These sources have estimated SO emissions ranging from about 30 kt yr to more than 4000 kt yr . Due to the coarse spatial resolution of OMI (relative to a point source) and the limited precision of individual SO column observations, data spanning a year are analyzed together using a wind rotation scheme to align the wind 19 19 vectors of all overpasses considered . The emissions were estimated by fitting OMI columns to a plume function 18,21 consisting of coordinates, wind speed and direction, with a single parameter representing the total mass . Other fitting parameters, including an effective lifetime (5.9 hours), are specified . Emissions are then calculated as the ratio of mass to lifetime, effectively assuming a steady-state. The operational OMI retrievals use an effective air mass factor of 0.36 for all locations. In the emission estimate algorithm, OMI data for each emission source were adjusted using an air mass factor calculated for the location based on its elevation, surface albedo and sun/ 18,21 viewing geometry to better represent OMI sensitivities to the local source . Wet removal of pollutants by the summer monsoon rainfall causes a strong seasonality in air pollution in India, especially for aerosols . OMI generally also observes smaller SO columns over India during summer months (see monthly maps at https://avdc.gsfc.nasa.gov/pub/data/satellite/Aura/OMI/V03/L3/OMSO2m/ Monthly_mean_jpeg/). This seasonality may be partially attributed to the washout effects and the shorter lifetime of SO in summer, but may also reflect reduced coverage by OMI due to increased cloud cover during the mon - soon. For OMI-based emission estimates, since pixels with a small cloud fraction from an entire year are analyzed, the collection of data used to derive emissions for a given source may be more representative of non-summer conditions. The impact of this seasonal change in sampling on estimated emissions is currently unclear, but it is unlikely to significantly ae ff ct their long-term trend. Aircraft measurements of ΔSO /ΔCO ratio. Between 8 May and 11 June 2016, a twin engine Y-12 2 2 research aircraft was flown on 11 missions over the heavily industrialized Hebei Province of China. A modified, commercial pulsed-fluorescence detector (TEI Model 43 C) was used to measure ambient SO . A Picarro cavity ring-down spectrometer (Model G2401-m) was used to measure CO . Profiles were flown from near the surface to the top of the planetary boundary layer, at ~1500 m above ground. The Δ SO /ΔCO ratio was determined 2 2 from the deviation from background in plumes. Only data with significant correlation between ∆SO and ∆CO 2 2 (R > 0.6) are included in Fig. S1. Population data and population exposure to SO pollution. Population data for 2005, 2010, and 2015 from the Gridded Population of the World, Version 4 (GPWv4) were used in this study. For each of the three years, the (30 arc seconds) GPWv4 population count and nation identifier data were used to calculate the counts of Chinese and Indian population for each grid cell in Fig. 1. An annual growth rate was then estimated for each grid cell between 2005 and 2010, and between 2010 and 2015, to interpolate population data to other years. With OMI SO (Ω) and population count (P) data now on the same grid, the population-weighted SO column 2 2 amount (Ω ) for the entire domain with n grid cells can be calculated as: in = PΩ ii i=1 Ω= , in = ∑ (1) i=1 th where P and Ω are population count and OMI SO column for the i grid cell, respectively. i i 2 Data Availability. Level 2 Principal Component Analysis SO data from OMI are available from the Goddard Earth Science Data and Information Service Center (http://disc.sci.gsfc.nasa.gov/). Wind reanalysis data are avail- able from ECMWF (http://apps.ecmwf.int/datasets/data/interim-full-daily). Derived SO emissions are available from the global SO monitoring website at NASA Goddard Space Flight Center (https://so2.gsfc.nasa.gov). The GPWv4 population data are available from the Socioeconomic Data and Applications Center (SEDAC) in NASA’s Earth Observing System Data and Information System and hosted by Center for International Earth Science Information Network at Columbia University (http://sedac.ciesin.columbia.edu/data/collection/gpw-v4). Aircra ft measurements acquired during the ARIAS campaign are available upon request from X. Ren (ren@umd.edu). References 1. BP P.L.C. BP Statistical Review of WorldEnergy 2016. London, UK: BP p.l.c. https://www.bp.com/content/dam/bp/pdf/energy- economics/statistical-review-2016/bp-statistical-review-of-world-energy-2016-full-report.pdf. Accessed 27 Mar. 2017 (2016). 2. Huang, R.-J. et al. High secondary aerosol contribution to particulate pollution during haze events in China. Nature 514, 218–222 (2014). SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 5 www.nature.com/scientificreports/ 3. Chakraborty, A., Gupta, T. & Tripathi, S. N. Chemical composition and characteristics of ambient aerosols and rainwater residues during Indian summer monsoon: Insight from aerosol mass spectrometry. Atmos. Environ. 136, 144–155 (2016). 4. Lelieveld, J., Evans., J. S., Fnais, M., Giannadaki, D. & Pozzer, A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature 525, 367–371 (2015). 5. Zhang, Q. et al. Transboundary health impacts of transported global air pollution and international trade. Nature 543, 705–709 (2017). 6. Wang, G. et al. Persistent sulfate formation from London Fog to Chinese haze. Proc. Natl. Acad. Sci. 113, 13630–13635 (2016). 7. Environmental Protection Agency. National Emissions Inventory. Research Triangle Park, NC: Environmental Protection Agency, https://www.epa.gov/air-emissions-inventories/national-emissions-inventory-nei. Accessed 15 Mar. 2017 (2014). 8. Zhang, Q., He, K. B. & Huo, H. Cleaning China’s air. Nature 484, 161–162 (2012). 9. Lu, Z., Streets, D. G., de Foy, B. & Krotkov, N. A. Ozone monitoring instrument observations of interannual increases in SO emissions from Indian coal-fired power plants during 2005-2012. Environ. Sci. Technol. 47, 13993–14000 (2013). 10. Streets, D. G. et al. Emissions estimation from satellite retrievals: A review of current capability. Atmos. Environ. 77, 1011–1042 (2013). 11. Martin, R. V. Satellite remote sensing of surface air quality. Atmos. Environ. 42, 7823–7843 (2008). 12. Beirle, S. et al. Megacity emissions and lifetimes of nitrogen oxides probed from space. Science 333, 1737–1739 (2011). 13. Lamsal, L. N. et al. Application of satellite observations for timely updates to global anthropogenic NOx emission inventories. Geophys. Res. Lett. 38, L05810 (2011). 14. Fioletov, V. E. et al. Application of OMI, SCIAMACHY, and GOME-2 satellite SO retrievals for detection of large emission sources. J. Geophys. Res. Atmos. 118, 11399–11418 (2013). 15. Li, C. et al. Recent large reduction in sulfur dioxide emissions from Chinese power plants observed by the Ozone Monitoring Instrument. Geophys. Res. Lett. 37, L08807 (2010). 16. Fioletov, V. E., McLinden, C. A., Krotkov, N., Moran, M. D. & Yang, K. Estimation of SO emissions using OMI retrievals. Geophys. Res. Lett. 38, L21811 (2011). 17. He, H. et al. Response of SO and particulate air pollution to local and regional emission controls: A case study in Maryland. Earth’s Future 4, 94–109 (2016). 18. McLinden, C. A. et al. Space-based detection of missing sulfur dioxide sources of global air pollution. Nature Geoscience 9, 496–500 (2016). 19. Fioletov, V. E., McLinden, C. A., Krotkov, N. & Li, C. Lifetimes and emissions of SO from point sources estimated from OMI. Geophys. Res. Lett. 42, 1969–1976 (2015). 20. Li, C., Joiner, J., Krotkov, N. A. & Bhartia, P. K. A fast and sensitive new satellite SO retrieval algorithm based on principal component analysis: application to the Ozone Monitoring Instrument. Geophys. Res. Lett. 40, 6314–6318 (2013). 21. Fioletov, V. E. et al. A global catalogue of large SO sources and emissions derived from the Ozone Monitoring Instrument. Atmos. Chem. Phys. 16, 11497–11519 (2016). 22. Carn, S. A., Fioletov, V. E., McLinden, C. A., Li, C. & Krotkov, N. A. A decade of global volcanic SO emissions measured from space. Scientific Reports 7, 44095 (2017). 23. Krotkov, N. A. et al. Aura OMI observations of regional SO and NO pollution changes from 2005 to 2015. Atmos. Chem. Phys. 16, 2 2 4605–4629 (2016). 24. Lu, Z. et al. Sulfur dioxide emissions in China and sulfur trends in East Asia since 2000. Atmos. Chem. Phys. 10, 6311–6331 (2010). 25. Klimont, Z., Smith, S. J. & Cofala, J. The last decade of global anthropogenic sulfur dioxide: 2000–2011 emissions. Environ. Res. Lett. 8, 014003 (2013). 26. Shekar Reddy, M. & Venkataraman, C. Inventory of aerosol and sulphur dioxide emissions from India: I—Fossil fuel combustion. Atmos. Environ. 36, 677–697 (2002). 27. Warner, J. X. et al. Increased atmospheric ammonia over the world’s major agricultural areas detected from space. Geophys. Res. Lett. 44, https://doi.org/10.1002/2016GL072305 (2017). 28. Wang, S. X. et al. Emission trends and mitigation options for air pollutants in East Asia. Atmos. Chem. Phys. 14, 6571–6603 (2014). 29. Song, C. et al. Health burden attributable to ambient PM in China. Environ. Pollution 223, 575–586 (2017). 2.5 30. Stocker, T. F et al. (eds)Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge Univ. Press, 2013). 31. Levelt, P. F. et al. The Ozone Monitoring Instrument. IEEE Trans. Geosci. Remote 44, 1093–1101 (2006). 32. Li, C. et al. New-generation NASA aura ozone monitoring instrument (OMI) volcanic SO dataset: Algorithm description, initial results, and continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS). Atmos. Meas. Tech. 10, 445–458 (2017). 33. Dee, D. P. et al. e Th ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 137, 553–597 (2011). 34. International Institute for Applied Systems Analysis. Emissions of Air Pollutants for the World Energy Outlook 2011 Energy Scenarios, Laxenburg, Austria: International Institute for Applied Systems Analysis, http://pure.iiasa.ac.at/9766/1/XO-11-028.pdf. Accessed 27 Mar. 2017 (2011). 35. Hyvärinen, A.-P. et al. Effect of the summer monsoon on aerosols at two measurement stations in Northern India – Part 1: PM and BC concentrations. Atmos. Chem. Phys. 11, 8271–8282 (2011). 36. Center for International Earth Science Information Network. Gridded Population of the World, Version 4 (GPWv4): Population Count. Palisades, NY: NASA Socioeconomic Data and Applications Center. http://dx.doi.org/10.7927/H4X63JVC. Accessed 16 Feb. 2017 (2016). 37. International Institute for Applied Systems Analysis. Implications of energy trajectories from the World Energy Outlook 2015 for India’s air pollution, Laxenburg, Austria: International Institute for Applied Systems Analysis, http://www.worldenergyoutlook.org/media/ weowebsite/2015/Air_pollution_emissions_impacts_India_WEO2015_IIASA.pdf. Accessed 27 Mar. 2017 (2015). 38. Klimont, Z. et al. Projections of SO , NO and carbonaceous aerosols emissions in Asia. Tellus B 61, 602–617 (2009). 2 x 39. Lu, Z., Zhang, Q. & Streets, D. G. Sulfur dioxide and primary carbonaceous aerosol emissions in China and India, 1996–2010. Atmos. Chem. Phys. 11, 9839–9864 (2011). 40. Xing, J. et al. Projections of air pollutant emissions and its impacts on regional air quality in China in 2020. Atmos. Chem. Phys. 11, 3119–3136 (2011). 41. Zhao, Y., Zhang, J. & Nielsen, C. P. The effects of energy paths and emission controls and standards on future trends in China’s emissions of primary air pollutants. Atmos. Chem. Phys. 14, 8849–8868 (2014). 42. Zhang, Q. et al. Asian emissions in 2006 for the NASA INTEX-B mission. Atmos. Chem. Phys. 9, 5131–5153 (2009). Acknowledgements We thank the NASA Earth Science Division (ESD) Aura Science Team program for funding of OMI SO2 product development and analysis (Grant # 80NSSC17K0240). The Dutch and Finnish built OMI instrument is part of the NASA’s Earth Observing System (EOS) Aura satellite payload. The OMI project is managed by the Royal Meteorological Institute of the Netherlands (KNMI) and the Netherlands Space Agency (NSO). We thank Dr. SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 6 www.nature.com/scientificreports/ Steve Fetter (University of Maryland) for helpful comments. We thank the National Science Foundations of China (91544217) and US (Grant # 1558259) for supporting aircraft measurements over China presented in this study. C.L. acknowledges partial support from NASA’s Earth Science New Investigator Program (Grant # NNX14AI02G). Author Contributions C.L., C.M., V.F. and S.C. planned the research. C.L., N.K. and J.J. provided satellite retrievals. C.L., C.M., V.F. and D.S. analyzed the satellite data. H.H., X.R., Z.L. and R.D. provided and analyzed aircraft data. C.L. wrote the paper. All authors commented on the manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-14639-8. Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 SciEntific REPO R tS | 7: 14304 | DOI:10.1038/s41598-017-14639-8 7

Journal

Scientific ReportsSpringer Journals

Published: Nov 9, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off