Integrated modeling of canopy photosynthesis, fluorescence, and the transfer of energy, mass, and momentum in the soil–plant–atmosphere continuum (STEMMUS–SCOPE v1.0.0)

Root water uptake by plants is a vital process that influences terrestrial energy, water, and carbon exchanges. At the soil, vegetation, and atmosphere interfaces, root water uptake and solar radiation predominantly regulate the dynamics and health of vegetation growth, which can be remotely monitored by satellites, using the soil–plant relationship proxy – solar-induced chlorophyll fluorescence. However, most current canopy photosynthesis and fluorescence models do not account for root water uptake, which compromises their applications under water-stressed conditions. To address this limitation, this study integrated photosynthesis, fluorescence emission, and transfer of energy, mass, and momentum in the soil–plant–atmosphere continuum system, via a simplified 1D root growth model and a resistance scheme linking soil, roots, leaves, and the atmosphere. The coupled model was evaluated with field measurements of maize and grass canopies. The results indicated that the simulation of land surface fluxes was significantly improved by the coupled model, especially when the canopy experienced moderate water stress. This finding highlights the importance of enhanced soil heat and moisture transfer, as well as dynamic root growth, on simulating ecosystem functioning.


Introduction
Root water uptake (RWU) by plants is a critical process controlling water and energy exchanges between the land surface and the atmosphere and, as a result, plant growth. The representation of RWU is an essential component of ecohydrological models that simulate terrestrial water, energy, and carbon fluxes (Seneviratne et al., 2010;Wang and Smith, 2004). However, most of these models consider the aboveground processes in much greater detail than belowground processes; therefore, they have a limited ability to represent the dynamic response of plant water uptake to water stress. A particular mechanism of importance for plants to mitigate water stress is the compensatory root water uptake (CRWU) which refers to the process by which water uptake from sparsely rooted but well-watered parts of the root zone compensates for stress in other parts (Jarvis, 2011). The failure to account for compensatory water uptake and the associated hydraulic lift from deep subsoil (Caldwell et al., 1998;Espeleta et al., 2004;Amenu and Kumar, 2007;Fu et al., 2016) can lead to significant uncertainties in simulating the plant growth and corresponding ecohydrological processes (Seneviratne et al., 2010).
Because the spatial (i.e., 1D vertical) pattern of RWU is determined by the spatial distribution of the root system, knowledge of the latter is essential for predicting the spatial distribution of water contents and water fluxes in soils. The distribution of roots and their growth are, in turn, sensitive to various physical, chemical, and biological factors, as well as to soil hydraulic properties that influence the availability of water for plants (Beaudoin et al., 2009). Many attempts have been made in the past to develop root growth models that account for the influence of various environmental factors such as temperature, aeration, soil water availability, and soil compaction. Existing root growth models range from complex, 3D root architecture models (Bingham and Wu, 2011;Leitner et al., 2010;Wu et al., 2005) to much simpler root growth models that are implemented within more complex models such as EPIC (Williams et al., 1989) and DSSAT (Robertson et al., 1993). Most of these models reproduce the measured rooting depth very well, but the distribution of new growth root is based on empirical functions rather than biophysical processes (Camargo and Kemanian, 2016; Table 1).
Modeling RWU requires the representation of above-and belowground processes, which can be realized considering the flow of water from soil through the plant to the atmosphere (i.e., the soil-plant-atmosphere continuum, SPAC model; Guo, 1992). The SPAC model represents a good compromise between simplicity (i.e., a small number of tuning parameters) and the ability to capture non-linear responses of RWU (and subsequently the ecosystem functioning) to drought events. Specifically, the SPAC model calculates the CRWU term using the gradient between the leaf water potential and the soil water potential of each soil layer. The most important parameters in the SPAC model include the leaf water potential, stomatal resistance, and the root resistance. Different from other macroscopic models using the root distribution function, the SPAC model explicitly needs the root length density at each soil layer to calculate the root resistance for each soil layer (Deng et al., 2017). The most practical method for obtaining the root length density is using a root growth model.
On other hand, remote sensing of solar-induced chlorophyll fluorescence (SIF) has been deployed to understand and monitor the ecosystem functioning under drought stress using models for vegetation photosynthesis and fluorescence (Zhang et al., 2018(Zhang et al., , 2020Mohammed et al., 2019;Shan et al., 2019). SCOPE (Soil Canopy Observation, Photochemistry, and Energy Fluxes) is one such model and simulates canopy reflectance and fluorescence spectra in the observation directions as well as photosynthesis and evapotranspiration as functions of leaf optical properties, canopy structure, and weather variables (Van der Tol et al., 2009). The SCOPE model provides a valuable means to study the link between remote sensing signals and ecosystem functioning; however, it does not consider the water budget in soil and vegetation.
As such, there is no explicit parametrization of the effects of soil moisture variations on the photosynthetic or stomatal parameters. Consequently, soil moisture effects are only "visible" in SCOPE if the lack of soil moisture affects the optical or thermal remote sensing signals (i.e., during water stress periods). The lack of such a link between soil moisture availability and remote sensing signals compromises the capacity of SCOPE to simulate and predict drought events on vegetation functioning.
The change in vegetation optical appearance as a result of soil moisture variations can only partially explain the soil moisture effect on ecosystem functioning (Bayat et al., 2018), which leads to considerably biased estimations of the gross primary productivity (GPP) and evapotranspiration (ET) under water-limited conditions. This presents a challenge with respect to using SCOPE for ecosystems in arid and semiarid areas, where water availability is the primary limiting factor for vegetation functioning. This challenge becomes even more relevant considering that soil moisture deficit or "ecological drought" is expected to increase in both frequency and severity in nearly all ecosystems around the world (Zhou et al., 2013). Bayat et al. (2019) incorporated the SPAC model into SCOPE to address water-stressed conditions at a grassland site, but the coupled model neglected the dynamic root distribution in different soil layers, and soil moisture only serves as a model input when it comes from measurements.
In this study, the modeling of aboveground photosynthesis, fluorescence emission, and energy fluxes in the vegetation layer by SCOPE will be fully coupled with a two-phase mass and heat transfer model -the STEMMUS model (Simultaneous Transfer of Energy, Mass and Momentum in Unsaturated Soil; a more detailed description of STEMMUS can be found in Sect. 2), by considering RWU based on a root growth model. The root growth model and the corresponding resistance scheme (from soil, through roots and leaves, to atmosphere) will be integrated for the dynamic modeling of water stress and the root system, enabling the seamless modeling of soil-water-plant energy, water, and carbon exchanges as well as SIF, thereby directly linking the vegetation dynamics (and its optical and thermal appearance) at the process level to soil moisture variability.
The rest of this article is structured as follows: Sect. 2 describes the coupling scheme between SCOPE and STEM-MUS and the data that were used to validate the coupled model; Sect. 3 verifies the coupled STEMMUS-SCOPE model using a maize agroecosystem and a grassland ecosystem located in semiarid regions and explores the dynamic responses of the leaf water potential and root length density to water stress; and the summary of this study and the further challenges are addressed in Sect. 4. Table 1. Comparison of land surface models (LSMs) and crop models in terms of sink term calculation of soil water balance. CRWU stands for compensatory root water uptake.

Model
Sink term calculation of soil water balance Root water uptake process Hydraulic redistribution (Richards and Caldwell, 1987) Compensatory uptake (Jarvis, 2011) Root distribution LSMs CLM5.0 Root length density of each soil layer; water stress is applied by the hydraulic conductance model (Lawrence et al., 2020) Extreme case of CRWU Following Darcy's law for porous media flow equations Empirical function depends on the plant functional type CLM4.5 Actual transpiration, root fraction of each soil layer, and soil integral soil water availability  The Ryel et al. (2002) function

Not considered
Empirical function CLM4.0 Actual transpiration, root fraction of each soil layer, and integral soil water availability (Couvreur et al., 2012, Sulis et al., 2019 HRWU scheme (RWU model based on hydraulic architecture)

HRWU scheme
Empirical function

CLM3 & IBIS2
Actual transpiration, physical root distribution, and the water availability in each layer (Zheng and Wang, 2007) The Ryel et al. (2002) function Dynamic root water uptake Empirical function

CoLM
Potential transpiration, root fraction in each layer, and water stress factor (Zhu et al., 2017) The Ryel et al. (2002) and the Amenu and Kumar (2007) function Empirical approach with a compensatory factor Empirical function JULES Potential transpiration, root fraction of each soil layer, and a weighted water stress in each layer (Eller et al., 2020) Not considered Not considered Exponential distribution with depth

Noah-MP
Based on the gradient in water potentials between root and soil, and root surface area (Niu et al., 2020) Extreme case of CRWU Following Darcy's law for porous media flow equations Process-based 1D root surface area growth model Table 1. Continued.

Model
Sink term calculation of soil water balance Root water uptake process Hydraulic redistribution (Richards and Caldwell, 1987) Compensatory uptake (Jarvis, 2011) Root distribution CABLE Based on the gradient in water potentials between the leaf, stem, and the weighted average of the soil (De Kauwe et al., 2020) Extreme case of CRWU Following Darcy's law for porous media flow equations

Empirical function
Crop models APSIM Potential transpiration and water supply factor but neglects root distribution (Keating et al., 2003) Not considered Not considered Empirical function CropSyst Difference in water potential between the soil and the leaf, and a total soil-root-shoot conductance (Stöckle et al., 2003) Not considered Considered by the leaf and soil water potential Linear decrease in soils with no limitations on root exploration DSSAT Water uptake per unit of root length is computed as an exponential function, and the actual RWU is the minimum of potential transpiration and the maximum capacity of root water uptake (Jones et al., 2003) Not considered Water uptake per unit of root length as a function of soil moisture Using an empirical function EPIC EPIC assumes that water is used preferentially from the top layers, and the potential water supply rate decreases exponentially downward (Williams et al., 2014) Not considered Not considered Not considered SWAP Based on the potential transpiration, root fraction, and an empiric stress factor relationship (van Dam, 2000) Not considered Based on soil water potential Function of relative rooting depth WOFOST The simplest one, it calculates water uptake as a function of the rooting depth and the water available at that rooting depth without regard for the soil water distribution with depth (Supit et al., 1994) Not considered Not considered Empirical function Table 1. Continued.

Model
Sink term calculation of soil water balance Root water uptake process Hydraulic redistribution (Richards and Caldwell, 1987) Compensatory uptake (Jarvis, 2011) Root distribution SPACSYS According to empirical root length density distribution in a soil layer, potential transpiration, and soil moisture (Wu et al., 2005) Not considered Not considered 1D (empirical function) or 3D root system (process based)

STICS
Based on the potential transpiration, root fraction, and soil water distribution, but not process based (Beaudoin et al., 2009) Not considered Not considered 1D root length density profile 2 Methodology and data 2.1 SCOPE and SCOPE_SM models SCOPE is a radiative transfer and energy balance model (Van der Tol et al., 2009). It simulates the transfer of optical, thermal, and fluorescent radiation in the vegetation canopy and computes ET using an energy balance routine. SCOPE includes a radiative transfer module for incident solar and sky radiation to calculate the top-of-canopy outgoing radiation spectrum, net radiation, and absorbed photosynthetically active radiation (aPAR); a radiative transfer module for thermal radiation emitted by soil and vegetation to calculate the top-of-canopy outgoing thermal radiation and net radiation; an energy balance module for latent heat, sensible heat, and soil heat flux; and a radiative module for chlorophyll fluorescence to calculate the top-of-canopy SIF (the observation zenith angle was set as 0 • in this study). Compared with other radiative transfer models that simplify the radiative transfer processes based on Beer's law, SCOPE has well-developed radiative transfer modules that consider the various leaf orientation and multiple scattering. SCOPE can provide detailed information about the net radiation of every leaf within the canopy. Furthermore, SCOPE incorporates an energy balance model that predicts not only the temperature of leaf but also the soil surface temperature (i.e., a vital boundary condition needed by STEM-MUS). In the original SCOPE model, soil is treated in a very simple way with several empirical functions describing the ground heat storage. Later, Bayat et al. (2019) extended the SCOPE model by including the moisture effects on the vegetation canopy, which resulted in the SCOPE_SM model. This model takes soil moisture as input and predicts the effects on several processes of the vegetation canopy using the SPAC concept. Appendix A1 lists the main equations for calculating the water stress factor within SCOPE (Bayat et al., 2019), and the reader is referred to Van der Tol et al. (2009) for a detailed formulation of SCOPE.
SCOPE_SM provides the basic framework to couple SCOPE with a soil process model. However, both SCOPE and SCOPE_SM ignored the soil heat and mass transfer processes and the dynamics of root growth. This can be overcome by introducing the STEMMUS model.

STEMMUS model
The STEMMUS model is a two-phase mass and heat transfer model with explicit consideration of the coupled liquid, vapor, dry air, and heat transfer in unsaturated soil (Zeng et al., 2011a, b;Zeng and Su, 2013;Yu et al., 2016Yu et al., , 2018. STEMMUS provides a comprehensive description of water and heat transfer in the unsaturated soil, which can compensate for what is currently neglected in SCOPE. In STEM-MUS, the soil layers can be set in a flexible manner, which is an improvement on the previous SPAC model that only considered the whole root zone soil water content as fixed layers (Williams et al., 1996). The water and heat transfer processes are vital for vegetation phenology development as well as freeze-thaw processes. The boundary condition needed by STEMMUS includes surface soil temperature, which is the output of SCOPE. In addition, STEMMUS already contains an empirical equation to calculate root water uptake and a simplified root growth module to calculate root fraction profile. As such, STEMMUS has an ideal model structure to be coupled with SCOPE. The main governing equations of STEMMUS are listed in Appendix A2.

Dynamic root growth and root water uptake
To obtain the root resistance of each soil layer, we incorporated a root growth module to simulate the root length density profile (see Appendix A3). The simulation of root growth refers to the root growth module in the INRA STICS crop growth model (Beaudoin et al., 2009), which includes the calculations of root front growth and root length growth. The root front growth is a function of temperature, with the depth of the root front beginning at the sowing depth for sown crops and at an initial value for transplanted crops or perennial crops (Beaudoin et al., 2009). The root length growth is calculated in each soil layer, considering the net assimilation rate and the allocation fraction of net assimilation to root, which is, in turn, a function of leaf area index (LAI) and root zone water content (Krinner et al., 2005). The root length density profile is then used to calculate the root resistance to water flow radially across the roots, soil hydraulic resistance, and plant axial resistance to flow from the soil to the leaves (see Appendix A4).

STEMMUS-SCOPE v1.0.0 coupling
The coupling starts with an initial soil moisture (SM) profile simulated by STEMMUS, which enables the calculation of the water stress factor as a reduction factor of the maximum carboxylation rate (V cmax ). SCOPE v1.73 is then used to calculate net photosynthesis (A n ) or gross primary productivity (GPP), soil respiration (R s ), energy fluxes (net radiation, R n ; latent heat, LE; sensible heat, H ; and soil heat flux, G), transpiration (T ), and SIF, which is passed to STEMMUS as the root water uptake (RWU). Then, the gross primary production (GPP) can be calculated based on A n . Surface soil moisture is also used in calculating soil surface resistance and then calculating soil evaporation (E) . Furthermore, SCOPE can calculate soil surface temperature (T s0 ) based on energy balance, which is subsequently used as the top boundary condition of STEMMUS, and leaf water potential (LWP), which is a parameter to reflect plant water status, can be calculated through iteration. Based on RWU, STEMMUS calculates the soil moisture in each layer at the end of the time step, and the new soil moisture profile will be the soil moisture at the beginning of next time step, which is repeated as such until the end of simulation period. The time step of STEMMUS-SCOPE is flexible, and the time step used in this study was 30 min. Figure 1 shows the coupling scheme of STEMMUS and SCOPE, and Table B1 in the Appendix shows all of the parameter values used in this study.

Evapotranspiration partitioning
Most studies in partitioning evapotranspiration (ET) use sap flow and microlysimeter data from in situ measurements. In this study, we used a simple and practical method to separate evaporation (E) and transpiration (T ) proposed by Zhou et al. (2016). Although the behavior of plant stomata is influenced by environmental factors, the potential water use efficiency (uWUE p , g C hPa 0.5 (kg H 2 O) −1 ) at the stomatal scale in the ecosystem with a homogeneous underlying surface is assumed to be nearly constant, and variations in actual uWUE (g C hPa 0.5 (kg H 2 O) −1 ) can be attributed to the soil evaporation (Zhou et al., 2016). Thus, the method can be used to estimate T and E with the quantities of ET, uWUE, and uWUE p . Another assumption of this method is that the ecosystem T is equal to ET at some growth stages, so uWUE p can be estimated using the upper bound of the ratio of GPP √ VPD to ET (here VPD refers to the vapor pressure deficit; Zhou et al., 2014Zhou et al., , 2016. Zhou et al. (2016) used the 95th quantile regression between GPP

√
VPD and ET to estimate uWUE p , and they showed that the 95th quantile regression for uWUE p at flux tower sites was consistent with the uWUE derived at the leaf scale for different ecosystems. In addition, the variability in seasonal and interannual uWUE p was relatively small for a homogeneous canopy. Therefore, the calculations of uWUE p , uWUE, and T at the ecosystem scale were as follows: The calculation of the VPD was based on air temperature and relative humidity data, and the method of gap-filling was the marginal distribution sampling (MDS) method proposed by Reichstein et al. (2005). To calculate GPP, the complete series of net ecosystem exchange (NEE) was partitioned into gross primary production (GPP) and respiration (Re) using the method proposed by Reichstein et al. (2005). Finally, ET was calculated using the latent heat flux and air temperature. Based on GPP, ET, and VPD data, T can be calculated using the method proposed by Zhou et al. (2016).

Study site and data description
To evaluate the performance of STEMMUS-SCOPE in modeling ecohydrological processes, simulation was conducted to compare STEMMUS-SCOPE with SCOPE, SCOPE_SM, and STEMMUS using observations over a C 4 cropland (summer maize: from 11 June to 10 October 2017) at the Yangling station (34 • 17 N, 108 • 04 E; 521 m a.s.l.) and a C 3 grassland at the Vaira Ranch (US-Var) FLUXNET site (38 • 25 N, 120 • 57 W; 129 m a.s.l.; annual grasses: from 1 June to 8 August 2004). The seasonal variation in precipitation, irrigation, and SM for these two sites are presented in Fig. 2, and the differences in soil surface resistance, water stress factor (WSF), ET, photosynthesis, soil surface temperature (T s0 ), root water  Table B1 in the Appendix. uptake (RWU), and leaf water potential (LWP) between these four models are presented in Table 2. In this study, the LAI data of the Vaira Ranch (US-Var) FLUXNET site were from the MODIS 8 d LAI product instead of the field-measured LAI used by Bayat et al. (2019). For the soil water content employed by SCOPE_SM, the averaged root zone soil moisture was used for Yangling station, and the soil moisture at 10 cm depth was used for the Vaira Ranch site. For more detailed descriptions of these sites and data, the reader is referred to Wang et al. (2019Wang et al. ( , 2020a and Bayat et al. (2018Bayat et al. ( , 2019.

Performance metrics
The metrics used to evaluate the performance of the coupled STEMMUS-SCOPE model include the (1) root-meansquare error (RMSE), (2) coefficient of determination (R 2 ), and (3) the index of agreement (d). They are calculated as follows: where P i is the ith predicted value, O i is the ith observed value,Ō is the average of the observed values, and n is the number of samples. The simulated soil moisture at 20 cm depth agreed with the observed values in terms of the seasonal pattern. Although a slight overestimation occurred at initial and late stages, the dynamics in soil moisture resulting from precipitation or irrigation were well captured. Per the nature of the two models, the coupling of SCOPE with STEMMUS is not expected to improve the simulation of soil moisture. However, compared with SCOPE_SM, which used soil moisture measurements as inputs, the coupled STEMMUS-SCOPE model improves the simulation of soil moisture dynamics as measured. The deviation between the model simulations and the measurements can be attributed to the following two potential reasons. First, the field observations contain errors to a certain extent, and the soil moisture sensors may be not well calibrated. Second, in this simulation, we assumed that the soil texture was homogeneous in the vertical profile, whereas, in reality, the soil properties (e.g., soil bulk density and saturated hydraulic conductivity) may vary with depth and at different growth stages due to field management practices. For example, the soil bulk density at 40 cm was much higher than that at 20 cm due to the mechanical tillage, especially in the early stage.

Soil temperature modeling
Similar to soil moisture, only soil temperatures (T s ) simulated by STEMMUS and STEMMUS-SCOPE at 20 and 40 cm depth at the Yangling site are shown in Fig. 4. In general, both models can capture the dynamics of soil temperature well. For the simulation of temperature at 20 cm, the RMSE values were 2.56 and 2.58 • C and the d values were 0.92 and 0.92 for STEMMUS and STEMMUS-SCOPE, respectively. For the simulation of temperature at 40 cm, the RMSE values were 2.06 and 2.07 • C and the d values were 0.93 and 0.93, respectively. These results indicate that both models can simulate soil temperature well. However, some differences also exist between the simulation and observations. The largest difference occurred on DOY (day of year) 202, when the field was irrigated using the flood irrigation method. This irrigation activity may lead to boundary condition errors (i.e., for soil surface temperature), which cannot be estimated well enough (e.g., there is no monitoring of water temperature from the irrigation). Meanwhile, the measurements may also have some errors during this period. The fact that the observed soil temperature at 20 and 40 cm decreased to almost the same level at the same time indicates a potential pathway for preferential flow in the field (see precipitation and irrigation on DOY 202 in Fig. 2), and the sensors captured this phenomenon. Nevertheless, the model captures the soil temperature dynamics.

Energy balance modeling
Comparisons of the modeled and observed 30 min net radiation (R n ), sensible heat flux (H ), latent heat flux (LE), and soil heat flux (G) using SCOPE, SCOPE_SM, and STEMMUS-SCOPE are presented in Fig. 5      predict ET with a relatively higher accuracy, especially when the maize was under water stress (DOY 183-202 at Yangling station and DOY 90-220 at the US-Var site), and STEMMUS-SCOPE and SCOPE_SM performed similarly well. It is noteworthy that although STEMMUS considered the effect of soil moisture on ET, the accuracy of STEM-MUS was lower than that of the coupled model (see Fig. 6). The possible reason for this is the better representation of transpiration in the SCOPE model (see Fig. 7), which separates the canopy into 60 layers, whereas STEMMUS only treats the canopy as one layer. Moreover, the coupled model performed better for the grassland than for maize cropland.
The reason for this is that the grassland simulation used the dynamic V cmax data, whereas the maize simulation used a constant V cmax data. The modeled and observed daily transpiration at the maize cropland are presented in Fig. 7, and the modeled transpiration at the grassland site is presented in Fig. 8. For Yangling station, the R 2 values between the simulated and observed transpiration were 0.82, 0.86, 0.79, and 0.86, and the RMSEs were 0.60, 0.50, 0.67, and 0.50 mm d −1 , for SCOPE, SCOPE_SM, STEMMUS, and STEMMUS-SCOPE, respectively. Because it ignored the effect of water stress on transpiration, SCOPE failed to simulate transpiration accurately when the vegetation experienced water stress. As shown in Fig. 6a, SCOPE overestimated transpiration for the maize cropland at Yangling station from DOY 183 to 202 during the water stress period. Compared with SCOPE, SCOPE_SM, STEMMUS, and STEMMUS-SCOPE can capture the reduction in transpiration during the dry period. The performance of STEMMUS-SCOPE and SCOPE_SM was also better than that of STEMMUS. The possible reason for this is the more processed-based consideration of the radiative transfer and energy balance at the leaf level in the coupled STEMMUS-SCOPE model (as in SCOPE_SM) and the more accurate root water uptake (compared with that in SCOPE_SM). Nevertheless, STEMMUS-SCOPE slightly underestimated transpiration when the plant was undergoing severe water stress and slightly overestimated it after the field was irrigated. This is mainly because the actual V cmax was not only influenced by drought but was also related to the leaf nitrogen content (Xu and Baldocchi, 2003), which was not considered in the maize cropland simulation. Although measured T at the grassland was not available, we compared modeled T from the four models (Fig. 7). During the wet season (before DOY 85), the modeled T values from SCOPE, SCOPE_SM, and STEMMUS-SCOPE were similar and were higher than that from STEMMUS from DOY 64 to 82. During the dry season (after DOY 85), due to the simplified consideration of soil processes, the modeled T values from SCOPE and SCOPE_SM were both much higher than those from STEMMUS and STEMMUS-SCOPE. The reason for the better performance of the coupled model for the grassland (Fig. 6b) is that it also considers the effect of the leaf chlorophyll content (C ab ) on V cmax , in addition to a more detailed consideration of water stress as discussed above for the maize cropland.
As shown in Fig. 9 for soil evaporation at Yangling station, the simulated values from STEMMUS-SCOPE are closer to the observations than those from other models. When using  SCOPE to simulate soil evaporation, the soil moisture is set as constant (i.e., 0.25 m 3 m −3 ). Therefore, SCOPE generally underestimates soil evaporation when soil moisture is higher than 0.25 and overestimates it when it is lower than 0.25. Here, we use the average soil moisture at the root zone simulated by STEMMUS-SCOPE as the input data for SCOPE and SCOPE_SM in order to calculate soil surface resistance and soil evaporation. Although STEMMUS can capture variation in soil evaporation reasonably well, it has a higher RMSE than STEMMUS-SCOPE. This is probably attributed to the comprehensive consideration of radiation transfer in SCOPE, which is lacking in STEMMUS. Consequently, the simulation of soil net radiation by the coupled model was more accurate than that from STEMMUS alone. The RMSE value for STEMMUS-SCOPE was 0.60 mm d −1 , which was lower than those from the other three models (0.67, 0.65, and 0.64 mm d −1 , respectively). For STEMMUS-SCOPE, the major differences between simulations and observations occurred on rainy or irrigation days (see Fig. 2a), which may be caused by errors in the estimated soil surface resistance during these periods or the uncertainty of the ET partitioning method. The uncertainty of the ET partitioning method (Zhou et al., 2016) was mainly caused by (1) the uncertainty in the partitioning of GPP (less than 10 %) and Re based on NEE, which would result in some uncertainty in uWUE; (2) due to the seasonal variation in the atmospheric CO 2 concentration -the assumption of uWUE p being constant would cause some uncertainty (less than 3 %); (3) the assumption of T being equal to ET sometimes during the growing season would cause some uncertainty when vegetation is sparse. Because the observed E at the US-Var site was not available, a comparison of only modeled E is shown in Fig. 8, in which SCOPE modeled unrealistic E during the dry season, whereas the modeled E values from SCOPE_SM, STEM-MUS, and STEMMUS-SCOPE were consistent due to the use the simulated surface SM as the input for soil evaporation calculation.

Daily GPP modeling
Simulated GPP from SCOPE, SCOPE_SM, and STEMMUS-SCOPE and observed GPP are presented in Fig. 10. As shown, similar to the simulation of transpiration, SCOPE cannot respond to water stress when simulating GPP. After introducing a soil water stress factor in STEMMUS-SCOPE and SCOPE_SM, the simulation of GPP was improved in both models. For Yangling station, the consistency between simulated and observed GPP at mid and late stages was higher than that at early and rapid growth stages. The difference usually occurred when soil moisture increased. For the US-Var site, STEMMUS-SCOPE simulated GPP well during the whole period, whereas SCOPE_SM slightly underestimated GPP around DOY 80 when this site transits from the wet season to the dry season. This indicates that only using the surface SM cannot reflect the actual root zone SM when the vegetation is experiencing moderate water stress. Under such conditions, the hydraulic redistribution (HR) and compensatory root water uptake (CRWU) process enable the vegetation to utilize the water in the deep soil layer. Only using the surface soil water content to calculate RWU in SCOPE_SM ignored the effect of the HR and CRWU process, and the effect of water stress was overestimated. However, the surface soil moisture can reflect root zone soil moisture well when the vegetation is not under water stress or severe water stress. A similar underestimation of GPP was also found by Bayat et al. (2019).
3.6 Simulation of leaf water potential (LWP), water stress factor (WSF), and root length density (RLD) The simulated 30 min leaf water potential and water stress factor at Yangling station are presented in Fig. 11. The leaf water potential was lower when vegetation was suffering water stress compared with other periods. The reason for this is that soil water potential is low due to the low soil moisture, and plants need to maintain an even lower leaf water poten- tial to suck water from the soil and transfer it to leaves. During mid and late stages, the leaf water potential was sensitive to transpiration demand due to the slowdown of root system growth. As continuous measurements of the leaf water potential are not available, we compared only the magnitude of simulated leaf water potential to measurements reported in the literature. Many studies have measured midday leaf water potential or dawn leaf water potential. Fan et al. (2015) reported that the leaf water potential of well-watered maize remained high at between −73 and −88 m and that leaf water potential would decrease when the soil water content was lower than 80 % of field capacity. Martineau et al. (2017) reported that the midday leaf water potential of well-watered maize was around −0.82 MPa (about 84.8 m in water pressure head; note that 0.1 MPa is equal to 10.339 m water pressure head) and that the midday leaf water potential decreased to −1.3 MPa (about 134.4 m in water head) when the maize was suffering water stress. Moreover, O'Toole and Cruz (1980) studied the response of leaf water potential to water stress in rice and concluded that the leaf water potential of rice can be lower than −80 to −120 m when the vegetation was under water stress and the leaves started curling, which was similar to the simulated leaf water potential of maize in this study. Aston and Lawlor (1979) revealed the relationship between transpiration, root water uptake, and leaf water potential of maize. These field studies found that leaf water potential was often very low and reached trough values at midday. Elfving et al. (1972) developed a water flux model based on the SPAC system, evaluated it for orange trees, and reported about −120 m for the trough value of leaf water potential under non-limiting environmental conditions, which was slightly lower than the simulation in this study.
In this study, the calculation of the water stress factor considered the effect of soil moisture and root distribution. The severe water stress occurred from DOY 183 to 202, and the coupled model performed very well during this period. Due to feedback, water stress can also influence root water uptake and root growth and, consequently, influence soil moisture and root dynamics in next time step. This indicates that the water stress equation used in this study can characterize the reduction in V cmax reasonably well. Root length density is another vital parameter in calculating root water uptake. As shown in Table 3, the simulated peak root length density and maximum rooting depth of maize at Yangling station was comparable to the measured values from other sites. Many previous studies have revealed that root length density is influenced by soil moisture, bulk density, tillage, and soil mineral nitrogen (Amato and Ritchie, 2002; Chassot et al., 2001;Schroder et al., 1996). In this study, as we assumed that the soil was homogenous. STEMMUS-SCOPE considered the effect of soil moisture but neglected the effect of bulk density and soil mineral nitrogen. Amato and Ritchie (2002) also found a similar result to this study with respect to the root length density in a maize field. Peng et al. (2012) studied temporal and spatial dynamics in the root length density of fieldgrown maize and found that 80 % root length density was distributed at 0-30 cm depth with peak values from 0.86 to 1.00 cm cm −3 . Ning et al. (2015) also reported a similar observation of root length density. Chassot et al. (2001) and Qin et al. (2006) reported that root length density can reach 1.59 cm cm −3 in the Swiss Midlands. In Stuttgart, Germany, Wiesler and Horst (1994) observed the root growth and nitrate utilization of maize under field conditions. The observed root length density was 2.45-2.80 cm cm −3 at 0-30 cm depth, which was much higher than in other studies, and decreased to 0.01 cm cm −3 at 120-150 cm depth, which was consistent with the observations of Oikeh et al. (1999) at Samaru, Nigeria. Zhuang et al. (2001b) proposed a scaling model to estimate the distribution of the root length density of field-grown maize. In their study, the measured root length density in Tokyo, Japan, decreased from 0.4 to 0.95 cm cm −3 in the top soil layer to about 0.1 cm cm −3 in the bottom layer. Zhuang et al. (2001a) observed that the root length density of maize was mainly distributed at 0-60 cm depth, and the maximum values were about 0.9 cm cm −3 . These studies indicated that the root length density values were quite variable when this parameter was observed at different sites; never-theless, the simulated root length density in our study was of an order of magnitude that was similar to the observations from previous studies (Table 3). Figure 12 shows the modeled and observed 30 min canopy transpiration (T ), gross primary production (GPP), solarinduced fluorescence (SIF), and leaf water potential (LWP) from DOY 183 to 202 at Yangling station. The simulations by STEMMUS-SCOPE and SCOPE_SM were consistent with observations, whereas the simulated values from SCOPE were much higher than observations. The performance of STEMMUS-SCOPE and SCOPE_SM was consistent with that of SCOPE in the early morning and late afternoon, when photosynthesis was mainly limited by incident radiation rather than by water stress, intercellular CO 2 concentration, and V cmax . At midday, with increasing incident radiation, photosynthesis was mainly limited by water stress and V cmax , during which time the simulations by STEMMUS-SCOPE and SCOPE_SM were much better than that by SCOPE. The diurnal variation in the observed and modeled GPP was similar to that of T . Due to the lack of observed SIF, only the simulated SIF values were presented. As shown in Fig. 12, the SIF values simulated by STEMMUS-SCOPE and SCOPE_SM were reduced when the vegetation was experiencing water stress, which indicated that both the simulated SIF from STEMMUS-SCOPE and SCOPE_SM can respond to water stress. However, the accuracy of the simulated SIF requires further validation with field observations. Figure 13 shows the relationship among 30 min GPP, SIF, and LWP on DOY 199 at Yangling station. There was a strong linear relationship between SIF and GPP when the maize was well-watered (Fig. 13a). However, SIF kept increasing, whereas GPP tended to saturate when the maize was suffering water stress. This result is consistent with the previous study conducted for cotton and tobacco leaves (Van  der Tol et al., 2014). Because SCOPE_SM used the averaged root zone SM and ignored vertical root and soil water distribution, it overestimated GPP and SIF. When the maize was experiencing drought, the LWP was maintained at a low level. With GPP and T increasing, the plant decreased LWP in order to extract enough water from the root zone.

Diurnal variation in T, GPP, SIF, and LWP
The SPAC system enabled STEMMUS-SCOPE to simulate 30 min LWP. To better detect the response of simulated SIF to simulated LWP, we chose a cloudless day (DOY 199), and a liner relationship between the simulated SIF and LWP was obtained (Fig. 13b). Sun et al. (2016) reported that the SIFsoil moisture-drought relationship depended on variations in both absorbed PAR and fluorescence yield in response to water stress, whereas the LWP can reflect both the effect of absorbed PAR and the soil moisture status. The strong correlation between GPP, LWP, and SIF indicates the potential for using SIF as an effective signal for characterizing the response of photosynthesis to water stress. In the future, more studies should focus on the measurement of SIF, GPP, and LWP simultaneously for different vegetation types across different environmental conditions (radiation, soil moisture, and CO 2 concentration) to reveal how the water stress affects these relationships.

Limitations that need to be overcome
The new coupled model notably improved simulations of carbon and water fluxes when vegetation was suffering water stress. However, this study mainly aimed to improve the response of SCOPE to drought by introducing the vertical soil water and root profile. Some critical processes were followed that existed in SCOPE_SM and STEMMUS. As with any model, some modules in STEMMUS-SCOPE, such as plant hydraulics and root growth, could be improved upon in future development. First, to date many LSMs (e.g., CLM 5, Noah-MP, JULES, and CABLE) have incorporated a state-of-the-art plant hy-  De Kauwe et al., 2015). Although STEMMUS-SCOPE integrated a 1D root growth model and a relatively novel RWU model, its hydraulics model followed that in SCOPE_SM and ignored the most exciting recent advances in our understanding of plant hydraulics: hydraulic failure due to the loss of hydraulic conductivity owing to embolism and refilling for recovery from xylem embolism . Because STEMMUS-SCOPE performed well in maize cropland and grassland, the influence of embolism and refilling on water transfer from the soil through vegetation to the atmosphere cannot be fully detected. The value of using plant water potential instead of soil water potential to constrain model predictions has been demonstrated in many case studies (De Kauwe et al., 2020;Niu et al., 2020;Medlyn et al., 2016;Xu et al., 2016;Williams et al., 1996). Niu et al. (2020) followed the plant hydraulic model developed by Xu et al. (2016) and represented the plant stomatal water stress factor as a function of the plant water storage. CLM 5.0 also introduced a new formulation for WSF, which is based on leaf water potential (ψ leaf ) instead of soil water potential (ψ soil ; Kennedy et al., 2019). These new formulations based on plant water potential could offer significant improvements for plant drought responses. Furthermore, STEMMUS-SCOPE presently does not account for plant water storage; this may result in underestimating morning LE and overestimating afternoon LE. Some field observations have shown that the plant do not immediately respond when soil moisture is enhanced (Mackay et al., 2019), instead there are long lags, which were ignored in this study, between soil water recovery from drought and plant responses to the recovery. The WSF in STEMMUS-SCOPE directly comes from soil moisture and cannot reflect true stomatal response when vegetation is experiencing drought. For example, in early morning, the low stomatal aperture was induced by low PAR rather than by SM. Consequently, STEMMUS-SCOPE needs to introduce advanced hydraulics after the model has been tested in a wide range of ecosystems, particularly for vegetation exposed to frequent drought cycles or prolonged periods of severe drought events. It is important, however, to note that explicit representations of plant hydraulics require additional model parameters and increase the parameterization burden. This is the most challenging limitation to STEMMUS-SCOPE with respect to incorporating these hydraulics models, and we have chosen a trade-off between mechanism and practicality.
Second, as mentioned above, STEMMUS-SCOPE adapted the macroscopic RWU model and a simplified 1D root growth model in order to save on computational costs, although it predicted maximum root depth well, which is the most critical factor when calculating WSF and RWU. Such a simplification would likely ease the migration of our model into larger-scale models, such as Earth system models. However, STEMMUS-SCOPE oversimplified metabolic processes of the roots, including root exudates, root maintenance respiration, root growth respiration, and root turnover, which are also critical and have been incorporated in Noah-MP (Niu et al., 2020). This simplification could result in uncertainties in modeling the root growth and root water uptake. Meanwhile, there was no validation of the seasonal vertical root length distribution based on in situ observations, which need to be validated in the next step. Furthermore, the model presently does not account for the feedback between hydraulic controls over carbon allocation and the role of root growth on soil-plant hydraulics, which could also be considered in future model development.

Conclusions
A fundamental understanding of coupled energy, water, and carbon flux is vital for obtaining information on ecohydrological processes and functioning under climate change. The coupled model, STEMMUS-SCOPE, integrating radiative transfer, photochemistry, energy balance, root system dynamics, and soil moisture and soil temperature dynamics, has been proven to be a practical model to simulate detailed land surface processes such as evapotranspiration and GPP. In the coupled model, STEMMUS could provide the root zone moisture profile to SCOPE, which was used to calculate the water stress factor. On the other hand, SCOPE could provide the net carbon assimilation and soil surface temperature to STEMMUS, which was subsequently used as the top boundary condition and as the input for root growth model. This study explores the role of dynamic root growth in affecting canopy photosynthesis activities, fluorescence emissions, and evapotranspiration, which has not been reported before. The coupled model has been successfully applied in a maize field and a grassland and can be used to describe ET partitioning, canopy photosynthesis, reflectance, and fluorescence emissions. The results show that by considering dynamic root growth and the associated root water uptake, the simulated SIF of the coupled STEMMUS-SCOPE model can respond to water stress, whereas this is not the case for SCOPE_SM.
Through the intercomparison of SCOPE, SCOPE_SM, STEMMUS, and STEMMUS-SCOPE, we concluded that the coupled STEMMUS-SCOPE model can be used to investigate vegetation states under water-stressed conditions and to simultaneously understand the dynamics of soil heat and mass transfer, as well as the root growth. By considering the vertical distribution of soil moisture and the root system, the simulation of water and carbon fluxes, especially when vegetation was suffering moderate water stress, was significantly improved. However, the need remains for further studies to enhance the capacity of STEMMUS-SCOPE with respect to understanding ecosystem functioning. First of all, the estimation of the soil boundary condition, especially during the irrigation period, which has a significant influence on the simulation of soil temperature, requires further consideration. Second, the realism of the present model in modeling the water-stressed SIF will be subject to further studies. Nevertheless, STEMMUS-SCOPE may be used as an effective forward simulator to simulate remote sensing signals and to assimilate remote sensing data, such as solar-induced chlorophyll fluorescence, in order to improve the estimation of water and carbon fluxes. STEMMUS-SCOPE could also be used to investigate regional or global land surface processes, especially in arid and semiarid regions, due to its sensitivity to water-stressed conditions.

Appendix A A1 Photosynthesis and evapotranspiration under water stress in SCOPE
The C 4 photosynthesis is calculated in the SCOPE model as the minimum of three processes (Collatz et al., 1991(Collatz et al., , 1992: (1) the carboxylation rate limited by ribulose bisphosphatecarboxylase-oxygenase activity (known as rubisco (enzyme) limited), V c , described in Eq. (A1); (2) the carboxylation rate limited by ribulose 1-5 bisphosphate regeneration rate (known as RuBP (electron transport/light) limited), V e , described in Eq. (A2); and (3) at low CO 2 concentrations, the carboxylation rate limited by intercellular CO 2 partial pressure (p i ), V s , described in Eq. (A3).
The C 3 photosynthesis is calculated in the SCOPE model as the minimum of two processes (Farquhar et al., 1980): (1) the carboxylation rate limited by ribulose bisphosphatecarboxylase-oxygenase activity (known as rubisco (enzyme) limited), V c , described in Eq. (A5); and (2) the carboxylation rate limited by ribulose 1-5 bisphosphate regeneration rate (known as RuBP (electron transport/light) limited), V e , described in Eq. (A6).
Here, V cmax is the maximum carboxylation rate (µmol m −2 s −1 ), p i is the intercellular CO 2 partial pressure (Pa), k p is a pseudo-first-order rate constant for PEP carboxylase with respect to C i , P is the atmospheric pressure, A n is the net photosynthesis (µmol m −2 s −1 ), WSF is the total water stress factor, J is the electron transport rate (µmol m −2 s −1 ), C i is the intercellular CO 2 concentration (µmol m −3 ), C a is CO 2 concentration in the boundary layer (µmol m −3 ), m is Ball-Berry parameter, and RH is relative humidity at the leaf surface (%). In addition, leaf stomatal resistance r c (s m −1 ) is calculated as where ρ a is specific mass of air (kg m −3 ), M a is molecular mass of dry air (g mol −1 ), and p is atmosphere pressure (hPa). The calculation of latent heat flux (LE) is as follows: where λ is vaporization heat of water (J kg −1 ), q i is the humidity in stomata or soil pores (kg m −3 ), q a is the humidity above the canopy (kg m −3 ), r c is stomatal or soil surface resistance (s m −1 ), and r a is aerodynamic resistance (s m −1 ). In the study of Bayat et al. (2019), the water stress factor was calculated based on the root zone soil moisture content neglecting the distribution of root length. In this study, the water stress factor considered both root length distribution and water content in root zone. We use a sigmoid formulation rather than the piecewise function by Bayat et al. (2019). The calculations are as follows: where θ w is the soil water content at wilting point, θ f is the soil water content at field capacity, θ sat is the saturated soil water content, WSF(i) is the water stress factor at each soil layer, RF(i) is the ratio of root length in soil layer i (its calculation can be found in the Appendix A4), and SM(i) is the soil moisture at each soil layer.

A2.1 Soil water conservation equation
The soil water conservation equation is as follows: where ρ L and ρ V (kg m −3 ) are the density of liquid water and water vapor, respectively; q L and q V (m 3 m −3 ) are the volumetric water content for liquid and water vapor, respectively; z(m) is the vertical space coordinate (positive upwards); S (cm s −1 ) is the sink term for the root water extraction; K (m s −1 ) is hydraulic conductivity; h (cm) is the pressure head; T s ( • C) is the soil temperature; P g (Pa) is the mixed pore-air pressure; γ w (kg m −2 s −2 ) is the specific weight of water; D TD (kg m −1 s −1 ( • C) −1 ) is the transport coefficient for adsorbed liquid flow due to the temperature gradient; D Vh (kg m −2 s −1 ) is the isothermal vapor conductivity; D VT (kg m −1 s −1 ( • C) −1 ) is the thermal vapor diffusion coefficient; D Va is the advective vapor transfer coefficient (Zeng et al., 2011a, b); q Lh , q LT , and q La (kg m −2 s −1 ) are the liquid water fluxes driven by the gradient of matric potential, temperature, and air pressure, respectively; and q Vh , q VT , and q Va (kg m −2 s −1 ) are the water vapor fluxes driven by the gradient of matric potential, temperature, and air pressure, respectively.

A2.2 Dry-air conservation equation
The dry-air conservation equation is as follows: where ε is the porosity, ρ da (kg m −3 ) is the density of dry air, S a (= 1 − S L ) is the degree of air saturation in the soil, S L (= θ L /ε) is the degree of saturation in the soil, H c is Henry's constant, D e (m 2 s −1 ) is the molecular diffusivity of water vapor in soil, K g (m 2 ) is the intrinsic air permeability, m a (kg m −2 s −1 ) is the air viscosity, q L (kg m −2 s −1 ) is the liquid water flux, θ a (= θ V ) is the volumetric fraction of dry air in the soil, and D Vg (m 2 s −1 ) is the gas-phase longitudinal dispersion coefficient (Zeng et al., 2011a,b).

A2.3 Energy balance equation
The energy balance equation is as follows: where C s , C L , C V , and C a (J kg −1 ( • C) −1 ) are the specific heat capacities of solids, liquid, water vapor, and dry air, respectively; ρ s (kg m −3 ), ρ L (kg m −3 ), ρ V (kg m −3 ), and ρ da (kg m −3 ) are the density of solids, liquid water, water vapor, and dry air, respectively; θ s is the volumetric fraction of solids in the soil; θ L , θ V , and θ a are the volumetric fraction of liquid water, water vapor, and dry air, respectively; T r ( • C) is the reference temperature; L 0 (J kg −1 ) is the latent heat of vaporization of water at temperature T r ; W (J kg −1 ) is the differential heat of wetting (the amount of heat released when a small amount of free water is added to the soil matrix); λ eff (W m −1 ( • C) −1 ) is the effective thermal conductivity of the soil; and q L , q V , and q a (kg m −2 s −1 ) are the liquid, vapor water and dry air flux, respectively.
A3 Dynamic root growth modeling

A3.1 Root front growth
The depth of the root front is firstly initialized either with the sowing depth for sown crops or with an initial value for transplanted crops or perennial crops. The root front growth stops when it reaches a certain depth of soil or a physical or chemical obstacle preventing root growth, but it also stops when the phenological stopping stage has been reached.
where Z is root front growth at the tth time step, D Z (cm) is the root zone depth, T air ( • C) is air temperature, T min ( • C) is the minimum temperature for root growth, T max ( • C) is the maximum temperature for root growth, and RGR (cm ( • C) −1 d −1 ) is the root growth rate of root front.

A3.2 Root length growth
In this study, the root distribution in the root zone was realized via simulating the root length growth in each soil layer.
where fr root is the allocation fraction of the net assimilation to root, fr root is assumed as a function of leaf area index (LAI) and root zone water content, A n is the net assimilation rate (µmol m −2 s −1 ), R C is the ratio of carbon to dry organic matter in root, R D is root density (g m −3 ), r root is the radius of the root, and Rl_tot (m m −3 ) is the total root length growth. The limiting factors for allocation are preliminarily computed, and they account for root zone soil moisture availability, A W , and light availability, A L .
where WSF is the averaged soil moisture stress factor in the root zone.
where K e = 0.15 is a constant light extinction coefficient.
where r min (= 0.15) is the minimum allocation coefficient to fine roots, and r 0 is a coefficient that indicates the theoretically unstressed allocation to fine roots.
where RF(i) is the allocation fraction of root growth length in layer i, and Rl(i) is the root growth length in layer i. For i = 1 to n − 1 (i = 1 refers to the top soil layer), For i = n, Here, Rl t i and Rl t−1 i are the root length of layer i at time step t and time step t − 1, respectively.
where Rl T is the total root length in the root zone, and Rl(i) is the root length in soil layer i. At the root front, the density is imposed and estimated by the parameter L v_front , and the growth in root length depends directly on the root front growth rate Z: (A26)

A4 Root water uptake
The equation to calculate root water uptake and transpiration is as follows: where ψ s,i is the soil water potential of layer i (pressure head, unit: m), ψ l is leaf water potential (m), r s,i is the soil hydraulic resistance (s m −1 ), r r,i is the root resistance to water flow radially across the roots (s m −1 ), and r x,i is the plant axial resistance to flow from the soil to the leaves (s m −1 ). e l and e a are vapor pressure of leaf and the atmosphere (hPa), respectively, and r a and r c are aerodynamic resistance and canopy resistance (s m −1 ), respectively. ρ da is the density of dry air (kg m −3 ). ρ V is the density of water vapor. P is the atmospheric pressure (Pa). The ratio of the molar mass of water to air is 0.622. ψ s,i is described as a function of soil moisture by Van Genuchten (1980), and the relevant parameters are shown in Table B1.
The r s is calculated by Reid and Huck (1990) as follows: where B is the root length activity factor, K is hydraulic conductivity of soil (m s −1 ), L v is root length density (m m −3 ), and d is the thickness of the soil layer (m). B is calculated as where r root is root radius (m). The r r is estimated as (Reid and Huck, 1990) follows: where P r is root radial resistivity (s m −1 ). The xylem resistance r x is estimated by Klepper et al. (1983): where P a is root axial resistivity (s m −3 ), Z mid is the depth of the midpoint of the soil layer, and f is a fraction defined for a specific depth as the number of roots that connect directly to the stem base to total roots crossing a horizontal plane at that depth. We can consider it equal to 0.22 based on Klepper et al. (1983). The updated root water uptake term is In contrast to other studies that need to calculate the compensatory water uptake and hydraulic redistribution after calculating the standard water uptake of each soil layer, the sink term in this study is calculated by a physically based model that contains the effect of root resistance and soil hydraulic resistance rather than only considering the root fraction; thus, the compensatory water uptake and hydraulic redistribution have been considered when calculating the sink term.
Appendix B   Mixed pore-air pressure Pa q L Liquid water flux kg m −2 s −1 q Lh Liquid water flux driven by the gradient of matric potential kg m −2 s −1 q LT Liquid water flux driven by the gradient of temperature kg m −2 s −1 q La Liquid water flux driven by the gradient of air pressure kg m −2 s −1 q V Water vapor flux kg m −2 s −1 q Vh Water vapor flux driven by the gradient of matric potential kg m −2 s −1 q VT Water vapor flux driven by the gradient of temperature kg m −2 s −1 q Va Water vapor flux driven by the gradient of air pressure kg m −2 s −1 q a Dry-air flux kg m −2 s −1 S Sink term for the root water extraction cm s −1 S a Degree of air saturation in the soil S L Degree of saturation in the soil SM(i) The soil moisture at a specific soil layer m 3 m −3 T s Soil temperature • C T s0 Soil surface temperature • C T r Reference temperature • C 20 20 W Differential heat of wetting J kg −1 1.001 × 10 3 1.001 × 10 3 WSF Total water stress factor WSF(i) Water stress factor at a specific soil layer Code and data availability. The development and validation of STEMMUS-SCOPE in this paper were conducted in MAT-LAB R2016a. The exact version of the model used to produce the results employed in this paper is archived on Zenodo (https://doi.org/10.5281/zenodo.3839092, Wang et al., 2020). The original source of the SCOPE model and STEMMUS model was obtained from Van der Tol et al. (2009) and Zeng et al. (2011a, b), respectively. The tower-based eddy-covariance measurements used for model validation were provided by the authors for the Yangling station, China (https://doi.org/10.4121/uuid:aa0ed483-701e-4ba0-b7b0-674695f5f7a7, Wang et al., 2019), and were obtained from the FLUXNET2015 Dataset and PLUMBER2 program for the Vaira Ranch (US-Var) FLUXNET site.
Author contributions. YW, YZ, HC, and ZS designed the study. YW developed the code, conducted the analysis, and wrote the paper. YW and HC collected and shared their eddy-covariance measurements for the purpose of model validation. All authors discussed, commented, and contributed to the revisions and final version of the article.
Competing interests. The authors declare that they have no conflict of interest.