A development of reduction scenarios of the short-lived climate pollutants (SLCPs) for mitigating global warming and environmental problems

This study presents the results of the ERTDF S-12 project for searching an optimum reduction scenario of the short-lived climate pollutants (SLCPs) to simultaneously mitigate the global warming and environmental problems. The study utilized REAS emission inventory, Asia-Pacific Integrated Model-Enduse (AIM/Enduse), MIROC6 climate model, NICAM non-hydrostatic atmospheric model, and models for estimating environmental damages to health, agriculture, and flood risks. Results of various scenario search indicate that it is difficult to attain simultaneous reduction of global warming and environmental damages, unless a significant reduction of CO2 is combined with carefully designed SLCP reductions for CH4, SO2, black carbon (BC), NOx, CO, and VOCs. In this scenario design, it is important to take into account the impact of small BC reduction to the surface air temperature and complex atmospheric chemical interactions such as negative feedback between CH4 and NOx reduction. We identified two scenarios, i.e., B2a and B1c scenarios which combine the 2D-scenario with SLCP mitigation measures using End-of-Pipe (EoP) and new mitigation technologies, as promising to simultaneously mitigate the temperature rise by about 0.33 °C by 2050 and air pollution in most of the globe for reducing damages in health, agriculture, and flood risk. In Asia and other heavy air pollution areas, health-care measures have to be enhanced in order to suppress the mortality increase due to high temperature in hot spot areas caused by a significant cut of particulate matter. For this situation, the B1b scenario is better to reduce hot spot areas and high-temperature damage to the public health.


Introduction
The concept of the short-lived climate pollutants (SLCPs) has been created in the development of global warming mitigation policies to identify anthropogenic atmospheric compositions that produce positive radiative forcing. Black carbon (BC), methane, and tropospheric ozone are key constituents among the SLCPs. UNEP and WMO (2011) proposed a scenario of future global warming mitigation by which a surface air temperature decrease of about 0.5°C can be attained by the introduction of SLCP reduction measures along with that of long-lived greenhouse gases (LLGHGs). An international initiative of the Climate and Clean Air Coalition (CCAC) was launched in 2012 with a scope of simultaneous mitigation of global warming and health problem through the reduction of SLCPs and wastes. Such a combined mitigation strategy of LLGHG and SLCP has become more important after the Paris Agreement of COP21 in 2015 set a goal of the global surface air temperature not exceeding 2°C and an effort goal of 1.5°C relative to the pre-industrial era. According to the IPCC special report on global warming of 1.5°C (IPCC 2018), these global warming mitigation goals require a significant reduction of LLGHG emission at a level of almost zero by 2050, which demands the maximum effort of the society. To meet the goals, SLCP reduction is useful to contribute to the mitigation effort, particularly given that a short lifetime of SLCP assures a quick reduction of surface temperature on the order of 0.5°C in less than 10 years right after the action of SLCP reduction made. On the other hand, the temperature reduction by LLGHGs takes several 10 years due to their long lifetime in the atmosphere.
The design of effective mitigation measures critically depends on our understanding of the complex climate effects of the short-lived atmospheric compositions including SLCPs and their precursors. This motivated various research efforts in the past taken to challenge the problem, e.g., UNEP ABC project (Ramanathan and Crutzen 2003;Nakajima et al. 2007), AeroCom (Schulz et al. 2006) and AeroCom-II , Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) (Shindell et al. 2013), European Union Seventh Framework Programme project-Evaluating the Climate, and Air Quality Impacts of Short-Lived Pollutants (ECLIPSE, Stohl et al. 2015), Asia Pacific Clean Air Partnership (APCAP) (Akimoto et al. 2015;UNEP 2018), NAPEX PDRMIP ), among others. Ministry of Environment of Japan carried out two SLCP-related strategic researches, i.e., S-7 and S-12 by the Environment Research and Technology Development Fund (ERTDF).
This paper summarizes the key results of the S-12 project that aimed at searching an optimum mitigation path of SLCP reduction, conducted in Japanese fiscal years of 2014-2018, through four research theme activities, i.e., (1) analysis of air quality change and development of emission inventories, (2) selection of mitigation technologies and scenario development, (3) investigation of SLCP impacts on the climate and environment, and (4) development for shared tools and a high-resolution atmospheric modeling to link a large variety of scales from regional to global. In this paper, depending on the context, we use the term SLCP to simply identify anthropogenic short-lived atmospheric constituents including scattering aerosols and their precursors as well as absorbing aerosols.

Methods/experimental
Evaluation of the SLCP effects on the surface air temperature One particular issue in designing the effective SLCP reduction scenario is the fact that anthropogenic atmospheric compositions of different species have different signs of the radiative forcing (abbreviated as RF, hereafter). Figure 1 shows direct RFs in China and Asia in a period of 1980-2010 and RF changes caused by 50% reduction of all the emissions. These are obtained from climate change simulations with the MIROC5/MIROC6 climate models (Watanabe et al. 2010;Tatebe et al. 2019), implemented with SPRINTARS aerosol model (Takemura et al. 2005) and CHASER atmospheric chemistry model (Sudo and Akimoto 2007). The figure indicates that a simple 50% reduction of the air pollutants produces a net positive forcing of + 10 mW m − 2 as a result of the cancelation of a negative forcing due to reduction in black carbon (BC), tropospheric ozone (TO 3 , hereafter), and methane (CH 4 ) by a positive forcing due to sulfate reduction, suggesting that a simple 50% reduction is not an optimal scenario to simultaneously mitigate the global warming and air pollution problem.
Another difficulty that complicates the scenario design is the presence of indirect effects between SLCPs and the climate system. Several key indirect effects have been identified in the past studies regarding sulfate, BC, and methane. Figure 2 shows contributions of methane and TO 3 to the net RFs as functions of emission change factors of nitrogen oxides (NO x ), carbon monoxide (CO), and volatile organic compounds (VOCs) including indirect RFs caused by feedbacks among atmospheric compositions. Reference concentrations for the emission factor were assumed to be the 2008 values of the HTAP2 multi-model study (Stjern et al. 2016). The figure indicates that a decrease in NO x emission produces a negative RF caused by indirect reduction of TO 3 and a positive RF due to indirect change of methane, resulting in a net positive RF that acts to warm the globe. This enhanced methane increase is caused by a prolonged  lifetime of methane due to the reduction of OH radical (Shindell et al. 2009;Akimoto et al. 2015). On the other hand, CO emission reduction can produce a negative forcing of methane. This experiment suggests that a careful design of reduced rates of NO x and CO is required to attain a reduction of TO 3 without increasing methane concentration. Another process of large uncertainty is the indirect climate effects of sulfate and BC aerosols, which exert a significant impact on the SLCP-induced climate response. Figure 3, cited from Takemura and Suzuki (2019), shows simulated equilibrium global surface air temperature changes as a function of the instantaneous direct RF, F i , of sulfate, and BC applied to the tropopause. The figure indicates that a sulfate reduction generates a significant temperature rise in a manner mostly proportional to the sulfate loading, whereas a BC reduction produces a much smaller temperature decrease with the sensitivity or the slope far smaller than that of sulfate. The large climate sensitivity of sulfate is considered to result from the superposition of direct RF and indirect RF due to induced cloud RF of same signs, whereas the small climate sensitivity of BC is caused by the cancelation of direct RF by indirect RF as also found by previous studies (Hansen et al. 2005;Hodnebrog et al. 2014;Stohl et al. 2015;Samset et al. 2016). Tables 1 and 2 and Fig. 4 compare the values of reported RFs and associated surface air temperature changes caused by CO 2 and SLCPs. Here, following Hansen et al. (2005), we define the instantaneous direct RF as the earth radiation budget perturbation at tropopause for RF computation, or simply referred to as the Top Of the Atmosphere (TOA), that occurs immediately when a climate change factor is added to the system, and the effective RF (ERF), F e , is defined as the net radiation imbalance with the Sea Surface Temperature (SST) fixed. The latter contains the rapid adjustment of the climate system and is considered to be a direct driver of the slow climate response to cause the global-mean surface air temperature change ΔT s . In this definition, the first and other aerosol indirect effects and semi-direct effects are included in the rapid adjustment. We performed 150-year run for equilibrium climate response using MIROC6 with ocean coupled to estimate ΔT s and feedback parameter β defined as, where the radiative imbalance at TOA (denoted by F) and T s are represented to depend on time (t). The β value thus defined is known to be similar to that derived from the regression analysis by the method of Gregory (2004) without a large dependence of β on the time period of the analysis. We therefore selected the equilibrium climate sensitivity parameter, λ = 1/β, and the forcing ratio, as useful parameters in our analysis of the SLCP impact for comparison among present and past studies that adopt different methods to evaluate F e with different time scales to separate rapid and slow climate responses. We compare the reported values of γ and λ in Fig. 4. Figure 4a shows that the magnitude of F e for sulfate is γ = 1.7 to 3 times larger than that of F i . The γ value for the total aerosol is in a range of γ = 2 to 4.5 and is noticeably larger than that for sulfate. This phenomenon can be explained by the fact that the magnitude of F i for the total aerosol is smaller than that for sulfate due to light absorption by absorbing materials externally and internally mixed in the total aerosol. The major contribution to the large γ value for sulfate is caused by the first indirect effect of sulfate aerosol acting as effective cloud condensation nuclei (CCN) that enhances lowlevel cloud albedo. On the other hand, the magnitude of F e of BC is significantly smaller by 25 to 75% than the magnitude of F i with an exception of GISS-E2-R model which produces a forcing larger than F i . Recent studies proposed that the low climate sensitivity of BC occurs as a result of the cancelation of negative RF due to BC reduction by a positive RF arising from the cloud system response generated by decreased atmospheric stability due to reduced heating by BC Suzuki and Takemura 2019). There is also a RF due to decreased semi-direct effect associated with decreasing solar radiation absorption by decreased BC. In this regard, it is important to recognize the existence of a large inter-model spread of the ERF induced by BC and the fact that the MIROC model, used in this study, belongs to the group of small ERF for BC (Shindell et al. 2013). This small γ value by the MIROC model is, however, similar to or even larger than the result of 14 km global simulations by NICAM (Tomita and Satoh 2004;Goto et al. 2020) implemented with more first-principle modeling of cloud formation and aerosol-cloud interaction processes. On the other hand, the γ value exceeding the unity in the GISS-ES-R model case also suggests that the BC-induced cloud response might be more complicated to cause cloud forcing that could amplify, rather than dampen, the initial forcing F i (Shindell et al. 2013).   Future studies are needed to conclude this issue, especially in the atmospheric moist convective adjustment process. Figure 4b shows that the λ value ranges from 0.6 to 0.8 for sulfate similar to that of doubling CO 2 ranging from 0.5 to 0.9 with the model spread smaller for sulfate than for CO 2 . On the other hand, the slow response to BC is large model dependent from 0.1 to 0.8. Although it is early to lead a conclusion before more future model comparisons, the figure suggests that CO 2 , sulfate and BC have different effective forcing sensitivities to the slow climate response with different forcing mechanisms, i.e., thermal infrared forcing by CO 2 , solar forcing to the earth's surface by sulfate, and solar absorption by BC in the atmosphere. Sulfate causes a unique large cooling of the surface through direct effect and aerosolcloud interaction to increase the cloud optical thickness. MIROC has a large value of the forcing ratio γ and small value of λ, so that the resulted temperature rise per instantaneous forcing is similar to those models with small γ and large λ values.
There are several issues to improve the large model uncertainties in the RF evaluation used in the model comparisons shown in Fig. 4. A well-known problem is a negative model bias of AOT and a part of which is attributed to the insufficient model treatment of missing nitrate and secondary organic aerosol (Shindell et al. 2013;Shrivastava et al. 2017). Another error source is the inter-model difference in the representation of the mixed state of BC with other atmospheric compositions Oshima and Koike 2013;Matsui et al. 2018) that significantly influences the RF estimate. Also, some of absorption by colored pollution aerosols, such as brown carbon (Kirchstetter et al. 2004) and combustion iron (Moteki et al. 2017), is approximated by BC absorption in a simplified model treatment. The BC vertical profile is also a key factor affecting the RF and ERF estimates and could be constrained by aircraft measurements through calibrating the wet removal process as claimed by Hodnebrog et al. (2014). Undefined BC sources in the common emission inventories are still a problem, for example, in India (e.g., Goto et al. 2011).
Additional numerical experiments would also be useful to better quantify regional responses to SLCP perturbations, including a shift of the tropical precipitation belt that occurs via both atmosphere and ocean energy transports (Zhao and Suzuki 2019) and modulations of the hydrological cycle over the tropical Asian monsoon regions (Takahashi et al. 2018). High-resolution modeling is also important. Sato et al. (2016) reported an increase in BC transportation to a polar region with their highresolution modeling, and there is also a recent argument regarding the anti-Twomy response of low clouds as revealed by active satellite remote sensing and high-resolution modeling (Michibata et al. 2016;Sato et al. 2018;Sato and Suzuki 2019). These studies suggest that we need high-resolution modeling that better represents the aerosol and cloud processes for more reliable estimates of RF and climate sensitivity due to the aerosolcloud interaction phenomenon in the scenario impact study. Figure 5 schematically illustrates various SLCP-related processes affecting the global air temperature and environmental problems (water resource, flooding, draught; health and agricultural damage). In the figure, solid and dashed lines respectively represent positive and negative interactions to produce increasing and decreasing outputs with increasing input. The thickness of the line indicates the strength of the process. As discussed in Figs. 3 and 4, a decrease in sulfate causes an increase in the global surface air temperature by both direct and cloud RF-related indirect effects, but at the same time, it decreases the public health damage, implying a need for a balanced cut of sulfate for reducing the health problem without accelerating the global warming. In contrast, a BC decrease causes a small net increase in the surface air temperature as a balance of positive direct and negative cloud RF-related indirect effects. It also decreases the environmental problems, so that a substantial cut of BC is desirable for simultaneous mitigation of global warming and environmental problems. Furthermore, a NO x decrease causes a decrease in TO 3 , but at the same time, it causes a negative indirect effect to increase methane as discussed in Fig. 2. In order to weaken the negative effect, we can combine CO and VOC reductions as a mitigation measure that cancel the NO x negative feedback and also serve to reduce TO 3 .

Construction of SLCP scenarios
The S-12 project adopted the REAS 2.1 emission inventory (Kurokawa et al. 2013) and the Asia-Pacific Integrated Model-Enduse (AIM/Enduse) to relate the mitigation technologies with emissions of SLCPs. Table 3 shows the codes and names of SLCP scenarios constructed by the present project using AIM/Enduse. Details of scenario design are given in Hanaoka and Masui (2020). For this work, the AIM/Enduse has been implemented with several SCLP processes to get a better consistency with the REAS inventory and to output scenario parameters for climate simulations. The reference scenario (Ref-scenario) is a modified version of the SSP2 scenario (O'Neill et al. 2015), in which the future mitigation policies and technologies are taken place with the current trends. The medium and maximum End-of-Pipe (A1 and A2 scenarios) scenarios are scenarios of enhancing EoP-diffusion both in developed and developing countries by 2050 for reducing SO 2 , NO x , BC, OC, PM2.5, and PM10. The 2D-scenario is designed for significant emission reductions of CO 2 and methane, by utilizing decarbonization mitigation measures toward the 2°C target of the Paris Agreement with reaching a carbon price at 400 US$/tCO 2 in 2050. The CCSscenario is one of 2D-scenarios with energy shift to coal and biomass power and CCS technology rather than using renewable energies. RES-and BLD-scenarios are 2D-scenarios with energy shift to renewables and with enhancing electrification in the building sector across the world, respectively. The TRT-scenario is for enhancing EV in the passenger transport sector across the world by 2050. We mainly use the abbreviated scenario codes in this paper, rather than full scenario names, for concise presentation.
Resulted time series of global emissions from 2010 to 2050 are shown in Fig. 6 and their relations in Fig. 7 between (a) CO 2 -SO 2 , (b) BC-SO 2 , (c) CO-NO x , and (d) VOC-NO x . An appendix also presents global spatial distributions of SLCP emission changes. First of all, we note that the simulated emissions at 2010 by the improved AIM model are consistent with those of EDGAR4.3 and HTAP without large gaps in the current condition at 2010. Secondly, Fig. 6 shows that Ref, A2, and A1 scenarios result in 43% increase of CO 2 emission and 53% increase in methane emission from 2010 to 2050. On the other hand, the introduction of 2D-scenarios reduces CO 2 emission by 50% and methane emission by 44% from 2010 to 2050. SLCP emission changes depend on how new mitigation technologies are selected among    Fig. 6 does not have much variety in the scenarios because the major sources of methane are from agriculture and waste sectors for which we do not pose many options for the new mitigation technologies. Future studies of larger emission cuts are needed to be sought by the innovation of new methane reduction technologies. A study of Fig. 7 for characteristic patterns of relations among SLCP compositions suggest B1a and B1c scenarios are advantageous for large BC cut and medium-cut of SO 2 , while B1c scenario is useful for a large cut of CO and VOC and medium-cut for NO x . These figures show that there is a large room for bringing different future air compositions by adopting different combinations of the new mitigation technologies. This means that the present study serves as a testbed for the search of the scenarios with results representing a part of the future full search which needs more simulations with climate and environmental models. It is also noted that the cost of introducing new mitigation technologies is larger than that of the EoP technology, that implies a careful discussion considering cost-benefit and trade-offs has to be included in the future scenario search. Figure 8 shows atmospheric concentrations of TO 3 in the Dobson Unit (DU) simulated with the constructed SLCP scenarios. The figure shows that the TO 3 concentration will increase in the Ref-scenario with increasing NO x as indicated in Fig. 6. It is interesting to find that A1 and A2 scenarios increase TO 3 even though NO x is decreased in these scenarios. This phenomenon is understood by increases in CO and VOC in these scenarios as shown in Fig. 6. Scenarios of new mitigation technologies are successful to reduce TO 3 , especially by 2.6 DU from 2010s to 2040s with B2a and B1c scenarios. This situation is similar to the indirect effect to the methane concentration as shown in Fig. 9. The methane emission is prescribed in the S-12 scenarios as given in Fig. 6, but the atmospheric concentration of methane is strongly affected by the indirect effect of SLCP compositions as discussed in the preceding section. Dotted lines in Fig. 9 indicate that the methane concentration will increase by 0.2 ppmV for the Ref-scenario more than simulated without the feedback process as presented by solid lines. We found more increase of methane by 0.3 ppmV for A2 scenario caused by a large decrease in OH radical with decreasing NO x in the scenario. On the other hand, the figure shows that the reduction of CO and VOC in the new technology scenarios successfully suppresses the methane increase due to NO x reduction. Figure 10 shows simulated aerosol optical thickness (AOT) and PM2.5 in the globe and East Asia (20°N-50°N, 100°E-150°E). The figure indicates that the A1 scenario does not improve AOT and PM2.5, but the A2 scenario can reduce AOT by 0.01 and PM2.5 by 0.3 μg m − 3 in the globe from the current state. It is important to find that the B1c scenario can perform good mitigation similar to B2a, i.e., AOT by 0.02 and PM2.5 by 0.9 μg m − 3 from the current state, even when the A1 scenario is adopted but combined with the new mitigation technologies of RES, BLD, and TRT. The East Asian region is more polluted than the global average, so that the role of EoP technology is more important. The Ref-scenario can maintain the present values of AOT and PM2.5, and A2 scenarios will decrease AOT by 0.14 and PM2.5 by 5.3 μg m − 3 , respectively, from 2010. And the B1c scenario can decrease AOT by 0.17 and PM2.5 by 7.3 μg m − 3 .
The analyses in this section are especially important to study the role of the new mitigation technologies to ease the problem that EoP-scenarios (A1 and A2) is beneficial for improving health problem but increases global warming. One conclusion from Figs. 8, 9, 10 is that a pollution control measure to introduce the maximum EoP-technology alone will not lead us to clean and warming mitigation goals, and therefore, we need to

SLCP impacts on the earth's climate and environment
Emission scenarios constructed by the AIM/Enduse were used for the MIROC6 climate model to perform seven ensemble runs of coupled atmosphere-ocean model simulation for each scenario with a grid resolution of 1.4°(T85) for standard simulations and with 0.56°resolution for high-resolution simulations. In detail, atmospheric concentrations of CO 2 and methane were simulated by the MAGICC6 model (Meinshausen et al.   Figure 11 shows the simulated time series of the global mean surface air temperature changes relative to Refscenario. The figure indicates that Ref-scenario will increase the global mean surface air temperature by about 1.0°C by 2050. The introduction of EoP-technology (A2 and A1 scenarios) alone increases temperatures more than Ref-scenario about 0.07°C with a large positive RF generated by a decrease of sulfate. This increase will be enhanced if we introduce the indirect effect of methane as discussed in Fig. 9. On the other hand, the introduction of 2D-measures and new mitigation technologies will produce a global mean temperature decrease of about 0.3°C from Ref-scenario by 2050. The best scenarios for warming mitigation are B1a and B1c which will decrease the temperature by 0.31°C and 0.33°C by 2050, respectively. These scenarios select the A1 scenario with less reduction of sulfate than A2 scenario indicating an effectiveness of a combination of technologies that simultaneously decreases SO 2 , BC, methane, NO x , CO, and VOC in a suitable amount so as not to produce a large warming by sulfate. Figures 12 presents global  other than hot spot areas in India, Myanmar, and China, where the scenarios cannot avoid a temperature increase reaching 0.3°C due to a large positive RF by sulfate reduction. Figure 14 shows global annual excess mortalities due to PM2.5 and high temperature for Ref2010 and seven scenarios in the 2040s. The excess mortality due to PM2.5 was estimated using the health impact assessment framework of Seposo et al. (2019) as follows: where ΔM PM2.5 is the excess mortality due to PM2.5, a product of population attributable fraction (PAF PM2.5 ) and the baseline mortality (N) for each grid; N was calculated by multiplying the population (Pop) by the baseline mortality rate (bmort); PAF PM2.5 is the proportion of the mortality that is attributable to PM2.5 and was calculated by using relative risks (RR PM2.5 ), where Z q is the scenario-specific PM2.5 concentration (annual mean) with varying theoretical minima Z 0 relative to the health outcomes, and an integrated exposureresponse function parameter of α, β, γ (Burnett et al. 2014). In this study, we estimated RRs for stroke, ischemic heart disease, chronic obstructive pulmonary disease, lung cancer, and acute respiratory infection. We applied the same framework in estimating the excess mortality due to a high temperature. To derive the heatrelated risks (RR temp ) above the optimum temperature (T op ), we applied a V-shaped temperature-mortality risk function as, where δ is the parameter of the risk function above T op ; T q is the scenario-specific daily maximum temperature. Based on the previous study (Honda et al. 2014), T op was defined as 84th percentile value of daily maximum temperature for each grid. We summed the daily excess mortalities through the year. We assumed no adaptation measures in the estimation. The figure indicates Ref and A1 scenarios cannot reduce the mortality due to PM2.5 from the situation in the 2010s, because the expected future population increase in polluted areas overcomes the air pollution mitigation effort in these scenarios. We need new mitigation technologies for a substantial decrease in the mortality. The introduction of B2a and B1c scenarios will prevent the mortality by 1.2 to 2 million persons per year in the 2040s relative to that of Ref2010.
On the other hand, although the magnitude is not as large as that of the PM2.5 impact, the annual mortality due to high temperature will increase from the present with all the scenarios, i.e., 0.022 million persons with B1b scenario to 0.037 million persons with Ref-scenario in the 2040s depending on the global mean temperature rises of 0.6°C and 1.0°C, respectively, as shown in Fig. 11. It is important to note that the mortality also increases in hot spot areas of temperature rise which remain in polluted areas, e.g., India, Myanmar, and China as shown in Fig. 13, caused by a large positive RF by sulfate reduction. This is the reason why B2a and B1c scenarios sit next to the B1b scenario for mitigation of the mortality, as understood by hot spots higher than those of B1b. An alternative solution of trade-off is to adopt B2a or B1c scenario for mitigation of health damage due to PM2.5 with enhanced anti-high temperature health care in the heavy air pollution areas in Asia where are expected to experience a hot spot-type temperature increase due to significant cut of particulate matter. Figure 15 presents changes in the rice production due to climate change and ozone estimated by a crop growth simulation model, MATCRO (Masutomi et al. 2016a(Masutomi et al. , 2016b, and a flux-based model for ozone damage on rice (Mills et al. 2018). Annual losses of the rice yield On the other hand, the climate change decreases the rice production by 0.21 ton ha − 1 to 0.80 ton ha − 1 with all the scenarios depending on the temperature rise. Different from the case of health impact due to high temperature, however, the least production loss is realized by B2a and B1c, not B1b, because the rice production is less affected by the hot spot in the polluted areas. Figure 16 shows the global population change under risk of flooding and precipitation change when BC and SO 2 concentration change from 0.5 to 10 times less or more than present (Yoshimura et al. 2018). Though the SLCP impacts on globally averaged land precipitation are rather straightforward as decreasing with increasing sulfate or BC, the impact on peak river discharge (i.e., flooding) is not simple because of various impact factors regarding hydrological processes over land including runoff, evaporation, snowmelt, river routing, etc. The figure indicates a V-shaped dependence of the risk with a minimum around the present state, suggesting that the human society has been adapted so as to minimize the disaster risk in the present climate condition. This result indicates the disaster prevention measures in the decreasing global mean precipitation should also be aware of. Although the magnitude of BC is small and noisy, we found that there is an increasing tendency of flood damage, partly because the precipitation increases in some parts of the land area. Figure 17 summarizes changes of CO 2 , SLCPs, PM2.5, AOT, and environmental parameters relative to Refscenario in the 2040s. The present study concludes that it is possible to create a CO 2 and SLCP reduction scenario that simultaneously mitigates global warming and damages on health, agriculture, and water risk by a welldesigned combination of various mitigation technologies as discussed in Figs. 12,13,14,15,and 16. Among the seven scenarios developed in this study, B2a and B1c scenarios are the best to simultaneously mitigate temperature rise, mortality due to PM2.5, and loss of rice production by ozone and climate change, whereas B1b scenario is the best to mitigate the mortality increase due to high temperature because the scenario eases high temperatures in the hot spots caused by positive RF due to sulfate decrease in the polluted area in Asia. Figure 18 summarizes changes of ERFs of CO 2 and SLCPs relative to the Ref-scenario in the 2040s. The figure indicates that the main contributors are emission cut of CO 2 and CH 4 by the 2D-scenario producing RFs of about − 0.87 Wm − 2 . The introduction of the EoPtechnology for air quality control will produce a small positive RFs of about + 0.05 Wm − 2 and will not help mitigation of the global warming. Furthermore, the introduction of new mitigation technologies of CCS, Fig. 16 Relationship between the proportion of the global population exposed to flood risk and changes in terrestrial precipitation due to BC (orange) and sulfate (blue) emissions changes (0.5 to 10 times, respectively). Error bars represent the standard deviation of interannual variation during the 20-year simulations Nakajima et al. Progress in Earth and Planetary Science (2020)  RES, BLD, and TRT reduces sulfate for improving air quality more than EoP-scenarios, and it produces positive ERF more than + 1 Wm -2 . The figure indicates, however, the advantage of the new mitigation technologies is to compensate for the positive sulfate ERF by negative ERFs by TO 3 and BC reductions to result in a larger net reduction of the total ERF about − 1.4 Wm -2 . The figure indicates the best scenario of reducing the total ERF is the B1c scenario with − 0.7 Wm -2 better than that of B2b. Also, it is important to note that the small value of BC ERF does not mean small climate impacts of BC, because of the reduced amount of ERF from instantaneous RF, F i /F e shown in Fig. 4a, is used to induce changes in the fast component of the climate system including global cloud and precipitation changes (Suzuki and Takemura 2019). Although these results suggest recommendable scenarios for simultaneous mitigation of global warming and environmental problems, we feel that the knowledge obtained by the present study is still limited to find the ultimate best scenario, and several undone works are left for future works. First, we should continue the effort of evaluating BC forcing because estimates of γ and λ values are largely scattered as shown in Fig. 4. Indirect effects on methane and TO 3 are also important processes in the scenario construction, because they increase scenario differences between EoP-scenarios and 2D-scenarios than presented in Figs. 17 and 18. Furthermore, emission change scenarios in Figs. 6 and 7 suggest that there is still a large room for bringing different future air compositions by adopting different combinations of the new mitigation technologies, especially for those of methane. The cost of introducing new mitigation technologies is larger than that of EoP technology, so that a careful discussion considering cost-benefit and the trade-off has to be included in the scenario search.

Conclusions
As found by Fig. 4, large model uncertainties of RFs still exist and can be an Achilles tendon for reliable scenario construction. Future studies are, therefore, needed for improving the model treatment of missing nitrate and secondary organic aerosol, representation of the mixed state of BC, and colored aerosols, BC vertical profile, and undefined BC sources in Asia, and other areas. Also, future detailed ensemble simulations with coupled atmosphere-ocean models and multi-model experiments are needed to investigate the BC impact in the situation of large model dependency as described in Section 2. Such additional experiments would also be useful to better quantify regional responses to SLCP perturbations.
Another important challenge for the effective SLCP mitigation is to cycle the scenario construction and their assessment by monitoring the air quality and climate state. One of significant problems met in the S-12 project is the long turnaround time of emission inventory construction by the bottom-up method based on various census data, which became a bottleneck in the project for understanding the emission gap between emission inventories and constructed scenarios. In order to ease this problem, we introduced NO x emissions in the Asian region retrieved from the satellite data (Yumimoto et al. 2015). Fig. 19 compares estimates of the annual NO x emission from the Chinese area using updated versions of REAS 2.2 and satellite retrieval (Kurokawa et al. 2017;Itahashi et al. 2019). The figure shows the emission inventory significantly exceeds the satellite-based estimate by about 2 Mton/year after 2011. We found the difference is caused by problems both in the inventory and the satellite methods. As for the bottom-up emission inventory REAS, introduction rates of denitrification equipment to large power plants were underestimated. On the other hand, the rapid increase in Chinese emissions caused negative biases in the satellite values. To overcome this problem, we developed a sequential update technique (Itahashi et al. 2019). The revised inventory and satellite values get closer to each other as illustrated in Fig. 19. This update helped reducing the gap between REAS emission and AIM/Enduse scenario in Asia as reported by Hanaoka and Masui (2020). This experience clearly indicates shortening the turn-around time of the scenario construction and air quality monitoring cycle is useful to healthy SLCP mitigation measures along with the global stock-take initiatives by UNFCCC.
Additional file 1: Figure A1 shows global spatial distributions of SLCP emissions in the 2010s (2010-2019) with the Ref-scenario (Ref2010) and emission changes with the seven scenarios in the 2040s (2040-2049) relative to Ref2010. It is found from the figure that the reduction rates of SLCP compositions depend on regions, so that large regional differences in the mixing state of SLCPs are expected resulting in different magnitudes and signs of RF. BLD-and TRT-technologies are effective to reduce BC and NO x emissions, especially in Asia, Africa, and South America. Large increase in NO x in Asia can be eased by combination of CCS-, BLD-, and TRT-technologies. Fig. A1(a). Global distributions of emissions (g m − 2 yr − 1 ) with Ref2010 and emission changes in the 2040s with the seven scenarios relative to Ref2010 for SO 2 . Fig. A1(b). Same as Fig. A1(a), but for NO x . Fig. A1(c). Same as Fig. A1(a), but for BC. Fig.  A1(d). Same as Fig. A1(a), but for OC.