Data assimilation of SMAP soil moisture into the PCR-GLOBWB hydrological model to improve discharge estimates via A Novel Local Particle Filter

Author(s):  
Zhenwu Wang ◽  
Rolf Hut ◽  
Natthachet Tangdamrongsub ◽  
Nick van de Giesen

<p>Assimilating surface soil moisture data or GRACE data, retrieved from satellite, into hydrological models has been proven to improve the accuracy of hydrological model estimations and predictions. For data assimilation applications in hydrology, the ensemble Kalm filter(EnKF) is the most commonly used data assimilation(DA) method. Particle filters are a type of non-Gaussian filter that doesn’t need the normality assumption that the EnKF needs. Adding localization defeats the curse of dimensionality that is a problem in normal particle filters. In the present study, we investigated our adaption of the local particle filter based on the Gamma test theory(LPF-GT) to improve discharge estimates by assimilating SMAP satellite soil moisture into the PCR-GLOBWB hydrological model. The study area is the Rhine river basin, driven by forcing data from April 2015 to December 2016. The improved discharge estimates are obtained by using DA to adjust the surface soil moisture in the model. The influence of DA to discharge is not direct but works through the dynamics of the hydrological model.  To explore the potential of LPF-GT, serval sensitivity experiments were conducted to figure out the impact of localization scales and the number of particles on DA's performance. The DA estimates were validated against in situ discharge measurements from gauge stations. To demonstrate the benefit of LPF-GT, EnKF was used as a benchmark in this research. Increases in Nash-Sutcliffe (0.05%– 38%) and decreases in normalized RMSE (0.02%–3.4%) validated the capability of LPF-GT. Results showed that localization scales' impact was substantial. The optimal value of the localization scale was obtained by tuning. LPF-GT achieved a satisfactory performance when only using a few particles, even with as little as five particles. The sample errors posed an adverse impact on the open-loop results. Further improvement could be achieved by considering reduce sample errors due to a small number of particles.</p>

2015 ◽  
Vol 19 (12) ◽  
pp. 4831-4844 ◽  
Author(s):  
C. Draper ◽  
R. Reichle

Abstract. A 9 year record of Advanced Microwave Scanning Radiometer – Earth Observing System (AMSR-E) soil moisture retrievals are assimilated into the Catchment land surface model at four locations in the US. The assimilation is evaluated using the unbiased mean square error (ubMSE) relative to watershed-scale in situ observations, with the ubMSE separated into contributions from the subseasonal (SMshort), mean seasonal (SMseas), and inter-annual (SMlong) soil moisture dynamics. For near-surface soil moisture, the average ubMSE for Catchment without assimilation was (1.8 × 10−3 m3 m−3)2, of which 19 % was in SMlong, 26 % in SMseas, and 55 % in SMshort. The AMSR-E assimilation significantly reduced the total ubMSE at every site, with an average reduction of 33 %. Of this ubMSE reduction, 37 % occurred in SMlong, 24 % in SMseas, and 38 % in SMshort. For root-zone soil moisture, in situ observations were available at one site only, and the near-surface and root-zone results were very similar at this site. These results suggest that, in addition to the well-reported improvements in SMshort, assimilating a sufficiently long soil moisture data record can also improve the model representation of important long-term events, such as droughts. The improved agreement between the modeled and in situ SMseas is harder to interpret, given that mean seasonal cycle errors are systematic, and systematic errors are not typically targeted by (bias-blind) data assimilation. Finally, the use of 1-year subsets of the AMSR-E and Catchment soil moisture for estimating the observation-bias correction (rescaling) parameters is investigated. It is concluded that when only 1 year of data are available, the associated uncertainty in the rescaling parameters should not greatly reduce the average benefit gained from data assimilation, although locally and in extreme years there is a risk of increased errors.


2014 ◽  
Vol 18 (1) ◽  
pp. 173-192 ◽  
Author(s):  
A. L. Barbu ◽  
J.-C. Calvet ◽  
J.-F. Mahfouf ◽  
S. Lafont

Abstract. The land monitoring service of the European Copernicus programme has developed a set of satellite-based biogeophysical products, including surface soil moisture (SSM) and leaf area index (LAI). This study investigates the impact of joint assimilation of remotely sensed SSM derived from Advanced Scatterometer (ASCAT) backscatter data and the Copernicus Global Land GEOV1 satellite-based LAI product into the the vegetation growth version of the Interactions between Soil Biosphere Atmosphere (ISBA-A-gs) land surface model within the the externalised surface model (SURFEX) modelling platform of Météo-France. The ASCAT data were bias corrected with respect to the model climatology by using a seasonal-based CDF (Cumulative Distribution Function) matching technique. A multivariate multi-scale land data assimilation system (LDAS) based on the extended Kalman Filter (EKF) is used for monitoring the soil moisture, terrestrial vegetation, surface carbon and energy fluxes across the domain of France at a spatial resolution of 8 km. Each model grid box is divided into a number of land covers, each having its own set of prognostic variables. The filter algorithm is designed to provide a distinct analysis for each land cover while using one observation per grid box. The updated values are aggregated by computing a weighted average. In this study, it is demonstrated that the assimilation scheme works effectively within the ISBA-A-gs model over a four-year period (2008–2011). The EKF is able to extract useful information from the data signal at the grid scale and distribute the root-zone soil moisture and LAI increments throughout the mosaic structure of the model. The impact of the assimilation on the vegetation phenology and on the water and carbon fluxes varies from one season to another. The spring drought of 2011 is an interesting case study of the potential of the assimilation to improve drought monitoring. A comparison between simulated and in situ soil moisture gathered at the twelve SMOSMANIA (Soil Moisture Observing System–Meteorological Automatic Network Integrated Application) stations shows improved anomaly correlations for eight stations.


2013 ◽  
Vol 10 (7) ◽  
pp. 9057-9103 ◽  
Author(s):  
A. L. Barbu ◽  
J.-C. Calvet ◽  
J.-F. Mahfouf ◽  
S. Lafont

Abstract. The land monitoring service of the European Copernicus programme has developed a set of satellite-based biogeophysical products, including surface soil moisture (SSM) and leaf area index (LAI). This study investigates the impact of joint assimilation of remotely sensed SSM derived from ASCAT backscatter data and the GEOV1 satellite-based LAI into the ISBA-A-gs land surface model within the SURFEX modelling platform of Meteo-France. The ASCAT data were bias corrected with respect to the model climatology by using a seasonal-based CDF (Cumulative Distribution Function) matching technique. A multivariate multi-scale land data assimilation system (LDAS) based on the Extended Kalman Filter (EKF) is used for monitoring the soil moisture, terrestrial vegetation, surface carbon and energy fluxes across the France domain at a spatial resolution of 8 km. Each model grid box is divided in a number of land covers, each having its own set of prognostic variables. The filter algorithm is designed to provide a distinct analysis for each land cover while using one observation per grid box. The updated values are aggregated by computing a weighted average. In this study, it is demonstrated that the assimilation scheme works effectively within the ISBA-A-gs model over a four-year period (2008–2011). The EKF is able to extract useful information from the data signal at the grid scale and to distribute the root-zone soil moisture and LAI increments among the mosaic structure of the model. The impact of the assimilation on the vegetation phenology and on the water and carbon fluxes varies from one season to another. The spring drought of 2011 is an interesting case study showing the potential of the assimilation to improve drought monitoring. A comparison between simulated and in situ soil moisture gathered at the twelve SMOSMANIA stations shows improved anomaly correlations for eight stations.


2006 ◽  
Vol 7 (3) ◽  
pp. 421-432 ◽  
Author(s):  
Wade T. Crow ◽  
Emiel Van Loon

Abstract Data assimilation approaches require some type of state forecast error covariance information in order to optimally merge model predictions with observations. The ensemble Kalman filter (EnKF) dynamically derives such information through a Monte Carlo approach and the introduction of random noise in model states, fluxes, and/or forcing data. However, in land data assimilation, relatively little guidance exists concerning strategies for selecting the appropriate magnitude and/or type of introduced model noise. In addition, little is known about the sensitivity of filter prediction accuracy to (potentially) inappropriate assumptions concerning the source and magnitude of modeling error. Using a series of synthetic identical twin experiments, this analysis explores the consequences of making incorrect assumptions concerning the source and magnitude of model error on the efficiency of assimilating surface soil moisture observations to constrain deeper root-zone soil moisture predictions made by a land surface model. Results suggest that inappropriate model error assumptions can lead to circumstances in which the assimilation of surface soil moisture observations actually degrades the performance of a land surface model (relative to open-loop assimilations that lack a data assimilation component). Prospects for diagnosing such circumstances and adaptively correcting the culpable model error assumptions using filter innovations are discussed. The dual assimilation of both runoff (from streamflow) and surface soil moisture observations appears to offer a more robust assimilation framework where incorrect model error assumptions are more readily diagnosed via filter innovations.


2017 ◽  
Vol 21 (4) ◽  
pp. 2015-2033 ◽  
Author(s):  
David Fairbairn ◽  
Alina Lavinia Barbu ◽  
Adrien Napoly ◽  
Clément Albergel ◽  
Jean-François Mahfouf ◽  
...  

Abstract. This study evaluates the impact of assimilating surface soil moisture (SSM) and leaf area index (LAI) observations into a land surface model using the SAFRAN–ISBA–MODCOU (SIM) hydrological suite. SIM consists of three stages: (1) an atmospheric reanalysis (SAFRAN) over France, which forces (2) the three-layer ISBA land surface model, which then provides drainage and runoff inputs to (3) the MODCOU hydro-geological model. The drainage and runoff outputs from ISBA are validated by comparing the simulated river discharge from MODCOU with over 500 river-gauge observations over France and with a subset of stations with low-anthropogenic influence, over several years. This study makes use of the A-gs version of ISBA that allows for physiological processes. The atmospheric forcing for the ISBA-A-gs model underestimates direct shortwave and long-wave radiation by approximately 5 % averaged over France. The ISBA-A-gs model also substantially underestimates the grassland LAI compared with satellite retrievals during winter dormancy. These differences result in an underestimation (overestimation) of evapotranspiration (drainage and runoff). The excess runoff flowing into the rivers and aquifers contributes to an overestimation of the SIM river discharge. Two experiments attempted to resolve these problems: (i) a correction of the minimum LAI model parameter for grasslands and (ii) a bias-correction of the model radiative forcing. Two data assimilation experiments were also performed, which are designed to correct random errors in the initial conditions: (iii) the assimilation of LAI observations and (iv) the assimilation of SSM and LAI observations. The data assimilation for (iii) and (iv) was done with a simplified extended Kalman filter (SEKF), which uses finite differences in the observation operator Jacobians to relate the observations to the model variables. Experiments (i) and (ii) improved the median SIM Nash scores by about 9 % and 18 % respectively. Experiment (iii) reduced the LAI phase errors in ISBA-A-gs but had little impact on the discharge Nash efficiency of SIM. In contrast, experiment (iv) resulted in spurious increases in drainage and runoff, which degraded the median discharge Nash efficiency by about 7 %. The poor performance of the SEKF originates from the observation operator Jacobians. These Jacobians are dampened when the soil is saturated and when the vegetation is dormant, which leads to positive biases in drainage and/or runoff and to insufficient corrections during winter, respectively. Possible ways to improve the model are discussed, including a new multi-layer diffusion model and a more realistic response of photosynthesis to temperature in mountainous regions. The data assimilation should be advanced by accounting for model and forcing uncertainties.


2007 ◽  
Vol 8 (2) ◽  
pp. 194-206 ◽  
Author(s):  
Joaquín Muñoz Sabater ◽  
Lionel Jarlan ◽  
Jean-Christophe Calvet ◽  
François Bouyssel ◽  
Patricia De Rosnay

Abstract Root-zone soil moisture constitutes an important variable for hydrological and weather forecast models. Microwave radiometers like the L-band instrument on board the European Space Agency’s (ESA) future Soil Moisture and Ocean Salinity (SMOS) mission are being designed to provide estimates of near-surface soil moisture (0–5 cm). This quantity is physically related to root-zone soil moisture through diffusion processes, and both surface and root-zone soil layers are commonly simulated by land surface models (LSMs). Observed time series of surface soil moisture may be used to analyze the root-zone soil moisture using data assimilation systems. In this paper, various assimilation techniques derived from Kalman filters (KFs) and variational methods (VAR) are implemented and tested. The objective is to correct the modeled root-zone soil moisture deficiencies of the newest version of the Interaction between Soil, Biosphere, and Atmosphere scheme (ISBA) LSM, using the observations of the surface soil moisture of the Surface Monitoring of the Soil Reservoir Experiment (SMOSREX) over a 4-yr period (2001–04). This time period includes contrasting climatic conditions. Among the different algorithms, the ensemble Kalman filter (EnKF) and a simplified one-dimensional variational data assimilation (1DVAR) show the best performances. The lower computational cost of the 1DVAR is an advantage for operational root-zone soil moisture analysis based on remotely sensed surface soil moisture observations at a global scale.


2017 ◽  
Vol 44 ◽  
pp. 89-100 ◽  
Author(s):  
Luca Cenci ◽  
Luca Pulvirenti ◽  
Giorgio Boni ◽  
Marco Chini ◽  
Patrick Matgen ◽  
...  

Abstract. The assimilation of satellite-derived soil moisture estimates (soil moisture–data assimilation, SM–DA) into hydrological models has the potential to reduce the uncertainty of streamflow simulations. The improved capacity to monitor the closeness to saturation of small catchments, such as those characterizing the Mediterranean region, can be exploited to enhance flash flood predictions. When compared to other microwave sensors that have been exploited for SM–DA in recent years (e.g. the Advanced SCATterometer – ASCAT), characterized by low spatial/high temporal resolution, the Sentinel 1 (S1) mission provides an excellent opportunity to monitor systematically soil moisture (SM) at high spatial resolution and moderate temporal resolution. The aim of this research was thus to evaluate the impact of S1-based SM–DA for enhancing flash flood predictions of a hydrological model (Continuum) that is currently exploited for civil protection applications in Italy. The analysis was carried out in a representative Mediterranean catchment prone to flash floods, located in north-western Italy, during the time period October 2014–February 2015. It provided some important findings: (i) revealing the potential provided by S1-based SM–DA for improving discharge predictions, especially for higher flows; (ii) suggesting a more appropriate pre-processing technique to be applied to S1 data before the assimilation; and (iii) highlighting that even though high spatial resolution does provide an important contribution in a SM–DA system, the temporal resolution has the most crucial role. S1-derived SM maps are still a relatively new product and, to our knowledge, this is the first work published in an international journal dealing with their assimilation within a hydrological model to improve continuous streamflow simulations and flash flood predictions. Even though the reported results were obtained by analysing a relatively short time period, and thus should be supported by further research activities, we believe this research is timely in order to enhance our understanding of the potential contribution of the S1 data within the SM–DA framework for flash flood risk mitigation.


Sign in / Sign up

Export Citation Format

Share Document