Prediction of plant growth based on statistical methods and remote sensing data

While data like HJ-1 CCD images have advantageous spatial characteristics for describing crop properties, the temporal resolution of the data is rather low, which can be easily made worse by cloud contamination. In contrast, although Moderate Resolution Imaging Spectroradiometer (MODIS) can only achieve a spatial resolution of 250 m in its normalised difference vegetation index (NDVI) product, it has a high temporal resolution, covering the Earth up to multiple times per day. To combine the high spatial resolution and high temporal resolution of different data sources, a new method (Spatial and Temporal Adaptive Vegetation index Fusion Model [STAVFM]) for blending NDVI of different spatial and temporal resolutions to produce high spatial–temporal resolution NDVI datasets was developed based on Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). STAVFM defines a time window according to the temporal variation of crops, takes crop phenophase into consideration and improves the temporal weighting algorithm. The result showed that the new method can combine the temporal information of MODIS NDVI and spatial difference information of HJ-1 CCD NDVI to generate an NDVI dataset with both high spatial and high temporal resolution. An application of the generated NDVI dataset in crop biomass estimation was provided. An average absolute error of 17.2% was achieved. The estimated winter wheat biomass correlated well with observed biomass (R 2 of 0.876). We conclude that the new dataset will improve the application of crop biomass estimation by describing the crop biomass accumulation in detail. There is potential to apply the approach in many other studies, including crop production estimation, crop growth monitoring and agricultural ecosystem carbon cycle research, which will contribute to the implementation of Digital Earth by describing land surface processes in detail.

Surface soil moisture has a wide application in climate change, agronomy, water resources, and in many other domain of science and engineering. Measurement of soil moisture at high spatial and temporal resolution at regional and global scale is needed for the prediction of flood, drought, planning and management of agricultural productivity to ensure food security. Recent advancement in microwave remote sensing, especially after the launch of Sentinel operational satellites has enabled the scientific community to estimate soil moisture at higher spatial and temporal resolution with greater accuracy. This study evaluates the potential of Sentinel-1A satellite images to estimate soil moisture in a semi-arid region. Exactly at the time when satellite passes over the study area, we have collected soil samples at 37 different locations and measured the soil moisture from 5 cm below the ground surface using ML3 theta probe. We processed the soil samples in laboratory to obtain volumetric soil moisture using the oven dry method. We found soil moisture measured from calibrated theta probe and oven dry method are in good agreement with Root Mean Square Error (RMSE) 0.025 m 3 /m 3 and coefficient of determination (R 2 ) 0.85. We then processed Sentinel-1A images and applied modified Dubois model to calculate relative permittivity of the soil from the backscatter values ( σ ∘ ). The volumetric soil moisture at each pixel is then calculated by applying the universal Topp’s model. Finally, we masked the pixels whose Normalised Difference Vegetation Index (NDVI) value is greater than 0.4 to generate soil moisture map as per the Dubois NDVI criterion. Our modelled soil moisture accord with the measured values with RMSE = 0.035 and R 2 = 0.75. We found a small bias in the modelled soil moisture ( 0.02 m 3 / m 3 ). However, this has reduced significantly ( 0.001 m 3 / m 3 ) after applying a bias correction based on Cumulative Distribution Function (CDF) matching. Our approach provides a first-order estimate of soil moisture from Sentinel-1A images in sparsely vegetated agricultural land.

Monitoring impacts of soil bund on spatial variation of teff and finger millet yield with Sentinel-2 and spectroradiometric data in Ethiopia

Satellite data holds considerable potential as a source of information on rice crop growth which can be used to inform agronomy. However, given the typical field sizes in many rice-growing countries such as China, data from coarse spatial resolution satellite systems such as the Moderate Resolution Imaging Spectroradiometer (MODIS) are inadequate for resolving crop growth variability at the field scale. Nevertheless, systems such as MODIS do provide images with sufficient frequency to be able to capture the detail of rice crop growth trajectories throughout a growing season. In order to generate high spatial and temporal resolution data suitable for mapping rice crop phenology, this study fused MODIS data with lower frequency, higher spatial resolution Landsat data. An overall workflow was developed which began with image preprocessing, calculation of multi-temporal normalized difference vegetation index (NDVI) images, and spatiotemporal fusion of data from the two sensors. The Spatial and Temporal Adaptive Reflectance Fusion Model was used to effectively downscale the MODIS data to deliver a time-series of 30 m spatial resolution NDVI data at 8-day intervals throughout the rice-growing season. Zonal statistical analysis was used to extract NDVI time-series for individual fields and signal filtering was applied to the time-series to generate rice phenology curves. The downscaled MODIS NDVI products were able to characterize the development of paddy rice at fine spatial and temporal resolutions, across wide spatial extents over multiple growing seasons. These data permitted the extraction of key crop seasonality parameters that quantified inter-annual growth variability for a whole agricultural region and enabled mapping of the variability in crop performance between and within fields. Hence, this approach can provide rice crop growth data that is suitable for informing agronomic policy and practice across a wide range of scales.

The impacts of global change, including extreme heat and water scarcity, are threatening an ever-growing urban world population. Evapotranspiration (ET) mitigates the urban heat island, reducing the effect of heat waves. It can also be used as a proxy for vegetation water use, making it a crucial tool to plan resilient green cities. To optimize the trade-off between urban greening and water security, reliable and up-to-date maps of ET for cities are urgently needed. Despite its importance, few studies have mapped urban ET accurately for an entire city in high spatial and temporal resolution. We mapped the ET of Berlin, Germany in high spatial (10-m) and temporal (hourly) resolution for the year of 2019. A novel machine learning (ML) approach combining Sentinel-2 time series, open geodata, and flux footprint modeling was applied. Two eddy flux towers with contrasting surrounding land cover provided the training and testing data. Flux footprint modeling allowed us to incorporate comprehensive land cover types in training the ML models. Open remote sensing and geodata used as model inputs included Normalized Difference Vegetation Index (NDVI) from Sentinel-2, building height, impervious surface fraction, vegetation fraction, and vegetation height. NDVI was used to indicate vegetation phenology and health, as plant transpiration contributes to the majority of terrestrial ET. Hourly reference ET (RET) was calculated and used as input to capture the temporal dynamics of the meteorological conditions. Predictions were carried out using Random Forest (RF) regression. Weighted averages extracted from hourly ET maps using flux footprints were compared to measured ET from the two flux towers. Validation showed that the approach is reliable for mapping urban ET, with a mean R2 of 0.76 and 0.56 and a mean RMSE of 0.0289 mm and 0.0171 mm at the more vegetated site and the city-center site, respectively. Lastly, the variation of ET between Local Climate Zones (LCZs) was analyzed to support urban planning. This study demonstrated the capacity to map urban ET at an unprecedented high spatial and temporal resolution with a novel methodology, which can be used to support the sustainable management of green infrastructure and water resources in an urbanizing world facing climate change.

Accurate long-term temperature and precipitation estimates at high spatial and temporal resolutions are vital for a wide variety of climatological studies. We have produced a new, publicly available, daily, gridded maximum temperature, minimum temperature, and precipitation dataset for China with a high spatial resolution of 1 km and over a long-term period (1961 to 2019). It has been named the HRLT and the dataset is publicly available at https://doi.org/10.1594/PANGAEA.941329 (Qin and Zhang, 2022). In this study, the daily gridded data were interpolated using comprehensive statistical analyses, which included machine learning, the generalized additive model, and thin plate splines. It is based on the 0.5° × 0.5° grid dataset from the China Meteorological Administration, together with covariates for elevation, aspect, slope, topographic wetness index, latitude, and longitude. The accuracy of the HRLT daily dataset was assessed using observation data from meteorological stations across China. The maximum and minimum temperature estimates were more accurate than the precipitation estimates. For maximum temperature, the mean absolute error (MAE), root mean square error (RMSE), Pearson’s correlation coefficient (Cor), coefficient of determination after adjustment (R2), and Nash-Sutcliffe modeling efficiency (NSE) were 1.07 °C, 1.62 °C, 0.99, 0.98, and 0.98, respectively. For minimum temperature, the MAE, RMSE, Cor, R2, and NSE were 1.08 °C, 1.53 °C, 0.99, 0.99, and 0.99, respectively. For precipitation, the MAE, RMSE, Cor, R2, and NSE were 1.30 mm, 4.78 mm, 0.84, 0.71, and 0.70, respectively. The accuracy of the HRLT was compared to those of the other three existing datasets and its accuracy was either greater than the others, especially for precipitation, or comparable in accuracy, but with higher spatial resolution or over a longer time period. In summary, the HRLT dataset, which has a high spatial resolution, covers a longer period of time and has reliable accuracy, is suitable for future environmental analyses, especially the effects of extreme weather.

Monitoring crop growth and estimating crop yield are essential for managing agricultural production, ensuring food security, and maintaining sustainable agricultural development. Combining the mechanistic framework of a crop growth model with remote sensing observations can provide a means of generating realistic and spatially detailed crop growth information that can facilitate accurate crop yield estimates at different scales. The main objective of this study was to develop a robust estimation methodology of within-field winter wheat yield at a high spatial resolution (20 m × 20 m) by combining a light use efficiency-based model and Sentinel-2 data. For this purpose, Sentinel-2 derived leaf area index (LAI) time series were assimilated into the Simple Algorithm for Yield Estimation (SAFY) model using an ensemble Kalman filter (EnKF). The study was conducted on rainfed winter wheat fields in southern Sweden. LAI was estimated using vegetation indices (VIs) derived from Sentinel-2 data with semi-empirical models. The enhanced two-band vegetation index (EVI2) was found to be a useful VI for LAI estimation, with a coefficient of determination (R2) and a root mean square error (RMSE) of 0.80 and 0.65 m2/m2, respectively. Our findings demonstrate that the assimilation of LAI derived from Sentinel-2 into the SAFY model using EnKF enhances the estimation of within-field spatial variability of winter wheat yield by 70% compared to the baseline simulation without the assimilation of remotely sensed data. Additionally, the assimilation of LAI improves the accuracy of winter wheat yield estimation by decreasing the RMSE by 53%. This study demonstrates an approach towards practical applications of freely accessible Sentinel-2 data and a crop growth model through data assimilation for fine-scale mapping of crop yield. Such information is critical for quantifying the yield gap at the field scale, and to aid the optimization of management practices to increase crop production.

Parameterization of air temperature in high temporal and spatial resolution from a combination of the SEVIRI and MODIS instruments

Application of Landsat-8 and Sentinel-2 for retrieval of chlorophyll-a in a shallow freshwater lake

To identify factors affecting and key environmental factors of the Oncomelania hupensis snail density in the Yangtze River Delta region using machine learning methods. Administrative village-level O. hupensis snail survey data in the Yangtze River Delta (including Shanghai Municipality, Jiangsu Province, Zhejiang Province and Anhui Province) from 2011 to 2021 were retrieved from the Information Management System for Parasitic Disease Control of Chinese Center for Disease Control and Prevention. Environmental factor data were captured from the Google Earth Engine platform, including elevation, slope, terrain, normalized difference vegetation index (NDVI), vegetation type, soil type, total petroleum hydrocarbon (TPH), ammonium nitrogen, inorganic nitrogen, dissolved oxygen, pH of water, chemical oxygen demand (COD) and inorganic phosphorus, and climatic factor data in the study region were retrieved from the Copernicus Climate Data Store, including annual precipitation, aridity index and annual mean temperature (AMT). O. hupensis snail survey data in the Yangtze River Delta region from 2011 to 2021 were randomly divided into a training set (70%) and a test set (30%), and five machine learning models were selected for machine learning model construction and comparative analysis of the O. hupensis snail density using the software R 4.3.0, including random forest (RF), eXtreme gradient boosting (XGBoost), support vector machine (SVM), gradient boosting machine (GBM) and neural network (NN). The XGBoost model was employed to construct a predictive model for the O. hupensis snail density, and the impact of each environmental factor on O. hupensis snail distribution was quantified. The SHapley Additive exPlanations (SHAPs) values were calculated to estimate the average contribution of each variable to the model prediction, and the core environmental factors affecting the O. hupensis snail population density were screened. Among the five machine learning models, the XGBoost model exhibited the optimal comprehensive performance, with the coefficient of determination (R2) of 0.855, mean squared error (MSE) of 0.188, root mean squared error (RMSE) of 0.434 and mean absolute error (MAE) of 0.155, respectively. Analysis of factors affecting the O. hupensis snail density with the XGBoost model showed that among the 16 environmental factors, the top four high-impact factors ranked by SHAPs values included annual precipitation, elevation, aridity index and NDVI, with cumulative SHAPs contributions of 75%, which was higher than that of other environmental factors. If NDVI was higher than 0.6, the O. hupensis snail density increased with NDVI and peaked if NDVI was 0.8 (1.60 snails/0.1 m2). The O. hupensis snail density increased with elevation if the elevation ranged from 14 to 40 m, and slowly rose if the annual precipitation ranged from 900 to 1 300 mm, and then increased rapidly to the peak (1.52 snails/0.1 m2) if the annual precipitation ranged from 1 300 to 1 500 mm. In addition, the O. hupensis snail density increased rapidly to the maximum (1.60 snails/0.1 m2) if the aridity index ranged from 0.8 to 1.1, and decreased gradually if the aridity index exceeded 1.1. The XGBoost model shows excellent performance in prediction of the O. hupensis snail density and identification of key environmental factors in the Yangtze River Delta region. Annual precipitation, elevation, aridity index and NDVI are key environmental factors affecting the distribution and density of O. hupensis snails in the Yangtze River Delta region.

Estimation of crop gross primary production (GPP): I. impact of MODIS observation footprint and impact of vegetation BRDF characteristics

Coral-spawn slicks: Reflectance spectra and detection using optical satellite data

Normalized difference vegetation index (NDVI) time series data, derived from optical images, play a crucial role for crop mapping and growth monitoring. Nevertheless, optical images frequently exhibit spatial and temporal discontinuities due to cloudy and rainy weather conditions. Existing algorithms for reconstructing NDVI time series using multi-source remote sensing data still face several challenges. In this study, we proposed a novel method, an enhanced gap-filling and Whittaker smoothing (EGF-WS), to reconstruct NDVI time series (EGF-NDVI) using Google Earth Engine. In EGF-WS, NDVI calculated from MODIS, Landsat-8, and Sentinel-2 satellites were combined to generate high-resolution and continuous NDVI time series data. The MODIS NDVI was employed as reference data to fill missing pixels in the Sentinel–Landsat NDVI (SL-NDVI) using the gap-filling method. Subsequently, the filled NDVI was smoothed using a Whittaker smoothing filter to reduce residual noise in the SL-NDVI time series. With reference to the all-round performance assessment (APA) metrics, the performance of EGF-WS was compared with the conventional gap-filling and Savitzky–Golay filter approach (GF-SG) in Fusui County of Guangxi Zhuang Autonomous Region. The experimental results have demonstrated that the EGF-WS can capture more accurate spatial details compared with GF-SG. Moreover, EGF-NDVI of Fusui County exhibited a low root mean square error (RMSE) and a high coefficient of determination (R2). In conclusion, EGF-WS holds significant promise in providing NDVI time series images with a spatial resolution of 10 m and a temporal resolution of 8 days, thereby benefiting crop mapping, land use change monitoring, and various ecosystems, among other applications.

Several studies have demonstrated the potential of actively heated fiber optics (AHFO) to measure soil water content (SWC) at high spatial and temporal resolutions. This study tested the feasibility of the AHFO technique to measure soil water in the surface soil of a crop grown field over a growing season using an in-situ calibration approach. Heat pulses of five minutes duration were applied at a rate of 7.28 W m−1 along eighteen fiber optic cable transects installed at three depths (0.05, 0.10 and 0.20 m) at six-hour intervals. Cumulative temperature increase (Tcum) during heat pulses was calculated at locations along the cable. While predicting commercial sensor measurements, the AHFO showed root mean square errors (RMSE) of 2.8, 3.7 and 3.7% for 0.05, 0.10 and 0.20 m depths, respectively. Further, the coefficients of determination (R2) for depth specific relationships were 0.87 (0.05 m depth), 0.46 (0.10 m depth), 0.86 (0.20 m depth) and 0.66 (all depths combined). This study showed a great potential of the AHFO technique to measure soil water at high spatial resolutions (<1 m) and to monitor soil water dynamics of surface soil in a crop grown field over a cropping season with a reasonable compromise between accuracy and practicality.

The precise mapping of coastal wetlands holds great significance in the context of monitoring carbon sequestration and storage within coastal ecosystems, particularly in light of climate change and human-induced activities. Time series and multi-source remote sensing data offers distinct advantage in the spatial and temporal mapping of land use, particularly in wetland systems, encompassing various types of vegetation. The primary aim of this study was to delineate the spatial distribution of land use within the Liao River Delta (LRD) wetland. This was achieved by employing a stacking ensemble model that integrated time series GaoFen-1 (GF-1) optical imagery, GaoFen-3 (GF-3) synthetic aperture radar (SAR) imagery, and Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) imagery. The first step involved the application of the enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM) to fuse the GF-1 NDVI and MODIS NDVI datasets, resulting in the production of time series NDVI data. Subsequently, the parameters pertaining to vegetation phenology were obtained by employing the threshold method on time series NDVI data. We compiled feature datasets that encompassed GF-1 spectral bands, spectral indices, phenological parameters, and GF-3 SAR features. In order to mitigate data redundancy, the Recursive Feature Elimination and Cross-Validation (RFECV) model was employed to identify and select significant features. Finally, the stacking ensemble model was constructed by combining five base models [K-Nearest Neighbors (KNN), Random Forest (RF), Adaptive Boosting (AdaBoost), Extreme Gradient Boosting (XGBoost), and Light Gradient Boosting Machine (LightGBM)] to perform wetland classification. The findings showed the following: (1) ESTARFM was able to successfully fuse GF-1 NDVI and MODIS NDVI data, resulting in a spatiotemporal fusion with a coefficient of determination (R2) of 0.85 and a root mean square error (RMSE) of 0.07. (2) In the process of wetland classification, the RFECV algorithm was employed to select relevant features. Specifically, 75 spectral band features, 89 spectral index features, 13 SAR features, and 7 phenological parameters were identified as significant for this task. (3) A stacking ensemble model was constructed using the aforementioned multi-source features. This model exhibited a robust and consistent performance in wetland classification, achieving the highest overall accuracy of 94.33%. Notably, this accuracy improvement ranged from approximately 0.09% to 10.02% when compared to the individual base models. Thus, the present study has the potential to be utilized in the context of fine-scale wetland monitoring, thereby offering valuable assistance to the field of wetland environmental research.