Soil Vegetation Atmosphere Transfer Research Articles

Extended HYDRUS-1D freezing module emphasizes thermal conductivity schemes for simulation of soil hydrothermal dynamics

Soil thermal conductivity (λ) is required to investigate coupled heat and water transport in disciplines such as agriculture, hydrology and engineering. Parameterization schemes or models of λ are also the critical input parameter for various numerical simulation programs like the widely used HYDRUS, one of the most commonly used models for mimicking water, heat, and solute transport. However, λ has not received enough attention in HYDRUS, and it remains unclear how different λ schemes affect the simulated soil water and thermal regimes. Thus, we programmed 24 λ schemes (including two built-in schemes) used in mainstream land surface, hydrological, and soil–vegetation–atmosphere transfer models into HYDRUS-1D (freezing module) to assess the effects of different λ schemes on soil temperature and water content simulations under freezing–thawing. The results showed that the 24 λ schemes performed differently in the simulation of soil temperature within 1 m below ground, with eight λ schemes, i.e., de Vries 1963 scheme/DV1963 (R=0.99), Camillo and Schmugge 1981/CS1981 scheme (R=0.98), Desborough and Pitman 1998/DP1998 scheme (R=0.97), Cass et al. 1984/CS1984 scheme (R=0.96), Shmakin 1998/SA1998 scheme (R=0.95), Dharssi et al. 2009/DI2009 scheme (R=0.95), Becker et al. 1992/BB1992 scheme (R=0.95), Hubrechts 1998/HL1998 scheme (R=0.94), performing superiorly to the built-in λ schemes. This study highlights the importance of choosing appropriate λ schemes in soil water and heat simulations.

Impact of vegetation composition and seasonality on sensitivity of modelled CO2 exchange in temperate raised bogs

Encroachment of vascular plants (VP) in temperate raised bogs, as a consequence of altered hydrological conditions and nutrient input, is widely observed. Effects of such vegetation shift on water and carbon cycles are, however, largely unknown and identification of responsible plant physiological traits is challenging. Process-based modelling offers the opportunity of gaining insights into ecosystem functioning beyond observations, and to infer decisive trait shifts of plant functional groups. We adapted the Soil–Vegetation–Atmosphere Transfer model pyAPES to a temperate raised bog site by calibration against measured peat temperature, water table and surface CO2 fluxes. We identified the most important traits determining CO2 fluxes by conducting Morris sensitivity analysis (MSA) under changing conditions throughout the year and simulated VP encroachment. We further investigated transferability of results to other sites by extending MSA to parameter ranges derived from literature review. We found highly variable intra-annual plant traits importance determining ecosystem CO2 fluxes, but only a partial shift of importance of photosynthetic processes from moss to VP during encroachment. Ecosystem respiration was dominated by peat respiration. Overall, carboxylation rate, base respiration rate and temperature sensitivity (Q10) were most important for determining bog CO2 balance and parameter ranking was robust even under the extended MSA.

Recent Developments to the SimSphere Land Surface Modelling Tool for the Study of Land–Atmosphere Interactions

Soil–Vegetation–Atmosphere Transfer (SVAT) models are a promising avenue towards gaining a better insight into land surface interactions and Earth’s system dynamics. One such model developed for the academic and research community is the SimSphere SVAT model, a popular software toolkit employed for simulating interactions among the layers of vegetation, soil, and atmosphere on the land surface. The aim of the present review is two-fold: (1) to deliver a critical assessment of the model’s usage by the scientific and wider community over the last 15 years, and (2) to provide information on current software developments implemented in the model. From the review conducted herein, it is clearly evident that from the models’ inception to current day, SimSphere has received notable interest worldwide, and the dissemination of the model has continuously grown over the years. SimSphere has been used so far in several applications to study land surface interactions. The validation of the model performed worldwide has shown that it is able to produce realistic estimates of land surface parameters that have been validated, whereas detailed sensitivity analysis experiments conducted with the model have further confirmed its structure and architectural coherence. Furthermore, the recent inclusion of novel functionalities in the model, as outlined in the present review, has clearly resulted in improving its capabilities and in opening up new opportunities for its use by the wider community. SimSphere developments are also ongoing in different aspects, and its use as a toolkit towards advancing our understanding of land surface interactions from both educational and research points of view is anticipated to grow in the coming years.

Implications of seasonal changes in photosynthetic traits and leaf area index for canopy CO2 and H2O fluxes in a Japanese cedar (Cryptomeria japonica D. Don) plantation

Recently, requests for forest management to conserve water resources and sequester carbon have increased in Japan, where approximately 27% and 12% of the total land area are covered by forest plantations and Japanese cedar (Cryptomeria japonica D. Don). To formulate better forest management practices that maximize forest functions by maintaining water and carbon cycles, an understanding of how leaf-level ecophysiological traits influence canopy-level CO2 and H2O exchanges in Japanese cedar plantations is required. In this study, observations of eddy covariance flux and leaf ecophysiology in a Japanese cedar plantation in southern Japan indicate that significantly greater productivity and seasonal variations tend to occur in canopy CO2 and H2O fluxes, with seasonal patterns of photosynthetic traits (maximum carboxylation rate normalized at 25 °C, Vcmax,25) and leaf area index (LAI). We also examined how biotic (e.g., Vcmax,25 and LAI) and abiotic (e.g., meteorological variables) factors govern canopy fluxes using a multi-layer soil–vegetation–atmosphere transfer (SVAT) model. The model validation suggests that seasonal variations in Vcmax,25 and LAI must be included to reproduce the fluxes observed in the plantation. Quantitative experiments using the validated SVAT model indicate that considering cold acclimation during the winter improved the reproducibility of the model with regard to the measured fluxes, constant high intra-annual LAI had the least impact on forest productivity, and a consistently high Vcmax,25 value increased forest productivity but considerably decreased the nitrogen use efficiency and required a larger nitrogen supply from the ecosystem during the winter. These findings highlight the necessity of decreasing photosynthetic ability during the winter, which safeguards the ecosystem's nitrogen resources to sustain the productivity of the Japanese cedar plantation throughout the year.

Does interception evaporation occur from both sides of leaves at the wet Japanese cypress canopy during and after rainfall?

AbstractThe interception by the lower (abaxial) side of needle leaves not only contributes to the forest evapotranspiration but will also cause a diffusion barrier of CO2 over the stomata and depresses the plant gas exchange. This study combined the eddy covariance (EC) technique and leaf wetness measurement with a soil–vegetation–atmosphere transfer (SVAT) multilayer model to examine the occurrence of interception evaporation from both sides of leaves at a Japanese cypress forest canopy relating to rainfall intensity and different wetting periods. We compared the measured latent heat flux (λE) with the simulated wet canopy λE with two models that interception evaporation only happens from the upper (adaxial) side of leaves and both sides of leaves. Both models showed a low λE during the rainfall as the EC data did. The simulated λE at the wet period after rainfall indicates that the interception evaporation from both sides of leaves is more likely to happen after heavy rainfall (>15 mm/12 h). However, for the most frequent small rainfall events (0–5 mm/12 h) at this site, interception evaporation is more likely to occur only from the adaxial side than both sides, which helps the wet leaves to maintain stomata opening and process CO2 uptake after rainfall.

Understanding root-zone soil moisture in agricultural regions of Central Mexico using the ensemble Kalman filter, satellite-derived information, and the THEXMEX-18 dataset

ABSTRACT An Ensemble Kalman Filter (EnKF)-based assimilation algorithm was implemented to estimate root-zone soil moisture (RZSM) using a Soil-Vegetation-Atmosphere Transfer (SVAT) model during a complete growing season of corn in Central Mexico. Synthetic and field soil moisture (SM) observations and NASA SMAP SM retrievals were used to understand the effect of vertically spatial updates and uncertainties in meteorological forcings on RZSM estimates. Assimilation of RZSM every 3 days using SM observations at 4 depths lowered the averaged standard deviation (ASD) and the root mean square error (RMSE) by 60 % and 50 %, respectively, compared to the open-loop ASD. The assimilation of synthetic SM at the top 0-5 cm obtained RZSM closer to observations compared to THEXMEX-18 SM measurements and SMAP SM retrievals. Differences between EnKF estimates and SM observations and SMAP SM retrievals are mainly due to misrepresentation of vegetation conditions. The results improved SM estimates up to 10-cm depth using SMAP SM retrievals; however, additional studies are needed to improve SM at deeper layers. The implemented methodology can estimate SM at the top 10 cm of the soil every 3 days to mitigate the impact of the climate change on agricultural production over rainfed areas, particularly in developing countries.

Simsphere: A Downloadable Soil–Vegetation–Atmosphere–Transfer (SVAT) Model for Teaching and Research

Abstract Simsphere, a soil–vegetation–atmosphere–transfer (SVAT) model developed at The Pennsylvania State University, can be downloaded from the web for use by students and researchers. In existence for several decades, Simsphere has figured in both the classroom and in research at several universities. As such, Simsphere has been supported by a knowledgeable group of academic users and has been applied in a variety of applications, such as in remote sensing of surface soil water content and in the assessment of water and ozone stresses on plants. This paper describes the model and how it can be downloaded and run.

Hydrological Functioning of Maize Crops in Southwest France Using Eddy Covariance Measurements and a Land Surface Model

The primary objective of this study is to evaluate the representation of the energy budget for irrigated maize crops in soil–vegetation–atmosphere transfer (SVAT) models. To this end, a comparison between the original version of the interactions between the soil–biosphere–atmosphere (ISBA) model based on a single-surface energy balance and the new ISBA-multi-energy balance (ISBA-MEB) option was carried out. The second objective is to analyze the intra- and inter-seasonal variability of the crop water budget by implementing ISBA and ISBA-MEB over six irrigated maize seasons between 2008 and 2019 in Lamasquère, southwest France. Seasonal dynamics of the convective fluxes were properly reproduced by both models with R2 ranging between 0.66 and 0.80 (RMSE less than 59 W m−2) for the sensible heat flux and between 0.77 and 0.88 (RMSE less than 59 W m−2) for the latent heat flux. Statistical metrics also showed that over the six crop seasons, for the turbulent fluxes, ISBA-MEB was consistently in better agreement with the in situ measurements with RMSE 8–30% lower than ISBA, particularly when the canopy was heterogeneous. The ability of both models to partition the evapotranspiration (ET) term between soil evaporation and plant transpiration was also acceptable as transpiration predictions compared very well with the available sap flow measurements during the summer of 2015; (ISBA-MEB had slightly better statistics than ISBA with R2 of 0.91 and a RMSE value of 0.07 mm h−1). Finally, the results from the analysis of the inter-annual variability of the crop water budget can be summarized as follows: (1) The partitioning of the ET revealed a strong year-to-year variability with transpiration ranging between 40% and 67% of total ET, while soil evaporation was dominant in 2008 and 2010 due to the late and poor canopy development; (2) drainage losses are close to null because of an impervious layer at 60 cm depth; and (3) this very specific condition limited the inter-annual variability of irrigation scheduling as crops can always extract water that is stored in the root zone.

Evapotranspiration partition using the multiple energy balance version of the ISBA-A-g<sub>s</sub> land surface model over two irrigated crops in a semi-arid Mediterranean region (Marrakech, Morocco)

Abstract. The main objective of this work is to question the representation of the energy budget in soil–vegetation–atmosphere transfer (SVAT) models for the prediction of the turbulent fluxes in the case of irrigated crops with a complex structure (row) and under strong transient hydric regimes due to irrigation. To this end, the Interaction between Soil, Biosphere, and Atmosphere (ISBA-A-gs) is evaluated at a complex open olive orchard and, for the purposes of comparison, on a winter wheat field taken as an example of a homogeneous canopy. The initial version of ISBA-A-gs, based on a composite energy budget (hereafter ISBA-1P for one patch), is compared to the new multiple energy balance (MEB) version of ISBA that represents a double source arising from the vegetation located above the soil layer. In addition, a patch representation corresponding to two adjacent, uncoupled source schemes (hereafter ISBA-2P for two patches) is also considered for the olive orchard. Continuous observations of evapotranspiration (ET), with an eddy covariance system and plant transpiration (Tr) with sap flow and isotopic methods were used to evaluate the three representations. A preliminary sensitivity analyses showed a strong sensitivity to the parameters related to turbulence in the canopy introduced in the new ISBA–MEB version. For wheat, the ability of the single- and dual-source configuration to reproduce the composite soil–vegetation heat fluxes was very similar; the root mean square error (RMSE) differences between ISBA-1P, ISBA-2P and ISBA–MEB did not exceed 10 W m−2 for the latent heat flux. These results showed that a composite energy balance in homogeneous covers is sufficient to reproduce the total convective fluxes. The two configurations are also fairly close to the isotopic observations of transpiration in spite of a light underestimation (overestimation) of ISBA-1P (ISBA–MEB). At the olive orchard, contrasting results are obtained. The dual-source configurations, including both the uncoupled (ISBA-2P) and the coupled (ISBA–MEB) representations, outperformed the single-source version (ISBA-1P), with slightly better results for ISBA–MEB in predicting both total heat fluxes and evapotranspiration partition. Concerning plant transpiration in particular, the coupled approach ISBA–MEB provides better results than ISBA-1P and, to a lesser extent, ISBA-2P with RMSEs of 1.60, 0.90, and 0.70 mm d−1 and R2 of 0.43, 0.69, and 0.70 for ISBA-1P, ISBA-2P and ISBA–MEB, respectively. In addition, it is shown that the acceptable predictions of composite convective fluxes by ISBA-2P for the olive orchard are obtained for the wrong reasons as neither of the two patches is in agreement with the observations because of a bad spatial distribution of the roots and a lack of incoming radiation screening for the bare soil patch. This work shows that composite convection fluxes predicted by the SURFace EXternalisée (SURFEX) platform and the partition of evapotranspiration in a highly transient regime due to irrigation is improved for moderately open tree canopies by the new coupled dual-source ISBA–MEB model. It also points out the need for further local-scale evaluations on different crops of various geometry (more open rainfed agriculture or a denser, intensive olive orchard) to provide adequate parameterisation to global database, such as ECOCLIMAP-II, in the view of a global application of the ISBA–MEB model.

Temporal interpolation of land surface fluxes derived from remote sensing – results with an unmanned aerial system

Abstract. Remote sensing imagery can provide snapshots of rapidly changing land surface variables, e.g. evapotranspiration (ET), land surface temperature (Ts), net radiation (Rn), soil moisture (θ), and gross primary productivity (GPP), for the time of sensor overpass. However, discontinuous data acquisitions limit the applicability of remote sensing for water resources and ecosystem management. Methods to interpolate between remote sensing snapshot data and to upscale them from an instantaneous to a daily timescale are needed. We developed a dynamic soil–vegetation–atmosphere transfer model to interpolate land surface state variables that change rapidly between remote sensing observations. The “Soil–Vegetation, Energy, water, and CO2 traNsfer” (SVEN) model, which combines the snapshot version of the remote sensing Priestley–Taylor Jet Propulsion Laboratory ET model and light use efficiency GPP models, now incorporates a dynamic component for the ground heat flux based on the “force-restore” method and a water balance “bucket” model to estimate θ and canopy wetness at a half-hourly time step. A case study was conducted to demonstrate the method using optical and thermal data from an unmanned aerial system at a willow plantation flux site (Risoe, Denmark). Based on model parameter calibration with the snapshots of land surface variables at the time of flight, SVEN interpolated UAS-based snapshots to continuous records of Ts, Rn, θ, ET, and GPP for the 2016 growing season with forcing from continuous climatic data and the normalized difference vegetation index (NDVI). Validation with eddy covariance and other in situ observations indicates that SVEN can estimate daily land surface fluxes between remote sensing acquisitions with normalized root mean square deviations of the simulated daily Ts, Rn, θ, LE, and GPP of 11.77 %, 6.65 %, 19.53 %, 14.77 %, and 12.97 % respectively. In this deciduous tree plantation, this study demonstrates that temporally sparse optical and thermal remote sensing observations can be used to calibrate soil and vegetation parameters of a simple land surface modelling scheme to estimate “low-persistence” or rapidly changing land surface variables with the use of few forcing variables. This approach can also be applied with remotely-sensed data from other platforms to fill temporal gaps, e.g. cloud-induced data gaps in satellite observations.

Using Statistical and Dynamical Downscaling to Assess Climate Change Impacts on Mine Reclamation Cover Water Balances

The oil sands industry in Canada uses soil–vegetation–atmosphere-transfer (SVAT) water balance models, calibrated against short-term (<≈ 10 years) field monitoring data, to evaluate long-term (≈60 years) reclamation cover design performance. These evaluations use long-term historical climate data; however, the effects of climate change should also be incorporated in these analyses. Although statistical downscaling of global climate change projections is commonly used to obtain local, site-specific climate, high resolution dynamical downscaling can also be used. The value of this latter approach to obtain local site-specific projections for mine reclamation covers has not been evaluated previously. This study explored the differences in key water balance components of three reclamation covers and three natural sites in northern Alberta, Canada, under future, site-specific, statistical, and dynamical climate change projections. Historical meteorological records were used to establish baseline periods. Temperature datasets were used to calculate potential evapotranspiration (PET) using the Hargreaves–Samani method. Statistical downscaling uses the Long Ashton Research Station Weather Generator (LARS-WG) and global circulation model (GCM) projections of temperature and precipitation. Dynamical climate change projections were generated on a 4 km grid using the weather research and forecasting (WRF) model. These climate projections were applied to a physically-based water balance model (i.e. Hydrus-1D) to simulate actual evapotranspiration (AET) and net percolation (NP) for the baseline and future periods. The key findings were: (a) LARS-WG outperformed WRF in simulating baseline temperatures and precipitation; (b) both downscaling methods showed similar directional shifts in the future temperatures and precipitation; (c) this, in turn, created similar directional shifts in future growing season median AET and NP, although the increase in future NP for LARS-WG was higher than that for WRF. The relative increases in future NP were much higher than the relative increases in future AET, particularly for the reclamation covers.

Limitations in Validating Derived Soil Water Content from Thermal/Optical Measurements Using the Simplified Triangle Method

We assess the validity of the surface moisture availability parameter (Mo) derived from satellite-based optical/thermal measurements using the simplified triangle method. First, we show that Mo values obtained from the simplified triangle method agree closely with those generated from a soil/vegetation/atmosphere/transfer (SVAT) model for scenes over a field site at the Allahabad district, India. Next, we compared Mo values from the simplified triangle method for these same overpass scenes with surface soil water content measured at depths of 5 and 15 cm at this field site. Although a very weak correlation exists between remotely sensed values of Mo for the full scenes and measured soil water content measured at both depths, correlations increasingly improve for the 5 cm samples (but not for the 15 cm samples) as pixels were limited to increasingly smaller vegetation fractions. We conclude that the simplified triangle method would yield reasonable values of Mo and demonstrate good agreement with ground measurements, provided that validation is limited to pixels with little or no vegetation and to soil depths of 5 cm or less.

Characterizing uncertainty in the hydraulic parameters of oil sands mine reclamation covers and its influence on water balance predictions

Abstract. One technique to evaluate the performance of oil sands reclamation covers is through the simulation of long-term water balance using calibrated soil–vegetation–atmosphere transfer models. Conventional practice has been to derive a single set of optimized hydraulic parameters through inverse modelling (IM) based on short-term (&lt;5–10 years) monitoring datasets. This approach is unable to characterize the impact of variability in the cover properties. This study utilizes IM to optimize the hydraulic properties for 12 soil cover designs, replicated in triplicate, at Syncrude's Aurora North mine site. The hydraulic parameters for three soil types (peat cover soil, coarse-textured subsoil, and lean oil sand substrate) were optimized at each monitoring site from 2013 to 2016. The resulting 155 optimized parameter values were used to define distributions for each parameter/soil type, while the progressive Latin hypercube sampling (PLHS) method was used to sample parameter values randomly from the optimized parameter distributions. Water balance models with the sampled parameter sets were used to evaluate variations in the maximum sustainable leaf area index (LAI) for five illustrative covers and quantify uncertainty associated with long-term water balance components and LAI values. Overall, the PLHS method was able to better capture broader variability in the water balance components than a discrete interval sampling method. The results also highlight that climate variability dominates the simulated variability in actual evapotranspiration and that climate and parameter uncertainty have a similar influence on the variability in net percolation.

Data assimilation of high-resolution thermal and radar remote sensing retrievals for soil moisture monitoring in a drip-irrigated vineyard

Efficient water use assessment and irrigation management is critical for the sustainability of irrigated agriculture, especially under changing climate conditions. Due to the impracticality of maintaining ground instrumentation over wide geographic areas, remote sensing and numerical model-based fine-scale mapping of soil water conditions have been applied for water resource applications at a range of spatial scales. Here, we present a prototype framework for integrating high-resolution thermal infrared (TIR) and synthetic aperture radar (SAR) remote sensing data into a soil-vegetation-atmosphere-transfer (SVAT) model with the aim of providing improved estimates of surface- and root-zone soil moisture that can support optimized irrigation management strategies. Specifically, remotely-sensed estimates of water stress (from TIR) and surface soil moisture retrievals (from SAR) are assimilated into a 30-m resolution SVAT model over a vineyard site in the Central Valley of California, U.S. The efficacy of our data assimilation algorithm is investigated via both the synthetic and real data experiments. Results demonstrate that a particle filtering approach is superior to an ensemble Kalman filter for handling the nonlinear relationship between model states and observations. In addition, biophysical conditions such as leaf area index are shown to impact the relationship between observations and states and must therefore be represented accurately in the assimilation model. Overall, both surface and root-zone soil moisture predicted via the SVAT model are enhanced through the assimilation of thermal and radar-based retrievals, suggesting the potential for improving irrigation management at the agricultural sub-field-scale using a data assimilation strategy.

Ability of a soil–vegetation–atmosphere transfer model and a two-source energy balance model to predict evapotranspiration for several crops and climate conditions

Abstract. The heterogeneity of Agroecosystems, in terms of hydric conditions, crop types and states, and meteorological forcing, is difficult to characterize precisely at the field scale over an agricultural landscape. This study aims to perform a sensitivity study with respect to the uncertain model inputs of two classical approaches used to map the evapotranspiration of agroecosystems: (1) a surface energy balance (SEB) model, the Two-Source Energy Balance (TSEB) model, forced with thermal infrared (TIR) data as a proxy for the crop hydric conditions, and (2) a soil–vegetation–atmosphere transfer (SVAT) model, the SEtHyS model, where hydric conditions are computed from a soil water budget. To this end, the models' skill was compared using a large and unique in situ database covering different crops and climate conditions, which was acquired over three experimental sites in southern France and Morocco. On average, the models provide 30 min estimations of latent heat flux (LE) with a RMSE of around 55 W m−2 for TSEB and 47 W m−2 for SEtHyS, and estimations of sensible heat flux (H) with a RMSE of around 29 W m−2 for TSEB and 38 W m−2 for SEtHyS. A sensitivity analysis based on realistic errors aimed to estimate the potential decrease in performance induced by the spatialization process. For the SVAT model, the multi-objective calibration iterative procedure (MCIP) is used to determine and test different sets of parameters. TSEB is run with only one set of parameters and provides acceptable performance for all crop stages apart from the early growing season (LAI &lt; 0.2 m2 m−2) and when hydric stress occurs. An in-depth study on the Priestley–Taylor key parameter highlights its marked diurnal cycle and the need to adjust its value to improve flux partitioning between the sensible and latent heat fluxes (1.5 and 1.25 for France and Morocco, respectively). Optimal values of 1.8–2 were highlighted under cloudy conditions, which is of particular interest due to the emergence of low-altitude drone acquisition. Under developed vegetation (LAI &gt; 0.8 m2 m−2) and unstressed conditions, using sets of parameters that only differentiate crop types is a valuable trade-off for SEtHyS. This study provides some scientific elements regarding the joint use of both approaches and TIR imagery, via the development of new data assimilation and calibration strategies.

Multiple-constraint inversion of SCOPE. Evaluating the potential of GPP and SIF for the retrieval of plant functional traits

The most recent efforts to provide remote sensing (RS) estimates of plant function rely on the combination of Radiative Transfer Models (RTM) and Soil-Vegetation-Atmosphere Transfer (SVAT) models, such as the Soil-Canopy Observation Photosynthesis and Energy fluxes (SCOPE) model. In this work we used ground spectroradiometric and chamber-based CO2 flux measurements in a nutrient manipulated Mediterranean grassland in order to: 1) develop a multiple-constraint inversion approach of SCOPE able to retrieve vegetation biochemical, structural as well as key functional traits, such as chlorophyll concentration (Cab), leaf area index (LAI), maximum carboxylation rate (Vcmax) and the Ball-Berry sensitivity parameter (m); and 2) compare the potential of the of gross primary production (GPP) and sun-induced fluorescence (SIF), together with up-welling Thermal Infrared (TIR) radiance and optical reflectance factors (RF), to estimate such parameters. The performance of the proposed inversion method as well as of the different sets of constraints was assessed with contemporary measurements of water and heat fluxes and leaf nitrogen content, using pattern-oriented model evaluation.The multiple-constraint inversion approach proposed together with the combination of optical RF and diel GPP and TIR data provided reliable estimates of parameters, and improved predicted water and heat fluxes. The addition of SIF to this scheme slightly improved the estimation of m. Parameter estimates were coherent with the variability imposed by the fertilization and the seasonality of the grassland. Results revealed that fertilization had an impact on Vcmax, while no significant differences were found for m. The combination of RF, SIF and diel TIR data weakly constrained functional traits. Approaches not including GPP failed to estimate LAI; however GPP overestimated Cab in the dry period. These problems might be related to the presence of high fractions of senescent leaves in the grassland. The proposed inversion approach together with pattern-oriented model evaluation open new perspectives for the retrieval of plant functional traits relevant for land surface models, and can be utilized at various research sites where hyperspectral remote sensing imagery and eddy covariance flux measurements are simultaneously taken.

An improved perspective in the spatial representation of soil moisture: potential added value of SMOS disaggregated 1 km resolution “all weather” product

Abstract. This study uses the synergy of multi-resolution soil moisture (SM) satellite estimates from the Soil Moisture Ocean Salinity (SMOS) mission, a dense network of ground-based SM measurements, and a soil–vegetation–atmosphere transfer (SVAT) model, SURFEX (externalized surface), module ISBA (interactions between soil, biosphere and atmosphere), to examine the benefits of the SMOS level 4 (SMOS-L4) version 3.0, or “all weather” high-resolution soil moisture disaggregated product (SMOS-L43.0; ∼1 km). The added value compared to SMOS level 3 (SMOS-L3; ∼25 km) and SMOS level 2 (SMOS-L2; ∼15 km) is investigated. In situ SM observations over the Valencia anchor station (VAS; SMOS calibration and validation – Cal/Val – site in Europe) are used for comparison. The SURFEX (ISBA) model is used to simulate point-scale surface SM (SSM) and, in combination with high-quality atmospheric information data, namely from the European Centre for Medium-Range Weather Forecasts (ECMWF) and the Système d'analyse fournissant des renseignements atmosphériques à la neige (SAFRAN) meteorological analysis system, to obtain a representative SSM mapping over the VAS. The sensitivity to realistic initialization with SMOS-L43.0 is assessed to simulate the spatial and temporal distribution of SSM. Results demonstrate the following: (a) All SMOS products correctly capture the temporal patterns, but the spatial patterns are not accurately reproduced by the coarser resolutions, probably in relation to the contrast with point-scale in situ measurements. (b) The potential of the SMOS-L43.0 product is pointed out to adequately characterize SM spatio-temporal variability, reflecting patterns consistent with intensive point-scale SSM samples on a daily timescale. The restricted temporal availability of this product dictated by the revisit period of the SMOS satellite compromises the averaged SSM representation for longer periods than a day. (c) A seasonal analysis points out improved consistency during December–January–February and September–October–November, in contrast to significantly worse correlations in March–April–May (in relation to the growing vegetation) and June–July–August (in relation to low SSM values &lt; 0.1 m3 m−3 and low spatial variability). (d) The combined use of the SURFEX (ISBA) SVAT model with the SAFRAN system, initialized with SMOS-L43.0 1 km disaggregated data, is proven to be a suitable tool for producing regional SM maps with high accuracy, which could be used as initial conditions for model simulations, flood forecasting, crop monitoring and crop development strategies, among others.

A Critical Evaluation on the Role of Aerodynamic and Canopy–Surface Conductance Parameterization in SEB and SVAT Models for Simulating Evapotranspiration: A Case Study in the Upper Biebrza National Park Wetland in Poland

Evapotranspiration (ET) estimation through the surface energy balance (SEB) and soil-vegetation-atmosphere-transfer (SVAT) models are uncertain due to the empirical parameterizations of the aerodynamic and canopy-substrate conductances (gA and gS) for heat and water vapor transfers. This study critically assessed the impact of conductance parameterizations on ET simulation using three structurally different SEB and SVAT models for an ecologically important North-Eastern European wetland, Upper Biebrza National Park (UBNP) in two consecutive years 2015 and 2016. A pronounced ET underestimation (mean bias −0.48 to −0.68 mm day−1) in SEBS (Surface Energy Balance System) was associated with an overestimation of gA due to uncertain parameterization of momentum roughness length and bare soil’s excess resistance to heat transfer (kB−1) under low vegetation cover. The systematic ET overestimation (0.65–0.80 mm day−1) in SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes) was attributed to the overestimation of both the conductances. Conductance parameterizations in SEBS and SCOPE appeared to be very sensitive to the general ecohydrological conditions, with a tendency of overestimating gA (gS) under humid (arid) conditions. Low ET bias in the analytical STIC (Surface Temperature Initiated Closure) model as compared to SEBS/SCOPE indicated the critical need for calibration-free conductance parameterizations for improved ET estimation.

The Impact of Climate Change on the Water Balance of Oil Sands Reclamation Covers and Natural Soil Profiles

Abstract The design of reclamation soil covers at oil sands mines in northern Alberta, Canada, has been conventionally based on the calibration of soil–vegetation–atmosphere transfer (SVAT) models against field monitoring observations collected over several years, followed by simulations of long-term performance using historical climate data. This paper evaluates the long-term water balances for reclamation covers on two oil sands landforms and three natural coarse-textured forest soil profiles using both historical climate data and future climate projections. Twenty-first century daily precipitation and temperature data from CanESM2 were downscaled based on three representative concentration pathways (RCPs) employing a stochastic weather generator [Long Ashton Research Station Weather Generator (LARS-WG)]. Relative humidity, wind speed, and net radiation were downscaled using the delta change method. Downscaled precipitation and estimated potential evapotranspiration were used as inputs to simulate soil water dynamics using physically based models. Probability distributions of growing season (April–October) actual evapotranspiration (AET) and net percolation (NP) for the baseline and future periods show that AET and NP at all sites are expected to increase throughout the twenty-first century regardless of RCP, time period, and soil profile. Greater increases in AET and NP are projected toward the end of the twenty-first century. The increases in future NP at the two reclamation covers are larger (as a percentage increase) than at most of the natural sites. Increases in NP will result in greater water yield to surface water and may accelerate the rate at which chemical constituents contained within mine waste are released to downstream receptors, suggesting these potential changes need to be considered in mine closure designs.

Climate change over the high-mountain versus plain areas: Effects on the land surface hydrologic budget in the Alpine area and northern Italy

Abstract. Climate change may intensify during the second half of the current century. Changes in temperature and precipitation can exert a significant impact on the regional hydrologic cycle. Because the land surface serves as the hub of interactions among the variables constituting the energy and water cycles, evaluating the land surface processes is essential to detail the future climate. In this study, we employ a trusted soil–vegetation–atmosphere transfer scheme, called the University of Torino model of land Processes Interaction with Atmosphere (UTOPIA), in offline simulations to quantify the changes in hydrologic components in the Alpine area and northern Italy, between the period of 1961–1990 and 2071–2100. The regional climate projections are obtained by the Regional Climate Model version 3 (RegCM3) via two emission scenarios – A2 and B2 from the Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios. The hydroclimate projections, especially from A2, indicate that evapotranspiration generally increases, especially over the plain areas, and consequently the surface soil moisture decreases during summer, falling below the wilting point threshold for an extra month. In the high-mountain areas, due to the earlier snowmelt, the land surface becomes snowless for an additional month. The annual mean number of dry (wet) days increases remarkably (slightly), thus increasing the risk of severe droughts, and slightly increasing the risk of floods coincidently. Our results have serious implications for human life, including agricultural production, water sustainability, and general infrastructures, over the Alpine and adjacent plain areas and can be used to plan the managements of water resources, floods, irrigation, forestry, hydropower, and many other relevant activities.