Comment on hess-2021-508

A spatial downscaling method for SMAP soil moisture through visible and shortwave-infrared remote sensing data

The Advanced Scatterometer (ASCAT), Soil Moisture Active Passive (SMAP), and European Space Agency-Climate Change Initiative (ESA–CCI) soil moisture (SM) products are widely used in agricultural drought monitoring, water resource management, and climate analysis applications. However, the performance of these SM products varies significantly across regions and environmental conditions, due to in sensor characteristics, retrieval algorithms, and the lack of localized calibration. This study proposes a multi-source long short-term memory (MLSTM) for improving ASCAT, ESA–CCI, and SMAP SM products by combining in-situ SM measurements and four key auxiliary variables: precipitation (PRE), land surface temperature (LST), fractional vegetation cover (FVC), and evapotranspiration (ET). First, the in-situ measured data from four in-situ observation networks were corrected using the LSTM method to match the grid sizes of ASCAT (0.1°), ESA–CCI (0.25°), and SMAP (0.1°) SM products. The RPE, LST, FVC, and ET were used as inputs to the LSTM to obtain loss data against in-situ SM measurements. Second, the ASCAT, ESA–CCI, and SMAP SM datasets were used as inputs to the LSTM to generate loss data, which were subsequently corrected using LSTM-derived loss data based on in-situ SM measurements. When the mean squared error (MSE) loss values were minimized, the improvement for ASCAT, ESA–CCI, and SMAP products was considered the best. Finally, the improved ASCAT, ESA–CCI, and SMAP were produced and evaluated by the correlation coefficient (R), root mean square error (RMSE), and standard deviation (SD). The results showed that the RMSE values of the improved ASCAT, ESA–CCI, and SMAP products against the corrected in-situ SM data in the OZNET network were lower, i.e., 0.014 cm3/cm3, 0.019 cm3/cm3, and 0.034 cm3/cm3, respectively. Compared with the ESA–CCI and SMAP products, the ASCAT product was greatly improved, e.g., in the SNOTEL network, the Root Mean-Square Deviation (RMSD) values of 0.1049 cm3/cm3 (ASCAT) and 0.0662 cm3/cm3 (improved ASCAT). Overall, the MLSTM-based algorithm has the potential to improve the global satellite SM product.

Abstract. Very high-resolution (∼10–100 m) surface soil moisture (SM) observations are important for applications in agriculture, among other purposes. This is the original goal of the S2MP (Sentinel-1/Sentinel-2-Derived Soil Moisture Product) algorithm, which was designed to retrieve surface SM at the agricultural plot scale by simultaneously using Sentinel-1 (S1) backscatter coefficients and Sentinel-2 (S2) NDVI (Normalized Difference Vegetation Index) as inputs to a neural network trained with Water Cloud Model simulations. However, for many applications, including hydrology and climate impact assessment at regional level, large maps with a high resolution (HR) of around 1 km are already a significant improvement with respect to most of the publicly available SM datasets, which have resolutions of about 25 km. In this study, the S2MP algorithm was adapted to work at 1 km resolution and extended from croplands to herbaceous vegetation types. A target resolution of 1 km also allows the evaluation of the interest in using NDVI derived from Sentinel-3 (S3) instead of S2. Two sets of SM maps at 1 km resolution were produced with S2MP over six regions of ∼104 km2 in Spain, Tunisia, North America, Australia, and the southwest and southeast regions of France for the whole year of 2019. The first set was derived from the combination of S1 and S2 data (S1 + S2 maps), while the second one was derived from the combination of S1 and S3 (S1 + S3 maps). S1 + S2 and S1 + S3 SM maps were compared to each other, to those of the 1 km resolution Copernicus Global Land Service (CGLS) SM and Soil Water Index (SWI) datasets, and to those of the Soil Moisture Active Passive (SMAP) + S1 product. The S2MP S1 + S2 and S1 + S3 SM maps are in very good agreement in terms of correlation (R≥0.9), bias (≤0.04 m3 m−3), and standard deviation of the difference (SDD≤0.03 m3 m−3) over the six domains investigated in this study. In a second step, the S1 + S3 S2MP maps were compared to the other HR maps. S1 + S3 SM maps are well correlated to the CGLS SM maps (R∼0.7–0.8), but the correlations with respect to the other HR maps (CGLS SWI and SMAP + S1) drop significantly over many areas of the six domains investigated in this study. The highest correlations between the HR maps were found over croplands and when the 1 km pixels have a very homogeneous land cover. The bias among the different maps was found to be significant over some areas of the six domains, reaching values of ±0.1 m3 m−3. The S1 + S3 maps show a lower SDD with respect to CGLS maps (≤0.06 m3 m−3) than with respect to the SMAP + S1 maps (≤0.1 m3 m−3) for all the six domains. Finally, all the HR datasets (S1 + S2, S1 + S3, CGLS, and SMAP + S1) were also compared to in situ measurements from five networks across five countries, along with coarse-resolution (CR) SM products from SMAP, SMOS, and the European Space Agency Climate Change Initiative (CCI). While all the CR and HR products show different bias and SDD, the HR products show lower correlations than the CR ones with respect to in situ measurements. The discrepancies in between the different HR datasets, except for the more simple land cover conditions (homogeneous pixels with croplands) and the lower performances with respect to in situ measurement than coarse-resolution datasets, show the remaining challenges for large-scale HR SM mapping.

Soil moisture plays a significant role in disciplines such as hydrology, meteorology and agriculture, and passive microwave remote sensing has become a widely used technique for global soil moisture estimation over the past three decades. Several satellite missions carrying radiometers have been launched over the past years. Among them are Japan Aerospace Exploration Agency's (JAXA's) Advanced Microwave Scanning Radiometer-EOS (AMSR-E) launched on NASA's Aqua satellite, European Space Agency's (ESA's) Soil Moisture and Ocean Salinity (SMOS) mission, JAXA's Advanced Microwave Scanning Radiometer 2 (AMSR2) onboard the GCOM-W satellite, and NASA's Soil Moisture Active Passive (SMAP) mission. Based on the availability of these four missions, there is an opportunity to develop a consistent inter-calibrated longterm soil moisture data record. This study focuses on the parametrization of the tau-omega model for soil moisture retrieval at L-, C- and X-band using brightness temperature observations from the four missions and in-situ soil moisture and soil temperature data from the SMAP core validation sites across various land cover types. The ancillary data sets used in the SMAP baseline algorithm are used for the retrievals at different frequencies. Simultaneous calibrations of the vegetation parameter $b$ and roughness parameter $h$ at both horizontal and vertical polarizations are performed. A set of model parameters to successfully retrieve soil moisture at different validation sites at L-, C- and X -band are presented. A preliminary comparison of SMAP and AMSR2 soil moisture retrievals against in situ observations at the Yanco (Australia) and TxSON (U.S.) sites showed the best accuracy and correlation at L-band (RMSD=0.03-0.05 m3/m3; R=0.92-0.94). The C-/X-band performance was not as satisfactory as L-band (RMSD=0.08-0.17 m3/m3; R=0.24-0.79). This also indicates that retrieval at higher frequencies can be very challenging when dense vegetation is present.

Global soil moisture (SM) mapping at high spatial and temporal resolution contributes significantly to various hydrologic and meteorological researches. This work presents an algorithm for combining fine-resolution Global Navigation Satellite System Reflectometry (GNSS-R) observations from the Cyclone Global Navigation Satellite System (CYGNSS) and coarse-resolution SM estimates from the SM active passive (SMAP) mission to estimate SM at 3 km resolution. In practice, the expression for downscaled 3 km SM is derived from the mathematical double-scale SM equations based on the linear assumption between SM and reflectivity, and the CYGNSS signal-to-noise ratio (SNR) data are leveraged to compensate the differences in heterogeneity between the two scales. Experimental validation over 150 in situ sites shows strong consistency between the SMAP/CYGNSS 3 km SM estimates and the in situ measurements, with a median correlation coefficient of 0.806 and a median unbiased root-mean-square error of 0.038 cm3/cm3. The contributions of this work are twofold: 1) introducing the normalized signal-to-noise ratio (NSNR) to account for the deviations in double-scale coefficients, which depend on the variations in vegetation and surface roughness between 36 and 3 km scales, without relying on any ongoing knowledge of ancillary data; and 2) achieving daily 3 km SM estimations at a quasi-global scale, and providing a new way for enhancing the temporal and spatial resolution of SMAP SM.

This current study explores satellite-based soil moisture downscaling approaches and applies them to common passive microwave retrievals. Three variations of a second-order polynomial regression were tested based on the surface temperature/greenness index concept and merged information from higher spatial resolution moderate-resolution imaging spectroradiometer with soil moisture active passive (SMAP) and soil moisture and ocean salinity (SMOS) products to obtain soil moisture estimates at higher resolutions (1 km). Downscaled products were evaluated at the Walnut Gulch Experimental Watershed (WGEW) in southeastern Arizona. Results show slight differences in performance among the three downscaling methods and little improvement between original low-resolution products and downscaled (1 km) products. Spatial analysis over WGEW demonstrates downscaled products were able to decipher small-scale heterogeneities in surface soil moisture, though spatial variability remains low compared to observations with a difference of only 0.06 m3/m3 in spatial standard deviation between observations and the mean between downscaling techniques. Results demonstrate the ability of both SMOS and SMAP to represent soil moisture accurately on the point scale without applying downscaling techniques in the region under study.

<p>Merging data from different instruments is required to construct long time data records of soil moisture (SM). This is the goal of projects such as the ESA Climate Change Initiative (CCI) for SM (Gruber et al., 2019), which uses both active and passive microwave sensors. Currently, the GLDAS v2.1 model is used as reference to re-scale active and passive time series by matching their Cumulative Density Function (CDF) to that of the model. Removing the dependency on models is important, in particular for data assimilation applications into hydrological or climate models, and it has been proposed (Van der Schalie et al., 2018) to use L-band data from one of the two instruments specifically designed to measure SM, ESA Soil Moisture and Ocean Salinity (SMOS) and NASA Soil Moisture Active Passive (SMAP) satellites, as reference to re-scale other time series.<br>To investigate this approach, AMSR-2 SM time series obtained from C1-, C2- and X-band observations using LPRM (Land Parameter Retrieval Model) were re-scaled by CDF-matching (Brocca et al., 2011) using different SMAP and SMOS official (SMAP L2 V005, SMOS L3 V300, SMOS NRT V100&V200) and research (SMOS IC V103) SM products as well as the SMAP and SMOS LPRM v6 SM data used by the ESA CCI. The time series re-scaled using L-band remote sensing data were compared to those re-scaled using GLDAS and were evaluated against in situ measurements at several hundred sites retrieved from the International Soil Moisture Network (Dorigo et al., 2011). The results were analyzed as a function of the land cover class and the Koppen-Geiger climate classification.<br>Overall, AMSR-2 time series re-scaled using SMAP L2, SMAP LPRM and SMOS IC data sets as reference gave the best correlations with respect to in situ measurements, similar to those obtained by the time series re-scaled using GLDAS and slightly better than those of the original AMSR-2 time series. These results imply that different SMAP and SMOS products could actually be used to replace GLDAS as reference for the re-scaling of other sensors time series within the ESA CCI. However, one must bear in mind that this study is limited to the re-scaling of AMSR-2 data at a few hundred sites.<br>For a more detailed assessment of the L-band data set to be used for a global re-scaling, it is necessary to investigate other effects such as the spatial coverage or the time series length. SMAP spatial coverage is better than that of SMOS in regions affected by radio frequency interference. In contrast, the length of SMAP time series can be too short to capture the long term SM variability for climate applications in some regions. The CDF of SMOS time series computed from the date of SMAP launch is significantly different to those of the full length SMOS time series in some regions of the Globe. Possible ways of using a coherent SMAP/SMOS L-band data set will be discussed.</p>

Soil moisture impacts the biosphere–atmosphere exchange of CO2 and CH4 and plays an important role in the terrestrial carbon cycle. A better representation of soil moisture would improve coupled carbon–water dynamics in terrestrial ecosystem models and could potentially improve model estimates of large-scale carbon fluxes and climate feedbacks. Here, we investigate using soil moisture observations from the Soil Moisture Active Passive (SMAP) satellite mission to inform simulated carbon fluxes in the global terrestrial ecosystem model LPJ-wsl. Results suggest that the direct insertion of SMAP reduces the bias in simulated soil moisture at in situ measurement sites by 40%, with a greater improvement at temperate sites. A wavelet analysis between the model and measurements from 26 FLUXNET sites suggests that the assimilated run modestly reduces the bias of simulated carbon fluxes for boreal and subtropical sites at 1–2-month time scales. At regional scales, SMAP soil moisture can improve the estimated responses of CO2 and CH4 fluxes to extreme events such as the 2018 European drought and the 2019 rainfall event in the Sudd (Southern Sudan) wetlands. The simulated improvements to land–surface carbon fluxes using the direct insertion of SMAP are shown across a variety of timescales, which suggests the potential of SMAP soil moisture in improving the model representation of carbon–water coupling.

In large river basins where in situ data were limited or absent, satellite-based soil moisture estimates can be used to supplement ground measurements for land and water resource management solutions. Consistent soil moisture estimation can aid in monitoring droughts, forecasting floods, monitoring crop productivity, and assisting weather forecasting. Satellite-based soil moisture estimates are readily available at the global scale but are provided at spatial scales that are relatively coarse for many hydrological modeling and decision-making purposes. Soil moisture data are obtained from NASA’s soil moisture active passive (SMAP) mission radiometer as an interpolated product at 9 km gridded resolution. This study implements a soil moisture downscaling algorithm that was developed based on the relationship between daily temperature change and average soil moisture under varying vegetation conditions. It applies a look-up table using global land data assimilation system (GLDAS) soil moisture and surface temperature data, and advanced very high resolution radiometer (AVHRR) and moderate resolution imaging spectroradiometer (MODIS) normalized difference vegetation index (NDVI) and land surface temperature (LST). MODIS LST and NDVI are used to obtain downscaled soil moisture estimates. These estimates are then used to enhance the spatial resolution of soil moisture estimates from SMAP 9 km to 1 km. Soil moisture estimates at 1 km resolution are able to provide detailed information on the spatial distribution and pattern over the regions being analyzed. Higher resolution soil moisture data are needed for practical applications and modelling in large watersheds with limited in situ data, like in the Lower Mekong River Basin (LMB) in Southeast Asia. The 1 km soil moisture estimates can be applied directly to improve flood prediction and assessment as well as drought monitoring and agricultural productivity predictions for large river basins.

Exploring SMAP and OCO-2 observations to monitor soil moisture control on photosynthetic activity of global drylands and croplands

Given the importance of soil moisture for hydrological applications, such as weather and flood forecasting, passive microwave remote sensing is a promising approach for retrieving soil moisture due to its high sensitivity to near-surface soil moisture, applicability to all weather conditions, direct relationship with the soil dielectric constant, and reduced effects from vegetation and roughness. However, passive microwave (radiometer) observations suffer from being relatively low spatial resolution, on the order of 36 km. It is proposed that this scale issue may be overcome by using active microwave (radar) observations, which have much higher resolution when using Synthetic Aperture Radar (SAR) techniques (<3km), and this is the approach being taken by NASA's Soil Moisture Active Passive (SMAP) mission, with a scheduled launch in late 2014. The rationale behind SMAP is to use the synergy between active and passive observations in a downscaling approach to overcome the individual limitations of each observation type, and ultimately provide a merged soil moisture data set at intermediate resolution (~9 km). The objective of this study is to test the proposed baseline downscaling approach for the SMAP mission using airborne data, thus assessing its viability for future application to SMAP data. The approach is based on the hypothesis of a near-linear relationship between radiometer brightness temperature (Tb) and SAR backscatter (σ), and has thus far received very limited testing. The experimental dataset used in this study was collected during the Soil Moisture Active Passive Experiment (SMAPEx) field campaigns over a study site in south-eastern Australia approximately 38km × 36km in size, equivalent to a SMAP radiometer pixel. This research focuses on the brightness temperature downscaling algorithm; according to the SMAP Algorithm Theoretical Basis Documents these downscaled brightness temperatures will subsequently be converted to soil moisture at fine resolution through the traditional passive microwave retrieval algorithm. The baseline downscaling algorithm was applied to high resolution data from SMAPEx, which include 1km resolution brightness temperature collected by the Polarimetric L-band Multibeam Radiometer (PLMR) and ~ 10m resolution backscatters collected by the Polarimetric L-band Imaging Synthetic aperture radar (PLIS). To minimize noise in the radar data and to approximate the SMAP radiometer/radar pixel ratios (36km Tb to 9km resolution using 3km σ) the PLIS data were aggregated to 250m resolution, so as to downscale 1km Tb to 250m resolution, thus keeping the same ratio of radiometer/SAR spatial resolution as the SMAP mission. Results showed that the Root-Mean-Square Error (RMSE) in Tb downscaled at 100m resolution was around 10K at h-polarization and 8K at v-polarization over a cropping area. This RMSE was reduced to 9K and 7K respectively when downscaling to 250m resolution, due to a decreased spatial heterogeneity during averaging. It was also noted that results at v-polarization were slightly better than those at h-polarization, since the backscatter is more linearly related to Tb at v-polarization than Tb at h-polarization. The accuracy of the downscaling over grassland sites was improved by approximately 3K with respect to the cropping area. This was attributed to the more heterogeneous conditions in cropping areas, compared to the relatively uniform conditions in the grassland area. However, one limitation of this study was the availability of only three days of data for estimating the linearity between radar and radiometer observations.

Long and consistent soil moisture time series at adequate spatial resolution are key to foster the application of soil moisture observations and remotely-sensed products in climate and numerical weather prediction models. The two L-band soil moisture satellite missions SMAP (Soil Moisture Active Passive) and SMOS (Soil Moisture and Ocean Salinity) are able to provide soil moisture estimates on global scales and in kilometer accuracy. However, the SMOS data record has an appropriate length of 7.5 years since late 2009, but with a coarse resolution of ∼25 km only. In contrast, a spatially-enhanced SMAP product is available at a higher resolution of 9 km, but for a shorter time period (since March 2015 only). Being the fundamental observable from passive microwave sensors, reliable brightness temperatures (Tbs) are a mandatory precondition for satellite-based soil moisture products. We therefore develop, evaluate and apply a copula-based data fusion approach for combining SMAP Enhanced (SMAP_E) and SMOS brightness Temperature (Tb) data. The approach exploits both linear and non-linear dependencies between the two satellite-based Tb products and allows one to generate conditional SMAP_E-like random samples during the pre-SMAP period. Our resulting global Copula-combined SMOS-SMAP_E (CoSMOP) Tbs are statistically consistent with SMAP_E brightness temperatures, have a spatial resolution of 9 km and cover the period from 2010 to 2018.&#13;\n&#13;\nCurrently, the data contains only horizontally polarised retrievals from SMOS's and SMAP's ascending and descending nodes, respectively.

The Soil Moisture Active Passive (SMAP) mission is designed to acquire L-band radiometer measurements for the estimation of soil moisture with 0.04 m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> /m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> volumetric accuracy in the top 5 cm for vegetation with water content of less than 5 kg/m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> . In regions near the coast or near inland bodies of water, the signal measured by the SMAP radiometer contains emissions from land and water, resulting in errors in the soil moisture estimation. In this article, the effort to extract the brightness temperature (TB) according to the land fraction or water fraction (depending on the center of the footprint location) from the affected SMAP measurements was addressed. A single pixel correction algorithm was applied and its performance was evaluated over simulated data. A data-driven approach for the estimation of land and water TB for data correction was developed. The correction algorithm was then applied to real data and its performance was assessed over the SMAP soil moisture retrievals. We showed that the single pixel algorithm is an effective and computationally efficient algorithm for removing land or water TB contamination from the SMAP data.

Root zone soil moisture (RZSM) is a vital variable for agricultural production, water resource management and runoff prediction. Satellites provide large-scale and long-term near-surface soil moisture retrievals, which can be used to estimate RZSM through various methods. In this study, we tested the utility of an exponential filter (ExpF) using in situ soil moisture by optimizing the optimal characteristic time length T_opt for different soil depths. Furthermore, the parameter analysis showed that T_opt correlated negatively with precipitation and had no significant correlation with selected soil properties. Two approaches were taken to obtain T_opt: (1) optimization of the Nash–Sutcliffe efficiency coefficient (NSE); (2) calculation based on annual average precipitation. The precipitation-based T_pre outperformed the station-specific T_opt and stations-averaged T_opt. To apply the ExpF on grid scale, the precipitation-based T_pre considering spatial variability was adopted in the ExpF to obtain RZSM from a new soil moisture dataset RF_SMAP_L3_P (Random Forest Soil Moisture Active Passive_L3_Passive) continuous in time and space over Huai River Basin. Finally, the performance of RF_SMAP_L3_P RZSM (0–100 cm) was evaluated using in situ measurements and compared with mainstream products, for instance, Soil Moisture Active Passive (SMAP) and Soil Moisture and Ocean Salinity Level 4 (SMOS L4) RZSM. The results indicated that RF_SMAP_L3_P RZSM could captured the temporal variation of measured RZSM best with R value of 0.586, followed by SMAP L4, which had the lowest bias value of 0.03, and SMOS L4 significantly underestimated the measured RZSM with bias value of −0.048 in the basin. Higher accuracy of RF_SMAP_L3_P RZSM was found in the flood period compared with the non-flood period, which indicates a better application for ExpF in wetter weather conditions.

Global-scale evaluation of SMAP, SMOS and ASCAT soil moisture products using triple collocation