ALBERT

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: We develop data assimilation techniques for non-linear dynamical systems modelled by moving mesh methods. Such techniques are valuable for explicitly tracking interfaces and boundaries in evolving systems. The unique aspect of these assimilation techniques is that both the states of the system and the positions of the mesh points are updated simultaneously using physical observations. Covariances between states and mesh points are generated either by a correlation structure function in a variational context or by ensemble methods. The application of the techniques is demonstrated on a one-dimensional model of a grounded shallow ice sheet. It is shown, using observations of surface elevation and/or surface ice velocities, that the techniques predict the evolution of the ice sheet margin and the ice thickness accurately and efficiently. This approach also allows the straightforward assimilation of observations of the position of the ice sheet margin.

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

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

Description: This paper introduces an Ensemble Square Root Filter (EnSRF), a deterministic Ensemble Kalman Filter, to the context of assimilating jointly observations of surface soil moisture (SSM) and 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), 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, which has been routinely operated in LDAS-Monde, in a real case over the well-studied 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. SEKF and EnSRF have a similar behaviour for LAI showing performances that are influenced by the vegetation type. For SSM, EnSRF estimates tend to be closer to observations than SEKF. The impact of assimilating SSM and LAI is also assessed 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 exhibited, as expected, showing seasonal patterns that vary geographically. Moderate correlation and anti-correlations are also noticed between LAI and soil moisture in spring, summer and autumn, their absolute value tending 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 and gross primary production (GPP) as well as measures of river discharges from gauging stations. The EnSRF shows a systematic albeit moderate improvement for evapotranspiration and GPP and a highly positive impact on river discharges, while the SEKF exhibits a more contrasting performance.

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