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.

[1] Near-surface soil moisture observations from the active microwave ASCAT and the passive microwave AMSR-E satellite instruments are assimilated, both separately and together, into the NASA Catchment land surface model over 3.5 years using an ensemble Kalman filter. The impact of each assimilation is evaluated using in situ soil moisture observations from 85 sites in the US and Australia, in terms of the anomaly time series correlation-coefficient, R. The skill gained by assimilating either ASCAT or AMSR-E was very similar, even when separated by land cover type. Over all sites, the mean root-zone R was significantly increased from 0.45 for an open-loop, to 0.55, 0.54, and 0.56 by the assimilation of ASCAT, AMSR-E, and both, respectively. Each assimilation also had a positive impact over each land cover type sampled. For maximum accuracy and coverage it is recommended that active and passive microwave observations be assimilated together. Citation: Draper, C. S., R. H. Reichle, G. J. M. De Lannoy, and Q. Liu (2012), Assimilation of passive and active microwave soil moisture retrievals, Geophys.

The Soil Moisture Active Passive (SMAP) mission Level-4 Surface and Root-Zone Soil Moisture (L4_SM) data product is generated by assimilating SMAP L-band brightness temperature observations into the NASA Catchment land surface model. The L4_SM product is available from 31 March 2015 to present (within 3 days from real time) and provides 3-hourly, global, 9-km resolution estimates of surface (0–5 cm) and root-zone (0–100 cm) soil moisture and land surface conditions. This study presents an overview of the L4_SM algorithm, validation approach, and product assessment versus in situ measurements. Core validation sites provide spatially averaged surface (root zone) soil moisture measurements for 43 (17) “reference pixels” at 9- and 36-km gridcell scales located in 17 (7) distinct watersheds. Sparse networks provide point-scale measurements of surface (root zone) soil moisture at 406 (311) locations. Core validation site results indicate that the L4_SM product meets its soil moisture accuracy requirement, specified as an unbiased RMSE (ubRMSE, or standard deviation of the error) of 0.04 m3 m−3 or better. The ubRMSE for L4_SM surface (root zone) soil moisture is 0.038 m3 m−3 (0.030 m3 m−3) at the 9-km scale and 0.035 m3 m−3 (0.026 m3 m−3) at the 36-km scale. The L4_SM estimates improve (significantly at the 5% level for surface soil moisture) over model-only estimates, which do not benefit from the assimilation of SMAP brightness temperature observations and have a 9-km surface (root zone) ubRMSE of 0.042 m3 m−3 (0.032 m3 m−3). Time series correlations exhibit similar relative performance. The sparse network results corroborate these findings over a greater variety of climate and land cover conditions.

[1] Eight years (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010) of Advanced Microwave Scanning Radiometer-EOS (AMSR-E) snow water equivalent (SWE) retrievals and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) observations are assimilated separately or jointly into the Noah land surface model over a domain in Northern Colorado. A multiscale ensemble Kalman filter (EnKF) is used, supplemented with a rule-based update. The satellite data are either left unscaled or are scaled for anomaly assimilation. The results are validated against in situ observations at 14 high-elevation Snowpack Telemetry (SNOTEL) sites with typically deep snow and at 4 lower-elevation Cooperative Observer Program (COOP) sites. Assimilation of coarse-scale AMSR-E SWE and fine-scale MODIS SCF observations both result in realistic spatial SWE patterns. At COOP sites with shallow snowpacks, AMSR-E SWE and MODIS SCF data assimilation are beneficial separately, and joint SWE and SCF assimilation yields significantly improved root-mean-square error and correlation values for scaled and unscaled data assimilation. In areas of deep snow where the SNOTEL sites are located, however, AMSR-E retrievals are typically biased low and assimilation without prior scaling leads to degraded SWE estimates. Anomaly SWE assimilation could not improve the interannual SWE variations in the assimilation results because the AMSR-E retrievals lack realistic interannual variability in deep snowpacks. SCF assimilation has only a marginal impact at the SNOTEL locations because these sites experience extended periods of near-complete snow cover. Across all sites, SCF assimilation improves the timing of the onset of the snow season but without a net improvement of SWE amounts.

This study explores the benefits of assimilating SMOS soil moisture retrievals for hydrologic modeling, with a focus on soil moisture and streamflow simulations in the Murray Darling Basin, Australia. In this basin, floods occur relatively frequently and initial catchment storage is known to be key to runoff generation. The land surface model is the Variable Infiltration Capacity (VIC) model. The model is calibrated using the available streamflow records of 169 gauge stations across the Murray Darling Basin. The VIC soil moisture forecast is sequentially updated with observations from the SMOS Level 3 CATDS (Centre Aval de Traitement des Données SMOS) soil moisture product using the Ensemble Kalman filter. The assimilation algorithm accounts for the spatial mismatch between the model (0.125°) and the SMOS observation (25 km) grids. Three widely-used methods for removing bias between model simulations and satellite observations of soil moisture are evaluated. These methods match the first, second and higher order moments of the soil moisture distributions, respectively. In this study, the first order bias correction, i.e. the rescaling of the long term mean, is the recommended method. Preserving the observational variability of the SMOS soil moisture data leads to improved soil moisture updates, particularly for dry and wet conditions, and enhances initial conditions for runoff generation. Second or higher order bias correction, which includes a rescaling of the variance, decreases the temporal variability of the assimilation results. In comparison with in situ measurements of OzNet, the assimilation with mean bias correction reduces the root mean square error (RMSE) of the modeled soil moisture from 0.058 m 3 /m 3 to 0.046 m 3 /m 3 and increases the correlation from 0.564 to 0.714. These improvements in antecedent wetness conditions further translate into improved predictions of associated water fluxes, particularly runoff peaks. In conclusion, the results of this study clearly demonstrate the merit of SMOS data assimilation for soil moisture and streamflow predictions at the large scale.

Land surface models are usually biased in at least a subset of the simulated variables even after calibration. Bias estimation may therefore be needed for data assimilation. Here, in situ soil moisture profile observations in a small agricultural field were merged with Community Land Model (CLM2.0) simulations using different algorithms for state and forecast bias estimation with and without bias correction feedback. Simple state updating with the conventional ensemble Kalman filter (EnKF) allows for some implicit forecast bias correction. It is possible to estimate the soil moisture bias explicitly and derive superior soil moisture estimates with a generalized EnKF that uses a simple persistence model for the bias and assumes that the a priori bias error covariance is proportional to the a priori state error covariance. For the case of bi‐weekly assimilation of the entire profile of soil moisture observations, bias estimation and correction typically reduces the RMSE in soil moisture (over the standard EnKF without bias correction) by around 60 percent. However, under the above assumptions, significant improvements are limited to state variables for which observations are available. Therefore, it is crucial to measure the state variables of interest. The best variant for state and bias estimation depends on the nature of the model bias and the output of interest to the user. In a model that is only biased for soil moisture, large and frequent increments for soil moisture updating may be required, which in turn may negatively impact the water balance and output fluxes. It is then better to post‐process the soil moisture with the bias analysis without updating the model state.

No abstract

[1] The Community Land Model (CLM2.0) has been used to simulate land surface processes in a small corn field. The subdivision of grid cells into patches in the CLM2.0 was explored for the generation of Monte Carlo simulations for use in calibration and ensemble generation. A distributed multiobjective calibration was developed for the optimal estimation of parameters and initial state variables for 36 soil moisture profiles. Since the resulting parameter and initial state values did not lead to perfect simulations for soil moisture, and in order to better understand the forecast uncertainty, ensemble runs were generated. The ensembles generated by CLM2.0 have been verified by several methods that are commonly used in meteorology. It was shown that the perfect model approach cannot be applied for bounded hydrological applications and that perturbation of parameters is a necessity to obtain a realistic assessment of the forecast error. Perturbation of forcings only captures more of the model uncertainty than perturbation of initial conditions only, but also causes a too limited spread in the ensembles. The generation of ensemble members through perturbation of the parameter set, found through calibration, does not necessarily result in ensembles that surround the calibrated deterministic control run for soil moisture. This is partially due the nonlinearity of the model in the parameters. It may also indicate that some parameter sets are not robust and not appropriate to perturb for ensemble generation. Consequently, the resulting ensemble mean may not represent the best forecast or a priori state estimation. During periods of extreme drought or precipitation, the ensemble probability density function (pdf) deviates far from normality and the model behaves very nonlinearly. For state estimation, methods like the ensemble Kalman filter are best suited for the propagation of the first moments to account for the nonlinear dynamics during crucial events for hydrological simulations. However, the a posteriori estimate for this technique will only be optimal in the limited class of linear filters, since the underlying pdfs cannot be assumed to be Gaussian.