Estimating the landscape and soil freeze-thaw (FT) dynamics in the Northern Hemisphere (NH) is crucial for understanding permafrost response to global warming and changes in regional and global carbon budgets. A new framework for surface FT-cycle retrievals using L-band microwave radiometry based on a deep convolutional autoencoder neural network is presented. This framework defines the landscape FT-cycle retrieval as a time-series anomaly detection problem, considering the frozen states as normal and the thawed states as anomalies. The autoencoder retrieves the FT-cycle probabilistically through supervised reconstruction of the brightness temperature (TB) time series using a contrastive loss function that minimizes (maximizes) the reconstruction error for the peak winter (summer). Using the data provided by the Soil Moisture Active Passive (SMAP) satellite, it is demonstrated that the framework learns to isolate the landscape FT states over different land surface types with varying complexities related to the radiometric characteristics of snow cover, lake-ice phenology, and vegetation canopy. The consistency of the retrievals is assessed over Alaska using in situ observations, demonstrating an 11% improvement in accuracy and reduced uncertainties compared to traditional methods that rely on thresholding the normalized polarization ratio (NPR).
The freezing front depth (z(ff)) of annual freeze-thaw cycles is critical for monitoring the dynamics of the cryosphere under climate change because z(ff) is a sensitive indicator of the heat balance over the atmosphere-cryosphere interface. Meanwhile, although it is very promising for acquiring global soil moisture distribution, the L-band microwave remote sensing products over seasonally frozen grounds and permafrost is much less than in wet soil. This study develops an algorithm, i.e., the brightness temperature inferred freezing front (BT-FF) model, for retrieving the interannual z(ff) with the diurnal amplitude variation of L-band brightness temperature (?T-B) during the freezing period. The new algorithm assumes first, the daily-scale solar radiation heating/cooling effect causes the daily surface thawing depth (z(tf)) variation, which leads further to ?T-B; second, ?T-B can be captured by an L-band radiometer; third, z(tf) and z(ff) are negatively linear correlated and their relation can be quantified using the Stefan equation. In this study, the modeled soil temperature profiles from the land surface model (STEMMUS-FT, i.e., simultaneous transfer of energy, mass, and momentum in unsaturated soil with freeze and thaw) and T-B observations from a tower-based L-band radiometer (ELBARA-III) at Maqu are used to validate the BT-FF model. It shows that, first, ?T-B can be precisely estimated from z(tf) during the daytime; second, the decreasing of z(tf) is linearly related to the increase of z(ff) with the Stefan equation; third, the accuracy of retrieved z(ff) is about 5-25 cm; fourth, the proposed model is applicable during the freezing period. The study is expected to extend the application of L-band T-B data in cryosphere/meteorology and construct global freezing depth dataset in the future.
The uncertainty of passive microwave retrievals of snowfall is notoriously high where high-frequency surface emissivity is significantly reduced and varies markedly in response to the changes in snowpack physical properties. Using the dense media radiative transfer theory, this article studies the potential effects of terrestrial snow-cover depth, density, and grain size on high-frequency channels 89 and 166 GHz of the radiometer onboard the Global Precipitation Measurement (GPM) core satellite, which are commonly used to capture snowfall scattering signals. Integrating the inference across all feasible grain sizes, ranges of snowpack density and depth are identified over which snowfall scattering signatures can be time-varying and potentially obscured. Using ten years of reanalysis data, the seasonal vulnerability of snowfall retrievals to the changes in snowpack emissivity in the Northern Hemisphere is mapped in a probabilistic sense and connections are made with the uncertainties of the GPM passive microwave snowfall retrievals. It is found that among different snow classes, relatively light Arctic tundra snow in fall, with a density below 260 kg m(-3), and shallow prairie snow during the winter, with a depth of less than 40 cm, can reduce the surface emissivity and obscure the snowfall passive microwave signatures. It is demonstrated that during winter, the highly vulnerable areas are over Kazakhstan and Mongolia with taiga and prairie snow. In the fall, these areas are largely over tundra and taiga snow in north of Russia and the Arctic Archipelagos as well as prairies in Canada and the Great Plains in the United States.
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.