<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article"><?xmltex \bartext{Research article}?>
  <front>
    <journal-meta><journal-id journal-id-type="publisher">ACP</journal-id><journal-title-group>
    <journal-title>Atmospheric Chemistry and Physics</journal-title>
    <abbrev-journal-title abbrev-type="publisher">ACP</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Atmos. Chem. Phys.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">1680-7324</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/acp-23-3311-2023</article-id><title-group><article-title>Global impact of the COVID-19 lockdown on surface concentration and health risk
of atmospheric benzene</article-title><alt-title>Ambient benzene variation during the COVID-19 period</alt-title>
      </title-group><?xmltex \runningtitle{Ambient benzene variation during the COVID-19 period}?><?xmltex \runningauthor{C.~Ling et al.}?>
      <contrib-group>
        <contrib contrib-type="author" equal-contrib="yes" corresp="no" rid="aff1">
          <name><surname>Ling</surname><given-names>Chaohao</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" equal-contrib="yes" corresp="no" rid="aff2">
          <name><surname>Cui</surname><given-names>Lulu</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="yes" rid="aff3 aff4">
          <name><surname>Li</surname><given-names>Rui</given-names></name>
          <email>rli@geo.ecnu.edu.cn</email>
        </contrib>
        <aff id="aff1"><label>1</label><institution>School of History and Geography, Minnan Normal University, Zhangzhou, 363000, China</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>State Key Joint Laboratory of Environment Simulation and Pollution Control, School of Environment, Tsinghua University, Beijing, 100084, China</institution>
        </aff>
        <aff id="aff3"><label>3</label><institution>Key Laboratory of Geographic Information Science of the Ministry of Education, <?xmltex \hack{\break}?>School of Geographic Sciences, East China Normal University, Shanghai, 200241, China</institution>
        </aff>
        <aff id="aff4"><label>4</label><institution>Institute of Eco-Chongming (IEC), 20 Cuiniao Road, Chenjia Town, <?xmltex \hack{\break}?>Chongming District, Shanghai, 202162, China</institution>
        </aff><author-comment content-type="econtrib"><p>These authors contributed equally to this work.</p></author-comment>
      </contrib-group>
      <author-notes><corresp id="corr1">Rui Li (rli@geo.ecnu.edu.cn)</corresp></author-notes><pub-date><day>16</day><month>March</month><year>2023</year></pub-date>
      
      <volume>23</volume>
      <issue>5</issue>
      <fpage>3311</fpage><lpage>3324</lpage>
      <history>
        <date date-type="received"><day>6</day><month>December</month><year>2022</year></date>
           <date date-type="rev-request"><day>2</day><month>January</month><year>2023</year></date>
           <date date-type="rev-recd"><day>24</day><month>February</month><year>2023</year></date>
           <date date-type="accepted"><day>28</day><month>February</month><year>2023</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2023 </copyright-statement>
        <copyright-year>2023</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://acp.copernicus.org/articles/.html">This article is available from https://acp.copernicus.org/articles/.html</self-uri><self-uri xlink:href="https://acp.copernicus.org/articles/.pdf">The full text article is available as a PDF file from https://acp.copernicus.org/articles/.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d1e127">To curb the spread of the COVID-19 pandemic, many countries around the world
imposed an unprecedented lockdown, producing reductions in pollutant
emissions. Unfortunately, the lockdown-driven global ambient benzene changes
still remain unknown. An ensemble machine-learning model coupled with chemical transport models (CTMs) was applied to estimate global
high-resolution ambient benzene levels. Afterwards, the extreme gradient boosting (XGBoost) algorithm
was employed to decouple the contributions of meteorology and emission
reduction to ambient benzene. The change ratio (<inline-formula><mml:math id="M1" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">dew</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) of the deweathered
benzene concentration from the pre-lockdown to lockdown period was in the order
of India (<inline-formula><mml:math id="M2" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.6</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M3" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> Europe (<inline-formula><mml:math id="M4" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">21.9</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M5" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United
States (<inline-formula><mml:math id="M6" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M7" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China (<inline-formula><mml:math id="M8" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">15.6</mml:mn></mml:mrow></mml:math></inline-formula> %). The detrended change (<inline-formula><mml:math id="M9" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>)
of the deweathered benzene level (change ratio in 2020 <inline-formula><mml:math id="M10" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula> change ratio in 2019)
followed the order of India (<inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>*</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M12" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> Europe (<inline-formula><mml:math id="M13" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">13.9</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M14" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China (<inline-formula><mml:math id="M15" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">13.3</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M16" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United States
(<inline-formula><mml:math id="M17" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6.00</mml:mn></mml:mrow></mml:math></inline-formula> %). Emission reductions derived from industrial activities
and transportation were major drivers for the benzene decrease during the
lockdown period. The highest decreasing ratio of ambient benzene in India
might be associated with local serious benzene pollution during the
business-as-usual period and restricted transportation after lockdown.
Substantial decreases in atmospheric benzene levels had significant health
benefits. The global average lifetime carcinogenic risk (LCR) and hazard
index (HI) decreased from <inline-formula><mml:math id="M18" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.89</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M19" display="inline"><mml:mrow><mml:mn mathvariant="normal">5.90</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.51</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M21" display="inline"><mml:mrow><mml:mn mathvariant="normal">5.40</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>,
respectively. China and India showed higher health benefits due to
benzene pollution mitigation compared with other countries, highlighting the
importance of benzene emission reduction.</p>
  </abstract>
    
<funding-group>
<award-group id="gs1">
<funding-source>National Natural Science Foundation of China</funding-source>
<award-id>42107113</award-id>
</award-group>
<award-group id="gs2">
<funding-source>Natural Science Foundation of Fujian Province</funding-source>
<award-id>2022J05178</award-id>
</award-group>
<award-group id="gs3">
<funding-source>Natural Science Foundation of Hunan Province</funding-source>
<award-id>2021JJ30555</award-id>
</award-group>
</funding-group>
</article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e393">Volatile organic compounds (VOCs) are an important class of organic
pollutants in the urban air and have aroused great attention (Kamal et al.,
2016; Koppmann, 2008; Mozaffar and Zhang, 2020). As one of the typical toxic
VOC species, benzene has a variety of negative impacts on human health
including respiratory irritation, asthma, and allergies (Cui et al., 2019;
Kim et al., 2013; Tang et al., 2007). Moreover, benzene has high chemical
reactivity and could participate in photochemical reactions in the
atmosphere, thereby leading to the formation of secondary organic aerosols
(SOAs) and ozone (Dumanoglu et al., 2014; Hsu et al., 2018; Li et al., 2019).
Given its high toxicity to human health and tremendous harm to air quality
(Dumanoglu et al., 2014; Lu et al., 2020), it is highly imperative to
decrease ambient<?pagebreak page3312?> benzene concentration. It has been well documented that
ambient benzene mainly originates from anthropogenic emissions (Mozaffar and
Zhang, 2020; Pakkattil et al., 2021). Therefore, understanding the response
of ambient benzene to anthropogenic emissions is favourable to evaluating the
effectiveness of abatement strategies and informing policy decisions.</p>
      <p id="d1e396">Recently, the ongoing global COVID-19 outbreak has resulted in paroxysmal
public health responses including travel restrictions, lockdown, curfews,
and quarantines around the world. These drastic lockdown measures have inevitably
triggered sweeping disruptions of social and economic activities and have
further affected the emissions and concentrations of some air pollutants
(Bauwens et al., 2020; Berg et al., 2021; Doumbia et al., 2021; P. Zheng et
al., 2021). The unexpected public health emergency provided us with an
unprecedented chance to assess the response of air pollutants to emission
reduction. Bauwens et al. (2020) observed that the average NO<inline-formula><mml:math id="M22" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula>
column in China during January–April 2020 decreased by about 40 % relative
to the same period in 2019 due to dramatic decreases in NO<inline-formula><mml:math id="M23" display="inline"><mml:msub><mml:mi/><mml:mi>x</mml:mi></mml:msub></mml:math></inline-formula>
emissions. Later on, Keller et al. (2021a) analysed the impact of the
COVID-19 lockdown on global NO<inline-formula><mml:math id="M24" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> concentrations and found that the
surface NO<inline-formula><mml:math id="M25" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> concentrations were 18 % lower than the business-as-usual scenario
from February 2020 onward. In addition, Hammer et al. (2021) estimated that
population-weighted mean PM<inline-formula><mml:math id="M26" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> concentrations in China, Europe, and
North America experienced changes of <inline-formula><mml:math id="M27" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">11</mml:mn></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M28" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">15</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M29" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula> to 1, and <inline-formula><mml:math id="M30" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula> to 1 <inline-formula><mml:math id="M31" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, respectively, during the COVID-19 lockdown period. Compared with
NO<inline-formula><mml:math id="M32" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> and PM<inline-formula><mml:math id="M33" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula>, ambient SO<inline-formula><mml:math id="M34" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> levels in China (<inline-formula><mml:math id="M35" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4.6</mml:mn></mml:mrow></mml:math></inline-formula> %) and
India (<inline-formula><mml:math id="M36" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">14</mml:mn></mml:mrow></mml:math></inline-formula> %) did not experience marked variations after lockdown (Zhang
et al., 2021; Zhao et al., 2020). It should be noted that the global O<inline-formula><mml:math id="M37" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:math></inline-formula>
concentration even increased by up to 50 % during this period (Keller et al.,
2021a). To date, most of the current studies have focused on regional or global
criteria pollutant (e.g. PM<inline-formula><mml:math id="M38" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula>, NO<inline-formula><mml:math id="M39" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula>, and O<inline-formula><mml:math id="M40" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:math></inline-formula>) concentration
changes after the COVID-19 outbreak, while few studies have assessed the impact
of the COVID-19 lockdown on ambient benzene levels.</p>
      <p id="d1e589">Mor et al. (2021) observed that the atmospheric
benzene level in Chandigarh, India, decreased by 27 % during the COVID-19
period. Afterwards, Pakkattil et al. (2021) demonstrated that the ambient
benzene levels in Delhi (<inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">93</mml:mn></mml:mrow></mml:math></inline-formula> %) and Mumbai (<inline-formula><mml:math id="M42" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">72</mml:mn></mml:mrow></mml:math></inline-formula> %) have undergone
drastic decreases after the COVID-19 lockdown. In China, Pei et al. (2022)
revealed that the VOC concentration in the Pearl River Delta (PRD) decreased by
19 % and the decrease rate of ambient benzene reached <inline-formula><mml:math id="M43" display="inline"><mml:mo>∼</mml:mo></mml:math></inline-formula> 40 % after lockdown. In Europe, Cai et al. (2022) revealed that the
ambient benzene level in Orléans even slightly increased after lockdown,
which might be associated with unfavourable meteorological conditions.
Although ground-level measurement could reflect regional ambient
benzene changes during the COVID-19 lockdown period to some extent, few
regions, especially in developing countries, have collected sufficient
observations for ambient benzene exposure assessment (Geddes et al., 2016;
Van Donkelaar et al., 2015). Moreover, the limited monitoring sites around
the world cannot accurately reflect the global benzene pollution because of
large spatial gaps and restricted spatial representativeness of these
ground-based sites (Shi et al., 2018). The health effect assessment based on
these scarce sites alone inevitably increases the probability of exposure
misclassification (Ling and Li, 2021). Fortunately, chemical transport
models (CTMs) give us an unparalleled chance to capture the full-coverage
ambient benzene level at the global scale. Although CTMs generally show
various biases owing to high uncertainties in initial conditions, input
variables, and parameterizations (Ivatt and Evans, 2020), the
machine-learning bias-correction method can significantly reduce bias in
air quality models (Bocquet et al., 2015). To date, some studies have
developed multiple machine-learning models to estimate the concentrations of
PM<inline-formula><mml:math id="M44" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> (Wei et al., 2021a, b), NO<inline-formula><mml:math id="M45" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> (Wei et al.,
2023), and O<inline-formula><mml:math id="M46" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:math></inline-formula> (Wei et al., 2022) around the world. Unfortunately, no
study has employed the ensemble technique to analyse the change in global
ambient benzene after the COVID-19 outbreak. Besides, nearly all of the current
studies have only used original observation data to assess the impact of the COVID-19
lockdown on the ambient benzene level (Pakkattil et al., 2021). Actually, the
concentrations of air pollutants are not only controlled by emissions but
also modulated by complex meteorological conditions (Hammer et al., 2021).
For instance, some pioneering studies have revealed that several severe haze
episodes still occurred even with the strict restrictions put in place in
China (Chang et al., 2020; Huang et al., 2021). Hence, it is necessary to
remove the effects of meteorological parameters and then to further quantify
the isolated contribution of emission reduction to the global ambient benzene
level and health risks during the COVID-19 lockdown period.</p>
      <p id="d1e647">In our study, a machine-learning model coupled with CTMs was applied to
estimate the global ambient benzene concentrations from 23 January to 30 June in 2019 and 2020. At first, the CTMs' output, emission inventory,
meteorological parameters, and many other geographical covariates were
integrated into the ensemble decision tree model to obtain global
full-coverage benzene concentrations in the atmosphere. Then, we also
examined the synergetic impacts from the anthropogenic emissions and
meteorological factors during the pre-lockdown and lockdown periods.
Finally, we estimated the emission-induced benzene concentrations before and
after the COVID-19 lockdown and quantified the benzene-related health benefits
due to the COVID-19 lockdown in major regions around the world. This study shows
important implications for developing control strategies to alleviate global
atmospheric benzene pollution.</p>
</sec>
<?pagebreak page3313?><sec id="Ch1.S2">
  <label>2</label><title>Data and methods</title>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>Data preparation</title>
<sec id="Ch1.S2.SS1.SSS1">
  <label>2.1.1</label><title>Ground-level benzene observation</title>
      <p id="d1e672">Our analysis was performed based on the recent development of unprecedented
public access to ground-level air quality observations. In our study, we
collected an air quality dataset of hourly surface benzene observations at
669 sites at the global scale during 23 January–30 June in 2019 and 2020
(Fig. S1). The start date of the COVID-19 lockdown in China was 23 January,
and the national lockdown lasted for about 1 month. The
lockdown end date in Wuhan was  8 April. The start and end dates of lockdown
in India were 25 March and  25 April, respectively. The lockdown
in the United States occurred firstly in California from  19 March, and then
the lockdown lasted for 1–2 months. The lockdown dates in most European
countries lasted from March to May. A detailed spatial distribution of
the monitoring sites in India, Europe, and the United States is depicted in Fig. S1. The surface benzene dataset in India was downloaded from the Central
Pollution Control Board (CPCB) online database, which has been widely
utilized in previous studies (Mahato et al., 2020; Mor et al., 2021; Sharma
et al., 2020). The CPCB database provides data quality assurance (QA) or
quality control (QC) programmes by developing strict procedures for sampling,
analysis, and calibration (Gurjar et al., 2016). The ground-level benzene
observations in Europe and the United States were collected from the air quality
data portal of the European Environment Agency (EEA) and United States
Environmental Protection Agency (EPA), respectively. Only days with more
than 12 h of available data are included in the analysis. All of the hourly
data were average to the daily scale.</p>
</sec>
<sec id="Ch1.S2.SS1.SSS2">
  <label>2.1.2</label><title>Independent variables</title>
      <p id="d1e683">The daily benzene concentrations at a global scale were simulated using
the GEOS-Chem (<uri>https://geos-chem.seas.harvard.edu/</uri>, last access: 10 March 2023) model (v12-01), which included the full gaseous
HO<inline-formula><mml:math id="M47" display="inline"><mml:msub><mml:mi/><mml:mi>x</mml:mi></mml:msub></mml:math></inline-formula>–O<inline-formula><mml:math id="M48" display="inline"><mml:msub><mml:mi/><mml:mi>x</mml:mi></mml:msub></mml:math></inline-formula>–NO<inline-formula><mml:math id="M49" display="inline"><mml:msub><mml:mi/><mml:mi>x</mml:mi></mml:msub></mml:math></inline-formula>–CO–non-methane VOC (NMVOC) chemistry and online aerosol
calculations. The simulation used assimilated meteorological observations
(GEOS MERRA-2) at 2<inline-formula><mml:math id="M50" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> <inline-formula><mml:math id="M51" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> 2.5<inline-formula><mml:math id="M52" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> horizontal resolution with
72 vertical levels for the years 2019 and 2020. The anthropogenic emission
inventory in 2019 was collected from the Community Emissions Data System (CEDS).
Then, the emission inventory in 2020 was calculated based on that in 2019
and on an updated adjustment factor proposed by Doumbia et al. (2021).</p>
      <p id="d1e742">The meteorological parameters were obtained from the NASA Goddard Earth
Observing System Composition Forecast (GEOS-CF) model (Keller et al.,
2021b). GEOS-CF integrates the GEOS-Chem atmospheric chemistry model into
the GEOS Earth system model (Hu et al., 2018; Long et al., 2015) and
provides global hourly analyses of meteorological variables at
0.25<inline-formula><mml:math id="M53" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> spatial resolution (Keller et al., 2021b). Meteorological
parameters including surface pressure (<inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>), relative humidity (RH), 2 m air
temperature (<inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>), total precipitation (TPREC), the 10 m latitudinal wind
component (<inline-formula><mml:math id="M56" display="inline"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn mathvariant="normal">10</mml:mn><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>), the 10 m longitudinal wind component (<inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mn mathvariant="normal">10</mml:mn><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>), and boundary
layer height (BLH) obtained from GEOS-CF were used to develop the model
(Fig. S2). In addition, cropland, forest, grassland, shrubland, and barren
land have also been integrated into the final model (Liu et al., 2020).</p>
      <p id="d1e810">All of the independent variables collected from multiple sources were
regridded to 0.25<inline-formula><mml:math id="M58" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> grids using spatial-interpolation algorithms.
During the process of model development, the most important procedure was to
remove some redundant explanatory variables and then to determine the
optimal variable group. The basic principle of the variable selection was to
eliminate the less important predictors. These variables generally suggested
that the <inline-formula><mml:math id="M59" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value of the submodel did not experience a significant
decrease or that it even experienced a slight increase when these redundant variables
were removed from the model. Finally, a total of 64 001 samples and 7
variables were utilized to predict the ambient benzene concentrations at the
global scale.</p>
</sec>
</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>Model development</title>
<sec id="Ch1.S2.SS2.SSS1">
  <label>2.2.1</label><title>The ensemble model development for atmospheric benzene estimates</title>
      <?pagebreak page3314?><p id="d1e849">In pioneering studies, random forest (RF), extreme gradient boosting
(XGBoost), and light gradient boosting machine (LightGBM) have exhibited good estimation accuracy (Li et al., 2021; Wei et al., 2021b). The RF model contains a great number of
decision trees; each one experiences an independent sampling procedure,
and all of these trees show the same distributions (Breiman, 2001). The RF model often
displays excellent prediction performance owing to the injected randomness.
The model accuracy is strongly dependent on the number of trees, splitting
features, and variable group. The detailed procedures are summarized as
follows:

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M60" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E1"><mml:mtd><mml:mtext>1</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>Z</mml:mi></mml:munderover><mml:msub><mml:mi>c</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mi>R</mml:mi><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>∈</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E2"><mml:mtd><mml:mtext>2</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msup><mml:msub><mml:mi>c</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mi mathvariant="normal">Δ</mml:mi></mml:msup><mml:mo>=</mml:mo><mml:mtext>average</mml:mtext><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mi mathvariant="normal">|</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>∈</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E3"><mml:mtd><mml:mtext>3</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mtext>BR</mml:mtext><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>(</mml:mo><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>q</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mo mathvariant="italic">{</mml:mo><mml:mi>X</mml:mi><mml:mi mathvariant="normal">|</mml:mi><mml:msub><mml:mi>X</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>≤</mml:mo><mml:mi>q</mml:mi><mml:mo mathvariant="italic">}</mml:mo><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="italic">&amp;</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:msub><mml:mtext>BR</mml:mtext><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>(</mml:mo><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>q</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mo mathvariant="italic">{</mml:mo><mml:mi>X</mml:mi><mml:mi mathvariant="normal">|</mml:mi><mml:msub><mml:mi>X</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>&gt;</mml:mo><mml:mi>q</mml:mi><mml:mo mathvariant="italic">}</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E4"><mml:mtd><mml:mtext>4</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:munder><mml:mo movablelimits="false">min⁡</mml:mo><mml:mrow><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>q</mml:mi></mml:mrow></mml:munder><mml:mfenced open="[" close="]"><mml:mrow><mml:mo movablelimits="false">min⁡</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>(</mml:mo><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>q</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:munder><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo movablelimits="false">min⁡</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>(</mml:mo><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>q</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:munder><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E5"><mml:mtd><mml:mtext>5</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:msub><mml:mi>c</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi mathvariant="normal">Δ</mml:mi></mml:msup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mtext>average</mml:mtext><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mi mathvariant="normal">|</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>∈</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>(</mml:mo><mml:mi>p</mml:mi><mml:mi>q</mml:mi><mml:mo>)</mml:mo><mml:mo>)</mml:mo><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mtext>and</mml:mtext></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:msub><mml:mi>c</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mi mathvariant="normal">Δ</mml:mi></mml:msup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mtext>average</mml:mtext><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mi mathvariant="normal">|</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>∈</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>(</mml:mo><mml:mi>p</mml:mi><mml:mi>q</mml:mi><mml:mo>)</mml:mo><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

              where (<inline-formula><mml:math id="M61" display="inline"><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M62" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) is the sample for <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>, 2, …, <inline-formula><mml:math id="M64" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> in
<inline-formula><mml:math id="M65" display="inline"><mml:mi>Q</mml:mi></mml:math></inline-formula> regions (<inline-formula><mml:math id="M66" display="inline"><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, …, <inline-formula><mml:math id="M68" display="inline"><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>); <inline-formula><mml:math id="M69" display="inline"><mml:mi>R</mml:mi></mml:math></inline-formula> denotes the weight
of each branch; BR represents the decision tree branch; <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msup><mml:msub><mml:mi>c</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mi mathvariant="normal">Δ</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> represents the optimal value; <inline-formula><mml:math id="M71" display="inline"><mml:mi>p</mml:mi></mml:math></inline-formula>
is the feature variable; <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is the average of the left branch; <inline-formula><mml:math id="M73" display="inline"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is
the average of the right branch; and <inline-formula><mml:math id="M74" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> represents the split point.</p>
      <p id="d1e1387">The XGBoost model is an improved algorithm of gradient boosting decision tree
(GBDT) models, and loss functions have been extended to the second-order
function. The detailed XGBoost algorithm is shown as the following formula
(Zhai and Chen, 2018):
              <disp-formula id="Ch1.E6" content-type="numbered"><label>6</label><mml:math id="M75" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:mo>[</mml:mo><mml:mi>l</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:msup><mml:mi>y</mml:mi><mml:mrow><mml:mi mathvariant="normal">Λ</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msup><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mrow><mml:msup><mml:mi>y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msup></mml:msub><mml:mi>l</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:msup><mml:mi>y</mml:mi><mml:mrow><mml:mi mathvariant="normal">Λ</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msup><mml:mo>)</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:msubsup><mml:mo>∂</mml:mo><mml:mrow><mml:msup><mml:mi>y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mi>l</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:msup><mml:mi>y</mml:mi><mml:mrow><mml:mi mathvariant="normal">Λ</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msup><mml:mo>)</mml:mo><mml:msubsup><mml:mi>f</mml:mi><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>]</mml:mo><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ε</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
            where <inline-formula><mml:math id="M76" display="inline"><mml:mrow><mml:msup><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> is the cost function at the <inline-formula><mml:math id="M77" display="inline"><mml:mi>t</mml:mi></mml:math></inline-formula>th period, <inline-formula><mml:math id="M78" display="inline"><mml:mo>∂</mml:mo></mml:math></inline-formula> represents the derivative of the original function, <inline-formula><mml:math id="M79" display="inline"><mml:mrow><mml:msubsup><mml:mo>∂</mml:mo><mml:mrow><mml:msup><mml:mi>y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> is the second derivative of the original function, <inline-formula><mml:math id="M80" display="inline"><mml:mi>l</mml:mi></mml:math></inline-formula> is the
differentiable convex loss function that reflects the difference of the predicted
value (<inline-formula><mml:math id="M81" display="inline"><mml:mrow><mml:msup><mml:mi>y</mml:mi><mml:mi mathvariant="normal">Λ</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula>) of the <inline-formula><mml:math id="M82" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>th instance at the <inline-formula><mml:math id="M83" display="inline"><mml:mi>t</mml:mi></mml:math></inline-formula>th period and
the target value (<inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>), <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> represents the increment, and
<inline-formula><mml:math id="M86" display="inline"><mml:mrow><mml:mi mathvariant="italic">ε</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the regularizer.</p>
      <p id="d1e1752">The LightGBM model is an updated version of the XGBoost method and significantly
improves the running speed of the modelling process. Moreover, this method could
decrease the cache miss by a large margin and further improve the
predictive accuracy. The detailed algorithms are as follows (Sun et al.,
2020):

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M87" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E7"><mml:mtd><mml:mtext>7</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mover accent="true"><mml:mi>f</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mi>arg⁡</mml:mi><mml:munder><mml:mo movablelimits="false">min⁡</mml:mo><mml:mi>f</mml:mi></mml:munder><mml:mi>L</mml:mi><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E8"><mml:mtd><mml:mtext>8</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>f</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>T</mml:mi></mml:munderover><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E9"><mml:mtd><mml:mtext>9</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi mathvariant="normal">Γ</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:mo>(</mml:mo><mml:msub><mml:mi>g</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msubsup><mml:mi>f</mml:mi><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

              where <inline-formula><mml:math id="M88" display="inline"><mml:mover accent="true"><mml:mi>f</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover></mml:math></inline-formula> is the lowest value of the cost function; <inline-formula><mml:math id="M89" display="inline"><mml:mrow><mml:mi>L</mml:mi><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>
is the cost function; <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> denotes the total regression trees;
<inline-formula><mml:math id="M91" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> represents each regression tree; and <inline-formula><mml:math id="M92" display="inline"><mml:mrow><mml:msub><mml:mi>g</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> represent
the first- and second-order gradient statistics of the cost function,
respectively.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1" specific-use="star"><?xmltex \currentcnt{1}?><?xmltex \def\figurename{Figure}?><label>Figure 1</label><caption><p id="d1e2024">The workflow of global atmospheric benzene modelling. CTM output
represents the simulated benzene concentration based on the GEOS-Chem model.
Meteo denotes the meteorological parameters derived from GEOS-CF reanalysis.
Emission represents the daily emissions of benzene. MOY and DOY are the month
of year and day of year, respectively. Simulated benzene represents the
predicted benzene concentrations based on the ensemble model. Deweathered
benzene denotes the benzene concentration after removing meteorological
effects.</p></caption>
            <?xmltex \igopts{width=398.338583pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/3311/2023/acp-23-3311-2023-f01.png"/>

          </fig>

      <p id="d1e2033">Although all of these models showed good performance in predicting air
pollutants, nearly all of the submodels still suffered from some
weaknesses in their prediction accuracy. Hence, it was necessary to collocate
these models using a back-propagation neutral network (BPNN) to further
simulate daily ambient benzene concentrations at the global scale. As
depicted in Fig. 1, three submodels including RF, XGBoost, and LightGBM
were stacked through the BPNN model to simulate the daily atmospheric benzene
levels at the global scale. Firstly, a 5-fold cross-validation method was
utilized to train each submodel to determine the optimal hyperparameter.
Then, the BPNN method was employed to further train the estimated
concentrations of three submodels against the observations (Fig. 1).
Lastly, the global ambient benzene concentrations were predicted on the
basis of the ensemble model.</p>
</sec>
<sec id="Ch1.S2.SS2.SSS2">
  <label>2.2.2</label><title>The meteorology-normalized benzene estimates</title>
      <p id="d1e2044">The ambient benzene concentration was influenced by both meteorological
parameters and emissions. To isolate the contribution of emissions, the
impacts of meteorological conditions should be removed. In our study, the
XGBoost approach was utilized to eliminate the impacts of meteorological
conditions. The simulated benzene concentration in each grid
(0.25<inline-formula><mml:math id="M94" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>) based on the method in Sect. 2.2.1 was treated as the
dependent variable. The daily benzene emissions, meteorological factors,
month of year (MOY), and day of year (DOY) in each grid were regarded as the
explanatory variables. The raw dataset was randomly classified into a
training dataset (90 % of input dataset) for developing the XGBoost model,
and the remaining samples were regarded as the test dataset. After the
development of the XGBoost model, the weather-normalized technique was
employed to predict the ambient benzene concentration at a specific time
point. The detailed deweathered algorithms were first introduced by Grange and
Carslaw (2019). The meteorology-normalized benzene level served as
the concentrations contributed by emissions alone. The differences in total
and deweathered benzene concentrations were regarded as the concentration
contributed by meteorology. In addition, the cross-validation (CV) <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value of the model used
for the separation of meteorology and emissions should also be higher than
0.50; otherwise the model could be considered unreliable.</p>
</sec>
</sec>
<sec id="Ch1.S2.SS3">
  <label>2.3</label><title>Health effect assessment</title>
      <?pagebreak page3315?><p id="d1e2076">In our study, the carcinogenic and non-carcinogenic risks of ambient benzene
were assessed based on the standard methodology of the United States
Environmental Protection Agency (EPA). The carcinogenic and
non-carcinogenic risks induced by benzene exposure were evaluated based
on the lifetime carcinogenic risk (LCR) and hazard index (HI). The formulas
for calculating benzene intake (BI), LCR, and HI are as follows (Table S1):

                <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M96" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E10"><mml:mtd><mml:mtext>10</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtext>BI</mml:mtext><mml:mo>=</mml:mo><mml:mo>(</mml:mo><mml:mi>C</mml:mi><mml:mo>×</mml:mo><mml:mtext>ET</mml:mtext><mml:mo>×</mml:mo><mml:mtext>EF</mml:mtext><mml:mo>×</mml:mo><mml:mtext>ED</mml:mtext><mml:mo>)</mml:mo><mml:mo>/</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">365</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">24</mml:mn><mml:mo>×</mml:mo><mml:mtext>AT</mml:mtext><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E11"><mml:mtd><mml:mtext>11</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtext>HI</mml:mtext><mml:mo>=</mml:mo><mml:mtext>BI</mml:mtext><mml:mo>/</mml:mo><mml:mtext>RfC</mml:mtext><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E12"><mml:mtd><mml:mtext>12</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtext>LCR</mml:mtext><mml:mo>=</mml:mo><mml:mtext>BI</mml:mtext><mml:mo>×</mml:mo><mml:mtext>IUR</mml:mtext><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            where <inline-formula><mml:math id="M97" display="inline"><mml:mi>C</mml:mi></mml:math></inline-formula> (<inline-formula><mml:math id="M98" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) denotes the concentration of the corresponding
ambient benzene; ET is the exposure time; EF represents the annual exposure
frequency (d a<inline-formula><mml:math id="M99" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>); ED is the exposure duration (a); AT<inline-formula><mml:math id="M100" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">nca</mml:mi></mml:msub></mml:math></inline-formula> and
AT<inline-formula><mml:math id="M101" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">ca</mml:mi></mml:msub></mml:math></inline-formula> denote the average exposure time for carcinogenic and
non-carcinogenic risks (a), respectively; BI means the benzene intake; RfC
represents the reference dose (<inline-formula><mml:math id="M102" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>); and IUR is the inhalation risk
(1/<inline-formula><mml:math id="M103" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>). The non-carcinogenic risk of the ambient benzene is
considered high when HI is above 1.0, whereas the health risk is not
obvious when HI is below 1.0. The carcinogenic risk was regarded as definite
risk when LCR was higher than <inline-formula><mml:math id="M104" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, while it was treated
as possible risk when this indicator was located between <inline-formula><mml:math id="M105" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>. The risk was treated as negligible
when the indicator was lower than <inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> (Dumanoglu et
al., 2014; Li et al., 2017).</p>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Results and discussion</title>
<sec id="Ch1.S3.SS1">
  <label>3.1</label><title>The model fitting and validation</title>
      <p id="d1e2355">The ensemble model was utilized to estimate the ambient benzene
concentrations at the global scale during 23 January–30 June in 2019 and
2020. The CV <inline-formula><mml:math id="M108" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value of the ensemble model (<inline-formula><mml:math id="M109" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.60</mml:mn></mml:mrow></mml:math></inline-formula>) was significantly higher than that of RF (0.52), XGBoost (0.53),
and LightGBM (0.55) (Fig. S3). Nevertheless, both the root-mean-square
error (RMSE) (1.18 <inline-formula><mml:math id="M110" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) and the mean absolute error (MAE) (0.59 <inline-formula><mml:math id="M111" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) of the ensemble model were significantly lower than those
of RF (RMSE and MAE: 1.41 and 0.72 <inline-formula><mml:math id="M112" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>), XGBoost (RMSE and MAE:
1.37 and 0.70 <inline-formula><mml:math id="M113" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>), and LightGBM (RMSE and MAE: 1.34 and 0.69 <inline-formula><mml:math id="M114" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>). The higher <inline-formula><mml:math id="M115" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value and the lower RMSE and MAE
suggested the higher accuracy of the ensemble model in air quality
simulation. In pioneering studies, Wolpert (1992) confirmed that the
joint use of multiple statistical models could decrease the probability of
overfitting and strengthen the predictive accuracy and transferability of
final models. Besides, our previous studies also demonstrated that the
stacking of various decision tree models could significantly outperform
an individual model because each decision tree model could suffer from some
weaknesses (Li et al., 2021). For instance, the dataset in the RF model
appeared to be overfitted when much noise existed in the training data of
regression problems (Breiman, 2001). Besides, the RF model might
underestimate/overestimate the extreme values of ambient benzene (Xue et
al., 2019), which could be neutralized by the XGBoost algorithm through the
boosting method (Li et al., 2020). For XGBoost algorithm, excessive leaf
nodes often showed low splitting gain, while the LightGBM model could make
up for this defect (Nemeth et al., 2019). Overall, the combination of these
decision tree models could overcome the weaknesses of the individual
models and enhance the robustness of the final model.</p>
      <p id="d1e2491">Although 10-fold CV has verified that the modelling performance of the ensemble
model was superior to that of the individual models, this method cannot examine the
spatial transferability of this model. In our study, many regions (except
India, Europe, and the United States) lacked monitoring sites for
ambient benzene. Fortunately, the CTMs' output provided a strong proxy to
predict the daily ambient benzene concentrations before and after the COVID-19
outbreak. In order to examine the spatial transferability of the
ensemble model, cross-validation was performed. In each round,
two-thirds of the benzene dataset in India, Europe, and the United States
was applied to train the model, and the remaining third was utilized to examine
the model (e.g. India <inline-formula><mml:math id="M116" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula> Europe for training and the United States for
testing). After three rounds, all of the simulated benzene concentrations
were compared with the corresponding observed values. As shown in Fig. S4,
the out-of-bag <inline-formula><mml:math id="M117" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value reached 0.58, which was slightly lower than the
<inline-formula><mml:math id="M118" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value (0.60) of the training model. In addition, the RMSE and MAE of the
fitting equation for the out-of-bag data were 1.18 and 0.62, respectively.
The results were in good agreement with those based on the CV database, indicating
the ensemble model showed satisfactory spatial generalization.</p>
      <?pagebreak page3316?><p id="d1e2523">The ensemble model can capture the spatiotemporal variation in ambient
benzene during the COVID-19 lockdown period, while the impact of the COVID-19
lockdown cannot be quantified because the contribution of meteorological
parameters cannot be removed based on this model alone. Therefore, it is
proposed to employ the XGBoost algorithm to isolate the contribution of
emission reduction to global atmospheric benzene. As depicted in Fig. S5,
the CV <inline-formula><mml:math id="M119" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value and slope of the fitting curve reached 0.65 and 0.62,
respectively. The results  suggest that the meteorology-normalized model was
robust because the CV <inline-formula><mml:math id="M120" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> value was much higher than 0.50.</p>
</sec>
<sec id="Ch1.S3.SS2">
  <label>3.2</label><title>The impact of the COVID-19 lockdown on the global atmospheric benzene level</title>
      <p id="d1e2556">The ensemble model was developed to expand the ground-observed benzene
measurement to the global scale and capture the global spatial variability
in ambient benzene. As shown in Fig. S6, the global simulated (total)
benzene concentration during  23 January–30 June in 2019 and 2020 ranged from
0.52 to 6.36 <inline-formula><mml:math id="M121" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, with an average value of <inline-formula><mml:math id="M122" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.92</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.23</mml:mn></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M123" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>. At the regional scale, the benzene concentration displayed
significant spatial variability. The benzene concentration followed the
order of India (<inline-formula><mml:math id="M124" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.44</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.14</mml:mn></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M125" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) <inline-formula><mml:math id="M126" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China (<inline-formula><mml:math id="M127" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.17</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.13</mml:mn></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M128" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) <inline-formula><mml:math id="M129" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> Europe (<inline-formula><mml:math id="M130" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.02</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.08</mml:mn></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M131" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) <inline-formula><mml:math id="M132" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United States (<inline-formula><mml:math id="M133" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.96</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.09</mml:mn></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M134" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>)
during  23 January–30 June in 2019 and 2020. Besides, the global simulated mean
benzene level underwent a slight decrease from <inline-formula><mml:math id="M135" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.93</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.06</mml:mn></mml:mrow></mml:math></inline-formula> in 2020
to <inline-formula><mml:math id="M136" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.90</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.06</mml:mn></mml:mrow></mml:math></inline-formula> in 2019. However, the inter-annual variation in ambient
benzene exhibited a remarkable spatial discrepancy at the global scale. As
depicted in Fig. S7, the change ratio of the simulated (total) benzene level
during the COVID-19 lockdown period (the difference in the benzene level
before the COVID-19 lockdown and that during the COVID-19 lockdown period) in 2020
was in the order of India (<inline-formula><mml:math id="M137" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">18.5</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M138" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> Europe (<inline-formula><mml:math id="M139" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.7</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M140" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China (<inline-formula><mml:math id="M141" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">11.7</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M142" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United States (<inline-formula><mml:math id="M143" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">11.5</mml:mn></mml:mrow></mml:math></inline-formula> %).
Compared with 2020, the change ratio of the benzene level during the same period
in 2019 followed the order of India (<inline-formula><mml:math id="M144" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.3</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M145" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> Europe
(<inline-formula><mml:math id="M146" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6.62</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M147" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United States (<inline-formula><mml:math id="M148" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6.46</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M149" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China
(<inline-formula><mml:math id="M150" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4.18</mml:mn></mml:mrow></mml:math></inline-formula> %). It should be noted that the simulated ambient benzene
concentration had a higher decreasing ratio in 2020 compared
with the same period in 2019 in nearly all of the major countries and regions around the
world, which might be associated with the local COVID-19 lockdown measures
in 2020.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2" specific-use="star"><?xmltex \currentcnt{2}?><?xmltex \def\figurename{Figure}?><label>Figure 2</label><caption><p id="d1e2907">The global average deweathered benzene concentrations in 2019 (23 January–30 June) <bold>(a)</bold> and 2020 (23 January–30 June) <bold>(b)</bold>. Panel <bold>(c)</bold> represents the
difference in deweathered benzene concentrations in 2020 and 2019 (unit:
<inline-formula><mml:math id="M151" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>).</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/3311/2023/acp-23-3311-2023-f02.png"/>

        </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3" specific-use="star"><?xmltex \currentcnt{3}?><?xmltex \def\figurename{Figure}?><label>Figure 3</label><caption><p id="d1e2946">The weekly variations in atmospheric benzene concentrations (<inline-formula><mml:math id="M152" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) in some major regions around the world during 23 January–30 June.
The red line and background denote mean values and standard deviation of
deweathered weekly benzene concentrations in 2020. The cyan line and
background denote mean values and standard deviation of deweathered weekly
benzene levels in 2019. The dashed vertical red line indicates COVID-19
restriction dates, and the dashed vertical black line indicates the beginning of easing
measures.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/3311/2023/acp-23-3311-2023-f03.png"/>

        </fig>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T1"><?xmltex \currentcnt{1}?><label>Table 1</label><caption><p id="d1e2978">The change ratio (%) of deweathered (<inline-formula><mml:math id="M153" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">dew</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and detrended
(<inline-formula><mml:math id="M154" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>) benzene concentrations in major regions around the world.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="5">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="right"/>
     <oasis:colspec colnum="3" colname="col3" align="right"/>
     <oasis:colspec colnum="4" colname="col4" align="right"/>
     <oasis:colspec colnum="5" colname="col5" align="right"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Change ratio</oasis:entry>
         <oasis:entry colname="col2">China</oasis:entry>
         <oasis:entry colname="col3">India</oasis:entry>
         <oasis:entry colname="col4">Europe</oasis:entry>
         <oasis:entry colname="col5">United States</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M155" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">dew</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in 2020</oasis:entry>
         <oasis:entry colname="col2"><inline-formula><mml:math id="M156" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">15.6</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M157" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.6</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"><inline-formula><mml:math id="M158" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">21.9</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col5"><inline-formula><mml:math id="M159" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M160" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">dew</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in 2019</oasis:entry>
         <oasis:entry colname="col2"><inline-formula><mml:math id="M161" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2.31</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M162" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7.40</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"><inline-formula><mml:math id="M163" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">8.04</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col5"><inline-formula><mml:math id="M164" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">10.2</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M165" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2"><inline-formula><mml:math id="M166" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">13.3</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M167" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"><inline-formula><mml:math id="M168" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">13.9</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col5"><inline-formula><mml:math id="M169" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6.00</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <p id="d1e3222">Due to the interference of meteorological conditions, we cannot quantify the
direct impact of the COVID-19 lockdown on ambient benzene based on the
comparison of simulated (total) benzene levels. Thus, the
meteorology-normalized method was employed to decouple the separated
contributions of emission reduction and meteorology to ambient benzene. In
our study, both the change ratio and the detrended change ratio were applied
to evaluate the impact of the COVID-19 lockdown on the global ambient benzene level.
The change ratio represents the variation in the ambient benzene level during the
lockdown period in 2020 compared with the pre-lockdown period in 2020. However,
the detrended change ratio reflects the difference in the change ratio in
2020 and the change ratio during the same period in 2019, which could avoid
the inter-annual system error and contingency of a single year. As
summarized in Figs. 2 and 3, the change ratio of the deweathered benzene
concentration from the pre-lockdown to lockdown period in 2020 was in the order
of India (<inline-formula><mml:math id="M170" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.6</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M171" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> Europe (<inline-formula><mml:math id="M172" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">21.9</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M173" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United
States (<inline-formula><mml:math id="M174" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M175" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China (<inline-formula><mml:math id="M176" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">15.6</mml:mn></mml:mrow></mml:math></inline-formula> %). Meanwhile, the change
ratio of the deweathered benzene concentration during the same time in 2019
followed the order of Europe (<inline-formula><mml:math id="M177" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">10.2</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M178" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> the United States
(<inline-formula><mml:math id="M179" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">8.04</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M180" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> India (<inline-formula><mml:math id="M181" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7.40</mml:mn></mml:mrow></mml:math></inline-formula> %) <inline-formula><mml:math id="M182" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> China (<inline-formula><mml:math id="M183" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2.31</mml:mn></mml:mrow></mml:math></inline-formula> %).
The large gap in the change ratio of the deweathered benzene level between 2019
and 2020 confirmed that the drastic and consequential quarantine measures
significantly decreased the ambient benzene concentrations in nearly all of
the regions with lockdown measures. Among all of the major countries and regions, India
experienced the most dramatic benzene decrease during 24 March–24 April 2020 (<inline-formula><mml:math id="M184" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.6</mml:mn></mml:mrow></mml:math></inline-formula> %) compared with the same period in 2019 (<inline-formula><mml:math id="M185" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7.4</mml:mn></mml:mrow></mml:math></inline-formula> %).
During this period, the prohibition of industrial activities and mass
transportation was proposed to curb the spread of the COVID-19 pandemic, leading
to a tremendous reduction in anthropogenic benzene emissions (Pathakoti et
al., 2021; Zhang et al., 2021). Sahu et al. (2022) revealed that the
substantial increase in OH radicals during the COVID-19 period also facilitated ambient benzene removal due to the photooxidation reaction. The decreasing
ratio of the deweathered benzene level in India was close to that of PM<inline-formula><mml:math id="M186" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula>
(<inline-formula><mml:math id="M187" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">26</mml:mn></mml:mrow></mml:math></inline-formula> %), while it is was markedly lower than that of NO<inline-formula><mml:math id="M188" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> (<inline-formula><mml:math id="M189" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">50</mml:mn></mml:mrow></mml:math></inline-formula> %)
(Zhang et al., 2021). Although both Europe and the United States also
performed stringent lockdown restrictions in some regions such as Italy,
Spain, and California (Guevara et al., 2021a; Keller et al., 2021a), the detrended change (<inline-formula><mml:math id="M190" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>: change ratio in 2020 <inline-formula><mml:math id="M191" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula> change ratio in 2019) for
deweathered benzene in Europe (<inline-formula><mml:math id="M192" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">13.9</mml:mn></mml:mrow></mml:math></inline-formula> %) and the United States (<inline-formula><mml:math id="M193" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:math></inline-formula> %) between 2020 and 2019 was still lower than that of India (<inline-formula><mml:math id="M194" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula> %) (Table 1). It was assumed that the absolute concentrations of
ambient benzene in Europe and the United States were much lower than that in
India. It should be noted that China displayed a relatively gently
decreasing ratio (<inline-formula><mml:math id="M195" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">15.6</mml:mn></mml:mrow></mml:math></inline-formula> %) after the COVID-19 outbreak, which was even lower
than the ratio in the United States. As the first epidemic epicentre
country, the Chinese government imposed a rapid lockdown measure in Wuhan and
other cities across China in an effort to prevent the spread of the COVID-19
pandemic (Wu et al., 2020). These restrictions interrupted a wide array of
economic activities and reduced primary air pollutant emissions and thus
resulted in remarkable decreases in deweathered NO<inline-formula><mml:math id="M196" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> (<inline-formula><mml:math id="M197" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">43.6</mml:mn></mml:mrow></mml:math></inline-formula> %) and
PM<inline-formula><mml:math id="M198" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> (<inline-formula><mml:math id="M199" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">22</mml:mn></mml:mrow></mml:math></inline-formula> %) (Dai et al., 2021). The gently decreasing ratio of
ambient benzene compared with other pollutants might be linked with the
source apportionment of atmospheric benzene. It was well known that
industrial sources (e.g. the chemical industry and solvent use) comprised the major
emission sector of benzene (Li et al., 2019). Although the contribution from
solvent use exhibited substantial decreases in some cities (Qi et al., 2021;
Wang<?pagebreak page3318?> et al., 2021), the chemical industry was not entirely interrupted even
during the COVID-19 lockdown period (Dai et al., 2021). B. Zheng et al. (2021)
also demonstrated that reduction in NMVOC emissions from the industrial sector was much less than for other
pollutants.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F4" specific-use="star"><?xmltex \currentcnt{4}?><?xmltex \def\figurename{Figure}?><label>Figure 4</label><caption><p id="d1e3529">The concentration difference for deweathered benzene between the
COVID-19 period in 2020 and the same period in 2019 in <bold>(a)</bold> East Asia, <bold>(b)</bold> South
Asia, <bold>(c)</bold> Europe, and <bold>(d)</bold> North America (difference <inline-formula><mml:math id="M200" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> deweathered benzene
concentration in 2020 <inline-formula><mml:math id="M201" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula> deweathered benzene concentration in 2019) (unit: <inline-formula><mml:math id="M202" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi><mml:mi mathvariant="normal">g</mml:mi><mml:mspace width="0.125em" linebreak="nobreak"/><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>).</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/3311/2023/acp-23-3311-2023-f04.png"/>

        </fig>

      <p id="d1e3584">Although the deweathered benzene concentrations in nearly all of the major
countries underwent obvious decreases during the COVID-19 lockdown period, the
change ratios of deweathered benzene in different regions of these countries
still showed large spatial variability. In China, most of the cities in eastern
China such as Beijing (<inline-formula><mml:math id="M203" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">30.6</mml:mn></mml:mrow></mml:math></inline-formula> %), Shanghai (<inline-formula><mml:math id="M204" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6.25</mml:mn></mml:mrow></mml:math></inline-formula> %), and Wuhan
(<inline-formula><mml:math id="M205" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">45.3</mml:mn></mml:mrow></mml:math></inline-formula> %) experienced dramatic decreases in deweathered benzene levels
(Fig. S8), which was mainly contributed by the simultaneous emission
reduction of the industrial and transportation sectors. Besides, enhanced
atmospheric oxidation capacity could have accelerated the benzene removal due to
the unbalanced decreases in VOC and NO<inline-formula><mml:math id="M206" display="inline"><mml:msub><mml:mi/><mml:mi>x</mml:mi></mml:msub></mml:math></inline-formula> emissions (Jensen et al.,
2021). However, the deweathered benzene concentrations in northeastern China
and Yunnan Province even exhibited slight increases after the COVID-19 outbreak
(Fig. 4). Dai et al. (2021) also found that the deweathered PM<inline-formula><mml:math id="M207" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula>
concentration in Kunming increased <inline-formula><mml:math id="M208" display="inline"><mml:mrow><mml:mo>∼</mml:mo><mml:mn mathvariant="normal">20</mml:mn></mml:mrow></mml:math></inline-formula> % after the COVID-19
outbreak. At first, the contribution of residential combustion sources
(62.1 %) to atmospheric benzene in Yunnan Province was higher than that of other
sectors (Guevara et al., 2021b; Kuenen et al., 2022). Moreover, the increase
in domestic emissions due to home quarantine regulations further increased the
ambient benzene concentration (10 %) in this province, which has been
demonstrated by the updated emission inventory in 2020 (Doumbia et al.,
2021). The slight increases in deweathered benzene levels in northeastern China
after the COVID-19 outbreak could be linked with earlier work
resumption (<uri>https://baijiahao.baidu.com/s?id=1658138056285012986&amp;wfr=spider&amp;for=pc</uri>, last access: 15 December 2022).
Based on the simulation result, the deweathered ambient benzene level in
northeastern China rebounded sharply after the third week and then returned
back to normal in late February. In India, the decreasing ratios of
deweathered benzene in Delhi, Mumbai, Kolkata, Bengaluru, Hyderabad,
Chennai, Ahmedabad, and Lucknow reached 21.6 %, 20.9 %, 73.7 %,
26.9 %, 38.0 %, 33.7 %, 25.1 %, and 33.3 %, respectively, during the COVID-19
lockdown period (Fig. S9). Among all of the major cities in
India, the ambient benzene level in Kolkata underwent the most drastic
decrease. It was assumed that Kolkata is a dusty city and filled
with vehicle emissions. Fortunately, the city experienced a complete stop of
vehicle movement and of burning of biomass and dust particles from construction works, which had been important sources for ambient benzene (Bera
et al., 2021; Kumar and Singh, 2003). In Europe, the deweathered benzene
levels in nearly all of the cities displayed marked decreases because most
European countries imposed lockdowns to combat the spread of the
COVID-19 pandemic (Guevara et al., 2021a). For example, private cars
and heavy goods vehicles (HGVs) on the road in London reduced by 80 % and
30 %–40 %, respectively (Shi et al., 2021). The drastic decrease in
transportation emissions triggered the <inline-formula><mml:math id="M209" display="inline"><mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mo>∗</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> value in London between 2020 and
2019 reaching <inline-formula><mml:math id="M210" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">43.6</mml:mn></mml:mrow></mml:math></inline-formula> %. In the United States, the decreasing ratios of
deweathered benzene levels in the cities of the eastern United States and
California were generally higher than those in the central United States, which
was in good agreement with the spatial variability in the PM<inline-formula><mml:math id="M211" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> decrease
(Hammer et al., 2021). It was closely associated with the length of the lockdown
period (<uri>https://en.wikipedia.org/wiki/COVID-19_lockdowns</uri>, last access: 14 January 2023).</p>
      <p id="d1e3682">In addition, ambient benzene levels were also strongly affected by
meteorological conditions that alter photochemical production, advection,
and depositional loss. Hence, we examined how meteorological parameters
influenced the temporal variability in ambient benzene during the COVID-19
lockdown period. In 2020, most of the major countries and regions including China
(3.9 %), Europe (5.2 %), and the United States (4.7 %) experienced
slightly unfavourable meteorological conditions, which is in good agreement
with the impact of meteorological conditions on ambient NO<inline-formula><mml:math id="M212" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula>
concentrations (Shi et al., 2021). Among all the meteorological parameters,
air temperature was the most important factor for ambient benzene in nearly
all of the regions around the world during the study period. Compared with
2019, the air temperatures in China, India, Europe, and the United States
increased by 0.4, 0.9, 0.4, and 0.2 <inline-formula><mml:math id="M213" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C during the same period in
2020, respectively. Jia and Xu (2014) demonstrated that the increased air
temperature generally suppressed benzene photooxidation and secondary
organic aerosol (SOA) formation. Thus, the increased air temperature was not
beneficial to the further reduction in ambient benzene. Apart from air
temperature, some other factors such as RH, rainfall amount, and wind speed
might have affected the ambient benzene level. For instance, the increased RH could have
been favourable to benzene oxidation and the higher rainfall amount could have
promoted benzene removal. However, in the machine-learning model, the
importance values of these variables were much lower than that of air
temperature. Overall, the results suggest that the unfavourable
meteorological conditions (air temperature) weakened the health benefits of reduced
ambient benzene due to drastic lockdown measures around the world.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F5" specific-use="star"><?xmltex \currentcnt{5}?><?xmltex \def\figurename{Figure}?><label>Figure 5</label><caption><p id="d1e3706">The carcinogenic risk differences (unit: 10<inline-formula><mml:math id="M214" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>) for atmospheric
benzene between the COVID-19 period in 2020 and the same period in 2019 in <bold>(a)</bold> East
Asia, <bold>(b)</bold> South Asia, <bold>(c)</bold> Europe, and <bold>(d)</bold> North America (difference <inline-formula><mml:math id="M215" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> benzene
concentration in 2020 <inline-formula><mml:math id="M216" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula> benzene concentration in 2019).</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/3311/2023/acp-23-3311-2023-f05.png"/>

        </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6" specific-use="star"><?xmltex \currentcnt{6}?><?xmltex \def\figurename{Figure}?><label>Figure 6</label><caption><p id="d1e3756">The non-carcinogenic risk differences (unit: 10<inline-formula><mml:math id="M217" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>) for
atmospheric benzene between the COVID-19 period in 2020 and the same period in
2019 in <bold>(a)</bold> East
Asia, <bold>(b)</bold> South Asia, <bold>(c)</bold> Europe, and <bold>(d)</bold> North America (difference <inline-formula><mml:math id="M218" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> benzene concentration in 2020 <inline-formula><mml:math id="M219" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula> benzene concentration in 2019).</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://acp.copernicus.org/articles/23/3311/2023/acp-23-3311-2023-f06.png"/>

        </fig>

</sec>
<sec id="Ch1.S3.SS3">
  <label>3.3</label><title>The effect of the COVID-19 lockdown on global health risks</title>
      <?pagebreak page3320?><p id="d1e3812">The global average LCR during 23 January–30 June in 2019 and 2020 was <inline-formula><mml:math id="M220" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.89</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M221" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.51</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, respectively, after removing the
contributions of meteorological conditions (Fig. S10).
Although the COVID-19 lockdown decreased the LCR value slightly, both of the
LCR values during the two periods were lower than the threshold level of
10<inline-formula><mml:math id="M222" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, suggesting that dwellings in most regions could avoid
carcinogenic risk through inhalation exposure to benzene (Li et al., 2017).
However, the LCR values showed significant spatial difference in different
regions. For instance, northern China often suffered from relatively high
benzene pollution, and the LCR value in this region decreased from <inline-formula><mml:math id="M223" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.03</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> (possible risk) during the study period in 2019 (the
same period to 2020) to <inline-formula><mml:math id="M224" display="inline"><mml:mrow><mml:mn mathvariant="normal">7.37</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> during the COVID-19 lockdown
period. The result verified that the stringent emission control measures
significantly decreased the health risk due to benzene exposure. The LCR
value across India only decreased from <inline-formula><mml:math id="M225" display="inline"><mml:mrow><mml:mn mathvariant="normal">6.55</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M226" display="inline"><mml:mrow><mml:mn mathvariant="normal">6.42</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> during the study period, whereas the value in the northern part of
India, including for example Bihar, decreased from <inline-formula><mml:math id="M227" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.14</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M228" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.09</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> due to the impact of the COVID-19 lockdown (Fig. 5). As
the most populous state of India, Bihar had more than 124 million
people at the time (<uri>http://kolkata.china-consulate.org/chn/lqgk/t1331638.htm</uri>, last access: 31 January 2023). The
results suggest that the COVID-19 lockdown certainly obtained remarkable
short-term health benefits through decreasing ambient benzene exposure.
The LCR values in Europe and the United States decreased from <inline-formula><mml:math id="M229" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.99</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M230" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.77</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M231" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.57</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M232" display="inline"><mml:mrow><mml:mn mathvariant="normal">4.63</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, respectively. Compared with China
and India, Europe and the United States faced relatively low
carcinogenic risk of benzene exposure even before the COVID-19 lockdown.
Although the COVID-19 lockdown further decreased the LCR values in these
regions, the overall carcinogenic risk was negligible.</p>
      <p id="d1e4048">Additionally, the non-carcinogenic risk around the world during the period
was also assessed based on HI. The average HI of ambient benzene exposure in
China, India, Europe, and the United States reduced from <inline-formula><mml:math id="M233" display="inline"><mml:mrow><mml:mn mathvariant="normal">8.92</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M234" display="inline"><mml:mrow><mml:mn mathvariant="normal">7.45</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M235" display="inline"><mml:mrow><mml:mn mathvariant="normal">6.32</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M236" display="inline"><mml:mrow><mml:mn mathvariant="normal">5.76</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> in 2019 to <inline-formula><mml:math id="M237" display="inline"><mml:mrow><mml:mn mathvariant="normal">8.53</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M238" display="inline"><mml:mrow><mml:mn mathvariant="normal">7.13</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M239" display="inline"><mml:mrow><mml:mn mathvariant="normal">5.81</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M240" display="inline"><mml:mrow><mml:mn mathvariant="normal">5.59</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, respectively, during the COVID-19 lockdown period in 2020 (Fig. 6). Although the
HI value in some regions including Bihar (<inline-formula><mml:math id="M241" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.52</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to
<inline-formula><mml:math id="M242" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.41</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) and Uttar Pradesh (<inline-formula><mml:math id="M243" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.04</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to
<inline-formula><mml:math id="M244" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.03</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) in India and Beijing–Tianjin–Hebei (BTH) (<inline-formula><mml:math id="M245" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.25</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M246" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.14</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) in China still
experienced decreases during the COVID-19 lockdown period, the HI values in
these regions were still significantly lower than the risk threshold (HI <inline-formula><mml:math id="M247" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 1). Therefore, the impact of the COVID-19 lockdown on the non-carcinogenic risk of
benzene exposure was insignificant.</p>
</sec>
</sec>
<sec id="Ch1.S4" sec-type="conclusions">
  <label>4</label><title>Conclusions and limitations</title>
      <p id="d1e4321">The drastic lockdown measures largely reduced the air pollutant emissions.
The meteorology-normalized ambient benzene concentrations in China
(<inline-formula><mml:math id="M248" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">15.6</mml:mn></mml:mrow></mml:math></inline-formula> %), India (<inline-formula><mml:math id="M249" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.6</mml:mn></mml:mrow></mml:math></inline-formula> %), Europe (<inline-formula><mml:math id="M250" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">21.9</mml:mn></mml:mrow></mml:math></inline-formula> %), and the United States
(<inline-formula><mml:math id="M251" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">16.2</mml:mn></mml:mrow></mml:math></inline-formula> %) experienced dramatic decreases after the COVID-19 outbreak.
Furthermore, the decreasing ratios in these major regions during the COVID-19
lockdown period were much higher than in the same period in 2019, indicating
the aggressive emission control measures efficiently decreased ambient
benzene<?pagebreak page3321?> concentrations. Emission reductions from industrial activities and
transportation were major drivers of the decrease in the ambient benzene
level during the lockdown period, while the relatively stable solvent use
emissions could have restricted a further decrease in benzene pollution. Besides,
the slight increase in domestic emissions during this period might be an
important reason for the benzene increase in some regions (e.g. Yunnan
Province). There is also an urgent need to control the household combustion
and solvent use emissions in addition to the emissions from the industrial and
transportation sectors.</p>
      <p id="d1e4364">Besides, substantial decreases in atmospheric benzene levels could have
sufficient health benefits. Dramatic decreases in benzene emissions in
Europe and the United States could not have brought about effective health benefits because
the ambient benzene levels in both of these regions during the business-as-usual
scenario were already significantly lower than the risk threshold. However, the
benzene decreases in the North China Plain (NCP), China, and Bihar, India, could
have had abundant health benefits because these regions often suffered from
severe atmospheric benzene pollution during the business-as-usual scenario.
Thus, more targeted abatement measures are needed to reduce the benzene
emissions in these areas. For instance, stricter industrial and vehicle
emission standards for VOC control should be implemented in China and India.
Moreover, some measures including limiting the number of coal-fired power
plants, adding environmentally friendly cars and clean fuels for vehicles
and vessels, and strengthening the labelling system for vehicles in use
should be pursued.</p>
      <p id="d1e4367">It should be noted that our study still suffered from some limitations.
First of all, the monitoring sites were not evenly distributed around the
world, and thus the simulation results might show higher uncertainty in
the regions lacking monitoring sites. Besides, the GEOS-Chem model still
suffered from some uncertainties due to an imperfect chemical mechanism and
inaccurate emission inventory. In future work, the model should be
further improved.</p>
</sec>

      
      </body>
    <back><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e4375">The CEDS emission inventory is available on the website of <ext-link xlink:href="https://doi.org/10.5281/zenodo.3754964" ext-link-type="DOI">10.5281/zenodo.3754964</ext-link> (McDuffie et al., 2020). The global meteorological parameters (reanalysis dataset) can be obtained from the website of are available here: <ext-link xlink:href="https://doi.org/10.5067/QBZ6MG944HW0" ext-link-type="DOI">10.5067/QBZ6MG944HW0</ext-link> (Global Modeling and Assimilation Office, 2015).</p>
  </notes><app-group>
        <supplementary-material position="anchor"><p id="d1e4384">The supplement related to this article is available online at: <inline-supplementary-material xlink:href="https://doi.org/10.5194/acp-23-3311-2023-supplement" xlink:title="pdf">https://doi.org/10.5194/acp-23-3311-2023-supplement</inline-supplementary-material>.</p></supplementary-material>
        </app-group><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e4393">CL and RL wrote the manuscript. RL and LC contributed to the
conceptualization of the study. CL and RL conducted the research and
visualized the results. LC revised the manuscript.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e4399">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d1e4405">Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
  </notes><ack><title>Acknowledgements</title><p id="d1e4411">This work was supported by the National Natural Science Foundation of China
(grant no. 42107113).</p></ack><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d1e4416">This research has been supported by the National Natural Science Foundation of China (grant no. 42107113), the Natural Science Foundation of Fujian Province (grant no. 2022J05178) and the Natural Science Foundation of Hunan Province (grant no. 2021JJ30555).</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e4422">This paper was edited by Zhanqing Li and reviewed by two anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><?label 1?><mixed-citation>Bauwens, M., Compernolle, S., Stavrakou, T., Müller, J. F., van Gent,
J., Eskes, H., Levelt, P. F., van der A, R., Veefkind, J. P., Vlietinck, J.,
Yu, H., and Zehner, C.: Impact of coronavirus outbreak on NO<inline-formula><mml:math id="M252" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> pollution
assessed using TROPOMI and OMI observations, Geophys. Res. Lett., 47,
e2020GL087978, <ext-link xlink:href="https://doi.org/10.1029/2020GL087978" ext-link-type="DOI">10.1029/2020GL087978</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><?label 1?><mixed-citation>
Bera, B., Bhattacharjee, S., Shit, P. K., Sengupta, N., and Saha, S.:
Significant impacts of COVID-19 lockdown on urban air pollution in Kolkata
(India) and amelioration of environmental health, Environ. Dev. Sustain.,
23, 6913–6940, 2021.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><?label 1?><mixed-citation>Berg, K., Romer Present, P., and Richardson, K.: Long-term air pollution and
other risk factors associated with COVID-19 at the census tract level in
Colorado, Environ. Pollut., 287, 117584, <ext-link xlink:href="https://doi.org/10.1016/j.envpol.2021.117584" ext-link-type="DOI">10.1016/j.envpol.2021.117584</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><?label 1?><mixed-citation>Bocquet, M., Elbern, H., Eskes, H., Hirtl, M., Žabkar, R., Carmichael, G. R., Flemming, J., Inness, A., Pagowski, M., Pérez Camaño, J. L., Saide, P. E., San Jose, R., Sofiev, M., Vira, J., Baklanov, A., Carnevale, C., Grell, G., and Seigneur, C.: Data assimilation in atmospheric chemistry models: current status and future prospects for coupled chemistry meteorology models, Atmos. Chem. Phys., 15, 5325–5358, <ext-link xlink:href="https://doi.org/10.5194/acp-15-5325-2015" ext-link-type="DOI">10.5194/acp-15-5325-2015</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><?label 1?><mixed-citation>
Breiman, L.: Random forests, Mach. Learn., 45, 5–32, 2001.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><?label 1?><mixed-citation>Cai, M., Yang, Y. G., Gibilisco, R. G., Grosselin, B., McGillen, M. R., Xue,
C. Y., Mellouki, A., and Daële, V.: Ambient BTEX Concentrations during
the COVID-19 Lockdown in a Peri-Urban Environment (Orléans, France),
Atmosphere, 13, 10, <ext-link xlink:href="https://doi.org/10.3390/atmos13010010" ext-link-type="DOI">10.3390/atmos13010010</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><?label 1?><mixed-citation>Chang, Y., Huang, R. J., Ge, X., Huang, X., Hu, J., Duan, Y., Zou, Z., Liu,
X., and Lehmann, M. F.: Puzzling haze events in China during the coronavirus
(COVID-19) shutdown, Geophys. Res. Lett., 47, e2020GL088533, <ext-link xlink:href="https://doi.org/10.1029/2020GL088533" ext-link-type="DOI">10.1029/2020GL088533</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><?label 1?><mixed-citation>
Cui, L., Li, R., Zhang, Y., Meng, Y., Zhao, Y., and Fu, H.: A geographically
and temporally weighted regression model for assessing intra-urban
variability of volatile organic compounds (VOCs) in Yangpu district,
Shanghai, Atmos. Environ., 213, 746–756, 2019.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><?label 1?><mixed-citation>Dai, Q., Hou, L., Liu, B., Zhang, Y., Song, C., Shi, Z., Hopke, P. K., and
Feng, Y.: Spring Festival and COVID-19 lockdown: disentangling PM sources in
major Chinese cities, Geophys. Res. Lett., 48, e2021GL093403, <ext-link xlink:href="https://doi.org/10.1029/2021GL093403" ext-link-type="DOI">10.1029/2021GL093403</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><?label 1?><mixed-citation>Doumbia, T., Granier, C., Elguindi, N., Bouarar, I., Darras, S., Brasseur, G., Gaubert, B., Liu, Y., Shi, X., Stavrakou, T., Tilmes, S., Lacey, F., Deroubaix, A., and Wang, T.: Changes in global air pollutant emissions during the COVID-19 pandemic: a dataset for atmospheric modeling, Earth Syst. Sci. Data, 13, 4191–4206, <ext-link xlink:href="https://doi.org/10.5194/essd-13-4191-2021" ext-link-type="DOI">10.5194/essd-13-4191-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><?label 1?><mixed-citation>
Dumanoglu, Y., Kara, M., Altiok, H., Odabasi, M., Elbir, T., and Bayram, A.:
Spatial and seasonal variation and source apportionment of volatile organic
compounds (VOCs) in a heavily industrialized region, Atmos. Environ., 98,
168–178, 2014.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><?label 1?><mixed-citation>Geddes, J. A., Martin, R. V., Boys, B. L., and van Donkelaar, A.: Long-term
trends worldwide in ambient NO<inline-formula><mml:math id="M253" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula> concentrations inferred from satellite
observations, Environ. Health Persp., 124, 281–289, 2016.</mixed-citation></ref>
      <ref id="bib1.bib13"><label>13</label><?label 1?><mixed-citation>Global Modeling and Assimilation Office (GMAO): MERRA-2 inst3_3d_asm_Np: 3d,3-Hourly,Instantaneous,Pressure-Level,Assimilation,Assimilated Meteorological Fields V5.12.4, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC) [data set], <ext-link xlink:href="https://doi.org/10.5067/QBZ6MG944HW0" ext-link-type="DOI">10.5067/QBZ6MG944HW0</ext-link> (last access: 28 February 2023), 2015.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><?label 1?><mixed-citation>
Grange, S. K. and Carslaw, D. C.: Using meteorological normalisation to
detect interventions in air quality time series, Sci. Total Environ., 653,
578–588, 2019.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><?label 1?><mixed-citation>Guevara, M., Jorba, O., Soret, A., Petetin, H., Bowdalo, D., Serradell, K., Tena, C., Denier van der Gon, H., Kuenen, J., Peuch, V.-H., and Pérez García-Pando, C.: Time-resolved emission reductions for atmospheric chemistry modelling in Europe during the COVID-19 lockdowns, Atmos. Chem. Phys., 21, 773–797, <ext-link xlink:href="https://doi.org/10.5194/acp-21-773-2021" ext-link-type="DOI">10.5194/acp-21-773-2021</ext-link>, 2021a.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><?label 1?><mixed-citation>Guevara, M., Jorba, O., Tena, C., Denier van der Gon, H., Kuenen, J., Elguindi, N., Darras, S., Granier, C., and Pérez García-Pando, C.: Copernicus Atmosphere Monitoring Service TEMPOral profiles (CAMS-TEMPO): global and European emission temporal profile maps for atmospheric chemistry modelling, Earth Syst. Sci. Data, 13, 367–404, <ext-link xlink:href="https://doi.org/10.5194/essd-13-367-2021" ext-link-type="DOI">10.5194/essd-13-367-2021</ext-link>, 2021b.</mixed-citation></ref>
      <ref id="bib1.bib17"><label>17</label><?label 1?><mixed-citation>
Gurjar, B. R., Ravindra, K., and Nagpure, A. S.: Air pollution trends over
Indian megacities and their local-to-global implications, Atmos. Environ.,
142, 475–495, 2016.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><?label 1?><mixed-citation>Hammer, M. S., van Donkelaar, A., Martin, R. V., McDuffie, E. E., Lyapustin,
A., Sayer, A. M., Hsu, N. C., Levy, R. C., Garay, M. J., and Kalashnikova,
O. V.: Effects of COVID-19 lockdowns on fine particulate matter
concentrations, Sci. Adv., 7, eabg7670, <ext-link xlink:href="https://doi.org/10.1126/sciadv.abg7670" ext-link-type="DOI">10.1126/sciadv.abg7670</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib19"><label>19</label><?label 1?><mixed-citation>
Hsu, C. Y., Chiang, H. C., Shie, R. H., Ku, C. H., Lin, T. Y., Chen, M. J.,
Chen, N. T., and Chen, Y. C.: Ambient VOCs in residential areas near a
large-scale petrochemical complex: Spatiotemporal variation, source
apportionment and health risk, Environ. Pollut., 240, 95–104, 2018.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><?label 1?><mixed-citation>Hu, L., Keller, C. A., Long, M. S., Sherwen, T., Auer, B., Da Silva, A., Nielsen, J. E., Pawson, S., Thompson, M. A., Trayanov, A. L., Travis, K. R., Grange, S. K., Evans, M. J., and Jacob, D. J.: Global simulation of tropospheric chemistry at 12.5 km resolution: performance and evaluation of the GEOS-Chem chemical module (v10-1) within the NASA GEOS Earth system model (GEOS-5 ESM), Geosci. Model Dev., 11, 4603–4620, <ext-link xlink:href="https://doi.org/10.5194/gmd-11-4603-2018" ext-link-type="DOI">10.5194/gmd-11-4603-2018</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib21"><label>21</label><?label 1?><mixed-citation>Huang, X., Ding, A., Gao, J., Zheng, B., Zhou, D., Qi, X., Tang, R., Wang,
J., Ren, C., and Nie, W.: Enhanced secondary pollution offset reduction of
primary emissions during COVID-19 lockdown in China, Natl. Sci. Rev., 8,
nwaa137, <ext-link xlink:href="https://doi.org/10.1093/nsr/nwaa137" ext-link-type="DOI">10.1093/nsr/nwaa137</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib22"><label>22</label><?label 1?><mixed-citation>Ivatt, P. D. and Evans, M. J.: Improving the prediction of an atmospheric chemistry transport model using gradient-boosted regression trees, Atmos. Chem. Phys., 20, 8063–8082, <ext-link xlink:href="https://doi.org/10.5194/acp-20-8063-2020" ext-link-type="DOI">10.5194/acp-20-8063-2020</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib23"><label>23</label><?label 1?><mixed-citation>Jensen, A., Liu, Z. Q., Tan, W., Dix, B., Chen, T. S., Koss, A., Zhu, L.,
Li, L., and Gouw, J. D.: Measurements of volatile organic compounds during the
COVID-19 lockdown in Changzhou, China, Geophys. Res. Lett., 48,
e2021GL095560, <ext-link xlink:href="https://doi.org/10.1029/2021GL095560" ext-link-type="DOI">10.1029/2021GL095560</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib24"><label>24</label><?label 1?><mixed-citation>
Jia, L. and Xu, Y. F.: Effects of relative humidity on ozone and
secondaryorganic aerosol formation from the photooxidationof benzene and
ethylbenzene, Aerosol Sci. Technol., 48, 1–12, 2014.</mixed-citation></ref>
      <ref id="bib1.bib25"><label>25</label><?label 1?><mixed-citation>
Kamal, M. S., Razzak, S. A., and Hossain, M. M.: Catalytic oxidation of
volatile organic compounds (VOCs) – A review, Atmos. Environ., 140, 117–134,
2016.</mixed-citation></ref>
      <ref id="bib1.bib26"><label>26</label><?label 1?><mixed-citation>Keller, C. A., Evans, M. J., Knowland, K. E., Hasenkopf, C. A., Modekurty, S., Lucchesi, R. A., Oda, T., Franca, B. B., Mandarino, F. C., Díaz Suárez, M. V., Ryan, R. G., Fakes, L. H., and Pawson, S.: Global impact of COVID-19 restrictions on the surface concentrations of nitrogen dioxide and ozone, Atmos. Chem. Phys., 21, 3555–3592, <ext-link xlink:href="https://doi.org/10.5194/acp-21-3555-2021" ext-link-type="DOI">10.5194/acp-21-3555-2021</ext-link>, 2021a.</mixed-citation></ref>
      <ref id="bib1.bib27"><label>27</label><?label 1?><mixed-citation>Keller, C. A., Knowland, K. E., Duncan, B. N., Liu, J., Anderson, D. C.,
Das, S., Lucchesi, R. A., Lundgren, E. W., Nicely, J. M., and Nielsen, E.:
Description of the NASA GEOS Composition Forecast Modeling System GEOS-CF
v1. 0, J. Adv. Model Earth Sy., 13, e2020MS002413, <ext-link xlink:href="https://doi.org/10.1029/2020MS002413" ext-link-type="DOI">10.1029/2020MS002413</ext-link>, 2021b.</mixed-citation></ref>
      <ref id="bib1.bib28"><label>28</label><?label 1?><mixed-citation>
Kim, K. H., Jahan, S. A., and Kabir, E.: A review on human health
perspective of air pollution with respect to allergies and asthma, Environ.
Int., 59, 41–52, 2013.</mixed-citation></ref>
      <ref id="bib1.bib29"><label>29</label><?label 1?><mixed-citation>Koppmann, R.: Volatile organic compounds in the atmosphere, John Wiley &amp;
Sons, Blackwell Pub Professional, <uri>https://book.douban.com/subject/2882453/</uri> (last access: 10 March 2023), 2008.</mixed-citation></ref>
      <ref id="bib1.bib30"><label>30</label><?label 1?><mixed-citation>Kuenen, J., Dellaert, S., Visschedijk, A., Jalkanen, J.-P., Super, I., and Denier van der Gon, H.: CAMS-REG-v4: a state-of-the-art high-resolution European emission inventory for air quality modelling, Earth Syst. Sci. Data, 14, 491–515, <ext-link xlink:href="https://doi.org/10.5194/essd-14-491-2022" ext-link-type="DOI">10.5194/essd-14-491-2022</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib31"><label>31</label><?label 1?><mixed-citation>
Kumar, B. and Singh, R. B.: Urban development and anthropogenic climate
change: experience in Indian metropolitan cities, Ltd, New Delhi, India, Manak Publication Pvt,
2003.</mixed-citation></ref>
      <ref id="bib1.bib32"><label>32</label><?label 1?><mixed-citation>
Li, B., Ho, S. S. H., Xue, Y., Huang, Y., Wang, L., Cheng, Y., Dai, W.,
Zhong, H., Cao, J., and Lee, S.: Characterizations of volatil<?pagebreak page3323?>e organic
compounds (VOCs) from vehicular emissions at roadside environment: The first
comprehensive study in Northwestern China, Atmos. Environ., 161, 1–12, 2017.</mixed-citation></ref>
      <ref id="bib1.bib33"><label>33</label><?label 1?><mixed-citation>Li, M., Zhang, Q., Zheng, B., Tong, D., Lei, Y., Liu, F., Hong, C., Kang, S., Yan, L., Zhang, Y., Bo, Y., Su, H., Cheng, Y., and He, K.: Persistent growth of anthropogenic non-methane volatile organic compound (NMVOC) emissions in China during 1990–2017: drivers, speciation and ozone formation potential, Atmos. Chem. Phys., 19, 8897–8913, <ext-link xlink:href="https://doi.org/10.5194/acp-19-8897-2019" ext-link-type="DOI">10.5194/acp-19-8897-2019</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib34"><label>34</label><?label 1?><mixed-citation>Li, R., Cui, L., Fu, H., Zhao, Y., Zhou, W., and Chen, J.: Satellite-based
estimates of wet ammonium (NH<inline-formula><mml:math id="M254" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">4</mml:mn></mml:msub></mml:math></inline-formula>-N) deposition fluxes across China during
2011–2016 using a space-time ensemble model, Environ. Sci. Technol., 54,
13419–13428, 2020.</mixed-citation></ref>
      <ref id="bib1.bib35"><label>35</label><?label 1?><mixed-citation>Li, R., Cui, L., Zhao, Y., Zhou, W., and Fu, H.: Long-term trends of ambient nitrate (NO3−) concentrations across China based on ensemble machine-learning models, Earth Syst. Sci. Data, 13, 2147–2163, <ext-link xlink:href="https://doi.org/10.5194/essd-13-2147-2021" ext-link-type="DOI">10.5194/essd-13-2147-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib36"><label>36</label><?label 1?><mixed-citation>Ling, C. and Li, Y.: Substantial changes of gaseous pollutants and health
effects during the COVID-19 lockdown period across China, GeoHealth, 5,
e2021GH000408, <ext-link xlink:href="https://doi.org/10.1029/2021GH000408" ext-link-type="DOI">10.1029/2021GH000408</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib37"><label>37</label><?label 1?><mixed-citation>Liu, H., Gong, P., Wang, J., Clinton, N., Bai, Y., and Liang, S.: Annual dynamics of global land cover and its long-term changes from 1982 to 2015, Earth Syst. Sci. Data, 12, 1217–1243, <ext-link xlink:href="https://doi.org/10.5194/essd-12-1217-2020" ext-link-type="DOI">10.5194/essd-12-1217-2020</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib38"><label>38</label><?label 1?><mixed-citation>Long, M. S., Yantosca, R., Nielsen, J. E., Keller, C. A., da Silva, A., Sulprizio, M. P., Pawson, S., and Jacob, D. J.: Development of a grid-independent GEOS-Chem chemical transport model (v9-02) as an atmospheric chemistry module for Earth system models, Geosci. Model Dev., 8, 595–602, <ext-link xlink:href="https://doi.org/10.5194/gmd-8-595-2015" ext-link-type="DOI">10.5194/gmd-8-595-2015</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib39"><label>39</label><?label 1?><mixed-citation>Lu, F., Li, S., Shen, B., Zhang, J., Liu, L., Shen, X., and Zhao, R.: The
emission characteristic of VOCs and the toxicity of BTEX from different
mosquito-repellent incenses, J. Hazard. Mater., 384, 121428, <ext-link xlink:href="https://doi.org/10.1016/j.jhazmat.2019.121428" ext-link-type="DOI">10.1016/j.jhazmat.2019.121428</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib40"><label>40</label><?label 1?><mixed-citation>Mahato, S., Pal, S., and Ghosh, K. G.: Effect of lockdown amid COVID-19
pandemic on air quality of the megacity Delhi, India, Sci. Total Environ.,
730, 139086, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2020.139086" ext-link-type="DOI">10.1016/j.scitotenv.2020.139086</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib41"><label>41</label><?label 1?><mixed-citation>McDuffie, E., Smith, S., O'Rourke, P., Tibrewal, K., Venkataraman, C., Marais, E., Zheng, B., Crippa, M., Brauer, M., and Martin, R.: CEDS_GBD-MAPS: Global Anthropogenic Emission Inventory of NO<inline-formula><mml:math id="M255" display="inline"><mml:msub><mml:mi/><mml:mi>x</mml:mi></mml:msub></mml:math></inline-formula>, SO<inline-formula><mml:math id="M256" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula>, CO, NH<inline-formula><mml:math id="M257" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:math></inline-formula>, NMVOCs, BC, and OC from 1970–2017 (2020_v1.0), Zenodo [data set], <ext-link xlink:href="https://doi.org/10.5281/zenodo.3754964" ext-link-type="DOI">10.5281/zenodo.3754964</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib42"><label>42</label><?label 1?><mixed-citation>Mor, S., Kumar, S., Singh, T., Dogra, S., Pandey, V., and Ravindra, K.:
Impact of COVID-19 lockdown on air quality in Chandigarh, India:
understanding the emission sources during controlled anthropogenic
activities, Chemosphere, 263, 127978, <ext-link xlink:href="https://doi.org/10.1016/j.chemosphere.2020.127978" ext-link-type="DOI">10.1016/j.chemosphere.2020.127978</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib43"><label>43</label><?label 1?><mixed-citation>Mozaffar, A. and Zhang, Y. L.: Atmospheric volatile organic compounds
(VOCs) in China: A review, Curr. Pollut. Rep., 31, 1–14, <ext-link xlink:href="https://doi.org/10.1007/s40726-020-00149-1" ext-link-type="DOI">10.1007/s40726-020-00149-1</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib44"><label>44</label><?label 1?><mixed-citation>
Nemeth, M., Borkin, D., and Michalconok, G.: The comparison of
machine-learning methods XGBoost and LightGBM to predict energy development,
Proceedings of the Computational Methods in Systems and Software, Springer,
208–215, 2019.</mixed-citation></ref>
      <ref id="bib1.bib45"><label>45</label><?label 1?><mixed-citation>Pakkattil, A., Muhsin, M., and Varma, M. R.: COVID-19 lockdown: Effects on
selected volatile organic compound (VOC) emissions over the major Indian
metro cities, Urban Clim., 37, 100838, <ext-link xlink:href="https://doi.org/10.1016/j.uclim.2021.100838" ext-link-type="DOI">10.1016/j.uclim.2021.100838</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib46"><label>46</label><?label 1?><mixed-citation>Pathakoti, M., Muppalla, A., Hazra, S., D. Venkata, M., A. Lakshmi, K., K. Sagar, V., Shekhar, R., Jella, S., M. V. Rama, S. S., and Vijayasundaram, U.: Measurement report: An assessment of the impact of a nationwide lockdown on air pollution – a remote sensing perspective over India, Atmos. Chem. Phys., 21, 9047–9064, <ext-link xlink:href="https://doi.org/10.5194/acp-21-9047-2021" ext-link-type="DOI">10.5194/acp-21-9047-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib47"><label>47</label><?label 1?><mixed-citation>Pei, C. L., Yang, W. Q., Zhang, Y. L., Song, W., Xiao, S. X., Wang, J.,
Zhang, J. P., Zhang, T., Chen, D. H., Wang, Y. J., Chen, Y. N., and Wang, X.
M.: Decrease in ambient volatile organic compounds during the COVID-19
lockdown period in the Pearl River Delta region, South China, Sci. Total
Environ., 823, 153720, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2022.153720" ext-link-type="DOI">10.1016/j.scitotenv.2022.153720</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib48"><label>48</label><?label 1?><mixed-citation>Qi, J., Mo, Z., Yuan, B., Huang, S., Huangfu, Y., Wang, Z., Li, X., Yang,
S., Wang, W., and Zhao, Y.: An observation approach in evaluation of ozone
production to precursor changes during the COVID-19 lockdown, Atmos.
Environ., 262, 118618, <ext-link xlink:href="https://doi.org/10.1016/j.atmosenv.2021.118618" ext-link-type="DOI">10.1016/j.atmosenv.2021.118618</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib49"><label>49</label><?label 1?><mixed-citation>Sahu, L. K., Tripathi, N., Gupta, M., Singh, V., Yadav, R., and Patel, K.:
Impact of COVID-19 pandemic lockdown in ambient concentrations of aromatic
volatile organic compounds in a metropolitan city of Western India, J.
Geophy. Res.-Atmos., 127, e2022JD036628, <ext-link xlink:href="https://doi.org/10.1029/2022JD036628" ext-link-type="DOI">10.1029/2022JD036628</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib50"><label>50</label><?label 1?><mixed-citation>Sharma, S., Zhang, M., Anshika, Gao, J., Zhang, H., and Kota, S. H.: Effect
of restricted emissions during COVID-19 on air quality in India, Sci. Total
Environ., 728, 138878, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2020.138878" ext-link-type="DOI">10.1016/j.scitotenv.2020.138878</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib51"><label>51</label><?label 1?><mixed-citation>Shi, X., Zhao, C., Jiang, J.H., Wang, C., Yang, X., and Yung, Y. L.: Spatial
representativeness of PM<inline-formula><mml:math id="M258" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> concentrations obtained using observations
from network stations, J. Geophy. Res.-Atmos., 123, 3145–3158, 2018.</mixed-citation></ref>
      <ref id="bib1.bib52"><label>52</label><?label 1?><mixed-citation>Shi, Z., Song, C., Liu, B., Lu, G., Xu, J., Van Vu, T., Elliott, R. J., Li,
W., Bloss, W. J., and Harrison, R. M.: Abrupt but smaller than expected
changes in surface air quality attributable to COVID-19 lockdowns, Sci.
Adv., 7, eabd6696, <ext-link xlink:href="https://doi.org/10.1126/sciadv.abd6696" ext-link-type="DOI">10.1126/sciadv.abd6696</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib53"><label>53</label><?label 1?><mixed-citation>Sun, X., Liu, M., and Sima, Z.: A novel cryptocurrency price trend
forecasting model based on LightGBM, Financ. Res. Lett., 32, 101084, <ext-link xlink:href="https://doi.org/10.1016/j.frl.2018.12.032" ext-link-type="DOI">10.1016/j.frl.2018.12.032</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib54"><label>54</label><?label 1?><mixed-citation>Tang, J., Chan, L., Chan, C., Li, Y. S., Chang, C., Liu, S., Wu, D., and Li,
Y.: Characteristics and diurnal variations of NMHCs at urban, suburban, and
rural sites in the Pearl River Delta and a remote site in South China,
Atmos. Environ., 41, 8620–8632, <ext-link xlink:href="https://doi.org/10.1016/j.atmosenv.2007.07.029" ext-link-type="DOI">10.1016/j.atmosenv.2007.07.029</ext-link>, 2007.</mixed-citation></ref>
      <ref id="bib1.bib55"><label>55</label><?label 1?><mixed-citation>
Van Donkelaar, A., Martin, R. V., Brauer, M., and Boys, B. L.: Use of
satellite observations for long-term exposure assessment of global
concentrations of fine particulate matter, Environ. Health Persp., 123,
135–143, 2015.</mixed-citation></ref>
      <ref id="bib1.bib56"><label>56</label><?label 1?><mixed-citation>Wang, M., Lu, S., Shao, M., Zeng, L., Zheng, J., Xie, F., Lin, H., Hu, K.,
and Lu, X.: Impact of COVID-19 lockdown on ambient levels and sources of
volatile organic compounds (VOCs) in Nanjing, China, Sci. Total Environ.,
757, 143823, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2020.143823" ext-link-type="DOI">10.1016/j.scitotenv.2020.143823</ext-link>, 2021.</mixed-citation></ref>
      <?pagebreak page3324?><ref id="bib1.bib57"><label>57</label><?label 1?><mixed-citation>Wei, J., Li, Z., Lyapustin, A., Sun, L., Peng, Y., Xue, W., Su, T., and
Cribb, M.: Reconstructing 1-km-resolution high-quality PM<inline-formula><mml:math id="M259" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> data
records from 2000 to 2018 in China: spatiotemporal variations and policy
implications, Remote Sens Environ., 252, 112136, <ext-link xlink:href="https://doi.org/10.1016/j.rse.2020.112136" ext-link-type="DOI">10.1016/j.rse.2020.112136</ext-link>, 2021a.</mixed-citation></ref>
      <ref id="bib1.bib58"><label>58</label><?label 1?><mixed-citation>Wei, J., Li, Z., Pinker, R. T., Wang, J., Sun, L., Xue, W., Li, R., and Cribb, M.: Himawari-8-derived diurnal variations in ground-level PM2.5 pollution across China using the fast space-time Light Gradient Boosting Machine (LightGBM), Atmos. Chem. Phys., 21, 7863–7880, <ext-link xlink:href="https://doi.org/10.5194/acp-21-7863-2021" ext-link-type="DOI">10.5194/acp-21-7863-2021</ext-link>, 2021b.</mixed-citation></ref>
      <ref id="bib1.bib59"><label>59</label><?label 1?><mixed-citation>Wei, J., Li, Z., Li, K., Dickerson, R., Pinker, R., Wang, J., Liu, X., Sun,
L., Xue, W., and Cribb, M.: Full-coverage mapping and spatiotemporal
variations of ground-level ozone (O<inline-formula><mml:math id="M260" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:math></inline-formula>) pollution from 2013 to 2020
across China, Remote Sens Environ., 270, 112775, <ext-link xlink:href="https://doi.org/10.1016/j.rse.2021.112775" ext-link-type="DOI">10.1016/j.rse.2021.112775</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib60"><label>60</label><?label 1?><mixed-citation>Wei, J., Li, Z., Wang, J., Li, C., Gupta, P., and Cribb, M.: Ground-level gaseous pollutants (NO<inline-formula><mml:math id="M261" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula>, SO<inline-formula><mml:math id="M262" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:math></inline-formula>, and CO) in China: daily seamless mapping and spatiotemporal variations, Atmos. Chem. Phys., 23, 1511–1532, <ext-link xlink:href="https://doi.org/10.5194/acp-23-1511-2023" ext-link-type="DOI">10.5194/acp-23-1511-2023</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib61"><label>61</label><?label 1?><mixed-citation>
Wolpert, D. H.: Stacked generalization, Neural Networks 5, 241–259, 1992.</mixed-citation></ref>
      <ref id="bib1.bib62"><label>62</label><?label 1?><mixed-citation>Wu, J., Gamber, M., and Sun, W.: Does wuhan need to be in lock-down during the Chinese lunar new year?, Int. J. Env. Res. Pub. He., 17, 3, <ext-link xlink:href="https://doi.org/10.3390/ijerph17031002" ext-link-type="DOI">10.3390/ijerph17031002</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib63"><label>63</label><?label 1?><mixed-citation>Xue, T., Zheng, Y., Tong, D., Zheng, B., Li, X., Zhu, T., and Zhang, Q.:
Spatiotemporal continuous estimates of PM<inline-formula><mml:math id="M263" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> concentrations in China,
2000–2016: A machine learning method with inputs from satellites, chemical
transport model, and ground observations, Environ. Int., 123, 345–357, 2019.
</mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bib64"><label>64</label><?label 1?><mixed-citation>Zhai, B. and Chen, J.: Development of a stacked ensemble model for
forecasting and analyzing daily average PM<inline-formula><mml:math id="M264" display="inline"><mml:msub><mml:mi/><mml:mn mathvariant="normal">2.5</mml:mn></mml:msub></mml:math></inline-formula> concentrations in
Beijing, China, Sci. Total Environ., 635, 644–658, 2018.</mixed-citation></ref>
      <ref id="bib1.bib65"><label>65</label><?label 1?><mixed-citation>Zhang, M., Katiyar, A., Zhu, S., Shen, J., Xia, M., Ma, J., Kota, S. H., Wang, P., and Zhang, H.: Impact of reduced anthropogenic emissions during COVID-19 on air quality in India, Atmos. Chem. Phys., 21, 4025–4037, <ext-link xlink:href="https://doi.org/10.5194/acp-21-4025-2021" ext-link-type="DOI">10.5194/acp-21-4025-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib66"><label>66</label><?label 1?><mixed-citation>
Zhao, Y. B., Zhang, K., Xu, X. T., Shen, H. Z., Zhu, X., Zhang, Y. X., Hu, Y. T., and Shen, G. H.: Substantial changes in nitrogen dioxide and ozone after excluding meteorological impacts during the Covid-19 outbreak in mainland China. Environ. Sci. Technol., 7, 402–408, 2020.</mixed-citation></ref>
      <ref id="bib1.bib67"><label>67</label><?label 1?><mixed-citation>Zheng, B., Zhang, Q., Geng, G., Chen, C., Shi, Q., Cui, M., Lei, Y., and He, K.: Changes in China's anthropogenic emissions and air quality during the COVID-19 pandemic in 2020, Earth Syst. Sci. Data, 13, 2895–2907, <ext-link xlink:href="https://doi.org/10.5194/essd-13-2895-2021" ext-link-type="DOI">10.5194/essd-13-2895-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib68"><label>68</label><?label 1?><mixed-citation>Zheng, P., Chen, Z., Liu, Y., Song, H., Wu, C. H., Li, B., Kraemer, M. U.
G., Tian, H., Yan, X., Zheng, Y., Stenseth, N. C., and Jia, G.: Association
between coronavirus disease 2019 (COVID-19) and long-term exposure to air
pollution: Evidence from the first epidemic wave in China, Environ. Pollut.,
276, 116682, <ext-link xlink:href="https://doi.org/10.1016/j.envpol.2021.116682" ext-link-type="DOI">10.1016/j.envpol.2021.116682</ext-link>, 2021.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Global impact of the COVID-19 lockdown on surface concentration and health risk of atmospheric benzene</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
      
Bauwens, M., Compernolle, S., Stavrakou, T., Müller, J. F., van Gent,
J., Eskes, H., Levelt, P. F., van der A, R., Veefkind, J. P., Vlietinck, J.,
Yu, H., and Zehner, C.: Impact of coronavirus outbreak on NO<sub>2</sub> pollution
assessed using TROPOMI and OMI observations, Geophys. Res. Lett., 47,
e2020GL087978, <a href="https://doi.org/10.1029/2020GL087978" target="_blank">https://doi.org/10.1029/2020GL087978</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
      
Bera, B., Bhattacharjee, S., Shit, P. K., Sengupta, N., and Saha, S.:
Significant impacts of COVID-19 lockdown on urban air pollution in Kolkata
(India) and amelioration of environmental health, Environ. Dev. Sustain.,
23, 6913–6940, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
      
Berg, K., Romer Present, P., and Richardson, K.: Long-term air pollution and
other risk factors associated with COVID-19 at the census tract level in
Colorado, Environ. Pollut., 287, 117584, <a href="https://doi.org/10.1016/j.envpol.2021.117584" target="_blank">https://doi.org/10.1016/j.envpol.2021.117584</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
      
Bocquet, M., Elbern, H., Eskes, H., Hirtl, M., Žabkar, R., Carmichael, G. R., Flemming, J., Inness, A., Pagowski, M., Pérez Camaño, J. L., Saide, P. E., San Jose, R., Sofiev, M., Vira, J., Baklanov, A., Carnevale, C., Grell, G., and Seigneur, C.: Data assimilation in atmospheric chemistry models: current status and future prospects for coupled chemistry meteorology models, Atmos. Chem. Phys., 15, 5325–5358, <a href="https://doi.org/10.5194/acp-15-5325-2015" target="_blank">https://doi.org/10.5194/acp-15-5325-2015</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
      
Breiman, L.: Random forests, Mach. Learn., 45, 5–32, 2001.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
      
Cai, M., Yang, Y. G., Gibilisco, R. G., Grosselin, B., McGillen, M. R., Xue,
C. Y., Mellouki, A., and Daële, V.: Ambient BTEX Concentrations during
the COVID-19 Lockdown in a Peri-Urban Environment (Orléans, France),
Atmosphere, 13, 10, <a href="https://doi.org/10.3390/atmos13010010" target="_blank">https://doi.org/10.3390/atmos13010010</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
      
Chang, Y., Huang, R. J., Ge, X., Huang, X., Hu, J., Duan, Y., Zou, Z., Liu,
X., and Lehmann, M. F.: Puzzling haze events in China during the coronavirus
(COVID-19) shutdown, Geophys. Res. Lett., 47, e2020GL088533, <a href="https://doi.org/10.1029/2020GL088533" target="_blank">https://doi.org/10.1029/2020GL088533</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
      
Cui, L., Li, R., Zhang, Y., Meng, Y., Zhao, Y., and Fu, H.: A geographically
and temporally weighted regression model for assessing intra-urban
variability of volatile organic compounds (VOCs) in Yangpu district,
Shanghai, Atmos. Environ., 213, 746–756, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
      
Dai, Q., Hou, L., Liu, B., Zhang, Y., Song, C., Shi, Z., Hopke, P. K., and
Feng, Y.: Spring Festival and COVID-19 lockdown: disentangling PM sources in
major Chinese cities, Geophys. Res. Lett., 48, e2021GL093403, <a href="https://doi.org/10.1029/2021GL093403" target="_blank">https://doi.org/10.1029/2021GL093403</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
      
Doumbia, T., Granier, C., Elguindi, N., Bouarar, I., Darras, S., Brasseur, G., Gaubert, B., Liu, Y., Shi, X., Stavrakou, T., Tilmes, S., Lacey, F., Deroubaix, A., and Wang, T.: Changes in global air pollutant emissions during the COVID-19 pandemic: a dataset for atmospheric modeling, Earth Syst. Sci. Data, 13, 4191–4206, <a href="https://doi.org/10.5194/essd-13-4191-2021" target="_blank">https://doi.org/10.5194/essd-13-4191-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
      
Dumanoglu, Y., Kara, M., Altiok, H., Odabasi, M., Elbir, T., and Bayram, A.:
Spatial and seasonal variation and source apportionment of volatile organic
compounds (VOCs) in a heavily industrialized region, Atmos. Environ., 98,
168–178, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
      
Geddes, J. A., Martin, R. V., Boys, B. L., and van Donkelaar, A.: Long-term
trends worldwide in ambient NO<sub>2</sub> concentrations inferred from satellite
observations, Environ. Health Persp., 124, 281–289, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
      
Global Modeling and Assimilation Office (GMAO): MERRA-2 inst3_3d_asm_Np: 3d,3-Hourly,Instantaneous,Pressure-Level,Assimilation,Assimilated Meteorological Fields V5.12.4, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC) [data set], <a href="https://doi.org/10.5067/QBZ6MG944HW0" target="_blank">https://doi.org/10.5067/QBZ6MG944HW0</a> (last access: 28 February 2023), 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
      
Grange, S. K. and Carslaw, D. C.: Using meteorological normalisation to
detect interventions in air quality time series, Sci. Total Environ., 653,
578–588, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
      
Guevara, M., Jorba, O., Soret, A., Petetin, H., Bowdalo, D., Serradell, K., Tena, C., Denier van der Gon, H., Kuenen, J., Peuch, V.-H., and Pérez García-Pando, C.: Time-resolved emission reductions for atmospheric chemistry modelling in Europe during the COVID-19 lockdowns, Atmos. Chem. Phys., 21, 773–797, <a href="https://doi.org/10.5194/acp-21-773-2021" target="_blank">https://doi.org/10.5194/acp-21-773-2021</a>, 2021a.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
      
Guevara, M., Jorba, O., Tena, C., Denier van der Gon, H., Kuenen, J., Elguindi, N., Darras, S., Granier, C., and Pérez García-Pando, C.: Copernicus Atmosphere Monitoring Service TEMPOral profiles (CAMS-TEMPO): global and European emission temporal profile maps for atmospheric chemistry modelling, Earth Syst. Sci. Data, 13, 367–404, <a href="https://doi.org/10.5194/essd-13-367-2021" target="_blank">https://doi.org/10.5194/essd-13-367-2021</a>, 2021b.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
      
Gurjar, B. R., Ravindra, K., and Nagpure, A. S.: Air pollution trends over
Indian megacities and their local-to-global implications, Atmos. Environ.,
142, 475–495, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
      
Hammer, M. S., van Donkelaar, A., Martin, R. V., McDuffie, E. E., Lyapustin,
A., Sayer, A. M., Hsu, N. C., Levy, R. C., Garay, M. J., and Kalashnikova,
O. V.: Effects of COVID-19 lockdowns on fine particulate matter
concentrations, Sci. Adv., 7, eabg7670, <a href="https://doi.org/10.1126/sciadv.abg7670" target="_blank">https://doi.org/10.1126/sciadv.abg7670</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
      
Hsu, C. Y., Chiang, H. C., Shie, R. H., Ku, C. H., Lin, T. Y., Chen, M. J.,
Chen, N. T., and Chen, Y. C.: Ambient VOCs in residential areas near a
large-scale petrochemical complex: Spatiotemporal variation, source
apportionment and health risk, Environ. Pollut., 240, 95–104, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
      
Hu, L., Keller, C. A., Long, M. S., Sherwen, T., Auer, B., Da Silva, A., Nielsen, J. E., Pawson, S., Thompson, M. A., Trayanov, A. L., Travis, K. R., Grange, S. K., Evans, M. J., and Jacob, D. J.: Global simulation of tropospheric chemistry at 12.5&thinsp;km resolution: performance and evaluation of the GEOS-Chem chemical module (v10-1) within the NASA GEOS Earth system model (GEOS-5 ESM), Geosci. Model Dev., 11, 4603–4620, <a href="https://doi.org/10.5194/gmd-11-4603-2018" target="_blank">https://doi.org/10.5194/gmd-11-4603-2018</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>21</label><mixed-citation>
      
Huang, X., Ding, A., Gao, J., Zheng, B., Zhou, D., Qi, X., Tang, R., Wang,
J., Ren, C., and Nie, W.: Enhanced secondary pollution offset reduction of
primary emissions during COVID-19 lockdown in China, Natl. Sci. Rev., 8,
nwaa137, <a href="https://doi.org/10.1093/nsr/nwaa137" target="_blank">https://doi.org/10.1093/nsr/nwaa137</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>22</label><mixed-citation>
      
Ivatt, P. D. and Evans, M. J.: Improving the prediction of an atmospheric chemistry transport model using gradient-boosted regression trees, Atmos. Chem. Phys., 20, 8063–8082, <a href="https://doi.org/10.5194/acp-20-8063-2020" target="_blank">https://doi.org/10.5194/acp-20-8063-2020</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>23</label><mixed-citation>
      
Jensen, A., Liu, Z. Q., Tan, W., Dix, B., Chen, T. S., Koss, A., Zhu, L.,
Li, L., and Gouw, J. D.: Measurements of volatile organic compounds during the
COVID-19 lockdown in Changzhou, China, Geophys. Res. Lett., 48,
e2021GL095560, <a href="https://doi.org/10.1029/2021GL095560" target="_blank">https://doi.org/10.1029/2021GL095560</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>24</label><mixed-citation>
      
Jia, L. and Xu, Y. F.: Effects of relative humidity on ozone and
secondaryorganic aerosol formation from the photooxidationof benzene and
ethylbenzene, Aerosol Sci. Technol., 48, 1–12, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>25</label><mixed-citation>
      
Kamal, M. S., Razzak, S. A., and Hossain, M. M.: Catalytic oxidation of
volatile organic compounds (VOCs) – A review, Atmos. Environ., 140, 117–134,
2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>26</label><mixed-citation>
      
Keller, C. A., Evans, M. J., Knowland, K. E., Hasenkopf, C. A., Modekurty, S., Lucchesi, R. A., Oda, T., Franca, B. B., Mandarino, F. C., Díaz Suárez, M. V., Ryan, R. G., Fakes, L. H., and Pawson, S.: Global impact of COVID-19 restrictions on the surface concentrations of nitrogen dioxide and ozone, Atmos. Chem. Phys., 21, 3555–3592, <a href="https://doi.org/10.5194/acp-21-3555-2021" target="_blank">https://doi.org/10.5194/acp-21-3555-2021</a>, 2021a.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>27</label><mixed-citation>
      
Keller, C. A., Knowland, K. E., Duncan, B. N., Liu, J., Anderson, D. C.,
Das, S., Lucchesi, R. A., Lundgren, E. W., Nicely, J. M., and Nielsen, E.:
Description of the NASA GEOS Composition Forecast Modeling System GEOS-CF
v1. 0, J. Adv. Model Earth Sy., 13, e2020MS002413, <a href="https://doi.org/10.1029/2020MS002413" target="_blank">https://doi.org/10.1029/2020MS002413</a>, 2021b.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>28</label><mixed-citation>
      
Kim, K. H., Jahan, S. A., and Kabir, E.: A review on human health
perspective of air pollution with respect to allergies and asthma, Environ.
Int., 59, 41–52, 2013.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>29</label><mixed-citation>
      
Koppmann, R.: Volatile organic compounds in the atmosphere, John Wiley &amp;
Sons, Blackwell Pub Professional, <a href="https://book.douban.com/subject/2882453/" target="_blank"/> (last access: 10 March 2023), 2008.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>30</label><mixed-citation>
      
Kuenen, J., Dellaert, S., Visschedijk, A., Jalkanen, J.-P., Super, I., and Denier van der Gon, H.: CAMS-REG-v4: a state-of-the-art high-resolution European emission inventory for air quality modelling, Earth Syst. Sci. Data, 14, 491–515, <a href="https://doi.org/10.5194/essd-14-491-2022" target="_blank">https://doi.org/10.5194/essd-14-491-2022</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib31"><label>31</label><mixed-citation>
      
Kumar, B. and Singh, R. B.: Urban development and anthropogenic climate
change: experience in Indian metropolitan cities, Ltd, New Delhi, India, Manak Publication Pvt,
2003.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib32"><label>32</label><mixed-citation>
      
Li, B., Ho, S. S. H., Xue, Y., Huang, Y., Wang, L., Cheng, Y., Dai, W.,
Zhong, H., Cao, J., and Lee, S.: Characterizations of volatile organic
compounds (VOCs) from vehicular emissions at roadside environment: The first
comprehensive study in Northwestern China, Atmos. Environ., 161, 1–12, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib33"><label>33</label><mixed-citation>
      
Li, M., Zhang, Q., Zheng, B., Tong, D., Lei, Y., Liu, F., Hong, C., Kang, S., Yan, L., Zhang, Y., Bo, Y., Su, H., Cheng, Y., and He, K.: Persistent growth of anthropogenic non-methane volatile organic compound (NMVOC) emissions in China during 1990–2017: drivers, speciation and ozone formation potential, Atmos. Chem. Phys., 19, 8897–8913, <a href="https://doi.org/10.5194/acp-19-8897-2019" target="_blank">https://doi.org/10.5194/acp-19-8897-2019</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib34"><label>34</label><mixed-citation>
      
Li, R., Cui, L., Fu, H., Zhao, Y., Zhou, W., and Chen, J.: Satellite-based
estimates of wet ammonium (NH<sub>4</sub>-N) deposition fluxes across China during
2011–2016 using a space-time ensemble model, Environ. Sci. Technol., 54,
13419–13428, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib35"><label>35</label><mixed-citation>
      
Li, R., Cui, L., Zhao, Y., Zhou, W., and Fu, H.: Long-term trends of ambient nitrate (NO3−) concentrations across China based on ensemble machine-learning models, Earth Syst. Sci. Data, 13, 2147–2163, <a href="https://doi.org/10.5194/essd-13-2147-2021" target="_blank">https://doi.org/10.5194/essd-13-2147-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib36"><label>36</label><mixed-citation>
      
Ling, C. and Li, Y.: Substantial changes of gaseous pollutants and health
effects during the COVID-19 lockdown period across China, GeoHealth, 5,
e2021GH000408, <a href="https://doi.org/10.1029/2021GH000408" target="_blank">https://doi.org/10.1029/2021GH000408</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib37"><label>37</label><mixed-citation>
      
Liu, H., Gong, P., Wang, J., Clinton, N., Bai, Y., and Liang, S.: Annual dynamics of global land cover and its long-term changes from 1982 to 2015, Earth Syst. Sci. Data, 12, 1217–1243, <a href="https://doi.org/10.5194/essd-12-1217-2020" target="_blank">https://doi.org/10.5194/essd-12-1217-2020</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib38"><label>38</label><mixed-citation>
      
Long, M. S., Yantosca, R., Nielsen, J. E., Keller, C. A., da Silva, A., Sulprizio, M. P., Pawson, S., and Jacob, D. J.: Development of a grid-independent GEOS-Chem chemical transport model (v9-02) as an atmospheric chemistry module for Earth system models, Geosci. Model Dev., 8, 595–602, <a href="https://doi.org/10.5194/gmd-8-595-2015" target="_blank">https://doi.org/10.5194/gmd-8-595-2015</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib39"><label>39</label><mixed-citation>
      
Lu, F., Li, S., Shen, B., Zhang, J., Liu, L., Shen, X., and Zhao, R.: The
emission characteristic of VOCs and the toxicity of BTEX from different
mosquito-repellent incenses, J. Hazard. Mater., 384, 121428, <a href="https://doi.org/10.1016/j.jhazmat.2019.121428" target="_blank">https://doi.org/10.1016/j.jhazmat.2019.121428</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib40"><label>40</label><mixed-citation>
      
Mahato, S., Pal, S., and Ghosh, K. G.: Effect of lockdown amid COVID-19
pandemic on air quality of the megacity Delhi, India, Sci. Total Environ.,
730, 139086, <a href="https://doi.org/10.1016/j.scitotenv.2020.139086" target="_blank">https://doi.org/10.1016/j.scitotenv.2020.139086</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib41"><label>41</label><mixed-citation>
      
McDuffie, E., Smith, S., O'Rourke, P., Tibrewal, K., Venkataraman, C., Marais, E., Zheng, B., Crippa, M., Brauer, M., and Martin, R.: CEDS_GBD-MAPS: Global Anthropogenic Emission Inventory of NO<sub><i>x</i></sub>, SO<sub>2</sub>, CO, NH<sub>3</sub>, NMVOCs, BC, and OC from 1970–2017 (2020_v1.0), Zenodo [data set], <a href="https://doi.org/10.5281/zenodo.3754964" target="_blank">https://doi.org/10.5281/zenodo.3754964</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib42"><label>42</label><mixed-citation>
      
Mor, S., Kumar, S., Singh, T., Dogra, S., Pandey, V., and Ravindra, K.:
Impact of COVID-19 lockdown on air quality in Chandigarh, India:
understanding the emission sources during controlled anthropogenic
activities, Chemosphere, 263, 127978, <a href="https://doi.org/10.1016/j.chemosphere.2020.127978" target="_blank">https://doi.org/10.1016/j.chemosphere.2020.127978</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib43"><label>43</label><mixed-citation>
      
Mozaffar, A. and Zhang, Y. L.: Atmospheric volatile organic compounds
(VOCs) in China: A review, Curr. Pollut. Rep., 31, 1–14, <a href="https://doi.org/10.1007/s40726-020-00149-1" target="_blank">https://doi.org/10.1007/s40726-020-00149-1</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib44"><label>44</label><mixed-citation>
      
Nemeth, M., Borkin, D., and Michalconok, G.: The comparison of
machine-learning methods XGBoost and LightGBM to predict energy development,
Proceedings of the Computational Methods in Systems and Software, Springer,
208–215, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib45"><label>45</label><mixed-citation>
      
Pakkattil, A., Muhsin, M., and Varma, M. R.: COVID-19 lockdown: Effects on
selected volatile organic compound (VOC) emissions over the major Indian
metro cities, Urban Clim., 37, 100838, <a href="https://doi.org/10.1016/j.uclim.2021.100838" target="_blank">https://doi.org/10.1016/j.uclim.2021.100838</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib46"><label>46</label><mixed-citation>
      
Pathakoti, M., Muppalla, A., Hazra, S., D. Venkata, M., A. Lakshmi, K., K. Sagar, V., Shekhar, R., Jella, S., M. V. Rama, S. S., and Vijayasundaram, U.: Measurement report: An assessment of the impact of a nationwide lockdown on air pollution – a remote sensing perspective over India, Atmos. Chem. Phys., 21, 9047–9064, <a href="https://doi.org/10.5194/acp-21-9047-2021" target="_blank">https://doi.org/10.5194/acp-21-9047-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib47"><label>47</label><mixed-citation>
      
Pei, C. L., Yang, W. Q., Zhang, Y. L., Song, W., Xiao, S. X., Wang, J.,
Zhang, J. P., Zhang, T., Chen, D. H., Wang, Y. J., Chen, Y. N., and Wang, X.
M.: Decrease in ambient volatile organic compounds during the COVID-19
lockdown period in the Pearl River Delta region, South China, Sci. Total
Environ., 823, 153720, <a href="https://doi.org/10.1016/j.scitotenv.2022.153720" target="_blank">https://doi.org/10.1016/j.scitotenv.2022.153720</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib48"><label>48</label><mixed-citation>
      
Qi, J., Mo, Z., Yuan, B., Huang, S., Huangfu, Y., Wang, Z., Li, X., Yang,
S., Wang, W., and Zhao, Y.: An observation approach in evaluation of ozone
production to precursor changes during the COVID-19 lockdown, Atmos.
Environ., 262, 118618, <a href="https://doi.org/10.1016/j.atmosenv.2021.118618" target="_blank">https://doi.org/10.1016/j.atmosenv.2021.118618</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib49"><label>49</label><mixed-citation>
      
Sahu, L. K., Tripathi, N., Gupta, M., Singh, V., Yadav, R., and Patel, K.:
Impact of COVID-19 pandemic lockdown in ambient concentrations of aromatic
volatile organic compounds in a metropolitan city of Western India, J.
Geophy. Res.-Atmos., 127, e2022JD036628, <a href="https://doi.org/10.1029/2022JD036628" target="_blank">https://doi.org/10.1029/2022JD036628</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib50"><label>50</label><mixed-citation>
      
Sharma, S., Zhang, M., Anshika, Gao, J., Zhang, H., and Kota, S. H.: Effect
of restricted emissions during COVID-19 on air quality in India, Sci. Total
Environ., 728, 138878, <a href="https://doi.org/10.1016/j.scitotenv.2020.138878" target="_blank">https://doi.org/10.1016/j.scitotenv.2020.138878</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib51"><label>51</label><mixed-citation>
      
Shi, X., Zhao, C., Jiang, J.H., Wang, C., Yang, X., and Yung, Y. L.: Spatial
representativeness of PM<sub>2.5</sub> concentrations obtained using observations
from network stations, J. Geophy. Res.-Atmos., 123, 3145–3158, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib52"><label>52</label><mixed-citation>
      
Shi, Z., Song, C., Liu, B., Lu, G., Xu, J., Van Vu, T., Elliott, R. J., Li,
W., Bloss, W. J., and Harrison, R. M.: Abrupt but smaller than expected
changes in surface air quality attributable to COVID-19 lockdowns, Sci.
Adv., 7, eabd6696, <a href="https://doi.org/10.1126/sciadv.abd6696" target="_blank">https://doi.org/10.1126/sciadv.abd6696</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib53"><label>53</label><mixed-citation>
      
Sun, X., Liu, M., and Sima, Z.: A novel cryptocurrency price trend
forecasting model based on LightGBM, Financ. Res. Lett., 32, 101084, <a href="https://doi.org/10.1016/j.frl.2018.12.032" target="_blank">https://doi.org/10.1016/j.frl.2018.12.032</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib54"><label>54</label><mixed-citation>
      
Tang, J., Chan, L., Chan, C., Li, Y. S., Chang, C., Liu, S., Wu, D., and Li,
Y.: Characteristics and diurnal variations of NMHCs at urban, suburban, and
rural sites in the Pearl River Delta and a remote site in South China,
Atmos. Environ., 41, 8620–8632, <a href="https://doi.org/10.1016/j.atmosenv.2007.07.029" target="_blank">https://doi.org/10.1016/j.atmosenv.2007.07.029</a>, 2007.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib55"><label>55</label><mixed-citation>
      
Van Donkelaar, A., Martin, R. V., Brauer, M., and Boys, B. L.: Use of
satellite observations for long-term exposure assessment of global
concentrations of fine particulate matter, Environ. Health Persp., 123,
135–143, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib56"><label>56</label><mixed-citation>
      
Wang, M., Lu, S., Shao, M., Zeng, L., Zheng, J., Xie, F., Lin, H., Hu, K.,
and Lu, X.: Impact of COVID-19 lockdown on ambient levels and sources of
volatile organic compounds (VOCs) in Nanjing, China, Sci. Total Environ.,
757, 143823, <a href="https://doi.org/10.1016/j.scitotenv.2020.143823" target="_blank">https://doi.org/10.1016/j.scitotenv.2020.143823</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib57"><label>57</label><mixed-citation>
      
Wei, J., Li, Z., Lyapustin, A., Sun, L., Peng, Y., Xue, W., Su, T., and
Cribb, M.: Reconstructing 1-km-resolution high-quality PM<sub>2.5</sub> data
records from 2000 to 2018 in China: spatiotemporal variations and policy
implications, Remote Sens Environ., 252, 112136, <a href="https://doi.org/10.1016/j.rse.2020.112136" target="_blank">https://doi.org/10.1016/j.rse.2020.112136</a>, 2021a.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib58"><label>58</label><mixed-citation>
      
Wei, J., Li, Z., Pinker, R. T., Wang, J., Sun, L., Xue, W., Li, R., and Cribb, M.: Himawari-8-derived diurnal variations in ground-level PM2.5 pollution across China using the fast space-time Light Gradient Boosting Machine (LightGBM), Atmos. Chem. Phys., 21, 7863–7880, <a href="https://doi.org/10.5194/acp-21-7863-2021" target="_blank">https://doi.org/10.5194/acp-21-7863-2021</a>, 2021b.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib59"><label>59</label><mixed-citation>
      
Wei, J., Li, Z., Li, K., Dickerson, R., Pinker, R., Wang, J., Liu, X., Sun,
L., Xue, W., and Cribb, M.: Full-coverage mapping and spatiotemporal
variations of ground-level ozone (O<sub>3</sub>) pollution from 2013 to 2020
across China, Remote Sens Environ., 270, 112775, <a href="https://doi.org/10.1016/j.rse.2021.112775" target="_blank">https://doi.org/10.1016/j.rse.2021.112775</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib60"><label>60</label><mixed-citation>
      
Wei, J., Li, Z., Wang, J., Li, C., Gupta, P., and Cribb, M.: Ground-level gaseous pollutants (NO<sub>2</sub>, SO<sub>2</sub>, and CO) in China: daily seamless mapping and spatiotemporal variations, Atmos. Chem. Phys., 23, 1511–1532, <a href="https://doi.org/10.5194/acp-23-1511-2023" target="_blank">https://doi.org/10.5194/acp-23-1511-2023</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib61"><label>61</label><mixed-citation>
      
Wolpert, D. H.: Stacked generalization, Neural Networks 5, 241–259, 1992.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib62"><label>62</label><mixed-citation>
      
Wu, J., Gamber, M., and Sun, W.: Does wuhan need to be in lock-down during the Chinese lunar new year?, Int. J. Env. Res. Pub. He., 17, 3, <a href="https://doi.org/10.3390/ijerph17031002" target="_blank">https://doi.org/10.3390/ijerph17031002</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib63"><label>63</label><mixed-citation>
      
Xue, T., Zheng, Y., Tong, D., Zheng, B., Li, X., Zhu, T., and Zhang, Q.:
Spatiotemporal continuous estimates of PM<sub>2.5</sub> concentrations in China,
2000–2016: A machine learning method with inputs from satellites, chemical
transport model, and ground observations, Environ. Int., 123, 345–357, 2019.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib64"><label>64</label><mixed-citation>
      
Zhai, B. and Chen, J.: Development of a stacked ensemble model for
forecasting and analyzing daily average PM<sub>2.5</sub> concentrations in
Beijing, China, Sci. Total Environ., 635, 644–658, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib65"><label>65</label><mixed-citation>
      
Zhang, M., Katiyar, A., Zhu, S., Shen, J., Xia, M., Ma, J., Kota, S. H., Wang, P., and Zhang, H.: Impact of reduced anthropogenic emissions during COVID-19 on air quality in India, Atmos. Chem. Phys., 21, 4025–4037, <a href="https://doi.org/10.5194/acp-21-4025-2021" target="_blank">https://doi.org/10.5194/acp-21-4025-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib66"><label>66</label><mixed-citation>
      
Zhao, Y. B., Zhang, K., Xu, X. T., Shen, H. Z., Zhu, X., Zhang, Y. X., Hu, Y. T., and Shen, G. H.: Substantial changes in nitrogen dioxide and ozone after excluding meteorological impacts during the Covid-19 outbreak in mainland China. Environ. Sci. Technol., 7, 402–408, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib67"><label>67</label><mixed-citation>
      
Zheng, B., Zhang, Q., Geng, G., Chen, C., Shi, Q., Cui, M., Lei, Y., and He, K.: Changes in China's anthropogenic emissions and air quality during the COVID-19 pandemic in 2020, Earth Syst. Sci. Data, 13, 2895–2907, <a href="https://doi.org/10.5194/essd-13-2895-2021" target="_blank">https://doi.org/10.5194/essd-13-2895-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib68"><label>68</label><mixed-citation>
      
Zheng, P., Chen, Z., Liu, Y., Song, H., Wu, C. H., Li, B., Kraemer, M. U.
G., Tian, H., Yan, X., Zheng, Y., Stenseth, N. C., and Jia, G.: Association
between coronavirus disease 2019 (COVID-19) and long-term exposure to air
pollution: Evidence from the first epidemic wave in China, Environ. Pollut.,
276, 116682, <a href="https://doi.org/10.1016/j.envpol.2021.116682" target="_blank">https://doi.org/10.1016/j.envpol.2021.116682</a>, 2021.

    </mixed-citation></ref-html>--></article>
