This paper describes the technical development and accuracy assessment of the most recent and improved version of the SoilGrids system at 250m resolution (June 2016 update). SoilGrids provides global predictions for standard numeric soil properties (organic carbon, bulk density, Cation Exchange Capacity (CEC), pH, soil texture fractions and coarse fragments) at seven standard depths (0, 5, 15, 30, 60, 100 and 200 cm), in addition to predictions of depth to bedrock and distribution of soil classes based on the World Reference Base (WRB) and USDA classification systems (ca. 280 raster layers in total). Predictions were based on ca. 150,000 soil profiles used for training and a stack of 158 remote sensing-based soil covariates (primarily derived from MODIS land products, SRTM DEM derivatives, climatic images and global landform and lithology maps), which were used to fit an ensemble of machine learning methods—random forest and gradient boosting and/or multinomial logistic regression—as implemented in the packages , , and . The results of 10–fold cross-validation show that the ensemble models explain between 56% (coarse fragments) and 83% (pH) of variation with an overall average of 61%. Improvements in the relative accuracy considering the amount of variation explained, in comparison to the previous version of SoilGrids at 1 km spatial resolution, range from 60 to 230%. Improvements can be attributed to: (1) the use of machine learning instead of linear regression, (2) to considerable investments in preparing finer resolution covariate layers and (3) to insertion of additional soil profiles. Further development of SoilGrids could include refinement of methods to incorporate input uncertainties and derivation of posterior probability distributions (per pixel), and further automation of spatial modeling so that soil maps can be generated for potentially hundreds of soil variables. Another area of future research is the development of methods for multiscale merging of SoilGrids predictions with local and/or national gridded soil products (e.g. up to 50 m spatial resolution) so that increasingly more accurate, complete and consistent global soil information can be produced. SoilGrids are available under the Open Data Base License.

Soil moisture is a key environmental variable, important to e.g. farmers, meteorologists, and disaster management units. Here, we present a method to retrieve Surface Soil Moisture (SSM) from the Sentinel-1 satellites, which carry C-band Synthetic Aperture Radar (S-1 CSAR) sensors that provide the richest freely available SAR data source so far, unprecedented in accuracy and coverage. Our SSM retrieval method, adapting well-established change detection algorithms, builds the first globally deployable soil moisture observation dataset with 1km resolution. The paper provides an algorithm formulation to be operated in data cube architectures and High Performance Computing (HPC) environments. It includes the novel Dynamic Gaussian Upscaling (DGU) method for spatial upscaling of SAR imagery, harnessing its field-scale information and successfully mitigating effects from the SAR's high signal complexity. Also, a new regression-based approach for estimating the radar slope is defined, coping with Sentinel-1's inhomogeneity in spatial coverage. We employ the S-1 SSM algorithm on a 3yr S-1 data cube over Italy, obtaining a consistent set of model parameters and product masks, unperturbed by coverage discontinuities. An evaluation of therefrom generated S-1 SSM data, involving a 1km Soil Water Balance Model (SWBM) over Umbria, yields high agreement over plains and agricultural areas, with low agreement over forests and strong topography. While positive biases during the growing season are detected, excellent capability to capture small-scale soil moisture changes as from rainfall or irrigation is evident. The S-1 SSM is currently in preparation towards operational product dissemination in the Copernicus Global Land Service (CGLS).

Crop monitoring is of great importance for e.g., yield prediction and increasing water use efficiency. The Copernicus Sentinel-1 mission operated by the European Space Agency provides the opportunity to monitor Earth’s surface using radar at high spatial and temporal resolution. Sentinel-1’s Synthetic Aperture Radar provides co- and cross-polarized backscatter, enabling the calculation of microwave indices. In this study, we assess the potential of Sentinel-1 VV and VH backscatter and their ratio VH/VV, the cross ratio (CR), to monitor crop conditions. A quantitative assessment is provided based on in situ reference data of vegetation variables for different crops under varying meteorological conditions. Vegetation Water Content (VWC), biomass, Leaf Area Index (LAI) and height are measured in situ for oilseed-rape, corn and winter cereals at different fields during two growing seasons. To quantify the sensitivity of backscatter and microwave indices to vegetation dynamics, linear and exponential models and machine learning methods have been applied to the Sentinel-1 data and in situ measurements. Using an exponential model, the CR can account for 87% and 63% of the variability in VWC for corn and winter cereals. In oilseed-rape, the coefficient of determination ( R 2 ) is lower ( R 2 = 0.34) due to the large difference in VWC between the two growing seasons and changes in vegetation structure that affect backscatter. Findings from the Random Forest analysis, which uses backscatter, microwave indices and soil moisture as input variables, show that CR is by and large the most important variable to estimate VWC. This study demonstrates, based on a quantitative analysis, the large potential of microwave indices for vegetation monitoring of VWC and phenology.

No abstract

Soil moisture is a key environmental variable, important to e.g., farmers, meteorologists, and disaster management units. We fuse surface soil moisture (SSM) estimates from spatio-temporally complementary radar sensors through temporal filtering of their joint signal and obtain a kilometre-scale, daily soil water content product named SCATSAR-SWI. With 25 km Metop ASCAT SSM and 1 km Sentinel-1 SSM serving as input, the SCATSAR-SWI is globally applicable and achieves daily full coverage over operated areas. We employ a near-real-time-capable SCATSAR-SWI algorithm on a fused 3 year ASCAT-Sentinel-1-SSM data cube over Italy, obtaining a consistent set of model parameters, unperturbed by coverage discontinuities. An evaluation of a therefrom generated SCATSAR-SWI dataset, involving a 1 km Soil Water Balance Model (SWBM) over Umbria, yields comprehensively high agreement with the reference data (median R = 0.61 vs. in situ; 0.71 vs. model; 0.83 vs. ASCAT SSM). While the Sentinel-1 signal is attenuated to some extent, the ASCAT's signal dynamics are fully transferred to the SCATSAR-SWI and benefit from the Sentinel-1 parametrisation. Using the SM2RAIN approach, the SCATSAR-SWI shows excellent capability to reproduce 5 day-accumulated rainfall over Italy, with R = 0.89 against observed rainfall. The SCATSAR-SWI is currently in preparation towards operational product dissemination in the Copernicus Global Land Service (CGLS).

Recently, the slope and curvature estimation of the backscatter-incidence angle relationship within the TU Wien retrieval algorithm has been improved. Where previously only climatologies of the slope and curvature parameters were available, i.e., one value for every day of year, slope and curvature are now calculated for every day. This enables the retrieval of time series of vegetation optical depth (τ a ) from backscatter observations. This study demonstrates the ability to detect interannual variability in vegetation dynamics using τ a derived from backscatter provided by the advanced scatterometer on-board Metop-A. τ a time series over Australia for the period 2007-2014 are compared to leaf area index (LAI) from SPOT-VEGETATION by calculating the rank correlation coefficient (r s ) for original time series and anomalies. High values for r s are found over bare soil and sparse vegetation in central Australia with median r s values of 0.78 and 0.58, respectively. Forests and ephemeral lakes and rivers impact the retrieval of τ a , and the negative values for r s are found in these areas. Looking at the annual averages of τ a , LAI, and surface soil moisture, significantly high values are found for the anomalously wet years 2010 and 2011. Patterns in the increased τ a correspond to regions with increased soil moisture and LAI. Values for τ a and LAI are anomalous especially in sparsely vegetated regions, where the flush of grasses increases τ a and LAI. Regions with enough precipitation and higher woody vegetation component show a smaller increase in 2010 and 2011. This study demonstrates the skill of τ a , and subsequently of scatterometers, to monitor the vegetation dynamics thanks to the multiincidence angle observation capability.

Vegetation products based on microwave remote sensing observations, such as Vegetation Optical Depth (VOD), are increasingly used in a variety of applications. One disadvantage is the often coarse spatial resolution of tens of kilometers of products retrieved from microwave observations from spaceborne radiometers and scatterometers. This can potentially be overcome byusing new high-resolution Synthetic Aperture Radar (SAR) observations from Sentinel-1. However, the sensitivity of Sentinel-1 backscatter to vegetation dynamics, or its use in radiative transfer models, such as the water cloud model, has only been tested at field to regional scale. In this study, we compared the cross-polarization ratio (CR) to vegetation dynamics as observed in microwave-based Vegetation Optical Depth from coarse-scale satellites over Europe. CR was obtained from Sentinel-1 VH and VV backscatter observations at 500 m sampling and resampled to the spatial resolution of VOD from the Advanced SCATterometer (ASCAT) on-board the Metop satellite series. Spatial patterns between median CR and ASCAT VOD correspond to each other and to vegetation patterns over Europe. Analysis of temporal correlation between CR and ASCAT VOD shows that high Pearson correlation coefficients (Rp) are found over croplands and grasslands (median Rp > 0.75). Over deciduous broadleaf forests, negative correlations are found. This is attributed to the effect of structural changes in the vegetation canopy which affect CR and ASCAT VOD in different ways. Additional analysis comparing CR to passive microwave-based VOD shows similar effects in deciduous broadleaf forests and high correlations over crop- and grasslands. Though the relationship between CR and VOD over deciduous forests is unclear, results suggest that CR is useful for monitoring vegetation dynamics over crop- and grassland and a potential path to high-resolution VOD.

Australia is frequently subject to droughts and floods. Its hydrology is strongly connected to oceanic and atmospheric oscillations (climate modes) such as the El Niño-Southern Oscillation (ENSO), Indian Ocean dipole (IOD), and southern annular mode (SAM). A global 32-yr dataset of remotely sensed surface soil moisture (SSM) was used to examine hydrological variations in mainland Australia for the period 1978-2010. Complex empirical orthogonal function (CEOF) analysis was applied to extract independent signals and to investigate their relationships to climate modes. The annual cycle signal represented 46.3% of the total variance and a low but highly significant connection with SAM was found. Two multiannual signals with a lesser share in total variance (6.3% and 4.2%) were identified. The first one had an unstable period of 2-5 yr and reflected an east-west pattern that can be associated with ENSO and SAM but not with IOD. The second one, a 1-to 5-yr oscillation, formed a dipole pattern between the west and north and can be linked to ENSO and IOD. As expected, relationships with ENSO were found throughout the year and are especially strong during southern spring and summer in the east and north. Somewhat unexpectedly, SAM impacts strongest in the north and east during summer and is proposed as the key driver of the annual SSM signal. The IOD explains SSM variations in the north, east, and southeast during spring and also in the west during winter.