North China Plain as a hot spot of ozone pollution exacerbated by extreme high temperatures

. A large population in China has been exposed to both severe ozone (O 3 ) pollution and extreme heat under global warming. Here, the spatiotemporal characteristics of coupled extremes in surface O 3 and heat (OPCs) over China are investigated using surface observations, a process-based chemical transport model (GEOS-Chem), and multi-model simulations from Phase 6 of the Coupled Model Intercomparison Project (CMIP6). North China Plain (NCP; 37–41 ◦ N; 114–120 ◦ E) is identiﬁed as a hot spot of OPCs, where more than half of the O 3 pollution days are accompanied by high temperature extremes. OPCs over NCP exceeded 40 d during 2014–2019, exhibiting an increasing trend. Both O 3 concentrations and temperatures are elevated during OPCs compared with O 3 pollution days occurring individually (OPIs). Therefore, OPCs impose more severe health impacts to humans than OPIs, but the stronger health effects are mainly driven by the higher temperatures. GEOS-Chem simulations further reveal that enhanced chemical production resulting from hot and stable atmospheric conditions under anomalous weather patterns primarily contributes to the exacerbated O 3 levels during OPCs. In the future, CMIP6 projections suggest increased occurrences of OPCs over NCP in the middle of this century, but by the end of this century, OPCs may decrease or increase depending on the pollutant emission scenarios. However, for all future scenarios, extreme high temperatures will play an increasingly important role in modulating O 3 pollution in a warming climate.


Introduction
With rapid economic development, car ownership, and fossil fuel consumption, China has been struck by severe air pollution in recent decades (Lu et al., 2018). Research and air quality controls have been prioritized to tackle the problem of particulate matter (PM 2.5 ; T. . Since the implementation of China's Action Plan on the Prevention and Control of Air Pollution Plan in 2013, anthropogenic emissions of many air pollutants and their precursor gases, including sulfur dioxide (SO 2 ), nitrogen oxide (NO x ), carbon monoxide (CO), black carbon (BC), and organic carbon (OC), decreased by 21 %-59 % between 2013 and 2017, despite an 11 % increase in anthropogenic emissions of nonmethane volatile organic compounds (NMVOCs; Zheng et al., 2018). Correspondingly, the annual average PM 2.5 concentrations decreased from 72 to 47 µg m −3 in 74 major cities in China (Huang et al., 2018). In contrast, ozone (O 3 ) con-4706 P. Wang et al.: North China Plain as a hot spot of ozone pollution centrations in China show an apparently increasing trend during 2013-2017, with the annual average O 3 concentrations in 74 key cities increasing from 140 to 160 µg m −3 (Huang et al., 2018). During the warm season (April to September) of the same period, the daily maximum 8 h average O 3 concentration (MDA8 O 3 ) increased at a rate of 3 % yr −1 , far exceeding the rates in many other countries, such as Japan, Korea, and some in Europe (Lu et al., 2018). Long-term exposure to high O 3 concentrations can seriously damage human health, agriculture, buildings, and ecology (Sharma et al., 2017;Yue et al., 2017). Therefore, the rising O 3 concentration in recent years has caused great public concerns in China.
With global warming, extreme high temperatures and heat events have become natural hazards in China in recent decades, with substantial effect on socioeconomics, ecosystems, and human health (Lau and Nath, 2014;Meehl and Tebaldi, 2004). For instance, southern China was hit by a widespread heat wave with a record-breaking maximum temperature of 43.2 • C during summer 2003. The extreme heat event lasted for more than 40 d and caused heightened levels of human mortality (Tan et al., 2007;Wang et al., 2018). Such disastrous high temperatures have become more frequent in China (Wang et al., 2016). Mideastern China experienced an excessively long heat wave over a wide-ranging area from mid-July to mid-August 2018. The local maximum temperatures exceeded 40 • C, and the spatial extent involved 18 provinces, resulting in record-breaking overloaded power grids in many areas (M. Lu et al., 2020). Recent studies found that extreme high temperatures and heat events have intensified in the past 60 years and are expected to become more frequent and severe in the coming decades (P. Wang et al., 2019a. Extreme high temperatures are conducive to O 3 pollution. Specifically, high temperatures can increase the production rate of surface O 3 in the presence of abundant O 3 precursors (Camalier et al., 2007;Lu et al., 2019a). As O 3 concentration increases nonlinearly with temperature, extreme high temperatures have disproportionate impacts on O 3 (Lin et al., 2020). Therefore, O 3 pollution often co-occurs with extreme heat (Schnell and Prather, 2017). Besides the direct impacts of air temperatures on O 3 production, the co-occurrence of extreme heat and O 3 pollution arise from their shared underlying drivers, i.e., persistent high pressure, strong solar radiation, low humidity, and weak wind speeds (P. Wang et al., , 2018Perkins, 2015). Hence, despite reductions in anthropogenic emissions of O 3 precursors in the U.S., Europe, and China, high O 3 episodes will likely continue in the future due to increasing heat waves under climate warming .
The coupled extremes in heat and O 3 pollution lead to higher mortality rates than O 3 pollution or hot extreme acting alone (Krug et al., 2019). While the impacts of extreme high temperatures on O 3 pollution have been investigated using case studies in China (Ma et al., 2019;Pu et al., 2017), there is a gap in understanding the spatiotemporal characteristics and underlying mechanisms of coupled extremes in high temperatures and O 3 pollution due to a lack of systematic analyses. Although extreme high temperatures are expected to be more frequent and intense in the future with accelerated warming, surface O 3 concentrations are expected to decrease because of curtailment in O 3 precursor emissions. Therefore, considerable uncertainties exist in the future changes of coupled extremes in heat and O 3 .
In this study, based on the available surface O 3 concentrations and air temperature observations during 2014-2019, we investigate the spatiotemporal characteristics of cooccurrences of extremes in air temperatures and surface O 3 in China, highlighting North China Plain (NCP; defined here as 37-41 • N; 114-120 • E; see Fig. 1) as a hot spot which has already suffered from the most severe O 3 pollution in recent years (K. . The underlying mechanisms governing the coupled extreme are examined using the global chemical transport model GEOS-Chem. The associated health burden during the coupled extreme days is also discussed. In addition, future projections of the coupled extremes in the warming climate are explored based on the latest multi-model simulations from Phase 6 of the Coupled Model Intercomparison Project (CMIP6).

Observed O 3 concentration and reanalysis data
Hourly O 3 concentrations for 2014-2019 were obtained from China National Environmental Monitoring Centre (CNEMC). The network covered 944 sites in 2014 that grew to about 1600 sites in 2019. The daily maximum air temperatures (Tmax) for more than 2000 observation sites during the same period were provided by the National Meteorological Information Centre of the China Meteorological Administration (CMA). The dataset has been quality controlled, homogenized (Li et al., 2004), and widely used in previous works (Wang et al., 2019b). Here in this study, we focus on the extreme high temperatures and surface O 3 of warm season during May to September. To unify the spatial resolutions of Tmax and O 3 concentration, the two observational datasets are mapped to 1 • × 1 • grid boxes, and the values in each box represent the averaged observations within that box. The spatial distributions of averaged daily Tmax and MDA8 O 3 over May to September during 2014-2019 are shown in Fig. S1 in the Supplement. We have also tested the grid size of 0.5 • and found that the different grid resolutions have a negligible influence on the results.
Meteorological conditions during extremes of O 3 and high temperatures are calculated using variables derived from the new Japanese 55-year Reanalysis (JRA-55) at 1.25 • × 1.25 • resolution (Ebita et al., 2011), including geopotential height (HGT), winds, relative humidity (RH), 2 m air temperature (T2m), surface soil moisture (SM), downward solar radiation where x i indicates the parameter x on day i, n is the total number of days during May to September for 2014-2019, s indicates the standard deviation of the daily time series and [x i ] is the standardized anomaly for parameter x on day i. The standardized meteorological variables enable a direct comparison among their magnitudes during extreme O 3 and/or high temperatures.

GEOS-Chem model
To explore the physical and chemical mechanisms related to the O 3 extremes, the 3-D global chemical transport model (GEOS-Chem, version 12.9.3) was utilized to simulate O 3 concentrations during May-September for 2014-2017, driven by assimilated meteorological data of version 2 of Modern-Era Retrospective analysis for Research and Applications (MERRA-2; Gelaro et al., 2017). The simulations were performed at a horizontal resolution of 2 • latitude × 2.5 • longitude with 47 vertical levels. By examining the simulations of surface O 3 over the U.S. with a regional climate model and the global GEOS-Chem model, Fiore et al. (2003) indicated that the ability to resolve local O 3 maxima is compromised, but the spatial correlation improves when the model resolution coarsens. The coarse-resolution global model can successfully capture the synoptic-scale pro-cesses modulating O 3 concentrations, whereas a finer spatial resolution may improve the representation of processes occurring on smaller scales. The anthropogenic emissions of O 3 precursor gases, including CO, NO x , and volatile organic compounds (VOCs) in China, were obtained from the MEIC emission inventory (http://meicmodel.org/, last access: 10 April 2020), which included emissions from industry, power, residential, and transportation sectors. Biogenic volatile organic compound (BVOC) emissions also play a vital role in modulating the formation of ozone and secondary organic aerosols (Ma et al., 2021;Gao et al., 2022).

CMIP6 data
We use O 3 and Tmax outputs from future projections of Scenario Model Intercomparison Project (ScenarioMIP) in the CMIP6 archive to determine how the coupled extremes will change in a warmer climate. ScenarioMIP is the primary activity within CMIP6 that provides multi-model climate projections driven by different scenarios of future emissions and land use changes (O'Neill et al., 2016), produced based on the shared socioeconomic pathways (SSPs) combining socioeconomic developments and the feedback of global climate changes (Z. Li et al., 2020). More details about the SSP scenarios can be found in O'Neill et al. (2016). Currently, four SSP scenarios in ScenarioMIP simulations provide hourly O 3 concentration and daily Tmax from the present day to the end of the 21st century (2015-2100), i.e., SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5 (combination of low, intermediate, relatively high, and high societal vulnerabilities and forcing levels, respectively). Among the four SSPs, SSP3-7.0 and SSP2-4.5 have the weakest and medium air pollution controls pathways, respectively, while strong air pollution controls are assumed in SSP1-2.6 and SSP5-8.5 (Gidden et al., 2019). Five global climate models (GCMs), i.e., MOHC.UKESM1-0-LL, CESM2-WACCM, GFDL-ESM4, MPI-ESM-1-2-HAM, and EC-Earth3-AerChem from ScenarioMIP under CMIP6 that provide both hourly O 3 and daily Tmax, were adopted in this work. The horizontal resolutions and institutions of the five GCMs are listed in Table S1 in the Supplement. Note that the numbers of available models vary across different scenarios (see Table S2 in the supplement for details). The results from the five GCMs are re-gridded to the observation boxes using linear interpolation to facilitate spatial comparison. In this study, 2015-2019 is regarded as the historical period and the overall performance of the CMIP6 simulations in reproducing the occurrences of coupled extremes is evalu-ated against the observations during 2015-2019. For the projection of coupled extremes, we focused on the two periods of 2046-2050 and 2096-2100 under different SSPs.

Identification of extremes in O 3 and temperature
Following Schnell and Prather (2017), in this study, we use the local-specific thresholds for each grid to identity the extreme cases of surface air temperatures and O 3 concentrations, specifically, the 90th percentile of daily Tmax and daily MDA8 O 3 from May to September for 2014-2019. The local-specific thresholds have been widely used in recent studies of ozone pollution (e.g., Schnell and Prather, 2017;Lin et al., 2019;Qin et al., 2021). Note that the 90th percentile of MDA8 O 3 over NCP, Yangtze River Delta, Sichuan Basin, and Pearl River Delta are 97.7, 84.4, 73.7, and 76.8 ppb, respectively, close to China's grade II air quality standard for MDA8 O 3 (around 80 ppb under standard atmospheric conditions).
To characterize the co-occurrences of extremes in high temperatures and surface O 3 and investigate the impacts of extreme high temperatures on O 3 pollution, the following extremes are defined: -Total O 3 pollution days (OPs): all days when daily MDA8 O 3 is above its threshold.
-Individual O 3 pollution days (OPIs): days when MDA8 O 3 is above its threshold while Tmax is lower than its threshold.
-Coupled extreme days (OPCs): days when both daily Tmax and daily MDA8 O 3 exceed their corresponding thresholds.
We use a co-occurrence frequency ratio (CF, in percent) to characterize the dependence of extreme high O 3 levels on extreme high temperatures. CF is defined as the ratio of the frequency of OPCs (days) to the frequency of OPs (days). Thus, a higher CF value indicates a higher dependence of O 3 pollution on extreme high temperatures:

Health impact of coupled extremes
Following Lee et al. (2017), in this study, we apply a ratio index to describe the combined human health impacts caused by O 3 and temperatures during OPCs, which represents the potential enhancement in mortality rates (MR) related to O 3 and temperature levels during OPC compared to those during OPIs. And the MR is defined as below: MR = Daily mortality during OPCs Daily mortality during OPIs where RR ozone,i (RR ozone,j ) and RR temperature,i (RR temperature,j ) are the relative risks due to O 3 concentration and temperature exceedance, respectively, on a coupled extreme day i (an individual O 3 pollution day j ), m is the total days of coupled extremes, and n is the total days of individual O 3 pollution day. MR ozone (MR temperature ) is the enhanced mortality rates attributed to O 3 concentration (temperature) changes, while MR is the combined effects from both O 3 and temperature changes.
Because China has higher air pollution levels and may also differ in terms of age structure, population sensitivity to air pollution/heat exposures, and components of air pollution mixture compared with developed countries (Chen et al., 2018), we use China-specific concentration and temperature response functions in the present study, as indicated in the recent nationwide studies (Yin et al., 2017;Huang et al., 2015). β 1 is the concentration response factor corresponding to a 0.24 % [95 % confidence interval: 0.13 %, 0.35 %] increase in daily mortality per 10 µg m −3 increase in MDA8 O 3 above C 0 (Yin et al., 2017). Following Huang et al. (2015) in calculating RR temperature in 66 Chinese communities, β 2 indicates a 1.09 % [95 % confidence interval: 0.72 % to 1.46 %] excess mortality per 1 • C increase in temperature above T 0 . Note that the algorithms here to calculate MR, MR ozone , and MR temperature do not consider the possible amplification/inhibition effect of combining O 3 and air temperature in affecting human health. Previous studies have claimed that O 3 -related mortality increases with higher temperatures, although several studies have presented contrasting results or inconsistent relationships for different regions (Chen et al., 2014;Jhun et al., 2014;Ren et al., 2008). By analyzing the total mortality rates associated with short-term O 3 exposure over East Asia among four seasons, Chen et al. (2014) found that the higher temperatures in summer significantly increased the O 3 -related mortality rates.

Spatial and temporal patterns of coupled extremes
The spatial patterns of OPCs and their ratio to the total O 3 pollutions days (CF values) during May to September for the recent 6 years (2014-2019) highly resemble each other (Fig. 1), with the highest values located over NCP which has suffered the most severe O 3 pollution in recent years (Fig. S1a in the Supplement). The highest OPCs ex-ceed 40 d over NCP and the corresponding CF is more than 56 % (Fig. 1). That means, the coupled extreme days account for more than half of the total O 3 pollution days, indicating a strong dependence of O 3 pollution on extreme high temperatures over NCP. It has been suggested that the dependence of O 3 concentration on high temperature increases with the O 3 levels (Lin et al., 2020). However, coupled extremes occur much less frequently over the Yangtze River Delta (YRD; 30-33 • N; 118-122 • E) compared with NCP, and the regional averaged CF in YRD is below 20 %, even though the MDA8 O 3 level and temperature in YRD are both as high as those in NCP (Fig. S1 in the Supplement). The distinctive relationships between extreme high temperature and O 3 concentration over NCP and YRD are driven by their different climatology during warm season. Southern China receives substantial monsoon rainfall during summer, accompanied by increased relative humidity and reduced radiation (Zhou and Yu, 2005), which can suppress surface O 3 levels (Han et al., 2020). Delineating the local daily maximum air temperatures (Tmax) and RH of all O 3 pollution days over NCP and YRD (Fig. S2 in the Supplement), OPCs occur more frequently over NCP than over YRD, and a higher fraction of the O 3 pollution days over NCP co-occur with extreme high temperatures and low-to-moderate RH (Fig. S2a in the Supplement). Humid environment dampens the occurrence of O 3 pollution over YRD and extreme O 3 pollution mostly occurs on days with relatively low RH when air temperatures are moderate (Fig. S2b), which explains the lower OPCs and CF in YRD compared with NCP. Therefore, we focus on the coupled extremes over NCP.
Daily variations of the occurrence of OPIs and OPCs over NCP during 2014-2019 are shown in Fig. 2. O 3 pollution days have appeared since 2015, but coupled extremes OPCs have only been observed since 2017, mostly during May to July (Fig. 2a). The abrupt increase in the occurrence of coupled extremes in 2017 is consistent with the significant increasing trends of both MDA8 O 3 and Tmax (95 % confidence level) over NCP in recent years (Fig. 2b). The strong increasing trends of MDA8 O 3 and air temperatures are consistent with previous results (K. Li et al., , 2020. As addressed previously (K. Li et al., 2020), the temperature trends during 2014-2019 reflect interannual climate variability rather than a long-term warming trend. Notably, daily MDA8 O 3 exhibits increasing sensitivity to Tmax from 2014 to 2019 (Fig. 2c), supporting the increase in OPCs during the same time period. Note that the linear regression slopes between daily MDA8 O 3 and Tmax are not strictly monotonically increasing. For example, the slopes are 3.96, 3.43, and 4.56 ppb • C in 2017, 2018, and 2019 (Fig. 2c) Both O 3 levels and air temperatures are higher during OPCs than during OPIs over the NCP region (Fig. 3a, b), with the regional mean anomalies of Tmax and MDA8 O 3 during OPCs reaching 3.36 • and 5.49 ppb, respectively, compared with those during OPIs. A north-south contrast in the MDA8 O 3 and Tmax difference between OPCs and OPIs is evident (Fig. 3b), suggesting that contrasting environments north and south of the Yangtze River during the summer monsoon may play a key role in the dependence of O 3 pollution on extreme Tmax in China. Figure 4 shows the composites of normalized anomalies (see Sect. 2) of meteorological fields during coupled extreme days over NCP for 2014-2019. During OPCs, anomalous high pressure and anticyclonic circulation dominate NCP and the surrounding region north of the Yangtze River in the midtroposphere (500 hPa), with anomalous easterlies prevailing over NCP (Fig. 4a). Associated with the anomalous highpressure system is clear sky with enhanced downward solar radiation (DSR) at the surface (Fig. 4c), leading to hotter near surface temperature (Fig. 4b), reduced RH and soil moisture (Fig. 4d, e), and enhanced surface sensible heat flux (Fig. 4f) that further intensifies the temperatures (Fig. 4b).

Weather patterns and ozone processes during coupled extremes
These anomalous conditions are all stronger during OPCs than OPIs over NCP (Fig. S3 in the Supplement) and more conducive to O 3 pollution . Among the meteorological factors, the intensification in surface temperatures is the strongest among different meteorological variables with the highest magnitudes ( Fig. S3b in the Supplement), supporting that air temperature is the most influential meteorological variable of surface O 3 over NCP (K. . The impacts of weather patterns on surface O 3 level can be understood via changes in physical and chemical processes, both sensitive to meteorology . The contributions of different chemical and physical processes to OPCs over NCP under the anomalous weather pattern of Fig. 4 are quantified by GEOS-Chem simulations of O 3 during May to September of 2014-2017. GEOS-Chem can reasonably capture the spatial pattern and magnitude of OPCs in observations during 2014-2017 (Text S1 and Fig. S4 in the Supplement). Four processes affecting O 3 levels are considered, including net chemical production, horizontal advection, vertical advection, and mixing (diffusion plus dry deposition), and are listed in Table 1. For both OPIs and OPCs, chemical production contributes the most to O 3 mass within the boundary layer. Compared with OPIs, the higher O 3 levels during OPCs (Fig. 3b) are contributed by stronger chemical production and weakened mixing, but vertical advection and horizontal advection tend to reduce the O 3 concentrations, with enhanced chemical production playing the dominant role. Therefore, we conclude that the hotter near surface temperature induced by anomalous weather patterns and amplified by land-atmosphere feedbacks during OPCs (Fig. 4) is the primary cause of the enhanced formation of O 3 and eventually a higher surface O 3 level than during OPIs. Besides meteorological effects, the O 3 precursor emissions should partially contribute to the spatiotemporal variations of OPCs over China. It is reported that surface O 3 pollution levels are strongly correlated with daytime surface temperatures, especially in highly polluted regions, with strong precursor emissions (Porter and Heald, 2019). NCP has the highest anthropogenic emissions compared with the other regions in China, which should benefit the higher correlations between surface O 3 and air temperatures, and thus the higher OPCs therein. Moreover, the increasing trend of OPCs over NCP in recent years may be associated with the continued anthropogenic increases in O 3 , as well as the unmitigated emissions of VOCs (K. Li et al., 2020), emphasizing the need for controlling anthropogenic emissions of VOCs. In addition, Fu et al. (2015) have indicated that the enhanced biogenic emissions and the accelerated photochemical reaction rates both act to increase surface ozone over the USA during 1988-2011. Thus, the increasing trend of biogenic emissions due to vegetation biomass variability over China (Cao et al., 2022) may also have potential impacts on the variations of OPCs.

Health impacts of coupled extremes
As both surface O 3 and air temperatures are amplified during coupled relative to individual O 3 pollution days (Fig. 3), we investigated the potential influences of OPCs on human health. The enhanced mortality rates for OPCs compared with OPIs during May to September for each year of 2017-2019 are illustrated in Fig. 5 and attributed to air temperature and/or O 3 concentration changes (MR temperature , MR ozone , and MR; see Sect. 2). It should be noted that coupled extreme days are only observed since 2017. MR, MR ozone , and MR Temperature are above 1.0 for all 3 years, indicating a harsher environment for people to survive during OPCs. Importantly, MR temperature is significantly higher than MR ozone for all years of 2017-2019, suggesting that extreme high temperature caused many more mortalities than extreme O 3 concentrations over NCP. The averaged MR ozone , MR Temperature , and MR for 2017-2019 are 1.003, 1.037, and 1.040, respectively. Compared with the individual O 3 pollution days OPIs, the daily mortality rate in NCP increases by 4.0 % during coupled extremes OPCs, the majority of which are attributed to the temperature increase, with less than one-tenth contributed by the O 3 concentration increase. That is, coupled extremes amplify health impacts compared with individual O 3 pollution days primarily because of the higher mortality risk associated with elevated air temperatures. Moreover, we estimate that around 100 daily excess deaths over NCP are attributable to the higher temperatures and O 3 level during OPCs than OPIs (See Text S3 in the Supplement).

Projected coupled extremes in future climate
As O 3 precursors (i.e., NO x and NMVOCs) are expected to keep declining due to the continued emission controls in China while extreme high temperatures will become more frequent and intense under global warming, uncertainties exist in the projection of the co-occurrences of extremes in high temperatures and O 3 pollution. Here, we investigate the projections of OPCs and CF values based on CMIP6 simulations under SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5. OPCs in the simulations are identified in the same way as for the observations (see Text S2 and Fig. S5 in the Supplement for details). We focus on the historical period of 2015-2019 (referred to as 2019) and the projected periods   S5 in the Supplement). The averaged OPCs over NCP under each SSP increase from the historical period to the mid-century (Fig. 6a), with a maximum increase under SSP5-8.5 (spatial distribution shown in Fig. S6 in the Supplement). From mid-century to the end of the century, OPCs decrease under SSP1-2.6, SSP2-4.5, and SSP5-8.5, but OPCs by 2100 obviously sur-   (Turnock et al., 2020), MDA8 O 3 in 2100 under this scenario is highest among the four SSPs (Fig. 7). In contrast, OPCs are substantially reduced to below 5 d by 2100 under SSP1-2.6 and SSP2-4.5, highlighting the benefit of strong actions in mitigating climate and reducing air pollutant emissions. In the future by 2050 and 2100, NCP will still be the most vulnerable region in China to the coupled extreme (Fig. S6 in the Supplement), while most other areas will be much less threatened by the coupled extremes by the end of the century under d,and h in the Supplement).
Unlike OPCs, CF over NCP obviously increases by the 2050 and 2100 scenarios compared with 2019 under all four SSPs (Fig. 6b). The projected increases of CF over NCP indicate the higher dependence of O 3 pollution on extreme high temperatures in the future, consistent with the increased sensitivity of MDA8 O 3 to Tmax at higher Tmax in historical periods (Fig. 2c). Spatially, the NCP region will still see the highest CF values in the future, especially under SSP1-2.6, . This means that regardless of the economic pathways, extreme high temperature will play an increasingly important role in modulating O 3 pollution in the warming climate. Therefore, besides the management strategies on pollutant emissions, global warming mitigations will undoubtedly benefit O 3 pollution control, especially for regions facing severe air quality issues. Note that for the future changes of OPCs, the influences of natural variability are less considered, whereas previous studies have emphasized the significant role of natural variability on altering the robustness of climate projections and their impacts on air quality (e.g, Garcia-Menendez et al., 2017). The detection of the anthropogenic-forced signal demands a larger model ensemble and a longer simulation length that deserves further explorations.

Discussion and conclusions
Climate change can impact local air quality. Higher temperatures associated with climate change can lead to an increase in surface O 3 , and high temperatures and surface O 3 are highly temporally correlated over many regions (Porter and Heald, 2019). A large population in China has been increasingly exposed to both severe O 3 pollution and extreme heat under global warming. With combined surface observations of air temperature and O 3 concentration, processbased model simulations, and multi-model projections, this study first presented a comprehensive analysis of the cooccurrences of extreme high temperatures and O 3 pollution in China. It is highlighted that NCP is a hot spot in China most threatened by the co-occurrence of extremes in heat and O 3 pollution. The higher co-occurrence over NCP than other regions in China is linked to their distinctive relationships to meteorological variables, as temperature is the top meteorological factor directly leading to O 3 pollution over NCP, whereas relative humidity is the most influential variable for O 3 pollution over southern China (Han et al., 2020). Recently, the compound extreme events (e.g., co-occurrence of two extreme weather events simultaneously) have been raised as a substantial concern to O 3 formation. For example, the co-occurrences of heat wave and air stagnation promote higher O 3 concentration compared with the single extreme events of heat wave or stagnation in the USA in the future relative to the present Gao et al., 2020).
The concurrent increasing trends in both surface O 3 and temperature over NCP in recent years account for the increasing coupled extremes in surface O 3 and heat in recent years. Besides, has been previously reported that the increasing trend of temperature is higher over northern China than southern China (P. Qian and Qin, 2006). The increase in air temperature can accelerate the O 3 production. Using a physically based model (GEOS-Chem), we have provided support for the dominant role of higher temperatures associated with stable atmospheric conditions under favorable weather patterns in amplifying O 3 pollution through enhanced chemical production during coupled extremes, compared with the individual ozone pollution days not accompanied by extreme temperatures. In addition, the increases in surface O 3 over NCP are much stronger than the other regions in recent years, which is also possibly linked to the stimulation effect from enhanced hydroperoxyl radicals (HO 2 ) due to a reduction in aerosol sink resulting from the decrease in PM 2.5 during this period (K. . Thus, the hot spot of co-occurrences of extremes in heat and O 3 over NCP could be attributed to the co-effects of stronger increasing trends of temperature and surface O 3 therein. It is a prevalent concept that the coupled extremes pose greater health impacts or risks to humans than the simply summed impacts of the single extremes acting alone (Smith et al., 2014). It is revealed here that both the O 3 concentration and air temperatures are higher during the coupled ex- tremes than the individual O 3 pollution, leading to an even heavier health burden to humans. This study also underscores the elevated air temperatures during the coupled extremes as the major driver for increased mortality rates, while the simultaneously elevated O 3 concentrations act as an additional stressor. However, as mentioned above, how the interactions between temperature and O 3 influence human health during coupled extremes is still an unanswered question that deserves future studies using more health-related data.
Currently, China has the highest emission of greenhouse gases, and the emission rates have increased significantly since the beginning of the 21st century (Friedlingstein et al., 2020). To prevent the dangerous climate change impacts, the Chinese government has declared an ambitious goal by pledging to peak emissions before 2030 and reach carbon neutrality before 2060. With global warming, hot extremes in China are projected to be more frequent, stronger, and longer lasting, which may present challenges for O 3 pollution control. Based on ScenarioMIP simulations from CMIP6, this study demonstrates that the coupled extremes over NCP are projected to be more frequent in the middle of this century, but their frequency decreases or increases by the end of the century under strong or weak air pollution control scenarios, respectively. Also, with higher sensitivity of O 3 concentration to higher temperatures, O 3 extremes will increasingly co-occur with extreme high temperatures over NCP as the climate warms, regardless of the economic pathways. Thus, our results further reinforce the notion that determined actions are vital to make our communities less vulnerable to climate change impacts already in progress. On the other hand, tropospheric O 3 levels are projected to increase in the near decades (see Fig. 7b; Turnock et al., 2020). As the third important anthropogenic greenhouse gas after CO 2 and CH 4 , higher tropospheric O 3 levels can cause temperature changes by altering the energy balance between the atmosphere and the earth (Dang and Liao, 2019), which may cause feedback on the air quality. Thus, potential co-benefits may be gained through O 3 pollution control and climate change managements in suppressing the occurrences of coupled extremes and tackling their consequences to air quality, human health, and climate.