Dry deposition to vegetation is a major sink of ground-level ozone and is responsible for about 20 % of the total tropospheric ozone loss. Its parameterization in atmospheric chemistry models represents a significant source of uncertainty for the global tropospheric ozone budget and might account for the mismatch with observations. The model used in this study, the Modular Earth Submodel System version 2 (MESSy2) linked to the fifth-generation European Centre Hamburg general circulation model (ECHAM5) as an atmospheric circulation model (EMAC), is no exception.
Like many global models, EMAC employs a “resistance in series” scheme with the major surface deposition via plant stomata which is hardly sensitive to meteorology, depending only on solar radiation. Unlike many global models, however, EMAC uses a simplified high resistance for non-stomatal deposition which makes this pathway negligible in the model. However, several studies have shown this process to be comparable in magnitude to the stomatal uptake, especially during the night over moist surfaces.
Hence, we present here a revised dry deposition in EMAC including meteorological adjustment factors for stomatal closure and an explicit cuticular pathway.
These modifications for the three stomatal stress functions have been included in the newly developed MESSy VERTEX submodel, i.e. a process model describing the vertical exchange in the atmospheric boundary layer, which will be evaluated for the first time here.
The scheme is limited by a small number of different surface types and generalized parameters.
The MESSy submodel describing the dry deposition of trace gases and aerosols (DDEP) has been revised accordingly.
The comparison of the simulation results with measurement data at four sites shows that the new scheme enables a more realistic representation of dry deposition. However, the representation is strongly limited by the local meteorology.
In total, the changes increase the dry deposition velocity of ozone up to a factor of 2 globally, whereby the highest impact arises from the inclusion of cuticular uptake, especially over moist surfaces.
This corresponds to a 6 % increase of global annual dry deposition loss of ozone resulting globally in a slight decrease of ground-level ozone but a regional decrease of up to 25 %. The change of ozone dry deposition is also reasoned by the altered loss of ozone precursors. Thus, the revision of the process parameterization as documented here has, among others, the potential to significantly reduce the overestimation of tropospheric ozone in global models.
Introduction
Ground-level ozone is a secondary air pollutant which is harmful for humans and ecosystems. Besides chemical destruction, a large fraction of it is removed by dry deposition which accounts for about 20 % of the total O3 loss .
The process description of dry deposition considers boundary-layer meteorology (e.g. turbulence), chemical properties of the trace gases and surface types. In most global models, dry deposition of trace gases is parameterized using the “resistance in series” analogy by . The largest deposition rates of ozone occur over dense vegetation where it mainly follows two pathways: through leaf openings (stomata) and to leaf waxes (cuticle) .
Thereby, stomatal uptake is commonly parameterized following the empirical multiplicative approach by which uses a predefined minimum resistance and multiple environmental response factors like in , and . More advanced formulations often used by land surface models are based on the CO2 assimilation by plants during photosynthesis . Both approaches rely on the choice and constraints of ecosystem-dependent parameters and have different advantages . A further role in coupling stomata to ecosystems is played by stomatal optimization models, whereas optimal stomatal activity with a maximum amount of carbon gain and a minimum loss of water is calculated based on ecophysiological processes e.g.. Of particular interest are stomatal optimization models which, based on ecophysiological processes, maximize carbon gain while minimizing water loss. According to , these models are promising in representing stomatal behaviour and improving carbon cycle modelling. Non-stomatal deposition has been less investigated by now; therefore, most models use predefined constant resistances or scale it with leaf area index e.g., while some apply an explicit parameterization based on the observational findings of enhance cuticular uptake under leaf surface wetness .
The different parameterizations of the (surface) resistances cause main model uncertainties in computing dry deposition fluxes of trace gases, which depend on the response to hydroclimate and land-type-specific properties . Thereby, it has been shown that the original Wesely-based parameterization generally captures well the seasonal and diurnal cycles of dry deposition velocity, whereas model–observation discrepancy at seasonal scales arises from biased land type and leaf area index input data . stated that discrepancies of up to 8 ppb in ground-level ozone arise from different parameterizations.
The current dry deposition scheme of EMAC uses six surface types, where the parameterized processes represent the forest canopy as a whole (big-leaf approach). Thereby, the uptake over vegetation relies on stomatal deposition as the only pathway determined by the photosynthetically active radiation .
According to and , the stomatal uptake in parameterizations often lacks the dependence on meteorological and environmental variables (leaf area index, temperature, vapour pressure deficit). Moreover, several studies (e.g. ) found the contribution of an additional process to dry deposition at the leaf covering of plants. firstly derived a parameterization from field studies which establishes the important link of this process to meteorology.
In general, findings by , and highlight the importance of considering the dry deposition–meteorology dependence in global models. Such an extension would realistically enhance the sensitivity of dry deposition to climate variability and would result in a more accurate prediction of ground-level ozone.
Given the importance of ozone as a major tropospheric oxidant, air pollutant and greenhouse gas, an accurate representation of dry deposition is desirable .
Additionally, the significance of a realistic representation of land–atmosphere feedbacks rises in light of the changing Earth's climate with the projected increase of extreme events' frequency and intensity .
Here, we present a revision of the existing Wesely-based dry deposition scheme in the Modular Earth Submodel System (MESSy), which has a very simplified representation of vegetation and soil. The modifications are done by well-established findings about the controls of stomatal and cuticular uptake of trace gases. The calculation of stomatal deposition fluxes is extended by including the vegetation density, two meteorological adjustment factors and an improved soil moisture availability function for plant stomata following the multiplicative algorithm by . For the first time in MESSy, a parameterization for cuticular dry deposition dependent on important meteorological and environmental variables is implemented explicitly .
In Sect. , a description of the model setup and the simulations is provided, whereas especially the transition to the new vertical exchange scheme is described in detail. Subsequently, the new VERTEX scheme is evaluated. In Sect. , the impact of the changes on ozone dry deposition is evaluated on daily and seasonal scales by comparison with measurements at four different sites. Here, advantages, uncertainties and missing processes in the revised scheme are identified. Next, the global impact on ground-level ozone is assessed by separating the effect of the different implemented parameterizations. Then, Sect. provides a description of the uncertainties in modelling stomatal conductance and Sect. comprises an investigation of the sensitivity to model resolution.
Section summarizes the main findings and the remaining process and model uncertainties which form the basis for the provided recommendations. Section describes planned future developments.
Model description
This study uses the ECHAM/MESSy atmospheric chemistry model.
MESSy v2.54 provides a flexible infrastructure for coupling processes to build comprehensive Earth system models (ESMs) and is utilized here with the fifth-generation European Centre Hamburg general circulation model (ECHAM5; ) as an atmospheric general circulation model. The dry deposition process of gases is calculated within the submodel DDEP . This is described in Sect. . It relies on the VERTEX vertical exchange submodel (Sect. ), former E5VDIFF, which contains the calculation of stomatal uptake (Eq. ) and soil moisture stress (Eq. ). The stomatal uptake parameterization is the base for the evapotranspiration scheme in VERTEX (Appendix ) which also incorporates the soil moisture stress.
The new VERTEX vertical exchange submodel
The VERTEX submodel represents land–atmosphere exchange and vertical diffusion as an alternative to the default E5VDIFF submodel in ECHAM5/MESSy. In 2016, Huug Ouwersloot branched VERTEX off from E5VDIFF. He optimized the code and applied bug fixes. This includes changes in calculation of the transfer coefficients for vertical diffusion, the latent heat vaporization, the convective transfer coefficient, the storage of the friction velocity, the roughness length over sea, the kinematic heat and moisture fluxes and the 2 and 10 m friction velocity.
A detailed description can be found in the Supplement.
Dry deposition over vegetation
Dry deposition of trace gases to vegetation is calculated according to the multiple resistance scheme by shown in Fig. . The scheme, originally designed for a regional model with 11 land types and five seasonal categories, is used here with six generalized land types . This was adapted by to the surface scheme of the ECHAM climate model . The vegetation canopy is represented as one system; i.e. the detailed structure and plant characteristics are neglected (one big-leaf approach). Only one assumption about the canopy structure is made: the leaves are horizontally oriented and the leaf density is uniformly vertically distributed . This is required in the formula for the calculation of stomatal resistance (Eq. ).
Dry deposition resistance analogy (adapted from ); modified resistors are marked with red boxes.
The resistances (in sm-1) in the big-leaf approach account for mass and energy transfer mainly exerted by the boundary layer turbulence (Ra), molecular diffusion via the quasi-laminar boundary layer (Rqbr) and heterogeneous losses at the surface (Rs) . With these, the dry deposition velocity vd of a trace gas X (in sm-1) is defined as follows:
vd(X)=1Ra+Rqbr(X)+Rs(X).
The dry deposition flux fd(X) (in moleculesm-2s-1) is determined by multiplying the dry deposition velocity with the trace gas concentration C(X) (in moleculesm-3):
fd(X)=-vd(X)⋅C(X).
The total resistance over land combines the resistances over snow, soil, vegetation (veg) and wet skin (ws) weighted by the respective land-covered fraction of a grid box . In the following, only the latter two are considered.
The resistances Ra and Rqbr are commonly parameterized with standard formulations from micrometeorology .
For the surface resistance over vegetation (Rs,veg), the parameterization according to is used:
1Rs,veg(X)=1Rcan+Rs,soil(X)+1Rcut(X)+1Rstom,corr(X)+Rmes(X)︸Rleaf(X),
which consists of the soil resistance (Rs,soil(X)), the in-canopy aerodynamic resistance (Rcan) (as in ) and the leaf resistance (Rleaf(X)).
The gas uptake by leaves (leaf) can be separated in two parallel pathways: the cuticular (cut) and the stomatal (stom) with its associated mesophilic pathway (mes), where the latter has negligible resistance for ozone and highly soluble species . In contrast to the default formulation in MESSy (Eq. ), the resistances in the updated scheme are provided at canopy scale in order to avoid linear scaling with the leaf area index (LAI, area of leaves [m2]/surface area [m2]). In fact, the linear scaling of resistances with LAI assumes that the leaves act in parallel and overestimates the uptake for high LAI values (>3–4) .
Furthermore, the quasi-laminar boundary resistance of individual leaves is included through the cuticular deposition scheme (see Sect. ), whereas Rqbr,veg is a separate term in the old formulation (Eq. ).
Due to the importance of stomatal and cuticular uptake for ozone dry deposition, their respective parameterizations are modified in this study (see Sect. and ). Also, ozone deposition to soil might be an important pathway but process understanding remains limited due to scant observational constraints . showed an exponential increase of soil resistance with surface relative humidity in three agricultural data sets which, however, varies much between different sites and contradicts previous findings . Models by, e.g. apply a linear dependence on soil water content for parameterizing soil resistance. These parameterizations rely on input variables like the minimum soil resistance which introduce an uncertainty due to measurement constraints. Also, the performance of a mechanistic model as proposed by depends on many input variables and parameters whose estimation is challenging and mostly biome dependent. Due to these uncertainties and limitations, the current parameterization of soil resistance in MESSy (see for details) was not modified in this study.
Uptake through plant stomata
The stomata are actively regulated openings between the plant cells. They are scattered mostly over the lower (hypostomatous) epidermis of leaves. They control the H2O and CO2 exchange by plants which is the essential coupling of vegetation to the atmosphere and therefore to weather and climate.
Here, the default parameterization of stomatal resistance (Eq. ) is extended by adding dependencies on meteorological variables according to the Simple Biosphere Model (SiB) by based on previous work by for temperature (T) and vapour pressure deficit (VPD):
Rstom,corr(X)=Rstom(PAR,LAI)f(Ws)⋅f(T)⋅f(VPD)⋅DH2OD(X).
The optimal stomatal resistance for water (Rstom(PAR,LAI)) is corrected with the ratio of the molecular diffusivity of the species (D(X)) and water (DH20). The optimal stomatal resistance depends on the photosynthetically active radiation (PAR) and LAI :
Rstom(PAR,LAI)=kcbdPARlndexp(kLAI)+1d+1-lnd+exp(-kLAI)d+1,
where k=0.9 is the extinction coefficient, c=100sm-1 is the minimum stomatal resistance, and a=5000Jm-3, b=10Wm-2 and d=a+b⋅cc⋅PAR are fitting parameters . For historical reasons, LAI was set to 1 in order to obtain the stomatal resistance at leaf level . This has been changed and the seasonal evolution of stomatal resistance now follows the LAI which, in our study, is based on a 5-year climatology of monthly normalized difference vegetation index (NDVI) satellite data .
First, the stomatal resistance is corrected by the inverse of the temperature stress factor (1/f(T)) derived by :
6f(T)=b3(T-Tl)(Th-T)b47b3=(T0-Tl)(Th-T0)-b48b4=(Th-T0)/(Th-Tl),
where the empirical parameters are Th=318.15 K, Tl=268.15 K and T0=298.15 K.
Secondly, following the analysis by , a stress factor dependent on vapour pressure deficit (1/f(VPD)) was added to the calculation of stomatal resistance in VERTEX:
9pH2O,sat(T)=0.61078exp17.1⋅T(pH2O)235+T(pH2O)10VPD=pH2O,sat(T)-pH2O=1-RH100pH2O,sat(T)11f(VPD)=VPD-12,
with T(pH2O) (in K) as the surface temperature, pH2O (in kPa) as the pressure of water vapour and pH2O(T) [kPa] the pressure of saturated air. The vapour pressure deficit is calculated according to .
While the stomatal resistance at canopy scale is actually calculated within the MESSy VERTEX submodel, the submodel DDEP uses it for the calculation of dry deposition fluxes. Thus, in DDEP, the user can choose between the old scheme based on and the new scheme actually using the stomatal resistance at canopy scale. The latter is activated by setting the DDEP &CTRL namelist parameter l_ganzeori to .FALSE..
How the stomatal resistance is calculated is chosen in VERTEX by the &CTRL namelist parameter irstom.
irstom=0 activates the original parameterization.
Separate modifications:
irstom=2: variable LAI,
irstom=3: T dependency and
irstom=4: VPD dependency, respectively.
irstom=5: all modifications.
irstom=1: stomatal resistance with variable LAI at leaf scale. Instead of choosing LAI of 1 in Eq. () to represent the stomatal resistance at leaf level, as is done by the original code, Eq. (5) is calculated at canopy level using the actual LAI and then multiplied by LAI to obtain the average stomatal resistance at leaf level. For this case, the DDEP namelist parameter l_ganzeori has to be set to .TRUE..
The stomatal activity of plants and the strength of surface–atmosphere coupling strongly depend on the parameterized plant–water stress . The soil water budget is represented by a “bucket scheme” where the soil water in a single layer is prescribed by a geographically varying predefined field capacity and soil wetness governed by transpiration, precipitation, runoff, snow melt and drainage . This scheme is used by so-called “first-generation” models. However, EMAC controls evapotranspiration through the stomatal resistance (Appendix ), which is the most important feature of biophysical (“second-generation”) land surface models. Thereby, the stomatal resistance is calculated often like the one described here (Eq. ) including temperature, VPD and soil moisture stress . The originally used plant–water stress function of and , however, relies on leaf water potential (f(ψ)) for different plant types, which is difficult to estimate. Hence, EMAC uses a plant–water stress function dependent on soil moisture (f(Ws). The default parameterization (Eq. , ifws=0 in VERTEX &CTRL), applies the permanent wilting point of plants (Wpwp, 35 % of field capacity
maximum amount of water the soil can hold against gravity over periods of several days
) as a lower threshold in the calculation of the soil moisture stress factor (f(Ws)). However, soil moisture is significantly underpredicted by the model in some regions and the calculated f(Ws) can be 0 for long periods. This is unrealistic and effectively shuts down dry deposition, e.g. during the dry season in the Amazon region. For this reason, f(Ws) is parameterized here according to the original formulation by by removing the lower limit:
f(Ws)=1Ws(t)>WcrWs(t)WcrWs(t)≤Wcr,
where Ws(t) is the surface soil wetness (in m). Wcr (in m) is defined as the critical soil moisture level (75 % of the field capacity) at which the transpiration of plants is reduced. The modified parameterization in Eq. () can be applied by setting the &CTRL parameter ifws= in the VERTEX namelist.
Cuticular deposition
According to several field studies (e.g. ), cuticular deposition is an important contributor to ozone uptake and should not be neglected in models. Therefore, an explicit parameterization of cuticular deposition as used in many North American air quality modelling studies has been implemented.
The gas uptake by leaf surfaces is based on two parallel routes, for which an analogy to ozone (highly reactive) and sulfur dioxide (very soluble) is used. The cuticular resistance is calculated as
Rcut(X)=Rcut,d(O3)10-5⋅H(X)+sreac(X),
where H(X) is the effective Henry's law coefficient as a measure of the solubility. The reactivity of a species is rated by the parameter sreac. For highly reactive species (sreac=1), the same property as for ozone is assumed (second term in Eq. ), while for less reactive species (sreac=0.1,0) the uptake is effectively reduced .
For soluble species, the uptake at wet skin is assumed to be similar to the one of sulfur dioxide and is calculated as
Rws(X)=1/3Rcut,w(SO2)+10-7⋅H(X)+sreac(X)Rcut,w(O3)-1,
where Rcut,w(SO2) and Rcut,w(O3) are the resistances of sulfur dioxide and ozone at wet surfaces, respectively.
The constant values of the default formulae (Eqs. , ) are replaced by parameterizations which account for the meteorological dependence of cuticular uptake according to :
15Rcut,d(O3/SO2)=Rcut,d0(O3/SO2)exp(0.03⋅RH)⋅LAI0.25⋅u*16Rcut,w(O3/SO2)=Rcut,w0(O3/SO2)LAI0.5⋅u*,
where the cuticular resistance of O3 and SO2, respectively, is distinguished for dry canopies (Rcut,d) and wet canopies (Rcut,w) depending on relative humidity (RH in %), LAI (in m2m-2) and friction velocity (u* in ms-1).
The input parameters are Rcut,d0(O3)=5000 sm-1, Rcut,w0(O3)=300sm-1 and Rcut,d0(SO2)=2000sm-1. For rain and dew conditions, values of 50 sm-1 and 100 sm-1 are prescribed for Rcut,w0(SO2).
In contrast to traditional approaches, these parameterizations also consider the aerodynamic and the quasi-laminar boundary resistances of individual leaves. For the usage in MESSy, this can be switched on via l_ganzeori=.FALSE. in the &CTRL namelist of DDEP.
Simulations
In order to answer the different research questions of this study, two different types of simulations have been performed (Table ).
The first kind were simulations to investigate dry deposition and the effect of the modifications in VERTEX:these simulations are based on the Chemistry-Climate Model Initiative (CCMI) setup . To allow for comparison with measurements, the model dynamics have been nudged towards realistic meteorology by the assimilation of data from the European Centre for Medium-Range Weather Forecasts (ECMWF) . Additionally, the QCTM mode is used; i.e. the chemistry does not feed back to the dynamics, resulting in the same meteorology for all simulations . All modifications for the dry deposition scheme are employed in a 7-year simulation (REV, 2009–2015). Additionally, a 1.5-year simulation covering the period 2017 to July 2018 (2017 as spin-up) has been performed to cover the measurement periods (Sect. ). For the same periods, simulations with the same configuration, except applying the default dry deposition scheme (DEF), have been conducted. The individual effects of the different modifications are investigated by two 2-year simulations employing the different namelist switches (Sect. ). Moreover, a free-running sensitivity simulation with an additional temperature and drought stress factor for evapotranspiration (Appendix ) has been performed aiming at an improved representation of local meteorology especially in the Amazon. The station simulation output and the global output are analysed in Sects. and , respectively. In addition, two 2-year simulations are realized for different horizontal resolutions (REST42, REST63) to investigate the resolution dependency of dry deposition (Sect. ). All these simulations use 31 model layers with the top at 10 hPa and take the first year of simulation as spin-up.
The second kind were simulations for the evaluation of VERTEX as the boundary layer scheme:two pure dynamical (i.e. without chemistry) 30-year simulations with the old (clim-E5) and the new boundary layer description (clim-VER), respectively, have been performed.
All simulations were performed at the Jülich Supercomputing Centre with the JURECA Cluster .
List of EMAC simulations
SimulationSpatial resolutionTime periodRemarks(1) Dry deposition mechanism: CCMI chemistry, nudged, no feedbacks (QCTM) REST42T42L31 (2.8∘× 2.8∘)2009/2010irstom=5, ifws=1, l_ganzeori=FREST63T63L31 (1.9∘× 1.9∘)2009/2010irstom=5, ifws=1, l_ganzeori=FREV (revised)T106L31 (1.1∘× 1.1∘)2009–2015, 2017–June 2018irstom=5, ifws=1, l_ganzeori=FDEF (default)T106L31 (1.1∘× 1.1∘)2009–2015, 2017–June 2018default ddep schemeREV-fwsT106L31 (1.1∘× 1.1∘)2009/2010irstom=5, ifws=0, l_ganzeori=FREV-fTfDT106L31 (1.1∘× 1.1∘)2009/2010irstom=2, ifws=1, l_ganzeori=FREV-NNTRT106L31 (1.1∘× 1.1∘)2014/2015free-running, all ddep modifications (as REV), all stress factors applied to evapotranspiration (izwet=1)(2) Climatology comparison: no chemistry, free-running clim-E5T42L90 (2.8∘× 2.8∘,up to 0.01 hPa)1979–2008E5VDIFF for vertical exchangeclim-VERT42L90 (2.8∘× 2.8∘,up to 0.01 hPa)1979–2008VERTEX for vertical exchangeVERTEX evaluation
In order to advise the usage of VERTEX (with the default settings) as the default vertical exchange submodel in MESSy, the dynamics produced by both submodels are compared. Therefore, two dynamical, free-running 30-year simulations have been performed using the E5VDIFF or the VERTEX submodels, respectively. To obtain a comparable radiative imbalance at TOA (top of the atmosphere) with VERTEX, the four cloud parameters have been tuned in advance according to . The tuning factors can be found in Table . The radiative imbalance at TOA is slightly positive at present-day conditions ; here, E5VDIFF gives a negative value. The difference between the tuned VERTEX and E5VDIFF is small and within the uncertainty range of ±0.4 Wm-2.
Overview of tuning parameter settings and global mean properties.
ParametersEMAC (E5VDIFF)EMAC (VERTEX)Cloud mass flux above level of non-buoyancy0.30.3Entrainment rate for shallow convection1×10-31×10-3Entrainment rate for deep convection1 ×10-41 ×10-4Conversion rate to rain in convective clouds1.5 ×10-41.6 ×10-4PropertiesObserved*EMAC (E5VDIFF)EMAC (VERTEX)Total cloud cover [%]67.1267.27Water vapour path [kgm-2]25.0324.83Liquid water path [kgm-2]0.0770.077Total precipitation [mm d-1]1.281.31Surface net shortwave [Wm-2]152–167158.27158.32Surface net longwave [Wm-2]-(40–57)-54.82-54.93Surface sensible heat flux [Wm-2]-(16–19)-18.75-19.65Surface latent heat flux [Wm-2]-(75–87)-87.45-88.73Planetary albedo [%]32.3832.37Shortwave net at TOA [Wm-2]238–244230.99231.00Longwave net at TOA [Wm-2]-(237–241)-232.46-232.55Radiation imbalance at TOA [Wm-2]-1.47-1.55
*
Additionally, global mean values of surface temperature, cloud liquid water, relative humidity and planetary boundary layer height of EMAC using E5VDIFF and EMAC using VERTEX with the respective uncertainty range for the period 1979–2008 are represented in Fig. . The results for cloud liquid water and planetary boundary height show no significant differences between the VERTEX and E5VDIFF simulations since the annual mean of each falls in the confidence interval of the other. This is not always the case for surface temperature and relative humidity. However, the 30-year means of surface temperature and relative humidity simulated by E5VDIFF and VERTEX are not significantly different.
Global mean properties and the uncertainty range (95 % confidence interval in shaded) of the climatology simulations with E5VDIFF (clim-E5) and with VERTEX (clim-VER) for the period 1979–2008.
Evaluation with deposition measurements
To assess the impact of the code revision/modifications on the variability of dry deposition, we compare the sensitivity simulations DEF, REV, REV-fTfVPD, REV-fws and REV-NNTR (see Table , all at T106L31 resolution) with dry deposition measurements at four field sites (listed in Table ). The chosen data sets are the best available of ozone dry deposition (flux data and ozone mixing ratio or velocity data) with the required temporal resolution and coverage of diverse biomes of the world. The analysis is aimed at covering the recent decade, which includes the most extreme drought and heat
events (where the stomatal stress factors are aimed). For the reason of uniqueness and importance of atmospheric processes in a remote and pristine forest like the Amazon, we included measurements from, among others, the Amazon Tall Tower Observatory (ATTO). Ozone dry deposition fluxes were measured with the eddy covariance and gradient method (Ontario). From these data, deposition velocities were calculated by the means of ozone concentration data. The eddy covariance technique determines a turbulent flux by the covariance of the measured vertical velocity and the gas concentration. Due to the stochastic nature of turbulence, these measurements have an uncertainty of 10 % to 20 % under typical observation conditions . For the gradient method used at the Borden forest research station, the dry deposition flux was estimated from concentration gradients below and above the canopy and the eddy diffusivity according to the Monin–Obukhov similarity theory. The estimated dry deposition velocities (Vd) show an uncertainty of ≈ 20 %, which is due to the assigned canopy, the inherent limitations of the algorithm and the measurement uncertainties in concentrations. However, results are in good agreement with other eddy covariance measurements .
Dry deposition measurements. In the description of vegetation/climate, the reported LAI (in m2m-2) is given in brackets; vdmod and vdobs are the average measured and modelled dry deposition velocity.
SiteVegetation/climateLocation (height)Time periodvdmod(vdobs)cms-1ReferenceHyytiälä, southern Finland (SMEAR II)Boreal forest, Scots Pine, (LAI of 3–4)/cold temperate61.85∘ N, 24.28∘ E (22 m/16 ma)2010–20120.29 (0.28)Lindcove research station, California, USACitrus orchard (LAI of 3)/Mediterranean36.35∘ N, 119.09∘ W (131 m)Oct 2009–Nov 20100.22 (0.49) and bBorden research station, Ontario, CanadaMixed forest (LAI of 4.6)/temperate44.19∘ N, 79.56∘ W (33 m)2010–20120.34 (0.47)Amazon Tall Tower(ATTO), Manaus, BrazilRainforest (LAI of 6)/tropical humid2.15∘ S, 59.01∘ W (41 m)Nov 2015, Apr/May 20180.18 (0.67), 0.33 (1,0)Available on request: Matthias Sörgel (m.soergel@mpic.de)
a Meteorological measurement height. b Ozone data are not available here.
Annual cycle of dry deposition
The annual cycle of dry deposition is mainly driven by the evolution of vegetation and is generally represented well in models .
Here, we use the long time series measured at Borden and Hyytiälä to identify the impact of the code modifications on the annual cycle of dry deposition velocity. The available micrometeorological data help to distinguish the different effects. From the hourly data, we calculated multiyear (2010–2012) monthly means. To explore the contribution of stomatal and cuticular uptake, the individual velocities are calculated for O3 according to the model calculations :
Gcut,d=(1-ws)⋅(1-cvs)⋅vegRcut,d(O3)17Gcut,w=ws⋅(1-cvs)Rcut,w(O3)18Gns=Gcut,d+Gcut,w19Gstom=(1-ws)⋅(1-cvs)⋅vegRstom,corr(O3)20vp=GpGstom+Gns⋅vd,
where G names the individual conductances (inverse of resistance) of stomata (stom), dry cuticle (cut,d), wet cuticle (cut,w) and non-stomata (ns). Here, veg, ws and cvs give the vegetation fraction, the wet skin fraction and the snow-covered fraction, respectively. Gp and vp are the individual conductance and the velocity of one pathway. Further terms are described in Sect. .
The multiyear (2010–2012) annual cycle of the simulated dry deposition velocity at Borden forest (Fig. a) captures the observed cycle well until June. The new scheme reproduces the observations better than the old scheme. This is a consequence of the increase in nighttime mean velocities due to the much larger cuticular contribution (Fig. a, b). However, due to the overestimated stomatal uptake in the default scheme (see Sect. ), only slight deviations from the new dry deposition scheme are visible in the daily mean shown in Fig. a. The mismatch of the simulated and measured Vd from August to October is a consequence of the underestimation of relative humidity leading to too-low simulated cuticular deposition (Fig. c, e). This effect exceeds the impact of the overestimation of relative humidity (only) in summer, because the LAI is higher in summer. In general, the cuticular uptake parameterization accounting for LAI, friction velocity, RH and surface wetness conditions performs, in our simulations, better than parameterizations without these dependencies as expected from the study of . Unfortunately, the cuticular uptake parameterization also introduces uncertainties to the modelled non-stomatal uptake. Moreover, accounting for biogenic volatile organic compounds (BVOCs) like in would enhance in-canopy loss of ozone, significantly increase non-stomatal dry deposition and lead to improved simulation results . The representation of in-canopy air chemistry is outside the scope of the present study but planned within a subsequent study.
In contrast, the amplitude of the annual cycle and the mean of dry deposition fluxes in Hyytiälä are overestimated by both schemes during spring and summer (Fig. b). For the default scheme, this is due to the oversimplification of the stomatal uptake that only accounts for a constant LAI of 1 m2m-2 (see Sect. ), which is far from the measured LAI of 3–4 m2m-2 during this period . Enabling the new scheme (REV), increases the dry deposition velocity which reproduces the measured values in autumn better. The contribution of non-stomatal dry deposition of 25 %–45 % during the day reported by is represented partly by that. However, the new scheme leads to an even higher overestimation by the model from April to July. The sensitivity simulation REV-fws (default f(Ws)) points to the increase of the soil moisture stress function (see Sect. , Eq. ) as one reason for the overestimation of Vd in summer (Fig. b, f). Moreover, the overestimation in June/July is partly (∼ 10 %) due to the too-high model LAI compared to the measured values of 3–4 (Fig. a). The remaining gap (Fig. f) can be explained by restricting the analysis to wet conditions (RH > 70 %) only and the analysis of the sensitivity simulation REF-fTfD (no f(T) and f(VPD)). This suggests that the overestimated Vd (Fig. c) in summer is due to the stress factors for stomatal uptake since the modelled and measured temperatures are a mismatch. VPD has been identified by as a strong driver of daytime total deposition velocity, which confirms the importance of inclusion of VPD dependence for stomatal uptake.
Measured (obs) and modelled (DEF, REV) multiyear mean (2010–2012) and REV-fws (2010) annual cycle. (a, c, e) Borden forest and (b, d, f) Hyytiälä; arrows indicate 1σ.
Importance of stress factors for the diurnal variation of deposition
The short-term measurements at Lindcove research station and at ATTO are used to assess the impact of the stress factors on the diurnal cycle of dry deposition velocity in spring and summer. Additionally, micrometeorological and additional flux data make possible to consider the stomatal resistance (∼ inverse of the velocity, calculations according to ) and the underlying meteorological conditions. Since the respective micrometeorological measurements are not available at ATTO, data extracted from the ERA5 global climate reanalysis at the 1000 hPa pressure level are used here.
The diurnal cycle of dry deposition velocity at the Lindcove research station follows the solar variation and is generally well reproduced by the model with the best match in spring (Fig. ). The revised dry deposition scheme reduces the underestimation of measured nighttime Vd due to the inclusion of cuticular uptake, which identified as an important ozone sink for exactly this measurement site. The measured dry deposition velocity increases at sunrise (around 15:00 UTC) and remains almost constant during the day. This is only reproduced by the revised dry deposition scheme. The comparison of the dry deposition velocity from the revised scheme (red line) and the velocity without stomatal T and VPD stress (gray line) in Fig. a illustrates the necessity of accounting for the stress factors. This is consistent with , who reported a high negative correlation of Vd(sto) with VPD and temperature and related it to stomatal stress. The direct comparison of the stomatal resistances calculated from measured and modelled variables (Fig. c) shows an improvement of the modelled resistances (comparing DEF and REV). However, the modelled daytime stomatal resistance is still too high compared to the measurements. This points to an underestimation of stomatal uptake by the model during the day. A small fraction can be explained by the direct effect of the stomatal soil moisture stress in the model which does not occur in reality since the citrus orchard was watered during the measurement campaign.
Contrastingly, in summer, the model underestimation of Vd is higher than in spring (Fig. b). As seen from the comparison of stomatal resistance values (Fig. d), the model underestimates the stomatal uptake. This is because the irrigation of the orchard leads to cooling sustained evapotranspiration and keeps f(T) low. Thus, in the model, a too-high temperature stress acts on the stomata. Moreover, neglecting the soil moisture stress on stomata would bring the stomatal resistance values closer since the irrigation at the site ensures a constant and high soil moisture. The irrigation of the citrus orchard during the day also enhances surface wetness and favours deposition at cuticles which cannot be captured by the model. estimate the stomatal contribution to only account for 20 %–45 % of the total daytime dry deposition flux during both seasons and point to soil deposition and reactions of ozone with NO and VOCs as major sinks at the citrus orchard, especially during flowering season. The contribution of these pathways is expected to be enhanced by the inclusion of further BVOCs within the chemical mechanism and the explicit parameterization of in-canopy residence and transport.
Diurnal cycles of measured (obs) and modelled (DEF, REV, REV-fTfD) ozone dry deposition velocity and stomatal resistance in spring and summer 2010 at Lindcove research station.
Diurnal cycles of measured (obs) and modelled (DEF, REV, REV-NNTR: free-running f(T) and f(VPD) for evapotranspiration) ozone dry deposition velocities in wet and dry seasons at ATTO (gray: standard deviation).
Multiyear (2010–2015) mean absolute values in boreal summer.
Multiyear (2010–2015) mean absolute changes in boreal summer: i.e. difference between the revised and default scheme (REV – DEF).
Mean changes (2010) of dry deposition velocity in boreal summer. (a)f(Ws) modification; (b) temperature and VPD stress.
Tropical forests are known to be effective O3 sinks with observed mean midday maximum dry deposition velocity of 2.3 cms-1 due to much higher LAI compared to other sites (e.g. Lindcove). The measured dry deposition velocity at ATTO shown in Fig. a and b is no exception but shows a high variability (standard deviation). The diurnal cycle follows the solar radiation with maximum Vd at 15:00 UTC and highest amplitude during the wet season (April–May 2018). The amplitude of the diurnal cycle is highly underestimated in both EMAC simulations, with the highest mismatch during daytime. This is similar for other models. In fact, report a general and large underestimation of dry deposition velocities by models over tropical forests with highest predicted values of 0.25 cms-1. Here, the simulation with the revised dry deposition scheme (REV) shows only a minor increase of Vd during the wet season. Since stomatal uptake is known to be an important daytime sink , the underestimation of the total dry deposition flux is partly attributed to a too-low simulated stomatal uptake caused by the overestimation of temperature and the underestimation of relative humidity (Fig. ). The increase of dry deposition velocity by the new scheme is mainly due to the lowered soil moisture stress on stomata (f(Ws)) shown in Fig. e. also links stomatal uptake to the efficiency of turbulent mixing in transporting ozone down to the canopy. In general, 10 % of the total ozone sink during daytime and 39 % during night are associated with in-canopy processes . and identified the oxidation of sesquiterpenes as an important contributor to the chemical nighttime sink. Cuticular deposition might also play a role in humid conditions during night , which is underestimated by the model due to the biased relative humidity (Fig. c).
The uncertainty introduced by the mismatching meteorology becomes even more obvious when comparing measurements and simulations for November 2015. This month was characterized by temperatures of 2 to 3 degreeC above average and unusually little rainfall (compared to usual conditions in this season) due to a strong El Niño event . The dryness is overestimated by the model with a too-high temperature (Δ=+5 to +8K), too-low relative humidity (Δ=-30 % to -40 %)) and too-dry soil. The lack of available soil moisture (f(Ws)=0) effectively shuts down stomatal deposition in the default simulation (DEF), whereas the modification of the soil moisture stress function (neglecting the artificial lower limit; see Eq. ) in the revised model (REV) allows for an increased deposition (Fig. b). The temperature and relative humidity biases result in corresponding mismatching stress factors for the stomata that are double the ones derived from reanalysis data (Fig. f). This mismatch leads to an underestimation of stomatal uptake. This result is confirmed by the sensitivity simulation REV-NNTR for which no meteorological nudging has been applied, and the stress factors f(T) and f(VPD) are also used for the calculation of evapotranspiration. The REV-NNTR simulation yields much more realistic results compared to the measurements, capturing at least 50 % of the measured Vd during the day (Fig. b). This improvement is partly due to the omission of nudging, as the latter can have a detrimental effect on precipitation and evaporation . The temperature bias of the model is associated with the missing soil moisture buffer simulated by the bucket scheme. Incorporating a 5-layer scheme has been shown to lead to a more realistic soil water storage capacity, especially in the Amazon, and to a removal of this bias . Nevertheless, the REV-NNTR simulation suggests that the stress factors f(T) and f(VPD) significantly contribute to buffer soil moisture and ameliorate the dryness bias.
Relative change [%] and absolute change [Tg yr-1] (numbers on bars) of annual global loss by dry deposition of O3, SO2, HNO3 and HCHO (REV – DEF).
Relative change of multiyear (2010–2015) mean at ground level (DEF – REV).
Relative change of multiyear (2010–2015) zonal mean (DEF – REV).
Relative change of multiyear (2010–2015) boreal summer mean (DEF – REV).
Global impact on ground-level ozone
Given the importance of dry deposition for ground-level ozone and the uncertainty of dry deposition parameterizations in models , the global impact of the implemented code changes is assessed in this section.
The global (boreal) summer mean distributions of deposition velocity and ground-level mixing ratio for O3 shown in Fig. a–b are generally in the same range as reported for global models (e.g. ). However, like most global models, EMAC overestimates tropospheric ozone in comparison to satellite observations . Applying the revised dry deposition scheme increases the mean summer Vd by up to 0.5 cms-1 (Fig. a). The highest fraction of this increase arises from the inclusion of cuticular uptake at wet surfaces (Vcut,w) (Fig. b). The effect is large over the most northern continental regions (Fig. d) and even more pronounced where LAI is high like in Scandinavia and eastern Canada (for LAI distribution, see Fig. a). Additionally, the uptake at dry surfaces (Vcut,d) is enhanced with up to 0.3 cms-1 higher dry deposition velocity (Fig. c). This is because the default scheme applies a very high constant resistance for this process.
Concerning the stomatal deposition, the impacts of three different stress factors are considered. First, over relatively dry soil, i.e. where soil moisture exceeds 35 % of field capacity (wilting point of plants), the soil moisture stress is reduced by the modified parameterization. Neglecting the plants' wilting point as the lower limit for soil moisture stress on stomata weakens the dependency on field capacity. Thus, dry deposition is enhanced by up to 0.32 cms-1, as illustrated in Fig. a. Second, the inclusion of temperature and (third) VPD adjustment factors, indeed, leads to a spatially varying impact of ±0.27 cms-1 change in Vd (Fig. b). In humid and cold temperate regions, like Siberia and Canada, no temperature stress appears and the VPD adjustment factor increases the stomatal uptake. In the eastern US, Kazakhstan and central Amazon during boreal summer, stomata are stressed by temperature and VPD. This effect is overpredicted by the model, as the humidity over the Amazon forest is probably too low in the model (see Fig. ). The stress factors are shown in Fig. c and d.
However, the overall decrease in ozone concentration dampens the impact of the change in dry deposition flux.
In total, the changes by the revised dry deposition scheme increase the multiyear mean (2010–2015) loss of ozone by dry deposition from 946 to 1001 Tgyr-1. Accordingly, (boreal) summer ground-level ozone over land is reduced by up to 12 ppb (24 %), peaking over Scandinavia, Asia, central Africa and eastern Canada (Fig. b). In the Northern Hemisphere, also the zonal mean of the tropospheric ozone mixing ratio show a noticeable reduction far from the ground compared to the default scheme (Fig. a). This has the potential to reduce the positive bias of tropospheric ozone on the Northern Hemisphere (20 %) reported by .
However, besides ozone, also other atmospheric tracer gases are affected by the change in dry deposition. The global annual dry deposition flux of odd oxygen (Ox)
, which includes many important tropospheric trace gases, increases from 978 to 1032 Tgyr-1 due to the revision. This is in good agreement with the reported numbers by and . In Fig. , we show additionally the absolute and relative change of the multiyear annual average dry deposition loss of SO2, NO2, HNO3 and HCHO. As a very soluble species, the loss of SO2 is increased by the revised dry deposition scheme, whereas the predefined low cuticular and wet skin resistance of HNO3 in the old scheme were replaced with the new mechanism, leading to an decrease in dry deposition. The altered loss of NO2 and HCHO and other ozone precursors at ground level, especially soluble oxygenated VOCs, contributes to the total change in ozone loss. NO2 is deposited almost 40 % more significantly, contributing to the net reduction in ozone production but is mostly counterbalanced by other processes. The change of HCHO dry deposition flux is small on a global and annual scale and only important regionally, mostly in (boreal) summer, when it decreases HCHO at ground level (Fig. b) by up to 25 %. Thereby, the change in wet uptake is highest but is partially counterbalanced by other effects.
This leads to lower HO2 production from HCHO photooxidation and lower NO-to-NO2 conversion and thus lower ozone production . These effects also impact the OH mixing ratio (Figs. b, b) which controls the methane lifetime predicted by the model. However, for a clearer effect, a longer simulated time period would be needed. A detailed analysis of the trace gas budgets is beyond the scope of this paper and will be investigated in a subsequent study.
Uncertainties in modelling stomatal conductance
Dry deposition is a highly uncertain term in modelling ozone pollution . Its representation is generally limited by a lack of measurements and process understanding but also to a large extent driven by the quality of land cover information . Although the dry deposition scheme by is commonly used in global and regional models (e.g. MOZART, GEOS-Chem), the approach has some constraints .
The disadvantage of the big-leaf approach used in MESSy is that a vertical variation of leaf properties, affecting, for instance, the attenuation of solar radiation, is not considered e.g..
Regarding stomatal uptake, we neglect the mesophyll resistance as reactions inside the leaf are commonly assumed to not limit stomatal ozone uptake, whereas, besides mostly supporting laboratory studies e.g., a few contradicting findings exist e.g..
The here-used empirical multiplicative algorithm by for stomatal modelling has one general drawback concerning that the environmental responses to stomata are treated clearly in contrast to experimental evidence . However, Jarvis-type models have been shown to be able to compete with the semi-mechanistic Anet-gs models which link stomatal uptake to the CO2 assimilation during plant photosynthesis . The critics in state that the Jarvis model cannot capture the afternoon depression of ozone dry deposition is due to the original used VPD stress factor which has been replaced here by a mechanistic one based on the optimized exchange of CO2 and water by plants .
Furthermore, a larger set of land cover types is expected to improve the vegetation-dependent variation of dry deposition. The parameters used to model dry deposition of stomata, cuticle and soil are biome dependent and using generalized ones like for the input cuticular resistance can lead to differences in dry deposition . Exemplary discrepancies for the stomatal conductance calculated with different parameter sets are shown in Fig. as the summer mean of 2010.
Thereby, the temperature stress factor has been calculated as in Eq. () using the obtained surface temperature by EMAC (Fig. a, c) and applied to the model (DEFAULT) stomatal conductance (Eq. ) with two different parameter sets for coniferous and mixed forest by
Used parameters: Tmin=0∘C, Topt=18∘C, Tmax=36∘C.
and
Used parameters: Tmin=-3∘C, Topt=21∘C, Tmax=42∘C.
. obtained the parameters from a set of measurements in mixed hardwood/coniferous forest in Washington.
In general, the parameters are related to measurements where the absolute values are influenced by multiple factors like genotype and local climatic conditions . So, for global modelling, mostly simplified parameters have to be used like in the European Monitoring and Evaluation Programme (EMEP) .
Absolute difference of stomatal conductance applied with the temperature stress factor calculated for two different parameter sets by (Simp) and (Zh) in a comparison with the here-used parameter set by (Jar).
Sensitivity to model resolution
The simulation of dry deposition depends on meteorology including boundary layer processes, radiation (cloud distribution and reflectivity) and ozone chemistry as well as on input fields like vegetation density (LAI) . Model horizontal resolution inherently affects the amplitude and distribution of (regridded) surface processes and the artificial dilution of ozone precursors that are emitted. This aspect is investigated here by analysing simulations at three different spatial resolutions: 2.8∘× 2.8∘, 1.9∘× 1.9∘ and 1.1∘× 1.1∘ (REST42, REST63, and REV (T106) in Table ).
In Fig. a, the resolution dependency is shown for the annual dry deposition flux of ozone on different continental regions. The annual dry deposition fluxes differ by up to 40 Tgyr-1 globally between the different resolutions, with highest dry deposition at high resolution (T106). For the Northern Hemisphere (and consequently globally), this difference is driven by the higher annual mean ground-level ozone compared to the lower resolutions (Fig. c). However, this effect cannot be disentangled from the effect of decreased dry deposition velocity on ground-level ozone. Globally, increasing differences in O3 are anti-correlated with relative humidity as shown in Fig. a (ρ=-0.8). The impact of humidity on ozone chemistry is considered to be relatively weak , but showed for the US that only dry deposition establishes the observed anti-correlation between ozone and relative humidity. A dominating positive correlation of the dry deposition flux with the velocity only occurs on the Southern Hemisphere extratropics (SH_exT), which is highest between T63 and T106 (Fig. c). This can be attributed to discrepancies in stomatal deposition (Fig. d) driven by differences in humidity which might be caused by different moisture cycles and transpiration.
Ozone and dry deposition at three different resolutions (T42: 2.8∘× 2.8∘, T63: 1.9∘× 1.9∘, T106: 1.1∘× 1.1∘) and the different regions: Northern Hemisphere extratropics (NH_exT: 90–30∘ N), tropics (30∘ N–30∘ S), Southern Hemisphere extratropics (SH_exT: 90–30∘ S) and the whole Earth (global).
Correlations of resolution dependent relative differences of ozone, dry deposition and meteorological variables for the whole Earth (global) and the Southern Hemisphere extratropics (SH_exT) for the four boreal seasons: spring (MAM), summer (JJA), autumn (SON) and winter (DJF).
Conclusion and recommendations
Dry deposition to the Earth's surface is a key process for the representation of ground-level ozone in global models. Its parameterizations constitutes a relevant part of the model uncertainty . Revising the dry deposition scheme of EMAC leads to an improved representation of surface ozone in regions with a positive model ozone bias (e.g. Europe). The highest increase in ozone dry deposition is due to the implementation of cuticular uptake whose contribution is important especially during night over moist surfaces. The extension of the stomatal uptake with temperature and VPD adjustment factors accounts for the desired link of plant activity to hydroclimate as recommended by . Especially in drought-stressed regions (e.g. citrus orchards), the dependence on vapour pressure deficit leads to a realistic depression of stomatal uptake at noon. Also the dependence of dry deposition on soil moisture has been modified since the current representation of soil moisture in the model is not satisfactory. Specifically, the model simulates a too-dry soil for the Amazon basin, causing stomatal closure and thus an underestimation of dry deposition (Sect. ). We have indications that the dry bias is a consequence of meteorological nudging in EMAC and also the missing representation of organized convection in the tropics . The sensitivity of the vegetation to droughts is comparably high in the Amazon region because the model soil cannot hold water in the catchment for a realistic time period and exhibits a memory effect . Deeper root zones or buffering of the soil moisture below the root zone would improve the water holding capacity . With an improved representation of soil moisture, the more realistic parameterization of the soil moisture stress on stomatal uptake could be re-enabled.
In general, the inclusion of the strong link between dry deposition and meteorology reveals some limitations of the dry deposition scheme associated with the inaccurate representation of local meteorology.
The results also indicate that an improved representation of important non-stomatal dry deposition like in-canopy reactions of ozone with volatile organic compounds (e.g. citrus orchards; Sect. ) would lower the positive model–observation discrepancy. This can be achieved with the inclusion of further BVOCs and an explicit parameterization of the transport dynamics in the boundary layer in model simulations . Explicit field measurements could foster further process understanding, which is required for a detailed process description within the models, especially over tropical rainforests. The seasonal variability of the simulated dry deposition velocity could be further improved by using as model input the time series of vegetation cover from imaging products which also capture land use changes and vegetation trend that are known to impact dry deposition significantly .
Outlook
The representation of gaseous dry deposition in MESSy will be further improved by using the MODIS time series of LAI which captures multi-annual vegetation changes.
As the next step of dry deposition modelling in MESSy, a biome-dependent dry deposition model coupled to CO2 assimilation (White et al., 2004) will be applied. Biome-dependent vegetation cover information, required for this scheme, is then provided by global input data which, however, represent only the annual cycle of vegetation. Coupling MESSy to the recently available dynamic vegetation model LPJ-GUESS, which provides detailed vegetation information with the temporal variability required for a climate model, could be a further improvement. By now, the one-way coupling of LPJ-GUESS as a MESSy submodel is only in the initial evaluation phase of the coupling with the atmospheric model (Forrest et al., 2020).
Default dry deposition scheme
The default dry deposition scheme of MESSy uses the following equations described in .
For surface resistance over vegetation (in sm-1),
1Rs,veg(X)=1Rcan+Rs,soil(X)+Rqbr,veg(X)+LAIrcut(X)+LAIrstom,corr(X)+rmes(X),
where Rcan(X), Rs,soil(X) and Rqbr,veg(X) are the in-canopy aerodynamic resistance, the soil resistance and the quasi-laminar boundary resistance at canopy scale (in sm-1). rcut(X), rstom,corr(X) and rmes(X) are the cuticular resistance, stomatal resistance and mesophyll resistance at leaf scale scaled with LAI (in m2m-2) to canopy scale.For stomatal resistance,
rstom,corr=rstom(PAR)fws⋅DH2OD(O3).
For soil moisture stress function,
f(Ws)=1Ws(t)≥Wcr(=75%)Ws(t)-WpwpWcr-WpwpWpwp<Ws(t)<Wcr0Ws(t)≤Wpwp(=35%).
For cuticular resistance,
rcut(X)=rcut(O3)10-5⋅H(O3)+sreac(O3),
where rcutO3=1×10-5sm-1, H(O3)=0.01 and sreac=1.For wet skin resistance,
Rws(O3)=1/3Rws(SO2)+10-7⋅H(O3)+sreac(O3)Rcut,w(O3)-1,
where Rws(O3)=2000sm-1 and Rws(SO2)=100sm-1.
Evapotranspiration
Plants play a key role in the water and energy cycle and thus contribute to the land–atmosphere coupling, which drives the global climate. In this context, transpiration is an important process, as plants lose water during the necessary CO2 uptake via their stomata. The amount depends on the aperture behaviour of the respective plant in the respective environmental conditions . Thus, the latent heat flux incorporates the canopy resistance. The formulation is based on the Monin–Obukhov stability theory:
E=ρCh|v|β(qa-hqs(Ts,ps))β=1+Ch|v|Rstomfws-1,
where ρ is the density of air, |v| is the absolute value of the horizontal wind speed and Ch is the transfer coefficient of heat, whereas ra=1/(Ch|v|).qs and qa are the saturation-specific humidity and the atmospheric specific humidity, whereas the relative humidity h at the surface limits the evapotranspiration from bare soil. β determines the ratio of transpiration between water-stressed plants (β<1) and well-watered plants (β=1) .
The formula for the canopy stomatal resistance Rstom is given in Eq. (). In order to adapt the transpiration to temperature and vapour pressure deficit, the T and VPD adjustment factors can be applied to Rstom inversely like in the new dry deposition scheme via izwet=1 in the VERTEX &CTRL namelist. The modification of the soil moisture stress function f(Ws) (old: Eq. ; new: Eq. ) affects evapotranspiration directly.
Measured and modelled (DEF, REV) annual cycles at Borden forest.
Measured (obs) and modelled (DEF, REV) multiyear (2010–2012) and REV-fTfD (2010) annual cycles at Hyytiälä.
Differences of meteorology between EMAC and ERA5 at ATTO.
Boreal summer mean vegetation and meteorological variables predicted by EMAC.
Code availability
The Modular Earth Submodel System (MESSy) is continuously further developed and applied by a consortium of institutions. The usage of MESSy and access to the source code is licensed to all affiliates of institutions which are members of the MESSy Consortium. Institutions can become a member of the MESSy Consortium by signing the MESSy Memorandum of Understanding. More information can be found on the MESSy Consortium Website http://www.messy-interface.org (last access: 17 August 2020). The code presented here has been based on MESSy version 2.54 and will be available in the next official release (version 2.55). The exact code version used to produce the results of this paper is archived in the MESSy code repository and can be made available to members of the MESSy community upon request.
Data availability
The measurement data from Ontario are freely available at http://data.ec.gc.ca/data/air/monitor/special-studies-of-atmospheric-gases-particles-and-precipitation-chemistry/borden-forest-ozone-and-sulphur-dioxide-dry-deposition-study with the “Open Government Licence-Canada” (https://open.canada.ca/en/open-government-licence-canada, last access: 14 November 2019). The measurement data from Hyytiälä (Creative Commons 4.0 Attribution (CC BY) license https://creativecommons.org/licenses/by/4.0/, last access: 5 May 2020) can be accessed at https://avaa.tdata.fi/web/smart/smear/download. The data from Lindcove station were provided by Silvano Fares . The dry deposition measurement data from the Amazon Tall Tower Observatory were provided by Matthias Sörgel and are available upon request. The used ERA5 global climate reanalysis by ECMWF is available through the Climate Data Store (https://cds.climate.copernicus.eu, ).
The supplement related to this article is available online at: https://doi.org/10.5194/gmd-14-495-2021-supplement.
Author contributions
DT (and AK) initiated and supervised the study. DT and TE discussed the model developments which were implemented by AK and TE. HO originally wrote the MESSy VERTEX vertical diffusion submodel. SF provided the measurement data from Lindcove and further related theoretical calculations. IM conducted the dry deposition measurements at Hyytiälä and gave related support. TE performed the EMAC simulations and the data analyses, prepared the figures and wrote the manuscript.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
The work described in this paper has received funding from the Initiative and Networking Fund of the Helmholtz Association through the project “Advanced Earth System Modelling Capacity (ESM)”. The content of this paper is the sole responsibility of the author(s) and it does not represent the opinion of the Helmholtz Association, and the Helmholtz Association is not responsible for any use that might be made of the information contained.
The author(s) acknowledge the Environment and Climate Change Canada and the United States Environmental Protection Agency for the provision of the dry deposition velocity data at the Borden forest research station. Moreover, the personnel at the SMEAR II station of INAR – Institute for Atmospheric and Earth System Research, University of Helsinki, Finland, is acknowledged.
Concerning the measurement data from Amazon Tall Tower, we thank the Instituto Nacional de Pesquisas da Amazônia (INPA) and the Max Planck Society for continuous support. We thank for the support by the German Federal Ministry of Education and Research (BMBF contracts 01LB1001A, 01LK1602B and 01LP1606B) and the Brazilian Ministério da Ciência, Tecnologia e Inovação (MCTI/FINEP contract 01.11.01248.00) as well as the Amazon State University (UEA), FAPEAM, LBA/INPA and SDS/CEUC/RDS-Uatumã.
The measurements were conducted by Matthias Sörgel, Anywhere Tsokankunku, Stefan Wolff and Rodrigo Souza. For the usage of data from the ERA5 global climate reanalysis (generated using Copernicus Atmosphere Monitoring Service Information, 2020), we acknowledge the Copernicus Climate Change and Atmosphere Monitoring Service (https://apps.ecmwf.int/datasets/licences/copernicus/, last access: 10 March 2020). Neither the European Commission nor ECMWF is responsible for any use that may be made of the Copernicus information or data it contains.
Financial support
The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
Review statement
This paper was edited by Jason Williams and reviewed by Dennis Baldocchi and two anonymous referees.
ReferencesAltimir, N., Kolari, P., Tuovinen, J.-P., Vesala, T., Bäck, J., Suni, T., Kulmala, M., and Hari, P.: Foliage surface ozone deposition: a role for surface moisture?, Biogeosciences, 3, 209–228, 10.5194/bg-3-209-2006, 2006.Andersson, C. and Engardt, M.: European ozone in a future climate: Importance
of changes in dry deposition and isoprene emissions, J. Geophys.
Res.-Atmos., 115, D02303, 10.1029/2008JD011690, 2010.
Baldocchi, D. D., Hicks, B. B., and Camara, P.: A canopy stomatal resistance
model for gaseous deposition to vegetated surfaces, Atmos. Environ., 21, 91–101, 1987.
Ball, J. T., Woodrow, I. E., and Berry, J. A.: A model predicting stomatal
conductance and its contribution to the control of photosynthesis under
different environmental conditions, in: Progress in photosynthesis research, Springer, Dordrecht, 221–224, 1987.Bourtsoukidis, E., Behrendt, T., Yañez-Serrano, A. M., Hellén, H.,
Diamantopoulos, E., Catão, E., Ashworth, K., Pozzer, A., Quesada, C.,
Martins, D., Sá, M., Araujo, A., Brito, J., Artaxo, P.,
Kesselmeier, J., Lelieveld, J., and Williams, J.: Strong sesquiterpene emissions from Amazonian soils,
Nat. Commun., 9, 1–11, 10.1038/s41467-018-04658-y, 2018.Clifton, O., Paulot, F., Fiore, A., Horowitz, L., Correa, G., Baublitz, C.,
Fares, S., Goded, I., Goldstein, A., 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, 10.1029/2020JD032398,
2020a.Clifton, O. E., Fiore, A. M., Munger, J., Malyshev, S., Horowitz, L.,
Shevliakova, E., Paulot, F., Murray, L., and Griffin, K.: Interannual
variability in ozone removal by a temperate deciduous forest, Geophys.
Res. Lett., 44, 542–552, 10.1002/2016GL070923, 2017.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, S. L., and Tai, A. P. K.: Dry deposition of ozone over land: processes, measurement, and modeling, Rev. Geophys., 58, e2019RG000670, 10.1029/2019RG000670,
2020b.
Collatz, G. J., Ribas-Carbo, M., and Berry, J.: Coupled photosynthesis-stomatal conductance model for leaves of C4 plants, Funct. Plant Biol., 19, 519–538, 1992.Combe, M., de Arellano, J. V.-G., Ouwersloot, H. G., and Peters, W.: Plant
water-stress parameterization determines the strength of land–atmosphere
coupling, Agr. Forest Meteorol., 217, 61–73,
10.1016/j.agrformet.2015.11.006, 2016.Copernicus: ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate, 10.24381/cds.adbb2d47, available at: https://cds.climate.copernicus.eu/cdsapp#!/home (last access: 8 April 2020),
2017.Coumou, D. and Rahmstorf, S.: A decade of weather extremes, Nature Clim. Change, 2, 491, 10.1038/NCLIMATE1452, 2012.
Cowan, I. and Farquhar, G.: Stomatal function in relation to leaf metabolism
and environment, Symp. Soc. Exp. Biol., 31, 471–505, 1977.Damour, G., Simonneau, T., Cochard, H., and Urban, L.: An overview of models of stomatal conductance at the leaf level, Plant Cell Environ., 33,
1419–1438, 10.1111/j.1365-3040.2010.02181.x, 2010.Deckert, R., Jöckel, P., Grewe, V., Gottschaldt, K.-D., and Hoor, P.: A quasi chemistry-transport model mode for EMAC, Geosci. Model Dev., 4, 195–206, 10.5194/gmd-4-195-2011, 2011.
Delworth, T. L. and Manabe, S.: The influence of potential evaporation on the
variabilities of simulated soil wetness and climate, J. Climate, 1,
523–547, 1988.
Emberson, L., Ashmore, M., Cambridge, H., Simpson, D., and Tuovinen, J.-P.:
Modelling stomatal ozone flux across Europe, Environ. Pollut., 109,
403–413, 2000.Fares, S., Savi, F., and Conte, A.: Measurement data at Lindcove Orange Orchard, FLUXNET2015 IT-Cp2 Castelporziano2, Dataset,
10.18140/FLX/1440233, 2012–2014.Fares, S., Weber, R., Park, J.-H., Gentner, D., Karlik, J., and Goldstein,
A. H.: Ozone deposition to an orange orchard: Partitioning between stomatal
and non-stomatal sinks, Environ. Pollut., 169, 258–266,
10.1016/j.envpol.2012.01.030, 2012.Fares, S., Matteucci, G., Mugnozza, G. S., 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,
10.1016/j.atmosenv.2012.11.007, 2013.Fisher, R., Williams, M., Da Costa, A. L., Malhi, Y., Da Costa, R., Almeida,
S., and Meir, P.: The response of an Eastern Amazonian rain forest to drought stress: results and modelling analyses from a throughfall exclusion
experiment, Glob. Change Biol., 13, 2361–2378,
10.1111/j.1365-2486.2007.01417.x, 2007.Freire, L., Gerken, T., Ruiz-Plancarte, J., Wei, D., Fuentes, J., Katul, G.,
Dias, N., Acevedo, O., and Chamecki, M.: Turbulent mixing and removal of
ozone within an Amazon rainforest canopy, J. Geophys. Res.-Atmos., 122, 2791–2811, 10.1002/2016JD026009, 2017.
Ganzeveld, L. and Lelieveld, J.: Dry deposition parameterization in a chemistry general circulation model and its influence on the distribution of reactive trace gases, J. Geophys. Res.-Atmos., 100,
20999–21012, 1995.
Ganzeveld, L., Lelieveld, J., and Roelofs, G.-J.: A dry deposition
parameterization for sulfur oxides in a chemistry and general circulation
model, J. Geophys. Res.-Atmos., 103, 5679–5694, 1998.Ganzeveld, L., Lelieveld, J., Dentener, F., Krol, M., Bouwman, A., and Roelofs, G.-J.: Global soil-biogenic NOx emissions and the role of canopy processes, J. Geophys. Res.-Atmos., 107, 4298, 10.1029/2001JD0012892002, 2002.
Giorgetta, M. A., Roeckner, E., Mauritsen, T., Bader, J., Crueger, T., Esch,
M., Rast, S., Kornblueh, L., Schmidt, H., Kinne, S., Hohenegger, C., Möbis, B., Krismer, T., Wieners, K.-H., and Stevens, B.: The atmospheric general circulation model ECHAM6-model description, Reports on Earth System Science, No.135, Max-Planck-Institut für Meteorologie, Hamburg, 2013.Hagemann, S. and Stacke, T.: Impact of the soil hydrology scheme on simulated
soil moisture memory in a GCM, Clim. Dynam., 44.7, 1731–1750, 10.1007/s00382-014-2221-6, 2013.Hardacre, C., Wild, O., and Emberson, L.: An evaluation of ozone dry deposition in global scale chemistry climate models, Atmos. Chem. Phys., 15, 6419–6436, 10.5194/acp-15-6419-2015, 2015.Hogg, A., Uddling, J., Ellsworth, D., Carroll, M. A., Pressley, S., Lamb, B.,
and Vogel, C.: Stomatal and non-stomatal fluxes of ozone to a northern mixed hardwood forest, Tellus B, 59, 514–525,
10.1111/j.1600-0889.2007.00269.x, 2007.Hoshika, Y., Osada, Y., De Marco, A., Penuelas, J., and Paoletti, E.: Global
diurnal and nocturnal parameters of stomatal conductance in woody plants and major crops, Global Ecol. Biogeogr., 27, 257–275,
10.1016/j.agrformet.2017.01.005, 2018.Hu, L., Jacob, D. J., Liu, X., Zhang, Y., Zhang, L., Kim, P. S., Sulprizio,
M. P., and Yantosca, R. M.: Global budget of tropospheric ozone: Evaluating
recent model advances with satellite (OMI), aircraft (IAGOS), and ozonesonde observations, Atmos. Environ., 167, 323–334,
10.1016/j.atmosenv.2017.08.036, 2017.Huang, L., McDonald-Buller, E. C., McGaughey, G., Kimura, Y., and Allen, D. T.: The impact of drought on ozone dry deposition over eastern Texas, Atmos. Environ., 127, 176–186, 10.1016/j.atmosenv.2015.12.022, 2016.Jacob, D. J. and Winner, D. A.: Effect of climate change on air quality,
Atmos. Environ., 43, 51–63, 10.1016/j.atmosenv.2008.09.051,
2009.
Jarvis, P.: The interpretation of the variations in leaf water potential and
stomatal conductance found in canopies in the field, Phil. Trans. R. Soc.
Lond. B, 273, 593–610, 1976.Jeuken, A. B. M., Siegmund, P. C., Heijboer, L. C., Feichter, J., and
Bengtsson, L.: On the potential of assimilating meteorological analyses in a global climate model for the purpose of model validation, J.
Geophys. Res.-Atmos., 101, 16939–16950, 10.1029/96JD01218, 1996.Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., Tost, H., Riede, H., Baumgaertner, A., Gromov, S., and Kern, B.: Development cycle 2 of the Modular Earth Submodel System (MESSy2), Geosci. Model Dev., 3, 717–752, 10.5194/gmd-3-717-2010, 2010.Jöckel, P., Tost, H., Pozzer, A., Kunze, M., Kirner, O., Brenninkmeijer, C. A. M., Brinkop, S., Cai, D. S., Dyroff, C., Eckstein, J., Frank, F., Garny, H., Gottschaldt, K.-D., Graf, P., Grewe, V., Kerkweg, A., Kern, B., Matthes, S., Mertens, M., Meul, S., Neumaier, M., Nützel, M., Oberländer-Hayn, S., Ruhnke, R., Runde, T., Sander, R., Scharffe, D., and Zahn, A.: Earth System Chemistry integrated Modelling (ESCiMo) with the Modular Earth Submodel System (MESSy) version 2.51, Geosci. Model Dev., 9, 1153–1200, 10.5194/gmd-9-1153-2016, 2016.
Jones, H.: Plants and Microclimate, Cambridge University Press,
Cambridge, 1992.Jülich Supercomputing Centre: JURECA: Modular supercomputer at
Jülich Supercomputing Centre, J. Large-Scale Res. Facilities, 4, p. 132, 10.17815/jlsrf-4-121-1, 2018.Katul, G. G., Palmroth, S., and Oren, R.: Leaf stomatal responses to vapour
pressure deficit under current and CO2-enriched atmosphere explained by the economics of gas exchange, Plant Cell Environ., 32, 968–979,
10.1111/j.1365-3040.2009.01977.x, 2009.Katul, G. G., Oren, R., Manzoni, S., Higgins, C., and Parlange, M. B.:
Evapotranspiration: a process driving mass transport and energy exchange in
the soil-plant-atmosphere-climate system, Rev. Geophys., 50, RG3002,
10.1029/2011RG000366, 2012.Kavassalis, S. C. and Murphy, J. G.: Understanding ozone-meteorology
correlations: A role for dry deposition, Geophys. Res. Lett., 44,
2922–2931, 10.1002/2016GL071791, 2017.Kerkweg, A., Buchholz, J., Ganzeveld, L., Pozzer, A., Tost, H., and Jöckel, P.: Technical Note: An implementation of the dry removal processes DRY DEPosition and SEDImentation in the Modular Earth Submodel System (MESSy), Atmos. Chem. Phys., 6, 4617–4632, 10.5194/acp-6-4617-2006, 2006.
Keronen, P., Reissell, A., Rannik, U., Pohja, T., Siivola, E., Hiltunen, V.,
Hari, P., Kulmala, M., and Vesala, T.: Ozone flux measurements over a Scots
pine forest using eddy covariance method: performance evaluation and
comparison with flux-profile method, Boreal Environ. Res., 8,
425–444, 2003.Kharol, S., Shephard, M., McLinden, C., Zhang, L., Sioris, C., O'Brien, J.,
Vet, R., Cady-Pereira, K., Hare, E., Siemons, J., and Krotkov, N. A.: Dry deposition of
reactive nitrogen from satellite observations of ammonia and nitrogen dioxide over North America, Geophys. Res. Lett., 45, 1157–1166,
10.1002/2017GL075832, 2018.
Klimarechenzentrum: The ECHAM3 atmospheric general circulation
model, Deutsches Klimarechenzentrum, Hamburg, Techn. Rep, 6, 1992.
Kraus, H.: Die Atmosphäre der Erde: Eine Einführung in die
Meteorologie, Springer-Verlag, Berlin, Heidelberg, 2007 (in German).
Lamaud, E., Carrara, A., Brunet, Y., Lopez, A., and Druilhet, A.: Ozone fluxes above and within a pine forest canopy in dry and wet conditions, Atmos. Environ., 36, 77–88, 2002.Lamaud, E., Loubet, B., Irvine, M., Stella, P., Personne, E., and Cellier, P.: Partitioning of ozone deposition over a developed maize crop between stomatal and non-stomatal uptakes, using eddy-covariance flux measurements and modelling, Agr. Forest Meteorol., 149, 1385–1396,
10.1016/j.agrformet.2009.03.017, 2009.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, 10.1029/2018GB006157, 2019.
Lu, Y.-S.: Propagation of land surface model uncertainties in simulated
terrestrial system states, PhD thesis, Bonner Meteorologische Abhandlungen
Heft 84, University of Bonn, Germany, 120 pp., 2018.Makar, P., Staebler, R., Akingunola, A., Zhang, J., McLinden, C., Kharol, S.,
Pabla, B., Cheung, P., and Zheng, Q.: The effects of forest canopy shading
and turbulence on boundary layer ozone, Nat. Commun., 8, 15243, 10.1038/ncomms15243,
2017.Mammarella, I., Rannik, Ü., and Launiainen, S.: SMEAR II Hyytiälä forest eddy covariance, Institute for Atmospheric and Earth System Research, available at: https://avaa.tdata.fi/web/smart/smear/download, last access: 12 June 2020.Mauritsen, T. and Stevens, B.: Missing iris effect as a possible cause of mutedhydrological change and high climate sensitivity in models, Nat,
Geosci., 8, 346–351, 10.1038/NGEO2414, 2015.Mauritsen, T., Stevens, B., Roeckner, E., Crueger, T., Esch, M., Giorgetta, M., Haak, H., Jungclaus, J., Klocke, D., Matei, D., Mikolajewicz, U., Notz, D., Pincus, R., Schmidt, H., and Tomassini, L.: Tuning the climate of a global model, J. Adv. Model. Earth Sy., 4, M00A01,
10.1029/2012MS000154, 2012.Mészáros, R., Horváth, L., Weidinger, T., Neftel, A., Nemitz, E., Dämmgen, U., Cellier, P., and Loubet, B.: Measurement and modelling ozone fluxes over a cut and fertilized grassland, Biogeosciences, 6, 1987–1999, 10.5194/bg-6-1987-2009, 2009.National Centers for Environmental Information: State of the Climate:
Global Climate Report for Annual 2015, Tech. rep., available at: https://www.ncdc.noaa.gov/sotc/global/201513 (last access: 3 March 2020), 2016.Ran, L., Pleim, J., Song, C., Band, L., Walker, J. T., and Binkowski, F. S.: A photosynthesis-based two-leaf canopy stomatal conductance model for
meteorology and air quality modeling with WRF/CMAQ PX LSM, J.
Geophys. Res.-Atmos., 122, 1930–1952,
10.1002/2016JD025583, 2017.Rannik, Ü., Altimir, N., Mammarella, I., Bäck, J., Rinne, J., Ruuskanen, T. M., Hari, P., Vesala, T., and Kulmala, M.: Ozone deposition into a boreal forest over a decade of observations: evaluating deposition partitioning and driving variables, Atmos. Chem. Phys., 12, 12165–12182, 10.5194/acp-12-12165-2012, 2012.Rannik, Ü., Peltola, O., and Mammarella, I.: Random uncertainties of flux measurements by the eddy covariance technique, Atmos. Meas. Tech., 9, 5163–5181, 10.5194/amt-9-5163-2016, 2016.Righi, M., Eyring, V., Gottschaldt, K.-D., Klinger, C., Frank, F., Jöckel, P., and Cionni, I.: Quantitative evaluation of ozone and selected climate parameters in a set of EMAC simulations, Geosci. Model Dev., 8, 733–768, 10.5194/gmd-8-733-2015, 2015.
Roeckner, E., Bäuml, G., Bonaventura, L., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kirchner, I., Kornblueh, L., Manzini, E., Rhodin, A., Schlese, U., Schulzweida, U., and Tompkins, A.: The atmospheric general circulation model ECHAM 5. PART I: Model description, report no. 349, Max Planck Institute for Meteorology, Hamburg, 2003.Rummel, U., Ammann, C., Kirkman, G. A., Moura, M. A. L., Foken, T., Andreae, M. O., and Meixner, F. X.: Seasonal variation of ozone deposition to a tropical rain forest in southwest Amazonia, Atmos. Chem. Phys., 7, 5415–5435, 10.5194/acp-7-5415-2007, 2007.
Schulz, J.-P., Dümenil, L., and Polcher, J.: On the land
surface–atmosphere coupling and its impact in a single-column atmospheric
model, J. Appl. Meteorol., 40, 642–663, 2001.
Schwede, D., Zhang, L., Vet, R., and Lear, G.: An intercomparison of the
deposition models used in the CASTNET and CAPMoN networks, Atmos.
Environ., 45, 1337–1346, 2011.
Seinfeld, J. H. and Pandis, S. N.: Atmospheric chemistry and physics: from air pollution to climate change, John Wiley & Sons, Hoboken, NJ, 2016.
Sellers, P., Mintz, Y., Sud, Y. E. A., and Dalcher, A.: A simple biosphere
model (SiB) for use within general circulation models, J.
Atmos. Sci., 43, 505–531, 1986.
Sellers, P., Dickinson, R. E., Randall, D., Betts, A., Hall, F., Berry, J.,
Collatz, G., Denning, A., Mooney, H., Nobre, C., Sato, N., Field, C. B., and Henderson-Sellers, A.: Modeling the
exchanges of energy, water, and carbon between continents and the atmosphere, Science, 275, 502–509, 1997.
Sellers, P. J.: Canopy reflectance, photosynthesis and transpiration,
Int. J. Remote Sens., 6, 1335–1372, 1985.Seneviratne, S. I., Corti, T., Davin, E. L., Hirschi, M., Jaeger, E. B.,
Lehner, I., Orlowsky, B., and Teuling, A. J.: Investigating soil
moisture–climate interactions in a changing climate: A review, Earth-Sci. Rev., 99, 125–161, 10.1016/j.earscirev.2010.02.004, 2010.
Silva, S. J. and Heald, C. L.: Investigating dry deposition of ozone to
vegetation, J. Geophys. Res.-Atmos., 123, 559–573, 2018.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, 10.5194/acp-12-7825-2012, 2012.Solberg, S., Hov, Ø., Søvde, A., Isaksen, I., Coddeville, P., De Backer, H., Forster, C., Orsolini, Y., and Uhse, K.: European surface ozone in the extreme summer 2003, J. Geophys. Res.-Atmos., 113, D07307,
10.1029/2007JD009098, 2008.
Stella, P., Loubet, B., Lamaud, E., Laville, P., and Cellier, P.: Ozone
deposition onto bare soil: a new parameterisation, Agr. Forest
Meteorol., 151, 669–681, 2011.
Stella, P., Loubet, B., de Berranger, C., Charrier, X., Ceschia, E., Gerosa,
G., Finco, A., Lamaud, E., Serça, D., George, C., and Ciuraru, R.: Soil ozone
deposition: Dependence of soil resistance to soil texture, Atmos.
Environ., 199, 202–209, 2019.Stephens, G. L., Li, J., Wild, M., Clayson, C. A., Loeb, N., Kato, S.,
L'ecuyer, T., Stackhouse, P. W., Lebsock, M., and Andrews, T.: An update on
Earth's energy balance in light of the latest global observations, Nat.
Geosci., 5, 691–696, 10.1038/NGEO1580, 2012.Stevens, B. and Schwartz, S. E.: Observing and modeling Earth’s energy flows,
Surv. Geophys., 33, 779–816, 10.1007/s10712-012-9184-0, 2012.Sulis, M., Langensiepen, M., Shrestha, P., Schickling, A., Simmer, C., and
Kollet, S. J.: Evaluating the influence of plant-specific physiological
parameterizations on the partitioning of land surface energy fluxes, J. Hydrometeorol., 16, 517–533, 10.1175/JHM-D-14-0153.1, 2015.Sun, S., Moravek, A., von der Heyden, L., Held, A., Sörgel, M., and Kesselmeier, J.: Twin-cuvette measurement technique for investigation of dry deposition of O3 and PAN to plant leaves under controlled humidity conditions, Atmos. Meas. Tech., 9, 599–617, 10.5194/amt-9-599-2016, 2016.Tuovinen, J.-P., Emberson, L., and Simpson, D.: Modelling ozone fluxes to
forests for risk assessment: status and prospects, Ann. Forest Sci., 66, 1–14, 10.1051/forest/2009024, 2009.Tuzet, A., Perrier, A., Loubet, B., and Cellier, P.: Modelling ozone deposition fluxes: The relative roles of deposition and detoxification processes, Agr. Forest Meteorol., 151, 480–492,
10.1016/j.agrformet.2010.12.004, 2011.Val Martin, M., Heald, C., and Arnold, S.: Coupling dry deposition to
vegetation phenology in the Community Earth System Model: Implications for
the simulation of surface O3, Geophys. Res. Lett., 41, 2988–2996,
10.1002/2014GL059651, 2014.
Van Pul, W. and Jacobs, A.: The conductance of a maize crop and the underlying soil to ozone under various environmental conditions, Bound.-Lay. Meteorol., 69, 83–99, 1994.Wang, Y., Sperry, J. S., Anderegg, W. R., Venturas, M. D., and Trugman, A. T.: A theoretical and empirical assessment of stomatal optimization modeling, New Phytol., 227, 311–325, 10.1111/nph.16572, 2020.
Wesely, M.: Parameterization of surface resistances to gaseous dry deposition
in regional-scale numerical models, Atmos. Environ., 23, 1293–1304, 1989.
Wesely, M. and Hicks, B.: Some factors that affect the deposition rates of
sulfur dioxide and similar gases on vegetation, J. Air Pollut. Con. Ass., 27, 1110–1116, 1977.
Wesely, M. and Hicks, B.: A review of the current status of knowledge on dry
deposition, Atmos. Environ., 34, 2261–2282, 2000.Wong, A. Y. H., Geddes, J. A., Tai, A. P. K., and Silva, S. J.: Importance of dry deposition parameterization choice in global simulations of surface ozone, Atmos. Chem. Phys., 19, 14365–14385, 10.5194/acp-19-14365-2019, 2019.Wu, Z., Staebler, R., Vet, R., and Zhang, L.: Dry deposition of O3 and SO2
estimated from gradient measurements above a temperate mixed forest,
Environ. Pollut., 210, 202–210, 10.1016/j.envpol.2015.11.052,
2016.Wu, Z., Schwede, D. B., Vet, R., Walker, J. T., Shaw, M., Staebler, R., and
Zhang, L.: Evaluation and intercomparison of five North American dry
deposition algorithms at a mixed forest site, J. Adv. Model. Earth Sy., 10, 1571–1586, 10.1029/2017MS001231, 2018.Young, P. J., Naik, V., Fiore, A. M., Gaudel, A., Guo, J., Lin, M., Neu, J.,
Parrish, D., Rieder, H., Schnell, J., Tilmes, S., Wild, O., Zhang, L., Ziemke, J., Brandt, J., Delcloo, A., Doherty, R. M., Geels, C., Hegglin, M. I., Hu, L., Im, U., Kumar, R., Luhar, A., Murray, L., Plummer, D., Rodriguez, J., Saiz-Lopez, A., Schultz, M. G., Woodhouse, M. T., and Zeng, G.: Tropospheric Ozone Assessment
Report: Assessment of global-scale model performance for global and regional ozone distributions, variability, and trends, Elem. Sci. Anth., 6, 10, 10.1525/elementa.265, 2018.
Zhang, L., Brook, J. R., and Vet, R.: On ozone dry deposition – with
emphasis on non-stomatal uptake and wet canopies, Atmos. Environ.,
36, 4787–4799, 2002.Zhang, L., Brook, J. R., and Vet, R.: A revised parameterization for gaseous dry deposition in air-quality models, Atmos. Chem. Phys., 3, 2067–2082, 10.5194/acp-3-2067-2003, 2003.