scholarly journals Dynamic Calibration with an Ensemble Kalman Filter Based Data Assimilation Approach for Root-Zone Moisture Predictions

2007 ◽  
Vol 8 (4) ◽  
pp. 910-921 ◽  
Author(s):  
Nicola Montaldo ◽  
John D. Albertson ◽  
Marco Mancini

Abstract In the presence of uncertain initial conditions and soil hydraulic properties, land surface model (LSM) performance can be significantly improved by the assimilation of periodic observations of certain state variables, such as the near-surface soil moisture (θg), as observed from a remote platform. In this paper the possibility of merging observations and the model optimally for providing robust predictions of root-zone soil moisture (θ2) is demonstrated. An assimilation approach that assimilates θg through the ensemble Kalman filter (EnKF) and provides a physics-based update of θ2 is developed. This approach, as with other common soil moisture assimilation approaches, may fail when a key LSM parameter, for example, the saturated hydraulic conductivity (ks), is estimated poorly. This leads to biased model errors producing a violation of a main assumption (model errors with zero mean) of the EnKF. For overcoming this model bias an innovative assimilation approach is developed that accepts this violation in the early model run times and dynamically calibrates all the components of the ks ensemble as a function of the persistent bias in root-zone soil moisture, allowing one to remove the model bias, restore the fidelity to the EnKF requirements, and reduce the model uncertainty. The robustness of the proposed approach is also examined in sensitivity analyses.

2012 ◽  
Vol 16 (1) ◽  
pp. 105-119 ◽  
Author(s):  
B. Li ◽  
D. Toll ◽  
X. Zhan ◽  
B. Cosgrove

Abstract. Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of near surface soil moisture in a footprint area, and can be used to reduce bias of model estimates (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals, as they are without any scaling or pre-processing, to improve the estimated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed.


2011 ◽  
Vol 8 (4) ◽  
pp. 8131-8171 ◽  
Author(s):  
B. Li ◽  
D. Toll ◽  
X. Zhan ◽  
B. Cosgrove

Abstract. Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of soil moisture in a footprint area, and can be used to reduce model bias (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce model bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate the actual value of Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals to improve the mean of simulated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed.


2011 ◽  
Vol 8 (1) ◽  
pp. 1433-1468
Author(s):  
S. Nie ◽  
J. Zhu ◽  
Y. Luo

Abstract. The performance of the ensemble Kalman filter (EnKF) in soil moisture assimilation applications is investigated in the context of simultaneous state-parameter estimation in the presence of uncertainties from model parameters, initial soil moisture condition and atmospheric forcing. A physically-based land surface model is used for this purpose. Using a series of idealized twin experiments, model generated near-surface soil moisture observations are assimilated to estimate soil moisture state and three hydraulic parameters (the saturated hydraulic conductivity, the saturated soil moisture suction and a soil texture empirical parameter) in the model. The single imperfect parameter can be successfully estimated using the EnKF. Results show that all the three estimated parameters converge toward their respective true values, while the root mean squared errors (RMSE) of soil moisture associated with these parameters is on average reduced by 54% and 53% comparing with the non-parameter-estimation benchmark RMSE for near-surface layer and root zone layer, respectively. The performance of simultaneous multi-parameter estimation is significant degraded, mainly because the inherent balance relationship of these parameters is broken and the degree of freedom increases in assimilation processes. By introducing constraints between estimated parameters, the performance of the constraint-based simultaneous multi-parameter estimations are as good as that of single-parameter cases even assimilating temporal-sparse observations. In terms of the relative root mean squared error (RRE), the constraint-based estimation cases can achieve 36% to 53% in near-surface layer and 25% to 50% in root zone layer for different assimilation intervals ranging from 1-day to 40-days. This result suggests that the greatest advantage of this method can be displayed with a proper temporal-sparse assimilation interval of 10-days as actual measurement interval of conventional in situ soil moisture observations. As these obtained constraints are mostly in statistical sense, this constraint-based simultaneous state-parameter estimation scheme is supposed to be suitable for other land surface models in soil moisture assimilation applications.


2016 ◽  
Vol 20 (12) ◽  
pp. 4895-4911 ◽  
Author(s):  
Gabriëlle J. M. De Lannoy ◽  
Rolf H. Reichle

Abstract. Three different data products from the Soil Moisture Ocean Salinity (SMOS) mission are assimilated separately into the Goddard Earth Observing System Model, version 5 (GEOS-5) to improve estimates of surface and root-zone soil moisture. The first product consists of multi-angle, dual-polarization brightness temperature (Tb) observations at the bottom of the atmosphere extracted from Level 1 data. The second product is a derived SMOS Tb product that mimics the data at a 40° incidence angle from the Soil Moisture Active Passive (SMAP) mission. The third product is the operational SMOS Level 2 surface soil moisture (SM) retrieval product. The assimilation system uses a spatially distributed ensemble Kalman filter (EnKF) with seasonally varying climatological bias mitigation for Tb assimilation, whereas a time-invariant cumulative density function matching is used for SM retrieval assimilation. All assimilation experiments improve the soil moisture estimates compared to model-only simulations in terms of unbiased root-mean-square differences and anomaly correlations during the period from 1 July 2010 to 1 May 2015 and for 187 sites across the US. Especially in areas where the satellite data are most sensitive to surface soil moisture, large skill improvements (e.g., an increase in the anomaly correlation by 0.1) are found in the surface soil moisture. The domain-average surface and root-zone skill metrics are similar among the various assimilation experiments, but large differences in skill are found locally. The observation-minus-forecast residuals and analysis increments reveal large differences in how the observations add value in the Tb and SM retrieval assimilation systems. The distinct patterns of these diagnostics in the two systems reflect observation and model errors patterns that are not well captured in the assigned EnKF error parameters. Consequently, a localized optimization of the EnKF error parameters is needed to further improve Tb or SM retrieval assimilation.


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.


2018 ◽  
Vol 22 (12) ◽  
pp. 6611-6626 ◽  
Author(s):  
Sara Sadri ◽  
Eric F. Wood ◽  
Ming Pan

Abstract. Since April 2015, NASA's Soil Moisture Active Passive (SMAP) mission has monitored near-surface soil moisture, mapping the globe (between 85.044∘ N/S) using an L-band (1.4 GHz) microwave radiometer in 2–3 days depending on location. Of particular interest to SMAP-based agricultural applications is a monitoring product that assesses the SMAP near-surface soil moisture in terms of probability percentiles for dry and wet conditions. However, the short SMAP record length poses a statistical challenge for meaningful assessment of its indices. This study presents initial insights about using SMAP for monitoring drought and pluvial regions with a first application over the contiguous United States (CONUS). SMAP soil moisture data from April 2015 to December 2017 at both near-surface (5 cm) SPL3SMP, or Level 3, at ∼36 km resolution, and root-zone SPL4SMAU, or Level 4, at ∼9 km resolution, were fitted to beta distributions and were used to construct probability distributions for warm (May–October) and cold (November–April) seasons. To assess the data adequacy and have confidence in using short-term SMAP for a drought index estimate, we analyzed individual grids by defining two filters and a combination of them, which could separate the 5815 grids covering CONUS into passed and failed grids. The two filters were (1) the Kolmogorov–Smirnov (KS) test for beta-fitted long-term and the short-term variable infiltration capacity (VIC) land surface model (LSM) with 95 % confidence and (2) good correlation (≥0.4) between beta-fitted VIC and beta-fitted SPL3SMP. To evaluate which filter is the best, we defined a mean distance (MD) metric, assuming a VIC index at 36 km resolution as the ground truth. For both warm and cold seasons, the union of the filters – which also gives the best coverage of the grids throughout CONUS – was chosen to be the most reliable filter. We visually compared our SMAP-based drought index maps with metrics such as the U.S. Drought Monitor (from D0–D4), 1-month Standard Precipitation Index (SPI) and near-surface VIC from Princeton University. The root-zone drought index maps were shown to be similar to those produced by the root-zone VIC, 3-month SPI, and the Gravity Recovery and Climate Experiment (GRACE). This study is a step forward towards building a national and international soil moisture monitoring system without which quantitative measures of drought and pluvial conditions will remain difficult to judge.


2018 ◽  
Author(s):  
Sara Sadri ◽  
Eric F. Wood ◽  
Ming Pan

Abstract. Since April 2015, NASA's Soil Moisture Active Passive (SMAP) mission has monitored near-surface soil moisture, mapping the globe between the latitude bands of 85.044° N/S in 2–3 days depending on location. SMAP Level 3 passive radiometer product (SPL3SMP) measures the amount of water in the top 5 cm of soil except for regions of heavy vegetation (vegetation water content >4.5 kg/m2) and frozen or snow covered locations. SPL3SMP retrievals are spatially and temporally discontinuous, so the 33 months offers a short SMAP record length and poses a statistical challenge for meaningful assessment of its indices. The SMAP SPL4SMAU data product provides global surface and root zone soil moisture at 9-km resolution based on assimilating the SPL3SMP product into the NASA Catchment land surface model. Of particular interest to SMAP-based agricultural applications is a monitoring product that assesses the SMAP near-surface soil moisture in terms of probability percentiles for dry and wet conditions. We describe here SMAP-based indices over the continental United States (CONUS) based on both near-surface and root zone soil moisture percentiles. The percentiles are based on fitting a Beta distribution to the retrieved moisture values. To assess the data adequacy, a statistical comparison is made between fitting the distribution to VIC soil moisture values for the days when SPL3SMP are available, versus fitting to a 1979–2017 VIC data record. For the cold season (November–April), 57 % of grids were deemed to be consistent between the periods, and 68 % in the warm season (May–October), based on a Kolmogorov–Smirnov statistical test. It is assumed that if grids passed the consistency test using VIC data, then the grid had sufficient SMAP data. Our near-surface and root zone drought index on maps are shown to be similar to those produced by the U.S. Drought Monitor (from D0-D4) and GRACE. In a similar manner, we extend the index to include pluvial conditions using indices W0-W4. This study is a step forward towards building a national and international soil moisture monitoring system, without which, quantitative measures of drought and pluvial conditions will remain difficult to judge.


2015 ◽  
Vol 12 (8) ◽  
pp. 7971-8004 ◽  
Author(s):  
C. Draper ◽  
R. Reichle

Abstract. Nine years 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 one 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 one year of data is available, the associated uncertainty in the rescaling parameters should not greatly reduce the average benefit gained from data assimilation, but locally and in extreme years there is a risk of increased errors.


2016 ◽  
Author(s):  
Gabriëlle J. M. De Lannoy ◽  
Rolf H. Reichle

Abstract. Three different data products from the Soil Moisture Ocean Salinity (SMOS) mission are assimilated separately into the Goddard Earth Observing System Model, version 5 (GEOS-5) to improve estimates of surface and root-zone soil moisture. The first product consists of multi-angle, dual-polarization brightness temperature (Tb) observations at the bottom of the atmosphere extracted from Level 1 data. The second product is a derived SMOS Tb product that mimics the data at 40° incidence angle from the Soil Moisture Active Passive mission. The third product is the operational SMOS Level 2 surface soil moisture (SM) retrieval product. The assimilation system uses a spatially distributed ensemble Kalman filter (EnKF) with seasonally varying climatological bias mitigation for Tb assimilation, whereas a time-invariant cumulative density function matching is used for SM retrieval assimilation. All assimilation experiments improve the soil moisture estimates compared to model-only simulations during the period 1 July 2010 to 1 May 2015 and for 187 sites across the United States. Especially in areas where the satellite data are most sensitive to surface soil moisture, large skill improvements (e.g. increase in anomaly correlation by 0.1) are found in the surface soil moisture. The domain-average surface and root-zone skill metrics are similar among the various assimilation experiments, but large differences in skill are found locally. The observation-minus-forecast residuals and analysis increments reveal large differences in how the observations add value in the Tb and SM retrieval assimilation systems. The distinct patterns of these diagnostics in the two systems reflect observation and model errors patterns that are not well captured in the assigned EnKF error parameters. Consequently, a localized optimization of the EnKF error parameters is needed to further improve Tb or SM retrieval assimilation.


Sign in / Sign up

Export Citation Format

Share Document