This study aims to assess the potential of the LDAS-Monde platform, a land data assimilation system developed by Météo-France, to monitor the impact on vegetation state of the 2018 summer heatwave over Western Europe. The LDAS-Monde is driven by ECMWF’s (i) ERA5 reanalysis, and (ii) the Integrated Forecasting System High Resolution operational analysis (IFS-HRES), used in conjunction with the assimilation of Copernicus Global Land Service (CGLS) satellite-derived products, namely the Surface Soil Moisture (SSM) and the Leaf Area Index (LAI). The study of long time series of satellite derived CGLS LAI (2000–2018) and SSM (2008–2018) highlights marked negative anomalies for July 2018 affecting large areas of northwestern Europe and reflects the impact of the heatwave. Such large anomalies spreading over a large part of the domain of interest have never been observed in the LAI product over this 19-year period. LDAS-Monde land surface reanalyses were produced at spatial resolutions of 0.25° × 0.25° (January 2008 to October 2018) and 0.10° × 0.10° (April 2016 to December 2018). Both configurations of LDAS-Monde forced by either ERA5 or HRES capture well the vegetation state in general and for this specific event, with HRES configuration exhibiting better monitoring skills than ERA5 configuration. The consistency of ERA5- and IFS HRES-driven simulations over the common period (April 2016 to October 2018) allowed to disentangle and appreciate the origin of improvements observed between the ERA5 and HRES. Another experiment, down-scaling ERA5 to HRES spatial resolutions, was performed. Results suggest that land surface spatial resolution is key (e.g., associated to a better representation of the land cover, topography) and using HRES forcing still enhances the skill. While there are advantages in using HRES, there is added value in down-scaling ERA5, which can provide consistent, long term, high resolution land reanalysis. If the improvement from LDAS-Monde analysis on control variables (soil moisture from layers 2 to 8 of the model representing the first meter of soil and LAI) from the assimilation of SSM and LAI was expected, other model variables benefit from the assimilation through biophysical processes and feedback in the model. Finally, we also found added value of initializing 8-day land surface HRES driven forecasts from LDAS-Monde analysis when compared with model-only initial conditions.

Land data assimilation system (LDAS)-Monde, an offline land data assimilation system with global capacity, is applied over the CONtiguous US (CONUS) domain to enhance monitoring accuracy for water and energy states and fluxes. LDAS-Monde ingests satellite-derived surface soil moisture (SSM) and leaf area index (LAI) estimates to constrain the interactions between soil, biosphere, and atmosphere (ISBA) land surface model (LSM) coupled with the CNRM (Centre National de Recherches Météorologiques) version of the total runoff integrating pathways (CTRIP) continental hydrological system (ISBA-CTRIP). LDAS-Monde is forced by the ERA-5 atmospheric reanalysis from the European Center for Medium Range Weather Forecast (ECMWF) from 2010 to 2016 leading to a seven-year, quarter degree spatial resolution offline reanalysis of land surface variables (LSVs) over CONUS. The impact of assimilating LAI and SSM into LDAS-Monde is assessed over North America, by comparison to satellite-driven model estimates of land evapotranspiration from the Global Land Evaporation Amsterdam Model (GLEAM) project, and upscaled ground-based observations of gross primary productivity from the FLUXCOM project. Taking advantage of the relatively dense data networks over CONUS, we have also evaluated the impact of the assimilation against in situ measurements of soil moisture from the USCRN (US Climate Reference Network), together with river discharges from the United States Geological Survey (USGS) and the Global Runoff Data Centre (GRDC). Those data sets highlight the added value of assimilating satellite derived observations compared with an open-loop simulation (i.e., no assimilation). It is shown that LDAS-Monde has the ability not only to monitor land surface variables but also to forecast them, by providing improved initial conditions, which impacts persist through time. LDAS-Monde reanalysis also has the potential to be used to monitor extreme events like agricultural drought. Finally, limitations related to LDAS-Monde and current satellite-derived observations are exposed as well as several insights on how to use alternative datasets to analyze soil moisture and vegetation state.

LDAS-Monde is a global offline land data assimilation system (LDAS) that jointly assimilates satellite-derived observations of surface soil moisture (SSM) and leaf area index (LAI) into the ISBA (Interaction between Soil Biosphere and Atmosphere) land surface model (LSM). This study demonstrates that LDAS-Monde is able to detect, monitor and forecast the impact of extreme weather on land surface states. Firstly, LDAS-Monde is run globally at 0.25^∘ spatial resolution over 2010–2018. It is forced by the state-of-the-art ERA5 reanalysis (LDAS_ERA5) from the European Centre for Medium Range Weather Forecasts (ECMWF). The behaviour of the assimilation system is evaluated by comparing the analysis with the assimilated observations. Then the land surface variables (LSVs) are validated with independent satellite datasets of evapotranspiration, gross primary production, sun-induced fluorescence and snow cover. Furthermore, in situ measurements of SSM, evapotranspiration and river discharge are employed for the validation. Secondly, the global analysis is used to (i) detect regions exposed to extreme weather such as droughts and heatwave events and (ii) address specific monitoring and forecasting requirements of LSVs for those regions. This is performed by computing anomalies of the land surface states. They display strong negative values for LAI and SSM in 2018 for two regions: north-western Europe and the Murray–Darling basin in south-eastern Australia. For those regions, LDAS-Monde is forced with the ECMWF Integrated Forecasting System (IFS) high-resolution operational analysis (LDAS_HRES, 0.10^∘ spatial resolution) over 2017–2018. Monitoring capacities are studied by comparing open-loop and analysis experiments, again against the assimilated observations. Forecasting abilities are assessed by initializing 4 and 8 d LDAS_HRES forecasts of the LSVs with the LDAS_HRES assimilation run compared to the open-loop experiment. The positive impact of initialization from an analysis in forecast mode is particularly visible for LAI that evolves at a slower pace than SSM and is more sensitive to initial conditions than to atmospheric forcing, even at an 8 d lead time. This highlights the impact of initial conditions on LSV forecasts and the value of jointly analysing soil moisture and vegetation states.

This paper introduces an ensemble square root filter (EnSRF) in the context of jointly assimilating observations of surface soil moisture (SSM) and the leaf area index (LAI) in the Land Data Assimilation System LDAS-Monde. By ingesting those satellite-derived products, LDAS-Monde constrains the Interaction between Soil, Biosphere and Atmosphere (ISBA) land surface model (LSM), coupled with the CNRM (Centre National de Recherches Météorologiques) version of the Total Runoff Integrating Pathways (CTRIP) model to improve the reanalysis of land surface variables (LSVs). To evaluate its ability to produce improved LSVs reanalyses, the EnSRF is compared with the simplified extended Kalman filter (SEKF), which has been well studied within the LDAS-Monde framework. The comparison is carried out over the Euro-Mediterranean region at a 0.25^∘ spatial resolution between 2008 and 2017. Both data assimilation approaches provide a positive impact on SSM and LAI estimates with respect to the model alone, putting them closer to assimilated observations. The SEKF and the EnSRF have a similar behaviour for LAI showing performance levels that are influenced by the vegetation type. For SSM, EnSRF estimates tend to be closer to observations than SEKF values. The comparison between the two data assimilation approaches is also carried out on unobserved soil moisture in the other layers of soil. Unobserved control variables are updated in the EnSRF through covariances and correlations sampled from the ensemble linking them to observed control variables. In our context, a strong correlation between SSM and soil moisture in deeper soil layers is found, as expected, showing seasonal patterns that vary geographically. Moderate correlation and anti-correlations are also noticed between LAI and soil moisture, varying in space and time. Their absolute value, reaching their maximum in summer and their minimum in winter, tends to be larger for soil moisture in root-zone areas, showing that assimilating LAI can have an influence on soil moisture. Finally an independent evaluation of both assimilation approaches is conducted using satellite estimates of evapotranspiration (ET) and gross primary production (GPP) as well as measures of river discharges from gauging stations. The EnSRF shows a systematic albeit moderate improvement of root mean square differences (RMSDs) and correlations for ET and GPP products, but its main improvement is observed on river discharges with a high positive impact on Nash–Sutcliffe efficiency scores. Compared to the EnSRF, the SEKF displays a more contrasting performance.

Root zone soil moisture (RZSM) is critical for water resource management, drought monitoring and sub-seasonal flood climate prediction. While RZSM is not directly observable from space, several RZSM products are available and widely used at global and continental scales. This study conducts a comprehensive and quantitative evaluation of eight RZSM products using observations from 58 in situ soil moisture stations over the Huai River basin (HRB) in China. Attention is drawn to the potential factors that contribute to the uncertainties of model-based RZSM, including the errors in atmospheric forcing, vegetation parameterizations, soil properties and spatial scale mismatch. The results show that the Global Land Data Assimilation System Catchment Land Surface Model (GLDAS_CLSM) outperforms the other RZSM products with the highest correlation coefficient (R= 0.69) and the lowest unbiased root mean square error (ubRMSE = 0.018 m³ m⁻³), while SMOS Level 4 (L4) RZSM shows the worst performance among eight RZSM products. The RZSM products based on land surface models generally perform better in the wet season than in the dry season due to the enhanced ability to capture of the temporal dynamics of in situ observations in the wet season and the inertia of remaining high soil moisture values even in the dry season, while the SMOS L4 RZSM product, derived from SMOS L3 surface moisture (SSM) combined with an exponential filter method, performs better in the dry season due to the attenuated ground microwave radiation signal caused by the increased water vapour absorption and scattering in the wet season. The underestimated SMOS L3 SSM triggers the underestimation of RZSM in SMOS L4. The overestimated RZSM products based on land surface models could be associated with the overestimated precipitation amounts and frequency, the underestimated air temperature, and the underestimated ratio of transpiration to the total terrestrial evapotranspiration. In addition, the biased soil properties and flawed vegetation parameterizations affect the hydrothermal transport processes represented in different land surface models (LSMs) and lead to inaccurate soil moisture simulation. The scale mismatch between point and footprint also introduces representative errors. The comparison of frequency of normalized soil moisture between RZSM products and in situ observations indicates that the LSMs should focus on reducing the frequency of wet soil moisture, increasing the frequency of dry soil moisture and the ability to capture the frequency peak of soil moisture. The study provides some insights into how to improve the ability of land surface models to simulate the land surface states and fluxes by taking into account the issues mentioned above. Finally, these results can be extrapolated to other regions located in similar climate zones, as they share similar precipitation patterns that dominate the terrestrial water cycle.

This paper investigates to what extent soil moisture and vegetation density information can be extracted from the Advanced Scatterometer (ASCAT) satellite-derived radar backscatter (σ°) in a data assimilation context. The impact of independent estimates of the surface soil moisture (SSM) and leaf area index (LAI) of diverse vegetation types on ASCAT σ° observations is simulated over southwestern France using the water cloud model (WCM). The LAI and SSM variables used by the WCM are derived from satellite observations and from the Interactions between Soil, Biosphere, and Atmosphere (ISBA) land surface model, respectively. They permit the calibration of the four parameters of the WCM describing static soil and vegetation characteristics. A seasonal analysis of the model scores shows that the WCM has shortcomings over karstic areas and wheat croplands. In the studied area, the Klaus windstorm in January 2009 damaged a large fraction of the Landes forest. The ability of the WCM to represent the impact of Klaus and to simulate ASCAT σ° observations in contrasting land-cover conditions is explored. The difference in σ° observations between the forest zone affected by the storm and the bordering agricultural areas presents a marked seasonality before the storm. The difference is small in the springtime (from March to May) and large in the autumn (September to November) and wintertime (December to February). After the storm, hardly any seasonality was observed over four years. This study shows that the WCM is able to simulate this extreme event. It is concluded that the WCM could be used as an observation operator for the assimilation of ASCAT σ° observations into the ISBA land surface model.