Spaceborne Microwave Radiometers Research Articles

Pan-tropical daily L-band microwave land surface emissivity retrieval from GNSS-R observations

ABSTRACT The uncertainties in microwave land surface emissivity (MLSE) measurements have long limited the use of spaceborne microwave radiometer data. As an emerging observation method, Global Navigation Satellite System Reflectometry (GNSS-R) has demonstrated great potential in several land and ocean applications. In this study, a method for obtaining daily MLSE dataset in the pan-tropical region from Cyclone GNSS (CYGNSS) observations is presented and evaluated. The CYGNSS observations are first aggregated into the Equal-Area-Scalable-Earth (EASE) 2.0 36 km grid by a combined weight function of distance, time, and signal-to-noise ratio variance. Then, the method employs a pixel-by-pixel regression algorithm to conduct the daily MLSE retrieval using reference emissivity derived from the Soil Moisture Active Passive (SMAP) brightness temperature. The CYGNSS MLSE shows good agreement with SMAP MLSE, delivering an overall root-mean-square error (RMSE) of 0.022 and 0.017 for horizontal and vertical polarization, respectively, during the training set spanning the whole year of 2018. Furthermore, on the test set from January 2019 to May 2019, the RMSE values amounted to 0.030 and 0.023 for horizontal and vertical polarization, respectively. Temperature records from the International Soil Moisture Network are employed to calculate the emissivity and for in-situ validation, which yield an RMSE of 0.034 and 0.026 for the two polarizations, respectively. The proposed algorithm provides an encouraging approach to obtain accurate daily MLSE dataset for microwave remote sensing. Compared to the SMAP MLSE, the CYGNSS MLSE has a remarkable improvement of 86% in temporal resolution, greatly complementing the existing microwave emissivity datasets.

Modeling Snow and Ice Microwave Emissions in the Arctic for a Multi‐Parameter Retrieval of Surface and Atmospheric Variables From Microwave Radiometer Satellite Data

AbstractMonitoring surface and atmospheric parameters—like water vapor—is challenging in the Arctic, despite the daily Arctic‐wide coverage of spaceborne microwave radiometer data. This is mainly due to the difficulties in characterizing the sea ice surface emission: sea ice and snow microwave emission is high and highly variable. There are very few data sets combining relevant in situ measurements with co‐located remote sensing data, which further complicates the development of accurate retrieval algorithms. Here, we present a multi‐parameter retrieval based on the inversion of a forward model for both, atmosphere and surface, for non‐melting conditions. The model consists of a layered microwave emission model of snow and ice. Since snow scattering and emission effects, as well as temperature gradients, are taken into account, a high variability in brightness temperatures can be simulated. For ocean regions and the atmosphere existing parameterized forward models are used. By using optimal estimation, the forward model can be inverted allowing for the simultaneous and consistent retrieval of nine variables: integrated water vapor, liquid water path, sea ice concentration, multi‐year ice fraction, snow depth, snow‐ice interface temperature and snow‐air interface temperature as well as sea‐surface temperature and wind speed (over open ocean). In addition, the method provides retrieval uncertainty estimates for each retrieved parameter. To evaluate the forward model as well as the retrieval, we use the extensive data sets acquired during the year‐long Arctic expedition Multidisciplinary drifting Observatory for the Study of Arctic Climate (2019–2020) as a reference.

Reviews and syntheses: Recent advances in microwave remote sensing in support of terrestrial carbon cycle science in Arctic–boreal regions

Abstract. Spaceborne microwave remote sensing (300 MHz–100 GHz) provides a valuable method for characterizing environmental changes, especially in Arctic–boreal regions (ABRs) where ground observations are generally spatially and temporally scarce. Although direct measurements of carbon fluxes are not feasible, spaceborne microwave radiometers and radar can monitor various important surface and near-surface variables that affect terrestrial carbon cycle processes such as respiratory carbon dioxide (CO2) fluxes; photosynthetic CO2 uptake; and processes related to net methane (CH4) exchange including CH4 production, transport and consumption. Examples of such controls include soil moisture and temperature, surface freeze–thaw cycles, vegetation water storage, snowpack properties and land cover. Microwave remote sensing also provides a means for independent aboveground biomass estimates that can be used to estimate aboveground carbon stocks. The microwave data record spans multiple decades going back to the 1970s with frequent (daily to weekly) global coverage independent of atmospheric conditions and solar illumination. Collectively, these advantages hold substantial untapped potential to monitor and better understand carbon cycle processes across ABRs. Given rapid climate warming across ABRs and the associated carbon cycle feedbacks to the global climate system, this review argues for the importance of rapid integration of microwave information into ABR terrestrial carbon cycle science.

The extraordinarily large vortex structure of Typhoon In-fa (2021), observed by spaceborne microwave radiometer and synthetic aperture radar

Tropical cyclone (TC) vortex structure can be estimated from wind radii (34, 50, and 64 kt). Accurate wind radii are essential for assessing the impact range of TCs. Combining observations from active and passive microwave remote sensing instruments can provide long-time series data for monitoring changes in TC wind structure. Here, the evolution of the wind radii for Typhoon In-fa (2021), from its genesis to its first landfall, is evaluated using data from spaceborne microwave radiometers and synthetic aperture radars (SARs). Our results show that the retrieved wind radii are relatively close to those from the Joint Typhoon Warning Center (JTWC), the coastal automatic weather stations in China, and the advanced scatterometer. The time series for the wind radii show significant increasing trends during the lifecycle of In-fa and a massive vortex structure before In-fa made landfall. The average wind radii observed before landfall are all likely in the 90th percentile of the 2001–2021 JTWC historical wind radii records, indicating that the vortex structure is extraordinarily large before In-fa made landfall and represents a historical extreme. The extremely large vortex structure may be attributable to the interaction between In-fa and an extremely strong monsoon gyre.

Delayed Antarctic melt season reduces albedo feedback.

Antarctica's response to climate change varies greatly both spatially and temporally. Surface melting impacts mass balance and also lowers surface albedo. We use a 43-year record (from 1978 to 2020) of Antarctic snow melt seasons from space-borne microwave radiometers with a machine-learning algorithm to show that both the onset and the end of the melt season are being delayed. Granger-causality analysis shows that melt end is delayed due to increased heat flux from the ocean to the atmosphere at minimum sea-ice extent from warming oceans. Melt onset is Granger-caused primarily by the turbulent heat flux from ocean to atmosphere that is in turn driven by sea-ice variability. Delayed snowmelt season leads to a net decrease in the absorption of solar irradiance, as a delayed summer means that higher albedo occurs after the period of maximum solar radiation, which changes Antarctica's radiation balance more than sea-ice cover.

A Cloud Detection Neural Network Approach for the Next Generation Microwave Sounder Aboard EPS MetOp-SG A1

This work presents an algorithm based on a neural network (NN) for cloud detection to detect clouds and their thermodynamic phase using spectral observations from spaceborne microwave radiometers. A standalone cloud detection algorithm over the ocean and land has been developed to distinguish clear sky versus ice and liquid clouds from microwave sounder (MWS) observations. The MWS instrument—scheduled to be onboard the first satellite of the Eumetsat Polar System Second-Generation (EPS-SG) series, MetOp-SG A1—has a direct inheritance from advanced microwave sounding unit A (AMSU-A) and the microwave humidity sounder (MHS) microwave instruments. Real observations from the MWS sensor are not currently available as its launch is foreseen in 2024. Thus, a simulated dataset of atmospheric states and associated MWS synthetic observations have been produced through radiative transfer calculations with ERA5 real atmospheric profiles and surface conditions. The developed algorithm has been validated using spectral observations from the AMSU-A and MHS sounders. While ERA5 atmospheric profiles serve as references for the model development and its validation, observations from AVHRR cloud mask products provide references for the AMSU-A/MHS model evaluation. The results clearly show the NN algorithm’s high skills to detect clear, ice and liquid cloud conditions against a benchmark. In terms of overall accuracy, the NN model features 92% (88%) on the ocean and 87% (85%) on land, for the MWS (AMSU-A/MHS)-simulated dataset, respectively.

Analysis of the Effect of Enhanced FOV and Sampling Strategy on the Spatial Resolution Enhancement of Spaceborne Microwave Radiometer

Analysis of the Effect of Enhanced FOV and Sampling Strategy on the Spatial Resolution Enhancement of Spaceborne Microwave Radiometer

Optimal Channel Selection of Spaceborne Microwave Radiometer for Surface Pressure Retrieval over Oceans

Abstract The O2-band channel configuration of existing microwave radiometers is not optimal for surface pressure retrieval, which limits the surface pressure retrieval accuracy. In this study, we present the results of theoretically what might be the optimal microwave channels for surface pressure retrieval. An improved iterative selection method is used to select the channels that contain the highest cumulative content of surface pressure information. The selected optimal channel set comprises 16 channels, among which 10 channels are centered at the 50–60 GHz oxygen absorption band and 6 channels are centered around the 118.75 GHz oxygen absorption line. Two representative spaceborne microwave radiometers are used for comparisons, the Advanced Technology Microwave Sounder (ATMS) on board the Suomi National Polar-Orbiting Partnership (SNPP) satellite and the Microwave Humidity and Temperature Sounder (MWHTS) on board the Chinese Fengyun-3C (FY-3C) satellite. The results of information content analysis show that the optimal channel set contains more surface pressure information than that of the combination of SNPP/ATMS and FY-3C/MWHTS (SNPP/ATMS+FY-3C/MWHTS) channels. A representative dataset from the ERA5 data is input into the plane-parallel Microwave Radiative Transfer model to obtain the simulated brightness temperature observations of the selected optimal channels and the SNPP/ATMS+FY-3C/MWHTS channels. Using the simulated observations, retrieval experiments are performed. Experimental results show that retrieval accuracies of the optimal channel set are 1.09 and 1.64 hPa for clear-sky and cloudy conditions, respectively. The retrieval accuracies are 0.60 and 0.65 hPa better than that of the SNPP/ATMS+FY-3C/MWHTS channels for clear-sky and cloudy conditions, respectively.

Improved Estimation of O-B Bias and Standard Deviation by an RFI Restoration Method for AMSR-2 C-Band Observations over North America

Spaceborne microwave radiometer observations play vital roles in surface parameter retrievals and data assimilation, but widespread radio-frequency interference (RFI) signals in the C-band channel result in a lack of valuable data over large areas. Establishing repaired data based on existing observation information is crucial. In this study, Advanced Microwave Scanning Radiometer (AMSR)-2 C-band data affected by RFI were accurately repaired through the iterative principal component analysis (PCA) method in 2016 over the U.S. land area. The standard deviation (STD) and bias characteristics of the brightness temperature in the C-band vertical polarization channel were compared and analyzed before and after the restoration to verify the assimilation application prospect of the repaired data. Not only was the spatial continuity of the microwave imager observations significantly improved following restoration; the STD and bias of the observation minus background (OMB) of the restored data were basically consistent with those of the RFI-free data. The STD of OMB exhibited obvious seasonal variations, which were approximately 4.0 K from January to May and 3.0 K from June to December, whereas the biases were near zero in winter but negative (approximately −2.0 K) in summer. The surface type and terrain height also critically affected the STD and bias. The STD decreased with increasing terrain height, whereas the bias exhibited the opposite trend. The STD was largest in low-vegetation areas (4.0 K) but only approximately 2.0–3.0 K in pine forest and brush areas. These results show that the restored data have a high prospect for retrieval application and assimilation, and the STD and bias estimation results also provide a reference for land-based AMSR-2 data assimilation.

Sea Surface Wind Retrieval under Rainy Conditions from Active and Passive Microwave Measurements

The space-borne microwave radiometers and scatterometers can effectively measure global sea surface winds under non-precipitation. However, the measurements in rainy conditions significantly degrade, which are usually flagged as poor quality or invalidated for some scientific purposes. This paper develops a combined active–passive wind vector retrieval model for rainy conditions based on the HY-2B radiometer and scatterometer measurements. In our model, the polarization ratio of brightness temperatures at 6.925 GHz (PR06) is used as an indicator to implicitly represent the rain effect. For wind speed retrieval, a statistical regression model is trained as a function of PR06 and brightness temperatures of the radiometer. Moreover, two new geophysical model functions, including rain effect, are developed for wind direction inversion. Comparisons between HY-2B retrieval results and ERA5 wind products indicate that the retrieval model performs well under all rainy conditions. The overall root mean squared errors (RMSEs) of wind speed and direction retrievals are 1.60 m/s and 20.60°, respectively. With an increase in the rain rate, the wind retrieval performance degrades slightly and still provides a reliable retrieval result.

Global Snowmelt Onset Reflects Climate Variability: Insights from Spaceborne Radiometer Observations

Abstract Snowmelt is a critical component in the cryosphere and has a direct impact on Earth’s energy and water budget. Here, a 40-yr integrated melt onset (MO) dataset over sea ice, ice sheets, and terrestrial snow is compiled from spaceborne microwave radiometers and ERA5, allowing an overall assessment of the cryosphere. Results suggest that MO in both hemispheres shows latitudinal and vertical zonalities. The global cryosphere presented a trend toward earlier MO (−2 days decade−1) with hotpots distributed at the Northern Hemisphere high latitudes where the warming rate is much higher than that at lower latitudes. Overall, variations in MO showed a similar pattern to that in near-surface temperature. The advance of MO has been slowing down since the 1990s and no significant trend was observed during the so-called warming hiatus period (1998–2012). Regionally, climatic linkage analyses suggest the local MO variations were associated with different climate indices. MO in the pan-Arctic region is related with the Arctic Oscillation and North Atlantic Oscillation, while that in the pan-Antarctic region is associated with El Niño–Southern Oscillation and the southern annular mode. Occasionally, abnormal MO occurs as a result of extreme weather conditions. In February 2018, abnormal early melt events that occurred in the Arctic Ocean are found to be linked with the warm southerly flow due to sudden stratospheric warming. These findings suggest the satellite-based MO allows examining the dynamics and extremes in the climate system, both regionally and globally.

Impact of dynamic snow density on GlobSnow snow water equivalent retrieval accuracy

Abstract. Snow water equivalent (SWE) is an important variable in describing global seasonal snow cover. Traditionally, SWE has been measured manually at snow transects or using observations from weather stations. However, these measurements have a poor spatial coverage, and a good alternative to in situ measurements is to use spaceborne passive microwave observations, which can provide global coverage at daily timescales. The reliability and accuracy of SWE estimates made using spaceborne microwave radiometer data can be improved by assimilating radiometer observations with weather station snow depth observations as done in the GlobSnow SWE retrieval methodology. However, one possible source of uncertainty in the GlobSnow SWE retrieval approach is the constant snow density used in modelling emission of snow. In this paper, three versions of spatially and temporally varying snow density fields were implemented using snow transect data from Eurasia and Canada and automated snow observations from the United States. Snow density fields were used to post-process the baseline GlobSnow v.3.0 SWE product. Decadal snow density information, i.e. fields where snow density for each day of the year was taken as the mean calculated for the corresponding day over 10 years, was found to produce the best results. Overall, post-processing GlobSnow SWE retrieval with dynamic snow density information improved overestimation of small SWE values and underestimation of large SWE values, though underestimation of SWE values larger than 175 mm was still significant.

Arctic sea ice cover data from spaceborne synthetic aperture radar by deep learning

Abstract. Widely used sea ice concentration and sea ice cover in polar regions are derived mainly from spaceborne microwave radiometer and scatterometer data, and the typical spatial resolution of these products ranges from several to dozens of kilometers. Due to dramatic changes in polar sea ice, high-resolution sea ice cover data are drawing increasing attention for polar navigation, environmental research, and offshore operations. In this paper, we focused on developing an approach for deriving a high-resolution sea ice cover product for the Arctic using Sentinel-1 (S1) dual-polarization (horizontal-horizontal, HH, and horizontal-vertical, HV) data in extra wide swath (EW) mode. The approach for discriminating sea ice from open water by synthetic aperture radar (SAR) data is based on a modified U-Net architecture, a deep learning network. By employing an integrated stacking model to combine multiple U-Net classifiers with diverse specializations, sea ice segmentation is achieved with superior accuracy over any individual classifier. We applied the proposed approach to over 28 000 S1 EW images acquired in 2019 to obtain sea ice cover products in a high spatial resolution of 400 m. The validation by 96 cases of visual interpretation results shows an overall accuracy of 96.10 %. The S1-derived sea ice cover was converted to concentration and then compared with Advanced Microwave Scanning Radiometer 2 (AMSR2) sea ice concentration data, showing an average absolute difference of 5.55 % with seasonal fluctuations. A direct comparison with Interactive Multisensor Snow and Ice Mapping System (IMS) daily sea ice cover data achieves an average accuracy of 93.98 %. These results show that the developed S1-derived sea ice cover results are comparable to the AMSR and IMS data in terms of overall accuracy but superior to these data in presenting detailed sea ice cover information, particularly in the marginal ice zone (MIZ). Data are available at https://doi.org/10.11922/sciencedb.00273 (Wang and Li, 2020).

Physical Retrieval of Rain Rate from Ground-Based Microwave Radiometry

Because of its clear physical meaning, physical methods are more often used for space-borne microwave radiometers to retrieve the rain rate, but they are rarely used for ground-based microwave radiometers that are very sensitive to rainfall. In this article, an opacity physical retrieval method is implemented to retrieve the rain rate (denoted as Opa-RR) using ground-based microwave radiometer data (21.4 and 31.5 GHz) of the tropospheric water radiometer (TROWARA) at Bern, Switzerland from 2005 to 2019. The Opa-RR firstly establishes a direct connection between the rain rate and the enhanced atmospheric opacity during rain, then iteratively adjusts the rain effective temperature to determine the rain opacity, based on the radiative transfer equation, and finally estimates the rain rate. These estimations are compared with the available simultaneous rain rate derived from rain gauge data and reanalysis data (ERA5). The results and the intercomparison demonstrate that during moderate rains and at the 31 GHz channel, the Opa-RR method was close to the actual situation and capable of the rain rate estimation. In addition, the Opa-RR method can well derive the changes in cumulative rain over time (day, month, and year), and the monthly rain rate estimation is superior, with the rain gauge validated R2 and the root-mean-square error value of 0.77 and 22.46 mm/month, respectively. Compared with ERA5, Opa-RR at 31GHz achieves a competitive performance.

The Diurnal Cycle of Precipitation according to Multiple Decades of Global Satellite Observations, Three CMIP6 Models, and the ECMWF Reanalysis

AbstractNASA Precipitation Measurement Mission observations are used to evaluate the diurnal cycle of precipitation from three CMIP6 models (NCAR-CESM2, CNRM-CM6.1, CNRM-ESM2.1) and the ERA5 reanalysis. NASA’s global-gridded IMERG product, which combines spaceborne microwave radiometer, infrared sensor, and ground-based gauge measurements, provides high-spatiotemporal-resolution (0.1° and half-hourly) estimates that are suitable for evaluating the diurnal cycle in models, as determined against the ground-based radar network over the conterminous United States. IMERG estimates are coarsened to the spatial and hourly resolution of the state-of-the-art CMIP6 and ERA5 products, and their diurnal cycles are compared across multiple decades of June–August in the 60°N–60°S domain (IMERG and ERA5: 2000–19; NCAR and CNRM: 1979–2008). Low-precipitation regions (and weak-amplitude regions when analyzing the diurnal phase) are excluded from analyses so as to assess only robust diurnal signals. Observations identify greater diurnal amplitudes over land (26%–134% of the precipitation mean; 5th–95th percentile) than over ocean (14%–66%). ERA5, NCAR, and CNRM underestimate amplitudes over ocean, and ERA5 overestimates over land. IMERG observes a distinct diurnal cycle only in certain regions, with precipitation peaking broadly between 1400 and 2100 LST over land (2100–0600 LST over mountainous and varying-terrain regions) and 0000 and 1200 LST over ocean. The simulated diurnal cycle is unrealistically early when compared with observations, particularly over land (NCAR-CESM2 AMIP: −1 h; ERA5: −2 h; CNRM-CM6.1 AMIP: −4 h on average) with nocturnal maxima not well represented over mountainous regions. Furthermore, ERA5’s representation of the diurnal cycle is too simplified, with less interannual variability in the time of maximum relative to observations over many regions.

Covariation of Passive–Active Microwave Measurements over Vegetated Surfaces: Case Studies at L-Band Passive and L-, C- and X-Band Active

The analysis of soil and land cover scattering properties and their connection with the parameters of microwave scattering is a longstanding research topic. Recently, the advent of modern space-borne microwave radiometers like SMAP in addition to the trend towards open data for scientific use fostered the development of enhanced models based on data fusion from different platforms permitting more accurate assessments. SMAP was designed to operate on an integrated combination of a radiometer and a radar, both operating in L-band. Unexpected failure of the radar component encouraged scientists to experiment various combination of data from the surviving radiometer with other sources of radar data, notably C-band Sentinel-1 data. In this work, we present a case study on a possible combination of SMAP radiometer data with X-band radar data from TerraSAR-X and COSMO-SkyMed, comparing results with those provided by NASA from their standard production procedures. The study was performed on two test sites, one at an agricultural site in Germany and one in the Brazilian Amazon, to explore very different vegetation conditions. This work is a part of a broader research effort addressing the combination of multiple sources of passive and active microwave sensing data. The research question defining this research effort is whether the use of data from multiple active sources affords either obtaining more accurate estimates of active–passive co-variation parameters for a given observation period, or shortening the minimum observation period by increasing the temporal density of active samples. In this framework, this paper addresses a preliminary comparison of fresh and past results obtained from C-, X-, and L-band active sensing data. The observed relations offer interesting clues on the impact of band selection on soil vegetation analysis.

Assessing the Spaceborne 183.31-GHz Radiometric Channel Geolocation Using High-Altitude Lakes, Ice Shelves, and SAR Imagery

The goal of this work is to perform the geolocation error assessment of the channel imagery at 183.31 GHz of the Special Sensor Microwave Imager/Sounder (SSMIS). The frequency around 183.31 GHz still represents the highest channel frequency of current spaceborne microwave and millimeter-wave radiometers. The latter will be extended to frequencies up to 664 GHz, as in the case of EUMETSAT Ice Cloud Imager (ICI). This use of submillimeter observations unfortunately prevents a straightforward geolocation error assessment using landmark-based techniques. We used SSMIS data at 183.31 GHz as a submillimeter proxy to identify the most suitable targets for geolocation error validation in very dry atmospheric conditions, as suggested by radiative transfer modeling. Using a yearly SSMIS data set, three candidates’ landmark targets are selected: 1) high-altitude lakes and high-latitude bays using a coastline reference database and 2) Antarctic ice shelves using coastlines derived from Sentinel-1 Synthetic Aperture Radar (SAR) imagery. Data processing is carried out by using spatial cross correlation methods in the spatial frequency domain and performing a numerical sensitivity analysis to contour displacement. Cloud masking, based on a fuzzy-logic approach, is applied to automatically selected clear-air days. The results show that the average geolocation error is about 6.2 km for mountainous lakes and sea bays and 5.4 km for ice shelves, with a standard deviation of about 2.7 and 2.0 km, respectively. The results are in line with SSMIS previous estimates, whereas annual clear-air days are about 10% for mountainous lakes and sea bays and 18% for ice shelves.

Raindrop Signature from Microwave Radiometer Over Deserts

AbstractRainfall estimates from spaceborne microwave radiometers form the foundation of global precipitation data sets. Since the beginning of the satellite microwave rainfall estimation era in the 1980s, the primary signature leveraged over land for these estimates has been the brightness temperature (TB) depression due to ice particle scattering. Contrary to this practice, time series analyses based on observations from two spaceborne radars and two spaceborne radiometers reveal a TB increase at H19 due to raindrop emission as the primary cloud particle signature over desert terrain. Low surface emissivity supports the use of liquid raindrop emission as the primary signature over desert surfaces. In these regions, the surface rain rate better correlates with the liquid raindrop emission signal than with the scattering induced by ice further aloft, suggesting a new potential for improving rainfall estimation over deserts by exploiting the liquid raindrop emission signature.

The global long-term microwave Vegetation Optical Depth Climate Archive (VODCA)

Abstract. Since the late 1970s, space-borne microwave radiometers have been providing measurements of radiation emitted by the Earth’s surface. From these measurements it is possible to derive vegetation optical depth (VOD), a model-based indicator related to the density, biomass, and water content of vegetation. Because of its high temporal resolution and long availability, VOD can be used to monitor short- to long-term changes in vegetation. However, studying long-term VOD dynamics is generally hampered by the relatively short time span covered by the individual microwave sensors. This can potentially be overcome by merging multiple VOD products into a single climate data record. However, combining multiple sensors into a single product is challenging as systematic differences between input products like biases, different temporal and spatial resolutions, and coverage need to be overcome. Here, we present a new series of long-term VOD products, the VOD Climate Archive (VODCA). VODCA combines VOD retrievals that have been derived from multiple sensors (SSM/I, TMI, AMSR-E, WindSat, and AMSR2) using the Land Parameter Retrieval Model. We produce separate VOD products for microwave observations in different spectral bands, namely the Ku-band (period 1987–2017), X-band (1997–2018), and C-band (2002–2018). In this way, our multi-band VOD products preserve the unique characteristics of each frequency with respect to the structural elements of the canopy. Our merging approach builds on an existing approach that is used to merge satellite products of surface soil moisture: first, the data sets are co-calibrated via cumulative distribution function matching using AMSR-E as the scaling reference. To do so, we apply a new matching technique that scales outliers more robustly than ordinary piecewise linear interpolation. Second, we aggregate the data sets by taking the arithmetic mean between temporally overlapping observations of the scaled data. The characteristics of VODCA are assessed for self-consistency and against other products. Using an autocorrelation analysis, we show that the merging of the multiple data sets successfully reduces the random error compared to the input data sets. Spatio-temporal patterns and anomalies of the merged products show consistency between frequencies and with leaf area index observations from the MODIS instrument as well as with Vegetation Continuous Fields from the AVHRR instruments. Long-term trends in Ku-band VODCA show that since 1987 there has been a decline in VOD in the tropics and in large parts of east-central and north Asia, while a substantial increase is observed in India, large parts of Australia, southern Africa, southeastern China, and central North America. In summary, VODCA shows vast potential for monitoring spatial–temporal ecosystem changes as it is sensitive to vegetation water content and unaffected by cloud cover or high sun zenith angles. As such, it complements existing long-term optical indices of greenness and leaf area. The VODCA products (Moesinger et al., 2019) are open access and available under Attribution 4.0 International at https://doi.org/10.5281/zenodo.2575599.

Change of snow and ice melting time in High Mountain Asia

High Mountain Asia (HMA) is very sensitive to climate changes. In HMA, air temperature and precipitation shifts or increases are reflected in the timing of snowmelt onset. In this study, a long-term series (1979–2018) of snow melt onset time is first derived using spaceborne microwave radiometer data, following which the long-term trend of snow and ice-melt time of Tienshan, Altay, Karakorom, Hindu Kush, and Pamir is analyzed. Previous studies proposed many algorithms for detecting snowmelt onset and freeze-up using microwave remote sensing. The previously proposed algorithms were extensively applied to polar regions (including the Antarctic, Siberia, Greenland, and sea ice surface), where the impact of topography or mixed pixel problem is relatively small. However, for high mountain regions, complex topography and a potential mixed pixel issue would result in a very complicated and noisy satellite-observed active and passive microwave remote sensing signal. This study thus proposes a recently developed snow and ice-melt detection algorithm in which, for passive microwave data, a median filter is first applied to the original signal (brightness temperature) to suppress random, small, or short-duration signal variations. The differential average derivative of a particular date is then calculated using the time series first-order derivative as the average first-order derivative of a specific count of observations after this particular date. Generally, the differential average derivative is an indicator of sudden changes in time series brightness temperature. Results show that the melt onset time in the majority of HMA occurs earlier, except in the case of the Karakorom Mountains and a part of the West Kunlun Mountains. Moreover, the melt onset time derived from satellite morning (6:00 local time) pass data shows that the melt onset time of Karakorom, the West Kulun Mountains, and southeast Tibet remains stable, or even occurs later. The trend of southeast Tibet is unique, with its earlier melt based on afternoon (18:00 local time) pass data and its later melt based on morning pass data. This suggests that the melt–refreeze period of southeast Tibet is increasing together with the increasing diurnal temperature difference. Then, a 2-m air temperature in ERA5 reanalysis data is used for comparison with melt onset time for validation and analysis. ERA5 is the latest climate reanalysis produced by ECMWF, providing hourly data on many atmospheric, land-surface, and sea-state parameters. A strong correlation exists between monthly average air temperatures and melt onset time, with the maximum linear fit R 2 of 0.76. This strong correlation indicates the good data quality of ERA5 reanalysis and melt onset time. The unique trend of southeast Tibet can also be explained using ERA5 reanalysis data, which show that the monthly mean daily maximum air temperature is increasing, but the monthly mean daily minimum air temperature is decreasing. This renders a slope of the linear fit of melt onset data in the period of 1988–2018 in the whole HMA region. Because the satellite overpasses time difference, data from only 1988 to 2018 is used for trend analysis. The analysis of the relationship between the melt onset time change rate and elevation shows that the areas with an earlier melt are almost all located in low-elevation regions, and the rate of melt time change is positively correlated with elevation. This suggests that low-elevation regions are more affected by climate changes. This study provides objective evidence of the impact of climate change on the cryospheric system in HMA.