High Resolution Surface Soil Moisture Research Articles

A novel land surface temperature reconstruction method and its application for downscaling surface soil moisture with machine learning

Downscaling of soil moisture data is important for high resolution hydrological modeling. Most downscaling studies in the literature have used spatially discontinuous land surface temperature (LST) maps as the main auxiliary parameter, which limits the creation of continuous soil moisture maps. The number of studies on soil moisture downscaling with machine learning that use gapless LST maps is limited. With this motivation, a hybrid reconstruction method has been proposed in this study to practically obtain continuous LST maps, which are then used to produce high resolution surface soil moisture (SSM) datasets. The proposed method is shown to have high mean performance with R2 and RMSE values of 0.94 and 1.84°K, respectively, for the period between 2019 and 2022. The developed reconstructed LST maps were then used to downscale original 9 km spatial resolution soil moisture datasets of SMAP L3 and SMAP L4 with Random Forest (RF) machine learning algorithm. The RF model were run with four different rainfall datasets, and the MSWEP rainfall dataset was found to produce the best results. The use of antecedent rainfall values as input variables in machine learning models has been shown to improve the performance of the models R2 0.76 to 0.93. The accuracy of the downscaled data was later evaluated for Western Anatolia Basins (WAB) in Türkiye with 31 in-situ stations. The downscaled SMAP L4 had good average statistical indicators R (0.815 ± 0.1), RMSE (0.09 ± 0.047 cm3/cm3), and ubRMSE (0.058 ± 0.025 cm3/cm3). Downscaled SMAP L3 was also validated with in-situ observations with satisfactory R (0.79 ± 0.074), RMSE (0.09 ± 0.043 cm3/cm3), and ubRMSE (0.06 ± 0.026 cm3/cm3) statistics. Furthermore, the performance of the downscaled SMAP L3 was also cross validated with SMAP + Sentinel 1 (L2) dataset between 2019 and 2022. The mean statistics of R (0.761 ± 0.11) and Root Mean Squared Difference (RMSD) (0.05 ± 0.014 cm3/cm3) between downscaled SMAP L3 and L2 data revealed that the new reconstruction method of LST used in the RF model for downscaling of soil moisture performed well to obtain high resolution soil moisture datasets. The proposed technique also overcame the difficulties associated with coastal regions where data was masked for quality considerations, by not only enhancing overall spatial resolution but also filling these data gaps and giving a complete SSM coverage.

Soil moisture estimates at 1 km resolution making a synergistic use of Sentinel data

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.

Assessment of High-Resolution Surface Soil Moisture Products over the Qinghai–Tibet Plateau for 2009–2017

The surface soil moisture over the Qinghai–Tibet Plateau (QTP) has an important impact on the weather and climate of East Asia. Under climate warming, the imbalance of solid–liquid water of QTP has become a research hotspot, but the surface soil moisture dynamics over QTP are not clear owing to the lack of precise measurements over a large scale. In this paper, the quality of gridded surface soil moisture products including CSSPv2 high-resolution (6 km) simulation, ESA CCI satellite retrieval, ERA5 land-atmosphere coupled reanalysis, and GLDAS2.1 land reanalysis products (Noah, Catchment, VIC) is analyzed over QTP by comparison with the in situ measurements at 140 stations during 2009–2017. We find that the CSSPv2 product shows a higher correlation than the global satellite and reanalysis products, with correlation increased by 7.7%–115.6%. The root mean squared error of the CSSPv2 product is lower than that of other products, with the error decreased by 13.4%–46.3%. The triple collocation analysis using high-resolution simulation, global reanalysis, and satellite retrieval products over the entire plateau shows that the error of CSSPv2 is the lowest, followed by ESA CCI, while ERA5 is the highest. The soil moisture products of ESA CCI, ERA5, and CSSPv2 all show an increasing trend from April to September of 2009 to 2017, with wetting in the west and drying in the east. This study indicates that the CSSPv2 high-resolution surface soil moisture product has better performance over QTP than other global products, and the global satellite and reanalysis products may overestimate the surface soil moisture dynamics.

SMPD: a soil moisture-based precipitation downscaling method for high-resolution daily satellite precipitation estimation

Abstract. As a key component in the water and energy cycle, estimates of precipitation with high resolution and accuracy is of great significance for hydrological, meteorological, and ecological studies. However, current satellite-based precipitation products have a coarse spatial resolution (from 10 to 50 km) not meeting the needs of several applications (e.g., flash floods and landslides). The implementation of spatial downscaling methods can be a suitable approach to overcome this shortcoming. In this study, we developed a soil moisture-based precipitation downscaling (SMPD) method for spatially downscaling the integrated multisatellite retrievals for global precipitation measurement (IMERG) V06B daily precipitation product over a complex topographic and climatic area in southwestern Europe (Iberian Peninsula) in the period 2016–2018. By exploiting the soil-water balance equation, high-resolution surface soil moisture (SSM) and normalized difference vegetation index (NDVI) products were used as auxiliary variables. The spatial resolution of the IMERG daily precipitation product was downscaled from 10 to 1 km. An evaluation using 1027 rain gauge stations highlighted the good performance of the downscaled 1 km IMERG product compared to the original 10 km product, with a correlation coefficient of 0.61, root mean square error (RMSE) of 4.83 mm and a relative bias of 5 %. Meanwhile, the 1 km downscaled results can also capture the typical temporal and spatial variation behaviors of precipitation in the study area during dry and wet seasons. Overall, the SMPD method greatly improves the spatial details of the original 10 km IMERG product also with a slight enhancement of accuracy. It shows good potential to be applied for the development of high-quality and high-resolution precipitation products in any region of interest.

An Exponential Filter Model-Based Root-Zone Soil Moisture Estimation Methodology from Multiple Datasets

Modern smart agriculture initiative presents more requests for soil moisture (SM) monitoring over large agricultural areas. Remote sensing techniques facilitate high-resolution surface SM (SSM) estimation at a large scale but lack root zone SM (RZSM) information. Establishing the deduction method of RZSM from the SSM has long been the focus of most attention. Data assimilation methods are promising techniques for RZSM estimation, developing numerous assimilated reanalysis datasets, e.g., ERA5 and the latest Soil Moisture Active and Passive (SMAP) L4 SM product. However, data latency and large computation during data collecting and processing often inhibits further applications. This work proposes a rapid estimation scheme for estimating RZSM with short latency and small computations, based on the Exponential Filter (EF) method. The EF model with single parameter T was firstly calibrated and validated using the SSM and RZSM of ERA5 reanalysis dataset, obtaining the optimum parameter T map for each grid. Then, the fast-updating SMAP L3 SSM product together with the scale-matched optimum T were adopted as inputs into the EF model to retrieve RZSM estimation of each grid. Specifically, such estimation scheme was tested over the central and eastern agricultural areas of China, using a dense monitoring network of 796 SM observation sites, which contains various land uses, as well as meteorological and hydrological conditions. The calibrated optimum parameter T presented an increasing trend with good physical explanations. Furthermore, all the estimated RZSMs were found to have good performances on capturing the temporal-spatial variations of RZSM and well reflecting seasonal RZSM changes. Overall, such an estimation scheme was proven to be a desirable alternative for estimating RZSM over large agricultural areas.

Estimating High-Resolution Soil Moisture Over Mountainous Regions Using Remotely-Sensed Multispectral and Topographic Data

A surface soil moisture (SM) condition at high spatiotemportal resolutions is required by regional Earth system applications. Here, we mapped daily 1-km SM in the Babao River Basin in the northwest of China during the summers from 2013 to 2015 using a random forest (RF) method by merging SM information retrieved from <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in situ</i> measurements, optical/thermal remote sensing, and topographical indices. Relative importance analysis was used to determine the optimal predictors for estimating high-resolution SM. A specific RF model (RF <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">VI</sub> +sup) was constructed using the optimal predictors including remote sensing albedo, apparent thermal inertia (ATI), normalized difference vegetation index, normalized difference infrared index 5, soil adjusted vegetation index, and topographical indices (aspect and elevation). The RF <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">VI+sup</sub> also accounted for missing observations of the thermal index (e.g., ATI) over the mountainous regions. In the comparison between the SM estimates using the new RF <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">VI+sup</sub> model and other RF models, the spatial coverage of available estimates increased from 14% to 64% over the study region, the correlation coefficient values were improved to 0.75, the unbiased root-mean-squared difference values decreased to 0.032 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> . Thus, the proposed RF method provided accurate SM estimates with high spatiotemporal resolution over the mountainous regions, by merging multiresource datasets from <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in situ</i> measurements, remotely-sensed, and topographical indices.

A Semi-Physical Approach for Downscaling Satellite Soil Moisture Data in a Typical Cold Alpine Area, Northwest China

Microwave remote sensing techniques provide a direct measurement of surface soil moisture (SM), with advantages for all-weather observations and solid physics. However, most satellite microwave soil moisture products fail to meet the requirements of land surface studies for high-resolution surface soil moisture data due to their coarse spatial resolutions. Although many approaches have been proposed to downscale the spatial resolution of satellite soil moisture products, most of them have been tested in flat areas where the surface is relatively homogeneous. Thus, those established approaches are often inapplicable for downscaling in cold alpine areas with complex terrain where multiple factors control the variations in surface soil moisture. In this work, we re-inferred and verified the mathematical assumption behind a semi-physical approach for downscaling satellite soil moisture data and extended this approach for cold alpine areas. Instead of directly deriving SM from proxy variables, this approach relies on a relationship between two standardized variables of SM and apparent thermal inertia (ATI), in which the sub grid standard deviation for SM is estimated by a physical hydraulic model taking soil texture data as input. The approach was applied to downscale the soil moisture active passive (SMAP) daily data in a typical cold alpine basin, i.e., the Babao River basin located in the Qilian Mountains of Northwest China. We observed good linearity between the computed ATI and SM observations on most wireless sensor network sites installed in the study basin, which justifies the underlying assumption. The sub grid standard deviations for the SMAP grid estimated through the Mualem-van Genuchten model can broadly represent the real characteristics. The downscaled 1-km resolution results correlated well with the in-situ SM observations, with an average correlation coefficient of 0.74 and a small root mean square error (0.096 cm3/cm3). The downscaled results show more and consistent textural details than the original SMAP data. After removal of biases in the original SMAP data even higher agreements with the observations can be achieved. These results demonstrate the adequacy of the proposed semi-physical approach for downscaling satellite soil moisture data in cold alpine areas, and the resultant fine-resolution data can serve as useful databases for land surface and hydrological studies in those areas.

Temporal Calibration of an Evaporation-Based Spatial Disaggregation Method of SMOS Soil Moisture Data

The resolution of current satellite surface soil moisture (SM) estimates is very low, of tens of kilometers, which proves to be insufficient for various agricultural and hydrological applications. Amongst the existing downscaling approaches of remotely sensed SM, DISPATCH (DISaggregation based on a Physical And Theoretical scale CHange) improves the resolution of SMOS (Soil Moisture and Ocean Salinity) soil moisture data using soil evaporative efficiency (SEE) estimates at high resolution (HR) and a SEE(SM) model implemented at low resolution (LR). Defined as the ratio of actual to potential soil evaporation, SEE can be derived from the remotely sensed land surface temperature (LST) and normalized difference vegetation index (NDVI). The current version of DISPATCH uses a linear SEE(SM) model. This study aims at improving the SEE(SM) model and testing different calibration strategies, to ultimately have more robust and better downscaled SM products. A nonlinear SEE(SM) model is introduced and its influence on the derived HR SM products is studied over a range of conditions. Each model, linear and nonlinear, is calibrated from remote sensing data on a daily and a multi-date basis. The approaches were tested over two mixed dry and irrigated areas in Catalonia, Spain, and over one dry area in Morocco. When using the linear model, better statistical results were generally obtained using a daily calibration (current version of DISPATCH), most notably over one Spanish site. However, the best results were systematically obtained for an annually calibrated nonlinear model, in terms of all metrics considered: correlation coefficient, slope of the linear regression, bias, unbiased root mean square error. In particular, when using the annually calibrated nonlinear SEE (SM) model, the temporal slope of the linear regression between disaggregated and in situ soil moisture increased to 1.16 and 0.75 for one Spanish site and for the Moroccan site (as opposed to 0.44 and 0.58, respectively, when using the linear model with a daily calibration). The temporal correlation coefficient increased to 0.47 and 0.54 over the Spanish sites (as opposed to 0.18 and 0.27, respectively, when using the linear model with a daily calibration). Those contrasted results indicate compensation effects between the model type and the calibration strategy. Taking into account studies that report the strong nonlinear behavior of the SEE with respect to SM, the introduction of the nonlinear SEE(SM) model in DISPATCH, combined with a multi-date calibration, is proven to perform significantly better under various conditions, leading to more robust disaggregated SM products. The SEE modeling based on the nonlinear SM model, with a multi-date calibration, could be integrated into the CATDS—Centre Aval de Traitement des Données SMOS as a future product, as well as into existing evapotranspiration models, which are based on a combination of thermal and microwave data.

New Downscaling Approach Using ESA CCI SM Products for Obtaining High Resolution Surface Soil Moisture

ESA CCI SM products have provided remotely-sensed surface soil moisture (SSM) content with the best spatial and temporal coverage thus far, although its output spatial resolution of 25 km is too coarse for many regional and local applications. The downscaling methodology presented in this paper improves ESA CCI SM spatial resolution to 1 km using two-step approach. The first step is used as a data engineering tool and its output is used as an input for the Random forest model in the second step. In addition to improvements in terms of spatial resolution, the approach also considers the problem of data gaps. The filling of these gaps is the initial step of the procedure, which in the end produces a continuous product in both temporal and spatial domains. The methodology uses combined active and passive ESA CCI SM products in addition to in situ soil moisture observations and the set of auxiliary downscaling predictors. The research tested several variants of Random forest models to determine the best combination of ESA CCI SM products. The conclusion is that synergic use of all ESA CCI SM products together with the auxiliary datasets in the downscaling procedure provides better results than using just one type of ESA CCI SM product alone. The methodology was applied for obtaining SSM maps for the area of California, USA during 2016. The accuracy of tested models was validated using five-fold cross-validation against in situ data and the best variation of model achieved RMSE, R2 and MAE of 0.0518 m3/m3, 0.7312 and 0.0374 m3/m3, respectively. The methodology proved to be useful for generating high-resolution SSM products, although additional improvements are necessary.

The SMAP and Copernicus Sentinel 1A/B microwave active-passive high resolution surface soil moisture product

Soil Moisture Active Passive (SMAP) mission of NASA was launched in January 2015. Currently, SMAP has an L-band radiometer and a defunct L-band radar with a rotating 6-m mesh reflector antenna. On July 7th, 2015, the SMAP radar malfunctioned and became inoperable. Consequently, the production of high-resolution active-passive soil moisture product got hampered, and only ~2.5 months (April 15th, 2015 to July 7th, 2015) of data remain available. Therefore, during the SMAP post-radar phase, many ways were examined to restart the high-resolution soil moisture product generation of the SMAP mission. One of the feasible approaches was to substitute the SMAP radar with other available SAR data. Sentinel-1A/Sentinel-1B SAR data was found most suitable for combining with the SMAP radiometer data because of its nearly similar orbit configuration that allows overlapping of their swaths with a minimal time difference, a key feature/requirement for the SMAP active-passive algorithm. The Sentinel interferometric wide swath (IW) mode acquisition also provides the co-polarized and cross-polarized observations required for the SMAP active-passive algorithm. However, some differences do exist between the SMAP and Sentinel SAR data. They are mainly: 1) Sentinel has a C-band SAR whereas SMAP operates at L-band; 2) Sentinel has multiple incidence angles within its swath, and SMAP has one single incidence angle; and 3) Sentinel 1A/B Interferometric Wide (IW) swath width is ~250 km as compared to SMAP with 1000 km swath width. On any given day, the narrow swath width of the Sentinel observations significantly reduces the overlap spatial coverage between SMAP and Sentinel as compared to the original SMAP radar and radiometer swath coverage. Hence, the temporal resolution (revisit interval) suffers due to narrow overlapped swath width and degrades from 3 days to 12 days. One advantage of using very high-resolution resolution Sentinel-1A/Sentinel-1B data in the SMAP active-passive algorithm is the potential of obtaining the disaggregated brightness temperature and thus soil moisture at a much finer spatial resolution of 3 km and 1 km at global extent. The assessment of high-resolution product at 3 km and 1 km using the soil moisture calibration and validations sites shows reasonable accuracy of ~0.05 m3/m3. The SMAP-Sentinel1 active-passive high-resolution product is now available to the public (new version released in October 2018) through NSIDC (NASA DAAC). The duration of this product is from April 2015 to current date.

Data-driven stochastic model for basin and sub-grid variability of SMAP satellite soil moisture

This study focused on the utility of coarse surface soil moisture observations for applications that require high resolution surface soil moisture information. This was accomplished by quantifying the information content of average soil moisture for three different spatial scales of 81 km2, 790 km2, and 4400 km2. In situ point observations of soil moisture from 31 stations in Iowa were used to develop a spatial stochastic model that assumes hillslope-scale model parameters are independent. Soil moisture dry-downs and wetting regimes were analyzed using rain gauge and soil moisture sensor data. The statistical nature of dry-downs were parameterized using a power-law decay, and soil moisture increases due to rainfall were parameterized using a non-dimensional logistic curve that is a function of soil moisture deficit. The resulting stochastic model is used to quantify the magnitude of the standard deviation, σ(θ), and skewness G(θ) as a function of the areal average. We show that the greatest information content (small spatial standard deviation) of average observation corresponded to values near the minimum or the maximum soil moisture with σ(θ)<5%, while average observations for intermediate soil moisture values had the lowest information content with σ(θ)>20%. The differences in information content as a function of the areal average were consistent with the statistical nature of soil moisture that can be interpreted as small range bounded variable. However, this study provides quantitative estimates for the magnitude of the sub-grid and basin scale variability, documenting the utility for applications that require high resolution information. These results form the basis for the investigation of spatial runoff production in response to rainfall and to inform plot scale agriculture applications.

Effect of vegetation index choice on soil moisture retrievals via the synergistic use of synthetic aperture radar and optical remote sensing

The recent launch of the Sentinel-1 A and Sentinel-1B synthetic aperture radar (SAR) satellite constellation has provided high-quality SAR data with fine spatial and temporal sampling characterizations (6˜12 revisit days at 10 m spatial resolution). When combined with high-resolution optical remote sensing, this data can potentially be used for high-resolution soil moisture retrieval over vegetated areas. However, the suitability of different vegetation index (VI) types for the parameterization of vegetation water content in SAR vegetation scattering models requires further investigation. In this study, the widely-used physical-based Advanced Integral Equation Model (AIEM) is coupled with the Water Cloud Model (WCM) for the retrieval of field-scale soil moisture. Three different VIs (NDVI, EVI, and LAI) produced by two different satellite sensors (Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat) are selected to examine their impact on the parameterization of vegetation water content, and subsequently, on soil moisture retrieval accuracy. Results indicate that, despite the different sensitivity of estimated surface roughness parameters to various VIs (i.e., this sensitivity is highest when utilizing MODIS EVI and lowest in the LAI-based model), the optimum roughness parameters derived from each VI exhibit no discernible difference. Consequently, the soil moisture retrieval accuracies show no noticeable sensitivity to the choice of a particular VI. Generally, meadow and grassland sites with small differences in VI-derived roughness parameters exhibit good performance in soil moisture estimation. With respect to the relative components in the coupled model, the vegetative contribution to the scattering signal exceeds that of soil at VI about 0.6∼0.8 [-] in NDVI-based models and 0.4∼0.6 [-] in EVI-based models. This study provides insight into the proper selection of vegetation indices during the use of SAR and optical imagery for the retrieval of high-resolution surface soil moisture.

Mapping Soil Moisture at a High Resolution over Mountainous Regions by Integrating In Situ Measurements, Topography Data, and MODIS Land Surface Temperatures

Hydro-agricultural applications often require surface soil moisture (SM) information at high spatial resolutions. In this study, daily spatial patterns of SM at a spatial resolution of 1 km over the Babao River Basin in northwestern China were mapped using a Bayesian-based upscaling algorithm, which upscaled point-scale measurements to the grid-scale (1 km) by retrieving SM information using Moderate Resolution Imaging Spectroradiometer (MODIS)-derived land surface temperature (LST) and topography data (including aspect and elevation data) and in situ measurements from a wireless sensor network (WSN). First, the time series of pixel-scale (1 km) representative SM information was retrieved from in situ measurements of SM, topography data, and LST. Second, Bayesian linear regression was used to calibrate the relationship between the representative SM and the WSN measurements. Last, the calibrated relationship was used to upscale a network of in situ measured SM to map spatially continuous SM at a high resolution. The upscaled SM data were evaluated against ground-based SM measurements with satisfactory accuracy—the overall correlation coefficient (r), slope, and unbiased root mean square difference (ubRMSD) values were 0.82, 0.61, and 0.025 m3/m3, respectively. Moreover, when accounting for topography, the proposed upscaling algorithm outperformed the algorithm based only on SM derived from LST (r = 0.80, slope = 0.31, and ubRMSD = 0.033 m3/m3). Notably, the proposed upscaling algorithm was able to capture the dynamics of SM under extreme dry and wet conditions. In conclusion, the proposed upscaled method can provide accurate high-resolution SM estimates for hydro-agricultural applications.

A Spatially Consistent Downscaling Approach for SMOS Using an Adaptive Moving Window

The European Space Agency (ESA)'s Soil Moisture and Ocean Salinity (SMOS) is the first spaceborne mission using L-band radiometry to monitor the Earth's global surface soil moisture (SM). After more than 7 years in orbit, many studies have contributed to improve the quality and applicability of SMOS-derived SM maps. In this research, a novel downscaling algorithm for SMOS is proposed to obtain high-resolution (HR) SM maps at 1 km (L4), from the ∼40 km native resolution of the instrument. This algorithm introduces the concept of a shape adaptive moving window as an improvement of the current semi-empirical downscaling approach at SMOS Barcelona Expert Center, based on the “universal triangle”. Its inputs are as follows: the SMOS SM (L3 at ∼40 km spatial resolution), the vertical and the horizontal SMOS brightness temperatures (L1C at ∼40 km), and the HR normalized difference vegetation index and land surface temperature from optical-based sensors. The proposed method has the following advantages: HR SM maps are obtained while maintaining the dynamic range from the original L3 product; energy is conserved, because differences between aggregated L4 and L3 SM maps are negligible; and it can be applied to continental areas, even when they integrate different climates. A comparison of SMOS L3 and L4 products with in situ data for networks allocated in Spain and Australia shows good agreement in terms of correlation and root mean square error. The proposed method is shown to capture 1-km SM spatial variability while preserving the quality of SMOS at its native resolution.

Spatial Downscaling of SMAP Soil Moisture Using MODIS Land Surface Temperature and NDVI During SMAPVEX15

The Soil Moisture Active Passive (SMAP) mission provides a global surface soil moisture (SM) product at 36-km resolution from its L-band radiometer. While the coarse resolution is satisfactory to many applications, there are also a lot of applications which would benefit from a higher resolution SM product. The SMAP radiometer-based SM product was downscaled to 1 km using Moderate Resolution Imaging Spectroradiometer (MODIS) data and validated against airborne data from the Passive Active L-band System instrument. The downscaling approach uses MODIS land surface temperature and normalized difference vegetation index to construct soil evaporative efficiency, which is used to downscale the SMAP SM. The algorithm was applied to one SMAP pixel during the SMAP Validation Experiment 2015 (SMAPVEX15) in a semiarid study area for validation of the approach. SMAPVEX15 offers a unique data set for testing SM downscaling algorithms. The results indicated reasonable skill (root-mean-square difference of 0.053 m3/m3 for 1-km resolution and 0.037 m3/m3 for 3-km resolution) in resolving high-resolution SM features within the coarse-scale pixel. The success benefits from the fact that the surface temperature in this region is controlled by soil evaporation, the topographical variation within the chosen pixel area is relatively moderate, and the vegetation density is relatively low over most parts of the pixel. The analysis showed that the combination of the SMAP and MODIS data under these conditions can result in a high-resolution SM product with an accuracy suitable for many applications.

Retrieving High-Resolution Surface Soil Moisture by Downscaling AMSR-E Brightness Temperature Using MODIS LST and NDVI Data

A method to retrieve soil moisture at high spatial resolution is presented in this paper. The method is based on soil moisture retrieval with passive brightness temperature. The method of retrieving land surface temperature with passive microwave is combined with the relationship between the microwave polarization difference index (MPDI) and normalized difference of vegetation index (NDVI) to obtain high-resolution microwave brightness temperature and soil moisture. Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) 18.7-GHz brightness temperature at 25-km resolution is downscaled to 1-km using high-resolution MODIS visible/infrared (VIS/IR) data. High-resolution soil moisture is retrieved with the downscaled microwave brightness temperature using a single-channel algorithm (SCA) and the Qp model to deal with the influence of roughness. The method is applied to an area in northwest of China. The downscaled high-resolution soil moisture is tested with ground data collected at three sites within the Maqu monitoring network from July 1, 2008 to June 30, 2009. The trend of the time series of the downscaled soil moisture is similar to the ground measurements during this period with root mean-square error (RMSE) less than 0.12. The results show that the method is more suitable to moderate to drier soil conditions with bare surface or covered by sparse vegetation.

A method for retrieving high-resolution surface soil moisture from hydros L-band radiometer and Radar observations

NASA's Earth System Science Pathfinder Hydrospheric States (Hydros) mission will provide the first global scale space-borne observations of Earth's soil moisture using both L-band microwave radiometer and radar technologies. In preparation for the Hydros mission, an observation system simulation experiment (OSSE) has been conducted. As a part of this OSSE, the potential for retrieving useful surface soil moisture at spatial resolutions of 9 and 3 km was explored. The approach involved optimally merging relatively accurate 36-km radiometer brightness temperature and relatively noisy 3-km radar backscatter cross section observations using a Bayesian method. Based on the Hydros OSSE data sets with low and high noises added to the simulated observations or model parameters, the Bayesian method performed better than direct inversion of either the brightness temperature or radar backscatter observations alone. The root-mean-square errors of 9-km soil moisture retrievals from the Bayesian merging method were reduced by 0.5 %vol/vol and 1.4 %vol/vol from the errors of direct radar inversions for the entire OSSE domain of all 34 consecutive days for the low and high noise data sets, respectively. Improvement in soil moisture estimates using the Bayesian merging method over the direct inversions of radar or radiometer data were even more significant for soil moisture retrieval at 3-km resolution. However, to address the representativeness of these results at the global and multiyear scales, further performance comparison studies are needed, particularly with actual field data.