Gaussian processes retrieval of crop traits in Google Earth Engine based on Sentinel-2 top-of-atmosphere data

Extraction and spatiotemporal changes of open-pit mines during 1985-2020 using Google Earth Engine: A case study of Qingzhou City, Shandong Province, China

The global use of mineral resources has increased exponentially for decades and will continue to grow for the foreseeable future, resulting in increasingly negative impacts on the surrounding environment. However, to date, there are a lack of historical and current spatial extent datasets with high accuracy for mining areas in many parts of the world, which has hindered a more comprehensive understanding of the environmental impacts of mining. Using the Google Earth Engine cloud platform and the Landsat normalized difference vegetation index (NDVI) datasets, the spatial extent data of open-pit mining areas for eight years (1985, 1990, 1995, 2000, 2005, 2010, 2015, and 2020) was extracted by the Otsu algorithm. The limestone mining areas in Qingzhou, Shandong Province, China, was selected as a case study. The annual maximum NDVI was first derived from the Landsat NDVI datasets, and then the Otsu algorithm was used to segment the annual maximum NDVI images to obtain the extent of the mining areas. Finally, the spatiotemporal characteristics of the mining areas in the study region were analyzed in reference to previous survey data. The results showed that the mining areas were primarily located in Shaozhuang Town, Wangfu Street and the northern part of Miaozi Town, and the proportion of mining areas within these three administrative areas has increased annually from 88% in 1985 to more than 98% in 2010. Moreover, the open-pit mining areas in ​​Qingzhou gradually expanded from a scattered, point-like distribution to a large, contiguous distribution. From 1985 to 2020, the open-pit mining area expanded to more than 10 times its original size at a rate of 0.5 km²/year. In 2015, this area reached its maximum size of 19.7 km² and slightly decreased in 2020. Furthermore, the expansion of the mining areas in Qingzhou went through three stages: a slow growth period before 1995, a rapid expansion period from 1995 to 2005, and a shutdown and remediation period after 2005. A quantitative accuracy assessment was performed by calculating the Intersection over Union (IoU) of the extraction results and the visual interpretation results from Gaofen-2 images with 1-m spatial resolution. The IoU reached 72%. The results showed that it was feasible to threshold the Landsat annual maximum NDVI data by the Otsu algorithm to extract the annual spatial extent of the open-pit mining areas. Our method will be easily transferable to other regions worldwide, enabling the monitoring of mine environments.

Quantifying Irrigated Winter Wheat LAI in Argentina Using Multiple Sentinel-1 Incidence Angles

Synthetic aperture radar (SAR) data provides an appealing opportunity for all-weather day or night Earth surface monitoring. The European constellation Sentinel-1 (S1) consisting of S1-A and S1-B satellites offers a suitable revisit time and spatial resolution for the observation of croplands from space. The C-band radar backscatter is sensitive to vegetation structure changes and phenology as well as soil moisture and roughness. It also varies depending on the local incidence angle (LIA) of the SAR acquisition's geometry. The LIA backscatter dependency could therefore be exploited to improve the retrieval of the crop biophysical variables. The availability of S1 radar time-series data at distinct observation angles holds the feasibility to retrieve leaf area index (LAI) evolution considering spatiotemporal coverage of intensively cultivated areas. Accordingly, this research presents a workflow merging multi-date S1 smoothed data acquired at distinct LIA with a Gaussian processes regression (GPR) and a cross-validation (CV) strategy to estimate cropland LAI of irrigated winter wheat. The GPR-S1-LAI model was tested against in situ data of the 2020 winter wheat campaign in the irrigated valley of Colorador river, South of Buenos Aires Province, Argentina. We achieved adequate validation results for LAI with R C V 2 = 0.67 and RMSE _CV = 0.88 m² m^-2. The trained model was further applied to a series of S1 stacked images, generating temporal LAI maps that well reflect the crop growth cycle. The robustness of the retrieval workflow is supported by the associated uncertainties along with the obtained maps. We found that processing S1 smoothed imagery with distinct acquisition geometries permits accurate radar-based LAI modeling throughout large irrigated areas and in consequence can support agricultural management practices in cloudprone agri-environments.

Retrieval of carbon content and biomass from hyperspectral imagery over cultivated areas

Spaceborne imaging spectroscopy is a highly promising data source for all agricultural management and research disciplines that require spatio-temporal information on crop properties. Recently launched science-driven missions, such as the Environmental Mapping and Analysis Program (EnMAP), deliver unprecedented data from the Earth's surface. This new kind of data should be explored to develop robust retrieval schemes for deriving crucial variables from future routine missions. Therefore, we present a workflow for inferring crop carbon content (Carea ), and aboveground dry and wet biomass (AGBdry , AGBfresh ) from EnMAP data. To achieve this, a hybrid workflow was generated, combining radiative transfer modeling (RTM) with machine learning regression algorithms. The key concept involves: (1) coupling the RTMs PROSPECT-PRO and 4SAIL for simulation of a wide range of vegetation states, (2) using dimensionality reduction to deal with collinearity, (3) applying a semi-supervised active learning technique against a 4-years campaign dataset, followed by (4) training of a Gaussian process regression (GPR) machine learning algorithm and (5) validation with an independent in situ dataset acquired during the ESA Hypersense experiment campaign at a German test site. Internal validation of the GPR-Carea and GPR-AGB models achieved coefficients of determination (R 2) of 0.80 for Carea and 0.80, 0.71 for AGBdry and AGBfresh , respectively. The mapping capability of these models was successfully demonstrated using airborne AVIRIS-NG hyperspectral imagery, which was spectrally resampled to EnMAP spectral properties. Plausible estimates were achieved over both bare and green fields after adding bare soil spectra to the training data. Validation over green winter wheat fields generated reliable estimates as suggested by low associated model uncertainties (< 40%). These results suggest a high degree of model reliability for cultivated areas during active growth phases at canopy closure. Overall, our proposed carbon and biomass models based on EnMAP spectral sampling demonstrate a promising path toward the inference of these crucial variables over cultivated areas from future spaceborne operational hyperspectral missions.

Seasonal Mapping of Irrigated Winter Wheat Traits in Argentina with a Hybrid Retrieval Workflow Using Sentinel-2 Imagery

Earth observation offers an unprecedented opportunity to monitor intensively cultivated areas providing key support to assess fertilizer needs and crop water uptake. Routinely, vegetation traits mapping can help farmers to monitor plant development along the crop's phenological cycle, which is particularly relevant for irrigated agricultural areas. The high spatial and temporal resolution of the Sentinel-2 (S2) multispectral instrument leverages the possibility to estimate leaf area index (LAI), canopy chlorophyll content (CCC), and vegetation water content (VWC) from space. Therefore, our study presents a hybrid retrieval workflow combining a physically-based strategy with a machine learning regression algorithm, i.e., Gaussian processes regression, and an active learning technique to estimate LAI, CCC and VWC of irrigated winter wheat. The established hybrid models of the three traits were validated against in-situ data of a wheat campaign in the Bonaerense valley, South of the Buenos Aires Province, Argentina, in the year 2020. We obtained good to highly accurate validation results with LAI: R² = 0.92, RMSE = 0.43 m² m^-2, CCC: R² = 0.80, RMSE = 0.27 g m^-2 and VWC: R² = 0.75, RMSE = 416 g m^-2. The retrieval models were also applied to a series of S2 images, producing time series along the seasonal cycle, which reflected the effects of fertilizer and irrigation on crop growth. The associated uncertainties along with the obtained maps underlined the robustness of the hybrid retrieval workflow. We conclude that processing S2 imagery with optimised hybrid models allows accurate space-based crop traits mapping over large irrigated areas and thus can support agricultural management decisions.