Global assessment of climatic responses to ozone–vegetation interactions

. The coupling between surface ozone (O 3 ) and vegetation signiﬁcantly inﬂuences the regional to global climate. O 3 uptake by plant stomata inhibits the photosynthetic rate and stomatal conductance, impacting evapotranspiration through land surface ecosystems. Using a climate–vegetation–chemistry coupled model (the NASA GISS ModelE2 coupled with the Yale Interactive terrestrial Biosphere, or ModelE2-YIBs), we assess the global climatic responses to O 3 –vegetation interactions during the boreal summer of the present day (2005– 2014). High O 3 pollution reduces stomatal conductance, resulting in warmer and drier conditions worldwide. The most signiﬁcant responses are found in the eastern US and eastern China, where the surface air temperature increases by + 0.33 ± 0.87 and + 0.56 ± 0.38 °C, respectively. These temperature increases are accompanied by decreased latent heat and increased sensible heat in both regions. The O 3 –vegetation interaction also affects atmospheric pollutants. The surface maximum daily 8 h average O 3 concentrations increase by + 1.46 ± 3.02 ppbv in eastern China and + 1.15 ± 1.77 ppbv in the eastern US due to the O 3 -induced inhibition of stomatal uptake. With reduced atmospheric stability following a warmer climate, increased cloud cover but decreased relative humidity jointly reduce aerosol optical depth by − 0.06 ± 0.01 ( − 14.67 ± 12.15 %) over eastern China. This study suggests that vegetation feedback should be considered for a more accurate assessment of climatic perturbations caused by tropospheric O 3 .

Experimental studies have shown that excessive O 3 exposure reduces both plant photosynthesis and stomatal conduc-tance (Ainsworth et al., 2012;Lombardozzi et al., 2013).The reduction rates are dependent on the O 3 stomatal fluxes and the damage sensitivities, which vary among different vegetation types (Nussbaum and Fuhrer, 2000;Karlsson et al., 2004;Pleijel et al., 2004).Several exposure-based indices, such as accumulated hourly O 3 concentrations over a threshold of 40 ppb (AOT40) and the sum of all hourly average concentrations (SUM00), are used to assess O 3 -induced vegetation damage (Fuhrer et al., 1997;Paoletti et al., 2007).In addition, the flux-related POD y method (phytotoxic O 3 dose above a threshold flux of y) is also widely applied to consider the dynamic adjustment of stomatal conductance (Buker et al., 2015;Sicard et al., 2016).Considering the variability of plant sensitivities, different O 3 damage schemes have been proposed to quantify the impacts of O 3 on land Published by Copernicus Publications on behalf of the European Geosciences Union.9924 X. Zhou et al.: Global assessment of climatic responses to ozone-vegetation interactions carbon assimilation at regional to global scales (Anav et al., 2011;Lam et al., 2023;Lei et al., 2020).For example, Sitch et al. (2007) calculated simultaneous damage to both photosynthesis and stomatal conductance on the basis of instantaneous O 3 stomatal uptake.In contrast, Lombardozzi et al. (2012) estimated decoupled reductions in plant photosynthesis and stomatal conductance via different response relationships to cumulative O 3 stomatal uptake.The application of different schemes has resulted in a wide range of reductions in gross primary productivity (GPP) of 2 %-12 % globally, with regional hotspots of up to 20 %-30 % (Lombardozzi et al., 2015;Unger et al., 2020;Zhou et al., 2024).
O 3 -induced inhibition of stomatal conductance decreases dry deposition and consequently enhances surface O 3 concentrations (Clifton et al., 2020;Wesely and Hicks, 2000;Zhang et al., 2006).Using the Sitch et al. ( 2007) scheme with high O 3 damage sensitivities in ModelE2-YIBs (NASA GISS ModelE2 coupled with the Yale Interactive terrestrial Biosphere model), Gong et al. (2020) revealed that O 3 -vegetation interactions increased regional O 3 concentrations by 1.8 ppbv in the eastern US, 1.3 ppbv in Europe, and 2.1 ppbv in eastern China in 2010.In comparison, Sadiq et al. (2017) reported consistently stronger feedbacks on O 3 concentrations in these polluted regions via the scheme of Lombardozzi et al. (2012) embedded in the Community Earth System Model (CESM).Moreover, the inclusion of online O 3 -vegetation interactions in numerical models will also result in a greater loss of simulated land carbon assimilation due to the feedbacks of both ecosystems and surface O 3 .This is attributable to several factors.On the one hand, O 3 damage to leaf photosynthesis inhibits plant growth and decreases the leaf area index (LAI), leading to a greater reduction in GPP than in simulations without LAI changes (Yue et al., 2020).On the other hand, O 3 enhancement due to vegetation feedback may cause additional vegetation damage and result in further GPP losses (Lei et al., 2021).As a result, O 3 -vegetation interactions should be considered in global estimates of O 3 damage to ecosystem functions.
In addition to affecting surface O 3 , O 3 -vegetation interactions can also alter water and energy exchanges between the land and atmosphere through the modulation of stomatal conductance.For example, Lombardozzi et al. (2015) used the Community Land Model (CLM) and estimated that the cumulative uptake of O 3 by leaves resulted in a reduction of 2.2 % in transpiration but an increase of 5.4 % in runoff globally.Arnold et al. (2018) used CESM and reported that plant exposure to O 3 could decrease land-air moisture fluxes and atmospheric humidity, which would further reduce shortwave cloud forcing in polluted regions and induce widespread surface warming up to +1.5 K. Two recent studies utilized the WRF-Chem model and revealed considerable warming and associated meteorological perturbations due to O 3 -vegetation interactions in China (Zhu et al., 2022;Jin et al., 2023).However, all these modeling studies applied the same O 3 vegetation damage scheme proposed by Lombar-dozzi et al. (2012).It is necessary to assess the climatic responses to O 3 -vegetation interactions via different schemes to explore robust responses and associated uncertainties.
In this study, we quantified the global impacts of O 3vegetation interactions on climatic conditions and surface air pollutants during the 2010s using ModelE2-YIBs (Yue and Unger, 2015).This fully coupled framework was implemented with the semi-mechanistic O 3 damage scheme proposed by Sitch et al. (2007), which calculates aggregated O 3 damage to photosynthesis on the basis of varied sensitivities to instantaneous stomatal O 3 uptake across eight plant functional types (PFTs).We performed sensitivity experiments to quantify the responses of surface air temperature and precipitation to O 3 -vegetation interactions.The feedbacks to aerosols and O 3 concentrations were also examined.

Model descriptions
ModelE2-YIBs is a fully coupled climate-carbon-chemistry model that combines the NASA GISS ModelE2 with the YIBs vegetation model.ModelE2 is a general circulation model with a horizontal resolution of 2°× 2.5°in latitude and longitude and 40 vertical layers up to 0.1 hPa.It dynamically simulates gas-phase chemistry (NO x , HO x , O x , CO, CH 4 , and non-methane volatile organic compounds (NMVOCs)), aerosols (sulfate, nitrate, black and organic carbon, dust, and sea salt), and their interactions (Menon and Rotstayn, 2006).Both the physical and chemical processes are calculated every 0.5 h, and the radiation module is called every 2.5 h.The radiation module includes direct and indirect aerosol radiative effects and accounts for the absorption of multiple greenhouse gases (GHGs).For cloud optical parameters, Mie scattering, ray tracing, and matrix theory are used (Schmidt et al., 2006).The model outperforms 20 other IPCC-class climate models in simulating surface solar radiation (Wild et al., 2013) and has been extensively validated for meteorological and hydrological variables against observations and reanalysis data (Schmidt et al., 2014).
The YIBs model employs the well-established Farquhar model for leaf photosynthesis and the Ball-Berry model for stomatal conductance (Farquhar et al., 1980;Ball et al., 1987) as follows: A tot = min(J c , J e , J s ). (1) Here, the total leaf photosynthesis, denoted as A tot (µmol m −2 [leaf] s −1 ), is calculated considering both C 3 (Collatz et al., 1991) and C 4 plants (Collatz et al., 1992).A tot is derived from the minimum value of the constraints.The rate of carboxylation limited by ribulose-1,5bisphosphate carboxylase (Rubisco) is J c : The carboxylation rate restricted by the availability of light is J e : (3) The export-limited rates for C 3 plants and the phosphoenolpyruvate carboxylase (PEPC)-limited rates of carboxylation for C 4 plants are represented by J s : In these functions, V cmax (µmol m −2 s −1 ) is the maximum carboxylation capacity.c i and O i (Pa) represent the internal leaf CO 2 and oxygen partial pressure, respectively.* (Pa) denotes the CO 2 compensation point, whereas K c and K o (Pa) are Michaelis-Menten constants for the carboxylation and oxygenation of Rubisco, respectively.The parameters * , K c , and K o vary with temperature on the basis of the sensitivity of the vegetation to temperature (Q 10 coefficient).PAR (µmol m −2 s −1 ) is the absorbed photosynthetically active radiation, a leaf is the leaf-specific light absorbance that considers sunlit and shaded leaves, and α is the quantum efficiency.P atm (Pa) represents the ambient pressure.K s is set to 4000 as a constant following Oleson et al. (2010) to limit photosynthesis to C 4 plants that become saturated at lower CO 2 concentrations.
The stomatal conductance (g s , mol [H 2 O] m −2 s −1 ) is linked to the variations in A tot , with several parameters, such as the dark respiration rate (R d , µmol m −2 s −1 ), relative humidity (RH), and CO 2 concentration at the leaf surface (c s ): The model simulates the biophysical processes of eight PFTs, including tundra, C 3 and/or C 4 grass, shrubland, deciduous broadleaf forest, evergreen broadleaf forest, evergreen needleleaf forest, and cropland.Different values are assigned to parameters m and b for each PFT (Table S1 in the Supplement).The carbon taken up by the leaf then accumulates and is allocated to different organs to support plant development, resulting in dynamic changes in the LAI and tree growth.

O 3 -vegetation damage scheme
The YIBs model employs a semi-mechanistic parameterization proposed by Sitch et al. (2007) to estimate the impact of O 3 on photosynthesis through stomatal uptake.The scheme applies an undamaged factor (F ) (nmol m −2 s −1 ) to calculate both A tot and g s as follows: where A totd and g sd are the unaffected photosynthesis and stomatal conductance, respectively.The factor F is defined as follows: where a h (mmol m −2 s −1 ) is the high O 3 sensitivity coefficient, calibrated by Sitch et al. ( 2007) on data from field observations by Karlsson et al. (2004) and Pleijel et al. (2004) to represent the "high" sensitivity of the relative species of each PFT.
) is the specific threshold for O 3 damage.Both a h and F O 3 ,crit vary with vegetation types (Table S1): where [O 3 ] represents surface O 3 concentrations, and R a (s m −1 ) represents aerodynamic resistance, which expresses the turbulent transport efficiency in transferring sensible heat and water vapor between the land surface and a reference height.The constant k O 3 = 1.67 is the ratio of stomatal resistance to O 3 , which is estimated on the basis of the theoretical stomatal resistance to water (Laisk et al., 1989).When plants are exposed to [O 3 ] (Eq.9), A tot and g s decrease (Eqs.6 and 7) if the excess O 3 enters the leaves (Eq.8).The increased stomatal resistance protects plants by reducing the O 3 uptake of stomata.Consequently, the damage scheme describes changes in both the photosynthetic rate and stomatal conductance.

Experiments
To explore the coupled O 3 -vegetation effect, we performed two simulations using the ModelE2-YIBs model.adopted from the land use harmonization (LUH2) dataset (Hurtt et al., 2020).For each time slice simulation, the model was run for 30 years with all the input data fixed, and the first 10 years were used as the spin-up period.We calculated the average of the last 20 years and focused on the boreal summer season (June-July-August, JJA), when the interaction of vegetation and surface O 3 reaches its annual maximum (Fig. S3).To show the uncertainty introduced by the internal variability of the model, all the related global and regional values are denoted as mean/sum ± standard deviation of the last 20 model years.We explored the climatic responses to O 3 -vegetation interactions as the differences between O3_online and O3_offline at the global scale, with a focus on hotspot regions, such as the eastern US (30-40°N, 80-90°W) and eastern China (22.5-38°N, 106-122°E).

Data for model evaluation
We evaluated the simulated air pollutants, carbon fluxes, and meteorological variables from the O3_offline simulation against observational and reanalysis datasets.(Kalnay et al., 1996).For precipitation, we used the monthly data averaged from 2005 to 2014 from the Global Precipitation Climatology Project (GPCP) (Huffman et al., 1997;Adler et al., 2018).All these datasets were interpolated to the same resolution as the ModelE2-YIBs model.The root-mean-square error (RMSE) and normalized mean bias (NMB) were applied to quantify the deviations of the simulations from the observations: Here, S i and O i represent the simulated and observed values, respectively, and n denotes the total number of grid points used in the comparisons.

Control simulation and model evaluations
We first evaluated the air pollutants simulated by the control simulation (O3_offline) of the ModelE2-YIBs model (Fig. 1).Over a total of 503 grids with site-level O 3 measurements (Fig. 1b), the model replicated both the magnitude and spatial distribution of MDA8 O 3 , with a correlation coefficient (r) of 0.59 and an NMB of −2.54 % (Fig. 1c).The simulated summertime surface MDA8 O 3 concentration was high in regions with high anthropogenic emissions, such as western Europe and eastern China (Ohara et al., 2007), as well as in central Africa, which has frequent fire emissions (van der Werf et al., 2017).At the global scale, the model yielded an average MDA8 O 3 concentration of 43.93 ppbv, and observations revealed an average of 44.72 ppbv over the same grids.However, the model overestimated the concentrations over the North China Plain and slightly underestimated them over the US, likely due to biases in the emission inventories and the predicted climate that drive O 3 production.The AOD simulated at 550 nm by O3_offline (Fig. 1d) showed a spatial pattern similar to that of the satellite retrievals (Fig. 1e), with r = 0.75 and an NMB of −7.35 % globally (Fig. 1f).Both simulations and observations revealed AOD hotspots over northern Africa and the Middle East, where dust emissions dominate, and in northern India and eastern China, where anthropogenic emissions are high (Feng et al., 2020).
We then evaluated the simulated GPP and LAI via a control experiment for the boreal summer period (Fig. 2).The observations revealed GPP hotspots over boreal forests, such as those in the eastern US, Eurasia, and East Asia, and over tropical forests, such as those in the Amazon, central Africa, and Indonesia (Fig. 2b).The seasonal total GPP was estimated to be 41.63 Pg [C], which accounted for 35 % of the annual amount.The simulations captured the observed GPP pattern at the global scale, with r = 0.64 and NMB = −7.81% across 2581 grids (Fig. 2c), with underestimations in the tundra area and slight overestimations in the tropical-rainforest and evergreen-forest regions.The model simulated a seasonal total GPP of 38.69 Pg [C], equivalent to 34 % of the annual amount.The simulated LAI showed  patterns similar to those of GPP (Fig. 2d) and resembled the observed LAI (Fig. 2e), with a spatial correlation of r = 0.79 and a low NMB = −5.43% across 4435 grids globally (Fig. 2f).

O 3 damage to terrestrial ecosystems
We assessed the damaging effects of surface O 3 on ecosystems due to online O 3 -vegetation interactions (Fig. 3).The impacts of O 3 on biospheric variables were mainly located in regions characterized by abundant vegetation cover and elevated O 3 concentrations.On the global scale, O 3 induced a GPP reduction of −1.80 ± 0.61 Pg C yr −1 (−4.69 ± 1.56 %, Fig. 3a).This deleterious effect was more pronounced in specific regions, notably eastern China and the eastern US, with significant GPP declines of −25.40 ± 1.90 % and −20.14 ± 5.02 %, respectively, under conditions of high O 3 sensitivity (Fig. 3a and Table S2).Moreover, stomatal conductance significantly decreased in the middle latitudes of the Northern Hemisphere (Fig. 3b).The most substantial relative change of −30.62 ± 4.30 % was observed in eastern China, followed by −25.65 ± 9.32 % in the eastern US (Fig. 3b and Table S2).Although there are positive responses in some regions, they are not dominant and are hardly significant.These values were stronger than those for GPP (Fig. 3a), likely because of the climatic feedback to O 3 -vegetation interactions.
The opening of the plant stoma plays a crucial role in regulating energy and water exchange between the land surface and the atmosphere.The inhibition of stomatal conductance by surface O 3 leads to a warmer (Fig. 4a) and drier (Fig. 4b) climate in those hotspot regions, resulting in even stronger inhibitory effects on stomatal conductance.Following the changes in GPP, the global LAI decreased by, on average, 0.01 ± 0.01 m 2 m −2 (−0.62 ± 0.84 %), with regional maxima of −4.53 ± 1.14 % in eastern China and −5.87 ± 3.11 % in the eastern US (Table S2).

Global climatic responses to O 3 -vegetation interactions
In response to the O 3 -induced inhibition of stomatal conductance, surface air temperature increased by 0.05 ± 0.20 °C (Fig. 4a), whereas precipitation decreased by −0.01 ± 0.03 mm d −1 (Fig. 4b) at the global scale.The most significant change was the warming of 0.56 ± 0.38 °C and the precipitation reduction of −0.79 ± 1.05 mm d −1 (−16.18 ± 20.38 %) in eastern China (Table S3), followed by the greatest inhibition of stomatal conductance (Fig. 3b).Such warming and rainfall deficits also appeared in the eastern US and western Europe, where the O 3 -vegetation interactions were notable.The O 3 -induced inhibition of stomatal conductance decreased the latent heat flux (Fig. 4e) and consequent precipitation (Fig. 4b) in those hotspot regions.Moreover, the reduced latent heat flux promoted higher surface air temperatures (Fig. 4a), resulting in an increase in the sensible heat flux (Fig. 4f).Such warming has also been reported in field experiments, where relatively high O 3 exposure resulted in noticeable increases in canopy temperature, along with reductions in transpiration (Bernacchi et al., 2011;VanLoocke et al., 2012).Globally, temperature and precipitation showed patchy responses with both positive and negative anomalies, suggesting that the regional hotspots of O 3induced meteorological changes propagate to surrounding areas through atmospheric perturbations.
We further examined the changes in air humidity and cloud cover.The surface relative humidity decreased by −0.18 ± 0.53 % globally, with a similar pattern to that of precipitation (Fig. 4c).The most significant reductions occurred over eastern China and the eastern US, where both warming (Fig. 4a) and rainfall deficits (Fig. 4b) contributed to drought.However, in adjacent regions, such as northern China and the central US, both rainfall and surface relative humidity increased.These changes were associated with the regional increase in cloud cover (Fig. 4d).The sensible heat flux increased by 6.3 ± 5.4 W m −2 (16.54 ± 15.59 %) and 7.12 ± 3.86 W m −2 (25.46 ± 14.71 %) in the eastern US and eastern China, respectively, suggesting the transfer of thermal energy from the land to the atmosphere via O 3vegetation interactions (Fig. 4f and Table S3).The warming effect further triggered anomalous updrafts in the lower troposphere, represented by changes in vertical velocity (Fig. 5), leading to enhanced convection; reduced atmospheric stability; and, consequently, an increase in low-level cloud cover (Fig. 4d).However, despite the usual cooling effect associated with increased cloud cover due to reductions in radiation, in regions predominantly influenced by O 3 -vegetation interactions, this cooling effect was outweighed by O 3induced warming through the inhibition of stomatal conductance.Therefore, temperatures exhibited an overall increase of 0.56 ± 0.38 °C in eastern China and 0.33 ± 0.87 °C in the eastern US (Table S3).

Changes in air pollution caused by O 3 -vegetation interactions
Changes in surface water and heat fluxes induced by O 3vegetation interactions could feed back to affect air pollutants, such as O 3 and aerosols.As shown in Fig. 6a and Table S4, surface MDA8 O 3 concentrations increased by 1.46 ± 3.02 ppbv in eastern China and 1.15 ± 1.77 ppbv in the eastern US due to decreased dry deposition following O 3 inhibition of stomatal conductance.This indicates that high contemporary O 3 pollution may worsen air quality through O 3 -vegetation interactions.However, negative O 3 changes were predicted in the central US and western China, where increased rainfall dampened O 3 through chemical reactions and wet deposition.On a global scale, sur-   are smaller than those induced by O 3 -vegetation interactions (Table S2).
Aerosols also exhibited evident changes in response to O 3 -vegetation interactions.The AOD significantly decreased over hotspot regions, such as eastern China and the eastern US (Fig. 6b).In the ModelE2-YIBs model, sulfate was especially sensitive to clouds, which could enhance aerosol scavenging through cloud water precipitation (Koch et al., 2006).The large increase in cloud cover removed sulfate more efficiently than the other aerosol species did, leading to an average decrease of −1.94 ± 1.67 µg m −3 (−8.52 ± 6.88 %) in the PM 2.5 loading over eastern China (Fig. S4 and Table S4).Moreover, the reduction in surface relative humidity (Fig. 4c) in the regions with strong O 3 -vegetation interactions limited the hygroscopic growth of aerosols, leading to a more noticeable decrease in AOD (Petters and Kreidenweis, 2007;Pitchford et al., 2007) of −0.06 ± 0.05 (−14.67 ± 16.75 %) in eastern China (Table S4).Similar aerosol changes were found in the eastern US but with smaller reductions in PM 2.5 and in AOD of −0.27 ± 0.36 µg m −3 (−6.01 ± 7.9 %) and −0.01 ± 0.01 (−8.15 ± 9.38 %), respectively (Table S4).In addition to the key O 3 -vegetation coupling regions, positive but insignificant changes in AOD were predicted, leading to moderate AOD changes at the global scale (Fig. 6b).

Discussion and conclusions
We examined O 3 -vegetation feedbacks to climate and air pollution in the 2010s via the fully coupled climate-carbonchemistry model ModelE2-YIBs.During boreal summer, surface O 3 resulted in strong damage to GPP and inhibited stomatal conductance, with regional hotspots over eastern China and the eastern US.Consequently, surface transpiration was weakened, leading to decreased latent heat fluxes and relative humidity but increased surface air temperature.Moreover, surface warming increased cloud cover by reducing atmospheric stability.However, the increase in cloud cover decreased the surface temperature and promoted precipitation outside the key regions with intense O 3 -vegetation interactions.The O 3 -induced inhibition of stomatal conductance resulted in a localized increase in O 3 concentration.In contrast, the increased cloud cover and decreased relative humidity jointly reduced the AOD in hotspot regions.At the global scale, the mean changes in both climate and air pollution were moderate because of the offset between the changes with opposite signs.
Our predictions of the changes in water and heat fluxes caused by O 3 -vegetation interactions were consistent with those of previous studies (Lombardozzi et al., 2015;Arnold et al., 2018;Gong et al., 2020).For example, simulations by Lombardozzi et al. (2015) revealed that surface O 3 reduces global GPP by 8 %-12 % and reduces transpiration by 2 %-2.4 %, with regional reductions of up to 20 % in GPP and 15 % in transpiration in eastern China and the US.These changes are generally consistent with our results, although we predicted greater reductions in transpiration than in GPP due to O 3 -vegetation interactions.Using the same scheme as Lombardozzi et al. (2015), Sadiq et al. (2017) reported that O 3 -vegetation coupling induced surface warming of 0.5-1 °C and O 3 enhancement of 4-6 ppbv in eastern China and the eastern US.The magnitude of these responses was much stronger than our predictions, likely because they considered the accumulation effect of O 3 .In contrast, regional simulations by Jin et al. (2023) revealed that O 3 -vegetation coupling led to increases in temperature and in surface O 3 of up to 0.16 °C and 0.6 ppbv, respectively, in eastern China, both of which were smaller than our predictions.The damage scheme they use, which depends on cumulative O 3 uptake, omits the difference in impact on sunlit or shaded leaves https://doi.org/10.5194/acp-24-9923-2024Atmos.Chem.Phys., 24, 9923-9937, 2024  (Cao et al., 2024).The discrepancies in O 3 -vegetation feedbacks when the same O 3 damage schemes were used revealed uncertainties in the climate and chemistry models.Our predictions were within the range of previous estimates for both climatic and O 3 changes.There were several limitations in our simulated O3vegetation interactions.First, the semi-mechanistic O 3 damage scheme we used in the present study linked damage to photosynthesis with damage to stomatal conductance (Sitch et al., 2007), leading to a greater percentage of inhibition of stomatal conductance than of photosynthesis considering O 3 -vegetation feedbacks.However, some observations have shown that damage to stomatal conductance occurs more slowly and might not be proportional to the decline in photosynthetic rates (Gregg et al., 2006;Lombardozzi et al., 2012).Second, observations have shown large variability in terms of plant sensitivity to O 3 damage.The Sitch et al. ( 2007) scheme employs low to high ranges of sensitivity to indicate interspecific variabilities.In this study, we employed only high O 3 sensitivity to explore the maximum responses.The possible uncertainties due to varied O 3 damage sensitivities deserve further investigation.Third, large-scale observations were not available to validate the simulated regional to global responses of climate and air pollutants.The O 3 vegetation damage scheme has been extensively validated against sitelevel measurements of both photosynthesis (Yue and Unger, 2018) and stomatal conductance (Yue et al., 2016).However, we were conservative about the derived global responses given that previous studies have shown large discrepancies when the same O 3 damage scheme is used while being implemented in different climate and/or chemistry models (Lombardozzi et al., 2015;Sadiq et al., 2017;Jin et al., 2023).Furthermore, the 2°× 2.5°resolution of the current version of the ModelE2-YIBs model has limitations due to the high computational demands.However, high-resolution models exhibit improved simulations of extreme events (Chang et al., 2020;Ban et al., 2021), which have certain effects on O 3vegetation interactions (Mills et al., 2016;Lin et al., 2020).While chemical transport models with relatively coarse resolutions can increase biases in simulated air pollutants, they still capture large-scale patterns similar to those of fineresolution results and compare reasonably well to observational data (Wang et al., 2013;Li et al., 2016;Lei et al., 2020).Moreover, we omit the slow climatic feedback caused by air-sea interactions in the simulations.Studies have revealed that these interactions may result in different climatic perturbations compared to those in simulations with fast responses by the land surface alone (Yue et al., 2011).A dynamic ocean model would enrich future research.Moreover, this study did not isolate the different impacts of aerosols, even though the radiation module included both direct and indirect radiative effects.We will investigate this further in the future by identifying the main processes involved.Despite these uncertainties, our simulations revealed considerable changes in both climate and air pollutants in response to O 3 -vegetation interactions.The most intense warming, dryness, and O 3 enhancement were predicted in eastern China and the eastern US, affecting the regional climate and threatening public health for these top two economic centers.In contrast, for the first time, we revealed a reduction in aerosol loading in hotspot regions, suggesting both positive and negative effects on air pollutants via O 3vegetation feedbacks.Such interactions should be considered in Earth system models to better project future changes in climate and air pollutants following anthropogenic interventions against both O 3 precursor emissions and ecosystem functions.
Data availability.The observational data and model outputs that support the findings in this study are available from the corresponding authors upon reasonable request.
Author contributions.XY conceived the project.XZ performed the model simulations, conducted the analysis, and wrote the draft of the paper.XY, CT, and XL assisted in the interpretation of the results and contributed to the discussion and improvement of the paper.
Competing interests.The contact author has declared that none of the authors has any competing interests.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper.While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Figure 1 .
Figure 1.Evaluation of the present-day boreal summertime (June-August) air pollutants simulated by the ModelE2-YIBs model.The MDA8 O 3 (a-c) and AOD (d-f) from the O3_offline simulation (a, d) and observations (b, e) are compared.The correlation coefficient (r), rootmean-square error (RMSE), normalized mean bias (NMB), and number of grid cells (n) for the comparisons are listed on the mean bias maps (c, f).

Figure 2 .
Figure 2. Same as in Fig. 1 but for gross primary productivity (GPP; a-c) and the leaf area index (LAI; d-f).

Figure 3 .
Figure 3. Changes in present-day boreal-summertime biospheric variables induced by O 3 -vegetation interactions.The results shown are the changes in (a) GPP, (b) canopy conductance, and (c) LAI between the O3_online and O3_offline simulations.The black dots denote areas with significant changes (p < 0.1).

Figure 4 .
Figure 4. Changes in present-day boreal-summertime meteorological fields caused by O 3 -vegetation interactions.The results shown are changes in (a) surface air temperature, (b) precipitation, (c) surface relative humidity, (d) low-level cloudiness, (e) latent heat flux, and (f) sensible heat flux between the O3_online and O3_offline simulations.For heat fluxes, positive values (shaded in red) indicate that the upward fluxes change.The black dots denote areas with significant changes (p < 0.1).

Figure 5 .
Figure 5. Vertical profile of the vertical velocity.The results shown are changes in the vertical velocity in (a) eastern China and (b) the eastern US between the O3_online and O3_offline simulations.The solid red line denotes 0. Please note the differences in the scales.

Figure 6 .
Figure 6.Changes in present-day summertime atmospheric pollution caused by O 3 -vegetation interactions.The results shown are the changes in (a) O 3 , (b) AOD, and (c) PM 2.5 between the O3_online and O3_offline simulations.The black dots denote areas with significant changes (p < 0.1).

Financial support .
This research has been supported by the National Key Research and Development Program of China (grant no.2023YFF0805404) and the National Natural Science Foundation of China (grant no.42293323).