Soil-Vegetation-Atmosphere Transfer Models (SVAT) and Crop Simulation Models describe physical and physiological processes occurring in crop canopies. Remote sensing data may be used through ...assimilation procedures for constraining or driving SVAT and crop models. These models provide continuous simulation of processes such as evapotranspiration and, thus, direct means for interpolating evapotranspiration between remote sensing data acquisitions (which is not the case for classical evapotranspiration mapping methods). They also give access to variables other than evapotranspiration, such as soil moisture and crop production. We developed the coupling between crop, SVAT and radiative transfer models in order to implement assimilation procedures in various wavelength domains (solar, thermal and microwave). Such coupling makes it possible to transfer information from one model to another and then to use remote sensing information for retrieving model parameters which are not directly related to remote sensing data (such as soil initial water content, plant growth parameters, physical properties of soil and so on). Simple assimilation tests are presented to illustrate the main techniques that may be used for monitoring crop processes and evapotranspiration. An application to a small agricultural area is also performed showing the potential of such techniques for retrieving evapotranspiration and information on irrigation practices over wheat fields.
The streamflow of the Yellow River (YR) is strongly affected by human
activities like irrigation and dam operation. Many attribution studies
have focused on the long-term trends of streamflows, yet ...the contributions of
these anthropogenic factors to streamflow fluctuations have not been well
quantified with fully mechanistic models. This study aims to (1) demonstrate whether the mechanistic global land surface model ORCHIDEE (ORganizing Carbon and Hydrology
in Dynamic EcosystEms) is
able to simulate the streamflows of this complex rivers with human
activities using a generic parameterization for human activities and (2) preliminarily quantify the roles of irrigation and dam operation in
monthly streamflow fluctuations of the YR from 1982 to 2014 with a newly
developed irrigation module and an offline dam operation model.
Validations with observed streamflows near the outlet of the YR
demonstrated that model performances improved notably with incrementally
considering irrigation (mean square error (MSE) decreased by 56.9 %) and
dam operation (MSE decreased by another 30.5 %). Irrigation withdrawals
were found to substantially reduce the river streamflows by approximately
242.8±27.8×108 m3 yr−1 in line with independent
census data (231.4±31.6×108 m3 yr−1). Dam operation
does not change the mean streamflows in our model, but it impacts
streamflow seasonality, more than the seasonal change of precipitation.
By only considering generic operation schemes, our dam model is able to
reproduce the water storage changes of the two large reservoirs,
LongYangXia and LiuJiaXia (correlation coefficient of
∼ 0.9). Moreover, other commonly neglected factors, such as
the large operation contribution from multiple medium/small
reservoirs, the dominance of large irrigation districts for
streamflows (e.g., the Hetao Plateau), and special management policies
during extreme years, are highlighted in this study. Related
processes should be integrated into models to better project future YR water
resources under climate change and optimize adaption strategies.
This study investigates the impact of topography on five snow cover fraction (SCF) parameterizations developed for global climate models (GCMs), including two novel ones. The parameterization skill ...is first assessed with the High Mountain Asia Snow Reanalysis (HMASR), and three of them are implemented in the ORCHIDEE land surface model (LSM) and tested in global land–atmosphere coupled simulations. HMASR includes snow depth (SD) uncertainties, which may be due to the elevation differences between in situ stations and HMASR grid cells. Nevertheless, the SCF–SD relationship varies greatly between mountainous and flat areas in HMASR, especially during the snow-melting period. The new parameterizations that include a dependency on the subgrid topography allow a significant SCF bias reduction, reaching 5 % to 10 % on average in the global simulations over mountainous areas, which in turn leads to a reduction of the surface cold bias from −1.8 ∘C to about −1 ∘C in High Mountain Asia (HMA). Furthermore, the seasonal hysteresis between SCF and SD found in HMASR is better captured in the parameterizations that split the accumulation and the depletion curves or that include a dependency on the snow density. The deep-learning SCF parameterization is promising but exhibits more resolution-dependent and region-dependent features. Persistent snow cover biases remain in global land–atmosphere experiments. This suggests that other model biases may be intertwined with the snow biases and points out the need to continue improving snow models and their calibration. Increasing the model resolution does not consistently reduce the simulated SCF biases, although biases get narrower around mountain areas. This study highlights the complexity of calibrating SCF parameterizations since they affect various land–atmosphere feedbacks. In summary, this research spots the importance of considering topography in SCF parameterizations and the challenges in accurately representing snow cover in mountainous regions. It calls for further efforts to improve the representation of subgrid-scale processes affecting snowpack in climate models.
Land-atmosphere feedbacks, which are particularly important over the Sahel during the West African Monsoon (WAM), partly depend on a large range of processes linked to the land surface hydrology and ...the vegetation heterogeneities. This study focuses on the evaluation of a new land surface hydrology within the Noah-WRF land-atmosphere-coupled mesoscale model over the Sahel. This new hydrology explicitly takes account for the Dunne runoff using topographic information, the Horton runoff using a Green-Ampt approximation, and land surface heterogeneities. The previous and new versions of Noah-WRF are compared against a unique observation dataset located over the Dantiandou Kori (Niger). This dataset includes dense rain gauge network, surfaces temperatures estimated from MSG/SEVIRI data, surface soil moisture mapping based on ASAR/ENVISAT C-band radar data and in situ observations of surface atmospheric and land surface energy budget variables. Generally, the WAM is reasonably reproduced by Noah-WRF even if some limitations appear throughout the comparison between simulations and observations. An appreciable improvement of the model results is also found when the new hydrology is used. This fact seems to emphasize the relative importance of the representation of the land surface hydrological processes on the WAM simulated by Noah-WRF over the Sahel.
Land surface models (LSMs) use the atmospheric grid as their basic spatial decomposition because their main objective is to provide the lower boundary conditions to the atmosphere. Lateral water ...flows at the surface on the other hand require a much higher spatial discretization as they are closely linked to topographic details. We propose here a methodology to automatically tile the atmospheric grid into hydrological coherent units which are connected through a graph. As water is transported on sub-grids of the LSM, land variables can easily be transferred to the routing network and advected if needed. This is demonstrated here for temperature. The quality of the river networks generated, as represented by the connected hydrological transfer units, are compared to the original data in order to quantify the degradation introduced by the discretization method. The conditions the sub-grid elements impose on the time step of the water transport scheme are evaluated, and a methodology is proposed to find an optimal value. Finally the scheme is applied in an off-line version of the ORCHIDEE (Organising Carbon and Hydrology In Dynamic Ecosystems) LSM over Europe to show that realistic river discharge and temperatures are predicted over the major catchments of the region. The simulated solutions are largely independent of the atmospheric grid used thanks to the proposed sub-grid approach.
Plant activity in semi-arid ecosystems is largely controlled by pulses of precipitation, making them particularly vulnerable to increased aridity that is expected with climate change. Simple ...bucket-model hydrology schemes in land surface models (LSMs) have had limited ability in accurately capturing semi-arid water stores and fluxes. Recent, more complex, LSM hydrology models have not been widely evaluated against semi-arid ecosystem in situ data. We hypothesize that the failure of older LSM versions to represent evapotranspiration, ET, in arid lands is because simple bucket models do not capture realistic fluctuations in upper-layer soil moisture. We therefore predict that including a discretized soil hydrology scheme based on a mechanistic description of moisture diffusion will result in an improvement in model ET when compared to data because the temporal variability of upper-layer soil moisture content better corresponds to that of precipitation inputs. To test this prediction, we compared ORCHIDEE LSM simulations from (1) a simple conceptual 2-layer bucket scheme with fixed hydraulic parameters and (2) an 11-layer discretized mechanistic scheme of moisture diffusion in unsaturated soil based on Richards equations, against daily and monthly soil moisture and ET observations, together with data-derived estimates of transpiration / evapotranspiration, T∕ET, ratios, from six semi-arid grass, shrub, and forest sites in the south-western USA. The 11-layer scheme also has modified calculations of surface runoff, water limitation, and resistance to bare soil evaporation, E, to be compatible with the more complex hydrology configuration. To diagnose remaining discrepancies in the 11-layer model, we tested two further configurations: (i) the addition of a term that captures bare soil evaporation resistance to dry soil; and (ii) reduced bare soil fractional vegetation cover. We found that the more mechanistic 11-layer model results in a better representation of the daily and monthly ET observations. We show that, as predicted, this is because of improved simulation of soil moisture in the upper layers of soil (top ∼ 10 cm). Some discrepancies between observed and modelled soil moisture and ET may allow us to prioritize future model development and the collection of additional data. Biases in winter and spring soil moisture at the forest sites could be explained by inaccurate soil moisture data during periods of soil freezing and/or underestimated snow forcing data. Although ET is generally well captured by the 11-layer model, modelled T∕ET ratios were generally lower than estimated values across all sites, particularly during the monsoon season. Adding a soil resistance term generally decreased simulated bare soil evaporation, E, and increased soil moisture content, thus increasing transpiration, T, and reducing the negative bias between modelled and estimated monsoon T∕ET ratios. This negative bias could also be accounted for at the low-elevation sites by decreasing the model bare soil fraction, thus increasing the amount of transpiring leaf area. However, adding the bare soil resistance term and decreasing the bare soil fraction both degraded the model fit to ET observations. Furthermore, remaining discrepancies in the timing of the transition from minimum T∕ET ratios during the hot, dry May–June period to high values at the start of the monsoon in July–August may also point towards incorrect modelling of leaf phenology and vegetation growth in response to monsoon rains. We conclude that a discretized soil hydrology scheme and associated developments improve estimates of ET by allowing the modelled upper-layer soil moisture to more closely match the pulse precipitation dynamics of these semi-arid ecosystems; however, the partitioning of T from E is not solved by this modification alone.
Soil moisture is a key variable of land surface hydrology, and its correct representation in land surface models is crucial for local to global climate predictions. The errors may come from the model ...itself (structure and parameterization) but also from the meteorological forcing used. In order to separate the two source of errors, four atmospheric forcing datasets, GSWP3 (Global Soil Wetness Project Phase 3), PGF (Princeton Global meteorological Forcing), CRU-NCEP (Climatic Research Unit-National Center for Environmental Prediction), and WFDEI (WATCH Forcing Data methodology applied to ERA-Interim reanalysis data), were used to drive simulations in China by the land surface model ORCHIDEE-MICT(ORganizing Carbon and Hydrology in Dynamic EcosystEms: aMeliorated Interactions between Carbon and Temperature). Simulated soil moisture was compared with in situ and satellite datasets at different spatial and temporal scales in order to (1) estimate the ability of ORCHIDEE-MICT to represent soil moisture dynamics in China; (2) demonstrate the most suitable forcing dataset for further hydrological studies in Yangtze and Yellow River basins; and (3) understand the discrepancies of simulated soil moisture among simulations. Results showed that ORCHIDEE-MICT can simulate reasonable soil moisture dynamics in China, but the quality varies with forcing data. Simulated soil moisture driven by GSWP3 and WFDEI shows the best performance according to the root mean square error (RMSE) and correlation coefficient, respectively, suggesting that both GSWP3 and WFDEI are good choices for further hydrological studies in the two catchments. The mismatch between simulated and observed soil moisture is mainly explained by the bias of magnitude, suggesting that the parameterization in ORCHIDEE-MICT should be revised for further simulations in China. Underestimated soil moisture in the North China Plain demonstrates possible significant impacts of human activities like irrigation on soil moisture variation, which was not considered in our simulations. Finally, the discrepancies of meteorological variables and simulated soil moisture among the four simulations are analyzed. The result shows that the discrepancy of soil moisture is mainly explained by differences in precipitation frequency and air humidity rather than differences in precipitation amount.
The existing medium-resolution land cover time series produced under the European Space Agency's Climate Change Initiative provides 29 years (1992–2020) of annual land cover maps at 300 m resolution, ...allowing for a detailed study of land change dynamics over the contemporary era. Because models need two-dimensional parameters rather than two-dimensional land cover information, the land cover classes must be converted into model-appropriate plant functional types (PFTs) to apply this time series to Earth system and land surface models. The first-generation cross-walking table that was presented with the land cover product prescribed pixel-level PFT fractional compositions that varied by land cover class but that lacked spatial variability. Here we describe a new ready-to-use data product for climate modelling: spatially explicit annual maps of PFT fractional composition at 300 m resolution for 1992–2020, created by fusing the 300 m medium-resolution land cover product with several existing high-resolution datasets using a globally consistent method. In the resulting data product, which has 14 layers for each of the 29 years, pixel values at 300 m resolution indicate the percentage cover (0 %–100 %) for each of 14 PFTs, with pixel-level PFT composition exhibiting significant intra-class spatial variability at the global scale. We additionally present an updated version of the user tool that allows users to modify the baseline product (e.g. re-mapping, re-projection, PFT conversion, and spatial sub-setting) to meet individual needs. Finally, these new PFT maps have been used in two land surface models – Organising Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) and the Joint UK Land Environment Simulator (JULES) – to demonstrate their benefit over the conventional maps based on a generic cross-walking table. Regional changes in the fractions of trees, short vegetation, and bare-soil cover induce changes in surface properties, such as the albedo, leading to significant changes in surface turbulent fluxes, temperature, and vegetation carbon stocks. The dataset is accessible at https://doi.org/10.5285/26a0f46c95ee4c29b5c650b129aab788 (Harper et al., 2023).
The contribution of forests to carbon storage and biodiversity conservation highlights the need for accurate forest height and biomass mapping and monitoring. In France, forests are managed mainly by ...private owners and divided into small stands, requiring 10 to 50 m spatial resolution data to be correctly separated. Further, 35 % of the French forest territory is covered by mountains and Mediterranean forests which are managed very extensively. In this work, we used a deep-learning model based on multi-stream remote-sensing measurements (NASA's Global Ecosystem Dynamics Investigation (GEDI) lidar mission and ESA's Copernicus Sentinel-1 and Sentinel-2 satellites) to create a 10 m resolution canopy height map of France for 2020 (FORMS-H). In a second step, with allometric equations fitted to the French National Forest Inventory (NFI) plot data, we created a 30 m resolution above-ground biomass density (AGBD) map (Mg ha−1) of France (FORMS-B). Extensive validation was conducted. First, independent datasets from airborne laser scanning (ALS) and NFI data from thousands of plots reveal a mean absolute error (MAE) of 2.94 m for FORMS-H, which outperforms existing canopy height models. Second, FORMS-B was validated using two independent forest inventory datasets from the Renecofor permanent forest plot network and from the GLORIE forest inventory with MAE of 59.6 and 19.6 Mg ha−1, respectively, providing greater performance than other AGBD products sampled over France. Finally, we compared FORMS-V (for volume) with wood volume estimations at the ecological region scale and obtained an R2 of 0.63 with an MAE of 30 m3 ha−1. These results highlight the importance of coupling remote-sensing technologies with recent advances in computer science to bring material insights to climate-efficient forest management policies. Additionally, our approach is based on open-access data having global coverage and a high spatial and temporal resolution, making the maps reproducible and easily scalable. FORMS products can be accessed from https://doi.org/10.5281/zenodo.7840108 (Schwartz et al., 2023).
This paper presents an original methodology to retrieve surface (<5 cm) soil moisture over low vegetated regions using the two active microwave instruments of ERS satellites. The developed algorithm ...takes advantage of the multi-angular configuration and high temporal resolution of the Wind Scatterometer (WSC) combined with the SAR high spatial resolution. As a result, a mixed target model is proposed. The WSC backscattered signal may be represented as a combination of the vegetation and bare soil contributions weighted by their respective fractional covers. Over our temperate regions and time periods of interest, the vegetation signal is assumed to be principally due to forests backscattered signal. Then, thanks to the high spatial resolution of the SAR instrument, the forest contribution may be quantified from the analysis of the SAR image, and then removed from the total WSC signal in order to estimate the soil contribution. Finally, the Integral Equation Model (IEM, IEEE Transactions on Geoscience and Remote Sensing, 30 (2), (1992) 356) is used to estimate the effect of surface roughness and to retrieve surface soil moisture from the WSC multi-angular measurements. This methodology has been developed and applied on ERS data acquired over three different Seine river watersheds in France, and for a 3-year time period. The soil moisture estimations are compared with in situ ground measurements. High correlations (
R
2 greater than 0.8) are observed for the three study watersheds with a root mean square (rms) error smaller than 4%.