ALBERT

Description: Soil moisture is a key parameter when it comes to understanding the processes related to the water cycle on continental surfaces (infiltration, evapotranspiration, runoff, etc [...]

Description: LDAS-Monde is a global land data assimilation system (LDAS) developed by Centre National de Recherches Météorologiques (CNRM) to monitor land surface variables (LSV) at various scales, from regional to global. With LDAS-Monde, it is possible to jointly assimilate satellite-derived observations of surface soil moisture (SSM) and leaf area index (LAI) into the interactions between soil biosphere and atmosphere (ISBA) land surface model (LSM) in order to analyze the soil moisture profile together with vegetation biomass. In this study, we investigate LDAS-Monde’s ability to predict LSV states up to two weeks in the future using atmospheric forecasts. In particular, the impact of the initialization, and the evolution of the forecasted variables in the LSM are addressed. LDAS-Monde is an offline system normally driven by atmospheric reanalysis, but in this study is forced by atmospheric forecasts from the European Centre for Medium-Range Weather Forecasts (ECMWF) for the 2017–2018 period over the contiguous United States (CONUS) at a 0.2° × 0.2° spatial resolution. These LSV forecasts are initialized either by the model alone (LDAS-Monde open-loop, without assimilation) or by the analysis (assimilation of SSM and LAI). These two forecasts are then evaluated using satellite-derived observations of SSM and LAI, evapotranspiration (ET) estimates, as well as in situ measurements of soil moisture from the U.S. Climate Reference Network (USCRN). Results indicate that for the three evaluation variables (SSM, LAI, and ET), LDAS-Monde provides reasonably accurate and consistent predictions two weeks in advance. Additionally, the initial conditions after assimilation are shown to make a positive impact with respect to LAI and ET. This impact persists in time for these two vegetation-related variables. Many model variables, such as SSM, root zone soil moisture (RZSM), LAI, ET, and drainage, remain relatively consistent as the forecast lead time increases, while runoff is highly variable.

Description: The high computational resources and the time-consuming IO (input/output) are major issues in offline ensemble-based high-dimensional data assimilation systems. Bearing these in mind, this study proposes a sophisticated dynamically running job scheme as well as an innovative parallel IO algorithm to reduce the time to solution of an offline framework for high-dimensional ensemble Kalman filters. The dynamically running job scheme runs as many tasks as possible within a single job to reduce the queuing time and minimize the overhead of starting and/or ending a job. The parallel IO algorithm reads or writes non-overlapping segments of multiple files with an identical structure to reduce the IO times by minimizing the IO competitions and maximizing the overlapping of the MPI (Message Passing Interface) communications with the IO operations. Results based on sensitive experiments show that the proposed parallel IO algorithm can significantly reduce the IO times and have a very good scalability, too. Based on these two advanced techniques, the offline and online modes of ensemble Kalman filters are built based on PDAF (Parallel Data Assimilation Framework) to comprehensively assess their efficiencies. It can be seen from the comparisons between the offline and online modes that the IO time only accounts for a small fraction of the total time with the proposed parallel IO algorithm. The queuing time might be less than the running time in a low-loaded supercomputer such as in an operational context, but the offline mode can be nearly as fast as, if not faster than, the online mode in terms of time to solution. However, the queuing time is dominant and several times larger than the running time in a high-loaded supercomputer. Thus, the offline mode is substantially faster than the online mode in terms of time to solution, especially for large-scale assimilation problems. From this point of view, results suggest that an offline ensemble Kalman filter with an efficient implementation and a high-performance parallel file system should be preferred over its online counterpart for intermittent data assimilation in many situations.

Description: 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.

Description: Earth observations were used to evaluate the representation of land surface temperature (LST) and vegetation coverage over Iberia in two state-of-the-art land surface models (LSMs) – the European Centre for Medium-Range Weather Forecasts (ECMWF) Carbon-Hydrology Tiled ECMWF Scheme for Surface Exchanges over Land (CHTESSEL) and the Météo-France Interaction between Soil Biosphere and Atmosphere model (ISBA) within the SURface EXternalisée modeling platform (SURFEX-ISBA) for the 2004–2015 period. The results showed that the daily maximum LST simulated by CHTESSEL over Iberia was affected by a large cold bias during summer months when compared against the Satellite Application Facility on Land Surface Analysis (LSA-SAF), reaching magnitudes larger than 10 ∘C over wide portions of central and southwestern Iberia. This error was shown to be tightly linked to a misrepresentation of the vegetation cover.  In contrast, SURFEX simulations did not display such a cold bias. We show that this was due to the better representation of vegetation cover in SURFEX, which uses an updated land cover dataset (ECOCLIMAP-II) and an interactive vegetation evolution, representing seasonality. The representation of vegetation over Iberia in CHTESSEL was improved by combining information from the European Space Agency Climate Change Initiative (ESA-CCI) land cover dataset with the Copernicus Global Land Service (CGLS) leaf area index (LAI) and fraction of vegetation coverage (FCOVER). The proposed improvement in vegetation also included a clumping approach that introduces seasonality to the vegetation cover. The results showed significant added value, removing the daily maximum LST summer cold bias completely, without reducing the accuracy of the simulated LST, regardless of season or time of the day. The striking performance differences between SURFEX and CHTESSEL were fundamental to guiding the developments in CHTESSEL highlighting the importance of using different models. This work has important implications: first, it takes advantage of LST, a key variable in surface–atmosphere energy and water exchanges, which is closely related to satellite top-of-atmosphere observations, to improve the model's representation of land surface processes. Second, CHTESSEL is the land surface model employed by ECMWF in the production of their weather forecasts and reanalysis; hence systematic errors in land surface variables and fluxes are then propagated into those products. Indeed, we showed that the summer daily maximum LST cold bias over Iberia in CHTESSEL is present in the widely used ECMWF fifth-generation reanalysis (ERA5). Finally, our results provided hints about the interaction between vegetation land–atmosphere exchanges, highlighting the relevance of the vegetation cover and respective seasonality in representing land surface temperature in both CHTESSEL and SURFEX. As a whole, this work demonstrated the added value of using multiple earth observation products for constraining and improving weather and climate simulations.

Description: 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.

Description: The aim of this study is to estimate surface soil moisture at a spatial resolution of 500 m and a temporal resolution of at least 6 days, by combining remote sensing data from Sentinel-1 and optical data from Sentinel-2 and MODIS (Moderate-Resolution Imaging Spectroradiometer). The proposed methodology is based on the change detection technique, applied to a series of measurements over a three-year period (2015 to 2018). The algorithm described here as “Soil Moisture Estimations from the Synergy of Sentinel-1 and optical sensors (SMES)” proposes different options, allowing information from vegetation densities and seasonal conditions to be taken into account. The output from this algorithm is a moisture index ranging between 0 and 1, with 0 corresponding to the driest soils and 1 to the wettest soils. This methodology has been tested at different test sites (South of France, Central Tunisia, Western Benin and Southwestern Niger), characterized by a wide range of different climatic conditions. The resulting surface soil moisture estimations are compared with in situ measurements and already existing satellite-derived soil moisture ASCAT (Advanced SCATterometer) products. They are found to be well correlated, for the African regions in particular (RMSE below 6 vol.%). This outcome indicates that the proposed algorithm can be used with confidence to estimate the surface soil moisture of a wide range of climatically different sites.

Description: The MARINE (Model of Anticipation of Runoff and INundations for Extreme events) hydrological model is a distributed model dedicated to flash flood simulation. Recent developments of the MARINE model are explored in this work. On one hand, transfers of water through the subsurface, formerly relying on water height, now take place in a homogeneous soil column based on the soil saturation degree (SSF model). On the other hand, the soil column is divided into two layers, which represent, respectively, the upper soil layer and the deep weathered rocks (SSF–DWF model). The aim of the present work is to assess the accuracy of these new representations for the simulation of soil moisture during flash flood events. An exploration of the various products available in the literature for soil moisture estimation is performed. The efficiency of the models for soil saturation degree simulation is estimated with respect to several products either at the local scale or spatially distributed: (i) the gridded soil moisture product provided by the operational modeling chain SAFRAN-ISBA-MODCOU; (ii) the gridded soil moisture product provided by the LDAS-Monde assimilation chain, which is based on the ISBA-A-gs land surface model and assimilating satellite derived data; (iii) the upper soil water content hourly measurements taken from the SMOSMANIA observation network; and (iv) the Soil Water Index provided by the Copernicus Global Land Service (CGLS), which is derived from Sentinel-1 C-SAR and ASCAT satellite data. The case study is performed over two French Mediterranean catchments impacted by flash flood events over the 2017–2019 period. The local comparison of the MARINE outputs with the SMOSMANIA measurements, as well as the comparison at the basin scale of the MARINE outputs with the gridded LDAS-Monde and CGLS data, lead to the following conclusion: both the dynamics and the amplitudes of the soil saturation degree simulated with the SSF and SSF–DWF models are better correlated with both the SMOSMANIA measurements and the LDAS-Monde data than the outputs of the base model. Finally, the soil saturation degree simulated by the two-layers model for the deep layer is compared to the soil saturation degree provided by the LDAS-Monde product at corresponding depths. In conclusion, the developments presented for the representation of subsurface flow in the MARINE model enhance the soil saturation degree simulation during flash floods with respect to both gridded data and local soil moisture measurements.