Articles | Volume 11, issue 3
Model description paper
08 Mar 2018
Model description paper |  | 08 Mar 2018

Description and validation of the Simple, Efficient, Dynamic, Global, Ecological Simulator (SEDGES v.1.0)

Pablo Paiewonsky and Oliver Elison Timm

In this paper, we present a simple dynamic global vegetation model whose primary intended use is auxiliary to the land–atmosphere coupling scheme of a climate model, particularly one of intermediate complexity. The model simulates and provides important ecological-only variables but also some hydrological and surface energy variables that are typically either simulated by land surface schemes or else used as boundary data input for these schemes. The model formulations and their derivations are presented here, in detail. The model includes some realistic and useful features for its level of complexity, including a photosynthetic dependency on light, full coupling of photosynthesis and transpiration through an interactive canopy resistance, and a soil organic carbon dependence for bare-soil albedo. We evaluate the model's performance by running it as part of a simple land surface scheme that is driven by reanalysis data. The evaluation against observational data includes net primary productivity, leaf area index, surface albedo, and diagnosed variables relevant for the closure of the hydrological cycle. In this setup, we find that the model gives an adequate to good simulation of basic large-scale ecological and hydrological variables. Of the variables analyzed in this paper, gross primary productivity is particularly well simulated. The results also reveal the current limitations of the model. The most significant deficiency is the excessive simulation of evapotranspiration in mid- to high northern latitudes during their winter to spring transition. The model has a relative advantage in situations that require some combination of computational efficiency, model transparency and tractability, and the simulation of the large-scale vegetation and land surface characteristics under non-present-day conditions.


Figure 1Variable dependencies and updating in SEDGES


1 Introduction

Simulation of the land surface is a critical component of models of the Earth climate system. In such models, heat, energy, and moisture transfer between the land and atmosphere or ocean are calculated with information from the hydrological and vegetation components of these models. The Simple, Efficient, Dynamic, Global, Ecological Simulator (SEDGES) simulates the gross properties of vegetation, as well as some large-scale land surface characteristics, which are important for the simulation of energy and moisture exchange with the atmosphere. These simulated variables are to be used by a land surface scheme in which SEDGES is embedded in a climate or Earth system model. Models of similar complexity to SEDGES exist that simulate the vegetative if not also the hydrological and energy-transferring aspects of the land surface, include the Vegetation Continuous Description model (VECODE) (Brovkin et al.1997, 2002) and the Efficient Numerical Terrestrial Scheme (ENTS) (Williamson et al.2006). The purpose of these kinds of models is to efficiently and reasonably simulate dynamic land surface characteristics and behavior and provide these to the atmospheric components of Earth system models of intermediate complexity. In such a framework, the need for sophisticated simulations of the land surface is obviated by the simplifications present in the other model components (which reduce the benefit of a more realistic representations of land surface processes as compared to when the other model components are complex and realistic), by a desire to more easily understand the processes underlying experimental results, and by the computational burden that goes along with the increased complexity.

SEDGES is a new model that is based on the original SIMulator for Biospheric Aspects (SimBA) model (Kleidon2006b), which is coupled to the Planet Simulator (PlaSim) general circulation model (GCM) (Lunkeit et al.2007), and it is even more strongly based on a later version of SimBA that the lead author of this paper developed (Lunkeit et al.2011) and which is also coupled to PlaSim. Neither version of SimBA has been thoroughly evaluated. However, when coupled to PlaSim, the earlier version's net primary productivity was shown to be broadly reasonable in Kleidon (2006b), and it has been used successfully in a number of studies (e.g., Kleidon2006a, b; Bowring et al.2014). The SimBA model is based on a light and water-use efficiency model for crops (Monteith et al.1989) that was later adapted and expanded to forest canopies (Dewar1997). This approach to vegetative productivity is maintained in SEDGES and helps form the core of the model. Other aspects of SEDGES's structure that are particular to SEDGES (i.e., are uncommon in or absent from other models) are carried over from the original SimBA, including the determination of soil water-holding capacity and moist-soil leaf cover fraction by the amount of vegetation biomass (Fig. 1).

SEDGES builds upon SimBA by improving most of its parameterizations. Some of the updated and expanded model features are taken or adapted from currently existing models. Although SEDGES makes many improvements upon the original SimBA, it is still a simple model; i.e., it is not of “intermediate” complexity (cf. Willeit and Ganopolski2016). Compared to the original SimBA, SEDGES (and the land surface framework that it presupposes) has four major increases in complexity: the separation of evapotranspiration into soil and vegetative components, the inclusion of aerodynamic conductance in the formulation for carbon uptake by vegetation, full coupling of photosynthesis and transpiration through interactive canopy conductance, and soil organic carbon-dependent soil albedo. Choosing the appropriate level of complexity of represented processes within a land surface or dynamic global vegetation model often depends on the context in which the model is to be used and involves subjective weighting of trade-offs between increased accuracy and realism on the one hand and, on the other hand, robustness to the chosen values of poorly constrained parameters and model reliability in a wide range of situations (Prentice et al.2015).

This paper presents the SEDGES model's structure, equations, and ability to simulate ecological and hydrological variables when used in conjunction with a simple soil hydrological scheme and parameterization for aerodynamic conductance and forced by reanalysis data.

2 Model description: SEDGES

2.1 Overview of SEDGES

SEDGES is a simple, dynamic global vegetation model (DGVM). Within the context of the history of land surface modeling, the SEDGES framework (defined as SEDGES and the type of land surface model that it presupposes that it forms a part of) combines aspects of first- and third-generation models (Sellers et al.1997; Pitman2003; Prentice et al.2015), which include, respectively, the use of a simple single-layer “bucket” model for soil hydrology and the full coupling of photosynthesis and transpiration through interactive canopy conductance.

The intended main use of SEDGES is for it to be embedded within a land surface scheme and to thus help the land surface scheme simulate large-scale properties and behavior of the surface in which vegetation and soil play a role. These large-scale properties and behavior include interactions between the land surface and the atmosphere through simulated CO2, moisture, energy, heat, and momentum exchanges. SEDGES either directly simulates or else provides the greater land surface scheme with the necessary variables. SEDGES also simulates ecological variables (e.g., biomass, soil organic carbon, forest cover fraction, leaf cover fraction) that are not directly used outside it by either the rest of the land surface scheme or by the atmosphere, but these ecological variables almost always have a downstream impact on the simulations of the aforementioned land–atmosphere exchanges. SEDGES has been designed, in particular, for coupling with the Planet Simulator, and thus presupposes that the land surface scheme that it forms a part of be similar to that of PlaSim in its hydrology and scheme for surface evaporation. This section presents, in detail, the SEDGES model formulations for its output ecological and land surface variables, including their derivations.

Overall, sub-grid-scale heterogeneity is treated similarly in SEDGES as it is in the SLAM-1T (Simple Land-Atmosphere Mosaic – 1 Tile) configuration of CHASM (CHAmeleon Surface Model) (Desborough1999) and in the original Met Office Surface Exchange Scheme (MOSES) land surface scheme (Cox et al.1999). Because of its need to couple with the Planet Simulator GCM, which uses a bulk aerodynamic formulation with a single tile (i.e., one surface type per grid cell), SEDGES handles sub-grid-scale heterogeneity according to what is called the “average parameter” method (Giorgi and Avissar1997) or the “composite” approach (Li and Arora2012; Melton and Arora2014). In this approach, land surface characteristics are aggregated to provide a single, representative value on the scale of the entire grid cell. The framework in SEDGES for evapotranspiration (ET) is special because it also qualifies as a simplified mosaic (or “mixed”, Li and Arora2012) approach, such that only surface conductance differs between the surface types (or tiles). This mosaic approach is also used to extend a “big-leaf” formulation for vegetative CO2 uptake (Appendix B of Raupach1998) to a mixed soil–vegetation surface, enabling us to isolate the fraction of total ET that is due to transpiration and thus formulate vegetative control of water loss and CO2 uptake. The SEDGES framework neglects evaporation from intercepted canopy water and thus only distinguishes between two tiles when snow is absent: bare soil and vegetation. When there is snow cover, there are essentially three tiles: snow-free exposed bare soil, snow-covered surface, and snow-free exposed vegetation (see Sect. 2.3.1 for more detail).

Figure 2Canopy resistance dependencies and updating in SEDGES


Table 1SEDGES variables in the paper. Notation for the rightmost column is as follows: “I”, “O”, and “M” indicate variables that are, respectively, input to, output by, or modified by the SEDGES model. A blank field in that column indicates that the variable could be output if desired, but some code changes would be needed. “EI” denotes a variable from ERA-Interim that is used by the simple land surface scheme that we use to drive SEDGES with in Sect. 4. We denote by “M” only those output variables that are output by SEDGES that are expected to be already existing in the land surface scheme in which SEDGES is embedded.

Download Print Version | Download XLSX

Franks et al. (2013)Pisek et al. (2010); He et al. (2012)Raupach (1998)van de Griend and Owe (1994)Knorr (2000)Jenkinson et al. (1990)Roesch et al. (2001)Roeckner et al. (2003)Moody et al. (2007)Oke (1987)(Sellers et al.1996b)

Table 2SEDGES parameters in the paper.

Download Print Version | Download XLSX

The core driving variable in SEDGES is gross primary productivity (GPP), which impacts the living biomass, which, in turn, heavily influences or determines the values of almost every other simulated land surface variable (listed in Sect. 3). Forest cover and leaf cover fractions and (implicitly) rooting depth are parameterized to increase with biomass using respective relationships that are fixed. SEDGES (as well as SimBA) uses a simple single-layer bucket for the soil hydrology and has essentially one plant functional type. Soil organic carbon is treated as a single homogeneous reservoir. Surface albedo depends on soil organic carbon content, snow depth, biomass, and leaf cover fraction. Surface roughness increases monotonically with biomass. Gross primary productivity and canopy resistance1 are coupled through the vegetation's effort to conjointly minimize water loss and satisfy photosynthetic demand for CO2, as per the original Dewar (1997) model (see Sect. 2.2.6). Table 1 lists all SEDGES variables, their units, and their relationship with the SEDGES vegetation model. Their dependencies and how those variables are updated is shown in Figs. 1 and 2. Table 2 lists all parameters that are used in the paper along with their values and their units.

With its main use and design already indicated, SEDGES can also be forced offline with external data. The offline mode of SEDGES was used in combination with a simple soil hydrological scheme and parameterization for aerodynamic conductance and run using reanalysis climate data for the SEDGES model evaluation (see Sect. 4).

In Sect. 2.2 and 2.3, we describe in detail the equations of the ecological and physical variables in the SEDGES model. Because SEDGES uses the original SimBA model (Kleidon2006b) (and its code) as a basis, we explicitly mention in the text when significant changes are introduced in SEDGES compared with SimBA. We provide extra detail on the original SimBA for the cases in which Kleidon (2006b) provides insufficient information for comparison with SEDGES. The broadest structural increases in complexity in SEDGES have already been identified in the introduction; our focuses are narrower in scope in these coming sections.

2.2 Equations for ecological variables

2.2.1 Vegetation biomass

The most important SEDGES variable is biomass of live vegetation. This quantity is used to directly or indirectly derive almost all the land surface variables; these are presented in the coming sections. The prognostic equation for biomass is as follows:

(1) dC veg d t = NPP - L ,

where Cveg is the carbon in living biomass (kg C m−2), L is the litterfall and equals Cvegτveg, τveg is the residence time of the vegetative carbon (see Sect. 2.2.9) and equals 10 years, and NPP (net primary productivity) is approximated as 0.5 GPP.

2.2.2 Net primary productivity and gross primary productivity

SEDGES uses a constant NPP  GPP = 0.5 approximation. This approximation is supported by the conservative nature of the ratio of mitochondrial respiration to gross photosynthesis and (hence) of the ratio of net photosynthesis to gross photosynthesis over a wide variety of conditions, on timescales of weeks or more (see the brief review in Van Oijen et al.2010). Since each model time step in SEDGES is much shorter than this, it thus might seem incorrect to hold NPP  GPP fixed for each time step. It is appropriate to do so, however, because the NPP-to-GPP ratio in SEDGES only impacts biomass changes, and the latter occur on very long timescales. Finally, meta-analyses of previous studies have found a robust (DeLucia et al.2007; Litton et al.2007) linear or proportional relationship between NPP and GPP, with slope or proportionality constant of around 0.5 (Gifford2003; DeLucia et al.2007; Litton et al.2007), albeit with considerable variation in NPP  GPP across field sites (Amthor and Baldocchi2001; DeLucia et al.2007; Litton et al.2007).

GPP is calculated as the minimum of a light-limited rate, GPPL, and a water-limited rate, GPPW (Monteith et al.1989; Dewar1997). That is, GPP =min(GPPL,GPPW.)

2.2.3 Light-limited gross primary productivity

GPPL uses a light use efficiency (LUE) formulation (e.g., Yuan et al.2007) whose overall structure is very similar to that of the original version of SimBA (Kleidon2006b). The equation is as follows:

(2) GPP L = ϵ max f 1 ( CO 2 ) f 2 ( T sfc ) f APAR SW ,

where ϵmax is a globally constant maximum LUE parameter (5.0×10-10 kg C J−1); f1(CO2) is a CO2 fertilization function (described below); f2(Tsfc) is a temperature limitation function (described below); fAPAR is the fraction of photosynthetically active radiation (PAR) that is absorbed by green vegetation (see below); and SW is the downward flux of shortwave radiation just above the canopy surface (in W m−2). (Instead of this term, the original version of SimBA (Kleidon2006b) uses the net surface shortwave radiation flux.)

In Eq. (2), the first term on the right-hand side, ϵmax, is the LUE with respect to the total shortwave broadband radiation that is absorbed by the photosynthetic parts of the vegetation. This constant term is the maximum efficiency with which incident shortwave radiation can be used to synthesize vegetative carbon at the reference (360 ppmv) CO2 level (also see explanation of fAPAR, below).

The second term in Eq. (2), f1(CO2), increases productivity with increasing CO2 and is taken directly from Eq. (5) of Franks et al. (2013). For CO2>CO2comp, we have

(3) f 1 ( CO 2 ) = CO 2 norm c a - CO 2 comp c a + 2 CO 2 comp ,

where ca is the atmospheric CO2 concentration (ppmv), CO2comp=40 ppmv is the light compensation point in the absence of dark respiration, and CO2norm=360+2CO2comp360-CO2comp. Otherwise, f1(CO2)=0. This fertilization function replaces an earlier “beta” factor approach used in SimBA (Lunkeit et al.2011). The latter does not saturate and therefore yields unrealistically large fertilization at high CO2 levels.

The third term in Eq. (2) is f2(Tsfc). f2(Tsfc) is a ramp function, reducing productivity linearly from a surface temperature of Tcrit≡20 to 0 C, with f2(Tsfc)=0 for below 0 C temperatures. The 20 C value is the critical temperature below which productivity drops. See Appendix Appendix A for more discussion.

The fourth term in Eq. (2) is fAPAR. “fAPAR” refers to the fraction of PAR that is absorbed by photosynthesizing parts (i.e., green leaves) of plants. PAR is the portion ( 50 %) of incoming solar radiation that is in wavelengths usable for photosynthesis. We approximate fAPAR as follows:

(4) f APAR = 1 - e - k veg Ω c LAI ,

where kveg is a light extinction coefficient (set to 1 for horizontal leaves), Ωc is the clumping index (or factor) and is set to 0.7, a near-mean value for natural land cover types (Pisek et al.2010; He et al.2012), and LAI is the leaf area index (see Sect. 2.2.8).

Equation (4) uses a simple Beer–Lambert approach, extended to canopies with leaves that are nonrandomly distributed in space (Nilson1971) (e.g., “clumped”). Here, we have used the common approximation (e.g., see Gower et al.1999) that fAPAR equals the fraction of intercepted PAR (fIPAR). Other assumptions and simplifications in Eq. (4) include an azimuthally symmetric leaf distribution (Gower et al.1999), leaf absorptivity of 1 with respect to PAR, and the neglect of the influence of non-photosynthesizing plant parts. A final and important assumption of horizontal leaves2 eliminates zenith angle dependency of the light extinction coefficient (Campbell and Norman1998), which makes the canopy gap fraction the same from any angle. This, in turn, eliminates the dependency of fIPAR on diffusive/direct radiative partitioning and on zenith angle, thus reducing fIPAR (and fAPAR) to the simple Beer–Lambert expression of Eq. (4), with a constant light extinction coefficient. Earlier versions of SimBA and the derivative model (Dewar1997) assume a constant value of 0.5 for kveg.

The simplifications made in the last paragraph allow us to interchange fAPAR and the leaf cover fraction, fleaf. Thus, we can rewrite Eq. (4) as follows:

(5) f leaf = 1 - e - k veg Ω c LAI ,

where fleaf is the areal fraction of view that is covered by photosynthesizing plant parts when looking directly down on the land surface. That is, fleaf=1- the gap fraction from the nadir.

Figure 3Zonal multiyear monthly means and zonal multiyear annual means of ET for SEDGES and the Mueller et al. (2013) reference dataset for the 1989–2005 period over non-glaciated land.


2.2.4 Water-limited gross primary productivity

Most land plants take up the CO2 they need for photosynthesis through tiny pores in their leaves called “stomata”. Water is also lost (transpired) through these same openings. Water loss in excess of water uptake can lead to cavitation, hydraulic system collapse (e.g., see discussion in Sperry et al.2002), and/or permanent reduction in photosynthetic capacity (Lawlor and Cornic2002). To prevent and mitigate such occurrences in both the present and future, most land plants adjust their stomatal openings in such a way as to balance the short- and long-term costs of transpiration with the photosynthetic gain from increased CO2 intake (Cowan and Farquhar1977; Medlyn et al.2011; Buckley and Schymanski2014; Prentice et al.2014). Thus, water limitation in most land plants is closely tied to CO2 limitation on photosynthesis. For this reason, this CO2-limited rate in SEDGES is referred to as the “water-limited rate of GPP”. This rate is given as follows:

(6) GPP W = ( co 2 conv ) 1 - c i c a c a ( f leaf ) ρ 1.6 r c + r a ,

where rc is the canopy resistance (see Sect. 2.2.6); ra is the aerodynamic resistance; ρpsfcRdTsfc is the surface air density, where psfc is the surface pressure and Rd is the gas constant for dry air; co2conv =4.15×10-7 kg C kg air−1 ppmv−1; ci represents a “bulk leaf” intercellular concentration of CO2 (ppmv); and cica is set to 0.80. This last value was lowered from the  0.86 value that was apparently used in Raupach (1998) because such a high value is less representative of typical daytime cica values (e.g., see Prentice et al.2014). SimBA uses a value of 0.7. The ramifications of using a fixed value of cica are discussed in Sect. 6.

Equation (6) is essentially the bulk aerodynamic formula for CO2 in Appendix B of Raupach (1998), except that we have multiplied their right-hand side by fleaf to generalize from their big-leaf approach to the simplified mosaic framework (as mentioned in Sect. 2.1), and we have multiplied by the ratio of the molar masses of carbon and CO2 to convert the flux into carbon units. In addition, we neglect the contribution of leaf mitochondrial respiration for simplicity. (It can be shown that including daytime leaf mitochondrial respiration in the formulation of GPPW multiplies its current form by a factor of 11-X, where X is the ratio of leaf mitochondrial respiration under light to gross photosynthesis.)

The above bulk aerodynamic formulation for GPPW in SEDGES is more realistic than the diffusive scheme used in earlier versions of SimBA and in the Dewar (1997) model. The latter scheme, at the canopy level, assumes that rc ≫ ra (which is noted, too, by Medlyn et al.2017). This condition is likely to be satisfied for aerodynamically rough canopies and when using daily-averaged environmental conditions, such as in Dewar (1997). However, early diagnostic tests of SEDGES coupled with PlaSim showed that the condition was often not satisfied, thus motivating our implementation of the bulk aerodynamic formulation for CO2 uptake. Implementing the new formulation led to higher daytime values of rc for moist areas (which is shown by Medlyn et al.2017, to follow from theory). These values were more in line with reported maximum surface conductances for various vegetation types (Kelliher et al.1995). Use of the bulk aerodynamic formulation for GPPW also unifies vegetation–atmosphere exchanges of CO2 and water under the same physical framework, which enables the process of transpiration to be properly tied into carbon uptake through interactive canopy conductance (see Sect. 2.2.6).

2.2.5 Evapotranspiration and transpiration

As was mentioned earlier in Sect. 2.1, the SEDGES framework for ET is very similar to the approaches taken by the SLAM-1T configuration of CHASM (Desborough1999) and by the original MOSES land surface scheme (Cox et al.1999). These approaches all expand on the SimBA framework (Appendix Appendix C) by separating evapotranspiration into soil and vegetative components and (in SEDGES and MOSES) by coupling photosynthesis and transpiration through interactive canopy conductance. We find that these changes greatly improve the simulation of ET (compare Fig. A2 with Fig. 3) when using the forcing setup described in Sect. 4.

By design, SEDGES does not directly calculate the evapotranspiration; rather, it provides a surface wetness factor to be used by the bulk aerodynamic formulation of the land surface scheme that is coupled to the atmospheric model, which performs the actual calculation of ET. An exception to this exporting of ET calculation to outside of SEDGES is the explicit calculation by SEDGES of evaporation from exposed soil when there is snow cover.

Although the (total) evapotranspiration is calculated outside of SEDGES, the way it is simulated is intimately tied into many parameterizations within the actual SEDGES model, which presuppose a specific ET framework (as discussed in Sects. 2.1 and 3). The ET framework is thus presented in this section, along with the SEDGES parameterizations that help comprise it. For this simplified mosaic approach (Sect. 2.1), we remind the reader that the aerodynamic conductance and saturation specific humidity of the air at the surface are homogeneous throughout the grid cell:3

(7) ET = ρ ρ w r a C w Δ q ,

where ET is the evapotranspiration in units of volumetric liquid water m3 m−2 s−1, ρw is the density of liquid water, Δqqsatsfc-q, qsatsfc is the temperature- (Tsfc) and pressure- (psfc) dependent saturation specific humidity of the air at the surface, q is the specific humidity of the air at the lowest atmospheric model level, and Cw is a surface “wetness” factor that reduces evapotranspiration from the potential rate by incorporating the effects of canopy resistance and soil resistance. Cw ranges from 0 to 1, where “1” gives the potential evapotranspiration (PET) rate:

(8) PET = ρ ρ w r a Δ q .

Although the above framework for ET (and PET) is required in order for SEDGES to operate as designed, there is no stipulation that the coupled land surface model calculate qsatsfc or ra in any preset way. That said, it is desirable for the calculation of ra to depend, in some way, on the surface roughness length that is output by SEDGES (Sect. 2.3.3).

Under snow-free conditions, SEDGES derives the key ET variable, the surface wetness factor, as follows:

(9) C w = f leaf 1 + r c g a + 1 - f leaf 1 + r ss g a ,

where ga is the aerodynamic conductance4 =1ra, and rss is the soil surface resistance, whose formulation is taken from the Joint UK Land Environment Simulator (JULES) model (Best et al.2011), except for the use of the lower minimum soil resistance from van de Griend and Owe (1994). This reduction was made to compensate for the absence of evaporation from ponded surface water in SEDGES. Soil surface resistance is formulated as follows:

(10) r ss = min r ss max , r ss min β ss ,

where rssmin is the minimum soil surface resistance of 10 s m−1 (van de Griend and Owe1994), rssmax is the maximum soil surface resistance (a very large number), and βss limits evaporation from the soil surface and is considered a “water stress factor”. This last factor is formulated according to Best et al. (2011), except that SEDGES uses the entire soil hydrological layer instead of just the top of several layers. The water stress factor for the soil surface in SEDGES is as follows:

(11) β ss = W frac 2 ,

where the soil wetness fraction, Wfrac, is given as follows:

(12) W frac = W soil W max ,

where Wmax is the biomass-dependent soil bucket depth and Wsoil is the water content (as depth) within the bucket. The SimBA model uses a general water stress factor that is a ramp function, reducing ET linearly from a soil wetness fraction of 0.25 to 0. In SimBA, Cw is equal to this water stress factor. In this way, no distinction is made between soil moisture's impact on soil evaporation and its impact on transpiration.

Examining Eq. (9), we see that Cw is comprised of a vegetation term and a bare-soil term (since 1−fleaf is the bare-soil fraction (neglecting non-photosynthesizing plant parts for simplicity)). The vegetation term involves loss of water only due to transpiration; canopy interception and evaporation are neglected for simplicity. Under snow-free conditions, ET is partitioned into transpiration and bare-soil evaporation, respectively, by replacing Cw in Eq. (7) with the vegetation term, fleaf1+rcga, or with the bare-soil term, 1-fleaf1+rssga. Doing so with the vegetation term yields the following equation for transpiration:

(13) T = ρ ρ w f leaf r a + r c Δ q .

From the above formulations, one can see that we weight the contributions to evapotranspiration from transpiration and bare-soil evaporation, respectively, by the fractional coverages of vegetation and bare soil.

In the presence of snow cover, the surface is treated as a mosaic of different types, with each type sharing the same aerodynamic conductance, soil hydrology, and surface temperature, as discussed before (Sect. 2.1). The parts of the grid cell that are snow-covered evaporate at the potential rate. The portion of the grid cell considered to be snow-free surface is comprised of vegetation parts that protrude above the snowpack and are snow-free, the exposed bare soil, and the exposed, low-lying leaf cover. In these snow-free portions of the grid cell, transpiration is 0 and bare-soil evaporation occurs (Eqs. 7, 9, and 10). The final equation for Cw in the presence of snow in the grid cell is as follows:


where ffor is the forest cover fraction (Eq. 22), fsnow flat is the fraction of the flat portion of the grid cell that is covered by snow (see Eq. 26 and nearby text), and fsnow for≡0.12 is the snow-covered fraction of the vegetation cover that protrudes above the snow pack (i.e., of the forest cover).5

From Eq. (14) and the relationship ET =CwPET (from Eqs. 7 and 8), one can deduce that bare-soil evaporation, Esoil, when there is snow cover is given as follows:

(15) E soil = ( 1 - f for ) ( 1 - f leaf ) ( 1 - f snow flat ) 1 + r ss g a PET .

In SimBA, Cw=1 when there is snow cover. The formulation was changed in SEDGES when it was discovered that it led to excessive ET when snow was present. The original MOSES and its progenitor UKMO (United Kingdom Met Office) scheme also set their equivalent Cw to 1 when there is snow cover, which also leads to excessive ET under those conditions (Essery1998; Cox et al.1999).

2.2.6 Canopy resistance

Canopy resistance, rc, is determined at each time step in accordance with the supply/demand principle in the original Dewar (1997) paper. That is, the plants attempt to adjust rc so that CO2 uptake exactly matches the light-limited rate of canopy photosynthesis, GPPL. (GPPL is the rate of gross carbon uptake in the absence of any CO2 limitation.) However, rc must also be sufficiently high to prevent transpiration from exceeding the maximum rate suppliable by the roots (the “supply rate”). If the rc value that would be needed to satisfy the CO2 demand caused the transpiration to exceed the supply rate, then water CO2 limitation occurs and rc takes on the value that would yield transpiration at the supply rate. In the derivation of the equation for canopy resistance, the unconstrained adjustment of rc is the first step. The second step is the calculation of the minimum possible value of rc that satisfies the supply constraint on transpiration. If needed, the unconstrained rc is raised to this minimum value. rc is also constrained to not exceed a maximum possible value, rcmax. Details on the mechanics of the formulation are given in Appendix Appendix B.

The original SimBA model does not explicitly consider canopy resistance. Although the incorporation of an interactive canopy resistance scheme within a coupled water CO2 exchange framework adds complexity to SEDGES, we feel that this is a critical feature for a vegetation model to have when it is used as part of a land surface model. In particular, significant climatic impacts of stomatal closure induced by increased CO2 are well-established in the literature and include reductions in ET, increased runoff, and increased surface temperature (Sellers et al.1996a; Betts et al.1997; Levis et al.2000; Betts et al.2007; Cao et al.2010).

2.2.7 Soil organic carbon

The overall modeling framework for soil organic carbon is the same as that of both the original SimBA (Kleidon2006b) and the land surface component (ENTS) of GENIE (Grid Enabled Integrated Earth System Model) (Williamson et al.2006). In this simple framework, soil organic carbon is treated as a single homogeneous reservoir (in which litter is included) that has a single, temperature-dependent residence time. Soil respiration (see below) depends only on temperature. The prognostic equation for soil organic carbon is given by

(16) dC soil d t = L - R soil ,

where Csoil is the soil organic carbon (kg C m−2), L is the litterfall (Eq. 1), and Rsoil is the rate of soil respiration (kg C m−2 s−1).

Soil respiration is modeled according to the following equation:

(17) R soil = c 8 C soil τ soil 1 + e c 9 T soil - 254.85 , T soil > 256.1 K , 0 , T soil 256.1 K ,

where Tsoil is the soil temperature at  0.20 m depth (i.e., the middle depth of the topsoil layer in PlaSim), τsoil is the residence time (42 years) of the soil organic carbon at 10 C, c81+e106283.15-254.85, and c9=106 K.

The above formulation is based on that of ENTS (Williamson et al.2006), except that we replace its temperature dependency (Lloyd and Taylor1994) with that from the RothC (Rothamsted Carbon Model) soil organic carbon model (Jenkinson et al.1990; Clark et al.2011), since the latter's temperature function reduces a negative soil organic carbon bias in the high Arctic. In addition, we use c8 as a normalizing constant to make the soil organic carbon residence time, CsoilRsoil, at 10 C the same as in the original SimBA formulation. Moreover, for strictly numerical purposes, 256.1 K is used as the cutoff temperature, below which soil respiration is set to 0. In SimBA, soil respiration is parameterized as follows: Q10Tsfc-1010Csoilτsoil. (Here, the Tsfc is surface temperature in C, and Q10 is set to 2). In the current version of SEDGES, it is tacitly assumed, for simplicity, that all respired soil organic carbon is emitted directly to the atmosphere.

2.2.8 Leaf area index and leaf cover fraction

LAI is typically defined as “the one-sided leaf area per unit ground area”, or else, for conifers, as the projected area of the needle leaves (Monson and Baldocchi2014, p. 246). In SEDGES (as well as in SimBA), LAI is based on a moist soil value that gets subsequently reduced for conditions of low soil moisture. We first describe the moist soil formulation.

For sufficiently moist soils, LAI is a simple function of biomass. In the version of SimBA that came out with version 15 of the Planet Simulator model (Lunkeit et al.2007), this moist soil LAI is a linear function of forest cover fraction (which, in turn, depends directly on biomass). This parameterization is discarded in the subsequent SimBA in version 16 of Planet Simulator, in favor of a simpler LAI functional dependency on biomass (Lunkeit et al.2011) to reduce the high number of multiple equilibria in the coupled climate-vegetation system that Dekker et al. (2010) found when using PlaSim version 15. This new parameterization is maintained in SEDGES, but its three parameters have updated values.

The equation for LAI in the absence of soil moisture stress is a monotonically increasing function of biomass with decreasing slope:

(18) LAI m = LAI min + 2 π ( LAI max - LAI min ) atan ( c 6 C veg ) ,

where LAImin and LAImax are the minimum and maximum LAIs (respectively) under moist soil conditions. LAImin is set to 0.05 and represents a seeding source with negligible mass. LAImax is set to 7.

While winter-deciduous phenology is not included in SEDGES, the model simulates drought-deciduous phenology using a crude parameterization from SimBA that came originally from the ECHAM3 model (Klimarechenzentrum1993). LAIm is converted into leaf cover fraction, fleafm, by substituting it into Eq. (5), giving the following:

(19) f leaf m = 1 - e - k veg Ω c LAI m ,

where fleafm is the leaf cover fraction under moist soil conditions.

Then, we set

(20) f leaf = min ( f leaf m , f leaf dry ) ,


(21) f leaf dry = 1 , W frac W frac crit , lai , W frac W frac crit , lai , 0 < W frac < W frac crit , lai , 0 , W frac = 0 .

Here, fleafdry represents the maximum leaf cover permitted by the soil wetness fraction; Wfraccrit,lai is the critical soil wetness fraction at which leaf cover begins to get restricted and is set to 0.05.

Although the final LAI is merely a diagnostic, it has a one-to-one relationship with fleaf because the fleaf in Eq. (19) is converted back into LAI by inverting the LAI–fleaf relationship in Eq. (5) as follows:

(22) LAI = - ln ( 1 - f leaf ) k veg Ω c .

2.2.9 Forest cover fraction

In SEDGES, forest cover fraction is defined as the fraction of the grid cell above ground that is covered by trees. However, woody shrubs are counted as partial trees. It has no seasonal phenology. Forest cover fraction helps to determine the amount of surface that is covered by snow for the calculation of evaporation (see Eq. 14). In SimBA, forest cover affects the albedo of snow-covered land; instead, SEDGES has a dependency on biomass (Eq. 25). As in the last version of SimBA (Lunkeit et al.2011), forest cover is parameterized to commence at a biomass of 1 kg C m−2. Such cover would represent woody shrubs. Only at a biomass threshold of 1.5 kg C m−2 is the woody cover considered tall enough so as to protrude above the winter snowpack and lower the albedo (Sect. 2.3.1). Forest cover fraction is formulated as follows:

(23) f for = max 0 , 1 - e - c 1 ( C veg - c 2 ) ,

where c1 is a shape parameter and c2 is the biomass threshold for forest cover.

The derivation of the c1 and c2 parameter values involved translating NPP from outside data into SEDGES biomass. If one assumes long timescales and steady-state conditions, then the translation is readily achieved by setting the left-hand side of Eq. (1) to 0 and solving for Cveg in terms of NPP. Doing so results in the following relationship: NPP =Cvegτveg. Using this relationship allows modeled annual NPP from Cramer et al. (1999) and McGuire et al. (1992) to be interpreted as SEDGES biomass equivalent. Then, comparing the real-world spatial distribution of boreal forest, boreal woodland, and tundra with the NPP data in those sources, including the NPP ranges for boreal forest, boreal woodland, and tundra ecosystems in McGuire et al. (1992), yields a rough threshold value of 1 kg C m−2 for forest cover commencement. The shape parameter, c1, was tuned so that forest cover fraction values in the Hagemann (2002) land cover dataset would match the NPP values in Cramer et al. (1999) for natural land cover types, while maintaining a fairly close likeness between the SEDGES-parameterized and Hagemann (2002) relationships between forest cover fraction and LAI. (This likeness was assessed by scatterplotting the forest cover fraction and maximum LAI for the moist ecosystems given in Hagemann (2002) along with SEDGES-formulated forest cover fractions and LAIs for a wide range of biomass values.) The original SimBA's (Kleidon2006b) dependency of forest cover fraction on biomass uses an S-shaped curve that is equivalent to ffor=max(0,atan(Cveg-3)-atan(-3))0.5π-atan(-3) and is similar to the formulation in the later SimBA version (Lunkeit et al.2011). The two previous SimBA formulations give higher forest cover fractions than SEDGES except at some low biomass values. Neither SimBA parameterization yields a good match with the aforementioned reference data in the absence of using additional fitting parameters.

2.3 Equations for remaining land surface variables

While the equations for the ecological variables were presented in Sect. 2.2, here, we provide the formulations for the land surface variables used by the PlaSim GCM as well as some purely diagnostic variables.

2.3.1 Surface albedo

A surface albedo with respect to broadband shortwave radiation is obtained by incorporating the effects of snow cover on the albedo for snow-free conditions. In obtaining the snow-free albedo (α0), we simplify by ignoring dependencies on solar zenith angle, diffuse and direct radiation partitioning, the spectral composition of the incident light, the effect of soil moisture on soil albedo, and any impact of leaf litter. Moreover, we neglect variations in leaf reflectivity that often occur between differing plant species and leaf development stages. The formulation of snow-free albedo in SEDGES is of the same form as that for SimBA, the ECHAM5 GCM (Rechid et al.2009) and the ENTS land surface scheme (Williamson et al.2006); it is as follows:

(24) α 0 = α veg f leaf + α soil ( 1 - f leaf ) ,

where αveg is the albedo of a completely leaf-covered surface and αsoil is the albedo of the bare soil.

Soil albedo decreases linearly with soil organic carbon to a minimum value that is reached at 9 kg C m−2 and stays constant with further carbon increase. The equation is as follows:

(25) α soil = α peat - α sand c 7 C soil + α sand , C soil c 7 , α peat , C soil > c 7 ,

where c7= 9 kg C m−2 and is the soil organic carbon level at which the soil albedo saturates, αsand is the soil albedo in the complete absence of soil organic carbon, and αpeat is the albedo of soil when saturated with soil organic carbon.

This soil albedo formulation is taken from ENTS (Williamson et al.2006), except for the following modifications: the peat albedo has been increased from 0.11 to 0.12, the sand albedo has been increased from 0.30 to 0.32, and c7 was decreased from 15 to 9 kg C m−2. The first change was made for the sake of model simplicity (since the albedo of a fully leaved surface is 0.12 in SEDGES); the increase in sand albedo was made to better match observed values in the Sahara and Arabian deserts (e.g., Knorr and Schnitzler2006). A saturation level for soil albedo as a function of soil column carbon content was estimated to be around 9 kg C m−2 by visually comparing the soil/litter surface albedo map from Houldcroft et al. (2009) with a global gridded dataset of soil organic carbon content (Wieder et al.2011). Using this estimate for saturation value in SEDGES reduces a positive albedo bias in the high Arctic in summer as compared to using the standard ENTS value of 15 kg C m−2.

Figure 4(a) SEDGES monthly mean climatologies of surface albedo for 2001–2010. (b) MODIS monthly mean climatologies of white-sky albedo for 2001–2010 (NASA LP DAAC2008b).


The soil albedo formulation in SEDGES is an advent to SimBA, which has always used a fixed soil albedo of 0.30. The new dynamic scheme was adopted, in part, due to the important role played by ground albedo in Sahelian/Saharan vegetation–precipitation feedbacks, as seen through its effect on low-frequency precipitation variability (Vamborg et al.2014) and increased greening in the mid-Holocene (Vamborg et al.2011). The new scheme also gives more realistic snow-free albedo values in the high latitudes and in the hot deserts (Fig. 4). While the above formulation for snow-free albedo, α0, neglects, for simplicity, the radiative impact of non-photosynthesizing plant parts (e.g., stems and branches), this impact becomes important (although implicit) in the snow-covered formulation.

A new snow albedo scheme for SEDGES replaces the version from SimBA that linearly combines the albedos from the forested and flat portions of the grid cell according to the forest cover fraction. The SimBA formulation did not simulate well the sharp, real-world transition in snow-covered albedo from tundra to boreal forest (Loranty et al.2014), which is why a new exponential decay scheme was developed.

In SEDGES, for snow depth > 0, the surface is treated as consisting of two components: (1) “flat” surface (consisting of exposed soil, prostrate vegetation, and dwarf shrubs) that can be partially or entirely buried by the snow cover; (2) “forest” that is covered by trees and shrubs of sufficient stature so as to protrude above the snowpack and mask it from the sun via their leaves, stems, and/or branches, regardless of accumulated snow depth on the ground. A nonlinear combination of the albedo of the flat portion and the snow-covered forest albedo yields the final surface albedo (α):

(26) α = α snow for + ( α snow flat - α snow for ) e - c 4 max ( 0 , C veg - c 5 ) ,

where c4 is a shape parameter, and c5 is the biomass threshold (1.5 kg C m−2) at which the woody vegetation is sufficiently tall as to begin to mask the snow.

The albedo of the flat portion of the grid cell, αsnow flat, is formulated according to the same temperature dependency (Roeckner et al.2003) and snow cover fraction (Roesch and Roeckner2006) schemes for non-forested areas in the ECHAM5 model, except that we maintain the melting snow albedo of 0.4 from PlaSim and ECHAM4 (Roesch et al.2001) and we neglect the impact of sloping terrain. We thus have the following:

(27) α snow flat = α 0 + ( α deep snow , flat - α 0 ) 0.95 f snow flat ,

where fsnow flat≡tanh(100 swe) is the fraction of the flat portion of the grid cell that is covered by snow (swe is the snow depth in liquid water equivalent (m3 m−2)),6 and αdeep snow,flat is the albedo of a deep and pure snowpack and is temperature dependent. This dependency has the form of a ramp function, with a maximum albedo, αmax deep snow flat≡0.8, for surface temperatures 5 C and a minimum albedo, αmin deep snow flat≡0.4, at melting point. These values are the same as in PlaSim and SimBA, except that they have maximum albedo attainment at 10 C.

The albedo of the forest component of the grid cell when there is snow cover is as follows:

(28) α snow for min ( α snow flat , α max snow for ) .

Here, αmax snow for is the maximum snow-covered albedo for forest and is set to 0.30, a midrange white-sky albedo value for the different forest types (Moody et al.2007). Restricting αsnow for to not exceed αsnow flat is done because the lower snow depth and warmer snow that reduce the albedo of the flat area also lower the albedo of the forested area to at least the level of the flat area.

2.3.2 Soil bucket depth

We first introduced the soil bucket depth, Wmax, in Sect. 2.2.5 in Eq. (12). SEDGES uses a simplified version of the parameterization in the JeDi (Jena Diversity) dynamic global vegetation model (Pavlick et al.2013). The JeDi model uses a pipe representation of the rooting system (Shinozaki et al.1964) in which a square root relationship emerges between Wmax and coarse root biomass as a result of an assumed constant density of fine roots in the rooting zone. To adapt this formulation to SEDGES, we further assume that coarse root biomass is a fixed fraction of total biomass and that the unit plant-available water capacity is spatially constant. Doing so gives the following relationship:

(29) W max = min W max min , c 12 C veg ,

where c12 is a tuning parameter set to 0.10 kg C1∕2. This value yields soil bucket depths that are a reasonably good fit with a plant-available water dataset based on optimal rooting depths (Kleidon and Heimann1998; Hall et al.2006; Kleidon2011). In that dataset, the soil bucket depths in the most sparsely vegetated regions range from around 0.05 m in the Canadian polar desert to < 0.003 m in the hyperarid Atacama and northeast Sahara deserts. SEDGES adopts a minimum soil bucket depth, Wmaxmin, of 0.05 m.

The versions of SimBA (Kleidon2006b; Lunkeit et al.2007, 2011) compute Wmax using the following formulation: Wmax=Z×Wmaxmax+(1-Z)Wmaxmin, where Z is the forest cover fraction formulation in the original SimBA (see Sect. 2.2.9) and Wmaxmax is some maximum value for Wmax. This Wmax formulation was updated in SEDGES because the new formulation has the advantage of being derived from first principles (Pavlick et al.2013).

2.3.3 Surface roughness

SEDGES has been designed for use in land surface schemes in which the surface roughness lengths for momentum and water vapor are the same (e.g., schemes that use parameterizations from Louis1979; Louis et al.1982). As such, SEDGES returns only a single value for surface roughness length, z0, for a grid cell. The (total) surface roughness is comprised of a surface roughness due to orography and a surface roughness due to vegetation:

(30) z 0 = z 0 veg 2 + z 0 oro 2 ,

where z0oro is the orographic surface roughness and z0veg is the roughness due to vegetation. z0oro is used by some GCMs (including Planet Simulator) to account for the enhancement of surface drag due to sub-grid-scale topographic variation (e.g., see Beljaars et al.2004). z0oro is set to 0 by default because small-scale orographic variation has little effect on land-to-atmosphere moisture transfer (Huntingford et al.1998; Blyth1999).

In the SimBA formulation, surface roughness due to vegetation is parameterized assuming that only forest cover increases it appreciably and that prostrate vegetative cover does not. As such, sparsely vegetated surfaces are assigned too high of a surface roughness. SimBA has minimum and maximum values of z0veg of 0.05 and 2 m, respectively (Kleidon2006b; Lunkeit et al.2007, 2011). SEDGES replaces SimBA's linear dependency of z0veg on forest cover fraction with the following logistic curve that depends on biomass rather than forest cover fraction:

(31) z 0 veg = c 17 1 + e - c 16 ( C veg - c 15 ) - z 0 const ,

where z0constc171+e-c16(-c15)-z0min is a constant (in meters) that forces z0veg to equal z0min at zero biomass. The new formulation for z0veg is tuned by adjusting c15 and c16 so as to roughly match reference values in the literature for different land cover types (especially those in Hagemann2002). z0min is set to a bare-soil value of 0.01 m (Oke1987). c17 represents the approximate maximum value of z0veg and is assigned a value of 2.5 m, which is representative of tropical rain forests (Sellers et al.1996b). z0veg values that lie between these two extremes have been constrained by using the NPP–biomass relationship described in Sect. 2.2.9 and by visually comparing reference NPP and GPP values (Cramer et al.1999; Jung et al.2011), biome/land cover distribution maps (Ramankutty and Foley1999; Olson et al.2001), and the aforementioned literature values of z0veg for those land cover types.

3 How to couple SEDGES

As stated in Sect. 2.1, SEDGES is a DGVM and is to be embedded within a land surface scheme of a climate or Earth system model. As such, the input and output variables for SEDGES are not the same as those for the land surface scheme that SEDGES forms a part of because SEDGES has its own set of exchanges with the rest of the land surface scheme in which it lies.

At minimum, the following time-varying input fields are required by SEDGES: aerodynamic conductance (ga), surface temperature (Tsfc), surface downwelling (broadband) shortwave radiation (SW), soil temperature at  0.20 m (Tsoil), (bulk aerodynamic) potential evapotranspiration (PET), soil moisture content (Wsoil), surface pressure (psurf), and snow depth water equivalent (swe). Most of these variables are either calculated by or made accessible to fully fledged land surface schemes.

Other variables that are used by SEDGES will usually be part of the land surface scheme and depend on additional input variables (not listed in the previous paragraph). This is also the case when using ERA-Interim (ECMWF Reanalysis – Interim) data to force SEDGES (Sect. 4). In addition, the following parameter values (see Table 2) are not specified within the actual SEDGES code, which means that they must be declared and assigned values either in outside modules (e.g., in the non-SEDGES part of the land surface scheme) or by modifying the current SEDGES code: ϵmax, atmospheric CO2 concentration (ca), the gas constant for dry air (Rd), maximum transpiration rate (trmax), and residence times for vegetative and soil organic carbon at 10 C (τveg and τsoil).

SEDGES returns the following output fields: live vegetative biomass (Cveg), soil organic carbon (Csoil), soil water-holding capacity (Wmax), surface wetness factor (Cw), transpiration (T) (as a diagnostic), evaporation from bare soil (only when snow is present), canopy resistance (rc), minimum canopy resistance (rcmin), surface albedo (α), total surface roughness (z0), forest cover fraction (ffor), leaf area index (LAI), leaf cover fraction (fleaf), net primary productivity (NPP), light-limited gross primary productivity (GPPL), temperature limitation function (f2(Tsfc)), water-limited gross primary productivity (GPPW), litterfall (L), and soil respiration (Rsoil).

In general, the land surface modeling framework that SEDGES presupposes must be reconciled with that of the land surface scheme of the Earth system or climate model that one is incorporating SEDGES into. In particular, in the absence of modification to either SEDGES or the original land surface scheme, the simulation of evapotranspiration according to the simplified mosaic and bulk aerodynamic formulations described in Sect. 2.1 and 2.2.5 is required. The requirement arises because the water and CO2 exchanges predicted through this framework are presupposed by the canopy resistance equations. The simplest and recommended scenario with regards to reconciling soil hydrology is that the land surface scheme use a single-layer bucket model because this is what SEDGES “sees” when calculating soil surface and canopy resistances. As an example, Eq. (31), used in our forcing of SEDGES with ERA-Interim data (Sect. 4), gives the standard hydrological formulation for the simple single-layer bucket model.

In the likely scenario that the original land surface scheme of the Earth system or climate model does not perfectly match the framework used by SEDGES, it should, in many cases, be easy to achieve consistency of frameworks by adapting SEDGES and/or the original land surface scheme. Such adaptation would be easiest for the case in which the land surface scheme uses a single or multiple tile/mosaic approach because this approach could be converted to the simplified mosaic framework of SEDGES by, at each time step, replacing tile surface temperatures and soil wetness fractions with their respective grid-cell-weighted averages and by computing an effective surface roughness for the grid cell (e.g., Mason1988) in the case of significant areal water body coverage.

Finally, many land surface schemes distinguish between a surface roughness length of water vapor and heat from a surface roughness length of momentum, which is done to correct the typical bulk aerodynamic formulation for its replacement of temperature at the height of the roughness length with surface temperature (e.g., see discussion and references in Chen et al.1997a). Although SEDGES only computes a roughness length of momentum, there are many ways to convert this roughness length to the roughness length for vapor/heat (e.g., see analyses and references in Chen et al.2010).

An advantage of SEDGES is that (at least in principle) it can be run at any horizontal and temporal resolution. However, we encourage users to run the model using a sub-daily time step so that the generally positive diurnal covariances of light, temperature, aerodynamic conductance, and surface-to-air specific humidity difference can be adequately captured. These covariances should increase ET, overall. On the other hand, if users choose a daily time step or longer, then they should also anticipate an increase in vegetative productivity in water-limited regions due to the concomitant reduction in water stress. This reduction in water stress can be offset by decreasing trmax, the maximum transpiration parameter (Appendix Appendix B).

4 Forcing SEDGES with ERA-Interim data

In this section, we describe the process by which we force SEDGES using an external forcing dataset, with the goal of evaluating the capability of SEDGES as part of a land surface model. Repeated 1981–2010 6-hourly data from ERA-Interim (Dee et al.2011) are used as the driver. The following required forcing variables (Sect. 3) were obtained from the reanalysis for direct input to SEDGES: surface (skin) temperature, surface downwelling shortwave radiation, temperature of the second-from-top soil layer, surface pressure, and snow depth water equivalent. Shortwave downwelling surface radiation is derived as follows: 3-hourly periods are isolated, and the two periods that straddle each given time (e.g., 12:00 Z straddled by the 09:00–12:00 Z and 12:00–15:00 Z periods) are averaged to give the flux at that time. A glacier mask is derived by considering any land point that has a daily mean snow depth always greater than zero to be glaciated. No vegetation occurs on glaciated grid points. Finally, all ERA-Interim data are interpolated to T62 resolution using climate data operators (CDOs) before being read into SEDGES.

In order to evaluate the performance of SEDGES, it is necessary to simulate soil hydrology in an interactive way rather than relying on soil moisture content from the external forcing dataset. (Early versions of SEDGES that were forced using reanalysis soil moisture resulted in unrealistic simulations of vegetation, especially in dry regions.) In order to do so, other land surface modeling components are combined with SEDGES, such that these components incorporate hydrologically relevant SEDGES output (namely, Wmax, Cw, bare-soil evaporation (in the presence of snow), and z0) and provide SEDGES with the variables needed: ga, PET, and Wsoil.

We simulate the aerodynamic conductance outside of SEDGES using these additional variables from ERA-Interim from the lowest atmospheric model level: the u and v wind components, the specific humidity, and temperature. The aerodynamic conductance and (bulk aerodynamic) PET are calculated according to the formulations of the Planet Simulator model (Lunkeit et al.2011). Soil moisture content is simulated outside of SEDGES as a prognostic variable using a simple hydrological scheme and three additional ERA-Interim variables: precipitation, snowfall, and snowmelt. As we recommend above in Sect. 3, we use a single-layer bucket model, with bucket depth determined by SEDGES. The bucket model is simple. As such, runoff is from the surface and only occurs when the soil moisture content exceeds the bucket depth. Infiltration of liquid water from rainfall and from snowmelt is unrestricted. We use the discretized version of the following formulation to simulate soil moisture:

(32) Δ W soil Δ t = P - S + M - ET soil ,

where P is precipitation, S is the rate of snowfall, and M is the rate of snowmelt (all in liquid water equivalent). Again, Wsoil cannot be greater than Wfrac, with any excess going into runoff. Here, ETsoil is the combined soil evaporation and transpiration (and thus not coming from the snow cover). Recall from Sect. 2.2.5 that when there is snow cover, soil evaporation (Esoil) is calculated by SEDGES. Under conditions of sufficient snow cover (see below for more), transpiration is 0 and ETsoil reduces from being the (total) ET to being Esoil.

Soil hydrology in the framework we have constructed for driving SEDGES requires careful treatment in the presence of non-zero snow depths. Such care is especially needed because of the Tiled ECMWF Scheme for Surface Exchanges over Land (TESSEL) that is used by ERA-Interim and the data interpolation to the coarse T62 grid. The underlying TESSEL tile scheme and the spatial interpolation imply that data for a given grid cell that are fed into SEDGES can be physically inconsistent in the SEDGES framework (Sect. 2.1). The inconsistency arises because a grid cell can have a non-zero snow depth and yet have a homogeneous (as required and seen by SEDGES) surface temperature that is well above 0 C. This situation, in the absence of any modification to the SEDGES code, gives rise to common occurrences in which there is substantial transpiration from snow-covered leaves and in which snow-covered vegetation and ground are evaporating as if their surface temperatures were well above freezing. These occurrences are not only unphysical; they give rise to higher simulated ET from a given snowy, above-freezing region than would occur over the same region if it were divided up into its original small, homogeneous but differing tiles.

A second issue that must be addressed when there is snow is the tendency for ERA-Interim to overpredict freezing rain at the expense of snowfall (Dutra et al.2011). Because we treat rain as entering the soil without interception by the snowpack (as is done in the TESSEL land surface model; Dutra et al.2010), the freezing-rain bias led our early simulations to have an unrealistically large recharge of soil water reservoirs in winter in many locations, even though surface temperatures were well below 0 C. The handling of the aforementioned snow-temperature inconsistency and the freezing-rain bias are discussed next.

In order to simply resolve the above issue of snow with above-freezing temperature, we separate possible conditions into two cases: (1) swe > snow thresh and (2) swe  snow thresh, where we define snow thresh as the swe threshold above which the surface is treated as snow-covered with respect to evaporation from and the physiology of vegetation and the threshold at or below which the surface is treated as snow-free. In case 1, we assume that transpiration is negligible (due to combined cold temperature and physical coverage of leaves by snow) and thus introduce a slight modification to Eq. (2), by setting f2(Tsfc) to 0 (its value at freezing), which, in turn, forces the productivity to 0. In case 2, we assume that evaporation from snow is negligible, and f2(Tsfc) is as normal. Finally, we address the excessive partitioning of precipitation into liquid form at subfreezing temperatures by further assuming for case 1 that soil moisture cannot increase unless there is some snowmelt during the same time step. Snowmelt indicates the presence of above-freezing temperatures and, thus, the presence of liquid water that does not freeze on contact with the surface and can thus infiltrate into the soil.

5 Model evaluation

In this section we evaluate SEDGES as part of the simple hydrological scheme described in Sect. 4, forced with ERA-Interim data. We emphasize the results from the simulation forced with historical CO2, using reference datasets of vegetative carbon, LAI, surface albedo, tree cover fraction, soil organic carbon, GPP, and evapotranspiration.

Four simulations are carried out: three equilibrium simulations and one transient simulation. In the three equilibrium simulations, atmospheric CO2 levels of 280, 360, and 560 ppm are held fixed through the simulations. The vegetation is spun up from a vegetative and soil organic carbon-free state. The spin-up includes an acceleration of the carbon cycle for the first 260 years (only). From simulation year 261 to year 1500, the carbon cycle is run normally. The model is very nearly at equilibrium by the last 30 years in all three runs, and these years are thus analyzed. The transient simulation begins with the end state of the 280 ppm equilibrium simulation and is subsequently forced using observed, transient atmospheric CO2 values from 1832 to 2010. The CO2 data are taken from the Mauna Loa dataset (Keeling et al.1976) for 1959 and onward and from 20-year smoothed ice core data (Etheridge et al.1998) for the prior years. Unless stated explicitly, the results and analyses always refer to the transient CO2 simulation. The model performs reasonably well, overall; GPP is simulated exceptionally well. In the evaluation process, we concentrate on the comparison of SEDGES with observation-based data products, but we sometimes supplement these comparisons with comparisons to land surface models and their performances, especially to the ENTS model (Williamson et al.2006) because it is of similar complexity to SEDGES.

5.1 Evaluation of productivity

GPP that is simulated by reanalysis-forced SEDGES is evaluated against two observation-based datasets: the Multi-Tree Ensemble (MTE) (Jung et al.2011) and CARBONES (CARBON fluxES and pools over Europe and the globe; In addition, we use the model intercomparisons of Sitch et al. (2008), Piao et al. (2013), and Anav et al. (2015) and modeling-based results of Hemming et al. (2013) and Holden et al. (2013) to assess the relative performance of SEDGES.

Figure 5Zonal multiyear annual mean of GPP for SEDGES and the two reference datasets, MTE and CARBONES, for 1990–2009 over non-glaciated land.


Mueller et al. (2013)(Jung et al.2011)Olson et al. (1985)(Wieder et al.2011)

Table 3RMSEs and correlations between multiyear annual means of SEDGES variables and those for reference datasets

Download Print Version | Download XLSX

For the 1982–2008 period, SEDGES simulates a global mean annual GPP of 126 Pg C yr−1 as compared to 120 Pg C yr−1 for the observation-based MTE dataset (Jung et al.2011). Zonal means of GPP, taken over the glacier-free land mask obtained from the ERA-Interim data and for 1990–2009, are shown for SEDGES, MTE, and the CARBONES data assimilation dataset in Fig. 5. SEDGES captures the large-scale patterns very well. Compared to the two reference datasets, SEDGES has positive productivity biases in most of the Southern Hemisphere and around 35 N, and it has negative biases near the equator and in the high northern latitudes. Moving to the full spatial field, we can see in Fig. 6 that the annual mean GPP, for 1990–2009, in SEDGES and MTE compares very well, overall, with those of MTE and CARBONES. Regionally, as compared to the two reference datasets, SEDGES simulates too high a GPP in western Argentina, in nondesertic tropical Africa, around the Korean Peninsula, and on the southwestern and eastern Australian coasts and too little GPP in almost all the equatorial tropics and Amazonia (see Sect. 5.6 for discussion), the Pacific northwest of the United States, northwestern Europe, and parts of north-central and far eastern Siberia. In Siberia, SEDGES simulates too low a GPP in a dry region northeast of the Kolyma Mountains. This underestimate may be partially due to SEDGES's neglect of leaf mitochondrial respiration in the light, which has been found to be a large fraction (up to  0.3) of gross photosynthesis in Arctic tundra plants (McLaughlin et al.2014; Heskel et al.2014). Indeed, a fraction of 0.3 would, under purely water-limited conditions, increase GPP by 43 %.

The spatial correlations between SEDGES and the reference datasets are as high as can be expected (Table 3). Spatial correlations between multiyear annual means of GPP from MTE and three offline land surface models (ORCHIDEE, ORganizing Carbon and Hydrology in Dynamic EcosystEms (Krinner et al., 2005); JULES; and CLM4CN, Community Land Model with coupled Carbon and Nitrogen cycles, Lee et al.2013) range from  0.87 to 0.95 (Anav et al.2015), whereas it is 0.92 for SEDGES. The same correlations between the three models and CARBONES ranges from  0.83 to 0.87 (Anav et al.2015), whereas it is 0.86 for SEDGES.

Figure 6Multiyear annual mean of GPP for SEDGES and the two reference datasets, MTE and CARBONES, for 1990–2009 over non-glaciated land.


Figure 7Zonal multiyear monthly means of GPP for SEDGES and the MTE reference dataset for 1982–2010 over non-glaciated land.


SEDGES simulates the multiyear monthly means of GPP well, as shown by comparison with those for the MTE dataset in Fig. 7 for different latitudinal bands. From this figure, we can see that SEDGES captures the seasonal progressions of GPP well. A noticeable departure of SEDGES from MTE, however, is a slight phase shift in the Northern Hemisphere extratropics, with a generally earlier spring productivity increase in SEDGES than in MTE. We suspect that the phase shift is due to the absence of both winter-deciduous leaf phenology and physiological constraints on the speed of acclimation to temperature (discussed more in Appendix Appendix A). Moving to the full spatial field, Fig. 8 shows that the correlations in the seasonal cycle between SEDGES and MTE are generally high in areas with a strong seasonality of productivity such as the mid- to high latitudes, India, and the dry tropics. Temporal correlations are weaker, but still generally positive, in semiarid regions and in Amazonia. Correlations are often negative in regions with low absolute seasonal variation in GPP such as equatorial Africa, the northern hemispheric deserts, and southern Australia. SEDGES's problems in the dry regions may be due to the simple single-layer bucket soil hydrology that was used. The negative correlation in equatorial Africa is due to insufficient moisture stress in that region in ERA-Interim-forced SEDGES. The lack of moisture stress in this region causes dry-season GPP to exceed that of the wet season due to the reduced cloud cover, as is seen in the wetter parts of the Amazon (Wu et al.2016). However, the real African equatorial forest experiences a drop in GPP during the dry season relative to the wet season (Guan et al.2015). Thus, the real-world behavior in equatorial Africa is anticorrelated with that in the model. In equatorial Africa, the lack of moisture stress in ERA-Interim-forced SEDGES is likely attributable to a pronounced positive precipitation bias in the ERA-Interim forcing data in that region (Lorenz and Kunstmann2012; Dolinar et al.2016).

Figure 8Pearson correlation coefficients of the seasonal cycle (i.e., of the multiyear monthly means) between SEDGES and the MTE reference dataset for the 1986–2005 period.


The interannual variability in global GPP for 1990–2009 in SEDGES is 1.79 Pg C yr−1, whereas it is 2.50 Pg C yr−1 for the CARBONES dataset. Anav et al. (2015) report values of 3.23, 4.4, and 2.87 Pg C yr−1 for the offline-driven land surface models ORCHIDEE, JULES, and CLM4CN, respectively. The interannual variability in GPP at each grid point is shown in Fig. 9 for SEDGES and CARBONES. SEDGES successfully captures the spatial pattern of the CARBONES reference dataset, although its magnitudes are smaller, overall. These smaller magnitudes are probably attributable to a number of factors. Figure 9 shows that SEDGES most severely underestimates interannual GPP variability in cropland regions (cf. Ramankutty et al.2008), which can differ greatly in form and function from the natural vegetation simulated by SEDGES. Secondly, the underestimate of mean annual GPP by SEDGES in some regions such as northern Asia (Fig. 6) likely contributes to the model's underestimate of GPP variability in those same regions. Finally, SEDGES applies the same temperature limitation function to productivity (Appendix Appendix A) in all situations. In so doing, the model underestimates physiological constraints on the speed with which vegetation can adapt to new thermal regimes. Although in real life these constraints apply mostly on daily timescales, interannual variation in the frequency of extreme temperature events probably has greater impact on annual GPP in the real world than in SEDGES.

Figure 9Interannual variability in GPP for SEDGES and the CARBONES reference dataset for the 1990–2009 period.


The trend in annual GPP for 1990–2009 is captured extremely well by SEDGES. The CARBONES reference dataset yields a global annual GPP trend of 0.086 Pg C yr−2, whereas SEDGES gives a value of 0.080 Pg C yr−2. The spatial pattern of the trend is seen in Fig. 10. SEDGES captures the spatial pattern of GPP trend in the CARBONES dataset better than ORCHIDEE, JULES, and CLM4CN, as reported by Anav et al. (2015). SEDGES generally outperforms the other models in Africa, Australia, the semiarid areas of eastern Europe, and central and eastern Asia. Notable areas where SEDGES (but not the other models) deviates from CARBONES include the Amazon and parts of northwestern North American, in which SEDGES simulates more negative trends in GPP than the reference dataset. The negative GPP trend in SEDGES in the Amazon is very possibly due, at least in part, to its simplifications affecting the CO2 fertilization of productivity, as is explained next.

The SEDGES treatment of CO2 fertilization is seen via Eqs. (3) and (6) and notably lacks a direct dependency on temperature. However, Long (1991), using the widely used Farquhar model of photosynthesis (Farquhar et al.1980), finds that, for a given cica, the proportional increase in leaf photosynthesis with elevated CO2 increases substantially with temperature under both Rubisco- and light-limited conditions. In line with that result, a modeling study on the scale of a global grid (Hickler et al.2008) which uses a photosynthesis model based on that of Farquhar, finds that the fertilization response of NPP to elevated CO2 increases substantially, too, when moving from boreal forest to temperate forest to tropical forest. Hence, if SEDGES were to use a more complex parameterization of photosynthesis that is based on the Farquhar model, one would expect the Amazon region to show a more positive GPP trend in Fig. 10 and thus be closer to the CARBONES data.

In spite of the simplifications affecting CO2 fertilization in SEDGES, its global productivity response to CO2 forcing is comparable to that of other models of the land surface and vegetation. The results of transient experiments in which 10 such models are forced offline with repeated early 20th century climate and with increasing CO2 values show relative increases in global NPP per ppm CO2 of 0.05 to 0.20 % (with mean of 0.16 %) for the 1980–2009 period (Piao et al.2013); it is 0.102 % for SEDGES. Similar transient simulations with five of the aforementioned models (but also including projected CO2 levels) (Sitch et al.2008) give global GPP increases7 from 280 to 360 ppm CO2 that range from approximately 14 to 19 %; they give global GPP increases from 360 to 560 ppm CO2 that range from approximately 23 to 36 %. In comparison, the transient run of SEDGES gives an increase in global GPP from 280 to 360 ppm CO2 of 20 % and an at-equilibrium global GPP increase from 360 to 560 ppm CO2 of 30 %. Finally, Hemming et al. (2013), approximating the purely physiological effect8 of CO2 on productivity using a GCM, finds an increase in equilibrium GPP of 75 % when doubling CO2 from near-preindustrial values. A most likely increase in global NPP of 40 to 60 % with CO2 doubling from preindustrial levels is found by Holden et al. (2013), using a combination of global atmospheric CO2 rise and estimated land use changes since the preindustrial period. Doubling preindustrial CO2 (280 to 560 ppm) in equilibrium simulations in SEDGES induces an NPP increase of 55 %.

Figure 10Linear regression trend in annual GPP for SEDGES and the CARBONES reference dataset for 1990–2009.


5.2 Evaluation of vegetative carbon

Vegetative carbon is simulated reasonably well by SEDGES. Unless stated otherwise, we compare the means of the last 30 years of the CO2-varying run with an older reference dataset (Olson et al.1985). The Olson dataset was chosen because its data sources reflect the potential natural vegetation better than the sources of the more recent NDP-017b (NDP – numeric data package) dataset (Gibbs2006). SEDGES has an equilibrium global vegetative carbon of 530 Pg C under preindustrial levels of CO2 and a mean of 615 Pg C for 1981–2010 of the transient CO2 run. In comparison, the Olson et al. (1985) dataset and the more recent Gibbs (2006) dataset give global vegetative carbon of 451 and 560 Pg C, respectively (Jiang et al.2015, citing Gibbs2006).

Figure 11Vegetative carbon for SEDGES: mean over 1981–2010 in the transient CO2 simulation; vegetative carbon from the Olson et al. (1985) reference dataset, which represents pre-Iron Age vegetation save for the most extreme anthropogenic land cover changes.


Figure 11 shows the distributions of mean vegetative carbon for SEDGES (1981–2010) and that for the Olson et al. (1985) dataset. The spatial pattern of Olson vegetative carbon is captured well by SEDGES (Table 3: R=0.57, R2=0.33). The RMSE (root mean squared error) between SEDGES and the Olson datasets is 3.92 kg C m−2. (We calculate RMSE in this study by taking the square root of an area-weighted mean of the squared differences between SEDGES and the reference at each grid point.) Among the discrepancies, two tendencies prevail: the almost-direct relationship between GPP biases and vegetative carbon biases and an overall tendency to overestimate vegetative carbon in the tropics and to underestimate it in the mid- to high latitudes. SEDGES has large positive vegetative carbon biases in most of the tropics (except the western Amazon), southeast North America, China, and northwestern Europe. One should note that many of these biases reflect land cover change from forests to croplands (Ramankutty and Foley1999). In addition, positive biases in the wet-and-dry and semiarid tropics are partially attributable to the absence of a fire module in SEDGES; fire can have a large impact on vegetative cover in these regions (Bond et al.2005). SEDGES has large negative biases in western Canada and central Siberia. The subarctic is generally negatively biased. These vegetative carbon biases are due to a combination of the aforementioned GPP biases and the globally uniform residence time for vegetative carbon in SEDGES (Table 2). The SEDGES 10-year residence time for vegetative carbon matches the observation-based estimate CvegNPP given by Jiang et al. (2015).

Compared to SEDGES, ENTS (Williamson et al.2006) simulates vegetative carbon slightly closer to the Olson reference, overall. ENTS has a strong negative biomass bias in eastern Australia, and it greatly overpredicts biomass in eastern Brazil and northeastern Canada; SEDGES predicts these regions fairly well. On the other hand, ENTS also predicts more biomass in central Siberia (i.e., has a weaker negative bias), simulates western Canada well, has a less severe negative biomass in far eastern Siberia and Alaska, and has a less severe positive biomass bias in the tropics.

Figure 12Soil organic carbon for SEDGES: mean over 1981–2010 of the transient CO2 simulation; soil organic carbon from the HWSD reference dataset (Wieder et al.2011). The HWSD dataset has values for the top meter of the profile, only.


5.3 Evaluation of soil organic carbon

At present, land surface models have great difficulty in simulating soil organic carbon well. Tian et al. (2015), using 10 offline-forced terrestrial biosphere models at 0.5 resolution, find highly diverging model estimates of soil organic carbon dynamics as well as systematic biases due to absent processes affecting high-latitude soil organic carbon stocks. The models in Tian et al. (2015) have Pearson correlation coefficients under 0.4 with respect to the Harmonized World Soil Database (HWSD) reference dataset (Wieder et al.2011). SEDGES performs better with a 0.58 correlation (at the coarser T62 resolution), although it should be kept in mind that model performance in simulating soil organic carbon has been found to improve dramatically when aggregating on large spatial scales (Todd-Brown et al.2013). As listed in Table 3, SEDGES has a 7.9 kg C m−2 RMSE with respect to the HWSD data.

Figure 12 shows the full spatial distribution of soil organic carbon for SEDGES and the HWSD reference dataset. Compared to HWSD, SEDGES simulates less soil organic carbon in semiarid and arid areas, in the more northern parts of the Arctic, and generally more soil organic carbon in the remaining areas. Compared to the ENTS model (Williamson et al.2006), SEDGES gives a broadly similar simulation of soil organic carbon, as expected since their soil respiration formulations are almost identical. Some notable differences between the two are that SEDGES generally simulates more soil organic carbon than ENTS; SEDGES predicts too much soil organic carbon in southeastern North America, southeastern Asia, and the tropics, whereas ENTS generally does not; in high northern latitudes, ENTS (only) simulates a pattern of very high values in far eastern Siberia, Alaska/Yukon, and northeastern Canada with lower values in between, but this pattern is not seen in the HWSD dataset.

Taking the globe as a whole, SEDGES has a residence time for soil organic carbon of 25.7 years and a time-averaged soil organic carbon of 1516 Pg C for 1981–2010. Respective estimates for these values based on HWSD data are 24 years for residence time (Todd-Brown et al.2013) and range from 891 to 1657 Pg C for total terrestrial soil organic carbon storage (Tian et al.2015). Köchy et al. (2015) estimate total terrestrial soil organic carbon at a higher value of  3000 Pg C because they include estimates of carbon in the deeper soil layers that HWSD does not.

Figure 13Maximum LAI of the multiyear monthly means for SEDGES and for the BNU MODIS-based reference dataset for the 2001–2010 period.


5.4 Evaluation of LAI

We use reprocessed Moderate Resolution Imaging Spectroradiometer (MODIS) data from 2001 to 2010 from Bejing Normal University (BNU) (Yuan et al.2011) as the reference set. We interpolated these data to T62 resolution and derived multiyear monthly means from them. Outside of the tropics, we compare only the maximum of the multiyear monthly means of LAI between SEDGES and the reference data because of the absence of cold deciduous phenology in SEDGES and because serious deficiencies with MODIS in capturing the seasonal cycle of LAI in boreal coniferous forests have been found (Serbin et al.2013).

Figure 13 shows the maximum multiyear monthly mean LAI in SEDGES and the BNU MODIS-based reference for the years 2001 to 2010. SEDGES has strong positive biases in southeastern and extreme northern South America, in parts of Africa, in parts of southern and eastern Asia, in northern coastal Australia, and in parts of Europe. Many of these LAI biases lie where positive GPP biases occur and/or where croplands have replaced natural forests. Strong negative biases occur in parts of northern Asia. Of these, the biases east and northeast of Lake Baikal occur in the absence of negative GPP biases. Overall, the spatial patterns of maximum LAI in SEDGES and the MODIS-based reference have a correlation of 0.793.

Figure 14Temporal correlation at each grid point in the tropics (−20 to 20) between the multiyear monthly means of LAI for 2001–2010 of SEDGES and of the BNU MODIS-based reference dataset. Some desert grid points do not have values because they have LAIs of 0 in one or both sets of data for all 12 months.


Figure 14 shows the temporal correlation in the tropics between the multiyear monthly means of SEDGES and the BNU MODIS-based reference (i.e., the correlations of the two seasonal cycles). In spite of the crude scheme for drought-deciduous phenology described in Sect. 2.2.8, SEDGES does a reasonable job in simulating the seasonal cycle of LAI in much of the tropics. Negative correlations (indicative of poor model performance) occur near the equator, especially in the Amazon region, where LAI varies modestly, seasonally, in the real world (e.g., Wu et al.2016). The wet parts of the Amazon generally experience high water availability, even in the dry season. As a result, LAI there is decoupled from soil-moisture-induced water stress (in contrast with the drier portions of the Amazon) and LAI, generally, increases during the dry season and decreases during the wet season (Wu et al.2016). SEDGES cannot capture this behavior.

5.5 Evaluation of surface albedo

Figure 4 shows monthly mean climatologies of SEDGES surface albedo and MODIS white-sky albedo (NASA LP DAAC2008b)9 for the months of January, April, July, and October for the years 2001 to 2010. From these figures, we can see that SEDGES captures the overall spatial albedo pattern fairly well. In fact, spatial correlations between SEDGES and the reference range from a low of 0.77 in July and August to a high of 0.91 in April. The high correlation in boreal spring indicates relative strength in simulating the winter melt season at mid- to high northern latitudes. The relative weakness of SEDGES in Northern Hemispheric summer is due to the absence of a strong snow signature on the land and relative weakness of the model in correctly capturing the second- and third-order determinants of (non-glacier) terrestrial albedo, since these determinants come to the fore during boreal summer. Three such determinants are now discussed.

The most significant model deficiency in albedo simulation that gains prominence during boreal summer is the simulation of bare-soil surface albedos. Recall from Sect. 2.3.1 that bare surface albedo in SEDGES is determined by soil organic carbon content and by the albedo of soil in the absence of soil organic carbon (αsand), which is assumed constant. In the real world, geological processes, the quantity and properties of dead biomass litter, and the albedo of the underlying bedrock greatly impact αsand as well (Knorr and Schnitzler2006; Vamborg et al.2011). The absence of these real-world processes in SEDGES gives rise to an underestimate of albedo in the Arabian and the Saharan deserts and an overestimate of albedo in the snow-free season in most of the world's remaining deserts (Fig. 4). MODIS values of snow-free albedo (NASA LP DAAC2008a) in the polar desert (not shown) vary generally from 0.12 to 0.22 in July. The negative soil organic carbon bias in SEDGES in the high Arctic (Sect. 5.3) causes the snow-free albedo to be too high in polar deserts of the high Arctic.

A second major model deficiency in boreal summer albedo simulation stems from SEDGES's inability to distinguish between low-albedo ( 0.09) evergreen coniferous forests and high-albedo ( 0.16) mid- or high-latitude broadleaf deciduous forests. This deficiency is seen most clearly in July (Fig. 4) in eastern North America, in which SEDGES fails to capture the south-to-north decrease in albedo from the temperate to subarctic regions.

The third significant model deficiency in the simulation of albedo during boreal summer is a negative bias in areas that are grassland or savanna in the real world, particularly those found in Mongolia, the US Great Plains, and tropical Africa. In the latter two regions, this negative albedo bias is associated with a positive bias in vegetative carbon (Sect. 5.2) and LAI (Fig. 13).

Outside of boreal summer, positive biases in vegetative carbon in the mid- to high latitudes also cause a lowering of SEDGES January albedo as compared to MODIS (Fig. 4). The most affected regions are in central Asia and the North American prairies and High Plains.10 In March (not shown), there is a positive albedo bias in northern Siberia (specifically, in transition areas between boreal forest and tundra) whose source cause is a negative GPP bias seen in Fig. 6. This positive bias carries over, albeit less strongly, into April (Fig. 4).

Finally, it must be noted that our albedo evaluations of SEDGES have been with respect to white-sky albedo and not actual albedo.

Figure 15Multiyear annual means of ET for SEDGES and the Mueller et al. (2013) reference dataset for the 1989–2005 period over non-glaciated land.


5.6 Evaluation of evapotranspiration and runoff

Evapotranspiration from SEDGES is compared to that of a multisource compilation reference dataset, Mueller et al. (2013), for the period from 1989 to 2005. Simulated evapotranspiration by SEDGES qualitatively follows that of the reference dataset, successfully capturing the first-order seasonal zonal mean pattern and annual pattern (Figs. 3 and 15). Some salient deficiencies exist, however.

Globally, the annual means in Figs. 3 and 15 show that ET is excessive, overall. Future versions of SEDGES will improve this bias, possibly through the inclusion of a more realistic treatment of the ratio of intercellular to atmospheric CO2 cica instead of the currently fixed value of 0.80, which is probably too high, in general. The important role of cica in transpiration is discussed more below and in Sect. 6. On the other hand, the percentage of global evapotranspiration coming from transpiration is 64 % in SEDGES, which compares very well with a literature-based estimate of 61 % (with 15 % standard deviation), derived from 81 studies (Schlesinger and Jasechko2014).

Figure 16July mean ET for SEDGES and the Mueller et al. (2013) reference dataset for the 1989–2005 period over non-glaciated land.


Regionally, the most severe discrepancy from the Mueller et al. (2013) reference dataset is the excessive ET simulated by SEDGES in the mid- to high northern latitudes, which is most pronounced in the transition season from late winter to early spring (Fig. 3). This ET bias occurs because the new parameterization (see Sect. 2.2.5) to reduce excessive sublimation from snow-covered surfaces in the original SimBA parameterization is only moderately effective. This deficiency is possibly the most crucial one to improve in future model versions, possibly by implementing a scheme to decrease surface roughness from snow-covered terrain (e.g., Cox et al.1999). Another major contributor to the mid- to high northern-latitude positive ET bias is the absence of winter-deciduous leaf phenology and the assumed instantaneous temperature acclimation in SEDGES (Appendix Appendix A), which together increase spring GPP (Sect. 5.1) and thus, concomitantly, increase transpiration. Additionally, Fig. 16, which shows July ET for SEDGES and the Mueller et al. (2013) dataset, suggests that snow-free bare-soil evaporation from moist areas is excessive, as indicated by the relatively large positive ET biases in the high northern latitudes, which have a large contribution to total ET from soil evaporation due to their moist soils and low LAI (Fig. 13) (and hence low leaf cover fraction).

Figure 17Zonal monthly mean and zonal annual mean climatologies of runoff for SEDGES and the UNH-GRDC reference dataset (Fekete et al.2002) over non-glaciated land.


In addition to the northern mid- to high-latitude positive bias, SEDGES suffers from a negative ET bias in the equatorial region (Fig. 3), which has at least three likely causes. A positive bias in annual cloud cover fraction in the ERA-Interim data (Dolinar et al.2016) results in a negative bias in the incident surface shortwave radiation field (Zhang et al.2016) on the order of  5 to 30 % throughout most of the equatorial tropics for at least 6 months of the year. This bias may result in a relative reduction in GPP of similar magnitude in SEDGES, which would help to explain the negative GPP bias that SEDGES simulates here (seen in Figs. 5 and 6), while also explaining a reduction in transpiration of comparable (but lesser) magnitude. The latter would result from the increased canopy resistance (Eq. 13) that would ensue from the decreased demand for CO2 uptake due to the lower light. The lack of an increase in light use efficiency under cloudy skies in SEDGES further exacerbates the negative GPP bias caused by the surface shortwave reduction (see Sect. 6 for more discussion). The second reason for low ET near the equator is the complete absence of canopy interception loss in SEDGES, since this source of surface evaporation has been found to comprise around 20 % of total ET in tropical rain forests (Shuttleworth1988; Da Rocha et al.2004; Czikowsky and Fitzjarrald2009; Miralles et al.2011).11 Finally, the low vapor pressure deficits in the wet tropics do not activate a mechanism in SEDGES that, in other models, would increase transpiration. Less transpiration is simulated by SEDGES in this region because SEDGES's simplifying assumption of a fixed cica promotes underestimation (overestimation) of transpiration when vapor pressure deficits are low (high), due to the increase in cica with decreasing vapor pressure deficit (Medlyn et al.2011, 2012; Prentice et al.2014). It then follows from Fick's law of diffusion (on the plant scale) that higher cica decreases the leaf-to-air CO2 gradient, thus requiring that stomatal conductance increase to ensure sufficient CO2 inflow for photosynthesis. The higher stomatal conductance increases transpiration. A similar process happens on the canopy scale.

Conspicuously high ET can be seen toward the end of the wet season at both the northern and southern edges of the tropics in Fig. 3. A very similar zonal mean temporal pattern is found in transpiration in the STEAM (Simple Terrestrial Evaporation to Atmosphere Model) land surface model (Fig. 6 of Wang-Erlandsson et al.2014). Because this model was evaluated using ERA-Interim forcing data, the high ET anomalies that we see in Fig. 3 may be attributable to this particular forcing; the ET anomalies in SEDGES may also be attributable to the lack of canopy interception loss, which declines in tandem with precipitation toward the end of the wet season in Wang-Erlandsson et al. (2014), thus reducing the magnitude of the ET aberrations.

On a smaller regional scale, areas of complex terrain and topography (the Andes and the Himalayas) have very strong positive annual mean ET biases (Fig. 15). We believe that this is an interpolation artifact stemming from cold, snow-covered grid cells lying adjacent to warm, snow-free grid cells. With the spatial interpolation of ERA-Interim to the coarse T62 resolution, such cases often result from deep snow being distributed over the entire area and having a relatively warm temperature due to the influence of the low-lying grid cells at the original high resolution. This setup results in strong evaporation.

With respect to other models of similar complexity, SEDGES simulates ET much better than the original SimBA (see Appendix Appendix C). However, at least when looking at the climatological annual mean, ET is considerably better-simulated by ENTS (Williamson et al.2006) than by SEDGES. In general, ENTS has a negative ET bias in dry regions and a positive bias in moist regions, especially in parts of the tropics. However, these deviations are generally no greater than the biases that SEDGES has.

Simulated runoff is compared to the composite gridded dataset from the University of New Hampshire (UNH) and World Meteorological Organization's Global Runoff Data Centre (GRDC) (Fekete et al.2002). SEDGES simulates runoff well except for excesses in the equatorial tropics and in the midlatitudes of the Southern Hemisphere (Fig. 17), especially in equatorial Africa and in the Andes (Fig. 18). Excessive runoff in the equatorial tropics is due to the aforementioned negative ET bias, but it is also attributable to a general overestimate of precipitation in these latitudes by ERA-Interim, especially in equatorial Africa (Lorenz and Kunstmann2012; Dolinar et al.2016). Global runoff for SEDGES, excluding areas of land ice, is about 44 ×1015 kg yr−1, which is slightly higher than the 38 ×1015 kg yr−1 for the UNH/GRDC dataset and the global (excluding Antarctica) runoff estimate of 37 ×1015 kg yr−1 by Dai and Trenberth (2002).

It should be noted that the parameterization for soil water-holding capacity (Sect. 2.3.2) gives generally much lower values (not shown) in the wet and dry tropics than in the reference dataset (Kleidon and Heimann1998; Hall et al.2006; Kleidon2011). Such lower values of soil water-holding capacity should lead one to expect a positive runoff bias in these areas in SEDGES. This is not the case, however, when we compare those regions of the wet and dry tropics where the ERA-Interim precipitation forcing has insignificant bias (and where biases in SEDGES runoff are thus unlikely to be attributable to precipitation deficiencies in the forcing). The lack of runoff bias in these regions suggests that the parameterization of soil water-holding capacity in SEDGES (Eq. 28) works well in tandem with its ET scheme and the simple soil bucket hydrology (Eq. 31).

Figure 18Climatological annual mean of runoff for SEDGES and the UNH-GRDC reference dataset (Fekete et al.2002) over non-glaciated land.


Figure 19Mean forest cover fraction for SEDGES and tree cover fraction from the reference ISLSCP II dataset (DeFries and Hansen2009).


Finally, while it might seem easy to attribute deficiencies in the SEDGES hydrology to our use of the single-layer soil bucket, this type of model has been found to perform well under many conditions: wet soils, when precipitation events tend to be large relative to the bucket depth, and under strong root compensation or hydraulic redistribution of soil moisture (Guswa et al.2002; Guswa2005). Even when the aforementioned conditions do not hold, distortions stemming from the use of the single-layer bucket model are comparable to those coming from other aspects of SEDGES hydrology.

5.7 Evaluation of forest cover fraction

The overall global pattern of tree cover fraction of the International Satellite Land Surface Climatology Project (ISLSCP II) reference dataset (DeFries and Hansen2009) is captured reasonably well by SEDGES, as seen from Fig. 19. However, SEDGES has almost double the amount of arboreal cover. In Europe, China, and India, this is at least partially due to anthropogenic land cover change of natural forests and woodland to croplands (Ramankutty and Foley1999). Aside from this deforestation, the model's lack of competing plant functional types (especially boreal trees, warm region trees, and grasses) in conjunction with the constant residence time of 10 years for all vegetative carbon, leads to an overestimation of arboreal vegetation by SEDGES in South America, Africa, and North America, as well as a slight underestimation in boreal forest regions. A noted limitation of the reference dataset is its underestimation of real-world tree cover in areas with the highest cover fraction. Thus, SEDGES's overestimates in these regions may partially be due to this phenomenon, as well. Finally, a less inclusive definition for arboreal cover in the reference dataset than in SEDGES contributes to the cover overestimate in the latter. SEDGES assumes that forest is vegetated land that protrudes above the winter snowpack, whereas the ISLSCP II dataset requires that trees be at least 5 m tall to count toward tree cover and thus excludes short trees and tall shrubs that do not get completely buried by snow.

5.8 Model evaluation: summary and conclusions

The performance of the new SEDGES model has been evaluated for the present day through comparison with numerous reference land datasets. To do this, we have forced SEDGES offline with ERA-Interim data, using our recommended simple single-layer bucket soil hydrology to simulate terrestrial water storage, incorporating a simple scheme to handle a snow-temperature mismatch in the forcing data, and adopting the Planet Simulator's formulation for aerodynamic conductance. Note that the particular forcing data and hydrological implementations impact the analyzed output of the SEDGES simulations. In particular, the spin-up using repeated climatology from 1981–2010 yields conditions that differ from a spin-up using purely historical data. Not only do forcing and hydrology affect the model evaluation, but so do the reference datasets against which it is compared, since even they are not perfect representations of reality. Given all of these particularities, one must view the presented evaluation as a conditional, non-definitive, but yet informative guide to the strengths and the weakness of the model.

The following output variables of SEDGES have been examined: GPP, vegetative carbon, soil organic carbon, LAI, surface albedo, ET, runoff, and forest cover fraction. In comparison to the respective reference datasets, SEDGES simulates each variable at least reasonably well. Relative model strengths lie in the simulation of GPP; relative weaknesses in ET, snow-free albedo, and forest cover fraction. Given the simplicity of its formulation, the strength in GPP is unexpected, particularly the ability of SEDGES to capture the spatial patterns of temporal trends in GPP. Simulated ET is too high in the mid- to high northern latitudes and is too low in the wet tropics. SEDGES captures albedo patterns fairly well outside of boreal summer, but subtler differences among the different desert types and among the nondesert areas are not well-captured during the snow-free season. Simulated forest cover fraction is generally too high.

6 Discussion and conclusions

A new simplified model for the representation of dynamic ecological processes for use in conjunction with climate models has been developed. This new model was combined with a simple hydrological scheme and forced with reanalysis data. In evaluating its performance, we concentrated on the comparison with present-day observation-based data products, while also referring to land surface models (especially ENTS) and their performances for a relative comparison of the strengths and weaknesses of SEDGES.

The quantitative comparison highlighted strengths and weaknesses of the model. A notable strength is that SEDGES's simulation of gross primary productivity is comparable to and sometimes better than that of state-of-the-art dynamic global vegetation models. Our evaluation has also shown that SEDGES performs well in a number of other metrics. Overall, the results show that SEDGES can be used to adequately simulate modern land surface characteristics, including input variables to land–atmosphere coupling schemes used in climate models, as well as key variables of the hydrological and terrestrial carbon cycles.

The most severe weaknesses of our offline-run SEDGES are the aforementioned strong positive ET bias that it yields in the mid-to-high latitudes in winter and early spring and its overprediction of arboreal cover. The latter results from the use of only one plant functional type, which carries over from the SimBA model, on which SEDGES is based. Even an expansion to just two plant functional types (e.g., tree and herbaceous) and simulation of competition between them would entail a significant increase in model complexity. However, SEDGES's use of only one plant functional type, along with a fixed relationship between biomass and both wet-soil LAI and rooting depth hydrology, excessively constrains large-scale vegetation structure by excluding the emergence of location-adapted landscape characteristics. As a result of these issues, simulated vegetation in water-limited regions and wetness-induced changes in biomass should be treated as more heuristic than definitive. Returning to ET, although less severe, the positive ET bias is also seen in the mid-to-high latitudes during the snow-free season. In contrast, the equatorial tropics have a strong negative ET bias due to high cloud cover and concomitant low insolation, which increases canopy resistance and thus reduces transpiration (as seen in Eq. 13). This negative bias occurs for various reasons, including SEDGES's lack of distinction between sunlit and shaded leaves (see below), lack of evaporation from canopy-intercepted water in SEDGES (see Sect. 5.6 for more), and a positive cloud cover and a negative surface insolation bias in the forcing dataset (see Sect. 5.6 for more). The last of these may also be associated with further reductions in transpiration by way of a more stable and humid boundary layer, which would decrease the specific humidity difference between the surface and the lower atmosphere (Δq) and increase the aerodynamic resistance (ra) in Eq. (13). Regardless, ET would be less inhibited by these cloudy and humid conditions if the ratio of intercellular to atmospheric CO2, cica, were given a more realistic treatment that is in line with recent theoretical development and measurements pertaining to optimized stomatal conductance (Medlyn et al.2011; Prentice et al.2014; Lin et al.2015); such treatment is already occurring in some land surface models (e.g., De Kauwe et al.2015; Kala et al.2015; Willeit and Ganopolski2016). The new theory formalizes how and why, at the leaf level, cica increases/decreases with decreased/increased vapor pressure deficit (VPD) (which is approximately proportional to Δq). Because CO2 uptake is proportional to 1-cica (Eq. 6), increased cica results in decreased canopy resistance (in order to maintain the same level of CO2 uptake, i.e., match the light-limited rate of GPP), which results in decreased transpiration (Eq. 13) (unless transpiration is occurring at the maximum rate).

The implementation of variable cica was not made in SEDGES because the simple relationships between cica and VPD that are derived from optimization principles (Medlyn et al.2011, 2012; Prentice et al.2014) are incompatible with the SEDGES framework. The derivations for optimized cica both assume a Fick's law of diffusion transmission between the leaf and the outside air. In order to have consistency in moisture and CO2 fluxes between the land surface and atmosphere, the land surface scheme also needs to use a diffusive exchange between outside air and leaves, or at the least, to reasonably approximate such diffusive exchange. Diffusive exchange of moisture and CO2 is not used in the SEDGES framework. Instead, transfer occurs from surface to atmosphere through a bulk aerodynamic formulation. This formulation only approximates the purely diffusive scheme when canopy resistance (rc) greatly exceeds aerodynamic resistance (ra). As is said in Sect. 2.2.4, in early versions of SEDGES (coupled to Planet Simulator), it was found that this condition (rcra) was frequently violated. That said, in light of the aforementioned ET biases, we intend to incorporate variable cica in future versions of SEDGES. Note that it is not enough to simply adjust the globally fixed cica in the current model to resolve the ET problems because doing so would improve the bias in the mid- to high latitudes while worsening it in the equatorial tropics (or vice versa).

The use of a fixed value of cica is likely to play a more important role in transpiration (by way of affecting the water-use efficiency) than in GPP under non-water-limiting conditions (i.e., when GPP = GPPL, where GPPL is the light-limited rate of GPP (Sect. 2.2.3)). This is because, when rcra (i.e., using the diffusive approximation for heuristic purposes), transpiration is proportional to ΔqGPPL/1-cica. As such, changes in cica (whose values are typically closer to 1 than 0 for C3 plants) have greater relative impact on transpiration than they do in the equations for either RuBP (ribulose 1,5-bisphosphate) regeneration-limited or Rubisco-limited photosynthesis (at least for non-extreme conditions) in the Farquhar model (Farquhar et al.1980). Furthermore, explicit dependence on cica is not critical when modeling GPP. In remote-sensing applications, light use efficiency models without explicit ci dependence have been widely successful in modeling GPP for many years (e.g., Yuan et al.2007, and references within). An inter-model comparison between the aforementioned light use efficiency models and the explicitly ci-dependent enzyme-kinetic-based approaches (Schaefer et al.2012) shows that including explicit ci variation (and hence cica variation) does not apparently improve the modeling of GPP. On the other hand, for arid regions in which water supply is often limiting, the decrease in cica and concomitant increase in water-use efficiency (GPP per amount of transpired water) affects GPP somewhat strongly.

As the SEDGES model was developed, many of the formulations were calibrated, both independently and dependently of model simulations. As noted by Foley et al. (2013), model calibration should ideally be done separately from model evaluation, but this happens only rarely in practice. In this paper, some of the datasets that were most heavily relied on for calibration were also used for model evaluation. These include the MODIS albedo dataset (NASA LP DAAC2008b), the MTE GPP dataset (Jung et al.2011) (multiyear annual mean), the Mueller et al. (2013) ET dataset (multiyear annual mean), and the BNU MODIS-based LAI dataset (Yuan et al.2011). The least well-constrained parameters and relationships that SEDGES is sensitive to are the ϵmax, the soil organic carbon value at which soil albedo saturation occurs (c7), the trmax, the CO2 fertilization function, the low-temperature limitation function, and c4 and c5 in the relationship between biomass and snow-covered albedo. The model is only moderately sensitive to the c12 parameter, which governs the relationship between biomass and the soil bucket depth. The “representative” value for cica is also somewhat unconstrained (e.g., see Prentice et al.2014).

Our offline simulations involving SEDGES and the subsequent comparisons with standard datasets should of course not be considered a complete model evaluation. SEDGES undoubtedly has limitations that cannot be seen from the simulations used in this study. However, from the results of this study and from the formulations that are used in SEDGES, we can anticipate at least some of the situations in which the current model might perform less realistically and/or suggest how to modify SEDGES for these cases.

In Sect. 2.2.3, for the calculation of the light-limited rate of GPP (Eq. 2), ϵmax is the same regardless of the diffuse or direct nature of the absorbed solar radiation. However, this lack of distinction is questionable, as shown by two observation-based studies on the relationship between the diffuse fraction of SW (shortwave) radiation at the top of the canopy and LUE that control for the negative correlation between VPD and diffuse SW fraction in their results (Alton et al.2007; Williams et al.2014). When going from conditions of predominantly direct solar radiation to predominantly diffuse solar radiation, Alton et al. (2007) finds an observed 6 to 33 % increase in LUE in three forests, whereas Williams et al. (2014) finds an  17 % increase in LUE in shrub tundra. The increase in LUE is apparently due to a more even distribution of PAR among the leaves, which reduces light saturation among the sunlit leaves. The distinction between sunlit and shade leaves is missing in our model's single big-leaf approach to canopy radiation, which tacitly assumes a spatially averaged light profile at each level of the canopy (de Pury and Farquhar1997; Monson and Baldocchi2014, p. 355). Not including a sunlit/shade leaf distinction implies that, in the absence of water limitation, our model underpredicts GPP at low sun angles and under cloudy conditions (and overpredicts it for opposite conditions). This pattern is, in fact, reflected in the regional negative GPP biases that are mentioned in Sect. 5.1, with negative productivity biases occurring in regions with reduced downwelling surface SW radiation for their latitude (in the reanalysis), especially in the equatorial tropics. Recent studies show that including a sunlit/shade leaf distinction reduces errors in the simulation of GPP (Yuan et al.2014; Wu et al.2015) (but cf. Schaefer et al.2012). Thus, in a future version of SEDGES, we hope to incorporate the sunlit/shade leaf distinction.

Although gross photosynthesis is somewhat robust to changes in elevation (i.e., atmospheric pressure), there is an appreciable decline at high altitudes and high temperature (Terashima et al.1995) (according to the widely used Farquhar model; Farquhar et al.1980), which should affect places like the Tibetan Plateau and the South American Altiplano. In these regions, the calculations by Terashima et al. (1995) suggest that SEDGES would overestimate the rate of its light-limited gross primary productivity (Eq. 2) by  20 % and by potentially more for areas in these regions that have low biomass (since these would have a stronger feedback between leaf cover fraction and productivity).

Also in Sect. 2.2.2, we stated that SEDGES approximates net primary productivity as half of gross primary productivity (NPP  GPP = 0.5), while noting that this approximation is only valid on the timescales of weeks or longer. On shorter timescales, the ratio may deviate substantially from 0.5, due to the strong short-term dependency of autotrophic respiration on temperature (e.g., Monson and Baldocchi2014), and, as such, carbon uptake by the land surface may not be well-simulated on these short timescales.

With respect to simulating terrestrial carbon pool changes, SEDGES's use of a single soil organic carbon reservoir and its fixed residence time for vegetative carbon may also need to be considered when drawing inferences from model results in transient simulations (e.g., see Friend et al.2014).

These last issues with regards to the carbon cycle are part of a group of model deficiencies in simulating some ecological dynamics. Of these simulation deficiencies, the most severe are probably of the phenological changes (GPP, NPP, transpiration, and especially LAI) associated with green-up in the Northern Hemisphere. More generally, however, it is likely that simplifying assumptions made by SEDGES (namely, constant NPP  GPP as well as its universal temperature limitation function; Appendix Appendix C) lead the model to overestimate the capacity of vegetation to adapt to rapidly changing conditions, especially on daily timescales. As such, SEDGES may underestimate the less positive impacts of extreme weather events on vegetation.

In paleoclimate simulations, periods before the Paleogene may require the modification of SEDGES due to differences in atmosphere oxygen concentration and to the effect of evolution on large-scale plant characteristics and behavior. SEDGES tacitly assumes an atmospheric oxygen concentration near its present-day value of  21 % in both the ϵmax and in the light compensation point in the absence of dark respiration CO2comp. However, CO2comp is proportional to the O2 concentration in the atmosphere (e.g., p. 102 of Monson and Baldocchi2014), which, in turn, affects the CO2 fertilization effect in SEDGES (Eq. 3). Higher/lower O2 would lead to greater/lesser fertilization. Next, for a given CO2 concentration, higher O2 increases photorespiration and decreases gross photosynthesis, both experimentally and according to the aforementioned Farquhar model (Beerling et al.1998) (although the effect was also found to be less pronounced for evolutionarily older taxa). These results suggest that ϵmax should be lowered/raised in SEDGES for higher/lower-than-present atmospheric O2 levels. (In addition, although optimization theory by Prentice et al. (2014) and experimental measurements show that the ratio of intercellular to atmosphere CO2cica increases/decreases with higher/lower atmospheric O2, the increase is found to be only  0.06 when going from the current (21 %) to the extreme high (35 %) O2 level over the last 400 million years (Beerling et al.1998). This cica difference is small in comparison with other uncertainties and inaccuracies in the model and can thus probably be neglected.)

Apart from the effects of changing atmospheric O2 concentration, vegetation characteristics important to the large scale have evolved and changed along with the Earth, and this fact has been recognized by previous researchers (e.g., by excluding plant functional types that have not yet evolved in paleoclimate simulations (e.g., Horton et al.2010; Zhou et al.2012)). We recommend making changes in line with a recent review of climate-relevant changes in plant physiology on geological timescales (Boyce and Lee2017). From the early Paleozoic to the Early Devonian, the landscape was dominated by nonvascular vegetation that lacked stomatal control over water loss, substantial roots, and tree-like aboveground structures. In order to simulate a landscape with this kind of vegetation, canopy resistance (rc) should be set to near 0, Wmax should be set to its minimum value (Wmaxmin) (regardless of biomass), and the biomass residence time τveg should be decreased as needed to keep the highest biomass from going above 1.0 kg C m−2. Starting with the landscape dominance by vascular plants in the Carboniferous, modern-like landscapes of forests, deep-rooted vegetation, and stomatal control over water loss were widespread, and the current version SEDGES should be able to adequately simulate the vegetation of these geological times (provided that the aforementioned effects of differing O2 are included). The Middle and Late Devonian saw the development and expansion of vascular plants and would require a more in-depth analysis on how to properly simulate with SEDGES because it represents a transition between the two aforementioned time periods.

It is important to realize that this paper has focused on SEDGES forced offline, i.e., not fully integrated into a climate or Earth system model. While a stand-alone evaluation against present-day observations or in comparison with other models is not sufficient to make firm statements on how SEDGES would perform in a coupled mode, it is beyond the scope of this paper to evaluate SEDGES within a coupled climate or Earth system model. Nevertheless, the performance of SEDGES in a coupled system will need careful evaluation. First results with the coupled PlaSim-SEDGES model are discussed in Paiewonsky (2017). In general, it is to be expected that coupling may exacerbate or reduce biases present in the offline mode of SEDGES that has been evaluated in this paper. For example, the aforementioned ET biases that SEDGES has (Sect. 5.6) are likely to be diminished through a long-known negative feedback (e.g., Sato et al.1989): excessive ET moistens and stabilizes the boundary layer which reduces the specific humidity difference between the surface and the lower atmosphere (lowers the Δq and increases the ra in Eq. 7), thus reducing ET and thus feeding back negatively onto the original positive ET bias. On the other hand, there are two notable positive feedbacks that may occur in coupled mode. First, the positive snow-free surface albedo bias in the high Arctic in SEDGES (e.g., see July in Fig. 4) may feed back positively in these regions and worsen the initial bias by cooling the surface, thus reducing productivity, thus reducing soil organic carbon, and thus increasing the surface albedo. Second, through another albedo mechanism in boreal spring (the forest–tundra snow albedo feedback, e.g., Bonan et al.1992; Foley et al.1994), the high albedo bias in northern Siberia'a forest–tundra woodland zones would be expected to lead to additional cooling in this region, thus further reducing productivity, thus decreasing biomass, and thus increasing the snow-covered surface albedo, resulting in an exacerbation of the initial positive albedo bias in this region.

For its level of complexity, SEDGES and the land surface framework which it presupposes have the advantages of a flexible time step, canopy control over transpiration, transpiration that is fully coupled to photosynthesis, and vegetative productivity that depends on light and not just temperature and moisture. The short time step option permits SEDGES to be incorporated into models that resolve the diurnal cycle. Because of these strengths, we expect SEDGES to have advantages over vegetation models of similar complexity in simulating vegetation, primary productivity, and transpiration in past geological warm periods in regions that receive little sunlight (e.g., the Eocene Arctic; Maxbauer et al.2014, and references within) and in simulating cloud–vegetation feedbacks in all eras. The SEDGES framework uses photosynthetic-driven stomatal control over evapotranspiration, which is critical when studying hydrological changes on land under altered atmospheric CO2 concentrations, including anthropogenically induced warming scenarios (Betts et al.2007).

The SEDGES framework has the advantages of being easier to understand, process-wise, than more complex vegetation and hydrological schemes. Increased vegetation model complexity can obfuscate model behavior, and it need not improve performance. In an inter-model comparison of the performance of 24 vegetation models at 39 eddy covariance flux tower sites, Schaefer et al. (2012) finds insignificant effects of nitrogen cycle inclusion and of having a light use efficiency or enzyme kinetic approach to productivity on a model's capacity to simulate GPP. Moreover, even when present-day biome distributions are well-simulated by complex vegetation models, model behaviors can diverge drastically from each other under non-present conditions due to differing physiological assumptions (e.g., climatic limitations within plant functional types) that are not well-grounded and that may very well change under novel climatic, atmospheric CO2, and nutrient conditions (Fisher et al.2015).

In conclusion, we feel that SEDGES provides a new viable and computationally efficient alternative to currently implemented terrestrial vegetation/ecological models, in particular inside Earth system models of intermediate complexity and for research on the global-scale interactions between the physical climate system and the terrestrial biosphere.

Code and data availability

Both the SEDGES model code and the code we used to drive SEDGES with external data are available from the following digital repository:

CDOs were obtained from and were used to process much of the data. Version 6.2.1 of NCL (NCAR Command Language) was used to generate the plots (

ERA-Interim data (Dee et al.2011) are available from The atmospheric CO2 data are available from for the Mauna Loa record and from for the ice core data. We used the “” file of the MTE GPP dataset (Jung et al.2011), which is available from The CARBONES GPP dataset was obtained from The vegetative carbon dataset (Olson et al.1985) is available from The HWSD soil organic carbon reference dataset (Wieder et al.2011) is available from The LAI data (Yuan et al.2011) are available from The MODIS white-sky albedo data (NASA LP DAAC2008b) and MODIS snow-free albedo data (NASA LP DAAC2008a) are available from The reference ET dataset (Mueller et al.2013) can, upon registration, be download from The UNH/GRDC runoff dataset (Fekete et al.2002) can be obtained from The ISLSCP II tree cover dataset (DeFries and Hansen2009) is available from The rooting zone plant-available soil water storage capacity dataset (Kleidon and Heimann1998; Hall et al.2006; Kleidon2011) is available from the following URL: Accessing these last two datasets requires registering and signing in.

Appendix A: Notes on the temperature limitation function

In this section, we discuss the rationale behind our particular temperature limitation multiplier, f2(Tsfc), on the light-limited rate of gross primary productivity (Eq. 2) in Sect. 2.2.3.

f2(Tsfc) increases linearly from 0 at Tsfc of 0 C to 1 at Tsfc of 20 C and then plateaus at 1. The linear increase is chosen for its simplicity and because it is analogous to the often-used “growing degree days” (GDDs) metric (with 0 C base) in agriculture and ecology (Kauppi and Posch1985; Prentice et al.1992; Kaplan2001) (save for the shorter timescale). SEDGES's critical temperature of 20 C at which productivity is no longer limited lies close to the middle of the range of optimum temperatures for C3 plants (5–39 C) that was found in a survey of 212 temperature response curves of photosynthesis in the literature for different species and growth temperatures (Way and Yamori2014; Yamori et al.2014). While the critical temperature for the Kleidon (2006b) SimBA model is 5 C, it otherwise uses the same ramp function formulation as SEDGES.

Figure A1Multiyear annual mean of GPP for the original SimBA model (Kleidon2006b) and the two reference datasets, MTE and CARBONES, for 1990–2009 over non-glaciated land.


Figure A2Zonal multiyear monthly means and zonal multiyear annual means of ET for the original SimBA model (Kleidon2006b) and the Mueller et al. (2013) reference dataset for the 1989–2005 period over non-glaciated land.


Temperature limitation on photosynthesis has been well studied (e.g., see reviews by Sage and Kubien2007; Yamori et al.2014), but many uncertainties remain. In land plants, photosynthesis is significantly limited by both high and low temperatures. Individual plants of different species and from different environments have differing optimal temperatures for photosynthesis, such that away from the optimal temperature, photosynthetic rate decreases (Yamori et al.2014). Around the temperature optimum, there exists a range of temperatures for which photosynthesis is nearly as high as at the optimum (Berry and Bjorkman1980; Yamori et al.2014). Both the temperature dependency and the maximum rate of photosynthesis (realized at the temperature optimum) can shift as a plant acclimates to a different set of environmental conditions (Berry and Bjorkman1980; Way and Yamori2014; Yamori et al.2014).

In most plants, temperature limitation on photosynthesis depends on the species, current light levels, the environmental conditions under which the plant grew, and the internal levels of CO2 in the leaves (Berry and Bjorkman1980; Sage and Kubien2007; Way and Yamori2014; Yamori et al.2014). In addition, the current temperature limitation depends on the extent of acclimation to recent environmental conditions (e.g., Dietze2014)

The framework in which we use f2(Tsfc) assumes perfect acclimation and adaptation of the vegetation to the current weather and light conditions. As such, the SEDGES temperature limitation assumes the presence of the most productive set of plants that could ideally grow under those conditions and also allows for their instantaneous adaptation and acclimation to those conditions. In other words, for a given surface temperature (Tsfc), f2(Tsfc) represents the universal physiological constraints on productivity exerted by that (isolated) given temperature.

For high temperatures, a conscious decision was made to not include a decline in productivity in SEDGES. Although it is well-established that high temperatures (i.e., temperatures above the optimum) limit productivity for individual plants, some of that limitation that has been found empirically may have been due to an indirect effect of temperature on increasing the leaf-to-air vapor pressure deficit (for constant ambient relative humidity), which causes the optimum temperature to be lower than when the vapor pressure deficit is held fixed (Lin et al.2012). An even more important point is that a high temperature limitation reflects more the inability of plants that are adapted to a given environment to perform equally well under all conditions that they might face in that environment and less an intrinsic barrier for plant life at that location to adapt to a given high temperature regime (via some combination of species migrations, natural selection, evolution, and/or acclimation). For example, although the Yamori et al. (2014) study found C3 plants to generally have temperature optima between 10 and 35 C, light-saturated and CO2-saturated photosynthesis for the C3 hot desert plant, Rhazya stricta, shows no significant decline with leaf temperature even up to 44 C (Lawson et al.2014). Moreover, its (ordinary) daily photosynthetic rates are nearly as high as those of common C3 agricultural crops, in spite of living under natural conditions of a much higher vapor pressure deficit and presumably much drier soil (Lawson et al.2014), which reduce productivity. Similar results as for Rhazya stricta have been found for the C3 desert shrub, Larrea divaricata (Mooney et al.1978). Thus, the large-scale ability of vascular plant life to adapt to high temperatures through the increased prevalence of heat-tolerant species and phenotypes should be properly accounted for when modeling vegetation and the carbon cycle in past warm climates (such as the Cretaceous) as well as in future climates.

Appendix B: Numerical parameterization of canopy resistance and maximum transpiration rate

Canopy resistance is the greater of an unconstrained canopy resistance (rcu) which is determined by the light-limited rate of canopy photosynthesis (GPPL) and a canopy resistance set by the maximum rate at which the rooting zone can supply water for transpiration. The equation for the unconstrained canopy resistance is essentially derived by incorporating the canopy resistance (rc) value from the previous time step, using the mismatch between GPPW and GPPL in the last time step, and accounting for updates to the variables that affect GPPW and GPPL. Although this formulation of rcu depends on values from the previous time step, this dependency is arbitrary. An expression for rcu that depends only on values in the current time step could be derived and used. (Such an rcu would depend on all the variables that comprise GPPW and GPPL, and these variables would be evaluated at the current time step.)

Our starting equation for deriving the unconstrained canopy resistance uses the formulations of GPPL and GPPW given in Eqs. (2), (4), (5), and (6) and is as follows:


Here, the equivalence statement reflects the first stomatal goal of parsimoniously meeting the light-driven demand for CO2. The future values of GPPW and GPPL are approximated (for simplicity) by neglecting changes in slowly varying variables (e.g., ρ). Equation (B1) represents the standard case for deriving rc. Special cases are discussed below (e.g., when fleaft-1, f2(Tsfct-1) or SWt-1 are zero.) Solving Eq. (B1) for rcu yields


Equation (B1) has provided us with rcu, the canopy resistance that would occur in the complete absence of physiological limitations on plant water loss. In the real world, plants must restrict how much their stomata open when water cannot be extracted from the soil (or from internal plant storage) and transported to the leaf stomata fast enough to balance the loss to the atmosphere through the stomata. Doing so helps to keep water potentials within the plant from falling to levels that would cause “runaway” cavitation (Tyree and Sperry1988) or “hydraulic failure” (McDowell et al.2008). Closing stomata increases stomatal resistances across the leaves, which, on the canopy scale, increases the canopy resistance. In this way, canopy resistance can be constrained by the supply rate of water from the soil to the leaves.

To simulate the maximum supply rate of water for transpiration, SEDGES adapts the simple model of Federer (1982). The original Federer maximum supply rate is directly proportional to the soil wetness fraction multiplied by a fixed constant parameter. This parameter is the absolute maximum transpiration value, which we denote here by trmax. In SEDGES, we extend the original Federer formulation from the limited case of a fully vegetated surface to the case of a mixed surface comprised of green leaves and exposed soil. This is achieved by simply multiplying the original expression by the leaf cover fraction, fleaf. Thus, the supply rate for transpiration, Str, is as follows in SEDGES:

(B3) S tr = f leaf β tr trmax ,

where βtr is a water stress factor that affects the roots' ability to supply water for transpiration and equals the soil wetness fraction, Wfrac, and trmax has a value of 2.78×10-7 m s−1, which is taken from the BETHY (Biosphere Energy-Transfer HYdrology; Knorr2000) vegetation model (Knorr2000) and is unfortunately not well-constrained (Knorr2000; Knauer et al.2015). This maximum transpiration value is achieved only for a fully leaved (fleaf=1), fully wet soil.

The above extension of the original Federer formulation for water supply rate requires some justification. A reader who is familiar with the Lund–Potsdam–Jena Dynamic Global Vegetation Model (LPJ; Sitch et al.2003) model might notice that SEDGES uses a similar formulation as LPJ (Gerten et al.2004); that is, both models multiply the original simple Federer expression (Str=βtr trmax) by a fractional term on the right-hand side: fveg for LPJ and fleaf for SEDGES. However, because the definitions of fveg in LPJ and fleaf in SEDGES differ, it is inappropriate to simply borrow the LPJ formulation and apply it in SEDGES. In LPJ, fveg is the vegetative cover fraction, which is (approximately) the portion of a grid cell in which vegetation resides. Although this portion appears to change to some extent with daily phenology (Gerten et al.2004), the basic idea in LPJ seems to be that the vegetated portion of the grid cell is (essentially) completely vegetated, i.e., having roots, stems, and leaves. As such, it makes sense in LPJ to multiply the original Federer supply rate formulation by fveg because this is the fraction of the grid cell that functions as in the Federer model. In contrast, SEDGES assumes that vegetation resides in the entire grid cell. In other words, SEDGES assumes that vegetation is spread throughout the grid cell uniformly, which means that plant roots and aboveground parts are also distributed uniformly on the large scale. Thus, the vegetative fraction in SEDGES, fveg, is assumed to be 1. In contrast, in SEDGES, fleaf is created by the covering of bare soil by green leaves when looking down from above the canopy. As such, it is not obvious why the simple supply rate formulation in Federer (1982) is multiplied by fleaf in SEDGES. The answer lies in a close examination of the trmax parameter of the original model.

The aforementioned simple supply rate model of Federer (1982) is tested, evaluated, and calibrated in the same paper against a more sophisticated “Type I” (Guswa et al.2002) water uptake model that is forced using site-specific atmospheric observational data. Doing so reveals that the maximum transpiration rate (i.e., the trmax parameter) in the simple model depends on the following input parameters for the Type I model (in decreasing order of strength): rooting density, root internal resistance, depth of the rooting zone, and vegetation height/surface roughness. Although the first three of these dependencies are found to be substantial, those dependencies are not included in the simple model, which takes the trmax parameter to be a constant. In addition, even the more sophisticated Type I water uptake model neglects soil-drying-induced embolisms' increase in root xylem resistance (e.g., Linton et al.1998; Domec et al.2004), which some authors (Domec et al.2009; Javaux et al.2013) feel plays a very significant or even dominant role in the whole plant conductivity of water. A significant increase in rooting zone xylem resistance decreases the maximum transpiration rate within the framework of the Type I model given in Federer (1982).

The missing dependencies in the original simple Federer formulation for maximum supply rate of transpiration presented in the last paragraph are addressed (though somewhat crudely) in SEDGES by multiplying the original formulation by the leaf cover fraction, fleaf. Here, multiplication by fleaf becomes a proxy for the effect of decreasing root embolism with increasing soil wetness fraction in dry soils because fleaf increases with soil wetness fraction for very low values of soil wetness fraction (see Eq. 20). At low biomass values ( 0.25 kg C m−2), fleaf is a proxy for rooting density because fleaf scales with biomass (through Eqs. 5 and 17) and so does rooting density (implicitly). At higher biomass values (> 0.25 kg C m−2), fleaf serves as a proxy for rooting depth since both scale with biomass.

The water supply rate in the preceding discussion, as noted above, directly impacts the minimum canopy resistance, for which we will now derive an expression. Starting from Eqs. (7) and (13), it follows that

(B4) T ET = f leaf ( 1 + r c g a ) C w .

Starting from Eq. (B2), we define rcmin as the value of canopy resistance that matches the supply rate of transpiration in the current time step. That is, we have

(B5) f leaf β tr trmax f leaf 1 + r c min g a ET C w ,

where we note that ETCw is the PET and ga=1ra. Solving Eq. (B3) for rcmin yields

(B6) r c min = PET β tr trmax - 1 r a .

The final canopy resistance is restricted as follows:

(B7) r c = min ( r c max , max ( r c min , r c u , r c min min ) ) ,

where rcminmin is the absolute minimum possible value for rc and is set to 0 for simplicity and rcmax is the maximum possible value for rc and is set to an extremely large value, purely for the sake of model simplicity and elegance.

In deriving above Eqs. (B1), (B3), and (B4), we have excluded all cases in which we would have had to divide by variables with a zero value. These cases are handled as follows: firstly, SEDGES checks to see if the light-limited rate of GPP is zero in the next time step or if the soil is completely dry. If so, then canopy resistance is set to the maximum possible value. In the former case, there is no carbon benefit to keeping stomata open, and in the latter case, the supply rate is zero. If, however, the first check comes back negative for both conditions, then SEDGES checks for the case of parched soil getting rewetted. If this is the case, canopy resistance is set to the minimum value, rcmin. If, however, the second check also comes back negative, then a final check is made to see if the light-limited rate of GPP is zero in the current time step. Under this last scenario, either the sun has risen or the surface temperature has risen above freezing in between the last and current time steps. Here, an expression for rcu can be derived by substituting into Eq. (B1) a time-specifying version of Eq. (2): GPPLt+1=ϵmaxf1(CO2)f2(Tsfc)fleafSW. We then solve for the unconstrained canopy resistance and get the following:


where the asterisk is used to denote our special case of GPPLt+1>0, GPPL=0, Wsoil>0, and fleaft-1>0. After its computation in Eq. (B6), r*cu is restricted in the same way as rcu is in Eq. (B5). In the case of the first, second, and final checks all being negative, we have the standard scenario described by Eqs. (B1), (B3), and (B4).

Appendix C: Output of the original SimBA

As a supplementary experiment, we followed the same procedure as we did in forcing SEDGES (Sect. 4) in running the original SimBA model (Kleidon2006b) (that found in version 15 of Planet Simulator) offline, except that the only equilibrium simulations used preindustrial CO2 and were run for 900 years without carbon cycle acceleration (because it was not needed). The results for GPP and ET are shown in Figs. A1 and A2. Excessively high ET in Fig. A2 results from the Manabe-style parameterization for evaporation from the land surface (Manabe1969), which is at the potential rate as soil wetness is reduced until a given threshold soil wetness fraction (0.25 for SimBA) is reached, after which the surface wetness factor, Cw (Eq. 7), is reduced linearly with the soil wetness fraction until they both reach 0. This scheme for land surface evaporation was used with the first-generation land surface models (e.g., Pitman2003) and is widely known to yield excessive evapotranspiration from the land surface (e.g., see the intermodel comparisons in Chen et al.1997b; Desborough1999) and extremely dry soils (Chen et al.1997b). Such extremely dry soils are also found in the offline-driven SimBA. This, in tandem with a higher soil wetness fraction threshold for commencement of leaf fall (0.25 vs. 0.05 in SEDGES), causes GPP to be unrealistically low throughout the world (Fig. A1).

In line with the findings by Dekker et al. (2010), in this offline version of SimBA, we found multiple steady states when the model was run using preindustrial levels of CO2: equilibrium biomass was slightly higher in some regions when the model was initialized with tropical forest-like biomass of 12 kg C m−2 as opposed to a bare desert initialization. Only the forest-initialized run was subsequently continued through 2010.

Competing interests

The authors declare that they have no conflict of interest.


Pablo Paiewonsky thanks Daniela Dalmonech for alerting us to severe deficiencies in the simulated evapotranspiration in an early version of SEDGES that ultimately led us to revise and much improve the model. Pablo Paiewonsky's work on the model's development was financially supported, in part, by the NSF-funded Alliance for Graduate Education and Professoriate (AGEP) program, by a Carson Carr Graduate Diversity Fellowship from the Office of Diversity, Equity, and Inclusion of the State University of New York, and by the Department of Atmospheric and Environmental Sciences. We would also like to thank the reviewers for their many constructive and useful criticisms that improved the paper.

Edited by: Min-Hui Lo
Reviewed by: Guillaume Marie, Ray Anderson, and three anonymous referees


Alton, P., North, P., and Los, S.: The impact of diffuse sunlight on canopy light-use efficiency, gross photosynthetic product and net ecosystem exchange in three forest biomes, Glob. Change Biol., 13, 776–787, 2007. a, b

Amthor, J. S. and Baldocchi, D. D.: Terrestrial higher plant respiration and net primary production, in: Terrestrial global productivity, chap. 3, edited by: Roy, J., Saugier, B., and Mooney, H. A., Academic Press, 33–59, 2001. a

Anav, A., Friedlingstein, P., Beer, C., Ciais, P., Harper, A., Jones, C., Murray-Tortarolo, G., Papale, D., Parazoo, N. C., Peylin, P., Piao, S., Sitch, S., Viovy, N., Wiltshire, A., and Zhao, M.: Spatiotemporal patterns of terrestrial gross primary production: A review, Rev. Geophys., 53, 785–818, 2015. a, b, c, d, e

Beerling, D., Woodward, F., Lomas, M., Wills, M., Quick, W., and Valdes, P.: The influence of Carboniferous palaeoatmospheres on plant function: an experimental and modelling assessment, Philos. T. R. Soc. B, 353, 131–140, 1998. a, b

Beljaars, A., Brown, A. R., and Wood, N.: A new parametrization of turbulent orographic form drag, Q. J. Roy. Meteor. Soc., 130, 1327–1347, 2004. a

Berry, J. and Bjorkman, O.: Photosynthetic response and adaptation to temperature in higher plants, Ann. Rev. Plant Physiol., 31, 491–543, 1980. a, b, c

Best, M. J., Pryor, M., Clark, D. B., Rooney, G. G., Essery, R. L. H., Ménard, C. B., Edwards, J. M., Hendry, M. A., Porson, A., Gedney, N., Mercado, L. M., Sitch, S., Blyth, E., Boucher, O., Cox, P. M., Grimmond, C. S. B., and Harding, R. J.: The Joint UK Land Environment Simulator (JULES), model description – Part 1: Energy and water fluxes, Geosci. Model Dev., 4, 677–699,, 2011. a, b

Betts, A. K. and Ball, J. H.: Albedo over the boreal forest, J. Geophys. Res.-Atmos., 102, 28901–28909, 1997. a

Betts, R. A., Cox, P. M., Lee, S. E., and Woodward, F. I.: Contrasting physiological and structural vegetation feedbacks in climate change simulations, Nature, 387, 796–799, 1997. a

Betts, R. A., Boucher, O., Collins, M., Cox, P. M., Falloon, P. D., Gedney, N., Hemming, D. L., Huntingford, C., Jones, C. D., Sexton, D. M., and Webb, M. J.: Projected increase in continental runoff due to plant responses to increasing carbon dioxide, Nature, 448, 1037–1041, 2007. a, b

Blyth, E.: Estimating potential evaporation over a hill, Bound.-Lay. Meteorol., 92, 185–193, 1999. a

Bonan, G. B., Pollard, D., and Thompson, S. L.: Effects of boreal forest vegetation on global climate, Nature, 359, 716–718, 1992. a

Bond, W. J., Woodward, F. I., and Midgley, G. F.: The global distribution of ecosystems in a world without fire, New Phytol., 165, 525–538, 2005. a

Bowring, S. P. K., Miller, L. M., Ganzeveld, L., and Kleidon, A.: Applying the concept of “energy return on investment” to desert greening of the Sahara/Sahel using a global climate model, Earth Syst. Dynam., 5, 43–53,, 2014. a

Boyce, C. K. and Lee, J.-E.: Plant Evolution and Climate Over Geological Timescales, Ann. Rev. Earth Pl. Sc., 45, 61–87,, 2017. a

Brovkin, V., Ganopolski, A., and Svirezhev, Y.: A continuous climate-vegetation classification for use in climate-biosphere studies, Ecol. Model., 101, 251–261, 1997. a

Brovkin, V., Bendtsen, J., Claussen, M., Ganopolski, A., Kubatzki, C., Petoukhov, V., and Andreev, A.: Carbon cycle, vegetation, and climate dynamics in the Holocene: Experiments with the CLIMBER-2 model, Global Biogeochem. Cy., 16, 1139,, 2002. a

Buckley, T. N. and Schymanski, S. J.: Stomatal optimisation in relation to atmospheric CO2, New Phytol., 201, 372–377, 2014. a

Campbell, G. and Norman, J.: An introduction to environmental biophysics, 2nd edn., Springer Science + Business Media Inc., New York, NY, USA, 286 pp., 1998. a

Cao, L., Bala, G., Caldeira, K., Nemani, R., and Ban-Weiss, G.: Importance of carbon dioxide physiological forcing to future climate change, P. Natl. Acad. Sci. USA, 107, 9513–9518, 2010. a

Chen, F., Janjić, Z., and Mitchell, K.: Impact of atmospheric surface-layer parameterizations in the new land-surface scheme of the NCEP mesoscale Eta model, Bound.-Lay. Meteorol., 85, 391–421, 1997a. a

Chen, T. H., Henderson-Sellers, A., Milly, P., et al.: Cabauw experimental results from the project for intercomparison of land-surface parameterization schemes, J. Climate, 10, 1194–1215, 1997b. a, b

Chen, Y., Yang, K., Zhou, D., Qin, J., and Guo, X.: Improving the Noah land surface model in arid regions with an appropriate parameterization of the thermal roughness length, J. Hydrometeorol., 11, 995–1006, 2010. a

Clark, D. B., Mercado, L. M., Sitch, S., Jones, C. D., Gedney, N., Best, M. J., Pryor, M., Rooney, G. G., Essery, R. L. H., Blyth, E., Boucher, O., Harding, R. J., Huntingford, C., and Cox, P. M.: The Joint UK Land Environment Simulator (JULES), model description – Part 2: Carbon fluxes and vegetation dynamics, Geosci. Model Dev., 4, 701–722,, 2011. a

Cowan, I. and Farquhar, G.: Stomatal function in relation to leaf metabolism and environment, in: Symposia of the Society for Experimental Biology, edited by: Jennings, D., chap. 3, Cambridge University Press, Cambridge, UK, 471–505, 1977. a

Cox, P., Betts, R., Bunton, C., Essery, R., Rowntree, P., and Smith, J.: The impact of new land surface physics on the GCM simulation of climate and climate sensitivity, Clim. Dynam., 15, 183–203, 1999. a, b, c, d

Cramer, W., Kicklighter, D., Bondeau III, A., B. M., Churkina, G., Nemry, B., Ruimy, A., Schloss, A., and the participants of the Potsdam NPP model intercomparison: Comparing global models of terrestrial net primary productivity (NPP): overview and key results, Glob. Change Biol., 5, 1–15, 1999. a, b, c

Czikowsky, M. J. and Fitzjarrald, D. R.: Detecting rainfall interception in an Amazonian rain forest with eddy flux measurements, J. Hydrology, 377, 92–105, 2009. a, b

Da Rocha, H. R., Goulden, M. L., Miller, S. D., Menton, M. C., Pinto, L. D., de Freitas, H. C., and Silva Figueira, A. M.: Seasonality of water and heat fluxes over a tropical forest in eastern Amazonia, Ecol. Appl., 14, 22–32, 2004. a, b

Dai, A. and Trenberth, K. E.: Estimates of freshwater discharge from continents: Latitudinal and seasonal variations, J. Hydrometeorol., 3, 660–687, 2002. a

De Kauwe, M. G., Kala, J., Lin, Y.-S., Pitman, A. J., Medlyn, B. E., Duursma, R. A., Abramowitz, G., Wang, Y.-P., and Miralles, D. G.: A test of an optimal stomatal conductance scheme within the CABLE land surface model, Geosci. Model Dev., 8, 431–452,, 2015. a

de Pury, D. and Farquhar, G.: Simple scaling of photosynthesis from leaves to canopies without the errors of big-leaf models, Plant Cell Environ., 20, 537–557, 1997. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, I., Biblot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Greer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Holm, E. V., Isaksen, L., Kallberg, P., Kohler, M., Matricardi, M., McNally, A. P., Mong-Sanz, B. M., Morcette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thepaut, J. N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. a, b

DeFries, R. and Hansen, M.: ISLSCP II Continuous Fields of Vegetation Cover, 1992–1993, Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, USA, (last access: February 2016), 2009. a, b, c

Dekker, S. C., de Boer, H. J., Brovkin, V., Fraedrich, K., Wassen, M. J., and Rietkerk, M.: Biogeophysical feedbacks trigger shifts in the modelled vegetation-atmosphere system at multiple scales, Biogeosciences, 7, 1237–1245,, 2010. a, b

DeLucia, E., Drake, J. E., Thomas, R. B., and Gonzalez-Meler, M.: Forest carbon use efficiency: is respiration a constant fraction of gross primary production?, Glob. Change Biol., 13, 1157–1167, 2007. a, b, c

Desborough, C.: Surface energy balance complexity in GCM land surface models, Clim. Dynam., 15, 389–403, 1999. a, b, c

Dewar, R. C.: A simple model of light and water use evaluated for Pinus radiata, Tree Physiol., 17, 259–265, 1997. a, b, c, d, e, f, g

Dietze, M. C.: Gaps in knowledge and data driving uncertainty in models of photosynthesis, Photosynth. Res., 119, 3–14, 2014. a

Dolinar, E. K., Dong, X., and Xi, B.: Evaluation and intercomparison of clouds, precipitation, and radiation budgets in recent reanalyses using satellite-surface observations, Clim. Dynam., 46, 2123–2144, 2016. a, b, c

Domec, J.-C., Warren, J., Meinzer, F., Brooks, J., and Coulombe, R.: Native root xylem embolism and stomatal closure in stands of Douglas-fir and ponderosa pine: mitigation by hydraulic redistribution, Oecologia, 141, 7–16, 2004. a

Domec, J.-C., Noormets, A., King, J. S., Sun, G., McNulty, S. G., Gavazzi, M. J., Boggs, J. L., and Treasure, E. A.: Decoupling the influence of leaf and root hydraulic conductances on stomatal conductance and its sensitivity to vapour pressure deficit as soil dries in a drained loblolly pine plantation, Plant Cell Environ., 32, 980–991, 2009. a

Dutra, E., Balsamo, G., Viterbo, P., Miranda, P. M., Beljaars, A., Schär, C., and Elder, K.: An improved snow scheme for the ECMWF land surface model: description and offline validation, J. Hydrometeorol., 11, 899–916, 2010. a

Dutra, E., Kotlarski, S., Viterbo, P., Balsamo, G., Miranda, P., Schär, C., Bissolli, P., and Jonas, T.: Snow cover sensitivity to horizontal resolution, parameterizations, and atmospheric forcing in a land surface model, J. Geophys. Res.-Atmos., 116, D21109,, 2011. a

Essery, R.: Boreal forests and snow in climate models, Hydrol. Proc., 12, 1561–1567, 1998. a

Etheridge, D., Steele, L., Langenfelds, R., Francey, R., Barnola, J.-M., and Morgan, V.: Historical CO2 records from the Law Dome DE08, DE08-2, and DSS ice cores, Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennesse, USA, 1998. a

Farquhar, G., von Caemmerer, S., and Berry, J.: A Biochemical Model of Photosynthetic CO2 Assimilation in Leaves of C3 Species, Planta, 149, 78–90, 1980. a, b, c

Federer, C. A.: Transpirational supply and demand: plant, soil, and atmospheric effects evaluated by simulation, Water Resour. Res., 18, 355–362, 1982. a, b, c, d

Fekete, B. M., Vörösmarty, C. J., and Grabs, W.: High-resolution fields of global runoff combining observed river discharge and simulated water balances, Global Biogeochem. Cy., 16, 1042,, 2002. a, b, c, d

Fisher, R. A., Muszala, S., Verteinstein, M., Lawrence, P., Xu, C., McDowell, N. G., Knox, R. G., Koven, C., Holm, J., Rogers, B. M., Spessa, A., Lawrence, D., and Bonan, G.: Taking off the training wheels: the properties of a dynamic vegetation model without climate envelopes, CLM4.5(ED), Geosci. Model Dev., 8, 3593–3619,, 2015. a

Foley, A. M., Dalmonech, D., Friend, A. D., Aires, F., Archibald, A. T., Bartlein, P., Bopp, L., Chappellaz, J., Cox, P., Edwards, N. R., Feulner, G., Friedlingstein, P., Harrison, S. P., Hopcroft, P. O., Jones, C. D., Kolassa, J., Levine, J. G., Prentice, I. C., Pyle, J., Vázquez Riveiros, N., Wolff, E. W., and Zaehle, S.: Evaluation of biospheric components in Earth system models using modern and palaeo-observations: the state-of-the-art, Biogeosciences, 10, 8305–8328,, 2013. a

Foley, J. A., Kutzbach, J. E., Coe, M. T., and Levis, S.: Feedbacks between climate and boreal forests during the Holocene epoch, Nature, 371, 52–54, 1994. a

Franks, P. J., Adams, M. A., Amthor, J. S., Barbour, M. M., Berry, J. A., Ellsworth, D. S., Farquhar, G. D., Ghannoum, O., Lloyd, J., McDowell, N., Norby, R. J., Tissue, D. T., and von Caemmerer, S.: Sensitivity of plants to changing atmospheric CO2 concentration: from the geological past to the next century, New Phytol., 197, 1077–1094, 2013. a, b

Friend, A. D., Lucht, W., Rademacher, T. T., Keribin, R., Betts, R., Cadule, P., Ciais, P., Clark, D. B., Dankers, R., Falloon, P. D., Ito, A., Kahana, R., Kleidon, A., Lomas, M. R., Nishina, K., Ostberg, S., Pavlick, R., Peylin, P., Schaphoff, S., Vuichard, N., Warszawski, L., Wiltshire, A., and Ian Woodward, F.: Carbon residence time dominates uncertainty in terrestrial vegetation responses to future climate and atmospheric CO2, P. Natl. Acad. Sci. USA, 111, 3280–3285, 2014. a

Gerten, D., Schaphoff, S., Haberlandt, U., Lucht, W., and Sitch, S.: Terrestrial vegetation and water balance-hydrological evaluation of a dynamic global vegetation model, J. Hydrology, 286, 249–270, 2004. a, b

Gibbs, H.: Olson's major world ecosystem complexes ranked by carbon in live vegetation: An updated database using the GLC2000 land cover product (NDP-017b), Carbon Dioxide Information Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA,, 2006. a, b, c

Gifford, R. M.: Plant respiration in productivity models: conceptualisation, representation and issues for global terrestrial carbon-cycle research, Funct. Plant Biol., 30, 171–186, 2003. a

Giorgi, F. and Avissar, R.: Representation of heterogeneity effects in earth system modeling: Experience from land surface modeling, Rev. Geophys., 35, 413–437, 1997. a

Gower, S. T., Kucharik, C. J., and Norman, J. M.: Direct and indirect estimation of leaf area index, fAPAR, and net primary production of terrestrial ecosystems, Remote Sens. Environ., 70, 29–51, 1999. a, b

Guan, K., Pan, M., Li, H., Wolf, A., Wu, J., Medvigy, D., Caylor, K. K., Sheffield, J., Wood, E. F., Malhi, Y., Liang, M., Kimball, J. S., Saleska, S. R., Berry, J., Joiner, J., and Lyapustin, A. I.: Photosynthetic seasonality of global tropical forests constrained by hydroclimate, Nat. Geosci., 8, 284–289, 2015. a

Guswa, A.: Soil-moisture limits on plant uptake: An upscaled relationship for water-limited ecosystems, Adv. Water Resour., 28, 543–552, 2005. a

Guswa, A. J., Celia, M., and Rodriguez-Iturbe, I.: Models of soil moisture dynamics in ecohydrology: A comparative study, Water Resour. Res., 38, 1166,, 2002. a, b

Hagemann, S.: An improved land surface parameter dataset for global and regional climate models, Tech. Rep. 336, Max-Planck-Institut für Meteorologie, available at:, 2002. a, b, c, d

Hall, F. G., Brown de Colstoun, E., Collatz, G. J., Landis, D., Dirmeyer, P., Betts, A., Huffman, G. J., Bounoua, L., and Meeson, B.: ISLSCP Initiative II global data sets: Surface boundary conditions and atmospheric forcings for land-atmosphere studies, J. Geophys. Res.-Atmos., 111, D22S01,, 2006. a, b, c

He, L., Chen, J. M., Pisek, J., Schaaf, C. B., and Strahler, A. H.: Global clumping index map derived from the MODIS BRDF product, Remote Sens. Environ., 119, 118–130, 2012. a, b

Hemming, D., Betts, R., and Collins, M.: Sensitivity and uncertainty of modelled terrestrial net primary productivity to doubled CO2 and associated climate change for a relatively large perturbed physics ensemble, Agr. Forest Meteorol., 170, 79–88, 2013. a, b

Heskel, M. A., Bitterman, D., Atkin, O. K., Turnbull, M. H., and Griffin, K. L.: Seasonality of foliar respiration in two dominant plant species from the Arctic tundra: response to long-term warming and short-term temperature variability, Funct. Plant Biol., 41, 287–300, 2014. a

Hickler, T., Smith, B., Prentice, I. C., Mjöfors, K., Miller, P., Arneth, A., and Sykes, M. T.: CO2 fertilization in temperate FACE experiments not representative of boreal and tropical forests, Glob. Change Biol., 14, 1531–1542, 2008. a

Hikosaka, K. and Hirose, T.: Leaf angle as a strategy for light competition: optimal and evolutionarily stable light-extinction coefficient within a leaf canopy, Ecoscience, 4, 501–507, 1997. a

Holden, P. B., Edwards, N. R., Gerten, D., and Schaphoff, S.: A model-based constraint on CO2 fertilisation, Biogeosciences, 10, 339–355,, 2013. a, b

Horton, D. E., Poulsen, C. J., and Pollard, D.: Influence of high-latitude vegetation feedbacks on late Palaeozoic glacial cycles, Nat. Geosci., 3, 572–577, 2010. a

Houldcroft, C. J., Grey, W. M., Barnsley, M., Taylor, C. M., Los, S. O., and North, P. R.: New vegetation albedo parameters and global fields of soil background albedo derived from MODIS for use in a climate model, J. Hydrometeorol., 10, 183–198, 2009. a

Huntingford, C., Blyth, E., Wood, N., Hewer, F., and Grant, A.: The effect of orography on evaporation, Bound.-Lay. Meteorol., 86, 487–504, 1998. a

Jarvis, P. G. and McNaughton, K.: Stomatal control of transpiration: scaling up from leaf to region, Adv. Ecol. Res., 15, 1–49, 1986. a

Javaux, M., Couvreur, V., Vanderborght, J., and Vereecken, H.: Root water uptake: from three-dimensional biophysical processes to macroscopic modeling approaches, Vadose Z. J., 12, 4,, 2013. a

Jenkinson, D., Andrew, S., Lynch, J., Goss, M., and Tinker, P.: The turnover of organic carbon and nitrogen in soil [and discussion], Philos. T. R. Soc. B, 329, 361–368, 1990. a, b

Jiang, L., Yan, Y., Hararuk, O., Mikle, N., Xia, J., Shi, Z., Tjiputra, J., Wu, T., and Luo, Y.: Scale-Dependent Performance of CMIP5 Earth System Models in Simulating Terrestrial Vegetation Carbon, J. Climate, 28, 5217–5232, 2015. a, b

Jung, M., Reichstein, M., Margolis, H. A., Cescatti, A., Richardson, A. D., Arain, M. A., Arneth, A., Bernhofer, C., Bonal, D., Chen, J., Gianelle, D., Gobron, N., Kiely, G., Kutsch, W., Lasslop, G.,Law, B. E., Lindroth, A., Merbold, L., Montagnani, L., Moors, E. J., Papale, D., Sottocornola, M., Vaccari, F., and Williams, C.: Global patterns of land-atmosphere fluxes of carbon dioxide, latent heat, and sensible heat derived from eddy covariance, satellite, and meteorological observations, J. Geophys. Res.-Biogeo., 116, G00J07,, 2011. a, b, c, d, e, f

Kala, J., De Kauwe, M. G., Pitman, A. J., Lorenz, R., Medlyn, B. E., Wang, Y.-P., Lin, Y.-S., and Abramowitz, G.: Implementation of an optimal stomatal conductance scheme in the Australian Community Climate Earth Systems Simulator (ACCESS1.3b), Geosci. Model Dev., 8, 3877–3889,, 2015. a

Kaplan, J. O.: Geophysical applications of vegetation modeling, PhD thesis, Lund University, 2001. a

Kauppi, P. and Posch, M.: Sensitivity of boreal forests to possible climatic warming, Clim. Change, 7, 45–54, 1985. a

Keeling, C. D., Bacastow, R. B., Bainbridge, A. E., Ekdahl, C. A., Guenther, P. R., Waterman, L. S., and Chin, J. F.: Atmospheric carbon dioxide variations at Mauna Loa observatory, Hawaii, Tellus, 28, 538–551, 1976. a

Kelliher, F., Leuning, R., Raupach, M., and Schulze, E.-D.: Maximum conductances for evaporation from global vegetation types, Agr. Forest Meteorol., 73, 1–16, 1995. a

Kleidon, A.: Quantifying the biologically possible range of steady-state soil and surface climates with climate model simulations, Biologia, 61, 234–239, 2006a. a

Kleidon, A.: The climate sensitivity to human appropriation of vegetation productivity and its thermodynamic characterization, Global Planet. Change, 54, 109–127, 2006b. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Kleidon, A.: ISLSCP II Total plant-available soil water storage capacity of the rooting zone, Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, USA,, 2011. a, b, c

Kleidon, A. and Heimann, M.: A method of determining rooting depth from a terrestrial biosphere model and its impacts on the global water and carbon cycle, Glob. Change Biol., 4, 275–286, 1998. a, b, c

Klimarechenzentrum, D.: The ECHAM3 Atmospheric General Circulation Model, Tech. Rep. 6, Max-Planck-Institut für Meteorologie, 1993. a

Knauer, J., Werner, C., and Zaehle, S.: Evaluating stomatal models and their atmospheric drought response in a land surface scheme: A multibiome analysis, J. Geophys. Res.-Biogeo., 120, 1894–1911, 2015. a

Knorr, W.: Annual and interannual CO2 exchanges of the terrestrial biosphere: Process-based simulations and uncertainties, Global Ecol. Biogeogr., 9, 225–252, 2000. a, b, c, d

Knorr, W. and Schnitzler, K.-G.: Enhanced albedo feedback in North Africa from possible combined vegetation and soil-formation processes, Climate Dynam., 26, 55–63, 2006. a, b

Köchy, M., Hiederer, R., and Freibauer, A.: Global distribution of soil organic carbon – Part 1: Masses and frequency distributions of SOC stocks for the tropics, permafrost regions, wetlands, and the world, SOIL, 1, 351–365,, 2015. a

Krinner, G., Viovy, N., Noblet-Ducoudré, N. de, Oge?e, J., Polcher, J., Friedlingstein, P., Ciais, P., Sitch, S., and Prentice, I. C.: A dynamic global vegetation model for studies of the coupled atmosphere-biosphere system, Global Biogeochem. Cy., 19, GB1015,, 2005. 

Lawlor, D. and Cornic, G.: Photosynthetic carbon assimilation and associated metabolism in relation to water deficits in higher plants, Plant Cell Environ., 25, 275–294, 2002. a

Lawson, T., Davey, P. A., Yates, S. A., Bechtold, U., Baeshen, M., Baeshen, N., Mutwakil, M. Z., Sabir, J., Baker, N. R., and Mullineaux, P. M.: C3 photosynthesis in the desert plant Rhazya stricta is fully functional at high temperatures and light intensities, New Phytol., 201, 862–873, 2014. a, b

Lee, E., Felzer, B. S., and Kothavala, Z.: Effects of nitrogen limitation on hydrological processes in CLM4-CN, J. Adv. Model. Earth Syst., 5, 741–754,, 2013. a

Levis, S., Foley, J. A., and Pollard, D.: Large-scale vegetation feedbacks on a doubled CO2 climate, J. Climate, 13, 1313–1325, 2000. a

Li, R. and Arora, V. K.: Effect of mosaic representation of vegetation in land surface schemes on simulated energy and carbon balances, Biogeosciences, 9, 593–605,, 2012. a, b

Lin, Y.-S., Medlyn, B. E., and Ellsworth, D. S.: Temperature responses of leaf net photosynthesis: the role of component processes, Tree Physiol., 32, 219–231, 2012. a

Lin, Y.-S., Medlyn, B. E., Duursma, R. A., et al.: Optimal stomatal behaviour around the world, Nat. Clim. Change, 5, 459–464, 2015. a

Linton, M., Sperry, J. S., and Williams, D.: Limits to water transport in Juniperus osteosperma and Pinus edulis: implications for drought tolerance and regulation of transpiration, Funct. Ecol., 12, 906–911, 1998. a

Litton, C. M., Raich, J. W., and Ryan, M. G.: Carbon allocation in forest ecosystems, Glob. Change Biol., 13, 2089–2109, 2007. a, b, c

Lloyd, J. and Taylor, J.: On the temperature dependence of soil respiration, Funct. Ecol., 8, 315–323, 1994. a

Long, S.: Modification of the response of photosynthetic productivity to rising temperature by atmospheric CO2 concentrations: has its importance been underestimated?, Plant Cell Environ., 14, 729–739, 1991. a

Loranty, M. M., Berner, L. T., Goetz, S. J., Jin, Y., and Randerson, J. T.: Vegetation controls on northern high latitude snow-albedo feedback: observations and CMIP5 model simulations, Glob. Change Biol., 20, 594–606, 2014. a

Lorenz, C. and Kunstmann, H.: The hydrological cycle in three state-of-the-art reanalyses: Intercomparison and performance analysis, J. Hydrometeorol., 13, 1397–1420, 2012. a, b

Louis, J., Tiedtke, M., and Geleyn, J.: A short history of the operational PBL-parameterization at ECMWF, in: Proc. Workshop on Planetary Boundary Layer Parameterization, ECMWF, Reading, UK, 59–79, 1982. a

Louis, J.-F.: A parametric model of vertical eddy fluxes in the atmosphere, Bound.-Lay. Meteorol., 17, 187–202, 1979. a

Lunkeit, F., Bottinger, M., Fraedrich, K., Jansen, H., Kirk, E., Kleidon, A., and Luksch, U.: Planet Simulator-reference manual, version 15.0, Tech. rep., Meteorologisches Institut, Universität Hamburg, 2007. a, b, c, d

Lunkeit, F., Borth, H., Bottinger, M., Fraedrich, K., Jansen, H., Kirk, E., Kleidon, A., Luksch, U., Paiewonsky, P., Schubert, S., Sielmann, S., and Wan, H.: Planet Simulator-reference manual, version 16, Tech. rep., Meteorologisches Institut, Universität Hamburg, 2011. a, b, c, d, e, f, g, h

Manabe, S.: Climate and the ocean circulation 1: i. The atmospheric circulation and the hydrology of the earth's surface, Mon. Weather Rev., 97, 739–774, 1969. a

Mason, P.: The formation of areally-averaged roughness lengths, Q. J. Roy. Meteor. Soc., 114, 399–420, 1988. a

Maxbauer, D. P., Royer, D. L., and LePage, B. A.: High Arctic forests during the middle Eocene supported by moderate levels of atmospheric CO2, Geology, 42, 1027–1030, 2014. a

McDowell, N., Pockman, W. T., Allen, C. D., Breshears, D. D., Cobb, N., Kolb, T., Plaut, J., Sperry, J., West, A., Williams, D. G., and Yepez, E. A.: Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought?, New Phytol., 178, 719–739, 2008. a

McGuire, A. D., Melillo, J., Joyce, L., Kicklighter, D., Grace, A., Moore, B. I., and Vorosmarty, C.: Interactions between carbon and nitrogen dynamics in estimating net primary productivity for potential vegetation in North America, Global Biogeochem. Cy., 6, 101–124, 1992. a, b

McLaughlin, B. C., Xu, C.-Y., Rastetter, E. B., and Griffin, K. L.: Predicting ecosystem carbon balance in a warming Arctic: the importance of long-term thermal acclimation potential and inhibitory effects of light on respiration, Glob. Change Biol., 20, 1901–1912, 2014. a

Medlyn, B. E., Duursma, R. A., Eamus, D., Ellsworth, D. S., Prentice, I. C., Barton, C. V., Crous, K. Y., de Angelis, P., Freeman, M., and Wingate, L.: Reconciling the optimal and empirical approaches to modelling stomatal conductance, Glob. Change Biol., 17, 2134–2144, 2011. a, b, c, d

Medlyn, B. E., Duursma, R. A., Eamus, D., Ellsworth, D. S., Colin Prentice, I., Barton, C. V., Crous, K. Y., Angelis, P., Freeman, M., and Wingate, L.: Reconciling the optimal and empirical approaches to modelling stomatal conductance, Glob. Change Biol., 18, 3476–3476, 2012. a, b

Medlyn, B. E., De Kauwe, M. G., Lin, Y.-S., Knauer, J., Duursma, R. A., Williams, C. A., Arneth, A., Clement, R., Isaac, P., Limousin, J.-M., Linderson, M.-L., Meir, P., St. Paul, N.-M., and Wingate, L.: How do leaf and ecosystem measures of water-use efficiency compare?, New Phytol., 216, 758–770,, 2017. a, b

Melton, J. R. and Arora, V. K.: Sub-grid scale representation of vegetation in global land surface schemes: implications for estimation of the terrestrial carbon sink, Biogeosciences, 11, 1021–1036,, 2014. a

Miralles, D. G., De Jeu, R. A. M., Gash, J. H., Holmes, T. R. H., and Dolman, A. J.: Magnitude and variability of land evaporation and its components at the global scale, Hydrol. Earth Syst. Sci., 15, 967–981,, 2011. a

Monson, R. and Baldocchi, D.: Terrestrial biosphere-atmosphere fluxes, Cambridge University Press, 507 pp., 2014. a, b, c, d

Monteith, J., Huda, A., and Midya, D.: RESCAP: a resource capture model for sorghum and pearl millet, in: Modelling the growth and development of sorghum and pearl millet, edited by: Virmani, S., Tandon, H., and Alagarswamy, G., Research Bulletin, ICRISAT, Patancheru, India, 12, 30–34, 1989. a, b

Moody, E. G., King, M. D., Schaaf, C. B., Hall, D. K., and Platnick, S.: Northern Hemisphere five-year average (2000–2004) spectral albedos of surfaces in the presence of snow: Statistics computed from Terra MODIS land products, Remote Sens. Environ., 111, 337–345, 2007. a, b

Mooney, H. A., Björkman, O., and Collatz, G. J.: Photosynthetic acclimation to temperature in the desert shrub, Larrea divaricata: I. Carbon dioxide exchange characteristics of intact leaves, Plant Physiol., 61, 406–410, 1978. a

Mueller, B., Hirschi, M., Jimenez, C., Ciais, P., Dirmeyer, P. A., Dolman, A. J., Fisher, J. B., Jung, M., Ludwig, F., Maignan, F., Miralles, D. G., McCabe, M. F., Reichstein, M., Sheffield, J., Wang, K., Wood, E. F., Zhang, Y., and Seneviratne, S. I.: Benchmark products for land evapotranspiration: LandFlux-EVAL multi-data set synthesis, Hydrol. Earth Syst. Sci., 17, 3707–3720,, 2013 a, b, c, d, e, f, g, h, i, j

NASA LP DAAC: MODIS Terra + Aqua BRDF/Albedo Snow-free Quality 16-Day L3 Global 0.05Deg CMG V005 (MCD43C2), NASA EOSDIS Land Processes DAAC, USGS Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota, available at: (last access: February 2016), 2008a. a, b

NASA LP DAAC: MODIS Terra + Aqua BRDF/Albedo 16-Day L3 Global 0.05Deg CMG V005 (MCD43C3), NASA EOSDIS Land Processes DAAC, USGS Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota, available at: (last access: December 2014), 2008b. a, b, c, d

Nilson, T.: A theoretical analysis of the frequency of gaps in plant stands, Agr. Meteorol., 8, 25–38, 1971. a

Oke, T.: Boundary layer climates, Methuen, London, 2nd edn., 464 pp., 1987. a, b

Olson, D. M., Dinerstein, E., Wikramanayake, E. D., Burgess, N. D., Powell, G. V., Underwood, E. C., D'amico, J. A., Itoua, I., Strand, H. E., Morrison, J. C., Loucks, C. J., Allnutt, T. F., Ricketts, T. H., Kura, Y., Lamoreux, J. F., Wettengel, W. W., Hedao, P., and Kassem, K. R.: Terrestrial Ecoregions of the World: A New Map of Life on Earth A new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity, BioScience, 51, 933–938, 2001. a

Olson, J., Watts, J., and Allison, L.: Major world ecosystem complexes ranked by carbon in live vegetation: A Database (NDP-017), Carbon Dioxide Information Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA, (last access: March 2016), 1985. a, b, c, d, e, f

Paiewonsky, P.: State dependency of the forest-tundra-short wave feedback: comparing the mid-Pliocene and pre-industrial eras using a newly-developed vegetation model, PhD thesis, State University of New York at Albany, 2017. a

Pavlick, R., Drewry, D. T., Bohn, K., Reu, B., and Kleidon, A.: The Jena Diversity-Dynamic Global Vegetation Model (JeDi-DGVM): a diverse approach to representing terrestrial biogeography and biogeochemistry based on plant functional trade-offs, Biogeosciences, 10, 4137–4177,, 2013. a, b

Piao, S., Sitch, S., Ciais, P., Friedlingstein, P., Peylin, P., Wang, X., Ahlström, A., Anav, A., Canadell, J. G., Cong, N., Huntingford, C., Jung, M., Levis, S., Levy, P. E., Li, J., Lin, X., Lomas, M. R., Lu, M., Luo, Y., Ma, Y., Myneni, R. B., Poulter, B., Sun, Z., Wang, T., Viovy, N., Zaehle, S., and Zeng, N.: Evaluation of terrestrial carbon cycle models for their response to climate variability and to CO2 trends, Glob. Change Biol., 19, 2117–2132, 2013. a, b

Pisek, J., Chen, J. M., Lacaze, R., Sonnentag, O., and Alikas, K.: Expanding global mapping of the foliage clumping index with multi-angular POLDER three measurements: Evaluation and topographic compensation, ISPRS J. Photogramm., 65, 341–346, 2010. a, b

Pisek, J., Sonnentag, O., Richardson, A. D., and Mõttus, M.: Is the spherical leaf inclination angle distribution a valid assumption for temperate and boreal broadleaf tree species?, Agr. Forest Meteorol., 169, 186–194, 2013. a

Pitman, A.: The evolution of, and revolution in, land surface schemes designed for climate models, International J. Climatol., 23, 479–510, 2003. a, b

Prentice, I. C., Cramer, W., Harrison, S. P., Leemans, R., Monserud, R. A., and Solomon, A. M.: Special Paper: A Global Biome Model Based on Plant Physiology and Dominance, Soil Properties and Climate, J. Biogeogr., 19, 117–134, 1992. a

Prentice, I. C., Dong, N., Gleason, S. M., Maire, V., and Wright, I. J.: Balancing the costs of carbon gain and water transport: testing a new theoretical framework for plant functional ecology, Ecol. Lett., 17, 82–91, 2014. a, b, c, d, e, f, g

Prentice, I. C., Liang, X., Medlyn, B. E., and Wang, Y.-P.: Reliable, robust and realistic: the three R's of next-generation land-surface modelling, Atmos. Chem. Phys., 15, 5987–6005,, 2015. a, b

Ramankutty, N. and Foley, J. A.: Estimating historical changes in global land cover: Croplands from 1700 to 1992, Global Biogeochem. Cy., 13, 997–1027, 1999. a, b, c, d

Ramankutty, N., Evan, A. T., Monfreda, C., and Foley, J. A.: Farming the planet: 1. Geographic distribution of global agricultural lands in the year 2000, Global Biogeochem. Cy., 22, GB1003,, 2008. a

Raupach, M.: Influences of local feedbacks on land–air exchanges of energy and carbon, Glob. Change Biol., 4, 477–494, 1998. a, b, c, d

Raupach, M. and Finnigan, J.: Single layer models of evaporation from plant canopies are incorrect but useful, whereas multilayer models are correct but useless: discuss, Australian J. Plant Physiol., 15, 705–716, 1988. a

Rechid, D., Raddatz, T. J., and Jacob, D.: Parameterization of snow-free land surface albedo as a function of vegetation phenology based on MODIS data and applied in climate modelling, Theor. Appl. Climatol., 95, 245–255, 2009. a

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, Tech. Rep. 349, Max-Planck-Institut für Meteorologie, 2003. a, b

Roesch, A. and Roeckner, E.: Assessment of snow cover and surface albedo in the ECHAM5 general circulation model, J. Climate, 19, 3828–3843, 2006. a, b

Roesch, A., Wild, M., Gilgen, H., and Ohmura, A.: A new snow cover fraction parametrization for the ECHAM4 GCM, Clim. Dynam., 17, 933–946, 2001. a, b, c

Sage, R. F. and Kubien, D. S.: The temperature response of C3 and C4 photosynthesis, Plant Cell Environ., 30, 1086–1106, 2007. a, b

Sato, N., Sellers, P., Randall, D., Schneider, E., Shukla, J., Kinter III, J., Hou, Y., and Albertazzi, E.: Effects of implementing the simple biosphere model in a general circulation model, J. Atmos. Sci., 46, 2757–2782, 1989. a

Schaefer, K., Schwalm, C. R., Williams, C., et al.: A model-data comparison of gross primary productivity: Results from the North American Carbon Program site synthesis, J. Geophys. Res.-Biogeo., 117, G03010,, 2012. a, b, c

Schlesinger, W. H. and Jasechko, S.: Transpiration in the global water cycle, Agr. Forest Meteorol., 189, 115–117, 2014. a

Sellers, P., Bounoua, L., Collatz, G., Randall, D., Dazlich, D., Los, S., Berry, J., Fung, I., Tucker, C., Field, C., and Jensen, T.: Comparison of Radiative and Physiological Effects of Doubled Atmospheric CO2 on Climate, Science, 271, 1402–1406, 1996a. a

Sellers, P. J., Tucker, C. J., Collatz, G. J., Los, S. O., Justice, C. O., Dazlich, D. A., and Randall, D. A.: A revised land surface parameterization (SiB2) for atmospheric GCMs – Part II: The generation of global fields of terrestrial biophysical parameters from satellite data, J. Climate, 9, 706–737, 1996b. a, b

Sellers, P., Dickinson, R., Randall, D., Betts, A., Hall, F., Berry, J., Collatz, G., Denning, A., Mooney, H., Nobre, C., Sato, N., Field, C., and Henderson-Sellers, A.: Modeling the exchanges of energy, water, and carbon between continents and the atmosphere, Science, 275, 502–509, 1997. a

Serbin, S. P., Ahl, D. E., and Gower, S. T.: Spatial and temporal validation of the MODIS LAI and FPAR products across a boreal forest wildfire chronosequence, Remote Sens. Environ., 133, 71–84, 2013. a

Shinozaki, K., Yoda, K., Hozumi, K., and Kira, T.: A quantitative analysis of plant form-the pipe model theory: I. Basic analyses, Jpn. J. Ecology, 14, 97–105, 1964. a

Shuttleworth, W. J.: Evaporation from Amazonian rainforest, Proc. R. Soc. Lon. Ser.-A, 233, 321–346, 1988. a

Sitch, S., Smith, B., Prentice, I. C., Arneth, A., Bondeau, A., Cramer, W., Kaplan, J., Levis, S., Lucht, W., Sykes, M., Thonicke, K., and Venevsky, S.: Evaluation of ecosystem dynamics, plant geography and terrestrial carbon cycling in the LPJ dynamic global vegetation model, Glob. Change Biol., 9, 161–185, 2003. a

Sitch, S., Huntingford, C., Gedney, N., Levy, P., Lomas, M., Piao, S., Betts, R., Ciais, P., Cox, P., Friedlingstein, P., Jones, C., Prentice, I., and Woodward, F.: Evaluation of the terrestrial carbon cycle, future plant geography and climate-carbon cycle feedbacks using five Dynamic Global Vegetation Models (DGVMs), Glob. Change Biol., 14, 2015–2039, 2008. a, b

Sperry, J., Hacke, U., Oren, R., and Comstock, J.: Water deficits and hydraulic limits to leaf water supply, Plant Cell Environ., 25, 251–263, 2002. a

Terashima, I., Masuzawa, T., Ohba, H., and Yokoi, Y.: Is photosynthesis suppressed at higher elevations due to low CO2 pressure?, Ecology, 76, 2663–2668, 1995. a, b

Tian, H., Lu, C., Yang, J., Banger, K., Huntzinger, D. N., Schwalm, C. R., Michalak, A. M., Cook, R., Ciais, P., Hayes, D., Huang, M., Ito, A., Jain, A. K., Lei, H., Mao, J., Pan, S., Post, W. M., Peng, S., Poulter, B., Ren, W., Ricciuto, D., Schaefer, K., Shi, X., Tao, B., Wang, W., Wei, Y., Yang, Q., Zhang, B., and Zeng, N.: Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions, Global Biogeochem. Cy., 29, 775–792, 2015. a, b, c

Todd-Brown, K. E. O., Randerson, J. T., Post, W. M., Hoffman, F. M., Tarnocai, C., Schuur, E. A. G., and Allison, S. D.: Causes of variation in soil carbon simulations from CMIP5 Earth system models and comparison with observations, Biogeosciences, 10, 1717–1736,, 2013. a, b

Tyree, M. T. and Sperry, J. S.: Do woody plants operate near the point of catastrophic xylem dysfunction caused by dynamic water stress? Answers from a model, Plant Physiol., 88, 574–580, 1988. a

Vamborg, F. S. E., Brovkin, V., and Claussen, M.: The effect of a dynamic background albedo scheme on Sahel/Sahara precipitation during the mid-Holocene, Clim. Past, 7, 117–131,, 2011. a, b

Vamborg, F. S. E., Brovkin, V., and Claussen, M.: Background albedo dynamics improve simulated precipitation variability in the Sahel region, Earth Syst. Dynam., 5, 89–101,, 2014. a

van de Griend, A. A. and Owe, M.: Bare soil surface resistance to evaporation by vapor diffusion under semiarid conditions, Water Resour. Res., 30, 181–188, 1994. a, b, c

Van Oijen, M., Schapendonk, A., and Höglind, M.: On the relative magnitudes of photosynthesis, respiration, growth and carbon storage in vegetation, Ann. Bot.-London, 105, 739–797, 2010. a

Wang-Erlandsson, L., van der Ent, R. J., Gordon, L. J., and Savenije, H. H. G.: Contrasting roles of interception and transpiration in the hydrological cycle – Part 1: Temporal characteristics over land, Earth Syst. Dynam., 5, 441–469,, 2014. a, b

Way, D. A. and Yamori, W.: Thermal acclimation of photosynthesis: on the importance of adjusting our definitions and accounting for thermal acclimation of respiration, Photosynth. Res., 119, 89–100, 2014. a, b, c

Wieder, W., Boehnert, J., Bonan, G., and Langseth, M.: Regridded Harmonized World Soil Database v1.2, Data set, Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, USA, (last access: February 2016), 2011. a, b, c, d, e

Willeit, M. and Ganopolski, A.: PALADYN v1.0, a comprehensive land surface–vegetation–carbon cycle model of intermediate complexity, Geosci. Model Dev., 9, 3817–3857,, 2016 a, b

Williams, M., Rastetter, E. B., Van der Pol, L., and Shaver, G. R.: Arctic canopy photosynthetic efficiency enhanced under diffuse light, linked to a reduction in the fraction of the canopy in deep shade, New Phytol., 202, 1267–1276, 2014. a, b

Williamson, M., Lenton, T., Shepherd, J., and Edwards, N.: An efficient numerical terrestrial scheme (ENTS) for Earth system modelling, Ecol. Model., 198, 362–374, 2006. a, b, c, d, e, f, g, h, i

Wu, J., Albert, L. P., Lopes, A. P., Restrepo-Coupe, N., Hayek, M., Wiedemann, K. T., Guan, K., Stark, S. C., Christoffersen, B., Prohaska, N., Tavares, J. V., Marostica, S., Kobayashi, H., Ferreira, M. L., Campos, K. S., da Silva, R., Brando, P. M., Dye, D. G., Huxman, T. E., Huete, A. R., Nelson, B. W., and Saleska, S. R.: Leaf development and demography explain photosynthetic seasonality in Amazon evergreen forests, Science, 351, 972–976, 2016.  a, b, c

Wu, X., Ju, W., Zhou, Y., et al.: Performance of linear and nonlinear two-leaf light use efficiency models at different temporal scales, Remote Sens., 7, 2238–2278, 2015. a

Yamori, W., Hikosaka, K., and Way, D. A.: Temperature response of photosynthesis in C3, C4, and CAM plants: temperature acclimation and temperature adaptation, Photosynth. Res., 119, 101–117, 2014. a, b, c, d, e, f, g

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, b, c

Yuan, W., Liu, S., Zhou, G., Zhou, G., Tieszen, L. L., Baldocchi, D., Bernhofer, C., Gholz, H., Goldstein, A. H., Goulden, M. L., Hollinger, D. Y., Hu, Y., Law, B. E., Stoy, P. C., Vesala, T., Wofsy, S. C., and other AmeriFlux collaborators: Deriving a light use efficiency model from eddy covariance flux data for predicting daily gross primary production across biomes, Agr. Forest Meteorol., 143, 189–207, 2007. a, b

Yuan, W., Cai, W., Xia, J., Chen, J., Liu, S., Dong, W., Merbold, L., Law, B., Arain, A., Beringer, J., Bernhofer, C., Black, A., Blanken, P. D., Cescatti, A., Chen, Y., Francois, L., Gianelle, D., Janssens, I. A., Jung, M., Kato, T., Kiely, G., Liu, D., Marcolla, B., Montagnani, L., Raschi, A., Roupsard, O., Varlagin, A., and Wohlfahrt, G.: Global comparison of light use efficiency models for simulating terrestrial vegetation gross primary production based on the LaThuile database, Agr. Forest Meteorol., 192, 108–120, 2014. a

Zhang, X., Liang, S., Wang, G., Yao, Y., Jiang, B., and Cheng, J.: Evaluation of the Reanalysis Surface Incident Shortwave Radiation Products from NCEP, ECMWF, GSFC, and JMA Using Satellite and Surface Observations, Remote Sens., 8, 225,, 2016. a

Zhou, J., Poulsen, C. J., Rosenbloom, N., Shields, C., and Briegleb, B.: Vegetation-climate interactions in the warm mid-Cretaceous, Clim. Past, 8, 565–576,, 2012. a


Canopy resistance and conductance are aggregated or “bulk” versions of, respectively, stomatal resistance and conductance, such that the entire canopy can be treated like a big leaf (e.g., see Jarvis and McNaughton1986; Raupach and Finnigan1988)


While spherical leaf angle distribution is commonly assumed for unmeasured canopies, measurements of 58 temperate and boreal broadleaf tree species showed planophile (tending toward horizontal) leaf angle distributions in 30 and spherical distributions in just 5 (Pisek et al.2013). Moreover, a modeling study (Hikosaka and Hirose1997) used game theoretical arguments to show how a horizontal leaf angle distribution is to be evolutionarily expected, in general (save for canopies with high LAI and high self-shading). Thus, globally applying a horizontal canopy leaf distribution is not unreasonable.


In this paper, we define ET and transpiration (T) as positive when there is a net water vapor flux from surface to atmosphere. In the actual SEDGES code, ET and T are defined as negative when the flux is from surface to atmosphere. In our exposition, we have reversed the sign convention to aid interpretation of the model equations.


An important requirement for the aerodynamic conductance formulation is that the surface roughness for moisture be the same as that for momentum. See Sect. 2.3.3 for more discussion.


The value of 0.12 is a crude estimate that is obtained by assuming an albedo transition for forest vegetation from snow-free vegetation albedo to deep and cold snow albedo (0.80) that is linear with fsnow for. Depending on whether one uses values of 0.30 and 0.12 for the values of snow-covered forest and snow-free vegetation, respectively, or values taken from Betts and Ball (1997), one ends up with fsnow for values ranging from  0.04 to 0.26.


Roesch and Roeckner (2006) include the 0.95 multiplier in their definition of the snow-covered fraction of the non-forested part of the grid cell, but we do not. Note, too, that this scheme for non-forested snow cover fraction supplants the ECHAM4 parameterization fsnowflat=sweswe+0.01 used by SimBA, which was found to poorly estimate satellite-based observations compared to the adopted ECHAM5 parameterization (Roesch et al.2001).


Visually estimated from their Fig. 9.


In addition to the direct impact of CO2 on GPP, some vegetation feedback onto GPP via climate is included here; also, vegetation distributions and LAI are held fixed.


The original MODIS data were processed by spatially interpolating to fill in missing data, taking multiyear means (i.e., means across multiple years) of the 23 16-day periods for each year, spatially interpolating these into T62 resolution, temporally interpolating the 16-day periods linearly to obtain 24 new such periods such that each month was covered by two periods, taking the mean of the two periods for each month, masking out albedo values for each month and region that had polar-night-caused missing values in its original pre-interpolated comprising data, and masking out water and land ice grid points in the ERA-Interim T62 data that were used to force SEDGES. White-sky albedo from MODIS was used as the reference dataset because SEDGES's albedo formulation lacks dependency on the zenith angle of the direct beam.


Even though these areas would be mostly grasslands rather than cropland if not for human land cover change (Ramankutty and Foley1999), the winter albedo would be similarly high in either case.


Here, we combine in situ precipitation and ET data from Da Rocha et al. (2004) with a ground-observed interception loss-to-precipitation ratio of 0.116 from nearby (as reported by Czikowsky and Fitzjarrald2009) to get an interception loss-to-ET ratio of 0.22.

Short summary
This paper presents a simple vegetation model for use as part of a climate or Earth system model. The model equations and their derivations are presented. We evaluate the model’s performance offline using near-present-day real-world datasets and deem it to be satisfactory. The model is useful because it is fast, easy to understand, and general in its formulations. It was developed to better simulate climate–vegetation feedbacks for answering paleoclimate research questions.