Next Article in Journal
Changes in Domestic Energy and Water Usage during the UK COVID-19 Lockdown Using High-Resolution Temporal Data
Next Article in Special Issue
Qualitative Study on the Observations of Emissions, Transport, and the Influence of Climatic Factors from Sugarcane Burning: A South African Perspective
Previous Article in Journal
Custom-Made Foot Orthoses as Non-Specific Chronic Low Back Pain and Pronated Foot Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regional Features of Long-Term Exposure to PM2.5 Air Quality over Asia under SSP Scenarios Based on CMIP6 Models

by
Sungbo Shim
1,*,
Hyunmin Sung
1,
Sanghoon Kwon
1,
Jisun Kim
1,
Jaehee Lee
1,
Minah Sun
1,
Jaeyoung Song
1,
Jongchul Ha
1,
Younghwa Byun
1,
Yeonhee Kim
1,
Steven T. Turnock
2,3,
David S. Stevenson
4,
Robert J. Allen
5,
Fiona M. O’Connor
2,
Joao C. Teixeira
2,
Jonny Williams
6,
Ben Johnson
2,
James Keeble
7,8,
Jane Mulcahy
2 and
Guang Zeng
6
1
Innovative Meteorological Research Department, National Institute of Meteorological Sciences, Seogwipo-si 63568, Jeju-do, Korea
2
Met Office Hadley Centre, Exeter EX1 3PB, UK
3
University of Leeds Met Office Strategic (LUMOS) Research Group, School of Earth and Environment, University of Leeds, Leeds LS2 9JT, UK
4
School of GeoSciences, The University of Edinburgh, Edinburgh EH9 3FF, UK
5
Department of Earth and Planetary Sciences, University of California Riverside, Riverside, CA 92521, USA
6
National Institute for Water and Atmospheric Research, Wellington 6022, New Zealand
7
Department of Chemistry, University of Cambridge, Cambridge CB2 1TN, UK
8
National Centre for Atmospheric Science, University of Cambridge, Cambridge CB2 1EW, UK
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(13), 6817; https://doi.org/10.3390/ijerph18136817
Submission received: 23 May 2021 / Revised: 22 June 2021 / Accepted: 22 June 2021 / Published: 25 June 2021
(This article belongs to the Special Issue Advances in Air Pollution Meteorology Research)

Abstract

:
This study investigates changes in fine particulate matter (PM2.5) concentration and air-quality index (AQI) in Asia using nine different Coupled Model Inter-Comparison Project 6 (CMIP6) climate model ensembles from historical and future scenarios under shared socioeconomic pathways (SSPs). The results indicated that the estimated present-day PM2.5 concentrations were comparable to satellite-derived data. Overall, the PM2.5 concentrations of the analyzed regions exceeded the WHO air-quality guidelines, particularly in East Asia and South Asia. In future SSP scenarios that consider the implementation of significant air-quality controls (SSP1-2.6, SSP5-8.5) and medium air-quality controls (SSP2-4.5), the annual PM2.5 levels were predicted to substantially reduce (by 46% to around 66% of the present-day levels) in East Asia, resulting in a significant improvement in the AQI values in the mid-future. Conversely, weak air pollution controls considered in the SSP3-7.0 scenario resulted in poor AQI values in China and India. Moreover, a predicted increase in the percentage of aged populations (>65 years) in these regions, coupled with high AQI values, may increase the risk of premature deaths in the future. This study also examined the regional impact of PM2.5 mitigations on downward shortwave energy and surface air temperature. Our results revealed that, although significant air pollution controls can reduce long-term exposure to PM2.5, it may also contribute to the warming of near- and mid-future climates.

1. Introduction

Air pollution has now emerged as a leading global environmental health risk factor. In particular, long-term exposure to fine particulate matter, referred to as aerosols, is associated with increased rates of chronic bronchitis, reduced lung function, and increased mortality from lung cancer and heart disease [1]. The World Health Organization (WHO) reports that 92% of the world’s population resides in areas where the air-quality levels exceed the WHO’s ambient air-quality guidelines for the annual mean of particulate matter with a diameter of less than 2.5 μm (PM2.5) [2]. Moreover, approximately three million premature deaths occur annually as a result of long-term exposure to ambient air pollution [2,3]. Significantly high anthropogenic aerosol concentrations, especially recently in Asian regions, are linked to population growth and rapid industrialization. A recent report by the Organization for Economic Co–operation and Development (OECD) determined that South Korea, China, and India will endure the most significant economic damage and highest premature death rates in 2060 due to air pollution [4]. In addition, aerosols influence global and regional climate characteristics by altering the Earth’s radiative balance [5,6,7]. Previous studies [8,9,10] estimated that the present-day effective radiative forcing (ERF) by total aerosols is approximately −1 Wm−2, masking a considerable fraction of the well-mixed greenhouse gases (GHGs) induced by warming since the pre-industrial period.
Extensive research in recent decades has confirmed that the increase in surface PM2.5 concentrations in Asian regions, most notably since the middle 20th century, has been caused by the emission of primary aerosols and their precursors from industrial sources [7,11,12,13]. However, recent substantial efforts by the Chinese and Korean governments to mitigate anthropogenic emissions have resulted in a gradual decrease in the PM2.5 concentrations over East Asia [14,15]. Despite these efforts, current air pollution levels in Asian countries remain substantially higher than those of the United States and Europe. According to previous studies [16,17,18,19,20,21], the main cause of premature death in East Asia in the future will be the long-term exposure to elevated concentrations of ambient PM2.5. Moreover, regional changes in PM2.5 concentrations are greatly affected by local emissions and the long/distant transport of aerosols [22,23,24]. In addition, significant temperature variations in Asia during the 20th century may have been related to large regional historic anthropogenic aerosol emissions and their radiative effects, which are more regionally confined than those of GHGs due to their relatively short lifetimes [6,13]. Because of the widespread environmental impacts such as visibility impairment and severe threats posed to human health, understanding the effects of fine particulate matter on human health and climate interactions have become critical to regional policymakers.
Recent chemistry–climate modeling studies have attempted to account for the geographical variations in aerosol emissions and transportation and consider the effect of these factors on the climate system [25,26,27,28]. In addition, experiments performed as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) [29] and the Aerosol Chemistry Model Intercomparison Project (AerChemMIP) [30] have contributed significantly to the multimodel (including interactive tropospheric chemistry and aerosols) evaluation of historical and future changes in air pollutants [10,31,32,33,34]. Aerosol concentrations are expected to reduce globally in the future, albeit at different paces over different regions, with reductions initially expected to occur in developed countries in Western Europe and North America. However, Asia, the world’s most populated and rapidly developing region, has been the largest source of aerosols and their precursors over the last few decades [35,36]. Therefore, the changes in regional aerosol concentrations caused by future air pollution mitigation could have important consequences for human health and climate change in the global and Asian countries.
This study aimed to assess the annual mean surface PM2.5 concentration and air quality in Asian regions (East Asia, South Asia, and Southeast Asia) and conduct future projections using simulations from CMIP6 multimodel ensembles. Future scenarios have been developed for the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6: the shared socioeconomic pathways (SSP), which were combined with the Representative Concentration Pathways (RCP) [29]. The qualitative and quantitative aspects of these new scenarios have been extensively discussed in the literature [37,38,39,40]. In the SSP scenarios, the level and speed of air pollution controls in each region are implemented differently. While SSP3 (regional rivalry) and SSP4 (inequality) assume the slowest deployment of air pollution controls, SSP1 (sustainability) and SSP5 (conventional development) adopt the fastest and widest implementations, as shown in Table 1 [41,42]. Moreover, the SSP2 scenario implements significant advances in air pollution control. Air pollution controls in the SSP scenarios depend on the income levels of the countries. Therefore, each country in the Asian region will implement air-pollution controls with various schedules according to their individual institutional, financial, and technological capacities. Detailed information on air pollution control per country according to the SSP scenarios can be found in Rao et al. [43]. This study was designed to explore the future changes in PM2.5 concentrations and regional effects on air quality and climate using the new CMIP6 scenarios. A total of 38 simulations were analyzed from nine CMIP6 climate models for historical and future SSP scenarios. The CMIP6 models, experiments, and methodology are described in Section 2, followed by an analysis of the model data in Section 3. Summaries and Discussions are presented in Section 4.

2. Data and Methodology

2.1. Simulation Data from CMIP6 Archive

In this study, we used 38 simulations from nine CMIP6 models (Table 2) for the first realizations that were available at the time of this study (up to February 2021) from the Earth System Grid Federation (ESGF). All the model outputs were bilinearly interpolated onto a grid with a horizontal resolution of 1.875° × 1.25°, which is also used by the UK Earth System Model (UKESM1-0-LL) [44]. The bilinear interpolation method tends to underestimate local maxima but does not create fictitious local peaks. The multimodel ensemble (MME) was examined by using the arithmetic mean of the model output, with the same weight assigned to each model. Six different experiments were selected for each CMIP6 model: historical simulation (Historical) was forced by the observed history of anthropogenic sources and natural forcing trends over a 164-year period, from 1850 to 2014.
The future changes in air pollutants were then projected under different SSP scenarios throughout the 21st century (Table 3). For CMIP6, the four SSPs were categorized as SSP1–2.6 (SSP126), SSP2–4.5 (SSP245), SSP3–7.0 (SSP370), and SSP5–8.5 (SSP585) for the period from 2015 to 2100. O’Neill et al. [37] provide a more detailed description of the CMIP6 simulation design for the future scenarios. CMIP6 model data for the SSP3-7.0-lowNTCF scenario (SSP370 with cleaner air-quality policies; SSP370-lowNTCF) over the 2015–2055 period, derived from AerChemMIP [30], was also used in this study. The baseline historical simulation period of 1995–2014 was selected as present day (PD) in this study, whereas the near-future period of 2021–2040 and mid-future period of 2041–2060 were selected for the analysis of future changes in air pollution. Overall, regional analyses of three Asian subregions (Figure 1) were selected for this study, following the domain of Iturbide et al. (2020) [53] prepared for the IPCC sixth assessment report (AR6).

2.2. Satellite Data

This study used satellite-derived global annual mean surface PM2.5 concentration grid data from the NASA Socioeconomic Data and Application Center (SEDAC) for evaluating the simulated surface PM2.5 concentrations. These data were estimated using a combined geophysical–statistical method with information from aerosol optical depth (AOD) retrievals from the NASA Moderate Resolution Imaging Spectro Radiometer (MODIS), Multiangle Imaging Spectro Radiometer (MISR), and the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) instruments with the GEOS–Chem chemical transport model. The data were then calibrated to global ground-based measurements of PM2.5 concentrations for large-scale health and environmental research [54,55]. This high-resolution (0.01° × 0.01°) dataset is provided from the SEDAC website (available at https://sedac.ciesin.columbia.edu (accessed on 5 April 2021)) and covers the global land surface from 70° north to 55° south over the 1998–2016 period. Surface PM2.5 observational data were re-gridded onto the same resolution as the CMIP6 multimodel ensemble product (1.875° × 1.25°) for evaluation purposes. We focused on the 1998–2014 period to match the present-day (1995–2014) CMIP6 historical simulations.

2.3. Methodology

This study explored the impact of fine particulate matter on the air quality in Asia using SSP scenarios from the CMIP6 archive. However, not enough models provided a direct calculation of PM2.5 concentrations and as such we had to use an approximation that accounts for all aerosol components that were consistent across all CMIP6 models. Therefore, we calculated the PM2.5 concentrations using the offline method, as shown in Equation (1). Here, BC, OC, NH4, SO4, NO3, DU, and SS represented black carbon (BC), total organic carbon (OC), ammonium (NH4), sulfate (SO4), nitrate (NO3), dust (DU), and sea-salt (SS) particles, respectively, from the lowest model level in the CMIP6 individual model following previous studies [21,26,28,33,35,56]. The factors 0.1 and 0.25 indicated the DU and SS PM2.5 size fraction. Note that only a few CMIP6 models include the simulation of ammonium and nitrate particles in their aerosol-chemistry schemes (CESM2-WACCM model included NH4 particle and the GFDL-ESM4, GISS-E2-1G models included both NH4 and NO3). The analysis area was divided into three subdomains: East Asia (EAS), South Asia (SAS), and South East Asia (SEA) (Figure 1). The analyses of the regional annual mean PM2.5 gridded data were conducted based on a multimodel ensemble (MME).
Estimated PM2.5 = BC + OC + NH4 + SO4 + NO3 + (0.1 × DU) + (0.25 × SS)
This study also attempted to estimate the effects of long-term exposure to fine particulate matter using the annual mean surface PM2.5 concentrations calculated from Equation (1). The air-quality levels in the six categories were classified based on the annual mean PM2.5 from the WHO air-quality guidelines and interim targets (Table 4), namely the AQG (air-quality guideline), IT-3 (interim target 3), IT-2 (interim target 2), IT-1 (interim target 1), OT (over target), and ST (significant target). The higher AQI values indicated a greater level of air pollution and therefore pose a higher risk to health, which could lead to the premature death of residents. For example, an AQI value of ‘AQG’ denoted a good air-quality level with little or no potential to affect public health. However, an AQI value of IT-1 suggests that air pollution may have contributed to an approximately 15% higher premature mortality risk relative to the AQG level [57]. Therefore, the WHO guidelines for outdoor particulate matter recommend that the annual average PM2.5 does not exceed 10 µg/m3. In this study, changes in the annual mean surface PM2.5 concentrations and AQI index, and their effects on regional future climates were analyzed according to the SSP scenarios from the CMIP6 simulated data.

3. Results

3.1. Evaluating the Estimated PM2.5 from CMIP6 Models in the Present–Day Period

The performance of fine particulate matter was evaluated from the nine CMIP6 models (Table 2). The annual mean surface PM2.5 concentrations for the globe and for the EAS, SAS, and SEA domain regions were estimated using historical simulations in the individual CMIP6 model and MME, and the results were compared with the available satellite-based PM2.5 data for validation. The CMIP6 MME mean PM2.5 (black horizontal line in Figure 2) showed good agreement with the satellite-derived PM2.5 (black open circle in Figure 2) for the present-day period. For both the CMIP6 model-simulated and satellite-derived data, annual mean surface PM2.5 concentrations in EAS and SAS were more than twice as high than worldwide and in the SEA region. The CMIP6 MME for the present-day annual mean surface PM2.5 was 9.6 ± 2.5 µg/m3 globally, 23.6 ± 5.3 µg/m3 in EAS, 23.7 ± 6.2 µg/m3 in SAS, and 9.6 ± 2.8 µg/m3 in SEA. The simulated PM2.5 concentrations in EAS and SAS were marginally underestimated compared to the satellite data. This might be because only a few CMIP6 models account for nitrate aerosols as described in Section 2. And the deviation between the individual CMIP6 models was also greater than those of the global and SEA regions. The larger model diversity over the EAS and SAS regions was consistent with a previous study [33] that demonstrated that the inter-model differences might be attributed to different simulations of historical changes in the anthropogenic aerosol components.
The regional spatial distribution of the annual mean PM2.5, estimated based on the CMIP6 MME, was compared to satellite-derived PM2.5 concentrations for the present day (Figure 3). The results indicate that the industrialized and highly populated countries of China and India are exposed to the highest regional PM2.5 concentrations. For this reason, the annual mean concentration of PM2.5 in recent decades was the highest in the Beijing and Zhangzhou regions of central China, with values exceeding 50 µg/m3 (Figure 3a). Moreover, significant variation was simulated between the CMIP6 individual models in areas with large anthropogenic emission sources, such as Eastern China and Northern India (Figure 3b). Conversely, the MME PM2.5 mean and standard deviation was relatively low in southern India, Japan, and Indonesia. Overall, the PM2.5 concentrations in most of the regions analyzed in this study, except for Japan, the Philippines, and eastern Indonesia, exceeded the WHO air-quality guidelines (10 µg/m3 annual mean).
In this study, we adopted the CMIP6 MME analysis to understand the uncertainty due to differences in the physical process of various climate models. The spatial distribution of the annual mean PM2.5 for the present-day period, calculated from CMIP6 historical simulations, showed similar distribution patterns to those of satellite-derived PM2.5 data (Figure 3a,c). However, the simulated-PM2.5 concentrations were underestimated compared to the satellite-derived PM2.5 concentrations obtained for eastern China, northern India, and Thailand (Figure 3d). Conversely, for western China, Pakistan, and Indonesia, the CMIP6 model concentrations were overestimated relative to the satellite data. Moreover, the diversity between the CMIP6 models was large over eastern China and northern India, which experience high levels of PM2.5. Despite model diversity and regional biases (Figure 3b,d), the simulated PM2.5 concentrations from the CMIP6 MME showed reasonable domain-averaged values in the present-day climate compared to the satellite-derived PM2.5 concentrations (Figure 2). Therefore, the historical CMIP6 MME simulations were used as a reference period (1995–2014) for near- and mid-future climate changes for the remainder of this study.

3.2. Future Changes in Simulated PM2.5 Concentrations and the Air-Quality Index

Figure 4 shows the projected global and regional changes in primary aerosols and their precursor emissions relative to the present day. This data was used as the input data for the CMIP6 models. Anthropogenic aerosol and precursor (organic carbon, black carbon, and sulfur dioxide) emissions used in each CMIP6 model were obtained from the same dataset. It is worth noting that the emissions of natural aerosol sources such as dust and sea-salt are different, depending on the physical configuration of the individual CMIP6 model (not shown).
Overall, except for the SSP3-7.0 scenario, the emissions of anthropogenic aerosols and precursors for EAS, SEA, and worldwide showed a decreasing tendency in the near and mid-future. In the SAS region, SO2 showed an increasing trend (up to 50%) in all the near-future scenarios and a continuous increase in the mid-future for the SSP5-8.5, SSP3-7.0, and SSP2-4.5 scenarios (Figure 4a). In addition, future OC and BC emissions in SAS are also expected to increase in the near future (Figure 4b,c) for all the scenarios except SSP1-2.6 and SSP5-8.5, which include the rapid implementation of air pollution controls. For the SSP3-7.0 scenario, which specifies inadequate air pollution controls, anthropogenic aerosols and precursor emissions are expected to increase or remain at the present-day levels in all regions.
Future changes in the annual mean surface PM2.5 concentrations in the Asian regions were examined using the different CMIP6 SSP scenarios (Figure 5). For the significant air-quality control scenarios, SSP1-2.6 and SSP5-8.5, a considerable decrease in the annual mean surface PM2.5 was predicted for the near future in EAS (Figure 5a). The results indicated a decrease in PM2.5 of more than 50% by the end of the 21st century compared to present-day levels. These changes are driven by the large emission controls of anthropogenic aerosols and their precursors (Figure 4). The decreasing trend in the SSP2-4.5 scenario was similar to that in the SSP5-8.5 scenario; however, higher PM2.5 concentrations were simulated in the near-future period. In the SSP3-7.0 scenario, which includes weak air-quality controls, the annual mean PM2.5 concentrations increased until the mid-21st century, and then subsequently decreased to levels similar to the present day by the end of the 21st century.
In the SAS region, future changes of the annual mean PM2.5 in the SSP1-2.6 and SSP3-7.0 scenarios were similar to those in EAS. However, significantly large diversities were simulated between the CMIP6 models (Figure 5b), particularly for the far future. Moreover, unlike the predictions for EAS, changes in the annual mean PM2.5 for the SSP2-4.5 and SSP5-8.5 scenarios were characterized by a marginal increase in the near-future period (green and red solid lines in Figure 5b). It was also found that the future changes in PM2.5 concentrations in SAS exceeded those in EAS, particularly for the SSP3-7.0 scenario, which includes the weak implementation of air-quality controls. In the SEA region, both the changes in the annual mean PM2.5 and the model diversity were smaller than in the EAS and SAS regions (Figure 5c). However, a continuous increase (approximately 30% in the far future) in the PM2.5 in the SEA region was found for the SSP3-7.0 scenario, whereas the PM2.5 concentrations simulated for the EAS and SAS regions decreased after the mid-future period.
In this study, we calculated the AQI index for the present-day and near-future periods using the simulated annual surface mean PM2.5 concentrations obtained from the CMIP6 MME mean (Figure 6). Higher AQI values are an indicator of increasing long-term exposure to PM2.5, which could pose a severe threat to human health and increase the risk of premature death (Table 4). The present-day AQI values for the historical simulations in the EAS and SAS regions were typically above the level of ’interim target 2′ (the cyan-colored areas in Figure 6a). In particular, the AQI values over northern India and eastern China, around mega cities such as Beijing, Shanghai, Delhi, and Kolkata were at ‘over target’ level (the orange-colored area in Figure 6a). Based on the simulated PM2.5 concentrations for the present day, we determined that the PM2.5 levels for Japan, western China, and eastern SEA, situated far from air pollution source regions, were below the specified WHO AQG annual mean level of 10 µg/m3 (the indigo-colored area in Figure 6a). These results are consistent with a WHO report published in 2016 on annual mean PM2.5 levels in the relevant Asian regions [2].
In the SSP1-2.6 scenario, the AQI values in EAS are expected to improve substantially in the near future. In particular, the AQI level in China was lowered to the interim target range, and the AQI values for the Korean peninsula were below the specified WHO AQG levels (Figure 6b). The ‘over target’ area in northern SAS also showed a decrease compared to present-day levels. The spatial distributions of the AQI values showed similar patterns for the SSP2-4.5 and SSP5-8.5 scenarios, with AQI levels showing slightly improvement in EAS and worse in SAS compared to the present day (Figure 6c,e). These results were attributed to increases in the annual mean PM2.5 concentrations in SAS due to higher anthropogenic aerosol and precursor concentrations in the near-future period (Figure 4 and Figure 5b). In the SSP3-7.0 scenario, the AQI values were expected to increase significantly in EAS, SAS, and northern SEA (Figure 6d). In particular, the AQI values over the northern SAS and eastern EAS regions were calculated as ‘significantly over target’ (the red-colored area in Figure 6d). We also estimated resident population in areas exceeding the WHO interim target using the total population dataset from the Inter-Sectoral Impact Model Inter-Comparison Project phase 2b (ISIMIP2b) [59]. In the SSP3 scenario, more than 1.6 billion people (around 1 billion in SAS and around 0.6 billion in EAS) will be exposed to air pollution exceeding the IT levels of AQI in the near future. These changes of AQI values in the mid-future are much more pronounced than for the near future (Figure 7). In the SSP1-2.6 scenario, the AQI values in EAS, excluding inland China, were expected to be IT-3 (#1) or the lowest levels (AQG, #0) of AQI for the mid-future period. Conversely, in the SSP3-7.0 scenario, it was found that the AQI values in more than 50% of SAS regions were OT or ST levels.
The annual mean PM2.5 concentrations and AQI levels in EAS were shown to improve overall in the future with the implementation of air-quality controls (Figure 5 and Figure 6). However, the predicted AQI values remained above the level of ‘interim target 2′ and ‘interim target 1′ in eastern China. Moreover, our results indicated that the AQI values for the northern region of SAS in the near future will exceed present-day levels. CMIP6 models simulated that most regions over China and India in the near-future period will be exposed to levels of PM2.5 concentrations that exceed the WHO air-quality guideline values. In addition, the proportion of populations over the age of 65 years old in China and India is projected to increase significantly under the SSP scenarios (Figure 8). This ratio could be as high as 35% by 2050, then increasing to up to 60% by 2100 (Figure 8a). Rapid population ageing coupled with AQI levels higher than the WHO guidelines in the future may increase the risk of premature mortality from the long-term exposure to PM2.5.

3.3. Regional Response to Future Air Pollution Mitigation

From the middle of the 20th century to the present day, the annual mean surface PM2.5 concentrations across Asia have increased considerably due to industrialization (Figure 9). The aerosol-radiative effects due to rapidly increasing anthropogenic aerosols have played a critical role in cooling in recent decades, partially offsetting the GHGs warming [13,60,61]. Aerosol can also affect the SW indirectly by altering cloud properties such as albedo and lifetime, but have more uncertain effects. For this reason, regional PM2.5 concentrations in CMIP6 historical simulations are closely correlated with changes in the clear-sky surface downwelling shortwave radiative flux (SW). Particularly, clear-sky SW in EAS showed a significant negative trend, decreasing from 260 to 240 Wm−2 (blue circles in Figure 9). Therefore, the implementation of air pollution controls can accelerate global and regional warming by recovering reduced surface radiation in the present day. According to previous study [62], the warming in response to reduced anthropogenic aerosols in China is likely already occurring recently.
To understand the warming effects by anthropogenic aerosol reductions, we examined the differences between the weak air-quality control simulations (SSP3-7.0) and the strong air quality control simulations (SSP3-7.0-lowNTCF), which assumes the implementation of stringent air-quality control policies under the same conditions as the SSP3-7.0 scenario (Table 3, Figure 10). The significant difference in annual mean PM2.5 concentrations between the SSP3-7.0 and SSP3-7.0-lowNTCF (the orange and cyan lines in Figure 5) scenarios contributed to changes in the total aerosol optical depth (AOD) and clear-sky SW in the middle of the 21st century (Figure 10a–f). The regional AODs were reduced and the surface SW clear-sky increased considerably in SSP3-7.0-lowNTCF compared to SSP3-7.0 simulations. In particular, these changes in the EAS and SAS regions, which recorded high levels of PM2.5, showed relatively large differences compared to those in SEA. Moreover, warming trends (approximately +0.5 K in 2050) were simulated in response to increases in the surface SW clear-sky (Figure 10g–i). However, significant model diversities were found for the surface air temperature changes. Despite large model uncertainties of surface air temperature in all regions, implementations of strong air-quality controls seem to contribute to regional warming in the future. Therefore, reducing long-term exposure to PM2.5 can impact the rate of regional warming, so should be considered together when implementing climate change mitigation initiatives. Detailed warming mechanisms involved in climate sensitivity, aerosol-cloud interaction, and the temperature advections related to large-scale circulation require further research to understand CMIP6 model uncertainty; however, these factors are beyond the scope of this study.

4. Conclusions

This study performed an initial analysis of the long-term changes in PM2.5, AQI, and climate responses for three major regions in Asia under different future climate and air pollution control scenarios. The regional PM2.5 concentrations were estimated using simulations from nine CMIP6 models. The model performance of the estimated PM2.5 concentrations from historical simulations was evaluated against satellite-derived data from NASA’s Socioeconomic Data and Application Center. Elevated annual mean PM2.5 concentrations (>50 µg/m3) were simulated for the EAS and SAS regions, which are heavily industrialized and densely populated. Conversely, simulated present-day annual mean PM2.5 concentrations in SEA were relatively lower than in EAS and SAS. Overall, the simulated annual mean PM2.5 concentrations for all the analyzed regions except for Japan, the Philippines, and Indonesia exceeded the WHO’s recommended air-quality guidelines (10 µg/m3) for the present day. Moreover, a comparison of the estimated PM2.5 levels from the CMIP6 models with satellite data revealed an underestimation in the data for the three Asian domains, and particularly for EAS and SAS. Despite model diversity and regional biases, the simulated PM2.5 showed reasonable domain averaged values in the present day compared to the satellite-derived PM2.5 concentrations.
Future changes in the simulated annual mean PM2.5 were also examined for the various SSP scenarios. Rapid decreases in SSP1-2.6 and SSP5-8.5 for the scenarios assuming strict air-quality control were shown across the Asian regions in the near future. Unlike for EAS, the scenarios predicted an increase in PM2.5 concentrations in SAS in the near future. The changes in the annual mean value and model diversity for SEA were smaller than in the EAS and SAS regions. For the slowest deployment of air pollution controls (SSP3-7.0 scenario), annual mean PM2.5 showed an increase across all regions until the mid-21st century. We also attempted to understand the long-term exposure risk to PM2.5 using the AQI based on the WHO air-quality guidelines and interim targets. The present-day AQI level in EAS and SAS typically exceeded ‘interim target 2′, and several regions exceeded the range of the WHO interim target. In the SSP1-2.6, the CMIP6 models predicted a substantial decrease in AQI values in EAS and SAS. Similar trends were simulated under the SSP2-4.5 and SSP5-8.5 for EAS, whereas the AQI values in SAS were predicted to increase due to increasing PM2.5 levels in the near future. In the SSP3-7.0, the number of regions beyond interim target is expected to increase significantly. Moreover, the risk of premature mortality from long-term exposure to PM2.5 could significantly increase in China and India due to the rapid population ageing in the future.
This study is among a few climate-modeling studies that explored the impact of future air quality in Asia under the SSP scenarios. Overall, our results indicated that the CMIP6 climate model simulations successfully reproduced the regional surface PM2.5 concentrations for the present day. The results in this study were consistent with a previous study [33] that demonstrated that the decreasing (increasing) trends of simulated regional PM2.5 are related to the strong (weak) implementation of future air-quality controls. In addition, this study analyzed the AQIs based on the WHO guidelines for the three major Asian regions and investigated the regional impact of air pollution controls implementations on PM2.5 levels. An additional finding of this study was the potential acceleration of global warming in Asia with future decreases in anthropogenic aerosol emissions reducing the radiative cooling effect of aerosols.
Despite the significant findings of this study, several limitations were noted. We adopted the CMIP6 multimodel ensemble analysis to understand the impact of fine PM on regional air quality and their uncertainties. However, the estimated PM2.5 concentrations from the CMIP6 models were underestimated compared to the satellite-derived data. Only a few models consider ammonium nitrate as aerosol components in their schemes, which may be one of the reasons for underestimation of PM2.5 concentrations. In addition, the large model diversity may be attributed to the differences in the aerosol microphysical processes among the CMIP6 individual models. Detailed causes of model uncertainties should be explored in further studies; however, they are beyond the scope of this study. This study also emphasized the effects of premature mortality due to increasing PM2.5 concentrations, coinciding with a rapidly ageing population. However, a more detailed risk analysis with the available data from the CMIP6 simulations is required to verify this conclusion. These shortcomings will be addressed in our future research.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijerph18136817/s1, Table S1: Description of CMIP6 models used in this study.

Author Contributions

S.S. wrote the manuscript and contributed to the graphics. H.S., S.K., J.K. (Jisun Kim), J.L., M.S., J.S., J.H., Y.B., Y.K., S.T.T., D.S.S., R.J.A., F.M.O., J.C.T., J.W., B.J., J.K. (James Keeble), J.M., and G.Z. contributed to the discussion and revision of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Korea Meteorological Administration Research and Development Program “Development and Assessment of Climate Change Scenario” under Grant (KMA2018-00321). F.M.O., J.C.T., B.J., and J.M. have been supported by the BEIS and DEFRA Met Office Hadley Centre Climate Programme (GA01101). S.T.T. has been supported by the UK–China Research and Innovation Partnership Fund through the Met Office Climate Science for Service Partnership (CSSP) China as part of the Newton Fund. J.K. was financially supported by NERC through NCAS (grant no. R8/H12/83/003). G.Z. was supported by the NZ Government’s Strategic Science Investment Fund (SSIF) through the NIWA programme (CACV).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. Datasets from CMIP6 simulations are available through the Earth System Grid Federation (ESGF).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Burden of Disease from Ambient and Household Air Pollution. 2016. Available online: https://who.int/phe/health_topics/outdoorair/databases/en/ (accessed on 2 August 2020).
  2. WHO. Ambient Air Pollution: A Global Assessment of Exposure and Burden of Disease. 2016. Available online: https://who.int/publications/air-pollution-global-asessment/en/ (accessed on 2 August 2020).
  3. Kuehn, B.M. WHO: More than 7 million air pollution deaths each year. JAMA 2014, 311, 1486. [Google Scholar] [CrossRef] [PubMed]
  4. OECD. The Economic Consequences of Outdoor Air Pollution; OECD Publishing: Paris, France, 2016. [Google Scholar]
  5. IPCC. Climate Change 2014: Anthropogenic and Natural Radiative Forcing. In Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; pp. 659–740. [Google Scholar]
  6. Kim, J.; Gu, Y.; Liou, K.N.; Park, R.; Song, C. Direct and semi-direct radiative effects of anthropogenic aerosols in the western United States: Seasonal and geographical variations according to regional climate characteristics. Clim. Chang. 2012, 111, 859–877. [Google Scholar] [CrossRef]
  7. Seo, J.; Shim, S.; Kwon, S.-H.; Boo, K.-O.; Byun, Y.-H.; O’Connor, F.; Johnson, B.; Dalvi, M.; Folberth, G.; Teixeira, J.; et al. The impacts of aerosol emissions on historical climate in UKESM1. Atmosphere 2020, 11, 1095. [Google Scholar] [CrossRef]
  8. Boucher, O.; Randall, P.; Artaxo, P.; Bretherton, C.; Feingold, G.; Forster, P.; Kerminen, V.-M.; Kondo, Y.; Liao, H.; Lohmann, U.; et al. Clouds and Aerosols. In Climate Change 2013: The Physical Sciences Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  9. Smith, C.J.; Kramer, R.J.; Myhre, G.; Alterskjaer, K.; Collins, W.; Sima, A.; Boucher, O.; Dufresne, J.-L.; Nabat, P.; Michou, M.; et al. Effective radiative forcing and adjustments in CMIP6 models. Atmos. Chem. Phys. 2020, 20, 9591–9618. [Google Scholar] [CrossRef]
  10. O’Connor, F.M.; Abraham, N.L.; Dalvi, M.; Folberth, G.A.; Griffiths, P.T.; Hardacre, C.; Johnson, B.T.; Kahana, R.; Keeble, J.; Kim, B.; et al. Assessment of pre-industrial to present-day anthropogenic climate forcing in UKESM1. Atmos. Chem. Phys. 2021, 21, 1211–1243. [Google Scholar] [CrossRef]
  11. Ramanathan, V.; Crutzen, P.J.; Kiehl, J.T.; Rosenfeld, D. Aerosols, climate, and the hydrological cycle. Science 2001, 294, 2119–2124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Smith, S.J.; Aardenne, J.; Klimont, Z.; Andres, R.J.; Volke, A.; Arias, S.D. Anthropogenic sulfur dioxide emission: 1850–2005. Atmos. Chem. Phys. 2011, 11, 1101–1116. [Google Scholar] [CrossRef] [Green Version]
  13. Shim, S.; Kim, J.; Yum, S.S.; Lee, H.; Boo, K.-O.; Byun, Y.-H. Effects of anthropogenic and natural forcing on the summer temperature variations in East Asia during the 20th century. Atmosphere 2019, 10, 690. [Google Scholar] [CrossRef] [Green Version]
  14. Hong, S.Y.; Lee, J.J.; Lee, J.Y.; Kim, Y.P. Comparison of the fine particle concentrations in Seoul and other foreign mega-cities. Par. Aerosol Res. 2008, 4, 1–7. [Google Scholar]
  15. Han, S.H.; Kim, Y.P. Long-term trends of the concentrations of mass and chemical composition in PM2.5 over Seoul. J. Korean Soc. Atmos. Environ. 2015, 31, 143–156. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, Y.-M. Premature deaths attributable to exposure to ambient fine particulate matter in the republic of Korea: Implication for future environmental health. J. Korean Med. Sci. 2018, 33, e243. [Google Scholar] [CrossRef]
  17. Kim, J.-H.; Oh, I.-H.; Park, J.-H.; Cheong, H.-K. Premature deaths attributable to long-term exposure to ambient fine particulate matter in the Republic of Korea. J. Korean Med. Sci. 2018, 33, e251. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, Q.; Wang, J.; Zhou, J.; Ban, J.; Li, T. Estimation of PM2.5-associatied disease burden in China in 2020 and 2030 using population and air quality scenarios: A modelling study. Lancet Planet Heath 2019, 3, 71–80. [Google Scholar] [CrossRef] [Green Version]
  19. Hwang, I.C. Health effects of aging and fine particle matter. Seoul Inst. Policy Rep. 2019, 287, 20. [Google Scholar]
  20. Park, S.; Allen, R.J.; Lim, C.H. A likely increase in fine particulate matter and premature mortality under future climate change. Air Qual. Atmos. Health 2020, 13, 143–151. [Google Scholar] [CrossRef] [Green Version]
  21. Silva, R.A.; West, J.J.; Lamarque, J.F.; Shindell, D.T.; Collins, W.J.; Faluvegi, G.; Folberth, G.A.; Horowitz, L.W.; Nagashima, T.; Naik, V.; et al. Future global mortality from changes in air pollution attributable to climate change. Nat. Clim. Chang. 2017, 7, 647–651. [Google Scholar] [CrossRef] [Green Version]
  22. Lee, S.; Ho, C.-H.; Lee, Y.G.; Choi, H.-J.; Song, C.-K. Influence of transboundary air pollutants from China on the high-PM10 episode in Seoul, Korea for the period October 16-20. Atmos. Environ. 2013, 77, 430–439. [Google Scholar] [CrossRef]
  23. Park, S.-S.; Jung, S.-A.; Gong, B.-J.; Cho, S.-Y.; Lee, S.-J. Characteristics of PM2.5 haze episodes revealed by highly time-resolved measurements at an air pollution monitoring supersite in Korea. Aerosol Air Qual. Res. 2013, 13, 957–976. [Google Scholar] [CrossRef]
  24. Koo, Y.-S.; Choi, D.-R.; Kwon, H.-Y.; Jang, Y.-K.; Han, J.-S. Improvement of PM10 prediction in East Asia using inverse modeling. Atmos. Environ. 2015, 106, 318–328. [Google Scholar] [CrossRef]
  25. Colette, A.; Bessagnet, B.; Vautard, R.; Szopa, S.; Rao, S.; Schucht, S.; Klimont, Z.; Menut, L.; Clain, G.; Meleux, F.; et al. European atmosphere in 2050, a regional air quality and climate perspective under CMIP5 scenarios. Atmos. Chem. Phys. 2013, 13, 7451–7471. [Google Scholar] [CrossRef] [Green Version]
  26. Fiore, A.M.; Naik, V.; Spracklen, D.V.; Steiner, A.; Unger, N.; Prather, M.; Bergmann, D.; Cameron-Smith, P.J.; Cionni, I.; Collins, W.J.; et al. Global air quality and climate. Chem. Soc. Rev. 2012, 41, 6663–6683. [Google Scholar] [CrossRef] [Green Version]
  27. Tai, A.P.K.; Mickley, L.J.; Jacob, D.J. Impact of 2000–2050 climate change on fine particulate matter (PM2.5) air quality inferred from a multi-model analysis of meteorological modes. Atmos. Chem. Phys. 2012, 12, 11329–11337. [Google Scholar] [CrossRef] [Green Version]
  28. Chowdhury, S.; Dey, S.; Smith, K.R. Ambient PM2.5 exposures and expected premature mortality to 2100 in India under climate change scenarios. Nat. Commun. 2018, 9, 318. [Google Scholar] [CrossRef] [Green Version]
  29. Eyring, V.; Bony, S.; Meehl, G.A.; Senior, C.A.; Stevens, B.; Stouffer, R.J.; Taylor, K.E. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev. 2016, 9, 1937–1958. [Google Scholar] [CrossRef] [Green Version]
  30. Collins, W.J.; Lamarque, J.-F.; Schulz, M.; Boucher, O.; Eyring, V.; Hegglin, M.I.; Maycock, A.; Myhre, G.; Prather, M.; Shindell, D.; et al. AerChemMIP: Quantifying the effects of chemistry and aerosols in CMIP6. Geosci. Model Dev. 2017, 10, 585–607. [Google Scholar] [CrossRef] [Green Version]
  31. Allen, R.J.; Turnock, S.; Nabat, P.; Neubauer, D.; Lohmann, U.; Olivie, D.; Oshima, N.; Michou, M.; Wu, T.; Zhang, J.; et al. Climate and air quality impacts due to mitigation of non-methane near-term climate forcers. Atmos. Chem. Phys. 2020, 20, 9641–9663. [Google Scholar] [CrossRef]
  32. Griffiths, P.T.; Murray, L.T.; Zeng, G.; Archibald, A.T.; Emmons, L.K.; Galbally, I.; Hassler, B.; Horowitz, L.W.; Keeble, J.; Liu, J.; et al. Tropospheric ozone in CMIP6 simulations. Amos. Chem. Phys. 2020. in review. [Google Scholar] [CrossRef] [Green Version]
  33. Turnock, S.T.; Allen, R.J.; Andrews, M.; Bauer, S.E.; Deushi, M.; Emmons, L.; Good, P.; Horowitz, L.; John, J.G.; Michou, M.; et al. Historical and future changes in air pollutants from CMIP6 models. Atmos. Chem. Phys. 2020, 20, 14547–14579. [Google Scholar] [CrossRef]
  34. Allen, R.J.; Horowitz, L.W.; Naik, V.; Oshima, N.; O’Connor, F.M.; Turnock, S.; Shim, S.; Le Sager, P.; van Noije, T.; Tsigaridis, K.; et al. Significant climate benefits from near-term climate forcer mitigation in spite of aerosol reductions. Environ. Res. Lett. 2021, 16, 034010. [Google Scholar] [CrossRef]
  35. Shim, S.; Seo, J.; Kwon, S.-H.; Lee, J.-H.; Sung, H.-M.; Boo, K.-O.; Byun, Y.-H.; Lim, Y.-J.; Kim, Y.-H. Impact of future air quality in East Asia under SSP scenarios. Atmosphere 2020, 30, 439–454. [Google Scholar] [CrossRef]
  36. Li, Z.; Lau, W.K.-M.; Ramanathan, V.; Wu, G.; Ding, Y.; Manoj, M.G.; Liu, J.; Qian, Y.; Li, J.; Zhou, T.; et al. Aerosol and monsoon climate interactions over Asia. Rev. Geophys. 2016, 54, 866–929. [Google Scholar] [CrossRef]
  37. O’Neill, B.C.; Kriegler, E.; Riahi, K.; Ebi, K.L.; Hallegatte, S.; Carter, T.R.; Mathur, R.; Vuuren, D.P. A new scenario framework for climate change research: The concept of shared socioeconomic pathways. Clim. Chang. 2014, 122, 387–400. [Google Scholar] [CrossRef] [Green Version]
  38. Dellink, R.; Chateau, J.; Lanzi, E.; Magne, B. Long-term economic growth projections in the Shared Socioeconomic Pathways. Glob. Environ. Chang. 2015, 42, 200–214. [Google Scholar] [CrossRef]
  39. KC, S.; Lutz, W. The human core of the shared socioeconomic pathways: Population scenarios by age, sex and level of education for all countries to 2100. Glob. Environ. Chang. 2017, 42, 181–192. [Google Scholar] [CrossRef] [Green Version]
  40. Gidden, M.J.; Riahi, K.; Smith, S.J.; Fujimori, S.; Luderer, G.; Kriegler, E.; Vuuren, D.P.; Berg, M.; Feng, L.; Klein, D.; et al. Global emissions pathways under different socioeconomic scenarios for use in CMIP6: A dataset of harmonized emissions trajectories through the end of the century. Geosci. Model Dev. 2019, 12, 1443–1475. [Google Scholar] [CrossRef] [Green Version]
  41. O’Neil, B.C.; Tebaldi, C.; Vuuren, D.P.; Eyring, V.; Friedlingstein, P.; Hurtt, G.; Knutti, R.; Kriegler, E.; Lamarque, J.-F.; Lowe, J.; et al. The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6. Geosci. Model Dev. 2016, 9, 3461–3482. [Google Scholar] [CrossRef] [Green Version]
  42. Reis, L.A.; Drouet, L.; Dingenen, R.V.; Emmerling, J. Future global air quality indices under different socio-economic and climate assumptions. Sustainability 2018, 10, 3645. [Google Scholar] [CrossRef] [Green Version]
  43. Rao, S.; Klimont, Z.; Smith, S.J.; Dingenen, R.V.; Dentener, F.; Bouwman, L.; Riahi, K.; Amann, M.; Bodirsky, B.L.; Vuuren, D.P.; et al. Future air pollution in the Shared Socio-economic Pathways. Glob. Environ. Chang. 2017, 42, 346–358. [Google Scholar] [CrossRef]
  44. Sellar, A.A.; Jones, C.G.; Mulcahy, J.; Tang, Y.; Yool, A.; Wiltshire, A.; O’Connor, F.M.; Stringer, M.; Hill, R.; Palmieri, J.; et al. UKESM1: Description and evaluation of the UK Earth System Model. J. Adv. Model. Earth Syst. 2019, 11, 4513–4558. [Google Scholar] [CrossRef] [Green Version]
  45. Dunne, J.P.; Horowitz, L.W.; Adcroft, A.J.; Ginoux, P.; Held, I.M.; John, J.G.; Krasting, J.P.; Malyshev, S.; Naik, V.; Paulot, F.; et al. The GFDL Earth System Model version 4.1 (GFDL-ESM4.1): Model description and simulation characteristics. J. Adv. Model. Earth Syst. 2020, 12, e2019MS002014. [Google Scholar] [CrossRef]
  46. Seland, O.; Bentsen, M.; Olivie, D.; Toniazzo, T.; Gjermundsen, A.; Graff, L.S.; Debernard, J.B.; Gupta, A.K.; He, Y.-C.; Kirkevag, A.; et al. Overview of the Norwegian Earth System Model (NorESM2) and key climate response of CMIP6 DECK, historical, and scenario simulations. Geosci. Model Dev. 2020, 13, 6165–6200. [Google Scholar] [CrossRef]
  47. Bauer, S.E.; Tsigaridis, K.; Faluvegi, G.; Kelly, M.; Lo, K.K.; Miller, R.L.; Nazarenko, L.; Schmidt, G.A.; Wu, J. Historical (1850–2014) aerosol evolution and role on climate forcing using the GISS ModelE2.1 contribution to CMIP6. J. Adv. Model. Earth Syst. 2020, 12, e2019MS001978. [Google Scholar] [CrossRef]
  48. Hajima, T.; Watanabe, M.; Yamamoto, A.; Tatebe, H.; Noguchi, M.A.; Abe, M.; Ohgaito, R.; Ito, A.; Yamazaki, D.; Okajima, H.; et al. Development of the MIROC-ES2L Earth system model and the evaluation of biogeochemical processes and feedbacks. Geosci. Model Dev. 2020, 13, 2197–2244. [Google Scholar] [CrossRef]
  49. Yukimoto, S.; Kawai, H.; Koshiro, T.; Oshima, N.; Yoshida, K.; Urakawa, S.; Tsujino, H.; Deushi, M.; Tanaka, T.; Hosaka, M.; et al. The meteorological research institute Earth system model version 2.0, MRI-ESM2.0: Description and basic evaluation of the physical component. J. Meteorol. Soc. Jpn. 2019, 97, 931–965. [Google Scholar] [CrossRef] [Green Version]
  50. Emmons, L.K.; Schwantes, R.H.; Orlando, J.J.; Tyndall, G.; Kinnison, D.; Lamarque, J.; Marsh, D.; Mills, M.J.; Tilmes, S.; Bardeen, C.; et al. The Chemistry Mechanism in the Community Earth System Model version 2 (CESM2). J. Adv. Model. Earth Syst. 2020, 12, 1–21. [Google Scholar] [CrossRef] [Green Version]
  51. Wu, T.; Zhang, F.; Zhang, J.; Jie, W.; Zhang, Y.; Wu, F.; Li, L.; Yan, J.; Liu, X.; Lu, X.; et al. Beijing Climate Center Earth System Model verison 1 (BCC-ESM1): Model description and evaluation of aerosol simulations. Geosci. Model Dev. 2020, 13, 977–1005. [Google Scholar] [CrossRef] [Green Version]
  52. Tegen, I.; Neubauer, D.; Ferrachat, S.; Drian, C.S.-L.; Bey, I.; Schutgens, N.; Stier, P.; Watson-Parris, D.; Stanelle, T.; Schmidt, H.; et al. The global aerosol-climate model ECHAM6.3-HAM2.3—Part 1: Aerosol evaluation. Geosci. Model Dev. 2019, 12, 1643–1677. [Google Scholar] [CrossRef] [Green Version]
  53. Iturbide, M.; Gutierrez, J.M.; Alves, L.M.; Bedia, J.; Climadevilla, E.; Cofino, A.S.; Cerezo-Mota, R.; Di Luca, A.; Faria, S.H.; Gorodetskaya, I.; et al. An update of IPCC climate reference regions for subcontinental analysis of climate model data: Definition and aggregated datasets. Earth Syst. Sci. Data Discuss. 2020. submitted for publication. [Google Scholar] [CrossRef]
  54. Van Donkelaar, A.; Martine, R.V.; Brauer, M.; Hsu, N.C.; Kahn, R.A.; Levy, R.C.; Lyapustin, A.; Sayer, A.M.; Winker, D.M. Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellite. Environ. Sci. Technol. 2016, 50, 3762–3772. [Google Scholar] [CrossRef]
  55. Van Donkelaar, A.; Martine, R.V.; Brauer, M.; Hsu, N.C.; Kahn, R.A.; Levy, R.C.; Lyapustin, A.; Sayer, A.M.; Winker, D.M. Global Annual PM2.5 Grids from MODIS, MISR and SeaWiFS Aerosol Optical Depth (AOD) with GWR, 1998–2016; NASA Socioeconomic Data and Applications Center (SEDAC): Palisades, NY, USA, 2018. [Google Scholar] [CrossRef]
  56. Silva, R.A.; West, J.J.; Lamarque, J.F.; Shindell, D.T.; Collins, W.J.; Dalsoren, S.; Faluvegi, G.; Folberth, G.; Horowitz, L.W.; Nagashima, T.; et al. The effect of future ambient air pollution on human premature mortality to 2100 using ouput from the ACCMIP model ensemble. Atmos. Chem. Phys. 2016, 16, 9847–9862. [Google Scholar] [CrossRef] [Green Version]
  57. WHO (World Health Organization); Occupational and Environmental Health Team. WHO Air Quality Guidelines for Particulate Matter, Ozone, Nitrogen Dioxide and Sulfur Dioxide: Global Update 2005: Summary of Risk Assessment; World Health Organization: Geneva, Switzerland, 2006. [Google Scholar]
  58. Pope, C.A., III; Krewski, D.; Gapstur, S.M.; Turner, M.C.; Jerrett, M.; Burnett, R.T. Fine particulate air pollution and mortality: Response to Enstrom’s re-analysis of the American Cancer Society Cancer Prevention Study II cohort. Dose Response 2017, 15, 1–5. [Google Scholar] [CrossRef]
  59. Frieler, K.; Lange, S.; Piontek, F.; Reyer, C.P.O.; Schewe, J.; Warszawski, L.; Zhao, F.; Chini, L.; Denvil, S.; Emanuel, K.; et al. Assessing the impacts of 1.5 °C global warming—simulation protocol of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2b). Geosci. Model Dev. 2017, 10, 4321–4345. [Google Scholar] [CrossRef] [Green Version]
  60. Mitchell, J.F.B.; Johns, T.C.; Gregory, J.M.; Tett, S.F.B. Climate response to increasing levels of greenhouse gases and sulphate aerosols. Nature 1995, 376, 501–504. [Google Scholar] [CrossRef]
  61. Menon, S.; Hansen, J.; Nazarenko, L.; Luo, Y. Climate effects of black carbon aerosols in China and India. Science 2002, 297, 2250–2253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zheng, Y.; Zhang, Q.; Tong, D.; David, S.J.; Caldeira, K. Climate effects of China’s efforts to improve its air quality. Environ. Res. Lett. 2020, 15, 104052. [Google Scholar] [CrossRef]
Figure 1. Three subregional Asian domains analyzed in this study, adopted from Iturbide et al. (2020) [53]. Shaded colors indicate three major regions: East Asia (blue; EAS), South Asia (green; SAS), and Southeast Asia (orange; SEA).
Figure 1. Three subregional Asian domains analyzed in this study, adopted from Iturbide et al. (2020) [53]. Shaded colors indicate three major regions: East Asia (blue; EAS), South Asia (green; SAS), and Southeast Asia (orange; SEA).
Ijerph 18 06817 g001
Figure 2. Comparison of the domain averaged annual mean surface PM2.5 concentrations from satellite data and CMIP6 historical simulations. The black circles represent the satellite data (1998–2014) and the remaining symbols indicate data obtained from the CMIP6 historical simulations for the present day (1995–2014). The thick horizontal black lines represent the ensemble mean of the CMIP6 models. Domain regions are defined in Figure 1.
Figure 2. Comparison of the domain averaged annual mean surface PM2.5 concentrations from satellite data and CMIP6 historical simulations. The black circles represent the satellite data (1998–2014) and the remaining symbols indicate data obtained from the CMIP6 historical simulations for the present day (1995–2014). The thick horizontal black lines represent the ensemble mean of the CMIP6 models. Domain regions are defined in Figure 1.
Ijerph 18 06817 g002
Figure 3. The spatial distribution of annual mean surface PM2.5 concentrations from the CMIP6 historical simulations for the present–day period (1995–2014) and satellite data (1998–2014). (a) Multimodel mean, (b) standard deviation of the multimodel mean, (c) satellite–derived annual surface PM2.5 concentrations, and (d) the difference between multimodel mean and satellite data.
Figure 3. The spatial distribution of annual mean surface PM2.5 concentrations from the CMIP6 historical simulations for the present–day period (1995–2014) and satellite data (1998–2014). (a) Multimodel mean, (b) standard deviation of the multimodel mean, (c) satellite–derived annual surface PM2.5 concentrations, and (d) the difference between multimodel mean and satellite data.
Ijerph 18 06817 g003
Figure 4. Future changes (unit: %) in the near future (2021–2040, left-top triangle; NF) and mid-future (2041–2060, right–bottom triangle; MF) relative to the present day (1995–2014) for annual mean total emissions of (a) SO2, (b) OC, and (c) BC worldwide and for the Asian regions (EAS, SAS, and SEA) analyzed in this study under the various CMIP6 SSP scenarios.
Figure 4. Future changes (unit: %) in the near future (2021–2040, left-top triangle; NF) and mid-future (2041–2060, right–bottom triangle; MF) relative to the present day (1995–2014) for annual mean total emissions of (a) SO2, (b) OC, and (c) BC worldwide and for the Asian regions (EAS, SAS, and SEA) analyzed in this study under the various CMIP6 SSP scenarios.
Ijerph 18 06817 g004
Figure 5. Future changes in the annual surface mean PM2.5 concentrations relative to the present-day (1995–2014) mean for various CMIP6 SSP scenarios for the three Asian regions: (a) EAS, (b) SAS, and (c) SEA. Each line represents the MME mean, and the shaded areas represent a ±1 standard deviation of the mean. The MME regional mean value and ±1 standard deviation for the present day is indicated in the bottom left corner.
Figure 5. Future changes in the annual surface mean PM2.5 concentrations relative to the present-day (1995–2014) mean for various CMIP6 SSP scenarios for the three Asian regions: (a) EAS, (b) SAS, and (c) SEA. Each line represents the MME mean, and the shaded areas represent a ±1 standard deviation of the mean. The MME regional mean value and ±1 standard deviation for the present day is indicated in the bottom left corner.
Ijerph 18 06817 g005
Figure 6. Spatial distribution of AQI based on WHO air-quality guidelines and interim targets for long-term exposure of PM2.5 for (a) the present day from CMIP6 historical simulations, and for the near future from the (b) SSP1-2.6, (c) SSP2-4.5, (d) SSP3-7.0, and (e) SSP5-8.5 scenarios. AQI value classifications are defined in Table 4.
Figure 6. Spatial distribution of AQI based on WHO air-quality guidelines and interim targets for long-term exposure of PM2.5 for (a) the present day from CMIP6 historical simulations, and for the near future from the (b) SSP1-2.6, (c) SSP2-4.5, (d) SSP3-7.0, and (e) SSP5-8.5 scenarios. AQI value classifications are defined in Table 4.
Ijerph 18 06817 g006
Figure 7. Spatial distribution of AQI based on WHO air-quality guidelines and interim targets for long-term exposure of PM2.5 for the mid-future (2041-2060) from the (a) SSP1-2.6, (b) SSP2-4.5, (c) SSP3-7.0, and (d) SSP5-8.5 scenarios. AQI value classifications are defined in Table 4.
Figure 7. Spatial distribution of AQI based on WHO air-quality guidelines and interim targets for long-term exposure of PM2.5 for the mid-future (2041-2060) from the (a) SSP1-2.6, (b) SSP2-4.5, (c) SSP3-7.0, and (d) SSP5-8.5 scenarios. AQI value classifications are defined in Table 4.
Ijerph 18 06817 g007
Figure 8. Time variation of total population (unit: million, solid lines) and ratio over 65 years to total population (unit: %, dotted lines) for the different SSP scenarios across (a) China and (b) India. Datasets were obtained from the International Institute for Applied System Analysis (IIASA) [39].
Figure 8. Time variation of total population (unit: million, solid lines) and ratio over 65 years to total population (unit: %, dotted lines) for the different SSP scenarios across (a) China and (b) India. Datasets were obtained from the International Institute for Applied System Analysis (IIASA) [39].
Ijerph 18 06817 g008
Figure 9. Scatter plot and regression line of simulated annual mean surface PM2.5 concentrations and clear-sky surface downwelling shortwave radiative fluxes from 1950 to 2014 for the Asian regions. The closed circles represent the ensemble mean values for the present-day period (1995–2014).
Figure 9. Scatter plot and regression line of simulated annual mean surface PM2.5 concentrations and clear-sky surface downwelling shortwave radiative fluxes from 1950 to 2014 for the Asian regions. The closed circles represent the ensemble mean values for the present-day period (1995–2014).
Ijerph 18 06817 g009
Figure 10. Future changes in annual mean total aerosol optical depths (ac), clear-sky surface downwelling shortwave fluxes (df), and surface air temperatures (gi) relative to the present-day (1995–2014) mean for the SSP3-7.0 and SSP3-7.0-lowNTCF scenarios for the three Asian regions (EAS, SAS, and SEA). Each line represents the MME mean and the shading represents the ±1 standard deviation of the mean. The MME regional mean value and ±1 standard deviation for the present day is indicated near the horizontal lines.
Figure 10. Future changes in annual mean total aerosol optical depths (ac), clear-sky surface downwelling shortwave fluxes (df), and surface air temperatures (gi) relative to the present-day (1995–2014) mean for the SSP3-7.0 and SSP3-7.0-lowNTCF scenarios for the three Asian regions (EAS, SAS, and SEA). Each line represents the MME mean and the shading represents the ±1 standard deviation of the mean. The MME regional mean value and ±1 standard deviation for the present day is indicated near the horizontal lines.
Ijerph 18 06817 g010
Table 1. Shared socioeconomic pathways’ (SSPs) air pollution control levels (Based on [42]).
Table 1. Shared socioeconomic pathways’ (SSPs) air pollution control levels (Based on [42]).
SSP ScenariosEmission Factors
SSP1 and SSP5Strong decrease (fastest and widest implementation of air pollution controls)
SSP2Medium decrease (significant advancement in pollution control, yet less than in SSP1 and SSP5)
SSP3 and SSP4Weak decrease (slowest deployment of air pollution controls)
Table 2. List of CMIP6 climate models used for historical and future scenario experiments in this study. Circles indicate that the model was used, but asterisks indicate a lack of surface PM2.5 concentration data. Detailed descriptions of the CMIP6 models are provided in Supplementary Table S1.
Table 2. List of CMIP6 climate models used for historical and future scenario experiments in this study. Circles indicate that the model was used, but asterisks indicate a lack of surface PM2.5 concentration data. Detailed descriptions of the CMIP6 models are provided in Supplementary Table S1.
Model NameHistoricalSSP1–2.6SSP2–4.5SSP3–7.0SSP3–7.0-lowNTCFSSP5–8.5
UKESM1-0-LL [44]
GFDL-ESM4 [45]
NorESM2-LM [46]
GISS-E2-1-G [47]
MIROC-ES2L [48]*
MRI-ESM2-0 [49] *
CESM2-WACCM [50]
BCC-ESM1 [51]
MPI-ESM1.2-HAM [52]
Total number of models955775
Table 3. List of CMIP6 experiments used in this study.
Table 3. List of CMIP6 experiments used in this study.
ExperimentInformation
Historical
(1850–2014)
The historical simulations use forcing due to both the natural causes and human factors over the period 1850 to 2014. These simulations were used to evaluate model performance.
SSP1-2.6
(2015–2100)
This scenario represents the low end of the range of plausible future pathways, and depicts the best-case future scenario from a sustainability perspective.
SSP2-4.5
(2015–2100)
This scenario represents the medium part of the range of plausible future pathways.
SSP3-7.0
(2015–2100)
This scenario represents the medium to high end of plausible future pathways.
SSP3-7.0-lowNTCF
(2015–2055)
This scenario represents the SSP3-7.0 scenario with the reduced near-term climate forcer (NTCF) emissions, including aerosols.
SSP5-8.5
(2015–2100)
This scenario represents the high end of plausible future pathways. SSP5 is the only SSP scenario with emissions high enough to produce the 8.5 Wm−2 level of forcing in the year 2100.
Table 4. List of CMIP6 experiments used in this study.
Table 4. List of CMIP6 experiments used in this study.
IndexPM2.5 (µg/m3)Basis for the Selected Level [57]
5Significantly over target (ST)53–Defined as a concentration that exceeds 150% of the interim target-1 level.
4Over target (OT)35–53Defined as a concentration higher than the interim target and less than 150% of the interim target-1 level.
3Interim target 1 (IT-1)25–35Approximately 15% higher long-term mortality risk relative to the air-quality guideline level.
2Interim target 2 (IT-2)15–25These levels lower the risk of premature mortality by approximately 6% relative to the IT-1 level.
1Interim target 3 (IT-3)10–15These levels reduce the mortality risk by approximately 6% relative to the IT-2 level.
0Air-Quality Guideline (AQG)0–10Lower end of the range of significant effects on survival in response to long-term exposure to PM2.5. [58]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shim, S.; Sung, H.; Kwon, S.; Kim, J.; Lee, J.; Sun, M.; Song, J.; Ha, J.; Byun, Y.; Kim, Y.; et al. Regional Features of Long-Term Exposure to PM2.5 Air Quality over Asia under SSP Scenarios Based on CMIP6 Models. Int. J. Environ. Res. Public Health 2021, 18, 6817. https://doi.org/10.3390/ijerph18136817

AMA Style

Shim S, Sung H, Kwon S, Kim J, Lee J, Sun M, Song J, Ha J, Byun Y, Kim Y, et al. Regional Features of Long-Term Exposure to PM2.5 Air Quality over Asia under SSP Scenarios Based on CMIP6 Models. International Journal of Environmental Research and Public Health. 2021; 18(13):6817. https://doi.org/10.3390/ijerph18136817

Chicago/Turabian Style

Shim, Sungbo, Hyunmin Sung, Sanghoon Kwon, Jisun Kim, Jaehee Lee, Minah Sun, Jaeyoung Song, Jongchul Ha, Younghwa Byun, Yeonhee Kim, and et al. 2021. "Regional Features of Long-Term Exposure to PM2.5 Air Quality over Asia under SSP Scenarios Based on CMIP6 Models" International Journal of Environmental Research and Public Health 18, no. 13: 6817. https://doi.org/10.3390/ijerph18136817

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop