Passive Microwave Brightness Temperature Research Articles

Cross-calibration between MWRI and AMSR2 to improve consistency of snow depth products

The Microwave Radiation Imager (MWRI), boarded on the FY-3 series satellites: FY-3B, FY-3C, and FY-3D, is the first satellite-based microwave radiometer in China, commencing passive microwave brightness temperature data acquisition since 2010. The Advanced Microwave Scanning Radiometer 2 (AMSR2) boarded on the Global Change Observation Mission 1st-Water (GCOM-W1), has been operational since 2012. Despite the FY-3 series satellites are equipped with the same MWRI and all MWRIs sharing comparable parameters and configurations as AMSR2, disparities in observation times and satellite platforms result in inconsistencies in the data obtained by different satellites, which further impacting the consistency of retrieved geophysical parameters. To improve the consistency of brightness temperatures from FY-3B, FY-3C, FY-3D/MWRI, and GCOM-W1/AMSR2, cross-calibrations were conducted among brightness temperatures at ten-channel from above four platforms. The consistency of derived snow depth from MWRIs and AMSR2 in China before and after the calibration were also analyzed. The results show that the correlation coefficients of brightness temperatures at all channels between sensors exceed 0.98. After cross-calibration, the RMSEs and biases of brightness temperatures at all frequencies and snow depth in China derived from them reduce to varying degrees. The consistencies in both brightness temperatures and snow depth of FY-3B/MWRI, FY-3D/MWRI, and AMSR2 are higher than those of FY-3C and others. These findings advocate for the utilization of cross-calibrated brightness temperatures from FY-3B/MWRI, FY-3D/MWRI, and AMSR2, which share similar satellite overpass time, to derived a long-term snow depth dataset.

Snow Depth Estimation and Spatial and Temporal Variation Analysis in Tuha Region Based on Multi-Source Data

In the modelling of hydrological processes on a regional scale, remote-sensing snow depth products with a high spatial and temporal resolution are essential for climate change studies and for scientific decision-making by management. The existing snow depth products have low spatial resolution and are mostly applicable to large-scale studies; however, they are insufficiently accurate for the estimation of snow depth on a regional scale, especially in shallow snow areas and mountainous regions. In this study, we coupled SSM/I, SSMIS, and AMSR2 passive microwave brightness temperature data and MODIS, TM, and Landsat 8 OLI fractional snow cover area (fSCA) data, based on Python, with 30 m spatially resolved fractional snow cover area (fSCA) data obtained by the spatio-temporal dynamic warping algorithm to invert the low-resolution passive microwave snow depths, and we developed a spatially downscaled snow depth inversion method suitable for the Turpan–Hami region. However, due to the long data-processing time and the insufficient arithmetical power of the hardware, this study had to set the spatial resolution of the result output to 250 m. As a result, a day-by-day 250 m spatial resolution snow depth dataset for 20 hydrological years (1 August 2000–31 July 2020) was generated, and the accuracy was evaluated using the measured snow depth data from the meteorological stations, with the results of r = 0.836 (p ≤ 0.01), MAE = 1.496 cm, and RMSE = 2.597 cm, which are relatively reliable and more applicable to the Turpan–Hami area. Based on the spatially downscaled snow depth data produced, this study found that the snow in the Turpan–Hami area is mainly distributed in the northern part of Turpan (Bogda Mountain), the northwestern part of Hami (Barkun Autonomous Prefecture), and the central part of the area (North Tianshan Mountain, Barkun Mountain, and Harlik Mountain). The average annual snow depth in the Turpan–Hami area is only 0.89 cm, and the average annual snow depth increases with elevation, in line with the obvious law of vertical progression. The annual mean snow depth in the Turpan–Hami area showed a “fluctuating decreasing” trend with a rate of 0.01 cm·a−1 over the 20 hydrological years in the Turpan–Hami area. Overall, the spatially downscaled snow depth inversion algorithm developed in this study not only solves the problem of coarse spatial resolution of microwave brightness temperature data and the difficulty of obtaining accurate shallow snow depth but also solves the problem of estimating the shallow snow depth on a regional scale, which is of great significance for gaining a further understanding of the snow accumulation information in the Tuha region and for promoting the investigation and management of water resources in arid zones.

Hydrometeor Identification Using GMI Passive Microwave Brightness Temperatures

Abstract Using passive microwave brightness temperatures Tb from the Global Precipitation Measurement (GPM) Microwave Imager (GMI) and hydrometeor identification (HID) data from dual-polarization ground radars, empirical lookup tables are developed for a multifrequency estimation of the likelihood a precipitation column includes certain hydrometeor types, as a function of Tb. Eight years of collocated Tb and HID data from the GPM Validation Network are used for development and testing of the GMI-based HID retrieval, with 2015–20 used for training and 2021–22 used for testing the GMI-based HID retrieval. The occurrence of profiles with hail and graupel are both slightly underpredicted by the lookup tables, but the percentage of profiles predicted is highly correlated with the percentage observed (0.98 correlation coefficient for hail and 0.99 for graupel). By having snow appear before rain in the hierarchy, the sample size for rain, without ice aloft, is fairly small, and the percentage of rain profiles is less than snow for all Tb.

An Innovative Approach: Sea Ice Types Classification Using Convolutional Neural Networks with DDDTDWT Filter

Studying sea ice and its interaction with climate change is crucial due to its significant impact on the environment, society, and global stability. The pressing need to address the underlying reasons for the rapid melting of Arctic and Antarctic sea ice is underscored by its adverse effects on the environment and society. In this proposed study, a Convolutional Neural Network (CNN) is utilized to predict ice types using data from the NSIDC DAAC Advanced Microwave Scanning Radiometer - Earth Observing System Sensor (AMSR-E) collection. This dataset contains parameters such as sea ice types and spans data products from June 2002, obtained from the NASA Data Centre. By employing hand-crafted features as input and a single layer of hidden nodes, the CNN used in this approach generates improved estimates of ice types, outperforming traditional image analysis methods. At each stage, ConvNets use diverse filter banks, feature extraction pooling layers, and fully connected layers with basic activation functions like Relu. This allows the network to build multifaceted hierarchies of features. The sea ice type estimates produced by the CNN are then compared with those obtained from passive microwave brightness temperature data using existing algorithms as well as a proposed CNN algorithm, resulting in an increased classification accuracy of 98.58%. This improvement is particularly notable in the reduction of the error rate, which has been effectively minimized from 3.01% without feature selection to 1.42% with infinite feature selection. When compared to existing algorithms, the CNN demonstrates superior performance. These findings underscore the impact of input patch size, CNN layer count, and input size on the model's performance.

Estimation of Lightning Flash Rate in Precipitation Features by Applying Shallow AI Neural Network Models to the GMI Passive Microwave Brightness Temperatures

AbstractThe satellite‐constellation passive‐microwave Brightness Temperature (TB) observations, with global coverage, and more additions from upcoming CubeSats, have been mainly used in surface precipitation retrievals. However, these observations can also be used to indicate the intensity of convective systems. This study attempts to relate Global Precipitation Mission (GPM) Microwave Imager (GMI) TBs to Geostationary Lightning Mapper (GLM) lightning flashes by using 4 years (02/2018–04/2022) of GPM Precipitation Feature (PF) database. GMI TBs are collocated to the GLM lightning counts, and to the ERA5 reanalysis 2‐m air temperature in PFs. Three Artificial Intelligent Neural Network Models (AI‐NN) are trained to classify PFs producing lightning in a 20‐min window respectively for land, ocean, or coast. The flash rate of the determined Lightning producing PFs (LPFs) is then quantified by three other AI‐NNs, each trained for one of the three regions. Though the models clearly capture the global geographical distribution of LPFs with a Probability Of Detection over 90%, high False Alarm Rates are found, ranging from 49.9% over land to 91.5% over the ocean. The importance of TB at each passive microwave channel varies regionally, corresponding to the different microphysical properties in various types of precipitation systems. The global lightning distribution is derived by applying the AI models to global PFs and is well compared to the lightning climatology from Lightning Imaging Sensor and Optical Transient Detector. This suggests that the use of passive microwave TBs can help to fill the gaps in lightning monitoring thanks to their global coverage.

Uncertainty Calibration of Passive Microwave Brightness Temperatures Predicted by Bayesian Deep Learning Models

Abstract Visible and infrared radiance products of geostationary orbiting platforms provide virtually continuous observations of Earth. In contrast, low-Earth orbiters observe passive microwave (PMW) radiances at any location much less frequently. Prior literature demonstrates the ability of a machine learning (ML) approach to build a link between these two complementary radiance spectra by predicting PMW observations using infrared and visible products collected from geostationary instruments, which could potentially deliver a highly desirable synthetic PMW product with nearly continuous spatiotemporal coverage. However, current ML models lack the ability to provide a measure of uncertainty of such a product, significantly limiting its applications. In this work, Bayesian deep learning is employed to generate synthetic Global Precipitation Measurement (GPM) Microwave Imager (GMI) data from Advanced Baseline Imager (ABI) observations with attached uncertainties over the ocean. The study first uses deterministic residual networks (ResNets) to generate synthetic GMI brightness temperatures with as little mean absolute error as 1.72 K at the ABI spatiotemporal resolution. Then, for the same task, we use three Bayesian ResNet models to produce a comparable amount of error while providing previously unavailable predictive variance (i.e., uncertainty) for each synthetic data point. We find that the Flipout configuration provides the most robust calibration between uncertainty and error across GMI frequencies, and then demonstrate how this additional information is useful for discarding high-error synthetic data points prior to use by downstream applications.

Improved cloudy-sky snow albedo estimates using passive microwave and VIIRS data

Land surface albedo (LSA) is an essential component of the surface radiation budget, and has been retrieved extensively as a basic remote sensing product; however, daily LSA products suffer from extensive data gaps primarily caused by cloud cover. Accordingly, several gap-filling methods were developed (e.g., spatiotemporal interpolation and data fusion with albedo climatology), although the traditional methods are limited by cloud scale and surface heterogeneity. Further, as the largest varying surface landscape feature, seasonal snow cover substantially influences LSA and represents a major uncertainty factor of gap recovery because previous studies failed to employ actual surface signals to capture such ephemeral but intense albedo changes under cloud cover. To address this issue, a three-step framework was proposed for estimating 1 km cloudy-sky LSA using passive microwave (PMW) data, albedo climatology, and Visible Infrared Imaging Radiometer Suite (VIIRS) clear-sky albedo: (1) All-sky snow albedo was estimated from PMW brightness temperatures using a statistical model, (2) Continuous albedo dynamics were generated by combining the all-sky snow albedo with snow-free climatological albedo, and (3) The 1 km cloudy-sky LSA was predicted after filtering 1 km VIIRS clear-sky LSA by the albedo dynamic series. PMW-derived snow albedo was assessed over the Contiguous US (CONUS), and the final 1 km cloudy-sky LSA was validated across 10 sites from SURFRAD and Core AmeriFlux in 2013. Based on the comparison with high-quality MODIS pixels, the estimated snow albedo yielded an overall RMSE of 0.064 over CONUS, with a bias of −0.010 (R2 = 0.845). The recovered 1 km cloudy-sky LSA produced RMSEs of 0.074 (0.137) for all (snow) samples, a significant improvement over the Global Land Surface Satellite (GLASS) gap-free albedo products especially on snow cases (p-value = 0.027). Corresponding RMSE in calculating surface net radiation was also decreased by 38.91 W·m−2; and anomalous snow samples were corrected as well. The temporal analysis and all-sky LSA mapping suggest that the recovered LSA has satisfactory spatiotemporal continuity, and successfully captured details of spatiotemporal variability, especially for ephemeral snow events. This study provides an innovative solution to recover gaps in LSA data, and considerably improves the LSA accuracy under cloud cover, which can inform snow melting modeling, hazard forecasting, and irrigation management.

Inference of Soil Freezing Front Depth During the Freezing Period From the L-Band Passive Microwave Brightness Temperature

The freezing front depth ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">zff</i> ) of annual freeze- thaw cycles is critical for monitoring the dynamics of the cryosphere under climate change because <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">zff</i> 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 seasonal 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 freezing front depth ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">zff</i> ) with the diurnal amplitude variation (DAV) of L band brightness temperature ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ΔTB</i> ) during the freezing period. The new algorithm assumes 1) the daily-scale solar radiation heating/cooling effect causes the daily surface thawing depth ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ztf</i> ) variation, which leads further to <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ΔTB</i> ; 2) Δ <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> can be captured by an L-band radiometer; 3) <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ztf</i> and <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">zff</i> 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 <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> observations from a tower- based L-band radiometer (ELBARA-III) at Maqu are used to validate the BT-FF model. It shows that: 1) <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ΔTB</i> can be precisely estimated from <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ztf</i> during the daytime; 2) the decreasing of <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">ztf</i> is linearly related to the increase of <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">zff</i> with the Stefan Equation; 3) the accuracy of retrieved freezing front depth is about 5-25 cm; 4) the proposed model is applicable during the freezing period. The study is expected to extend the application of L -band <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">TB</i> data in cryosphere/meteorology and construct global freezing depth data set in the future.

Impact of profile-averaged soil ice fraction on passive microwave brightness temperature Diurnal Amplitude Variations (DAV) at L-band

The dynamic change of frozen soil is crucial to land-surface modeling, carbon feedback studies, ground engineering (e.g., constructions), and microwave remote sensing. L-Band satellite missions Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) are currently exploited to characterize soil into freeze or thaw (FT) states. However, brightness temperatures (TB) at L-band contain more information besides the FT state, particularly over permafrost or seasonally frozen soil, which has not been explored via current retrieval algorithms. To examine the potential for L-band TB observations, we define an index called Profile-Averaged Frozen Soil Fraction (Ff) related to Diurnal Amplitude Variation (DAV) of TB (i.e., ΔTB) based on the optical depth of the frozen soil column. We evaluated Ff inferred from the 0th-order microwave transfer model with the SMAP L1c product, the ground-based European Space Agency L-Band Radiometer III (ELBARA-III) TB observations, and temperature profiles collected at the Maqu station in northeastern Tibet. While there is a clear relationship between Ff and ΔTB, no apparent link exists with the ice content fraction (Ffi) within a fixed-depth soil column. The proposed model certifies that the profile-averaged soil ice content Ff relates to the dynamic microwave penetration depth by math and field measurement. The model reproduces well ΔTB in Period Freezing but has problems in Period Thawing when melted surface water obstruct the microwave signals. Our findings can be used to exploit ΔTB between 6 am and 6 pm, as a typically overpassing time by the SMOS and SMAP satellites, for estimating Ff, which can be further applied in weather/climate forecasting and for improving land-surface modeling.

A 41-year (1979–2019) passive-microwave-derived lake ice phenology data record of the Northern Hemisphere

Abstract. Seasonal ice cover is one of the important attributes of lakes in middle- and high-latitude regions. The annual freeze-up and breakup dates as well as the duration of ice cover (i.e., lake ice phenology) are sensitive to the weather and climate; hence, they can be used as an indicator of climate variability and change. In addition to optical, active microwave, and raw passive microwave data that can provide daily observations, the Calibrated Enhanced-Resolution Brightness Temperature (CETB) dataset available from the National Snow and Ice Data Center (NSIDC) provides an alternate source of passive microwave brightness temperature (TB) measurements for the determination of lake ice phenology on a 3.125 km grid. This study used Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave/Imager (SSM/I), and Special Sensor Microwave Imager/Sounder (SSMIS) data from the CETB dataset to extract the ice phenology for 56 lakes across the Northern Hemisphere from 1979 to 2019. According to the differences in TB between lake ice and open water, a threshold algorithm based on the moving t test method was applied to determine the lake ice status for grids located at least 6.25 km away from the lake shore, and the ice phenology dates for each lake were then extracted. When ice phenology could be extracted from more than one satellite over overlapping periods, results from the satellite offering the largest number of observations were prioritized. The lake ice phenology results showed strong agreement with an existing product derived from Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) and Advanced Microwave Scanning Radiometer 2 (AMSR2) data (2002 to 2015), with mean absolute errors of ice dates ranging from 2 to 4 d. Compared with near-shore in situ observations, the lake ice results, while different in terms of spatial coverage, still showed overall consistency. The produced lake ice record also displayed significant consistency when compared to a historical record of annual maximum ice cover of the Laurentian Great Lakes of North America. From 1979 to 2019, the average complete freezing duration and ice cover duration for lakes forming a complete ice cover on an annual basis were 153 and 161 d, respectively. The lake ice phenology dataset – a new climate data record (CDR) – will provide valuable information to the user community about the changing ice cover of lakes over the last 4 decades. The dataset is available at https://doi.org/10.1594/PANGAEA.937904 (Cai et al., 2021).

Prediction of Pan-Arctic Sea Ice Using Attention-Based LSTM Neural Networks

Within the rapidly changing Arctic region, accurate sea ice forecasts are of crucial importance for navigation activities, such as the planning of shipping routes. Numerical climate models have been widely used to generate Arctic sea ice forecasts at different time scales, but they are highly dependent on the initial conditions and are computationally expensive. Recently, with the increasing availability of geoscience data and the advances in deep learning algorithms, the use of artificial intelligence (AI)-based sea ice prediction methods has gained significant attention. In this study, we propose a supervised deep learning approach, namely attention-based long short-term memory networks (LSTMs), to forecast pan-Arctic sea ice at monthly time scales. Our method makes use of historical sea ice concentration (SIC) observations during 1979–2020, from passive microwave brightness temperatures. Based on the persistence of SIC anomalies, which is known as one of the dominant sources of sea ice predictability, our approach exploits the temporal relationships of sea ice conditions across different time windows of the training period. We demonstrate that the attention-based LSTM is able to learn the variations of the Arctic sea ice and can skillfully forecast pan-Arctic SIC on monthly time scale. By designing the loss function and utilizing the attention mechanism, our approach generally improves the accuracy of sea ice forecasts compared to traditional LSTM networks. Moreover, it outperforms forecasts with the climatology and persistence based empirical models, as well as two dynamical models from the Copernicus Climate Change Service (C3S) datastore. This approach shows great promise in enhancing forecasts of Arctic sea ice using AI methods.

Deep Neural Network High Spatiotemporal Resolution Precipitation Estimation (Deep-STEP) Using Passive Microwave and Infrared Data

Abstract Recent developments in “headline-making” deep neural networks (DNNs), specifically convolutional neural networks (CNNs), along with advancements in computational power, open great opportunities to integrate massive amounts of real-time observations to characterize spatiotemporal structures of surface precipitation. This study aims to develop a CNN algorithm, named Deep Neural Network High Spatiotemporal Resolution Precipitation Estimation (Deep-STEP), that ingests direct satellite passive microwave (PMW) brightness temperatures (Tbs) at emission and scattering frequencies combined with infrared (IR) Tbs from geostationary satellites and surface information to automatically extract geospatial features related to the precipitable clouds. These features allow the end-to-end Deep-STEP algorithm to instantaneously map surface precipitation intensities with a spatial resolution of 4 km. The main advantages of Deep-STEP, as compared to current state-of-the-art techniques, are 1) it learns and estimates complex precipitation systems directly from raw measurements in near–real time, 2) it uses the automatic spatial neighborhood feature extraction approach, and 3) it fuses coarse-resolution PMW footprints with IR images to reliably retrieve surface precipitation at a high spatial resolution. We anticipate our proposed DNN algorithm to be a starting point for more sophisticated and efficient precipitation retrieval systems in terms of accuracy, fine spatial pattern detection skills, and computational costs.

Synergistic Evaluation of Passive Microwave and Optical/IR Data for Modelling Vegetation Transmissivity towards Improved Soil Moisture Retrieval.

Vegetation cover and soil surface roughness are vital parameters in the soil moisture retrieval algorithms. Due to the high sensitivity of passive microwave and optical observations to Vegetation Water Content (VWC), this study assesses the integration of these two types of data to approximate the effect of vegetation on passive microwave Brightness Temperature (BT) to obtain the vegetation transmissivity parameter. For this purpose, a newly introduced index named Passive microwave and Optical Vegetation Index (POVI) was developed to improve the representativeness of VWC and converted into vegetation transmissivity through linear and nonlinear modelling approaches. The modified vegetation transmissivity is then applied in the Simultaneous Land Parameters Retrieval Model (SLPRM), which is an error minimization method for better retrieval of BT. Afterwards, the Volumetric Soil Moisture (VSM), Land Surface Temperature (LST) as well as canopy temperature (TC) were retrieved through this method in a central region of Iran (300 × 130 km2) from November 2015 to August 2016. The algorithm validation returned promising results, with a 20% improvement in soil moisture retrieval.

Global Snow Depth Retrieval From Passive Microwave Brightness Temperature With Machine Learning Approach

Current global snow retrieval algorithms based on spaceborne microwave measurements inherit noticeable biases and uncertainties regarding spatial distribution and temporal variations. In this article, we present an improved spatiotemporally dynamic global snow depth retrieval algorithm to account for the heterogeneity of snowpacks in different seasons worldwide. The proposed model adopts nonlinear machine learning to retrieve snow depths from passive microwave measurements and other auxiliary information. We indirectly characterized the variation in snow grain size using the daily profiles of the temperature gradient within the snowpack. In addition, a zoning and multitemporal modeling strategy was employed to reduce the bias and uncertainty caused by snow heterogeneity across different ecoregions and seasons. The proposed model was implemented to retrieve the global daily snow depth from 2001 to 2010. The results were validated by <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in situ</i> observations and compared with the NASA Advanced Microwave Scanning Radiometer for EOS (AMSR-E) snow water equivalent product (AE_DySno). Satisfactory accuracy was achieved for different ecoregions with regard to daily, monthly, and yearly validations (the root-mean-square error (RMSE) varied from ~7.5 to ~12 cm; the Pearson correlation coefficient <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$R$ </tex-math></inline-formula> ranged from 0.75 to 0.85). The results of ten trials indicated the promising stability of the proposed model in different ecoregions with small variations in RMSE and <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$R$ </tex-math></inline-formula> values. Compared with the AE_DySno products, the estimation results did not exhibit the overestimation problem and provided snow depth patterns with greater spatial heterogeneity, showing RMSEs ~5 cm lower and <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$R$ </tex-math></inline-formula> values ~0.3 higher than those of the AE_DySno products.

A Fine-Resolution Snow Depth Retrieval Algorithm From Enhanced-Resolution Passive Microwave Brightness Temperature Using Machine Learning in Northeast China

As one of the major components in the hydrological system, seasonal snow cover in Northeast China has drawn much attention recently. Because of the coarse spatial resolution of the passive microwave (PMW), heterogeneity of snowpack, and forest cover, it is difficult for existing snow products to achieve high precision snow parameters (e.g. snow depth (SD) or snow water equivalent (SWE)) assessment and hydrological research in fine scale. In this study, a novel SD retrieval algorithm that considered both the spatiotemporal dynamic of snow characteristics and forest attenuation was developed by combining the Calibrated Enhanced Resolution Brightness Temperature (TB) data and other auxiliary information, and produced a fine resolution (i.e., 6.25 km &#x00D7; 6.25 km) and high accuracy SD data in Northeast China. Instead of complex physical models, the machine learning was used to untangle the nonlinear complex relationship between SD and the enhanced resolution TB, forest fraction (FF), and snow characteristics. The verification results at ground weather stations showed that the retrieved SD by the proposed algorithm had high consistency with the observed SD, its RMSE, bias, and correlation coefficient (R) of 6.32 cm, -0.23 cm, and 0.63, respectively. Compared with the existing SD products (WESTDC and AMSR2), the developed model greatly improved both in spatial resolution and retrieval accuracy. In general, the fine-resolution SD inversion model achieved satisfactory accuracy and stability, and it will be used to generate long-term SD dataset service for climate change and hydrological research in the future.

Estimating Terrestrial Snow Mass via Multi‐Sensor Assimilation of Synthetic AMSR‐E Brightness Temperature Spectral Differences and Synthetic GRACE Terrestrial Water Storage Retrievals

AbstractThis study explores multi‐sensor data assimilation (DA) using synthetic Advanced Microwave Scanning Radiometer (AMSR‐E) passive microwave brightness temperature spectral differences () and synthetic Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage (TWS) retrievals to improve estimates of snow water equivalent (SWE), subsurface water storage, and TWS across snow‐covered terrain. Results show that multi‐sensor DA improves SWE estimates by reducing the RMSE by 14.1% relative to a model‐only simulation. Multi‐sensor assimilation also yields the smallest TWS RMSE (reduced by 13.0% relative to a model‐only simulation). However, multi‐sensor DA does not always yield complementary updates, and can sometimes lead to conflicting changes to SWE, where the assimilation of synthetic generates positive SWE increments while the assimilation of synthetic TWS removes SWE, which can ultimately degrade the posterior SWE estimates. This synthetic experiment provides useful insight for future DA experiments using real‐world AMSR‐E/AMSR‐2 observations and GRACE/GRACE‐FO TWS retrievals to better characterize terrestrial freshwater storage across regional scales.

An investigation on microwave transmissivity at frequencies of 18.7 and 36.5 GHz for diverse forest types during snow season

ABSTRACT Forests have invariably been considered as an obstacle in retrieving land surface parameters from spaceborne passive microwave brightness temperature (TB) observations. For quantifying the effect of forests on microwave signals, several models have been developed. However, these models rarely reveal the dependence of microwave radiation on forest types, which can hardly meet the needs of high-accuracy retrieval of terrestrial parameters in forested regions. A ground-based microwave radiometric observation experiment was designed to investigate the dependence of microwave radiation on frequency, polarization, and forest type. Downward TB at 18.7- and 36.5-GHz for horizontal- and vertical-polarization from the forest canopy was measured at 14 sample plots in Northeast China, along with snowpack and forest structural parameters. By providing fits to experimental data, new empirical transmissivity models for three forest types were developed, as a function of woody stem volume and depending on the frequency/polarization. The proposed models give diverse asymptotic transmissivity saturation levels and the corresponding saturation point of woody stem volume for different forest types. Root-mean-square error results between TB simulations and Advanced Microwave Scanning Radiometer-2 observations are approximately 3–6 K. This study provides an experimental and theoretical reference for further development of inversion models for snow parameters in forested areas.

Surface melting over the Greenland ice sheet derived from enhanced resolution passive microwave brightness temperatures (1979–2019)

Abstract. Surface melting is a major component of the Greenland ice sheet surface mass balance, and it affects sea level rise through direct runoff and the modulation of ice dynamics and hydrological processes, supraglacially, englacially and subglacially. Passive microwave (PMW) brightness temperature observations are of paramount importance in studying the spatial and temporal evolution of surface melting due to their long temporal coverage (1979–present) and high temporal resolution (daily). However, a major limitation of PMW datasets has been the relatively coarse spatial resolution, which has historically been of the order of tens of kilometers. Here, we use a newly released PMW dataset (37 GHz, horizontal polarization) made available through a NASA “Making Earth System Data Records for Use in Research Environments” (MeASUREs) program to study the spatiotemporal evolution of surface melting over the Greenland ice sheet at an enhanced spatial resolution of 3.125 km. We assess the outputs of different detection algorithms using data collected by automatic weather stations (AWSs) and the outputs of the Modèle Atmosphérique Régional (MAR) regional climate model. We found that sporadic melting is well captured using a dynamic algorithm based on the outputs of the Microwave Emission Model of Layered Snowpack (MEMLS), whereas a fixed threshold of 245 K is capable of detecting persistent melt. Our results indicate that, during the reference period from 1979 to 2019 (from 1988 to 2019), surface melting over the ice sheet increased in terms of both duration, up to 4.5 (2.9) d per decade, and extension, up to 6.9 % (3.6 %) of the entire ice sheet surface extent per decade, according to the MEMLS algorithm. Furthermore, the melting season started up to 4.0 (2.5) d earlier and ended 7.0 (3.9) d later per decade. We also explored the information content of the enhanced-resolution dataset with respect to the one at 25 km and MAR outputs using a semi-variogram approach. We found that the enhanced product is more sensitive to local-scale processes, thereby confirming the potential of this new enhanced product for monitoring surface melting over Greenland at a higher spatial resolution than the historical products and for monitoring its impact on sea level rise. This offers the opportunity to improve our understanding of the processes driving melting, to validate modeled melt extent at high resolution and, potentially, to assimilate these data in climate models.

Snow Depth Fusion Based on Machine Learning Methods for the Northern Hemisphere

In this study, a machine learning algorithm was introduced to fuse gridded snow depth datasets. The input variables of the machine learning method included geolocation (latitude and longitude), topographic data (elevation), gridded snow depth datasets and in situ observations. A total of 29,565 in situ observations were used to train and optimize the machine learning algorithm. A total of five gridded snow depth datasets—Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) snow depth, Global Snow Monitoring for Climate Research (GlobSnow) snow depth, Long time series of daily snow depth over the Northern Hemisphere (NHSD) snow depth, ERA-Interim snow depth and Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) snow depth—were used as input variables. The first three snow depth datasets are retrieved from passive microwave brightness temperature or assimilation with in situ observations, while the last two are snow depth datasets obtained from meteorological reanalysis data with a land surface model and data assimilation system. Then, three machine learning methods, i.e., Artificial Neural Networks (ANN), Support Vector Regression (SVR), and Random Forest Regression (RFR), were used to produce a fused snow depth dataset from 2002 to 2004. The RFR model performed best and was thus used to produce a new snow depth product from the fusion of the five snow depth datasets and auxiliary data over the Northern Hemisphere from 2002 to 2011. The fused snow-depth product was verified at five well-known snow observation sites. The R2 of Sodankylä, Old Aspen, and Reynolds Mountains East were 0.88, 0.69, and 0.63, respectively. At the Swamp Angel Study Plot and Weissfluhjoch observation sites, which have an average snow depth exceeding 200 cm, the fused snow depth did not perform well. The spatial patterns of the average snow depth were analyzed seasonally, and the average snow depths of autumn, winter, and spring were 5.7, 25.8, and 21.5 cm, respectively. In the future, random forest regression will be used to produce a long time series of a fused snow depth dataset over the Northern Hemisphere or other specific regions.

Diurnal Cycle of Passive Microwave Brightness Temperatures over Land at a Global Scale

Satellite-borne passive microwave radiometers provide brightness temperature (TB) measurements in a large spectral range which includes a number of frequency channels and generally two polarizations: horizontal and vertical. These TBs are widely used to retrieve several atmospheric and surface variables and parameters such as precipitation, soil moisture, water vapor, air temperature profile, and land surface emissivity. Since TBs are measured at different microwave frequencies with various instruments and at various incidence angles, spatial resolutions, and radiometric characteristics, a mere direct integration of them from different microwave sensors would not necessarily provide consistency. However, when appropriately harmonized, they can provide a complete dataset to estimate the diurnal cycle. This study first constructs the diurnal cycle of land TBs using the non-sun-synchronous Global Precipitation Measurement (GPM) Microwave Imager (GMI) observations by utilizing a cubic spline fit. The acquisition times of GMI vary from day to day and, therefore, the shape (amplitude and phase) of the diurnal cycle for each month is obtained by merging several days of measurements. This diurnal pattern is used as a point of reference when intercalibrated TBs from other passive microwave sensors with daily fixed acquisition times (e.g., Special Sensor Microwave Imager/Sounder, and Advanced Microwave Scanning Radiometer 2) are used to modify and tune the monthly diurnal cycle to daily diurnal cycle at a global scale. Since the GMI does not cover polar regions, the proposed method estimates a consistent diurnal cycle of land TBs at global scale. Results show that the shape and peak of the constructed TB diurnal cycle is approximately similar to the diurnal cycle of land surface temperature. The diurnal brightness temperature range for different land cover types has also been explored using the derived diurnal cycle of TBs. In general, a large diurnal TB range of more than 15 K has been observed for the grassland, shrubland, and tundra land cover types, whereas it is less than 5K over forests. Furthermore, seasonal variations in the diurnal TB range for different land cover types show a more consistent result over the Southern Hemisphere than over the Northern Hemisphere. The calibrated TB diurnal cycle may then be used to consistently estimate the diurnal cycle of land surface emissivity. Moreover, since changes in land surface emissivity are related to moisture change and freeze–thaw (FT) transitions in high-latitude regions, the results of this study enhance temporal detection of FT state, particularly during the transition times when multiple FT changes may occur within a day.