Reconciling Flagging Strategies for Multi-Sensor Satellite Soil Moisture Climate Data Records

Reliable soil moisture retrievals from passive microwave satellite sensors are limited during certain conditions, e.g., snow coverage, radio-frequency interference, and dense vegetation. In these cases, the retrievals can be masked using flagging algorithms. Currently available single- and multi-sensor soil moisture products utilize different flagging approaches. However, a clear overview and comparison of these approaches and their impact on soil moisture data are still lacking. For long-term climate records such as the soil moisture products of the European Space Agency (ESA) Climate Change Initiative (CCI), the effect of any flagging inconsistency resulting from combining multiple sensor datasets is not yet understood. Therefore, the first objective of this study is to review the data flagging system that is used within multi-sensor ESA CCI soil moisture products as well as the flagging systems of two other soil moisture datasets from sensors that are also used for the ESA CCI soil moisture products: The level 3 Soil Moisture and Ocean Salinity (SMOS) and the Soil Moisture Active/Passive (SMAP). The SMOS and SMAP soil moisture flagging systems differ substantially in number and type of conditions considered, critical flags, and data source dependencies. The impact on the data availability of the different flagging systems were compared for the SMOS and SMAP soil moisture datasets. Major differences in data availability were observed globally, especially for northern high latitudes, mountainous regions, and equatorial latitudes (up to 37%, 33%, and 32% respectively) with large seasonal variability. These results highlight the importance of a consistent and well-performing approach that is applicable to all individual products used in long-term soil moisture data records. Consequently, the second objective of the present study is to design a consistent and model-independent flagging strategy to improve soil moisture climate records such as the ESA CCI products. As snow cover, ice, and frozen conditions were demonstrated to have the biggest impact on data availability, a uniform satellite driven flagging strategy was designed for these conditions and evaluated against two ground observation networks. The new flagging strategy demonstrated to be a robust flagging alternative when compared to the individual flagging strategies adopted by the SMOS and SMAP soil moisture datasets with a similar performance, but with the applicability to the entire ESA CCI time record without the use of modelled approximations.

Evaluation of a Microwave Emissivity Module for Snow Covered Area with CMEM in the ECMWF Integrated Forecasting System

The Community Microwave Emission Modelling platform (CMEM) has been developed by the European Centre for Medium-Range Weather Forecasts (ECMWF) as the forward operator radiative transfer model for low frequency passive microwave brightness temperatures (TB). It is used at ECMWF for L-band TB monitoring over snow free areas. In this paper, upgrades to CMEM are presented in order to explore forward modelling in snow-covered areas for coupled land-atmosphere numerical weather prediction systems. The upgrades enable to use CMEM on an extended range of frequencies and the Helsinki University of Technology multi-layer snow emission model is implemented. Offline CMEM experiments are evaluated against AMSR2 (Advanced Microwave Scanning Radiometer 2) observations showing that simulated TB is improved when using a multi-layer snow scheme, compared to a single-layer scheme. The improvements mainly result from a better representation of snow characteristics in the multi-layer snowpack model. CMEM is also evaluated in the Integrated Forecasting System and coupled to RTTOV (Radiative Transfer for TOVS). The numerical results show improved simulated TB at low frequency V polarization over snow-covered area compared to a configuration using emissivity atlas. Degradations at frequencies higher than 20 GHz indicate that further improvements are required in the emissivity and snowpack properties modelling.

Mapping of Snow Depth by Blending Satellite and In-Situ Data Using Two-Dimensional Optimal Interpolation—Application to AMSR2

The development of a snow depth product over North America is investigated by applying two-dimensional optimal interpolation to passive microwave satellite-derived and in-situ measured snow depth. At each snow-covered satellite footprint, the technique computes a snow depth increment as the weighted average of data increments, and updates the satellite-derived snow depth accordingly. Data increments are computed as the difference between the in-situ-measured and satellite snow depth at station locations surrounding the satellite footprint. Calculation of optimal weights is based on spatial lag autocorrelation of snow depth increments, modelled as functions of horizontal distance and elevation difference between pairs of observations. The technique is applied to Advanced Microwave Scanning Radiometer 2 (AMSR2) snow depth and in-situ snow depth obtained from the Global Historical Climatology Network. The results over North America during January–February 2017 indicate that the technique greatly enhances the performance of the satellite estimates, especially over mountain terrain, albeit with an accuracy inferior to that over low elevation areas. Moreover, the technique generates more accurate output compared to that from NOAA’s Global Forecast System, with implications for improving the utilization of satellite data in snow assessments and numerical weather prediction.

A Prognostic Nested k-Nearest Approach for Microwave Precipitation Phase Detection over Snow Cover

Monitoring changes of precipitation phase from space is important for understanding the mass balance of Earth's cryosphere in a changing climate. This paper examines a Bayesian nearest neighbor approach for prognostic detection of precipitation and its phase using passive microwave observations from the Global Precipitation Measurement (GPM) satellite. The method uses the weighted Euclidean distance metric to search through an a priori database populated with coincident GPM radiometer and radar observations as well as ancillary snow-cover data. The algorithm performance is evaluated using data from GPM official precipitation products, ground-based radars, and high-fidelity simulations from the Weather Research and Forecasting Model. Using the presented approach, we demonstrate that the hit probability of terrestrial precipitation detection can reach to 0.80, while the probability of false alarm remains below 0.11. The algorithm demonstrates higher skill in detecting snowfall than rainfall, on average by 10%. In particular, the probability of precipitation detection and its solid phase increases by 11% and 8%, over dry snow cover, when compared to other surface types. The main reason is found to be related to the ability of the algorithm in capturing the signal of increased liquid water content in snowy clouds over radiometrically cold snow-covered surfaces.

Comparison of commonly-used microwave radiative transfer models for snow remote sensing

This paper reviews four commonly-used microwave radiative transfer models that take different electromagnetic approaches to simulate snow brightness temperature (T_B): the Dense Media Radiative Transfer - Multi-Layer model (DMRT-ML), the Dense Media Radiative Transfer - Quasi-Crystalline Approximation Mie scattering of Sticky spheres (DMRT-QMS), the Helsinki University of Technology n-Layers model (HUT-nlayers) and the Microwave Emission Model of Layered Snowpacks (MEMLS). Using the same extensively measured physical snowpack properties, we compared the simulated T_B at 11, 19 and 37 GHz from these four models. The analysis focuses on the impact of using different types of measured snow microstructure metrics in the simulations. In addition to density, snow microstructure is defined for each snow layer by grain optical diameter (Do) and stickiness for DMRT-ML and DMRT-QMS, mean grain geometrical maximum extent (D_max) for HUT n-layers and the exponential correlation length for MEMLS. These metrics were derived from either in-situ measurements of snow specific surface area (SSA) or macrophotos of grain sizes (D_max), assuming non-sticky spheres for the DMRT models. Simulated T_B sensitivity analysis using the same inputs shows relatively consistent T_B behavior as a function of Do and density variations for the vertical polarization (maximum deviation of 18 K and 27 K, respectively), while some divergences appear in simulated variations for the polarization ratio (PR). Comparisons with ground-based radiometric measurements show that the simulations based on snow SSA measurements have to be scaled with a model-specific factor of Do in order to minimize the root mean square error (RMSE) between measured and simulated T_B. Results using in-situ grain size measurements (SSA or D_max, depending on the model) give a mean T_B RMSE (19 and 37 GHz) of the order of 16-26 K, which is similar for all models when the snow microstructure metrics are scaled. However, the MEMLS model converges to better results when driven by the correlation length estimated from in-situ SSA measurements rather than D_max measurements. On a practical level, this paper shows that the SSA parameter, a snow property that is easy to retrieve in-situ, appears to be the most relevant parameter for characterizing snow microstructure, despite the need for a scaling factor.