Snow depth variability in the Northern Hemisphere mountains observed from space.

Accurate snow depth observations are critical to assess water resources. More than a billion people rely on water from snow, most of which originates in the Northern Hemisphere mountain ranges. Yet, remote sensing observations of mountain snow depth are still lacking at the large scale. Here, we show the ability of Sentinel-1 to map snow depth in the Northern Hemisphere mountains at 1 km² resolution using an empirical change detection approach. An evaluation with measurements from ~4000 sites and reanalysis data demonstrates that the Sentinel-1 retrievals capture the spatial variability between and within mountain ranges, as well as their inter-annual differences. This is showcased with the contrasting snow depths between 2017 and 2018 in the US Sierra Nevada and European Alps. With Sentinel-1 continuity ensured until 2030 and likely beyond, these findings lay a foundation for quantifying the long-term vulnerability of mountain snow-water resources to climate change.

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach using SMOS L-Band Microwave Data.

The uncertainty of surface soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several established and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters. The second part of this paper evaluates the ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM, and the relationship with time series metrics obtained by comparing SM retrievals with in situ data. As expected for bounded variables, the total spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term SM bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread shows higher values when SM deviates further from its long-term mean SM. This reveals higher-order biases (e.g. differences in variances) in the retrieval error, which should be considered when characterizing retrieval error. The time-average anomaly ensemble SM spread of 0.037 m3/m3 approximates the actual time series unbiased root-mean-square-difference of 0.042 m3/m3 between ensemble mean retrievals and in situ data across the reference sites.

Assimilation of MODIS Snow Cover Fraction Observations into the NASA Catchment Land Surface Model.

The NASA Catchment land surface model (CLSM) is the land model component used for the Modern-Era Retrospective Analysis for Research and Applications (MERRA). Here, the CLSM versions of MERRA and MERRA-Land are evaluated using snow cover fraction (SCF) observations from the Moderate Resolution Imaging Spectroradiometer (MODIS). Moreover, a computationally-efficient empirical scheme is designed to improve CLSM estimates of SCF, snow depth, and snow water equivalent (SWE) through the assimilation of MODIS SCF observations. Results show that data assimilation (DA) improved SCF estimates compared to the open-loop model without assimilation (OL), especially in areas with ephemeral snow cover and mountainous regions. A comparison of the SCF estimates from DA against snow cover estimates from the NOAA Interactive Multisensor Snow and Ice Mapping System showed an improvement in the probability of detection of up to 28% and a reduction in false alarms by up to 6% (relative to OL). A comparison of the model snow depth estimates against Canadian Meteorological Centre analyses showed that DA successfully improved the model seasonal bias from -0.017 m for OL to -0.007 m for DA, although there was no significant change in root-mean-square differences (RMSD) (0.095 m for OL, 0.093 m for DA). The time-average of the spatial correlation coefficient also improved from 0.61 for OL to 0.63 for DA. A comparison against in situ SWE measurements also showed improvements from assimilation. The correlation increased from 0.44 for OL to 0.49 for DA, the bias improved from -0.111 m for OL to -0.100 m for DA, and the RMSD decreased from 0.186 m for OL to 0.180 m for DA.

Benefits and Pitfalls of GRACE Data Assimilation: a Case Study of Terrestrial Water Storage Depletion in India.

This study investigates some of the benefits and drawbacks of assimilating Terrestrial Water Storage (TWS) observations from the Gravity Recovery and Climate Experiment (GRACE) into a land surface model over India. GRACE observes TWS depletion associated with anthropogenic groundwater extraction in northwest India. The model, however, does not represent anthropogenic groundwater withdrawals and is not skillful in reproducing the interannual variability of groundwater. Assimilation of GRACE TWS introduces long-term trends and improves the interannual variability in groundwater. But the assimilation also introduces a negative trend in simulated evapotranspiration whereas in reality evapotranspiration is likely enhanced by irrigation, which is also unmodeled. Moreover, in situ measurements of shallow groundwater show no trend, suggesting that the trends are erroneously introduced by the assimilation into the modeled shallow groundwater, when in reality the groundwater is depleted in deeper aquifers. The results emphasize the importance of representing anthropogenic processes in land surface modeling and data assimilation systems.

Global Assessment of the SMAP Level-4 Surface and Root-Zone Soil Moisture Product Using Assimilation Diagnostics.

The Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and related land surface variables from 31 March 2015 to present with ~2.5day latency. The ensemble-based L4_SM algorithm assimilates SMAP brightness temperature (Tb) observations into the Catchment land surface model. This study describes the spatially distributed L4_SM analysis and assesses the observation-minus-forecast (O-F) Tb residuals and the soil moisture and temperature analysis increments. Owing to the climatological rescaling of the Tb observations prior to assimilation, the analysis is essentially unbiased, with global mean values of ~0.37 K for the O-F Tb residuals and practically zero for the soil moisture and temperature increments. There are, however, modest regional (absolute) biases in the O-F residuals (under ~3 K), the soil moisture increments (under ~0.01 m3 m-3), and the surface soil temperature increments (under ~1 K). Typical instantaneous values are ~6 K for O-F residuals, ~0.01 (~0.003) m3 m-3 for surface (root-zone) soil moisture increments, and ~0.6 K for surface soil temperature increments. The O-F diagnostics indicate that the actual errors in the system are overestimated in deserts and densely vegetated regions and underestimated in agricultural regions and transition zones between dry and wet climates. The O-F auto-correlations suggest that the SMAP observations are used efficiently in western North America, the Sahel, and Australia, but not in many forested regions and the high northern latitudes. A case study in Australia demonstrates that assimilating SMAP observations successfully corrects short-term errors in the L4_SM rainfall forcing.

Precipitation Estimation Using L-Band and C-Band Soil Moisture Retrievals.

An established methodology for estimating precipitation amounts from satellite-based soil moisture retrievals is applied to L-band products from the Soil Moisture Active Passive (SMAP) and Soil Moisture and Ocean Salinity (SMOS) satellite missions and to a C-band product from the Advanced Scatterometer (ASCAT) mission. The precipitation estimates so obtained are evaluated against in situ (gauge-based) precipitation observations from across the globe. The precipitation estimation skill achieved using the L-band SMAP and SMOS datasets is higher than that obtained with the C-band product, as might be expected given that L-band is sensitive to a thicker layer of soil and thereby provides more information on the response of soil moisture to precipitation. The square of the correlation coefficient between the SMAP-based precipitation estimates and the observations (for aggregations to ~100 km and 5 days) is on average about 0.6 in areas of high rain gauge density. Satellite missions specifically designed to monitor soil moisture thus do provide significant information on precipitation variability, information that could contribute to efforts in global precipitation estimation.