Assimilation of Terrestrial Water Storage from GRACE in a Snow-Dominated Basin

Terrestrial water storage (TWS) information derived from Gravity Recovery and Climate Experiment (GRACE) measurements is assimilated into a land surface model over the Mackenzie River basin located in northwest Canada. Assimilation is conducted using an ensemble Kalman smoother (EnKS). Model estimates with and without assimilation are compared against independent observational data sets of snow water equivalent (SWE) and runoff. For SWE, modest improvements in mean difference (MD) and root mean squared difference (RMSD) are achieved as a result of the assimilation. No significant differences in temporal correlations of SWE resulted. Runoff statistics of MD remain relatively unchanged while RMSD statistics, in general, are improved in most of the sub-basins. Temporal correlations are degraded within the most upstream sub-basin, but are, in general, improved at the downstream locations, which are more representative of an integrated basin response. GRACE assimilation using an EnKS offers improvements in hydrologic state/flux estimation, though comparisons with observed runoff would be enhanced by the use of river routing and lake storage routines within the prognostic land surface model. Further, GRACE hydrology products would benefit from the inclusion of better constrained models of post-glacial rebound, which significantly affects GRACE estimates of interannual hydrologic variability in the Mackenzie River basin

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems

Assimilation of Passive and Active Microwave Soil Moisture Retrievals

Root-zone soil moisture is an important control over the partition of land surface energy and moisture, and the assimilation of remotely sensed near-surface soil moisture has been shown to improve model profile soil moisture [1]. To date, efforts to assimilate remotely sensed near-surface soil moisture at large scales have focused on soil moisture derived from the passive microwave Advanced Microwave Scanning Radiometer (AMSR-E) and the active Advanced Scatterometer (ASCAT; together with its predecessor on the European Remote Sensing satellites (ERS. The assimilation of passive and active microwave soil moisture observations has not yet been directly compared, and so this study compares the impact of assimilating ASCAT and AMSR-E soil moisture data, both separately and together. Since the soil moisture retrieval skill from active and passive microwave data is thought to differ according to surface characteristics [2], the impact of each assimilation on the model soil moisture skill is assessed according to land cover type, by comparison to in situ soil moisture observations

L-band Microwave Remote Sensing and Land Data Assimilation Improve the Representation of Prestorm Soil Moisture Conditions for Hydrologic Forecasting

Recent advances in remote sensing and land data assimilation purport to improve the quality of antecedent soil moisture information available for operational hydrologic forecasting. We objectively validate this claim by calculating the strength of the relationship between storm-scale runoff ratio (i.e., total stream flow divided by total rainfall accumulation in depth units) and pre-storm surface soil moisture estimates from a range of surface soil moisture data products. Results demonstrate that both satellite-based, L-band microwave radiometry and the application of land data assimilation techniques have significantly improved the utility of surface soil moisture data sets for forecasting stream flow response to future rainfall events

Exploiting Soil Moisture, Precipitation and Streamflow Observations to Evaluate Soil Moisture/Runoff Coupling in Land Surface Models

Accurate partitioning of precipitation into infiltration and runoff is a fundamental objective of land surface models tasked with characterizing the surface water and energy balance. Temporal variability in this partitioning is due, in part, to changes in prestorm soil moisture, which determine soil infiltration capacity and unsaturated storage. Utilizing the National Aeronautics and Space Administration Soil Moisture Active Passive Level4 soil moisture product in combination with streamflow and precipitation observations, we demonstrate that land surface models (LSMs) generally underestimate the strength of the positive rank correlation between prestorm soil moisture and event runoff coefficients (i.e., the fraction of rainfall accumulation volume converted into stormflow runoff during a storm event). Underestimation is largest for LSMs employing an infiltrationexcess approach for stormflow runoff generation. More accurate coupling strength is found in LSMs that explicitly represent subsurface stormflow or saturationexcess runoff generation processes

Plant Water Uptake Thresholds Inferred From Satellite Soil Moisture

Empirical functions are widely used in hydrological, agricultural, and Earth system models to parameterize plant water uptake. We infer soil water potentials at which uptake is downregulated from its well‐watered rate and at which uptake ceases, in biomes with <60% woody vegetation at 36‐km grid resolution. We estimate thresholds through Bayesian inference using a stochastic soil water balance framework to construct theoretical soil moisture probability distributions consistent with empirical distributions derived from satellite soil moisture observations. The global median Nash–Sutcliffe efficiency between empirical soil moisture distributions and theoretical distributions using reference constants, inferred median parameters per biome, and spatially variable inferred parameters are 0.38, 0.59, and 0.8, respectively. Spatially variable thresholds capture location‐specific vegetation and climate characteristics and can be connected to biome‐level water uptake strategies. Results demonstrate that satellite soil moisture probability distributions encode information, valuable to understanding biome‐level ecohydrological adaptation and resistance to climate variability

The SMAP Level 4 Carbon PRODUCT for Monitoring Terrestrial Ecosystem-Atmosphere CO2 Exchange

The NASA Soil Moisture Active Passive (SMAP) mission Level 4 Carbon (L4_C) product provides model estimates of Net Ecosystem CO2 exchange (NEE) incorporating SMAP soil moisture information as a primary driver. The L4_C product provides NEE, computed as total respiration less gross photosynthesis, at a daily time step and approximate 14-day latency posted to a 9-km global grid summarized by plant functional type. The L4_C product includes component carbon fluxes, surface soil organic carbon stocks, underlying environmental constraints, and detailed uncertainty metrics. The L4_C model is driven by the SMAP Level 4 Soil Moisture (L4_SM) data assimilation product, with additional inputs from the Goddard Earth Observing System, Version 5 (GEOS-5) weather analysis and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. The L4_C data record extends from March 2015 to present with ongoing production. Initial comparisons against global CO2 eddy flux tower measurements, satellite Solar Induced Canopy Florescence (SIF) and other independent observation benchmarks show favorable L4_C performance and accuracy, capturing the dynamic biosphere response to recent weather anomalies and demonstrating the value of SMAP observations for monitoring of global terrestrial water and carbon cycle linkages

The International Mass Loading Service

The International Mass Loading Service computes four loadings: a) atmospheric pressure loading; b) land water storage loading; c) oceanic tidal loading; and d) non-tidal oceanic loading. The service provides to users the mass loading time series in three forms: 1) pre-computed time series for a list of 849 space geodesy stations; 2) pre-computed time series on the global 1deg x 1deg grid; and 3) on-demand Internet service for a list of stations and a time range specified by the user. The loading displacements are provided for the time period from 1979.01.01 through present, updated on an hourly basis, and have latencies 8-20 hours.Comment: 8 pages, 3 figures, to appear in the Proceedings of the Reference Frames for Applications in Geosciences Simposium, held in Luxemboug in October 201