key: cord-0949581-8shuzk7e
authors: Itahashi, Syuichi; Yamamura, Yuki; Wang, Zhe; Uno, Itsushi
title: Returning long-range PM(2.5) transport into the leeward of East Asia in 2021 after Chinese economic recovery from the COVID-19 pandemic
date: 2022-04-01
journal: Sci Rep
DOI: 10.1038/s41598-022-09388-2
sha: 3533b4ba9aa0c7f981ae2f3965d37bf21436fc19
doc_id: 949581
cord_uid: 8shuzk7e

Changes in the aerosol composition of sulfate (SO(4)(2−)) and nitrate (NO(3)(−)) from 2012 to 2019 have been captured as a paradigm shift in the region downwind of China. Specifically, SO(4)(2−) dramatically decreased and NO(3)(−) dramatically increased over downwind locations such as western Japan due to the faster reduction of SO(2) emissions than NO(x) emissions and the almost constant trend of NH(3) emissions from China. Emissions from China sharply decreased during COVID-19 lockdowns in February–March 2020, after which China’s economic situation seemed to recover going into 2021. Given this substantial change in Chinese emissions, it is necessary to clarify the impact of long-range PM(2.5) transport into the leeward of East Asia. In this study, ground-based aerosol compositions observed at three sites in western Japan were analysed. The concentrations of PM(2.5), SO(4)(2−) and NO(3)(−) decreased in 2020 (during COVID-19) compared with 2018–2019 (before COVID-19). In 2021 (after COVID-19), PM(2.5) and NO(3)(−) increased and SO(4)(2−) was unchanged. This suggests the returning long-range PM(2.5) transport in 2021. From numerical simulations, the status of Chinese emissions during COVID-19 did not explain this returning impact in 2021. This study shows that the status of Chinese emissions in 2021 recovered to that before COVID-19.

We analysed aerosol compositions for February and March in each year from 2018 to 2021. For simplicity in the following discussion, 2018, for example, stands for the average in February and March in 2018. These two months (February and March) corresponded to COVID-19 lockdown periods in China in 2020 16 and so are considered as the most suitable period for evaluating the change in PM 2.5 pollution before, during and after COVID-19. Herein, we refer to the two-year average for 2018 and 2019 as before COVID-19, 2020 as during COVID-19 and 2021 as after COVID-19. Figure 1 shows the NO 2 column over East Asia as measured by the Ozone Monitoring Instrument (OMI) (see "Methods") before, during, and after COVID-19. A high NO 2 column was found over north-eastern China before COVID-19 (Fig. 1a) . Compared with its before-COVID-19 status, the NO 2 column decreased dramatically during COVID-19 over the whole of eastern China (Fig. 1b) . Then, although some parts showed decreases, almost all parts of China showed an increased NO 2 column after COVID-19 compared to during COVID-19 (Fig. 1c) . These satellite results suggest substantial deceases in emissions over China in 2020 and then increased emissions along with the economic recovery in 2021. Figure 2 shows the observational results of mean, 75 percentile (3Q), median, and 25 percentile (1Q) for changes in PM 2.5 and its composition at three remote sites in Goto, Tsushima, and Oki located in western Japan (see the lower-right corner of Fig. 1c ) from 2018 to 2021. At each of these three sites, an aerosol chemical speciation analyszer (ACSA) measured the PM 2.5 , SO 4 2− and NO 3 − concentrations (see "Methods"). Because the NO 3 − concentration had a high peak in a short time and was otherwise generally close to zero, the mean value is sometimes higher than the 3Q value. In addition to the concentration analyses in Fig. 2 , changes in the concentration of PM 2.5 and its composition (denoted as C ) during and after COVID-19 are further illustrated in Fig. 3 based on the following equations:

The observational results (Figs. 2 and 3) show that the PM 2.5 , SO 4 2− and NO 3 − concentrations decreased substantially during compared with before COVID-19. The SO 4 2− concentration decreased gradually from 2018 to 2020, whereas the NO 3 − concentration increased from 2018 to 2019 and then decreased from 2019 to 2020. The increase in NO 3 − concentration in 2019 found at Goto and Oki was consistent with our previous study that reported the change in aerosol composition as a paradigm shift 4 . In 2021, increases in PM 2.5 and NO 3 − and almost the same level of SO 4 2− have been measured. Considering the expected economic recovery suggested from the NO 2 column (Fig. 1c) , the returning impact of long-range transport can be assumed. However, because PM 2.5 concentration and composition can also be affected by the meteorological field and emission sources other than www.nature.com/scientificreports/ anthropogenic sources, further discussion is required, and we conducted modeling simulations to clarify this point.

To identify a potential reason for the variation in air pollutants found in 2021, we conducted numerical simulations for 2018-2021. The experiments are described in Table 1 (see "Methods" for details about the modeling configurations). In the names of these experiments, 'M' and 'E' denote meteorology and Chinese anthropogenic emissions, respectively (e.g., M18E18 indicates the meteorological field for 2018 and the emission inventory for 2018). From 2018 to 2020, the meteorological field and available emissions corresponding to the simulation periods were used. Through comparison with ACSA observations for 2018-2020 (Figs. 2 and 3), we confirmed that the model generally captured the observed variations of PM 2.5 , SO 4 2− and NO 3 − concentrations. Statistical analyses were conducted to evaluate the modeling performance based on the correlation coefficient (R), the normalised mean bias (NMB) and the normalised mean error (NME), and the modeling performance was judged based on metrics from a review of modeling applications (see "Methods" for definitions). The results of the statistical analyses are given in Table 2 . For PM 2.5 , the model showed moderate correlation with observations but tended to underestimate the observed concentration; in some cases, the criteria for acceptable model performance were judged as not having been met. This underestimation stemmed mainly from the underestimation of organic aerosol in comparison with a model intercomparison study 23 . For SO 4 2− and NO 3 − , the model also tended to underestimate the observed concentrations, but in most cases it was judged to have performed well within the criteria. The model's underestimation for SO 4 2− was due to the insufficient oxidation process in the current modeling 24 . Considering our previous findings from a model intercomparison study in Japan 25 , one of the possible reasons for the model's underestimation for NO 3 − was the uncertainty of NH 3 emissions. In addition to emissions, partitioning between fine-and coarse-mode NO 3 − , and also the thermodynamic equilibrium between fine-mode NO 3 − and gas-phase nitric acid (HNO 3 ), would be an important aspect to refine the modeling performance on fine-mode NO 3 −26 . For 2021, two experiments were conducted because of the lack of emission information for that year (Table 1) . In these experiments, anthropogenic emissions from China have been perturbed. One was M21E20, using the meteorological field for 2021 but the Chinese emissions for 2020 (as during with the natural emissions for 2021; the other was M21E19, also using the meteorological field for 2021 but now the Chinese emissions for 2019 (as before COVID-19) with the natural emissions for 2021. Based on these two experiments, the impact www.nature.com/scientificreports/ of supposing emission changes in 2021 can be estimated. From comparing the simulation results and the ACSA observations for Tsushima and Oki in 2021 ( Fig. 2) , the observations showed increases in PM 2.5 and NO 3 − but almost none in SO 4 2− , whereas M21E20 showed decreases except for a slight increase in NO 3 − at Oki. The results from M21E20 did not explain the variation observed in 2021, whereas M21E19 gave almost flat variations in PM 2.5 and SO 4 2− and an increase in NO 3 − . The results from M21E19 were much closer to the observational results. The relative changes in the PM 2.5 , SO 4 2− and NO 3 − concentrations after COVID-19 ( Fig. 3 ) also clarified the differences in modeling performance for 2021. Experiment M21E20 showed negative changes after COVID-19 except for NO 3 − at Oki, results that are not consistent with the observations, whereas M21E19 showed better agreement with the observations. These analyses show that M21E19 outperformed M21E20 in capturing the changes observed after COVID-19. In terms of statistical scores (Table 2) , M21E19 also outperformed M21E20. These results clarify that the emission status in 2021 was not the same as during COVID-19 and was close to the status before COVID-19. Chinese economic activity and anthropogenic emissions recovered in 2021 and impacted the air quality in downwind regions through long-range PM 2.5 transport.

Note that experiment M21E20 gave negative changes after COVID-19. Based on the differences between M20E20 and M21E20 (i.e., the M21E20 results minus the M20E20 ones), the impact due to the meteorology and changes in natural emission in 2021 can be investigated. The impacts on air quality caused by the meteorological field have been reported with regard to chemical processes 27, 28 and transport process 29 , but a negative change in SO 4 2− will lead to a positive change in NO 3 − as examined previously as a paradigm shift in aerosol composition. Therefore, the negative changes seen for PM 2.5 , SO 4 2− and NO 3 − can be attributed to changes in the meteorological field. To clarify this point, several meteorological fields are shown in Fig. S1 as Table 2 . Statistical analyses of modeling performance. Notes: for PM 2.5 and SO 4 2− , the model performance goals (marked by **) are R > 0.70, NMB < ± 10% and NME < + 35%, and the model performance criteria (marked by *) are R > 0.4, NMB < ± 30% and NME < + 50%; for NO 3 − , the model performance goals (marked by **) are NMB < ± 15% and NME < + 65%, and the model performance criteria (marked by *) are NMB < ± 65% and NME < + 115% 64 . Where M21E19 is superior to M21E20 in a statistical score, the value is given in bold font. www.nature.com/scientificreports/ the difference between 2021 and 2020 (i.e., the M21 meteorological fields minus the M20 ones) overlaid with the observational results. The positive relative humidity over the East China Sea may enhance the aqueous-phase reaction of SO 4 2and did not relate to negative changes in PM 2.5 and SO 4 2− (Fig. S1a) . The change in precipitation was inhomogeneous over the East China Sea to western Japan (Fig. S1b) . Lower planetary boundary layer (PBL) height found over oceans will lead to accumulated pollutants and did not explain the negative concentrations (Fig. S1c) . The change in the 2-m temperature partly explains the NO 3 − change as understood by thermodynamic NO 3 − production. The positive and negative temperatures over western and eastern Japan (Fig. S1d ) are generally consistent with the negative and positive changes in NO 3 − concentration over these areas (Figs. 4c and S1d). The positive change in the 10-m wind speed could contribute to the negative changes in concentrations because of the faster transport of polluted air mass (Fig. S1e) . The most plausible explanation comes from the change in the 10-m wind direction (Fig. S1f.) , which was positive over eastern China to the East China Sea and western Japan. Further detailed illustrations of wind field are shown in Fig. 4d and e. During winter, the seasonal wind pattern is generally in the northwest direction 5 , and a northwest wind direction was simulated over the East China Sea to western Japan in 2020 (Fig. 4d) ; however, a north wind direction dominated in 2021 over western Japan (Fig. 4e) . Therefore, a positive difference in wind direction with positive wind speed was calculated over western Japan (Figs. 4f and S1e and f). The results observed at Tsushima and Oki are also displayed (Fig. S2) , and the observations at Tsushima show clearly the change in wind direction from NNW in 2020 to N in 2021. The results of HYSPLIT backward trajectory 32 analyses from Goto are also presented (Fig. S3) . These show the role of long-range transport from eastern China into western Japan during wintertime, and some trajectories were originated from northeastern China or eastern Japan in 2021; indicating the change of transport pattern in 2021 as also shown by simulated meteorological field (Figs. 4d and e and f) . These findings suggest that the meteorological conditions in 2021-especially the wind direction-prevented long-range PM 2.5 transport. Nevertheless, the observations showed increases in PM 2.5 and NO 3 − and almost no change in SO 4 2− in 2021. This point supports the assertion that the changes in concentrations in 2021 cannot be explained without the returning impact of long-range PM 2.5 transport associated with the growth in Chinese anthropogenic emissions after COVID-19.

Finally, as an extension of the analysis in our previous study 4 , the SO 4 2− and NO 3 − variations before, during and after COVID-19 are summarised in Fig. 5 . The scatterplot of SO 4 2− /NO 3 − at Goto from before to during COVID-19 has an inverted V shape. In this figure, the variations of SO 4 2− and NO 3 − reported as a paradigm shift in our previous study 4 www.nature.com/scientificreports/ were measured using a quadruple-type aerosol chemical speciation monitor (Q-ACSM) for fine aerosols of PM 1 , whereas the present study analysed PM 2.5 components using an ACSA. The variations of SO 4 2− decrease and NO 3 − increase were also found based on ACSA observations up to 2019, and the dramatic decreases in both SO 4 2− and NO 3 − during COVID-19 were furthermore clarified in the present study. The modeling simulation also reproduced well this inverted V-shaped variation. At Tsushima, the variations from during COVID-19 to after COVID-19 were relatively small. From the observational results, SO 4 2− did not change but NO 3 − increased slightly, and the model results also show small variations. Note that experiment M21E20 (grey closed circles) gave deviations from these small variations in both SO 4 2− and NO 3 − . These results also show that assuming the emission status from during COVID-19 for 2021 was not appropriate. For further improvement of the modeling for Tsushima, the updated emission status over the Republic of Korea should be used. At Oki, the scatterplot www.nature.com/scientificreports/ of SO 4 2− /NO 3 − from before to during COVID-19 also has an inverted V shape as found at Goto. The variation from during COVID-19 to after COVID-19 shows no change in SO 4 2− but an increase in NO 3 − . Experiment M21E19 reproduced this variation, but M21E20 did not. These analyses of the scatterplot of SO 4 2− /NO 3 − show substantial decreases in both SO 4 2− and NO 3 − during COVID-19 but suggest the return to the situation of aerosol composition change regarding the importance of NO 3 − after COVID-19 over the leeward of East Asia. The status of aerosol composition change should be monitored continuously because the nitrogen burden over East Asia is an issue of concern both currently and for the future 33, 34 . The atmospheric input of nitrogen via dry and wet deposition processes over East Asia has been estimated in previous studies, and its significance has been clarified [35] [36] [37] [38] . Targeting the proper management of nitrogen, the nitrogen burden during and after COVID-19 should be clarified in future work.

Observations. To obtain the spatial distribution of the NO 2 column, which can be regarded as a proxy for NO x emissions 39 , space-borne OMI measurements using the level-3 daily global nitrogen dioxide product (OMNO2d version 3.0) gridded into 0.25° × 0.25°4 0 was used. We analysed the tropospheric column with clouds screened under the condition of a cloud fraction of less than 30%.

For ground-based observation of PM 2.5 and its components, ACSA observation datasets for 10 sites in Japan were available 41 . The Japanese Ministry of the Environment (MOE) created this observation network in April 2017, and two remote observation sites at Goto (128.65°E, 32.60°N) and Oki (133.20°E, 36.28°N), located in western Japan, were analysed in the present study. The ACSA measured PM 2.5 and the secondary inorganic components of SO 4 2− and NO 3 − with an hourly temporal resolution; this high temporal resolution minimises the effect of the volatilisation of NO 3 −42 . Observations at the Goto site were unavailable from the autumn of 2020 onwards because of a typhoon; to compensate for these missing observations at Goto, observations at the Tsushima (129.29°E, 34.24°N) site were used for 2020 and 2021. The locations of the three sites at Goto, Tsushima and Oki are shown in Fig. 1c .

The lockdown in Wuhan began on 23 January 2020 and lasted until 8 April 2020, and those in other Chinese cities occurred within this period 43 . In Japan, a state of emergency was proclaimed for seven prefectures including Tokyo on 7 April 2020, and this was extended nationwide on 16 April 2020; it ended on 14 May 2020 for 39 prefectures except for mega-cities such Osaka and Tokyo, and on 25 May 2020 for all prefectures in Japan 44 . Emission changes over Japan can be supposed to have been low from February to March. As also seen from the Stringency Index released from One World Data 45 , the countermeasures against COVID-19 were the strictest in China from late January to the end of March, whereas other countries in East Asia had relatively loose measures in place during this period. Given these differences between China and Japan, the period analysed in the present study was chosen as February to March to focus on the variation of long-range PM 2.5 transport according to Chinese anthropogenic emissions.

Numerical modeling simulations. The numerical modeling simulations were performed using the regional chemical transport model of the Comprehensive Air quality Model with eXtensions (CAMx) version 6.4 46 . The simulation domain covering East Asia was configured by 215 × 120 grid points with 36-km horizontal resolution, centred at 116.2°E and 37.8°N on a Lambert conformal projection. The vertical resolution was 44 non-uniformly spaced layers from the surface to 50 hPa 37 . The meteorological fields to drive CAMx were calculated by the Weather Research and Forecasting (WRF) model version 4.1.1 47 . Gas-phase and aerosol chemistry were respectively adopted for SAPRC07TC 48, 49 and the CF scheme in CAMx. Based on previous findings 50, 51 , the NO 2 aqueous-phase reaction in the SO 4 2− oxidation process was added. In addition, to enhance the production of SO 4 2− during winter haze events, the heterogeneous reaction of SO 2 on aerosol surfaces was introduced following the work on Chinese haze events 52 . The initial and lateral boundary conditions were prepared from winter-averaged values of the extended CMAQ model over the northern hemisphere 53 .

In the present study, five modeling experiments were conducted, and the simulation setups are described in Table 1 . In the names of these experiments, 'M' and 'E' stand for meteorology and anthropogenic emissions, respectively. For each year, the corresponding meteorological field generated by WRF was used. The emissions inventories were prepared as follows. Anthropogenic emissions except those in Japan and China were based on Hemispheric Transport of Air Pollution (HTAP) version 2.2 54 . For Japan, the latest available domestic emission inventory for 2015 developed by the MOE 55 was used. For China, the emission projections were based on the Multi-resolution Emission Inventory in China (MEIC). Trends in Chinese emissions up to 2017 have been reported 56 , and their variation during the COVID-19 pandemic in 2020 have also been described 12 . Neither of those two references covered 2018, so instead the average data for 2017 and 2019 were calculated. From the MEIC estimations, the reduction in NO x emissions in 2020 compared to those in 2019 during February and March was 36% and 14%, respectively, and the reduction in SO 2 emissions was 27% and 11%, respectively. NH 3 emissions were not provided for COVID-19 12 , but because their main source is agriculture, we assumed that NH 3 emissions were unchanged during COVID-19 57 and we fixed them to the levels for 2017. Biogenic emissions were prepared from the Model of Emissions of Gases and Aerosols from Nature (MEGAN) 58 using the WRF-simulated meteorological field for the corresponding period. Biomass-burning emissions were taken from the beta version of the Global Fire Emissions Database (GFED) version 4.1 59 . SO 2 emissions from 15 volcanoes in Japan were prepared from a measurement report by the Japan Meteorological Agency (JMA) 60 . Sulphur emissions from ships could influence air quality, especially in areas close to the sea 61 . Because of the limited ship emission inventory, HTAP version 2.2 was used for 2018 and 2019, and a reduction in SO 2 emissions of 77% for 2020 and 2021 was applied based on a report by the International Maritime Organization (IMO) 62 . This reduction rate was consistent with a high-resolution emission inventory developed for Japan 63 www.nature.com/scientificreports/ in 2021-regarded as after the COVID-19 pandemic-is not yet available, so to consider the emission status in the region downwind of China from the observational facts, two experiments were conducted: (i) M21E20, which assumed the emission levels during COVID-19 for those in 2021, and (ii) M21E19, which assumed the emission levels before COVID-19 for those in 2021. By comparing these two experiments, the emission status in 2021 was investigated. The model performance was evaluated statistically using the metrics of R, NMB and NME defined as

where N is the total number of paired observations (O) and models (M), and these averages are denoted as O and M , respectively. The recommended metrics based on a literature review for the United States 64 are given in Table 2 .

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 20 January 2022; Accepted: 22 March 2022 

Estimated contributions of emissions controls, meteorological factors, population growth, and changes in baseline mortality to reductions in ambient PM 2.5 and PM 2.5 -related mortality in China

Separating contributions from anthropogenic emissions and meteorology

Drivers of improved PM 2.5 air quality in China from 2013 to 2017

Paradigm shift in aerosol chemical composition over regions downwind of China

Chinese province-scale source apportionments for sulfate aerosol in 2005 evaluated by the tagged tracer method

Toward synchronous evaluation of source apportionments for atmospheric concentration and deposition of sulfate aerosol over East Asia

Modeling investigation of controlling factors in the increasing ratio of nitrate to nonseasalt sulfate in precipitation over Japan

A 15-year record (2001-2015) of the ratio of nitrate to non-sea-salt sulfate in precipitation over East Asia

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak

An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China

Changes in China's anthropogenic emissions and air quality during the COVID-19 pandemic in 2020

Abrupt decline in tropospheric nitrogen dioxide over China after the outbreak of COVID-19

Non-negligible impacts of clean air regulations on the reduction of tropospheric NO 2 over East China during the COVID-19 pandemic observed by OMI and TROPOMI

Long-range air pollution transport in East Asia during the first week of the COVID-19 lockdown in China

Impacts of COVID-19 lockdown, Spring Festival and meteorology on the NO 2 variations in early 2020 over China based on in-situ observations, satellite retrievals and model simulations

Spatiotemporal mapping and assessment of daily ground NO 2 concentrations in China using high-resolution TRO-POMI retrievals

Changes in air quality related to the control of coronavirus in China: Implications for traffic and industrial emissions

Does lockdown reduce air pollution? Evidence from 44 cities in northern China

Comprehensive analyses of source sensitivities and apportionments of PM 2.5 and ozone over Japan via multiple numerical techniques

Source-receptor relationship of transboundary particulate matter pollution between China, South Korea and Japan: Approaches, current understanding and limitations

World economic outlook: Recovery during a pandemic-Health concerns, supply disruptions, price pressures

Advantages of continuous monitoring of hourly PM 2.5 component concentrations in Japan for model validation and source sensitivity analyses

Multiphase oxidation of sulfur dioxide in aerosol particles: Implications for sulfate formation in polluted environments

Model performance differences in fine-mode nitrate aerosol during wintertime over Japan in the J-STREAM model inter-comparison study

Seasonal variation of fine-and coarse-mode nitrates and related aerosols over East Asia: Synergetic observations and chemical transport model analysis

Unexpected air pollution with marked emission reductions during the COVID-19 outbreak in China

Chemistry of atmospheric fine particles during the COVID-19 pandemic in a megacity of eastern China

Trans-regional transport of haze particles from the North China Plain to Yangtze River Delta during winter

Model estimation of sulfate aerosol sources collected at Cape Hedo during an intensive campaign in

Contribution of volcanic SO 2 emission to high concentration PM 2.5 under summertime Pacific High condition

N deposition as a threat to the World's protected area under the Convention of Biological Diversity

The shape of nitrogen to come

Importance of coarse-mode nitrate produced via sea salt as atmospheric input to East Asian oceans

MICS-Asia III: overview of model intercomparison and evaluation of acid deposition over Asia

Nitrogen burden from atmospheric deposition in East Asian oceans in 2010 based on high-resolution regional numerical modeling

Model Inter-Comparison Study for Asia (MICS-Asia) phase III: Multimodel comparison of reactive nitrogen deposition over China

Inverse estimation of NO x emissions over China and India 2005-2016: contrasting recent trends and future perspectives

OMI/Aura NO 2 Cloud-Screened Total and Tropospheric Column L3 Global Gridded 0.25 degree × 0.25 degree V3

Nitrate transboundary heavy pollution over East Asia in winter

Effect of non-pharmaceutical interventions to contain COVID-19 in China

COVID-19 information and resources

COVID-19: Stringency index

Ramboll Environment and Health. User's guide, comprehensive air quality model with extensions version 6

A description of the advanced research WRF version 4. NCAR Tech. Note, NCAR/TN

Development of the SAPRC-07 chemical mechanism

Interpreting predictions from the SAPRC07 mechanism based on regional and continental simulations

Model performance differences in sulfate aerosol in winter over Japan based on regional chemical transport models of CMAQ and CAMx

Year-round modeling of sulfate aerosol over Asia through updates of aqueousphase oxidation and gas-phase reactions with stabilized Criegee intermediates

Heterogeneous chemistry: a mechanism missing in current models to explain secondary inorganic aerosol formation during the January 2013 haze episode in North China

Hemispheric CMAQ model version 5.3beta output data-2016 seasonally averaged 108km for N. hemisphere, UNC Dataverse, V1

2: A mosaic of regional and global emission grid maps for 2008 and 2010 to study hemispheric transport of air pollution

Trends in China's anthropogenic emissions since 2010 as the consequence of clean air actions

Ministory of Environment. PM 2.5 emission inveotory

Current and future global climate impacts resulting from COVID-19

The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1). An extended and updated framework for modeling biogenic emissions

Global fire emissions estimates during 1997-2016

Long-term trends of satellite-based fine-mode aerosol optical depth over the Seto Inland Sea, Japan, over two decades

Air pollution and energy efficiency: Study on the effects of the entry into force of the global 0.5% fuel oil sulphur content limit on human health

Development of air pollutants emission inventories for ships around Japan on a high geographical resolution

Recommendations on statistics and benchmarks to assess photochemical model performance

This work was supported by MEXT/JSPS KAKENHI grant number JP18H03359 and the Environment Research and Technology Development Fund (JPMEERF20195003) of the Environmental Restoration and Conservation Agency of Japan.

S.I. conducted the model simulations and analysed the model results, surface observations and satellite measurements. Y.Y. conducted the surface observations at the Tsushima site and analysed the observation dataset. Z.W. and I.U. conducted the analyses of the surface observations. S.I. and I.U. designed the study and also wrote the manuscript with valuable input from all co-authors. All authors reviewed the manuscript.

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-09388-2.Correspondence and requests for materials should be addressed to S.I.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.