Improving Simulated Soil Moisture Fields Through Assimilation of AMSR-E Soil Moisture Retrievals with an Ensemble Kalman Filter and a Mass Conservation Constraint

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

Data Assimilation of Remote Sensing Soil Moisture Retrievals with Local Ensemble Transform Kalman Filter Scheme

Department of Urban and Environmental Engineering (Environmental Science and Engineering)This thesis proposes a development of land data assimilation system to produce realistic land surface states, which is performed with diverse remote sensing retrievals using advanced land surface model (LSM) and data assimilation techniques. Remote sensed soil moisture retrievals with high-temporal and ???spatial resolution is recently available. For instance, the Advanced Scatterometer (ASCAT), a C-band active microwave remote sensing instrument, and the NASA Soil Moisture Active Passive (SMAP) L-band passive remote sensing retrievals provide global near-surface soil moisture condition in real-time. Bias corrected observation datasets are used in the assimilation based on cumulative distribution function fitting because there is a large discrepancy of soil moisture contents between each retrieval and LSM offline simulation by difference sensed layer depth and algorithms, and characteristics of model physics. This study performs the soil moisture data assimilation using these bias corrected satellite retrievals with Local Ensemble Transform Kalman Filter (LETKF) scheme. The impact of the soil moisture assimilation is evaluated with ground based in-situ soil moisture measurement network over the globe. The result reveals that each satellite retrieval provides significant added value in the data assimilation. The impact of the assimilation tends to be better improved when active and passive satellite retrievals are simultaneously used. The temporal correlation of the assimilated soil moisture increases surface soil moisture skills by ??R???0.12 over the continental U.S., and the improvement at root-zone is ??R???0.1. The result is explained by ???Assimilation Gain???, where the quality of assimilated satellite data and the number of assimilated observations strongly contribute to the skill improvement of assimilated soil moisture estimates. The skill improvement through the multi-retrieval assimilation is mostly significant in transitional climate regime where land-atmosphere interaction is strong, and the impact of soil moisture initialization is clearly shown in the forecast model. Furthermore, the skill improvement is also significant in other validated regions (e.g. western Europe, and central Tibetan Plateau). The magnitude of the skill improvement through the assimilation is large when the quality of satellite retrievals tends to be better than that of the open loop. The assimilated soil moisture estimates are widely used in understanding land surface physical process of hydrology cycle and the land-atmosphere interaction. Furthermore, the realistic land surface condition gives the better information in land surface monitoring system.clos

Assimilation of GRACE Terrestrial Water Storage into a Land Surface Model: Evaluation 1 and Potential Value for Drought Monitoring in Western and Central Europe

A land surface model s ability to simulate states (e.g., soil moisture) and fluxes (e.g., runoff) is limited by uncertainties in meteorological forcing and parameter inputs as well as inadequacies in model physics. In this study, anomalies of terrestrial water storage (TWS) observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission were assimilated into the NASA Catchment land surface model in western and central Europe for a 7-year period, using a previously developed ensemble Kalman smoother. GRACE data assimilation led to improved runoff correlations with gauge data in 17 out of 18 hydrological basins, even in basins smaller than the effective resolution of GRACE. Improvements in root zone soil moisture were less conclusive, partly due to the shortness of the in situ data record. In addition to improving temporal correlations, GRACE data assimilation also reduced increasing trends in simulated monthly TWS and runoff associated with increasing rates of precipitation. GRACE assimilated root zone soil moisture and TWS fields exhibited significant changes in their dryness rankings relative to those without data assimilation, suggesting that GRACE data assimilation could have a substantial impact on drought monitoring. Signals of drought in GRACE TWS correlated well with MODIS Normalized Difference Vegetation Index (NDVI) data in most areas. Although they detected the same droughts during warm seasons, drought signatures in GRACE derived TWS exhibited greater persistence than those in NDVI throughout all seasons, in part due to limitations associated with the seasonality of vegetation

Identifying and Addressing Land Surface Model Deficiencies with Data Assimilation

Land surface models (LSMs) encapsulate our understanding of terrestrial water and energy cycle physics and provide estimates of land surface states and fluxes when and where measurement gaps exist. Gaps in our understanding of the physics are a different issue. Data assimilation can address that issue both directly, through updating of prognostic model variables, or indirectly, when the simulated world conflicts with observation, necessitating adjustment of the model. Here we will focus on the latter case and present several examples, including (1) depth to bedrock adjustment to accommodate assimilated GRACE terrestrial water storage data; (2) steps to prevent immediate melting of assimilated snow cover; (3) irrigation's contribution to evapotranspiration; (4) lessons learned from soil moisture data assimilation; (5) the potential impact of satellite based runoff observatio

Assimilation of GRACE terrestrial water storage into a land surface model: Evaluation and potential value for drought monitoring in western and central Europe

A land surface model’s ability to simulate states (e.g., soil moisture) and fluxes (e.g., runoff) is limited by uncertainties in meteorological forcing and parameter inputs as well as inadequacies in model physics. In this study, anomalies of terrestrial water storage (TWS) observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission were assimilated into the NASA Catchment land surface model in western and central Europe for a 7-year period, using a previously developed ensemble Kalman smoother. GRACE data assimilation led to improved runoff estimates (in temporal correlation and root mean square error) in 17 out of 18 hydrological basins, even in basins smaller than the effective resolution of GRACE. Improvements in root zone soil moisture were less conclusive, partly due to the shortness of the in situ data record. GRACE data assimilation also had significant impacts in groundwater estimates including trend and seasonality. In addition to improving temporal correlations, GRACE data assimilation also reduced increasing trends in simulated monthly TWS and runoff associated with increasing rates of precipitation. The assimilation downscaled (in space and time) and disaggregated GRACE data into finer scale components of TWS which exhibited significant changes in their dryness rankings relative to those without data assimilation, suggesting that GRACE data assimilation could have a substantial impact on drought monitoring. Signals of drought in GRACE TWS correlated well with MODIS Normalized Difference Vegetation Index (NDVI) data in most areas. Although they detected the same droughts during warm seasons, drought signatures in GRACE derived TWS exhibited greater persistence than those in NDVI throughout all seasons, in part due to limitations associated with the seasonality of vegetation. Mass imbalances associated with GRACE data assimilation and challenges of using GRACE data for drought monitoring are discussed

To improve model soil moisture estimation in arid/semi-arid region using in situ and remote sensing information

Soil moisture plays a key role in water and energy exchange in the land hydrologic process. Effective soil moisture information can be used for many applications in weather and hydrological forecasting, water resources, and irrigation system management and planning. However, to accurate modeling of soil moisture variation in the soil layer is still very challenging. In this study, in situ and remote sensing information of near-surface soil moisture is assimilated into the Noah land surface model (LSM) to estimate deep-layer soil moisture variation. The sequential Monte Carlo-Particle Filter technique, being well known for capability of modeling high nonlinear and non-Gaussian processes, is applied to assimilate surface soil moisture measurement to the deep layers. The experiments were carried out over several locations over the semi-arid region of the US. Comparing with in situ observations, the assimilation runs show much improved from the control (non-assimilation) runs for estimating both soil moisture and temperature at 5-, 20-, and 50-cm soil depths in the Noah LSM. © 2012 Springer-Verlag

Evaluation of the Land Surface Water Budget in NCEP/NCAR and NCEP/DOE Reanalyses using an Off-line Hydrologic Model

The ability of the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) reanalysis (NRA1) and the follow-up NCEP/Department of Energy (DOE) reanalysis (NRA2), to reproduce the hydrologic budgets over the Mississippi River basin is evaluated using a macroscale hydrology model. This diagnosis is aided by a relatively unconstrained global climate simulation using the NCEP global spectral model, and a more highly constrained regional climate simulation using the NCEP regional spectral model, both employing the same land surface parameterization (LSP) as the reanalyses. The hydrology model is the variable infiltration capacity (VIC) model, which is forced by gridded observed precipitation and temperature. It reproduces observed streamflow, and by closure is constrained to balance other terms in the surface water and energy budgets. The VIC-simulated surface fluxes therefore provide a benchmark for evaluating the predictions from the reanalyses and the climate models. The comparisons, conducted for the 10-year period 1988–1997, show the well-known overestimation of summer precipitation in the southeastern Mississippi River basin, a consistent overestimation of evapotranspiration, and an underprediction of snow in NRA1. These biases are generally lower in NRA2, though a large overprediction of snow water equivalent exists. NRA1 is subject to errors in the surface water budget due to nudging of modeled soil moisture to an assumed climatology. The nudging and precipitation bias alone do not explain the consistent overprediction of evapotranspiration throughout the basin. Another source of error is the gravitational drainage term in the NCEP LSP, which produces the majority of the model\u27s reported runoff. This may contribute to an overprediction of persistence of surface water anomalies in much of the basin. Residual evapotranspiration inferred from an atmospheric balance of NRA1, which is more directly related to observed atmospheric variables, matches the VIC prediction much more closely than the coupled models. However, the persistence of the residual evapotranspiration is much less than is predicted by the hydrological model or the climate models

Assessment and enhancement of MERRA land surface hydrology estimates

The Modern-Era Retrospective Analysis for Research and Applications (MERRA) is a state-of-the-art reanalysis that provides, in addition to atmospheric fields, global estimates of soil moisture, latent heat flux, snow, and runoff for 1979 present. This study introduces a supplemental and improved set of land surface hydrological fields ("MERRA-Land") generated by rerunning a revised version of the land component of the MERRA system. Specifically, the MERRA-Land estimates benefit from corrections to the precipitation forcing with the Global Precipitation Climatology Project pentad product (version 2.1) and from revised parameter values in the rainfall interception model, changes that effectively correct for known limitations in the MERRA surface meteorological forcings. The skill (defined as the correlation coefficient of the anomaly time series) in land surface hydrological fields from MERRA and MERRA-Land is assessed here against observations and compared to the skill of the state-of-the-art ECMWF Re-Analysis-Interim (ERA-I). MERRA-Land and ERA-I root zone soil moisture skills (against in situ observations at 85 U.S. stations) are comparable and significantly greater than that of MERRA. Throughout the Northern Hemisphere, MERRA and MERRA-Land agree reasonably well with in situ snow depth measurements (from 583 stations) and with snow water equivalent from an independent analysis. Runoff skill (against naturalized stream flow observations from 18 U.S. basins) of MERRA and MERRA-Land is typically higher than that of ERA-I. With a few exceptions, the MERRA-Land data appear more accurate than the original MERRA estimates and are thus recommended for those interested in using MERRA output for land surface hydrological studies

Assessment of assimilating SMOS soil moisture information into a distributed hydrologic model

The role that soil moisture plays in terms of modulating hydrologic processes including infiltration and runoff generation makes it an essential component to capture for hydrologic modeling. This work aims to leverage information gained from SMOS to improve surface soil moisture simulations in the Russian River Basin (California, U.S.A). The basin's complex terrain offers a rigorous testing ground for SMOS soil moisture products. Data from seven in situ observation sites are used to assess model performance after assimilating SMOS-based soil saturation ratios. For a comparison of "best case" scenarios, the in situ observations themselves are assimilated. Results show that SMOS assimilated simulations shows modest improvement at most in situ locations. Despite the observed decrease in model performance at some locations, overall performance of simulations assimilated with SMOS-based saturation ratios remains high. Findings suggest that even in a complex environment, useful information may be extracted from SMOS estimates for hydrologic modeling