Quantifying Long-Term Land Surface and Root Zone Soil Moisture over Tibetan Plateau

Retrieving Soil Physical Properties by Assimilating SMAP Brightness Temperature Observations into the Community Land Model

This paper coupled a unified passive and active microwave observation operator-namely, an enhanced, physically-based, discrete emission-scattering model-with the community land model (CLM) in a data assimilation (DA) system. By implementing the system default local ensemble transform Kalman filter (LETKF) algorithm, the Soil Moisture Active and Passive (SMAP) brightness temperature TBp (p = Horizontal or Vertical polarization) assimilations for only soil property retrieval and both soil properties and soil moisture estimates were investigated with the aid of in situ observations at the Maqu site. The results indicate improved estimates of soil properties of the topmost layer in comparison to measurements, as well as of the profile. Specifically, both assimilations of TBH lead to over a 48% reduction in root mean square errors (RMSEs) for the retrieved clay fraction from the background compared to the top layer measurements. Both assimilations of TBV reduce RMSEs by 36% for the sand fraction and by 28% for the clay fraction. However, the DA estimated soil moisture and land surface fluxes still exhibit discrepancies when compared to the measurements. The retrieved accurate soil properties alone are inadequate to improve those estimates. The discussed uncertainties (e.g., fixed PTF structures) in the CLM model structures should be mitigated.

Global long term daily 1 km surface soil moisture dataset with physics informed machine learning

Although soil moisture is a key factor of hydrologic and climate applications, global continuous high resolution soil moisture datasets are still limited. Here we use physics-informed machine learning to generate a global, long-term, spatially continuous high resolution dataset of surface soil moisture, using International Soil Moisture Network (ISMN), remote sensing and meteorological data, guided with the knowledge of physical processes impacting soil moisture dynamics. Global Surface Soil Moisture (GSSM1 km) provides surface soil moisture (0-5 cm) at 1 km spatial and daily temporal resolution over the period 2000-2020. The performance of the GSSM1 km dataset is evaluated with testing and validation datasets, and via inter-comparisons with existing soil moisture products. The root mean square error of GSSM1 km in testing set is 0.05 cm3/cm3, and correlation coefficient is 0.9. In terms of the feature importance, Antecedent Precipitation Evaporation Index (APEI) is the most important significant predictor among 18 predictors, followed by evaporation and longitude. GSSM1 km product can support the investigation of large-scale climate extremes and long-term trend analysis.

Predicting root zone soil moisture using observations at 2121 sites across China

Root zone soil moisture (RZSM) is particularly useful for understanding hydrological processes, plant-land-atmosphere exchanges, and agriculture- and climate-related research. This study aims to estimate RZSM across China by using a one-parameter (T) exponential filter method (EF method) together with a random forest (RF) regionalization approach and by using a large dataset containing in situ observations collected at 2121 sites across China. First, at each site, T is optimized at each of four soil layers (10-20 cm, 20-30 cm, 30-40 cm and 40-50 cm) by using 0-10-cm soil layer observations and the corresponding calibration layers. Second, an RF classifier is built for each layer according to the calibrated T values and 14 soil, climate and vegetation parameters across 2121 sites. Third, the calibrated T at each soil layer is regionalized with an established RF classifier. Spatial T maps are given for each soil layer across China. Our results show that the EF method performs reasonably well in predicting RZSM at the 10-20-cm, 20-30-cm, 30-40-cm and 40-50-cm layers, with Nash-Sutcliffe efficiency (NSE) medians of 0.73, 0.52, 0.38 and 0.27, respectively, between the observations and estimations. The T parameter shows a spatial pattern in each soil layer and is largely controlled by climate regimes. This study offers an improved RZSM estimation method using a large dataset containing in situ observations; the proposed method also has the potential to be used in other parts of the world.

In Situ Observation-Constrained Global Surface Soil Moisture Using Random Forest Model

The inherent biases of different long-term gridded surface soil moisture (SSM) products, unconstrained by the in situ observations, implies different spatio-temporal patterns. In this study, the Random Forest (RF) model was trained to predict SSM from relevant land surface feature variables (i.e., land surface temperature, vegetation indices, soil texture, and geographical information) and precipitation, based on the in situ soil moisture data of the International Soil Moisture Network (ISMN.). The results of the RF model show an RMSE of 0.05 m3 m−3 and a correlation coefficient of 0.9. The calculated impurity-based feature importance indicates that the Antecedent Precipitation Index affects most of the predicted soil moisture. The geographical coordinates also significantly influence the prediction (i.e., RMSE was reduced to 0.03 m3 m−3 after considering geographical coordinates), followed by land surface temperature, vegetation indices, and soil texture. The spatio-temporal pattern of RF predicted SSM was compared with the European Space Agency Climate Change Initiative (ESA-CCI) soil moisture product, using both time-longitude and latitude diagrams. The results indicate that the RF SSM captures the spatial distribution and the daily, seasonal, and annual variabilities globally.

Monitoring Water and Energy Cycles at Climate Scale in the Third Pole Environment (CLIMATE-TPE)

A better understanding of the water and energy cycles at climate scale in the Third Pole Environment is essential for assessing and understanding the causes of changes in the cryosphere and hydrosphere in relation to changes of plateau atmosphere in the Asian monsoon system and for predicting the possible changes in water resources in South and East Asia. This paper reports the following results: (1) A platform of in situ observation stations is briefly described for quantifying the interactions in hydrosphere-pedosphere-atmosphere-cryosphere-biosphere over the Tibetan Plateau. (2) A multiyear in situ L-Band microwave radiometry of land surface processes is used to develop a new microwave radiative transfer modeling system. This new system improves the modeling of brightness temperature in both horizontal and vertical polarization. (3) A multiyear (2001–2018) monthly terrestrial actual evapotranspiration and its spatial distribution on the Tibetan Plateau is generated using the surface energy balance system (SEBS) forced by a combination of meteorological and satellite data. (4) A comparison of four large scale soil moisture products to in situ measurements is presented. (5) The trajectory of water vapor transport in the canyon area of Southeast Tibet in different seasons is analyzed, and (6) the vertical water vapor exchange between the upper troposphere and the lower stratosphere in different seasons is presented.

Enhanced Estimation of Root Zone Soil Moisture at 1 km Resolution Using SMAR Model and MODIS-Based Downscaled AMSR2 Soil Moisture Data

Root zone soil moisture (RZSM) is an essential variable for weather and hydrological prediction models. Satellite-based microwave observations have been frequently utilized for the estimation of surface soil moisture (SSM) at various spatio-temporal resolutions. Moreover, previous studies have shown that satellite-based SSM products, coupled with the soil moisture analytical relationship (SMAR) can estimate RZSM variations. However, satellite-based SSM products are of low-resolution, rendering the application of the above-mentioned approach for local and pointwise applications problematic. This study initially attempted to estimate SSM at a finer resolution (1 km) using a downscaling technique based on a linear equation between AMSR2 SM data (25 km) with three MODIS parameters (NDVI, LST, and Albedo); then used the downscaled SSM in the SMAR model to monitor the RZSM for Rafsanjan Plain (RP), Iran. The performance of the proposed method was evaluated by measuring the soil moisture profile at ten stations in RP. The results of this study revealed that the downscaled AMSR2 SM data had a higher accuracy in relation to the ground-based SSM data in terms of MAE (↓0.021), RMSE (↓0.02), and R (↑0.199) metrics. Moreover, the SMAR model was run using three different SSM input data with different spatial resolution: (a) ground-based SSM, (b) conventional AMSR2, and (c) downscaled AMSR2 products. The results showed that while the SMAR model itself was capable of estimating RZSM from the variation of ground-based SSM data, its performance increased when using downscaled SSM data suggesting the potential benefits of proposed method in different hydrological applications.

A Novel Scheme for Merging Active and Passive Satellite Soil Moisture Retrievals Based on Maximizing the Signal to Noise Ratio

In this research, we developed and evaluated a new scheme for merging soil moisture (SM) retrievals from both passive and active microwave satellite estimates, based on maximized signal-to-noise ratios, in order to produce improved SM products using least-squares theory. The fractional mean-squared-error (fMSE) derived from the triple collocation method (TCM) was used for this purpose. The proposed scheme was applied by using a threshold between signal and noise at fMSE equal to 0.5 to maintain the high-quality SM observations. In the regions where TCM is unreliable, we propose four scenarios based on the determinations of correlations between all three SM products of TCM at significance levels (i.e., p-values). The proposed scheme was applied to combine SM retrievals from Soil Moisture Active Passive (SMAP), Advanced Scatterometer (ASCAT), and Advanced Microwave Scanning Radiometer 2 (AMSR2) to produce SMAP+ASCAT and AMSR2+ASCAT SM datasets at a global scale for the period from June 2015 to December 2017. The merged SM dataset performance was assessed against SM data from ground measurements of international soil moisture network (ISMN), Global Land Data Assimilation System-Noah (GLDAS-Noah) and ERA5. The results show that the two merged SM datasets showed significant improvement over their parent products in the high average temporal correlation coefficients (R) and the lowest root mean squared difference (RMSE), compared with in-situ measurements over different networks of ISMN. Moreover, these datasets outperformed their parent products over different land cover types in most regions of the world, with a high overall average temporal R and the lowest overall average RMSE value with GLDAS and ERA5. In addition, the suggested scenarios improved SM performance in the regions with unreliable TCMs.

Multiyear in-situ L-band microwave radiometry of land surface processes on the Tibetan Plateau

We report a unique multiyear L-band microwave radiometry dataset collected at the Maqu site on the eastern Tibetan Plateau and demonstrate its utilities in advancing our understandings of microwave observations of land surface processes. The presented dataset contains measurements of L-band brightness temperature by an ELBARA-III microwave radiometer in horizontal and vertical polarization, profile soil moisture and soil temperature, turbulent heat fluxes, and meteorological data from the beginning of 2016 till August 2019, while the experiment is still continuing. Auxiliary vegetation and soil texture information collected in dedicated campaigns are also reported. This dataset can be used to validate the Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite based observations and retrievals, verify radiative transfer model assumptions and validate land surface model and reanalysis outputs, retrieve soil properties, as well as to quantify land-atmosphere exchanges of energy, water and carbon and help to reduce discrepancies and uncertainties in current Earth System Models (ESM) parameterizations. Measurement cases in winter, pre-monsoon, monsoon and post-monsoon periods are presented.

Evaluation of Drought Stress in Cereal through Probabilistic Modelling of Soil Moisture Dynamics

The early and accurate detection of drought episodes is crucial for managing agricultural yield losses and planning adequate policy responses. This study aimed to evaluate the potential of two novel indices, static and dynamic plant water stress, for drought detection and yield prediction. The study was conducted in SW Spain (Córdoba province), covering a 13-year period (2001–2014). The calculation of static and dynamic drought indices was derived from previous ecohydrological work but using a probabilistic simulation of soil moisture content, based on a bucket-type soil water balance, and measured climate data. The results show that both indices satisfactorily detected drought periods occurring in 2005, 2006 and 2012. Both their frequency and length correlated well with annual precipitation, declining exponentially and increasing linearly, respectively. Static and dynamic drought stresses were shown to be highly sensitive to soil depth and annual precipitation, with a complex response, as stress can either increase or decrease as a function of soil depth, depending on the annual precipitation. Finally, the results show that both static and dynamic drought stresses outperform traditional indicators such as the Standardized Precipitation Index (SPI)-3 as predictors of crop yield, and the R2 values are around 0.70, compared to 0.40 for the latter. The results from this study highlight the potential of these new indicators for agricultural drought monitoring and management (e.g., as early warning systems, insurance schemes or water management tools).

Influence of Spatial Resolution on Remote Sensing-Based Irrigation Performance Assessment Using WaPOR Data

This paper analyses the effect of the spatial assessment scale on irrigation performance indicators in small and medium-scale agriculture. Three performance indicators—adequacy (i.e., sufficiency of water use to meet the crop water requirement), equity (i.e., fairness of irrigation distribution), and productivity (i.e., unit of physical crop production/yield per unit water consumption)—are evaluated in five irrigation schemes for three spatial resolutions—250 m, 100 m, and 30 m. Each scheme has varying plot sizes and distributions, with average plot sizes ranging from 0.2 ha to 13 ha. The datasets are derived from the United Nations Food and Agricultural Organization (FAO) water productivity through open access of remotely sensed–derived data (the Water Productivity Open Access Portal—WaPOR) database. Irrigation indicators performed differently in different aspects; for adequacy, all three resolutions show similar spatial trends for relative evapotranspiration (ET) across levels for all years. However, the estimation of relative ET is often higher at higher resolution. In terms of equity, all resolutions show similar inter-annual trends in the coefficient of variation (CV); higher resolutions usually have a higher CV of the annual evapotranspiration and interception (ETIa) while capturing more spatial variability. For productivity, higher resolutions show lower crop water productivity (CWP) due to higher aboveground biomass productivity (AGBP) estimations in lower resolutions; they always have a higher CV of CWP. We find all resolutions of 250 m, 100 m, and 30 m suitable for inter-annual and inter-scheme assessments regardless of plot size. While each resolution shows consistent temporal trends, the magnitude of the trend in both space and time is smoothed by the 100 m and 250 m resolution datasets. This frequently results in substantial differences in the irrigation performance assessment criteria for inter-plot comparisons; therefore, 250 m and 100 m are not recommended for inter-plot comparison for all plot sizes, particularly small plots (<2 ha). Our findings highlight the importance of selecting the spatial resolution appropriate to scheme characteristics when undertaking irrigation performance assessment using remote sensing.

An Integrative Information Aqueduct to Close the Gaps between Satellite Observation of Water Cycle and Local Sustainable Management of Water Resources

The past decades have seen rapid advancements in space-based monitoring of essential water cycle variables, providing products related to precipitation, evapotranspiration, and soil moisture, often at tens of kilometer scales. Whilst these data effectively characterize water cycle variability at regional to global scales, they are less suitable for sustainable management of local water resources, which needs detailed information to represent the spatial heterogeneity of soil and vegetation. The following questions are critical to effectively exploit information from remotely sensed and in situ Earth observations (EOs): How to downscale the global water cycle products to the local scale using multiple sources and scales of EO data? How to explore and apply the downscaled information at the management level for a better understanding of soil-water-vegetation-energy processes? How can such fine-scale information be used to improve the management of soil and water resources? An integrative information flow (i.e., iAqueduct theoretical framework) is developed to close the gaps between satellite water cycle products and local information necessary for sustainable management of water resources. The integrated iAqueduct framework aims to address the abovementioned scientific questions by combining medium-resolution (10 m–1 km) Copernicus satellite data with high-resolution (cm) unmanned aerial system (UAS) data, in situ observations, analytical- and physical-based models, as well as big-data analytics with machine learning algorithms. This paper provides a general overview of the iAqueduct theoretical framework and introduces some preliminary results.