An Approach for Downscaling SMAP Soil Moisture by Combining Sentinel-1 SAR and MODIS Data

Estimation of 100 m root zone soil moisture by downscaling 1 km soil water index with machine learning and multiple geodata

Root zone soil moisture (RZSM) is crucial for agricultural water management and land surface processes. The 1 km soil water index (SWI) dataset from Copernicus Global Land services, with eight fixed characteristic time lengths (T), requires root zone depth optimization (Topt) and is limited in use due to its low spatial resolution. To estimate RZSM at 100-m resolution, we integrate the depth specificity of SWI and employed random forest (RF) downscaling. Topographic synthetic aperture radar (SAR) and optical datasets were utilized to develop three RF models (RF1: SAR, RF2: optical, RF3: SAR + optical). At the DEMMIN experimental site in northeastern Germany, Topt (in days) varies from 20 to 60 for depths of 10 to 30 cm, increasing to 100 for 40-60 cm. RF3 outperformed other models with 1 km test data. Following residual correction, all high-resolution predictions exhibited strong spatial accuracy (R ≥ 0.94). Both products (1 km and 100 m) agreed well with observed RZSM during summer but overestimated in winter. Mean R between observed RZSM and 1 km (100 m; RF1, RF2, and RF3) SWI ranges from 0.74 (0.67, 0.76, and 0.68) to 0.90 (0.88, 0.81, and 0.82), with the lowest and highest R achieved at 10 cm and 30 cm depths, respectively. The average RMSE using 1 km (100 m; RF1, RF2, and RF3) SWI increased from 2.20 Vol.% (2.28, 2.28, and 2.35) at 30 cm to 3.40 Vol.% (3.50, 3.70, and 3.60) at 60 cm. These negligible accuracy differences underpin the potential of the proposed method to estimate RZSM for precise local applications, e.g., irrigation management.

Assimilation of SMAP Products for Improving Streamflow Simulations over Tropical Climate Region—Is Spatial Information More Important Than Temporal Information?

Streamflow is one of the key variables in the hydrological cycle. Simulation and forecasting of streamflow are challenging tasks for hydrologists, especially in sparsely gauged areas. Coarse spatial resolution remote sensing soil moisture products (equal to or larger than 9 km) are often assimilated into hydrological models to improve streamflow simulation in large catchments. This study uses the Ensemble Kalman Filter (EnKF) technique to assimilate SMAP soil moisture products at the coarse spatial resolution of 9 km (SMAP 9 km), and downscaled SMAP soil moisture product at the higher spatial resolution of 1 km (SMAP 1 km), into the Soil and Water Assessment Tool (SWAT) to investigate the usefulness of different spatial and temporal resolutions of remotely sensed soil moisture products in streamflow simulation and forecasting. The experiment was set up for eight catchments across the tropical climate of Vietnam, with varying catchment areas from 267 to 6430 km2 during the period 2017–2019. We comprehensively evaluated the EnKF-based SWAT model in simulating streamflow at low, average, and high flow. Our results indicated that high-spatial resolution of downscaled SMAP 1 km is more beneficial in the data assimilation framework in aiding the accuracy of streamflow simulation, as compared to that of SMAP 9 km, especially for the small catchments. Our analysis on the impact of observation resolution also indicates that the improvement in the streamflow simulation with data assimilation is more significant at catchments where downscaled SMAP 1 km has fewer missing observations. This study is helpful for adding more understanding of performances of soil moisture data assimilation based hydrological modelling over the tropical climate region, and exhibits the potential use of remote sensing data assimilation in hydrology.

Enhancing Spatial Resolution of SMAP Soil Moisture Products through Spatial Downscaling over a Large Watershed: A Case Study for the Susquehanna River Basin in the Northeastern United States

Soil moisture (SM) with a high spatial resolution plays a paramount role in many local and regional hydrological and agricultural applications. The advent of L-band passive microwave satellites allowed for it to be possible to measure near-surface SM at a global scale compared to in situ measurements. However, their use is often limited because of their coarse spatial resolution. Aiming to address this limitation, random forest (RF) models are adopted to downscale the SMAP level-3 (L3SMP, 36 km) and SMAP enhanced (L3SMP_E, 9 km) SM to 1 km. A suite of predictors derived from the Sentinel-1 C-band SAR and MODIS is used in the downscaling process. The RF models are separately trained and verified at both spatial scales (i.e., 36 and 9 km) considering two experiments: (1) using predictors derived from the MODIS and Sentinel-1 along with other predictors such as elevation and brightness temperature and (2) using all predictors of the first experiment except for the Sentinel-1 predictors. Only dates when the Sentinel-1 images were available are considered for the comparison of the two experiments. The comparison of the results of the two experiments indicates that the removal of Sentinel-1 predictors from the second experiment only reduces the R value from 0.84 to 0.83 and from 0.91 to 0.86 for 36 and 9 km spatial scales, respectively. Among the predictors used in the downscaling, the brightness temperature in VV polarization is identified as the most important predictor, followed by NDVI, surface albedo and API. On the contrary, the Sentinel-1 predictors play a less important role with no marked contribution in enhancing the predictive accuracy of RF models. In general, the two experiments have limitation, such as a small sample size for the training of the RF model because of the scarcity of Sentinel-1 images (i.e., revisit time of 12 days). Therefore, based on this limitation, a third experiment is proposed, in which the Sentinel-1 predictors are not considered at all in the training of the RF models. The results of the third experiment show a good agreement between the downscaled L3SMP and L3SMP_E SM, and in situ SM measurements at both spatial scales. In addition, the temporal availability of the downscaled SM increased. Moreover, the downscaled SM from both SMAP products presented greater spatial detail while preserving the spatial patterns found in their original products. The use of the two SMAP SM products as background fields for the downscaling process does not show marked differences. Overall, this study demonstrates encouraging results in the downscaling of SMAP SM products over humid climate with warm summers dominated by vegetation.

Evaluation of the RF-Based Downscaled SMAP and SMOS Products Using Multi-Source Data over an Alpine Mountains Basin, Northwest China

Passive microwave surface soil moisture (SSM) products tend to have very low resolution, which massively limits their application and validation in regional or local-scale areas. Many climate and hydrological studies are urgently needed to evaluate the suitability of satellite SSM products, especially in alpine mountain areas where soil moisture plays a key role in terrestrial atmospheric exchanges. Aiming to overcome this limitation, a downscaling method based on random forest (RF) was proposed to disaggregate satellite SSM products. We compared the ability of the downscaled soil moisture active passive (SMAP) SSM and soil moisture and ocean salinity satellite (SMOS) SSM products to capture soil moisture information in upstream of the Heihe River Basin by using in situ measurements, the triple collocation (TC) method and temperature vegetation dryness index (TVDI). The results showed that the RF downscaling method has strong applicability in the study area, and the downscaled results of the two products after residual correction have more details, which can better represent the spatial distribution of soil moisture. The validation with the in situ SSM measurements indicates that the correlation between downscaled SMAP and in situ SSM is better than downscaled SMOS at both point and watershed scales in the Babaohe River Basin. From the TC method, the root mean square error (RMSE) of the CLDAS (CMA land data assimilation system), downscaled SMAP and downscaled SMOS were 0.0265, 0.0255 and 0.0317, respectively, indicating that the downscaled SMAP has smaller errors in the study area than others. However, the soil moisture distribution in the study area shown by the SMOS downscaled results is closer than the downscaled SMAP to the degree of drought reflected by TVDI. Overall, this study suggests that the proposed RF-based downscaling method can capture the variation of SSM well, and the downscaled SMAP products perform significantly better than the downscaled SMOS products after the accuracy verification and error analysis of the downscaled results, and it should be helpful to facilitate applications for satellite SSM products at small scales.

Downscaling of MODIS NDVI by Using a Convolutional Neural Network-Based Model with Higher Resolution SAR Data

The normalized difference vegetation index (NDVI) is a simple but powerful indicator, that can be used to observe green live vegetation efficiently. Since its introduction in the 1970s, NDVI has been used widely for land management, food security, and physical models. For these applications, acquiring NDVI in both high spatial resolution and high temporal resolution is preferable. However, there is generally a trade-off between temporal and spatial resolution when using satellite images. To relieve this problem, a convolutional neural network (CNN) based downscaling model was proposed in this research. This model is capable of estimating 10-m high resolution NDVI from MODIS (Moderate Resolution Imaging Spectroradiometer) 250-m resolution NDVI by using Sentinel-1 10-m resolution synthetic aperture radar (SAR) data. First, this downscaling model was trained to estimate Sentinel-2 10-m resolution NDVI from a combination of upscaled 250-m resolution Sentinel-2 NDVI and 10-m resolution Sentinel-1 SAR data, by using data acquired in 2019 in the target area. Then, the generality of this model was validated by applying it to test data acquired in 2020, with the result that the model predicted the NDVI with reasonable accuracy (MAE = 0.090, ρ = 0.734 on average). Next, 250-m NDVI from MODIS data was used as input to confirm this model under conditions replicating an actual application case. Although there were mismatch in the original MODIS and Sentinel-2 NDVI data, the model predicted NDVI with acceptable accuracy (MAE = 0.108, ρ = 0.650 on average). Finally, this model was applied to predict high spatial resolution NDVI using MODIS and Sentinel-1 data acquired in target area from 1 January 2020~31 December 2020. In this experiment, double cropping of cabbage, which was not observable at the original MODIS resolution, was observed by enhanced temporal resolution of high spatial resolution NDVI images (approximately ×2.5). The proposed method enables the production of 10-m resolution NDVI data with acceptable accuracy when cloudless MODIS NDVI and Sentinel-1 SAR data is available, and can enhance the temporal resolution of high resolution 10-m NDVI data.

Evaporative Fluxes and Surface Soil Moisture Retrievals in a Mediterranean Setting from Sentinel-3 and the “Simplified Triangle”

Earth Observation (EO) makes it possible to obtain information on key parameters characterizing interactions among Earth’s system components, such as evaporative fraction (EF) and surface soil moisture (SSM). Notably, techniques utilizing EO data of land surface temperature (Ts) and vegetation index (VI) have shown promise in this regard. The present study investigates, for the first time, the accuracy of one such technique, known as the “simplified triangle”, using Sentinel-3 EO data, acquired for 44 days in 2018 at three savannah FLUXNET sites in Spain. The technique was found to be able to predict both EF and SSM with reasonable accuracy when compared to collocated ground measurements. Comparisons performed for all days together showed relatively low Root Mean square Difference (RMSD) for both EF (0.191) and SSM (0.012 cm3 cm−3) and good correlation coefficients (R) of 0.721 and 0.577, respectively. Both EF and SSM were also largely in agreement with land cover and seasonal variability. The present study comprises the first detailed assessment of the “simplified triangle”, in this case, using Sentinel-3 data and in a Mediterranean setting. Findings, albeit preliminary, are of significant value regarding the use of the investigated technique as a tool of environmental management, and towards ongoing, worldwide efforts aiming at developing operationally relevant products based on the Ts/VI feature space and EO data based on new satellites such as Sentinel-3.

Applications of Remote Sensing in Precision Agriculture: A Review

Agriculture provides for the most basic needs of humankind: food and fiber. The introduction of new farming techniques in the past century (e.g., during the Green Revolution) has helped agriculture keep pace with growing demands for food and other agricultural products. However, further increases in food demand, a growing population, and rising income levels are likely to put additional strain on natural resources. With growing recognition of the negative impacts of agriculture on the environment, new techniques and approaches should be able to meet future food demands while maintaining or reducing the environmental footprint of agriculture. Emerging technologies, such as geospatial technologies, Internet of Things (IoT), Big Data analysis, and artificial intelligence (AI), could be utilized to make informed management decisions aimed to increase crop production. Precision agriculture (PA) entails the application of a suite of such technologies to optimize agricultural inputs to increase agricultural production and reduce input losses. Use of remote sensing technologies for PA has increased rapidly during the past few decades. The unprecedented availability of high resolution (spatial, spectral and temporal) satellite images has promoted the use of remote sensing in many PA applications, including crop monitoring, irrigation management, nutrient application, disease and pest management, and yield prediction. In this paper, we provide an overview of remote sensing systems, techniques, and vegetation indices along with their recent (2015–2020) applications in PA. Remote-sensing-based PA technologies such as variable fertilizer rate application technology in Green Seeker and Crop Circle have already been incorporated in commercial agriculture. Use of unmanned aerial vehicles (UAVs) has increased tremendously during the last decade due to their cost-effectiveness and flexibility in obtaining the high-resolution (cm-scale) images needed for PA applications. At the same time, the availability of a large amount of satellite data has prompted researchers to explore advanced data storage and processing techniques such as cloud computing and machine learning. Given the complexity of image processing and the amount of technical knowledge and expertise needed, it is critical to explore and develop a simple yet reliable workflow for the real-time application of remote sensing in PA. Development of accurate yet easy to use, user-friendly systems is likely to result in broader adoption of remote sensing technologies in commercial and non-commercial PA applications.

A New Approach for Soil Moisture Downscaling in the Presence of Seasonal Difference

The variation of soil moisture (SM) is a complex and synthetic process, which is impacted by numerous factors. The effects of these factors on soil moisture are dynamic. As a result, the relationship between soil moisture and explanatory variables varies with time and season. This kind of change should be considered in obtaining fine spatial resolution soil moisture products. We chose a study area with four distinct seasons in the temperate monsoon region. In this research, we established seasonal downscaling models to avoid the influence of seasonal differences. Precipitation, land surface temperature, evapotranspiration, vegetation index, land cover, elevation, slope, aspect and soil texture were taken as explanatory variables to produce fine spatial resolution SM. SM products derived from Advanced Microwave Scanning Radiometer–Earth Observing System (AMSR-E) and Advanced Microwave Scanning Radiometer 2 (AMSR2) were downscaled with the help of machine learning algorithms. We compared three machine learning algorithms of random forest (RF), support vector machine (SVM), and K-nearest neighbors (KNN) to determine the most suitable algorithm for this study. The results show that season-based downscaling is even better than continuous time series. In the analysis of seasonal differences, precipitation plays a dominant role, but its contribution rate is different in each season. Moreover, the influence of vegetation is more prominent in winter, while the influence of terrain is more important in the other three seasons. It could be noted that the accuracy of the RF model is the best among three machine learning algorithms, and the RF-downscaled products have superior matching performance to both AMSR (AMSR-E and AMSR2) SM products and in-situ measurements. The analysis indicates considering seasonal difference and the application of machine learning has high potential for spatial downscaling in remote sensing applications.

Comparison of spectral and spatial-based approaches for mapping the local variation of soil moisture in a semi-arid mountainous area

Accurate information on soil moisture (SM) is critical in various applications including agriculture, climate, hydrology, soil and drought. In this paper, various predictive relationships including regression (Multiple Linear Regression, MLR), machine learning (Random Forest, RF; Triangular regression, Tr) and spatial modeling (Inverse Distance Weighing, IDW and Ordinary kriging, OK) approaches were compared to estimate SM in a semi-arid mountainous watershed. In developing predictive relationship, Remote Sensing datasets including Landsat 8 satellite imagery derived surface biophysical characteristic, ASTER digital elevation model (DEM) derived surface topographical characteristic, climatic data recorded at the synoptic station and in situ SM data measured at Landsat 8 overpass time were utilized, while in spatial modeling, point-based SM measurements were interpolated. While 70%(calibration set) of the measured SM data were used for modeling, 30%(validation set) were used to evaluate modeling accuracy. Finally, the SM uncertainty maps were created for different models based on a bootstrapping approach. Among the environmental parameter sets, land surface temperature (LST) showed the highest impact on the spatial distribution of SM in the region at all dates. Mean R²(RMSE) between measured and modeled SM on three dates obtained from the MLR, RF, IDW, OK, and Tr models were 0.70(1.97%), 0.72(1.92%), 0.59(2.38%), 0.59(2.27%) and 0.71(1.99%), respectively. The results showed that RF and IDW produced the highest and lowest performance in SM modeling, respectively. Generally, the performance of RS-based models was higher than interpolation models for estimating SM due to the influence from combination of topographic parameters and surface biophysical characteristics. Modeled SM uncertainty with different models varies in the study area. The highest uncertainty in SM modeling was observed at the north part of the study area where the surface heterogeneity is high. Using RS data increased the accuracy of SM modeling because they can capture the surface biophysical characteristics and topographical properties heterogeneity.