Articles | Volume 23, issue 2
Atmos. Chem. Phys., 23, 1043–1071, 2023
Atmos. Chem. Phys., 23, 1043–1071, 2023
Research article
20 Jan 2023
Research article | 20 Jan 2023

Biogenic isoprene emissions, dry deposition velocity, and surface ozone concentration during summer droughts, heatwaves, and normal conditions in southwestern Europe

Biogenic isoprene emissions, dry deposition velocity, and surface ozone concentration during summer droughts, heatwaves, and normal conditions in southwestern Europe
Antoine Guion1,a, Solène Turquety1,b, Arineh Cholakian2, Jan Polcher2, Antoine Ehret1, and Juliette Lathière3 Antoine Guion et al.
  • 1LMD/IPSL, Sorbonne Université, ENS, PSL Université, École Polytechnique, Institut Polytechnique de Paris, CNRS, 75005 Paris, France
  • 2LMD/IPSL, École Polytechnique, Institut Polytechnique de Paris, ENS, PSL Université, Sorbonne Université, CNRS, 91128 Palaiseau, France
  • 3LSCE/IPSL, CEA-CNRS-UVSQ, Université Paris-Saclay, 91190 Gif-sur-Yvette, France
  • anow at: INERIS, 60550 Verneuil-en-Halatte, France​​​​​​​
  • bnow at: LATMOS/IPSL, Sorbonne Université, UVSQ, CNRS, 75005 Paris, France

Correspondence: Antoine Guion (


At high concentrations, tropospheric ozone (O3) deteriorates air quality, inducing adverse effects on human and ecosystem health. Meteorological conditions are key to understanding the variability in O3 concentration, especially during extreme weather events. In addition to modifying photochemistry and atmospheric transport, droughts and heatwaves affect the state of vegetation and thus the biosphere–troposphere interactions that control atmospheric chemistry, namely biogenic emissions of precursors and gas dry deposition. A major source of uncertainty and inaccuracy in the simulation of surface O3 during droughts and heatwaves is the poor representation of such interactions. This publication aims at quantifying the isolated and combined impacts of both extremes on biogenic isoprene (C5H8) emissions, O3 dry deposition, and surface O3 in southwestern Europe.

First, the sensitivity of biogenic C5H8 emissions, O3 dry deposition, and surface O3 to two specific effects of droughts, the decrease in soil moisture and in biomass, is analysed for the extremely dry summer 2012 using the biogenic emission model MEGANv2.1 and the chemistry transport model CHIMEREv2020r1. Despite a significant decrease in biogenic C5H8 emissions and O3 dry deposition velocity, characterized by a large spatial variability, the combined effect on surface O3 concentration remains limited (between +0.5 % and +3 % over the continent).

The variations in simulated biogenic C5H8 emissions, O3 dry deposition, and surface O3 during the heatwaves and agricultural droughts are then analysed for summer 2012 (warm and dry), 2013 (warm), and 2014 (relatively wet and cool). We compare the results with large observational data sets, namely O3 concentrations from Air Quality (AQ) e-Reporting (2000–2016) and total columns of formaldehyde (HCHO, which is used as a proxy for biogenic emissions of volatile organic compounds) from the Ozone Monitoring Instrument (OMI) of the Aura satellite (2005–2016).

Based on a cluster approach using the percentile limit anomalies indicator, we find that C5H8 emissions increase by +33 % during heatwaves compared to normal conditions, do not vary significantly during all droughts (either accompanied or not by a heatwave), and decrease by 16 % during isolated droughts. OMI data confirm an average increase in HCHO during heatwaves (between +15 % and +31 % depending on the product used) and decrease in HCHO (between 2 % and 6 %) during isolated droughts over the 2005–2016 summers. Simulated O3 dry deposition velocity decreases by 25 % during heatwaves and 35 % during all droughts. Simulated O3 concentrations increase by +7 % during heatwaves and by +3 % during all droughts. Compared to observations, CHIMERE tends to underestimate the daily maximum O3. However, similar sensitivity to droughts and heatwaves are obtained. The analysis of the AQ e-Reporting data set shows an average increase of +14 % during heatwaves and +7 % during all droughts over the 2000–2016 summers (for an average daily concentration value of 69 µg m−3 under normal conditions). This suggests that identifying the presence of combined heatwaves is fundamental to the study of droughts on surface–atmosphere interactions and O3 concentration.

1 Introduction

Tropospheric ozone (O3) plays a critical role in maintaining the oxidative capacity of the troposphere. However, as a high oxidant, it also deteriorates air quality at high concentrations, inducing adverse effects on human and ecosystem health. Both short- and long-term O3 exposure significantly increase the risk of morbidity and mortality from cardiovascular or respiratory causes (e.g. Jerrett et al.2009; Nuvolone et al.2018). It also causes visible damage on plant epidermis and photosynthesis inhibition, thereby seriously threatening vegetation growth (Anand et al.2021) and crop yields (De Andrés et al.2012). Tai et al. (2014) estimate that O3 pollution enhanced by global warming (Representative Concentration Pathway 8.5) could lead to a global crop yield reduction of 3.6 % in 2050.

Tropospheric O3 is a secondary air pollutant formed by photochemical reaction chains initiated by the oxidation of volatile organic compounds (VOCs) in the presence of nitrogen oxides (NOx). The VOC / NOx emission ratio determines the chemical regime of O3 production at the local scale in a highly nonlinear relationship (Jacob1999). Methane (CH4) and carbon monoxide (CO) are also O3 precursors in background conditions.

Meteorological conditions are a key driver for O3 formation and transport (Mertens et al.2020). Indeed, O3 production is favoured by warm and sunny conditions, and O3 peaks therefore occur mainly during summer. The overall objective of this paper is to quantify the variation in surface O3 concentrations over the Euro–Mediterranean area during extreme weather events of heatwaves and droughts. The Euro–Mediterranean region is often affected by both heatwaves (Russo et al.2015) and droughts (Spinoni et al.2015). The intensity and frequency of such events have increased over the last few decades in the region. According to projections based on climate models, these trends could last and worsen over the 21st century (Perkins-Kirkpatrick and Gibson2017; Spinoni et al.2018).

Heatwaves are often associated with high O3 concentration. For example, Vautard et al. (2005) show that the persistent heatwave of August 2003 led to high O3 concentration over western Europe for several days (daily maximum above 150 µg m−3; spatially averaged). Jaén et al. (2021) report hourly values of O3 exceeding 250 µg m−3 during a shorter heatwave (28–29 June 2019) at a local scale (Barcelona). These values are well above the European Union (EU) recommendation of an 8 h average maximum of 120 µg m−3.

Porter and Heald (2019) highlight that, despite a strong correlation between surface temperature and O3 concentration identified in Europe, this relationship is nonlinear and complex, as it depends on several mechanisms including precursor emissions, reaction rates, and lifetimes in the atmosphere. High temperatures favour emissions of biogenic VOCs (BVOCs), which can induce an increasing formation of O3 in the case of a NOx-limited regime, which is found in many areas in the Mediterranean Basin (Richards et al.2013). Persistent and intense solar radiation also allows long photochemical episodes (Jaén et al.2021). Since heatwaves in the Mediterranean are often related to a blocking situation of atmospheric systems, stagnant anticyclonic conditions lead to the accumulation of both precursors and O3 (Otero et al.2022). Moreover, extreme temperature and high O3 concentration can cause stomatal closure, thus reducing the dry deposition of O3, as highlighted by Gong et al. (2020) with the O3–vegetation feedback.

Droughts also significantly impact atmospheric chemistry through biosphere and atmospheric cycle perturbation (Wang et al.2017). Drought impact on air pollutants and, more specifically, O3 is difficult to investigate for several reasons. It is, firstly, related to their characterization. Due to their multiscalar character, droughts can last on timescales ranging from days to months, making their extent difficult to assess (Vicente-Serrano et al.2010). Moreover, drought definition depends on the type considered (Svoboda and Fuchs2016). Meteorological droughts correspond to a rainfall deficit or an excess of evapotranspiration, agricultural ones to soil water shortage for plant growth (soil dryness), and hydrological ones to surface and/or underground flow decrease. Another difficulty is that droughts affect not only the atmosphere but also the land biosphere through soil dryness and plant activity decline (e.g. Vicente-Serrano2007; Vicente-Serrano et al.2013). The variability in O3 concentration is generally not well represented in chemistry transport models (CTMs) compared to observations, partially because of the lack of interactions between the meteorology, terrestrial biosphere, and atmospheric chemistry (Wang et al.2017). Uncertainties related to BVOC emission modelling could lead to a BVOC-derived O3 uncertainty estimated at about 50 % in Europe (Curci et al.2010). Another major cause of uncertainties concerns the simulation of meteorological conditions.

Meteorological conditions specific to droughts are low relative humidity and low or absent precipitation and cloud cover. The latter leads to an increase in downward solar radiation. Temperature can also be higher during droughts than during normal conditions. Wang et al. (2017) found that drought periods in the United States over 1990–2014 were characterized by temperatures up to 2 C hotter and radiation 12.4 W m−2 higher. Their analysis of surface measurements shows a mean enhancement in O3 concentration of 3.5 ppbv (parts per billion by volume; 8 %), which is explained by enhanced precursor emissions and photochemistry. Finally, large number of O3 precursors emitted during biomass burning enhanced by droughts and heatwaves can contribute to O3 pollution peaks (e.g. Hodnebrog et al.2012).

The development of droughts and heatwaves can be linked (Miralles et al.2019). For example, through the soil moisture–temperature feedback, droughts can increase heatwave intensity due to a decrease in evapotranspiration and an increase in sensible heat (e.g. Zampieri et al.2009). It is therefore important to integrate such interactions for the simulation of droughts and heatwaves.

Through soil moisture deficit, droughts can induce considerable adverse effects on the biosphere, leading not only to a decrease in the overall biomass but also of the BVOC emission activity. Guion et al. (2021) report an averaged decrease of  10 % in leaf area index (LAI) observed by satellite during droughts in southern Europe. Areas of low altitude and vegetation with short root systems are more sensitive. Based on simulations from the Model of Emissions of Gases and Aerosols from Nature (MEGAN; Guenther et al.2006, 2012) and using MODIS (Moderate Resolution Imaging Spectroradiometer) LAI between 2003 and 2018, Cao et al. (2022) find that vegetation biomass variability in China is a major controlling factor of the interannual variation in BVOC emissions. Furthermore, Emmerson et al. (2019) simulated (also with MEGAN) a reduction in C5H8 emissions of 24 %–52 % over southeastern Australia due to soil dryness in summer.

Dry deposition velocity in the Mediterranean depends on both stomatal and non-stomatal conductance (Lin et al.2020; Sun et al.2022). Stomatal conductance that can be considerably affected under hot and dry conditions depends on several parameters. As identified by Kavassalis and Murphy (2017), the leaf-to-air vapour pressure deficit that depends on both temperature and relative humidity is a key variable regulating the stomatal conductance. Extreme hot and dry conditions can cause stomatal closure to reduce transpiration. Moreover, the lack of rainfall prevents wet deposition, and plant water stress slows dry deposition. However, there are still large uncertainties about the relationship between the soil moisture and the gas dry deposition (e.g. Clifton et al.2020b).

Considering the difficulty for models in simulating soil moisture accurately (Cheng et al.2017) and the uncertainties about plant and BVOC species response to water stress (Bonn et al.2019), the impact of soil moisture is often dismissed in biogenic emissions and dry deposition schemes used in CTMs. Therefore, knowledge about the respective effect on O3 remains limited. To evaluate the effects of droughts and heatwaves on surface O3 in southwestern Europe, the variation in both BVOC emissions and dry deposition as a function of meteorological and hydrological conditions are analysed in this paper including the soil dryness effect.

Agricultural droughts and heatwaves are identified based on the coupled meteorological and land surface vegetation regional model RegIPSL (regional Earth system model of Institute Pierre Simon Laplace) for the 1979–2016 period (Guion et al.2021). Summers are less affected by drought events (33 % of them) than by heatwaves (97 %). Nevertheless, drought periods are generally much longer than heatwaves (22.3 d vs. 5.8 d). Those two types of extreme weather events can be considerably correlated over the western Mediterranean (temporal correlation R between 0.4 and 0.5).

Atmospheric chemistry is simulated using the regional CHIMERE CTM (Menut et al.2021) over three selected summers (2012, 2013, and 2014), including the online calculation of biogenic emissions from the MEGANv2.1 model. Simulations are analysed in conjunction with surface O3 observations from the European surface network Air Quality e-Reporting (, last access: 30 July 2021). Concentrations of formaldehyde (HCHO) are used as a proxy for BVOC emissions. HCHO is produced in high yield during the oxidation of hydrocarbons. As an intermediate product of the oxidation of most VOC and being characterized by a short lifetime (about a few hours), HCHO has been used in many studies to infer VOC emissions (e.g. Millet et al.2006, 2008; Curci et al.2010). Observations of HCHO total column by the Ozone Monitoring Instrument (OMI) aboard the Aura satellite (Levelt et al.2006) are used.

The observations and models used for this research are presented in Sects. 2 and 3, respectively. The RegIPSL, CHIMERE, and MEGAN models, as well as the different modelling experiments undertaken, are detailed. Validation of surface O3, NO2, and 2 m temperature are presented in Sect. 4. The results (see Sect. 5) are divided into three parts. First, the sensitivity of C5H8 emissions, O3 dry deposition, and surface O3 concentration to drought effects (biomass decrease and soil dryness) is assessed in Sect. 5.1. Second, the variation in observed and simulated O3, integrating drought and heatwave effects, is analysed with respect to the variation in troposphere–atmosphere interactions during extreme weather events (see Sect. 5.2). Finally, the number and distribution of days exceeding the EU air quality standard are presented (see Sect. 5.3).

2 Observations

2.1 European surface network of O3 and NO2 concentration

The in situ measurements of surface O3 and NO2 provided by the European Environment Agency are used (EEA; The Air Quality (AQ) e-Reporting project gathers air quality data provided by EEA member countries as well as European collaborating countries. Merging with the statistics of the AirBase v8 project (2000–2012), AQ e-Reporting offers a multiannual time series of air quality measurement data until 2021, categorized by station types and zones across Europe. Only measurements from stations classified as background and rural are considered for this study.

Thunis et al. (2013) quantify the various sources of uncertainty for O3 measurements (e.g. linear calibration and ultraviolet photometry) and estimated a total uncertainty of 15 %, regardless of concentration level. There are also considerable uncertainties in the measurements of NO2. Lamsal et al. (2008) emphasize that the chemiluminescence analyser, the measurement technique primarily found in air quality stations, is subject to significant interference from other reactive species containing oxidized nitrogen (e.g. PAN and HNO3). This can lead to an overestimation of measured NO2 concentrations.

2.2 Satellite observations of HCHO total column

The HCHO total column retrieved from the observations of the backscatter solar radiation in the ultraviolet (306 and 360 nm) by the OMI/Aura instrument are used as a quantitative proxy for BVOC emissions. OMI has an Equator-crossing time at 13:45 LT (local time)​​​​​​​ on the ascending node, about 14.5 revolutions per day, and a swath width of 2600 km, allowing daily global coverage (Levelt et al.2006, 2018).

Uncertainties in the HCHO retrieval are significant, varying between 30 % (pixels with high concentrations) and 100 % (with low concentrations), mainly due to cloud and aerosol scattering along the field of view (González Abad et al.2015). Satellite retrievals of HCHO have been shown to have systematic low mean bias (20 %–51 %) compared to aircraft observations (Zhu et al.2016).

For this reason, we chose to focus on the relative difference in total HCHO during droughts and heatwaves (no correction factor has been applied). In order to quantify the uncertainty in the HCHO anomalies obtained, the analysis has been performed using two products. We primarily use the OMHCHOd level 3 product (Chance2019), which provides the HCHO total column with a spatial resolution of 0.1× 0.1. For comparison, the level 2 retrieval by the Belgian Institute for Space Aeronomy (BIRA) is also used (De Smedt et al.2015, hereafter referred to as OMI-BIRA). Both products are included in the intercomparison conducted by Zhu et al. (2016).

Only observations with a cloud fraction less than 0.3, a solar zenith angle less than 70, and a vertical column density within the range of 0.8 × 1015 and 7.6 × 1015 molec. cm−2 are selected in order to minimize OMI row anomalies, following Zhu et al. (2017). For a suitable comparison with model simulations, HCHO data were regridded on the Med-CORDEX (Coordinated Regional Climate Downscaling Experiment over the Mediterranean) domain using a bilinear method.

2.3 In situ observations of 2 m temperature, C5H8 surface concentration, and O3 dry deposition flux

To complete the validation of our simulations, several observational data sets are added to our study. First, in situ observations of temperature at 2 m above the surface (T2 m) from the E-OBS data set (Cornes et al.2018) are used to validate the simulated temperature by the Weather Research and Forecasting (WRF) model (Skamarock et al.2008). The E-OBS gridded product, with a resolution of 0.25× 0.25, was regridded to the Med-CORDEX domain (see Sect. 3.1.1).

Second, flux measurements are used for comparison with the different modelling experiments carried out to analyse the sensitivity of biogenic C5H8 emissions and O3 deposition to the effects of biomass decrease and soil dryness. However, there are few flux measurements available that cover at least several weeks during summer 2012. The Ersa station (FR0033R), located at Cape Corsica in France (42.97 N, 9.38 E), is used to assess the surface concentration of C5H8 measured by a steel canister instrument at 4.0 m above the surface. The data are provided by the EBAS infrastructure (, last access: 14 April 2022).

Dry deposition flux measurements of O3 are also used at one station. This is the Castel Porziano station (IT-Cp2), located in the Lazio region of Italy (41.70 N, 12.36 E), measuring at 14.9 m above the surface (eddy covariance technique). A full description of the measurement data is provided in Fares et al. (2013). The data were downloaded from the European Fluxes Database Cluster (EFDC;, last access: 16 October 2022).

3 Models


3.1.1 Model description

The regional WRF-ORCHIDEE model (RegIPSL) couples the land surface model of ORCHIDEE (ORganising Carbon and Hydrology In Dynamic EcosystEms; Maignan et al.2011) and the meteorological model WRF. ORCHIDEE is composed of three modules, namely SECHIBA, which simulates the water and energy cycle, STOMATE, which resolves the processes of the carbon cycle, allowing therefore an interactive phenology, and, finally, LPJ, which computes the competition between the plant functional types (PFTs). However, the LPJ module is not used for this research.

The RegIPSL model has been used following the recommendation of the Med-CORDEX initiative (Ruti et al.2016). On a Lambert conformal projection, the domain covers the Euro–Mediterranean region with a spatial resolution of 20 km. The full description of the configuration used with RegIPSL is presented in Guion et al. (2021).

The performed Med-CORDEX simulations over the 1979–2016 period have been evaluated against surface observations of temperature and precipitation, demonstrating that it is suitable for research on droughts and heatwaves (Guion et al.2021).

3.1.2 Indicators of drought and heatwave

Based on simulations performed by the RegIPSL model, we identified heatwaves and agricultural droughts over the Euro–Mediterranean region using indicators of extreme weather events (Guion2022; data access:, last access: 14 April 2022).

Following the approach of Lhotka and Kyselý (2015), we computed the percentile limit anomalies (PLAs) of temperature 2 m above the surface (PLAT2 m), for heatwave detection, and of soil dryness (PLASD), for agricultural drought detection. The soil dryness is computed as the complement of the soil wetness index, which is described in ORCHIDEE as the ratio between the soil moisture and the accessible water content.

For each cell of longitude (i) and latitude (j) of the domain, the monthly distribution of the 2 m temperature/soil dryness is normalized, and the 75th percentile is computed (p=0.75). The daily PLAT2m/SD indicator (Eq. 1) is equal to the daily deviation (dXi,j,d) between the reference variable (Xi,j,d) and the corresponding monthly percentile (Xi,j,mp). A heatwave/drought is identified when the daily PLAT2m/SD is positive for 3 consecutive days.

(1) d X i , j , d = X i , j , d - X i , j , m p .

The chosen percentile for the detection of extreme weather events is generally higher, between the 80th and 95th percentile (e.g. Stéfanon et al.2012). We chose the 75th percentile in this study for a larger population. As a result, events detected here cover not only extreme events but also periods that are significantly drier and warmer than normal conditions.

PLAT2 m was calculated based on 2 m temperature observations (E-OBS data set) in Guion et al. (2021). Although the intensity of the heatwaves was slightly overestimated with the Med-CORDEX simulations (+0.16 C for the mean bias and +0.50 C for the maximum bias), their temporal correlation with observations was high (R coefficient of about 0.9 over the whole Mediterranean).

3.1.3 Soil water stress for plants

The soil water stress (SWS) calculated by RegIPSL is used to describe the water stress for plant transpiration and photosynthesis (de Rosnay et al.2002). It is an index varying between 0 (plants fully stressed) and 1 (not stressed) that takes into account the water transfer between the soil layers along root profiles (11 layers in the version of ORCHIDEE used here). For each PFT (p), soil class (s), and soil layer (v), it is computed as follows:

(2) SWS p , s = SM p , s , v - SMw p , s , v SMnostress p , s , v - SMw p , s , v × nroot p , s , v ,

where SM is the soil moisture (liquid phase), SMw is the level at wilting point, SMnostress is the level at which SWS reaches 1, and nroot is the normalized root length fraction in each soil layer (between 0 and 1). There are three soil classes to avoid competition for water uptake, namely bare soil, soil with short root systems, and soil with long ones.


3.2.1 Model description and configuration

The CHIMERE v2020r1 Eulerian 3-dimensional regional CTM (Menut et al.2021) computes gas-phase chemistry and aerosol formation, transport, and deposition. It can be guided by precalculated meteorology (offline simulation) or used online with the WRF regional meteorological model (Skamarock et al.2008), including the option with the direct and indirect effects of aerosols. In this study, we chose to use the online version of CHIMERE but without any aerosol effects in order to reduce the calculation time and to isolate possible feedbacks.

Here, the reduced scheme MELCHIOR2 (Derognat2003) is used for gas-phase chemistry (44 species; almost 120 reactions). The aerosol module includes aerosol microphysics, secondary aerosol formation mechanisms, aerosol thermodynamics, and deposition, as detailed in Couvidat et al. (2018) and Menut et al. (2021). The photolysis rates are calculated online using the Fast-JX module version 7.0b (Bian and Prather2002), which accounts for the radiative impact of aerosols.

The WRF model used with CHIMERE, which has a different configuration from the RegIPSL model (Guion et al.2021), has 15 vertical layers, from 998 up to 300 hPa. The physics interface for the surface is provided by the Noah land surface model. The aerosol-aware Thompson scheme (Thompson and Eidhammer2014) is used for the microphysics parameterization. Horizontal and vertical advection are based on the scheme of Van Leer (1977). The reanalysis meteorological data for the initial and boundary conditions are provided by the National Centers of Environmental Prediction (NCEP). The physical and chemical time steps are, respectively, 30 and 10 min.

The simulations have been performed for June, July, and August (JJA, with a 5 d spinup) of the years 2012, 2013, and 2014 over two nested simulation domains (Fig. 1). The large domain covers northern Africa and Europe at 60 km horizontal resolution (EUMED60; 164 × 120 cells), and the small domain is close to that of Med-CORDEX but smaller; it covers the western Mediterranean region at 20 km horizontal resolution (Med-CORDEX; 222 × 93 cells). The area studied in more detail in the following is southwestern Europe (35–46 N, 10 W–20 E).

Figure 1Domains used for the CHIMERE simulations. EUMED60 at 60 km resolution is the large domain, and the smaller nested domain is Med-CORDEX at 20 km resolution. The southwestern Europe region studied here is shown in red.

Chemical boundary conditions for the larger domain are provided by the climatology from the global CTM LMDZ4-INCA3 (Hauglustaine et al.2014), for trace gases and non-dust aerosols, and from the GOCART model (Chin et al.2002), for dust.

The biogenic emissions are calculated online by the MEGAN model (see Sect. 3.2.2). Sea salts and dimethyl sulfide marine emissions are calculated online using the scheme of Monahan (1986) and Liss and Mervilat (1986), respectively. Mineral dust emissions are also calculated online (Marticorena and Bergametti1995; Alfaro and Gomes2001). Consistent with the WRF model, the land cover classification used for the calculation of biogenic and dust emissions is from the United States Geological Survey (USGS) land cover database (, last access: 14 April 2022). The biomass burning emissions from the APIFLAME v2.0 model (Turquety et al.2020) are used. The anthropogenic surface emissions are derived from the European Monitoring and Evaluation Programme (EMEP) database at 0.1× 0.1 resolution (, last access: 14 April 2022) for the year 2014. The inventory was preprocessed by the emiSURF program in order to calculate hourly emissions fluxes on the CHIMERE grid and the period of simulations (, last access: 14 April 2022).

3.2.2 MEGAN and the soil moisture factor

BVOC emissions were computed using the MEGAN model v2.04 (Guenther et al.2006), including several updates from the version 2.1 (Guenther et al.2012). It is used online in the CHIMERE model at an hourly time step. Emission fluxes are calculated based on emission factors at canopy standard conditions that are defined as 303 K for the air temperature (T), 1500 µmol m−2 s−1 for the photosynthetic photon flux density (PPFD) at the top of the canopy, 5 m2 m−2 for the LAI, and a canopy composed of 80 % mature, 10 % growing, and 10 % old foliage. Environmental conditions are taken into account using activity factors (γ) that represent deviations from the canopy standard conditions, so that the emission rate (ER; µg m−2 h−1) for a given species i is calculated in each model grid cell as follows:

(3) ER i = EF i γ LAI γ PPFD , i γ T , i γ Age , i γ SM , i γ CO 2 , i ρ i ,

where EF is the emission factor (µg m−2 h−1) provided by the MEGAN model, and ρ is the production/loss term within canopy that is assumed to be equal to 1.

γLAI is the activity factor based on LAI observations from the MODIS MOD15A2H product (Myneni et al.2015) improved by Yuan et al. (2011) (, last access: 14 April 2022). This improvement is undertaken with a two-step integrated method, where (1) the modified temporal spatial filter is used to fill the gaps and replace low-quality data with consistent ones, and (2) the TIMESAT (a program for analysing the time series of satellite sensor data) Savitzky–Golay filter is applied to smooth the final product. The temporal resolution is 8 d.

γPPFD, γT, and γAge are the activity factors accounting for light, temperature, and leaf age, respectively. They are calculated in MEGAN using meteorological variables from the WRF model. The activity factors accounting for the soil moisture (γSM) and CO2 concentration (γCO2) are applied on C5H8 species only. γCO2 is fixed (CO2 concentration at 395 ppm – parts per million), and γSM depends on the deviation between the soil wetness and a fixed wilting point, which is computed as follows:

(4) γ SM = 1 if   θ > θ 1 γ SM = ( θ - θ w ) Δ θ 1 if   θ w < θ < θ 1 γ SM = 0 if   θ < θ w ,

where θ is the soil wetness (volumetric water content; in m3 m−3), θw is the wilting point (level from which plants cannot extract water from the soil; in m3 m−3), Δθ1 is equal to 0.04 (empirical parameter), and θ1 is the sum of θw and Δθ1. θ values (third layer; 60 cm depth) are provided by the WRF-Noah land surface model (Ek et al.2003; Greve et al.2013). θw values are computed over the domain (Fig. S1 in the Supplement) using tabulated values from Chen and Dudhia (2001) that are soil type specific and parameterized on Noah soil wetness values. These values are spatialized using the soil texture map provided by the USGS (STATSGO-FAO product). However, emission response to drought is very sensitive to the wilting point (Müller et al.2008). Y. Wang et al. (2021) computed a relative difference in biogenic emissions varying between 50 % and 90 %, depending on the wilting point values.

Based on a comprehensive data set of in situ measurements of C5H8 emissions, Bonn et al. (2019) computed an activity factor of soil moisture as a function of the soil water availability. The fitted function is in relatively good agreement with the γSM algorithm in MEGAN.

As an alternative, Jiang et al. (2018) propose using an activity factor of soil moisture that integrates the soil water stress for plant photosynthesis. Below a critical value (fixed here at 0.5), the activity factor using SWS from RegIPSL (γSWS) decreases linearly as follows:

(5) γ SWS = 2 SWS if   SWS 0.5 γ SWS = 1 if   SWS > 0.5 .

The different soil wetness functions presented here will be subject to dedicated simulations for comparison.

Soil NO emissions are also included in the different simulations. Emission factors are based on the European soil emission inventory (Stohl et al.1996), which includes contributions from both forests and agricultural soils. Due to the strong temperature dependence, activity factors of soil NO are processed as for biogenic emissions in the MEGAN model. However, no dependence on soil moisture was parameterized.

3.2.3 Dry deposition scheme

Dry deposition in CHIMERE is resolved online using a canopy resistance approach, based on the scheme of the EMEP model (Wesely1989; Emberson et al.2000; Simpson et al.2003, 2012). Canopy resistance is calculated (Eq. 6) from stomatal conductance (gsto), which increases proportionally with LAI, and the bulk non-stomatal conductance (Gns). gsto varies between a maximum (gmax) and minimum (gmin) daytime stomatal conductance (mmol O3 m−2 s−1), depending on a series of meteorological conditions represented by factors of relative conductance (f). The parameters used to calculate these factors vary with the land cover (lc) and the seasonality. Because the land cover classes used in the dry deposition scheme are different from the one used in CHIMERE (USGS), a matrix of conversion is applied.

(6) g sto = g max , lc f phen , lc f light , lc max g min , lc , ( f temp , lc f VPD , lc f SWS , lc ) ,

where fX are factors of relative conductance determined by the leaf/needle age (X=phen), irradiance (X=light), temperature (X=temp), leaf-to-air vapour pressure deficit (X=VPD), and the soil moisture (X=SWS). Environmental variables required for fX factors are from the WRF model, except for fSWS that is based on SWS from RegIPSL. fSWS, for which the calculation is the same as γSWS, follows the parameterization suggested by Simpson et al. (2012) for the stomatal conductance.

Gns depends on three resistances, namely the external leaf uptake, the ground surface, and the in-canopy resistance. The external leaf uptake resistance is fixed at 2500 s m−1, the ground surface resistance is estimated from tabulated values (PFT specific), and the in-canopy resistance varies with the surface area index which is expressed in terms of LAI. The LAI used in the dry deposition scheme is parameterized (with no interannual variation). Information on phenology and biomass variation for each land cover type was collected from several studies in Emberson et al. (2000).

3.3 Modelling experiments

The years 2012, 2013, and 2014 were selected for the large diversity of hydroclimatic conditions and the availability of O3 observations over the study area. During summer 2012, southwestern Europe was affected by severe droughts (PLASD>0.1 of soil dryness) and heatwaves (PLAT2 m>1.5C; Fig. S2). The summer of 2013 was wet in the Iberian Peninsula, while summer 2014 was wet in Italy and the Balkans. A few heatwaves occurred in summer 2013 over southern France, northern Italy, and northern Spain. Summer 2014 was colder (PLAT2 m of 4.0 C) for the entire study area. The analysis method we propose here is divided into two steps.

Table 1Name, description, and aim of the simulations launched with the CHIMERE model on the nested Med-CORDEX domain for the summers of 2012, 2013, and 2014.

Download Print Version | Download XLSX

First, the focus is on drought effects (biomass decrease and soil dryness) that are generally poorly represented in modelling experiments over summer 2012. Several experiments are conducted on the Med-CORDEX domain in order to analyse the sensitivity of simulated C5H8 emissions, O3 dry deposition, and surface O3 concentration to drought effects. Table 1 summarizes the CHIMERE simulations conducted. The “Reference (R.)” simulations include all the emissions as defined in Sect. 3.2.1 but without accounting for the soil moisture factor in the emission scheme of MEGAN and deposition one of CHIMERE. A simulation with no biogenic emissions is also performed (“NoBio-emiss”) to estimate their relative contribution to O3 concentration.

The LAI used in the emission scheme of MEGAN is year dependent (MODIS observations). To evaluate the effect of biomass decrease by droughts, a simulation with LAI corresponding to a wet summer was used (“HighLAI-emiss”). Summer 2012, which was affected by an important biomass decrease over most of the study area (Fig. S3), has been simulated with the LAI of the wet summer of 2014 (higher than the 2012–2014 mean).

The LAI database used in the dry deposition scheme of CHIMERE does not vary with the year, so the biomass decrease during the summer 2012 is not reflected. In order to evaluate the importance of this effect, a simulation (“LAIdecr-dep”) with a LAI decrease of 5 % for forests and 20 % for grass and crops (close to the variations observed in 2012; Fig. S3) has been conducted.

The effect of soil dryness for the BVOC emissions is evaluated using the following three different approaches:

  • γSM from WRF-Noah in MEGAN (“γSM-emiss” experiment),

  • γSWS from RegIPSL in MEGAN (“γSWS-emiss” experiment), and

  • γSWSfit from RegIPSL in MEGAN (“γSWSfit-emiss” experiment). The fitted function of Bonn et al. (2019) is applied on the SWS as follows:

    (7) γ SWSfit = exp - exp 0.056 exp ( 1 ) - 2.3 - SWS + 1 .

The soil dryness effect on the dry deposition is evaluated using the fSWS from RegIPSL in CHIMERE (“fSWS-dep” experiment).

Second, based on a simulation configuration integrating all drought and heatwave effects (γSWSfit-emiss and LAIdecr/fSWS-dep in the same experiment, hereafter called the “all-emiss-dep” experiment), the variation in C5H8 emissions, O3 dry deposition velocity, and surface O3 concentration is analysed over the summers of 2012, 2013, and 2014. In order to compare isolated and combined extreme events with normal conditions, we used a cluster approach based on the PLAT2 m and PLASD indicators.

4 Validation of surface O3, NO2, and T2 m

Simulated surface O3 and NO2 concentrations by CHIMERE were compared to the EEA observations, while simulated T2 m by the WRF model was compared to the E-OBS observations (Cornes et al.2018). Table 2 presents the mean validation scores of the reference simulations over southwestern Europe for summer 2012. Validation scores of O3 and NO2 are very similar for the summers of 2013 and 2014.

Table 2Comparisons between observed and simulated surface O3 concentration, NO2 concentration, and 2 m temperature, averaged over southwestern Europe during summer 2012, for the reference CHIMERE simulation.

Download Print Version | Download XLSX

The simulated surface O3 is slightly higher than the observations for hourly and daily mean values (observation–model, or obs.–mod, bias of 1.93 and 0.13 µg m−3) and lower for the daily maximum values (+10.94 µg m−3). The simulated diurnal cycle is lower than the observed one because CHIMERE overestimates the daily minimum and underestimates the maximum. The temporal correlation (R, Pearson's coefficient) varies between 0.53 in the daily mean and 0.57 in the hourly mean. The spatial distribution of bias and the correlation scores of the daily maximum values for each summer is shown in the Supplement (Fig. S4). The highest bias (>+25 µg m−3) and lowest correlations (R of about 0.2) are obtained near large urban areas (e.g. Madrid and Milan), which may be less well represented in the model, even if the stations are classified as rural background. There are substantial uncertainties related to the NOx/ VOC emissions in (peri-)urban areas.

These O3 scores show an overall agreement with those calculated by Gaubert et al. (2014) and Menut et al. (2021) over Europe. The averaged daily maximum bias and correlation coefficients computed are slightly higher and lower, respectively, in this study mainly because the area is limited to southwestern Europe. Due to the multiple sources of O3 precursors and favourable temperature and light conditions, the Mediterranean is identified as a region of important uncertainty for modelling O3 concentration (Richards et al.2013).

Validation scores are also computed for the surface NO2 using the EEA observations, with a different distribution of stations that for O3. In agreement with the results of Menut et al. (2021), CHIMERE underestimates NO2 concentrations (+4.95 µg m−3 in the daily mean and +8.57 µg m−3 in the daily maximum) compared to the EEA observations. Interference from other oxidized nitrogen species in the chemiluminescence analyser (e.g. Lamsal et al.2008) may explain the constant underestimation compared to observations. The mean temporal correlation coefficients do not exceed 0.4. These low validation scores may also be linked to the emission inventory used in addition to the meteorological conditions (boundary layer height representation in particular). As a major O3 precursor, part of the uncertainty in this study is directly related to NO2 emissions and concentrations.

Finally, as a critical variable for simulating diurnal and seasonal O3 concentrations, T2 m values have been compared to the E-OBS observations. Averaged over southwestern Europe, the bias is close to 0 C ,while the RMSE is significant (8.19 C in the daily mean and 12.63 C in the daily maximum). The spatial distribution of the bias presents large variations (Fig. S4). The daily maximum T2 m is overestimated in the southern Mediterranean (up to 5 C), compared to observations, while it is underestimated in the northern part (up to 5 C). Nevertheless, the averaged scores of the temporal correlation are high (around 0.75). Even if such validation scores are close to those found in the scientific literature (e.g. Panthou et al.2018), the temperature uncertainties significantly contribute to those of the O3 simulated by CHIMERE.

5 Results

5.1 Sensitivity to soil dryness and biomass decrease effects

The sensitivity analysis to soil dryness and biomass decrease effects is focused on summer 2012, as droughts mainly occurred during this year (over the 2012–2014 period). Most of southwestern Europe was affected by a biomass decrease in 2012 (compared to the wet summer of 2014), except for eastern Spain. The LAI difference between summer 2012 and 2014 can reach 30 % over central Italy (Apennines region), corresponding to a mean decrease of 0.5 m2 m−2 (Fig. S3).

C5H8 is the main contributor to the total mass of BVOCs emitted (70 %) in our study area, followed by the model species APINEN (13 %, including, e.g., α-pinene and sabinene).

5.1.1 C5H8 emissions

Figure 2a shows the daily mean C5H8 emissions in southwestern Europe for summer 2012 simulated by the MEGAN model (reference simulation). The spatial distribution of C5H8 over the Mediterranean is heterogeneous, ranging from areas of zero or low emissions (1×10-3 g m−2 h−1) to high-emitting areas (1.5×10-3 g m−2 h−1), such as in the Balkans, Apennines, Sierra Morena, Sardinia, and central Portugal. Similar spatial distributions were found for APINEN model species. Differences in spatial patterns among species depend on the variation in emission factors over the land cover classes.

Figure 2Summer mean (June–August, JJA) of the daily mean C5H8 emissions (g m−2 h−1) for 2012 from the reference simulation (a). Mean relative difference (%) in C5H8 emissions due to the biomass decrease and soil dryness. The relative difference is computed between the reference and HighLAI-emiss simulation to quantify the biomass decrease effect (b) and the γSM-emiss (c), γSWS-emiss (d), and γSWSfit-emiss (e) simulation to quantify the soil dryness effect (for each configuration). For diagnostic purposes, areas of interest are designated within the dashed rectangles in panel (b), which are the Balkans (A), Po Valley (B), and central Spain (C).

The effect of biomass decrease and soil dryness on C5H8 emissions during the dry summer 2012 is also shown in Fig. 2. The biomass decrease averaged over June, July, and August is characterized by negative differences in emissions over most of the southern part of the region, reaching 20 % in northern Spain and 25 % in northern Italy (Fig. 2b).

The impact of soil dryness is assessed with the γSM-emiss, γSWS-emiss, and γSWSfit-emiss MEGAN simulations (Fig. 2c, d, and e, respectively). The quantified impact can vary considerably between them. In the γSM-emiss experience, C5H8 emissions stay constant as long as the soil moisture is above the wilting point (θw). Once this point is reached (water is missing), it is assumed that plant cannot synthesize C5H8 anymore, and the emissions decrease steeply. θw is constant and depends only on the soil type. For instance, silty soil is characterized by lower θw than clayey soil, so that soil dryness is more quickly reached for the latter. Due to the high spatial variability in soil type, the difference induced by γSM-emiss is also characterized by large variability. Within a same region (e.g. central Italy), the difference sharply varies between 50 % and 0 %, while the region is affected all over by an agricultural drought, according to the PLASD indicator (Fig. S2).

For the alternative approach based on the dynamical SWS function, the relative difference in C5H8 emissions (reference vs. γSWS-emiss experiment) is more spatially homogeneous than γSM-emiss, and the overall reduction is larger. The strongest stress values are located in plains and for plants with short root systems, which is in agreement with the sensitivity analysis performed by Vicente-Serrano (2007) about drought effects on vegetation. Semi-arid regions (e.g. Andalusia) are strongly affected, as the water budget is almost permanently in deficit (high solar radiation and low/no precipitation). Having adapted to recurrent droughts, some plant species (e.g. Arundo donax in a Moroccan ecotype) can reduce the C5H8 synthesis as the result of pressure selection to preserve their viability (Haworth et al.2017). However, the emission reduction could be overestimated, as irrigation is largely used in many semi-arid areas (e.g. García-Vila et al.2008). A specific analysis should be undertaken to cover the diversity of responses to drought stress in both natural and human-influenced systems within such regions.

The third experiment, γSWSfit-emiss, presents the same areas of C5H8 decrease as γSWS-emiss but to a lesser extent. C5H8 emissions decrease when SWS values are the lowest. Semi-arid regions are also more affected (50 %) than others (e.g. 20 % in central Italy).

Figure 3Daily mean C5H8 emissions (g m−2 h−1) during summer 2012, spatially averaged over the Balkans (a), Po Valley (b), and central Spain (c), for the different MEGAN experiments. The coloured bands highlight periods of droughts and heatwaves.


Three areas of high emissions are analysed more specifically (see the boxes in Fig. 2b), namely the Balkans, the Po Valley, and central Spain. Their temporal evolution of C5H8 emissions during summer 2012 for the different experiments is shown in Fig. 3. Droughts and heatwaves occurred over the three regions. Droughts do not negate the dependence on light and/or temperature, and emission peaks are driven by heatwaves with higher 2 m temperature and solar radiation (Fig. S6). However, the peak values are reduced if LAI decrease and soil dryness are accounted for. When averaged over the three regions, the biomass decrease induces a 3 % decrease in the total summer amount of C5H8 emitted. The soil dryness parameter induces difference of 12 % for γSM-emiss, 39 % for γSWS-emiss, and 13 % for γSWSfit-emiss. Based on experimental measurements over three summers (2012, 2013, and 2014) in a Mediterranean environment (Observatoire de Haute-Provence; Quercus pubescens plant species), Saunier et al. (2017) identified a 35 % emission decrease due to severe droughts. In California, Demetillo et al. (2019) measured a C5H8 concentration reduction of 50 % during severe droughts (2014 and 2015). These values are close to, and even above, the range of the simulated experiments. The γSWS-emiss experiment is considered here to be the upper limit of the reduction range of C5H8 emissions due to dry conditions, with SWS being included in the calculation of the emissions rate as soon as soil water stress for plant stress is below 0.5.

Figure 4Observed surface C5H8 concentration (pmol mol−1) during summer 2012 at the Ersa station (FR0033R; Cape Corsica) from the EBAS data set, compared to the different simulated experiments undertaken by the MEGAN-CHIMERE model. The shaded curve in black represents the precision range of the measurements. The reference and γSM-emiss experiments here have equal values. Maximum PLAT2 m and PLASD are +3.2 C and 0.02 of soil dryness index, respectively.


In addition, the different simulated experiments from the MEGAN-CHIMERE model have been compared to observations of surface C5H8 concentrations from EBAS data set (Fig. 4). The Ersa station, having an almost complete time series for summer 2012, is kept. Based on the PLA indicator, northern Corsica was affected by two heatwaves but no drought. However, we detected soil dryness conditions close to the drought limit (mean PLASD of 0.09). The temporal correlation between simulated and observed C5H8 concentrations is similar for each experiment (R coefficient around 0.68). Averaged over the summer, the lowest mean bias with the observations is obtained with γSWSfit-emiss (28 pmol mol−1) and the largest one with γSM-emiss (151 pmol mol−1). Over July and August, the γSWS-emiss experiment presents the lowest mean bias (+56 pmol mol−1). The reference and γSM-emiss experiments have equal values, as the soil wetness from WRF-Noah is above the local wilting point. However, this comparison with observations is made at a single station. It would be valuable to carry out the same exercise at several locations (including areas representative of a semi-arid environment) and over several summers.

Finally, biogenic emission models (such as MEGAN) consider a reduction in emissions in general when drought episodes occur only for isoprene (C5H8). However, there is a lot of uncertainty regarding how plant activity reacts to water stress. In the case of C5H8 species, the response to water stress could occur in two phases, i.e. a state of increasing emissions due to leaf temperature stimulation during the drought onset, followed by a state of C5H8 synthase limitation (Potosnak et al.2014). For monoterpenes, the response to water stress may be different, depending on the species (Bonn et al.2019). For instance, sabinene emissions (experimentally measured) from the European beech (Fagus sylvatica) decrease strongly with decreasing soil water availability, while trans-β-ocimene emissions (from the same plant species) remain constant.

5.1.2 O3 dry deposition velocity

Figure 5a presents the daily mean dry deposition velocity ranging from 1.0 to 1.8 cm s−1 within our areas of interest. The mean effect of biomass decrease (relative difference between the reference and the LAIdecr-dep experiment; Fig. 5b) over summer 2012 does not exceed 8 % and is lower over the forested areas (as it was prescribed in our experiment).

Figure 5Summer mean (JJA) of the daily mean O3 dry deposition velocity (cm s−1) for 2012 from the reference simulation (a). The mean relative difference (%) of the O3 dry deposition velocity is due to biomass decrease and soil dryness. The relative difference is computed between the reference and LAIdecr-dep simulation to quantify the biomass decrease effect (b) and LAIdecr-dep/fSWS-dep to quantify the combined biomass decrease and soil dryness effect (c).

The cumulative effects of biomass decrease and soil dryness are included in the LAIdecr-dep/fSWS-dep experiment (Fig. 5c). The soil dryness effect is much larger than that of the biomass decrease. Values range from 0 % to 35 %. Central Italy and the Iberian Peninsula are the most affected regions.

Figure 6Daily mean O3 dry deposition velocity (cm s−1) during summer 2012, spatially averaged over the Balkans (a), Po Valley (b), and central Spain (c) for the different CHIMERE experiments.


The temporal evolution of O3 dry deposition velocity for the three areas of interest is shown in Fig. 6. The biomass decrease effect is larger (e.g. 9 % over the Po Valley) during the first half of the summer. It is close to 1 % at the end of the summer. This is explained by the LAI decrease prescribed from mid-July in the dry deposition scheme. In addition to that, the dry deposition velocity is characterized by a decreasing trend over summer imposed by the fixed phenology factor (fphen).

Regarding the soil dryness, its effect is almost constant during the summer (e.g. 12 % over the Po Valley). Based on measurements over central Italy, Lin et al. (2020) computed a relative difference of about 50 % between August 2004 (wet summer;  0.8 cm s−1) and August 2003 (dry summer;  0.4 cm s−1). However, summer 2003 was characterized by considerable heatwaves in Italy (positive PLAT2 m), which might intensify the decrease. We emphasize that dry deposition velocity generally decreases during heatwaves (Fig. 6), as the near-surface temperature is above the optimal temperature of stomatal conductance. Both effects related to droughts and (the most intense) heatwaves lead to a reduction in O3 deposition. The sensitivity of the dry deposition velocity to soil moisture can be considerably different from one deposition scheme to another. Using the same deposition scheme as in the present study, Anav et al. (2018) calculated an average decrease of 10 % over Europe in dry O3 deposition.

Figure 7Observed O3 deposition flux (g cm−2) during summer 2012 at the Castel Porziano station (IT-Cp2; Lazio region) from the EFDC database, compared to the different simulated experiments undertaken by the CHIMERE model. The maximum PLAT2 m and PLASD are +3.88 C and +0.09 of soil dryness index, respectively.


To the best of our knowledge, there is no study yet in the scientific literature that fully assess the O3 dry deposition of CHIMERE against observations or its sensitivity to different meteorological forcings. Therefore, we have compared our results to the measured data available for summer 2012 at the Castel Porziano station from the European Eddy Fluxes Database Cluster (Fig. 7). Based on the PLA indicator, the Lazio region was affected by a severe drought throughout the summer (mean PLASD of +0.05). All simulated experiments overestimate O3 deposition compared to observations. The LAIdecr-dep/fSWS-dep experiment has the smallest average bias (0.19 × 10−6 g cm2) and the highest correlation (R coefficient of 0.40). Such overestimation has also been calculated for models whose gas deposition scheme is based on the description in Wesely (1989) of the big leaf parameterization (e.g. Michou et al.2005; Huang et al.2022). As the canopy conductance increases proportionally with the prescribed LAI (Emberson et al.2000), this could be explained by an overestimation of the LAI that is almost 2 times larger than the mean LAI reconstructed from MODIS over this area. The importance of representing processes dynamically (as opposed to fixed parameters, especially for the non-stomatal conductance) is also highlighted in order to better simulate the diurnal deposition cycle and the daily average values (Huang et al.2022).

Finally, identifying soil dryness as a major driver of interannual variability in O3 deposition velocity, Lin et al. (2019) emphasize that the error in modelling the deposition may considerably rely on the ability of the models to accurately simulate the precipitation distribution. They assessed the sensitivity of their deposition velocity scheme (Geophysical Fluid Dynamics Laboratory, LM3.0/LM4.0) to two different meteorological forcings and found a factor of 2 over northern Europe.

5.1.3 O3 surface concentration

Figure 8a maps the mean surface O3 concentration simulated by CHIMERE (reference experiment) over summer 2012. High concentrations (between 100 and 110 µg m−3) are located in the eastern part of the study area (e.g. the Balkans and Italy). The spatial distribution of the mean surface O3 shows similar patterns to the distribution of PLAT2 m (Fig. S2), highlighting the critical role of the temperature. Surface O3 remains high above the sea due to transport and the absence of deposition in the canopy (e.g. 120 µg m−3 over the Adriatic sea). The contribution of BVOC emissions on surface O3 is also presented (Fig. 8b). BVOC emissions are known to be significant precursors of O3 production in the Mediterranean region. It varies between 4 % (e.g. Iberian Peninsula) and 22 % (e.g. northern Italy) in our study, and that is included in the range of values reported in the scientific literature (e.g. Mertens et al.2020). Finally, our areas of interest are located in rural regions mostly characterized by a low-NOx regime (Fig. S7).

Figure 8Summer mean (JJA) of daily mean O3 surface concentration (µg m−3) for 2012 from the reference simulation (a). The mean contribution of biogenic emissions to O3 surface concentration is based on the NoBio-emiss experiment (b). The mean (c) and maximum absolute (d) relative difference (%) of the O3 surface concentration due to the biomass decrease and soil dryness. The relative difference is computed between the reference and γSWS-emiss and LAIdecr-dep/fSWS-dep simulation from CHIMERE.

A simulation including drought effects on both C5H8 emissions and O3 dry deposition has been conducted with MEGAN-CHIMERE (γSWS-emiss and LAIdecr/fSWS-dep experiment). The resulting effect on surface O3 is shown in Fig. 8. On average, during the summer (Fig. 8c), the O3 concentration is slightly higher over the continent (between +0.5 % and +3.0 %) due to the decrease in O3 deposition as a dominant effect, while it is slightly lower (between 1 % and 1.5 %) over the sea and ocean due to the lower transport of O3 precursors, compared to the reference experiment. However, the O3 increase over the Iberian Peninsula may be overestimated, as the LAI reduction we applied to the whole domain (e.g. 20 % for grass PFTs) is larger compared to the variation in MODIS LAI in this specific region (Fig. S3).

Drought effects on surface O3 induced by the C5H8 emission reduction is not constant. It is largest during combined heatwaves (i.e. when the biogenic contribution is high). As a result, drought effects on C5H8 emissions can be dominant (compared to the effects on O3 deposition) during simultaneous droughts and heatwaves, inducing a decrease in O3 peaks by a few micrograms per cubic metre (µg m−3) over both the continent and sea/ocean. The maximum absolute relative difference (Fig. 8d) reaches 5 % over the Po Valley and 14 % along the Strait of Gibraltar. O3 formation over the latter is favoured by large NOx shipping emissions.

Conducting a similar modelling experiment based on the γSM from MEGANv3, Jiang et al. (2018) simulated a maximum absolute relative difference in surface O3 of about 4 % in August 2010. O3 reduction (10 % on average) due to severe droughts was also measured in California over the period 2002–2015 (Demetillo et al.2019). This was identified as being related to a steep decrease in C5H8 concentrations.

Figure 9Daily mean O3 surface concentration (µg m−3) during summer 2012 that is spatially averaged over the Po Valley from the EEA observations and the different CHIMERE experiments.


For each area of interest, the temporal evolution of surface O3 based on the different CHIMERE experiments and the EEA observations (AQ e-Reporting) is presented (Fig. 9 for the Po Valley and Fig. S8 in the Supplement for the Balkans and central Spain). Since the resulting effects of biomass decrease and soil dryness on surface O3 are less than the bias between model and observations in our study (see Sect. 4), no simulation can be designated with certainty as the best fit.

5.2 Statistical variation during droughts and heatwaves

In this result section, a statistical analysis is performed using CHIMERE simulations and several observational data sets. The all-emiss-dep experiment has been chosen for simulating summers 2012–2014 because it includes drought and heatwave effects in the most comprehensive way. Moreover, the C5H8 emission approach γSWSfit-emiss has shown a good performance compared to observations (Fig. 4), remaining rather conservative and not being in the upper limit of the C5H8 reduction range.

Clusters of droughts and heatwaves (isolated or combined) are constructed based on the PLAT2 m and PLASD indicators, allowing the analysis of the statistical variation in simulated C5H8 emissions, O3 stomatal conductance, and O3 surface concentration for the summers of 2012–2014. We performed the same analysis on observations of HCHO total column for the summers of 2005–2016 and O3 the summers of 2000–2016. Using those indicators, the following conditions were defined and grouped into the following clusters: (1) heatwaves, (2) droughts, (3) heatwaves or droughts, (4) heatwaves and droughts, (5) heatwaves and not droughts, and (6) droughts and not heatwaves. Normal conditions are defined as (7) no drought and no heatwave.

5.2.1 C5H8 emissions

Figure 10 shows the distribution of C5H8 emission rates over the southwestern Mediterranean for clusters of extreme weather events. On average, the daily maximum C5H8 emission is significantly higher (t test; p<0.01) during isolated and combined heatwaves than normal conditions (no heatwave nor drought), with a mean value of 0.16×10-2 against 0.12×10-2 g m−2 h−1. However, emissions during droughts have the same mean value as in normal conditions. During isolated droughts (cluster drought and not heatwave), the mean daily maximum C5H8 emission rate is lower than normal conditions (-0.02×10-2 g m−2 h−1; non-significant difference). These results are in general agreement with the observed HCHO total column by satellite instrument, used as proxy for BVOC emissions variation (see Sect. 5.2.2). Weather conditions are considerably different between the droughts combined with heatwaves and those which are not (Fig. S9). The combined (respectively, isolated) droughts are characterized by a 2 m temperature of 23.9 C (respectively, 22.6 C), a shortwave radiation of 342.2 W m−2 (respectively, 296.0 W m−2), and a cloud cover of 1.8 % (respectively, 3.1 %). Those weather variables are used for the computation of the activity factors γP and γT, thus directly affecting emissions. During isolated droughts, γLAI (0.48) and γSWS (0.76) are smaller than for normal conditions (0.51 and 0.89, respectively, with a significant difference for both). Nevertheless, this negative signal is not large enough for a significant variation in the emission rates (compared to normal conditions). Gathering worldwide data from experimental measurements of biogenic emissions under different climate drivers, the scientific review presented by Feng et al. (2019, based on 74 articles) estimated a +53 % emission change in C5H8 during warm conditions and a 15 % change during dry conditions. Those variations are close to what we simulated (+35 % and 13 %, respectively).

Figure 10Daily maximum C5H8 emission rate (g m−2 h−1) simulated by the MEGAN model (γSWSfit-emiss experiment) over southwestern Europe, clustered by identified extreme weather events (PLAT2 m and PLASD indicators from RegIPSL). The number of days is indicated in parentheses. The analysed period is June–August in 2012, 2013 ,and 2014, covering a total of 276 d. The red squares show the mean of the distribution, and the black circles are the outliers. The blue dotted line indicates the mean value of the normal conditions (no hw. nor dr. cluster), and the blue stars show if the mean value of the considered cluster is significantly different (t test; at least p<0.1) from the normal conditions. The box covers the interquartile range (IQR) between Q1 (25th percentile) and Q3 (75th percentile). The lower whisker is limited to a statistical minimum (Q1  1.5 × IQR) and the upper one to a statistical maximum (Q3 + 1.5 × IQR).


5.2.2 HCHO total column

HCHO is a product of the oxidation of VOCs. Based on the reference and NoBio-emiss CHIMERE simulations, we computed that biogenic emissions contribute between 60 % and 80 % of the HCHO concentration over southwestern Europe. Variations in the HCHO concentration may therefore be used as indicator of BVOC emission variations during droughts and heatwaves. This allows us to use satellite observations of HCHO, which is particularly interesting due to the lack of in situ data. Observations of the HCHO total column from the OMI instrument are used.

Table 3Variation in HCHO total atmospheric column (molec. cm−2) (ΔHCHO) due to heatwaves, droughts, and isolated droughts in comparison to normal conditions (norm. cdt.; no heatwaves nor droughts) for the summers (JJA) between 2005–2016 (measurements at 13:00 LT, local time). Summer 2012 is compared with CHIMERE simulations. Results are computed for each pixel and averaged over the following areas of interest: the Balkans, central Italy, and central Spain. The asterisks mean that the difference with normal conditions is statistically significant (t test; p<0.1).

Download Print Version | Download XLSX

Table 3 presents the average difference in HCHO (ΔHCHO) during extreme events compared to normal conditions (no heatwave nor drought), based on the OMHCHOd product. Over the summers of 2005 to 2016, HCHO is significantly higher during heatwaves for the Balkans, Po Valley, and central Spain (+15 % on average). ΔHCHO is also positive during droughts but to a lesser extent (+3 % on average, non-significant for central Spain only). However, isolated droughts induce a significant decrease for the Balkans (7 %), Po Valley (6 %), and central Spain (6 %). Those results are consistent with the variation in C5H8 emissions simulated by MEGAN when a soil moisture parameter is considered (see Sect. 5.2.1). P. Wang et al. (2021) also report a significant HCHO column decrease (up to 30 %) induced by a prolonged drought in a forested area in China.

The observed and simulated HCHO total columns have been compared for the γSWSfit-emiss experiment over summer 2012. The time series showing the HCHO evolution for the three areas of interest are shown in the Supplement (Fig. S11). Simulated HCHO columns are generally higher than those from the observations, especially for the Balkans. As mentioned in Sect. 2.2, uncertainty in the observations is large (between 30 % and 100 %), with a large spatial variability (about a factor of 2 larger than in CHIMERE). HCHO products from OMI also present high uncertainty and, in particular, a systematic low mean bias (20 %–51 %; Zhu et al.2016).

HCHO differences could be due to incorrect specifications of land cover and, thus, of the EF of C5H8 (e.g. Curci et al.2010). Temperate tree PFTs are characterized by high C5H8 emission rates (respectively, 10 000 and 600 µg m−2 h−1 for broadleaf and needleleaf types) compared to grassland (800 µg m−2 h−1) or cropland (1 µg m−2 h−1Guenther et al.2012). After the aggregation of the USGS land cover classes, the vegetation type assumed in CHIMERE in the Balkans, for instance, is 57 % forest cover, 9 % grassland, and 33 % cropland. Using the MODIS MCD12 product (Friedl et al.2010), we find a different distribution, with 30 % forest cover, 25 % grassland, and 31 % cropland (Fig. S10). The choice and the temporal evolution of the land cover database are crucial for the calculation of C5H8 emissions (Chen et al.2018). Despite the use of satellite data for land cover, significant uncertainties remain in the calculation of C5H8 emissions due to the classification of vegetation types and species (Opacka et al.2021).

The comparison of the HCHO total column variations for CHIMERE and OMI during droughts and heatwaves relies on a limited number of cases. It is therefore difficult to support conclusions with a sufficient level of certainty. Nevertheless, results suggest that CHIMERE simulations are more sensitive to temperature than the OMHCHOd OMI observations (+52 % and +28 % during heatwaves, respectively, averaged over the three areas of interest). Since the summer 2012 was affected by long agricultural droughts, the inclusion of a soil dryness parameter in CHIMERE (γSWSfit-emiss experiment) reduces the simulated HCHO peaks (Fig. S11). However, the mean simulated ΔHCHO during heatwaves and droughts is similar in the γSWSfit-emiss and reference experiments (decrease of 2 % and 1 %, respectively). The temporal correlation does not vary significantly either (e.g. R coefficient around 0.5 for the Balkans).

Finally, the observed variations in total HCHO over the summers of 2005–2006 and the three areas considered here were also computed with OMI-BIRA retrieval. During heatwaves, both show a significant increase, i.e. +31 % for OMI-BIRA (and +15 % for OMHCHOd). The increase is lowers by +13 % (and +3 %) during drought. The variation becomes slightly negative during isolated droughts, i.e. 2 % (and 6 %). This comparison not only highlights the uncertainty in the satellite observations but also that the general behaviour is consistent in both retrievals and similar to what was obtained for 2012 using CHIMERE.

5.2.3 O3 stomatal conductance

Surface weather conditions are critical for the stomatal conductance and therefore influence the dry deposition velocity. Figure 11 shows the maximum daily stomatal conductance of O3 (LAIdecr/fSWS-dep CHIMERE experiment) clustered by simulated extreme weather events and averaged over the western Mediterranean. The same analysis has been performed on the dry deposition velocity, and the signals induced by extreme weather events are similar.

Figure 11Same as Fig. 10 with the simulated O3 stomatal conductance (µmol O3 m−2 s−1) by the CHIMERE model (LAIdecr/fSWS-dep experiment).


Droughts and heatwaves (isolated or combined) induce a significant decrease in the O3 stomatal conductance, quantified at 25 % for heatwaves and 35 % for droughts, compared to normal conditions. The activity factors mainly affected by droughts and heatwaves are ftemp, fVPD, and fSWS.

The variation in ftemp during heatwaves depends on the magnitude of the events and also on their location, since the percentile of the PLA indicator is defined for each grid cell. Over the Po Valley and Balkans, for instance, most heatwaves are characterized by temperatures close to optimal values of stomatal conductance that are fixed around 24 C (Fig. S6). However, for those occurring in central Spain (between 30 and 32 C), ftemp decreases by 7 % compared to normal conditions. The temperature limit before complete stomatal closure is set at 40 C. Therefore, exceptional heatwaves occurring in southern Spain, for instance, could quickly lead to an accumulation of O3 at the surface. fVPD that depends both on temperature and relative humidity significantly decreases during droughts and heatwaves (e.g. 6 % averaged over southwestern Europe). Finally, fSWS is the factor dominating the signal of variation in stomatal conductance. At the scale of southwestern Europe, this factor is the lowest during isolated droughts, with a mean decrease of 35 %.

5.2.4 O3 surface concentration

Figure 12 shows the distribution of the observed (summers 2000–2016 and 2012–2014) and simulated (summers 2012–2014) daily maximum surface O3 concentrations over southwestern Europe for each cluster of extreme events. Observed O3 (Fig. 12a) is significantly (t test; p<0.01) higher during heatwaves (+18 µg m−3) and droughts (+9 µg m−3) than during normal conditions. Considering all droughts over the United States of America, Wang et al. (2017) also computed a mean increase in surface O3 (+17 % compared to the average). During isolated droughts in our study area, the daily maximum O3 is larger (+4 µg m−3), but the difference is non-significant.

Figure 12Same as Fig. 10 with the observed surface O3 concentration (µg m−3) over the summers of 2000–2016 (a) and 2012–2014 (b) and with the simulated surface O3 (all-emiss-dep experiment) over the summers of 2012–2014 (c).


The distribution of the simulated surface O3 (all-emiss-dep experiment) by extreme events over the period 2012–2014 (Fig. 12c) presents similar signals, but of lower magnitude, with +9 µg m−3 during heatwaves, +3 µg m−3 during droughts, and a non-significant difference during isolated droughts compared to normal conditions. Based on the results discussed above, the difference between the heatwave and isolated drought cluster could be explained by the different conditions of biogenic emissions, dry deposition, temperature, and light. However, observations over the same period (Fig. 12b) present a significant increase in the daily maximum O3 during isolated droughts (+9 µg m−3), unlike what we simulated. The C5H8 emission reduction effect during such an extreme event could be counterbalanced to a larger extent by the O3 dry deposition decrease. It could also be explained by an underestimated impact of the enhanced photochemistry in the simulations, as we simulated favourable weather conditions during both combined and isolated droughts.

In summary, the variation in canopy–troposphere interactions (simulated by the MEGAN and CHIMERE models) during droughts and heatwaves is characterized by a consistent signal with respect to O3 observations (except for the isolated droughts over summers 2012–2014). Meteorological conditions are critical for the O3 budget, especially during summer droughts and heatwaves. In addition to uncertainties in the modelling of precursor emissions and O3 deposition (as mentioned above), differences between observations and simulations may also rely on meteorological uncertainties, such as the diurnal temperature cycle (see Sect. 4) and the planetary boundary layer height (PBLH). The nighttime representation of the PBLH is often misrepresented in the WRF model (e.g. Chu et al.2019).

5.3 Threshold level exceedance of O3

In the European Union, the air quality standard for O3 exposure is set at a daily maximum concentration of 120 µg m−3 (8 h average; EEA2020). Figure 13 shows the number of days above this threshold based on the AQ e-Reporting observations and the CHIMERE all-emiss-dep simulations for the summers of 2012, 2013, and 2014. The summers of 2012 and 2013 present a large number of days exceeding the standard concentration. Even if exceedances occur in similar regions for observations and simulations, the number of days is generally larger in the observations. This is due to the overall underestimated daily maximum in CHIMERE compared to observations (see Sect. 4). For example, in the Po Valley that is the most affected region in southwestern Europe, around 60 exceeding days were observed and 50 d were simulated over summer 2012. This region is known for its highly polluted air (O3 and its precursors) due to high anthropogenic emissions and unfavourable topographic and meteorological conditions for pollutants dispersion (e.g. Bigi et al.2012). Other regions affected by O3 peaks can be highlighted, southeastern France, central Spain, and central Italy, for the summers of both 2012 and 2013.

Figure 13Number of days above the air quality threshold value from the European Union for surface O3 (maximum daily 8 h mean of 120 µg m−3) over the summers of 2012, 2013, and 2014 (JJA; 92 d in total). EEA observations (left column) are compared to CHIMERE (all-emiss-dep) simulations (right column).

Table 4Average percentage of stations with at least 1 exceedance day (first row) regarding EU standard during summer (JJA). Considering only those stations, the second row shows the average number of exceeding days (exc. days) per station and the lower rows the average distribution of days above the EU standard as a function of extreme weather events. The mean exceedance concentration is indicated in parenthesis.

Download Print Version | Download XLSX

Table 4 presents the distribution characteristics of stations with at least 1 d above the EU standard during summer period. On average, over 2000–2016, that concerns 54 % of EEA stations over the western Mediterranean. The average number of exceedance days per station is 27 d (almost a third of the summer period). Summer 2012 is above the 2000–2016 average, with 61 % of stations affected and 28 d on average per station. For the same year, co-located values from CHIMERE simulations (all-emiss-dep experiment) are lower, i.e. 57 % of stations and 18 d on average.

Over all the summers between 2000 and 2016, 34 % of the exceeding days occurred during heatwaves (with a mean exceeding value of 24 µg m−3) and 27 % during droughts (+18 µg m−3). The number of days decreases by 13 % if we consider only the isolated droughts (14 %; +15 µg m−3). Summer 2012 was affected by exceptional droughts and heatwaves, resulting in high O3 pollution. Around 80 % of the days above the EU threshold occurred during heatwaves or droughts for both observations and simulations.

6 Discussion and conclusions

The analyses presented in this study were organized around two main objectives. The first one was to assess the sensitivity of biogenic emissions, O3 dry deposition, and surface O3 to the biomass decrease and soil dryness effect in a CTM model. The extremely dry summer of 2012 was chosen, and simulations were performed using the MEGAN v2.1 and CHIMERE v2020r1 model. We showed that the soil dryness parameter is critical during drought events, decreasing the C5H8 emissions and O3 dry deposition velocity considerably. This effect has a larger impact than the biomass decrease. However, the resulting effect on surface O3 remains limited.

In addition to the soil moisture activity factor used in MEGAN v2.1 that is mainly based on the wilting point (γSM), we proposed an innovative activity factor based on a soil water stress function (γSWS) simulated by the land surface and vegetation model ORCHIDEE. The latter induces a larger reduction in C5H8 emissions with more homogeneous patterns that follow the drought indicator. Furthermore, we adjusted this factor with a function fitted (γSWSfit) from experimental measurements in Bonn et al. (2019). By comparing the simulated surface concentration of C5H8 with observations at the Ersa measurement site (Corsica), γSWSfit showed promising results. However, such an evaluation should be carried out over several sites and several years in order to determine the added value of this approach with more certainty.

The second objective of this study was to quantify the variation in surface O3 over southwestern Europe during agricultural droughts when combined, or not, with heatwaves. Those extreme weather events were identified based on the RegIPSL model using the PLA indicator. During the 2000–2016 period, 59 % of summer drought days were not accompanied by heatwaves (isolated droughts). For the summers of 2012–2014, analysed more specifically in this study, the 2 m temperature is on average 5.5 % lower during isolated droughts compared to all droughts, the shortwave radiation is 13.5 % lower, and the cloud fraction is 42 % higher. As a result, not only surface O3 but also BVOC emissions and O3 dry deposition velocity are substantially different if the drought considered is accompanied by a heatwave or not.

Based on a cluster approach using the PLA indicator, we showed that observed surface O3 (summers of 2000–2016) is larger by +18 µg m−3 in the daily maximum during heatwaves and by +9 µg m−3 during droughts, compared to normal conditions. Despite a difference of several micrograms per cubic metre (µg m−3), CHIMERE correctly simulates the variations in O3 concentration between the clusters of extreme events. The overall increase in surface O3 during both heatwaves and droughts would be explained by the O3 precursor emission enhancement, O3 dry deposition decrease, and favourable weather conditions, so that all of these mechanisms lead to an accumulation of O3. However, we simulated a decrease in C5H8 emissions during isolated droughts, resulting in a non-significant difference in surface O3 compared to normal conditions (from both observations and simulations). Despite a significant bias between the HCHO total columns simulated by CHIMERE and observed by OMI, the satellite data confirm an average increase in HCHO (+3 % to +13 %, depending of the product considered) over our three regions of interest for all droughts and a decrease (2 % to 6 %) for isolated droughts.

Finally, almost half of the summer days (2000–2016 period) exceeding the EU standard of O3 for air quality in southwestern Europe occurred during droughts or heatwaves. However, this percentage can increase (up to 80 %) for exceptionally dry and hot summers, as in 2012. Only 14 % of the exceedance days occurred during isolated droughts.

The implementation of the dynamical effects of droughts in the MEGAN-CHIMERE model contributes to a better representation of biosphere–atmosphere interactions. However, comparisons between simulated and observed surface O3 still show large discrepancies. Important uncertainties appear to be related to BVOC emissions (especially due to the land cover classification), to NOx concentrations, for which CHIMERE presents limited performance scores of validation, to O3 deposition, and, finally, to meteorological conditions (e.g. temperature and PBLH). For instance, the simulated daily maximum temperature that is underestimated in the northern Mediterranean, compared to observations (see Sect. 4), may induced a decrease in O3 precursors, especially as the emission–temperature relationship of BVOCs is exponential, with an optimal temperature for C5H8 species (Guenther et al.1993). Such uncertainties need to be addressed to improve the simulation of O3 during the summer, especially over southwestern Europe.

In addition, the critical role of soil NO emissions in O3 production is increasingly studied (e.g. Romer et al.2018), especially in rural areas (e.g. Sha et al.2021). As hydroclimatic conditions are critical for soil NO emissions, the O3 budget during droughts and heatwaves is likely to be significantly influenced by soil NO emissions. In this study, the dependence on soil dryness for NO emissions is not included. Emission pulses can occur when rain follows a drought, and emission factors are higher with dry soil than with wet soil (Steinkamp and Lawrence2011; Weng et al.2020).

Some recent knowledge on fundamental processes that allow a better representation of surface–atmosphere interactions during extreme weather events is not yet integrated by much, if not all, of the modelling community. Among the many examples is the increased emission of monoterpenes and sesquiterpenes during the development of drought (e.g. Bonn et al.2019; Peron et al.2021) or the in-canopy chemistry that is ignored or approximated by the big leaf model approach (e.g. Clifton et al.2020a). Several studies show that non-stomatal conductance in the Mediterranean counts as much as stomatal conductance in the O3 sink budget (e.g. Gerosa et al.2009; Sun et al.2022). Moreover, non-stomatal conductance also seems to be significantly affected during droughts and heatwaves. In contrast to stomatal conductance, droughts and heatwaves should have opposite effects, i.e. an increase in non-stomatal conductance during heatwaves and a decrease during droughts (Wong et al.2022). However, such changes are not taken into account in the EMEP deposition scheme implemented in CHIMERE. There is a real need for a better representation of stomatal conductance in deposition schemes.

In conclusion, we provide in this paper a detailed analysis of the drought and heatwave effects on modelled biosphere–troposphere interactions controlling surface O3 concentration, and our findings are supported by several observational data sets. Heatwaves, and droughts to a lesser extent, induce a significant increase in the surface O3. Soil dryness and biomass decrease, as specific effects of droughts, are critical for the variation in C5H8 emissions and O3 dry deposition over southwestern Europe. We emphasize the need for a more dynamical representation of interactions between vegetation, meteorology, and atmospheric chemistry in models in order to improve the simulation of summer O3.

Data availability

The data set of indicators of heatwaves and agricultural droughts (PLAT2m/SD) is freely available at (last access: 14 April 2022; Guion2022). The data from AQ e-Reporting can be found at (last access: 14 April 2022; EEA2022), from the MODIS and OMI instruments at (last access: 14 April 2022; MODIS2022) and (Chance2019), from NCEP at (last access: 14 April 2022; NCEP-NCAR2022), from EBAS at (last access: 14 April 2022; EBAS2022), and from EFDC at (last access: 16 October 2022, login required; European Fluxes Database Cluster2022).


The supplement related to this article is available online at:

Author contributions

AG, ST, and JP conceptualized the paper. AG, ST, AC, JP, AE, and JL developed the methodology. AG carried out the analysis and designed the figures. AG and ST wrote the original draft. AG, ST, AC, JP, AE, and JL reviewed and edited the paper. All authors have read and agreed to the final version of the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We acknowledge the EEA, for providing the AQ e-Reporting data set, and Guillaume Siour, for the data extraction and preparation. The MCD153AH and OMHCHOd products were retrieved from the NASA EOSDIS Land Processes Distributed Active Archive Center (LP DAAC) and USGS Earth Resources Observation and Science (EROS) Center. Finally, we acknowledge the NCEP for the reanalysis meteorological data and the Norwegian Institute for Air Research (NILU) for the EBAS database.

Financial support

This work has been granted access to the HPC resources of TGCC under an allocation (no. 10274) made by GENCI (Grand Équipement National de Calcul Intensif) and funded by the Sorbonne Université (SU) and the Centre National d'Études Spatiales (CNES).

Review statement

This paper was edited by Leiming Zhang and reviewed by three anonymous referees.


Alfaro, S. C. and Gomes, L.: Modeling mineral aerosol production by wind erosion: Emission intensities and aerosol size distributions in source areas, J. Geophys. Res.-Atmos., 106, 18075–18084,, 2001. a

Anand, J. S., Anav, A., Vitale, M., Peano, D., Unger, N., Yue, X., Parker, R. J., and Boesch, H.: Ozone-induced gross primary productivity reductions over European forests inferred from satellite observations, Biogeosciences Discuss. [preprint],, 2021. a

Anav, A., Proietti, C., Menut, L., Carnicelli, S., De Marco, A., and Paoletti, E.: Sensitivity of stomatal conductance to soil moisture: implications for tropospheric ozone, Atmos. Chem. Phys., 18, 5747–5763,, 2018. a

Bian, H. and Prather, M. J.: Fast-J2: Accurate Simulation of Stratospheric Photolysis in Global Chemical Models, J. Atmos. Chem., 41, 281–296,, 2002. a

Bigi, A., Ghermandi, G., and Harrison, R. M.: Analysis of the air pollution climate at a background site in the Po valley, J. Environ. Monit., 14, 552–563,, 2012. a

Bonn, B., Magh, R.-K., Rombach, J., and Kreuzwieser, J.: Biogenic isoprenoid emissions under drought stress: different responses for isoprene and terpenes, Biogeosciences, 16, 4627–4645,, 2019. a, b, c, d, e, f

Cao, J., Situ, S., Hao, Y., Xie, S., and Li, L.: Enhanced summertime ozone and SOA from biogenic volatile organic compound (BVOC) emissions due to vegetation biomass variability during 1981–2018 in China, Atmos. Chem. Phys., 22, 2351–2364,, 2022. a

Chance, K.: OMI/Aura Formaldehyde (HCHO) Total Column Daily L3 Weighted Mean Global 0.1deg Lat/Lon Grid V003, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC) [data set],, 2019. a

Chance, K.: OMI/Aura Formaldehyde (HCHO) Total Column Daily L3 Weighted Mean Global 0.1deg Lat/Lon Grid V003, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2019. a

Chen, F. and Dudhia, J.: Coupling an Advanced Land Surface–Hydrology Model with the Penn State–NCAR MM5 Modeling System. Part I: Model Implementation and Sensitivity, Mon. Weather Rev., 129, 569–585,<0569:CAALSH>2.0.CO;2, 2001. a

Chen, W. H., Guenther, A. B., Wang, X. M., Chen, Y. H., Gu, D. S., Chang, M., Zhou, S. Z., Wu, L. L., and Zhang, Y. Q.: Regional to Global Biogenic Isoprene Emission Responses to Changes in Vegetation From 2000 to 2015, J. Geophys. Res., 123, 3757–3771,, 2018. a

Cheng, S., Huang, J., Ji, F., and Lin, L.: Uncertainties of soil moisture in historical simulations and future projections, J. Geophys. Res.-Atmos., 122, 2239–2253,, 2017. a

Chin, M., Ginoux, P., Kinne, S., Torres, O., Holben, B. N., Duncan, B. N., Martin, R. V., Logan, J. A., Higurashi, A., and Nakajima, T.: Tropospheric Aerosol Optical Thickness from the GOCART Model and Comparisons with Satellite and Sun Photometer Measurements, J. Atmos. Sci., 59, 461–483,<0461:TAOTFT>2.0.CO;2, 2002. a

Chu, Y., Li, J., Li, C., Tan, W., Su, T., and Li, J.: Seasonal and diurnal variability of planetary boundary layer height in Beijing: Intercomparison between MPL and WRF results, Atmos. Res., 227, 1–13,, 2019. a

Clifton, O. E., Fiore, A. M., Massman, W. J., Baublitz, C. B., Coyle, M., Emberson, L., Fares, S., Farmer, D. K., Gentine, P., Gerosa, G., Guenther, A. B., Helmig, D., Lombardozzi, D. L., Munger, J. W., Patton, E. G., Pusede, S. E., Schwede, D. B., Silva, S. J., Sörgel, M., Steiner, A. L., and Tai, A. P. K.: Dry Deposition of Ozone Over Land: Processes, Measurement, and Modeling, Rev. Geophys., 58, e2019RG000670,, 2020a. a

Clifton, O. E., Paulot, F., Fiore, A. M., Horowitz, L. W., Correa, G., Baublitz, C. B., Fares, S., Goded, I., Goldstein, A. H., Gruening, C., Hogg, A. J., Loubet, B., Mammarella, I., Munger, J. W., Neil, L., Stella, P., Uddling, J., Vesala, T., and Weng, E.: Influence of Dynamic Ozone Dry Deposition on Ozone Pollution, J. Geophys. Res.-Atmos., 125, e2020JD032398,, 2020b. a

Cornes, R. C., van der Schrier, G., van den Besselaar, E. J., and Jones, P. D.: An ensemble version of the E-OBS temperature and precipitation data sets, J. Geophys. Res.-Atmos., 123, 9391–9409,, 2018. a, b

Couvidat, F., Bessagnet, B., Garcia-Vivanco, M., Real, E., Menut, L., and Colette, A.: Development of an inorganic and organic aerosol model (CHIMERE 2017β v1.0): seasonal and spatial evaluation over Europe, Geosci. Model Dev., 11, 165–194,, 2018. a

Curci, G., Palmer, P. I., Kurosu, T. P., Chance, K., and Visconti, G.: Estimating European volatile organic compound emissions using satellite observations of formaldehyde from the Ozone Monitoring Instrument, Atmos. Chem. Phys., 10, 11501–11517,, 2010. a, b, c

De Andrés, J. M., Borge, R., de la Paz, D., Lumbreras, J., and Rodríguez, E.: Implementation of a module for risk of ozone impacts assessment to vegetation in the Integrated Assessment Modelling system for the Iberian Peninsula. Evaluation for wheat and Holm oak, Environ. Pollut., 165, 25–37,, 2012. a

Demetillo, M. A. G., Anderson, J. F., Geddes, J. A., Yang, X., Najacht, E. Y., Herrera, S. A., Kabasares, K. M., Kotsakis, A. E., Lerdau, M. T., and Pusede, S. E.: Observing Severe Drought Influences on Ozone Air Pollution in California, Environ. Sci. Technol., 53, 4695–4706,, 2019. a, b

Derognat, C.: Effect of biogenic volatile organic compound emissions on tropospheric chemistry during the Atmospheric Pollution Over the Paris Area (ESQUIF) campaign in the Ile-de-France region, J. Geophys. Res., 108, 8560,, 2003. a

de Rosnay, P., Polcher, J., Bruen, M., and Laval, K.: Impact of a physically based soil water flow and soil-plant interaction representation for modeling large-scale land surface processes, J. Geophys. Res.-Atmos., 107, ACL 3-1–ACL 3-19,, 2002. a

De Smedt, I., Stavrakou, T., Hendrick, F., Danckaert, T., Vlemmix, T., Pinardi, G., Theys, N., Lerot, C., Gielen, C., Vigouroux, C., Hermans, C., Fayt, C., Veefkind, P., Müller, J.-F., and Van Roozendael, M.: Diurnal, seasonal and long-term variations of global formaldehyde columns inferred from combined OMI and GOME-2 observations, Atmos. Chem. Phys., 15, 12519–12545,, 2015. a

EBAS: Station FR0033R – Ersa (Cape Corsica) – steel_canister – FR07L_SilcoCan+Silonite – isoprene – 2012-06-05 – 2014-06-27, NILU [data set],, last access: 14 April 2022. a

EEA: Air quality in Europe – 2020 report, European Environment Agency, (last access: 14 April 2022), 2020. a

EEA: Air Quality e-Reporting (AQ e-Reporting), European Environment Agency [data set],, last access: 14 April 2022. a

Ek, M. B., Mitchell, K. E., Lin, Y., Rogers, E., Grunmann, P., Koren, V., Gayno, G., and Tarpley, J. D.: Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model, J. Geophys. Res.-Atmos., 108, 2002JD003296,, 2003. a

Emberson, L., Simpson, D., Tuovinen, J.-P., Ashmore, M., and Cambridge, H.: Towards a model of ozone deposition and stomatal uptake over Europe, Research Note No. 42 – ISSN 0332-9879, Norwegian Meteorological Institute, emep/msc-w edn., 2000. a, b, c

Emmerson, K. M., Palmer, P. I., Thatcher, M., Haverd, V., and Guenther, A. B.: Sensitivity of isoprene emissions to drought over south-eastern Australia: Integrating models and satellite observations of soil moisture, Atmos. Environ., 209, 112–124,, 2019. a

European Fluxes Database Cluster: Station IT-CP2 – Castel Porziano (Lazio) – Eddy Covariance – O3 flux, European Fluxes Database Cluster [data set],, last access: 16 October 2022. a

Fares, S., Matteucci, G., Scarascia Mugnozza, G., Morani, A., Calfapietra, C., Salvatori, E., Fusaro, L., Manes, F., and Loreto, F.: Testing of models of stomatal ozone fluxes with field measurements in a mixed Mediterranean forest, Atmos. Environ., 67, 242–251,, 2013. a

Feng, Z., Yuan, X., Fares, S., Loreto, F., Li, P., Hoshika, Y., and Paoletti, E.: Isoprene is more affected by climate drivers than onoterpenes: A meta-analytic review on plant isoprenoid emissions, Plant Cell Environ., 42, 1939–1949,, 2019. a

Friedl, M. A., Sulla-Menashe, D., Tan, B., Schneider, A., Ramankutty, N., Sibley, A., and Huang, X.: MODIS Collection 5 global land cover: Algorithm refinements and characterization of new datasets, Remote Sens. Environ., 114, 168–182,, 2010. a

García-Vila, M., Lorite, I., Soriano, M., and Fereres, E.: Management trends and responses to water scarcity in an irrigation scheme of Southern Spain, Agr. Water Manage., 95, 458–468,, 2008. a

Gaubert, B., Coman, A., Foret, G., Meleux, F., Ung, A., Rouil, L., Ionescu, A., Candau, Y., and Beekmann, M.: Regional scale ozone data assimilation using an ensemble Kalman filter and the CHIMERE chemical transport model, Geosci. Model Dev., 7, 283–302,, 2014. a

Gerosa, G., Finco, A., Mereu, S., Vitale, M., Manes, F., and Denti, A. B.: Comparison of seasonal variations of ozone exposure and fluxes in a Mediterranean Holm oak forest between the exceptionally dry 2003 and the following year, Environ. Pollut., 157, 1737–1744,, 2009. a

Gong, C., Lei, Y., Ma, Y., Yue, X., and Liao, H.: Ozone–vegetation feedback through dry deposition and isoprene emissions in a global chemistry–carbon–climate model, Atmos. Chem. Phys., 20, 3841–3857,, 2020. a

González Abad, G., Liu, X., Chance, K., Wang, H., Kurosu, T. P., and Suleiman, R.: Updated Smithsonian Astrophysical Observatory Ozone Monitoring Instrument (SAO OMI) formaldehyde retrieval, Atmos. Meas. Tech., 8, 19–32,, 2015. a

Greve, P., Warrach-Sagi, K., and Wulfmeyer, V.: Evaluating Soil Water Content in a WRF-Noah Downscaling Experiment, J. Appl. Meteorol. Clim., 52, 2312–2327,, 2013. a

Guenther, A., Karl, T., Harley, P., Wiedinmyer, C., Palmer, P. I., and Geron, C.: Estimates of global terrestrial isoprene emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature), Atmos. Chem. Phys., 6, 3181–3210,, 2006. a, b

Guenther, A. B., Zimmerman, P. R., Harley, P. C., Monson, R. K., and Fall, R.: Isoprene and monoterpene emission rate variability: Model evaluations and sensitivity analyses, J. Geophys. Res., 98, 12609,, 1993. a

Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492,, 2012. a, b, c

Guion, A.: Agricultural droughts and heatwaves in the Euro-Mediterranean, daily indicator over the 1979–2016 period based on a coupled land surface and weather model (ORCHIDEE-WRF), IPSL Data Catalog [data set],, last access: 14 April 2022. a, b

Guion, A., Turquety, S., Polcher, J., Pennel, R., Bastin, S., and Arsouze, T.: Droughts and heatwaves in the Western Mediterranean: impact on vegetation and wildfires using the coupled WRF-ORCHIDEE regional model (RegIPSL), Clim. Dynam., 58, 2881–2903,, 2021.​​​​​​​ a, b, c, d, e, f

Hauglustaine, D. A., Balkanski, Y., and Schulz, M.: A global model simulation of present and future nitrate aerosols and their direct radiative forcing of climate, Atmos. Chem. Phys., 14, 11031–11063,, 2014. a

Haworth, M., Catola, S., Marino, G., Brunetti, C., Michelozzi, M., Riggi, E., Avola, G., Cosentino, S. L., Loreto, F., and Centritto, M.: Moderate Drought Stress Induces Increased Foliar Dimethylsulphoniopropionate (DMSP) Concentration and Isoprene Emission in Two Contrasting Ecotypes of Arundo donax, Front. Plant Sci., 8, 1016,, 2017. a

Hodnebrog, Ø., Solberg, S., Stordal, F., Svendby, T. M., Simpson, D., Gauss, M., Hilboll, A., Pfister, G. G., Turquety, S., Richter, A., Burrows, J. P., and Denier van der Gon, H. A. C.: Impact of forest fires, biogenic emissions and high temperatures on the elevated Eastern Mediterranean ozone levels during the hot summer of 2007, Atmos. Chem. Phys., 12, 8727–8750,, 2012. a

Huang, M., Crawford, J. H., Carmichael, G. R., Bowman, K. W., Kumar, S. V., and Sweeney, C.: Satellite soil moisture data assimilation impacts on modeling weather variables and ozone in the southeastern US – Part 2: Sensitivity to dry-deposition parameterizations, Atmos. Chem. Phys., 22, 7461–7487,, 2022. a, b

Jacob, D. J.: Introduction to atmospheric chemistry, Princeton university press, ISBN 9780691001852, 1999. a

Jaén, C., Udina, M., and Bech, J.: Analysis of two heat wave driven ozone episodes in Barcelona and surrounding region: Meteorological and photochemical modeling, Atmos. Environ., 246, 118037,, 2021. a, b

Jerrett, M., Burnett, R. T., Pope, C. A., Ito, K., Thurston, G., Krewski, D., Shi, Y., Calle, E., and Thun, M.: Long-Term Ozone Exposure and Mortality, New Engl. J. Med., 360, 1085–1095,, 2009. a

Jiang, X., Guenther, A., Potosnak, M., Geron, C., Seco, R., Karl, T., Kim, S., Gu, L., and Pallardy, S.: Isoprene emission response to drought and the impact on global atmospheric chemistry, Atmos. Environ., 183, 69–83,, 2018. a, b

Kavassalis, S. C. and Murphy, J. G.: Understanding ozone-meteorology correlations: A role for dry deposition: Ozone-Meteorology Correlations: Dry Dep, Geophys. Res. Lett., 44, 2922–2931,, 2017. a

Lamsal, L. N., Martin, R. V., van Donkelaar, A., Steinbacher, M., Celarier, E. A., Bucsela, E., Dunlea, E. J., and Pinto, J. P.: Ground-level nitrogen dioxide concentrations inferred from the satellite-borne Ozone Monitoring Instrument, J. Geophys. Res., 113, D16308,, 2008. a, b

Levelt, P., van den Oord, G., Dobber, M., Malkki, A., Huib Visser, Johan de Vries, Stammes, P., Lundell, J., and Saari, H.: The ozone monitoring instrument, IEEE T. Geosci. Remote, 44, 1093–1101,, 2006. a, b

Levelt, P. F., Joiner, J., Tamminen, J., Veefkind, J. P., Bhartia, P. K., Stein Zweers, D. C., Duncan, B. N., Streets, D. G., Eskes, H., van der A, R., McLinden, C., Fioletov, V., Carn, S., de Laat, J., DeLand, M., Marchenko, S., McPeters, R., Ziemke, J., Fu, D., Liu, X., Pickering, K., Apituley, A., González Abad, G., Arola, A., Boersma, F., Chan Miller, C., Chance, K., de Graaf, M., Hakkarainen, J., Hassinen, S., Ialongo, I., Kleipool, Q., Krotkov, N., Li, C., Lamsal, L., Newman, P., Nowlan, C., Suleiman, R., Tilstra, L. G., Torres, O., Wang, H., and Wargan, K.: The Ozone Monitoring Instrument: overview of 14 years in space, Atmos. Chem. Phys., 18, 5699–5745,, 2018. a

Lhotka, O. and Kyselý, J.: Characterizing joint effects of spatial extent, temperature magnitude and duration of heat waves and cold spells over Central Europe, Int. J. Climatol., 35, 1232–1244,, 2015. a

Lin, M., Malyshev, S., Shevliakova, E., Paulot, F., Horowitz, L. W., Fares, S., Mikkelsen, T. N., and Zhang, L.: Sensitivity of Ozone Dry Deposition to Ecosystem-Atmosphere Interactions: A Critical Appraisal of Observations and Simulations, Global Biogeochem. Cy., 33, 1264–1288,, 2019. a

Lin, M., Horowitz, L., Xie, Y., Paulot, F., Malyshev, S., Shevliakova, E., Finco, A., Gerosa, G., Kubistin, D., and Pilegaard, K.: Vegetation feedbacks during drought exacerbate ozone air pollution extremes in Europe, Nat. Clim. Change, 10, 444–451,, 2020. a, b

Liss, P. S. and Mervilat, L.: Air-Sea Gas Exchange Rates: Introduction and Synthesis, Springer Netherlands, Dordrecht, 113–127,, 1986. a

Maignan, F., Bréon, F.-M., Chevallier, F., Viovy, N., Ciais, P., Garrec, C., Trules, J., and Mancip, M.: Evaluation of a Global Vegetation Model using time series of satellite vegetation indices, Geosci. Model Dev., 4, 1103–1114,, 2011. a

Marticorena, B. and Bergametti, G.: Modeling the atmospheric dust cycle: 1. Design of a soil-derived dust emission scheme, J. Geophys. Res., 100, 16415,, 1995. a

Menut, L., Bessagnet, B., Briant, R., Cholakian, A., Couvidat, F., Mailler, S., Pennel, R., Siour, G., Tuccella, P., Turquety, S., and Valari, M.: The CHIMERE v2020r1 online chemistry-transport model, Geosci. Model Dev., 14, 6781–6811,, 2021. a, b, c, d, e

Mertens, M., Kerkweg, A., Grewe, V., Jöckel, P., and Sausen, R.: Attributing ozone and its precursors to land transport emissions in Europe and Germany, Atmos. Chem. Phys., 20, 7843–7873,, 2020. a, b

Michou, M., Laville, P., Serça, D., Fotiadi, A., Bouchou, P., and Peuch, V.-H.: Measured and modeled dry deposition velocities over the ESCOMPTE area, Atmos. Res., 74, 89–116,, 2005. a

Millet, D. B., Jacob, D. J., Turquety, S., Hudman, R. C., Wu, S., Fried, A., Walega, J., Heikes, B. G., Blake, D. R., Singh, H. B., Anderson, B. E., and Clarke, A. D.: Formaldehyde distribution over North America: Implications for satellite retrievals of formaldehyde columns and isoprene emission, J. Geophys. Res., 111, D24S02,, 2006. a

Millet, D. B., Jacob, D. J., Boersma, K. F., Fu, T.-M., Kurosu, T. P., Chance, K., Heald, C. L., and Guenther, A.: Spatial distribution of isoprene emissions from North America derived from formaldehyde column measurements by the OMI satellite sensor, J. Geophys. Res., 113, D02307,, 2008. a

Miralles, D. G., Gentine, P., Seneviratne, S. I., and Teuling, A. J.: Land-atmospheric feedbacks during droughts and heatwaves: state of the science and current challenges, Ann. NY Acad. Sci., 1436, 19–35,, 2019. a

MODIS: MODIS Leaf Area Index/FPAR, NASA [data set],, last access: 14 April 2022. a

Monahan, E.: In The Role of Air-Sea Exchange in Geochemical Cyclong, chap. The ocean as a source of atmospheric particles, Kluwer Academic Publishers, Dordrecht, Holland, 129–163, ISBN 978-94-010-8606-6, 1986. a

Müller, J.-F., Stavrakou, T., Wallens, S., De Smedt, I., Van Roozendael, M., Potosnak, M. J., Rinne, J., Munger, B., Goldstein, A., and Guenther, A. B.: Global isoprene emissions estimated using MEGAN, ECMWF analyses and a detailed canopy environment model, Atmos. Chem. Phys., 8, 1329–1341,, 2008. a

Myneni, R., Knyazikhin, Y., and Park, T.: MOD15A2H MODIS/terra leaf area index/FPAR 8-day L4 global 500 m SIN grid V006, NASA EOSDIS Land Processes DAAC, (last access: 14 April 2022), 2015. a

NCEP-NCAR: NCEP-NCAR Reanalysis 1, Physical Sciences Laboratory NOOA [data set],, last access: 14 April 2022. a

Nuvolone, D., Petri, D., and Voller, F.: The effects of ozone on human heatlh, Environ. Sci. Pollut. Res., 25, 8074–8088,,2018. a

Opacka, B., Müller, J.-F., Stavrakou, T., Bauwens, M., Sindelarova, K., Markova, J., and Guenther, A. B.: Global and regional impacts of land cover changes on isoprene emissions derived from spaceborne data and the MEGAN model, Atmos. Chem. Phys., 21, 8413–8436,, 2021. a

Otero, N., Jurado, O. E., Butler, T., and Rust, H. W.: The impact of atmospheric blocking on the compounding effect of ozone pollution and temperature: a copula-based approach, Atmos. Chem. Phys., 22, 1905–1919,, 2022. a

Panthou, G., Vrac, M., Drobinski, P., Bastin, S., and Li, L.: Impact of model resolution and Mediterranean sea coupling on hydrometeorological extremes in RCMs in the frame of HyMeX and MED-CORDEX, Clim. Dynam., 51, 915–932,, 2018. a

Perkins-Kirkpatrick, S. E. and Gibson, P. B.: Changes in regional heatwave characteristics as a function of increasing global temperature, Sci. Rep., 7, 12256,, 2017. a

Peron, A., Kaser, L., Fitzky, A. C., Graus, M., Halbwirth, H., Greiner, J., Wohlfahrt, G., Rewald, B., Sandén, H., and Karl, T.: Combined effects of ozone and drought stress on the emission of biogenic volatile organic compounds from Quercus robur L., Biogeosciences, 18, 535–556,, 2021. a

Porter, W. C. and Heald, C. L.: The mechanisms and meteorological drivers of the summertime ozone–temperature relationship, Atmos. Chem. Phys., 19, 13367–13381,, 2019. a

Potosnak, M. J., LeStourgeon, L., Pallardy, S. G., Hosman, K. P., Gu, L., Karl, T., Geron, C., and Guenther, A. B.: Observed and modeled ecosystem isoprene fluxes from an oak-dominated temperate forest and the influence of drought stress, Atmos. Environ., 84, 314–322,, 2014. a

Richards, N. A. D., Arnold, S. R., Chipperfield, M. P., Miles, G., Rap, A., Siddans, R., Monks, S. A., and Hollaway, M. J.: The Mediterranean summertime ozone maximum: global emission sensitivities and radiative impacts, Atmos. Chem. Phys., 13, 2331–2345,, 2013. a, b

Romer, P. S., Duffey, K. C., Wooldridge, P. J., Edgerton, E., Baumann, K., Feiner, P. A., Miller, D. O., Brune, W. H., Koss, A. R., de Gouw, J. A., Misztal, P. K., Goldstein, A. H., and Cohen, R. C.: Effects of temperature-dependent NOx emissions on continental ozone production, Atmos. Chem. Phys., 18, 2601–2614,, 2018. a

Russo, S., Sillmann, J., and Fischer, E. M.: Top ten European heatwaves since 1950 and their occurrence in the coming decades, Environ. Res. Lett., 10, 124003,, 2015. a

Ruti, P. M., Somot, S., Giorgi, F., Dubois, C., Flaounas, E., Obermann, A., Dell'Aquila, A., Pisacane, G., Harzallah, A., Lombardi, E., Ahrens, B., Akhtar, N., Alias, A., Arsouze, T., Aznar, R., Bastin, S., Bartholy, J., Béranger, K., Beuvier, J., Bouffies-Cloché, S., Brauch, J., Cabos, W., Calmanti, S., Calvet, J.-C., Carillo, A., Conte, D., Coppola, E., Djurdjevic, V., Drobinski, P., Elizalde-Arellano, A., Gaertner, M., Galàn, P., Gallardo, C., Gualdi, S., Goncalves, M., Jorba, O., Jordà, G., L'Heveder, B., Lebeaupin-Brossier, C., Li, L., Liguori, G., Lionello, P., Maciàs, D., Nabat, P., Önol, B., Raikovic, B., Ramage, K., Sevault, F., Sannino, G., Struglia, M. V., Sanna, A., Torma, C., and Vervatis, V.: Med-CORDEX Initiative for Mediterranean Climate Studies, B. Am. Meteorol. Soc., 97, 1187–1208,, 2016. a

Saunier, A., Ormeño, E., Boissard, C., Wortham, H., Temime-Roussel, B., Lecareux, C., Armengaud, A., and Fernandez, C.: Effect of mid-term drought on Quercus pubescens BVOCs' emission seasonality and their dependency on light and/or temperature, Atmos. Chem. Phys., 17, 7555–7566,, 2017. a

Sha, T., Ma, X., Zhang, H., Janechek, N., Wang, Y., Wang, Y., Castro García, L., Jenerette, G. D., and Wang, J.: Impacts of Soil NOx Emission on O3 Air Quality in Rural California, Environ. Sci. Technol., 55, 7113–7122,, 2021. a

Simpson, D., Fagerli, H., Jonson, J. E., Tsyro, S., Wind, P., and Tuovien, J. P.: Transboundary Acidification, Eutrophication and Ground Level Ozone in Europe. Part I: Unified EMEP Model Description, Technical report, Norwegian Meteorological Institute, ISSN 1504-6192, 2003. a

Simpson, D., Benedictow, A., Berge, H., Bergström, R., Emberson, L. D., Fagerli, H., Flechard, C. R., Hayman, G. D., Gauss, M., Jonson, J. E., Jenkin, M. E., Nyíri, A., Richter, C., Semeena, V. S., Tsyro, S., Tuovinen, J.-P., Valdebenito, Á., and Wind, P.: The EMEP MSC-W chemical transport model – technical description, Atmos. Chem. Phys., 12, 7825–7865,, 2012. a, b

Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Barker, D. M., Wang, W., and Powers, J. G.: A description of the Advanced Research WRF version 3. NCAR/TN-475+STR, Tech. rep., University Corporation for Atmospheric Research,, 2008. a, b

Spinoni, J., Naumann, G., Vogt, J. V., and Barbosa, P.: The biggest drought events in Europe from 1950 to 2012, J. Hydrol., 3, 509–524,, 2015. a

Spinoni, J., Vogt, J. V., Naumann, G., Barbosa, P., and Dosio, A.: Will drought events become more frequent and severe in Europe?, Int. J. Climatol., 38, 1718–1736,, 2018. a

Stéfanon, M., D'Andrea, F., and Drobinski, P.: Heatwave classification over Europe and the Mediterranean region, Environ. Res. Lett., 7, 014023,, 2012. a

Steinkamp, J. and Lawrence, M. G.: Improvement and evaluation of simulated global biogenic soil NO emissions in an AC-GCM, Atmos. Chem. Phys., 11, 6063–6082,, 2011. a

Stohl, A., Williams, E., Wotawa, G., and Kromp-Kolb, H.: A European inventory of soil nitric oxide emissions and the effect of these emissions on the photochemical formation of ozone, Atmos. Environ., 30, 3741–3755,, 1996. a

Sun, S., Tai, A. P. K., Yung, D. H. Y., Wong, A. Y. H., Ducker, J. A., and Holmes, C. D.: Influence of plant ecophysiology on ozone dry deposition: comparing between multiplicative and photosynthesis-based dry deposition schemes and their responses to rising CO2 level, Biogeosciences, 19, 1753–1776,, 2022. a, b

Svoboda, M. D. and Fuchs, B. A.: Handbook of drought indicators and indices., Geneva, World Meteorological Organization and global water partnership edn., oCLC: 958433349, ISBN 978-92-63-11173-9, 2016. a

Tai, A. P. K., Martin, M. V., and Heald, C. L.: Threat to future global food security from climate change and ozone air pollution, Nat. Clim. Change, 4, 817–821,, 2014. a

Thompson, G. and Eidhammer, T.: A Study of Aerosol Impacts on Clouds and Precipitation Development in a Large Winter Cyclone, J. Atmos. Sci., 71, 3636–3658,, 2014. a

Thunis, P., Pernigotti, D., and Gerboles, M.: Model quality objectives based on measurement uncertainty. Part I: Ozone, Atmos. Environ., 79, 861–868,, 2013. a

Turquety, S., Menut, L., Siour, G., Mailler, S., Hadji-Lazaro, J., George, M., Clerbaux, C., Hurtmans, D., and Coheur, P.-F.: APIFLAME v2.0 biomass burning emissions model: impact of refined input parameters on atmospheric concentration in Portugal in summer 2016, Geosci. Model Dev., 13, 2981–3009,, 2020. a

Van Leer, B.: Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection, J. Comput. Phys., 23, 276–299,, 1977. a

Vautard, R., Honore, C., Beekmann, M., and Rouil, L.: Simulation of ozone during the August 2003 heat wave and emission control scenarios, Atmos. Environ., 39, 2957–2967,, 2005. a

Vicente-Serrano, S. M.: Evaluating the impact of drought using remote sensing in a Mediterranean, semi-arid region, Nat. Hazards, 40, 173–208,, 2007. a, b

Vicente-Serrano, S. M., Beguería, S., and López-Moreno, J. I.: A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index, J. Climate, 23, 1696–1718,, 2010. a

Vicente-Serrano, S. M., Gouveia, C., Camarero, J. J., Beguería, S., Trigo, R., López-Moreno, J. I., Azorín-Molina, C., Pasho, E., Lorenzo-Lacruz, J., Revuelto, J., Morán-Tejeda, E., and Sanchez-Lorenzo, A.: Response of vegetation to drought time-scales across global land biomes, P. Natl. Acad. Sci. USA, 110, 52–57,, 2013. a

Wang, P., Liu, Y., Dai, J., Fu, X., Wang, X., Guenther, A., and Wang, T.: Isoprene Emissions Response to Drought and the Impacts on Ozone and SOA in China, J. Geophys. Res.-Atmos., 126, e2020JD033263,, 2021. a

Wang, Y., Xie, Y., Dong, W., Ming, Y., Wang, J., and Shen, L.: Adverse effects of increasing drought on air quality via natural processes, Atmos. Chem. Phys., 17, 12827–12843,, 2017. a, b, c, d

Wang, Y., Tan, X., Huang, L., Wang, Q., Li, H., Zhang, H., Zhang, K., Liu, Z., Traore, D., Yaluk, E., Fu, J. S., and Li, L.: The impact of biogenic emissions on ozone formation in the Yangtze River Delta region based on MEGANv3.1, Air Qual. Atmos. Hlth., 14, 763–774,, 2021. a

Weng, H., Lin, J., Martin, R., Millet, D. B., Jaeglé, L., Ridley, D., Keller, C., Li, C., Du, M., and Meng, J.: Global high-resolution emissions of soil NOx, sea salt aerosols, and biogenic volatile organic compounds, Scientific Data, 7, 148,, 2020.​​​​​​​ a

Wesely, M.: Parameterization of surface resistances to gaseous dry deposition in regional-scale numerical models, Atmos. Environ., 41, 52–63,, 1989. a, b

Wong, A. Y. H., Geddes, J. A., Ducker, J. A., Holmes, C. D., Fares, S., Goldstein, A. H., Mammarella, I., and Munger, J. W.: New Evidence for the Importance of Non-Stomatal Pathways in Ozone Deposition During Extreme Heat and Dry Anomalies, Geophys. Res. Lett., 49, e2021GL095717,, 2022.  a

Yuan, H., Dai, Y., Xiao, Z., Ji, D., and Shangguan, W.: Reprocessing the MODIS Leaf Area Index products for land surface and climate modelling, Remote Sens. Environ., 115, 1171–1187,, 2011. a

Zampieri, M., D'Andrea, F., Vautard, R., Ciais, P., de Noblet-Ducoudré, N., and Yiou, P.: Hot European summers and the role of soil moisture in the propagation of Mediterranean drought, J. Climate, 22, 4747–4758,, 2009. a

Zhu, L., Jacob, D. J., Kim, P. S., Fisher, J. A., Yu, K., Travis, K. R., Mickley, L. J., Yantosca, R. M., Sulprizio, M. P., De Smedt, I., González Abad, G., Chance, K., Li, C., Ferrare, R., Fried, A., Hair, J. W., Hanisco, T. F., Richter, D., Jo Scarino, A., Walega, J., Weibring, P., and Wolfe, G. M.: Observing atmospheric formaldehyde (HCHO) from space: validation and intercomparison of six retrievals from four satellites (OMI, GOME2A, GOME2B, OMPS) with SEAC4RS aircraft observations over the southeast US, Atmos. Chem. Phys., 16, 13477–13490,, 2016. a, b, c

Zhu, L., Mickley, L. J., Jacob, D. J., Marais, E. A., Sheng, J., Hu, L., Abad, G. G., and Chance, K.: Long-term (2005–2014) trends in formaldehyde (HCHO) columns across North America as seen by the OMI satellite instrument: Evidence of changing emissions of volatile organic compounds: HCHO Trend Across North America, Geophys. Res. Lett., 44, 7079–7086,, 2017. a

Short summary
At high concentrations, tropospheric ozone (O3) deteriorates air quality. Weather conditions are key to understanding the variability in O3 concentration, especially during extremes. We suggest that identifying the presence of combined heatwaves is essential to the study of droughts in canopy–troposphere interactions and O3 concentration. Even so, they are associated, on average, with an increase in O3, partly explained by an increase in precursor emissions and a decrease in dry deposition.
Final-revised paper