Search

Large-scale crop phenology monitoring is essential for agro-ecosystem policy. Remote sensing helps track crop development but requires high-temporal and spatial resolutions. While datasets with both attributes are now available, their large-scale applications require significant resources. Medium-resolution data offer daily observations but lack detail for smaller plots. This study generated crop-specific phenomaps for mainland France (2016–2020) using PROBA-V data. A spatial disaggregation method reconstructed NDVI time series for individual crops within mixed pixels. Then, phenometrics were extracted from disaggregated PROBA-V and Sentinel-2 separately and compared to observed phenological stages. Results showed that PROBA-V-based phenomaps closely matched observations at regional level, with moderate accuracy at municipal level. PROBA-V demonstrated a higher detection rate than Sentinel-2, especially in cloudy periods, and successfully generated phenomaps before Sentinel-2B's launch. The study highlights PROBA-V's potential for operational crop monitoring, i.e., wheat heading and oilseed rape flowering, with performance comparable to Sentinel-2. PROBA-V outputs complement Sentinel-2: phenometrics cannot be generated at plot level but are efficiently produced at regional or national scales to study phenological gradients more easily than with Sentinel-2 and with similar accuracy. This approach could be extended to MODIS or SPOT-VGT, to generate historical phenological data, providing that a crop map is available. [ABSTRACT FROM AUTHOR]

Accurate identification of crop phenology timing is crucial for agriculture. While remote sensing tracks vegetation changes, linking these to ground-measured crop growth stages remains challenging. Existing methods offer broad overviews but fail to capture detailed phenological changes, which can be partially related to the temporal resolution of the remote sensing datasets used. The availability of higher-frequency observations, obtained by combining sensors and gap-filling, offers the possibility to capture more subtle changes in crop development, some of which can be relevant for management decisions. One such dataset is Planet Fusion, daily analysis-ready data obtained by integrating PlanetScope imagery with public satellite sensor sources such as Sentinel-2 and Landsat. This study introduces a novel method utilizing Dynamic Time Warping applied to Planet Fusion imagery for maize phenology detection, to evaluate its effectiveness across 70 micro-stages. Unlike singular template approaches, this method preserves critical data patterns, enhancing prediction accuracy and mitigating labeling issues. During the experiments, eight commonly employed spectral indices were investigated as inputs. The method achieves high prediction accuracy, with 90% of predictions falling within a 10-day error margin, evaluated based on over 3200 observations from 208 fields. To understand the potential advantage of Planet Fusion, a comparative analysis was performed using Harmonized Landsat Sentinel-2 data. Planet Fusion outperforms Harmonized Landsat Sentinel-2, with significant improvements observed in key phenological stages such as V4, R1, and late R5. Finally, this study showcases the method's transferability across continents and years, although additional field data are required for further validation. [ABSTRACT FROM AUTHOR]

Crop identification at the field level using remote sensing data is a very important task. However, the use of multispectral data for the construction of vegetation indices is sometimes impossible or limited. For such situations, solutions based on the use of time series of synthetic aperture radar (SAR) indices are promising, eliminating the problems associated with cloudiness and providing an assessment of crop development characteristics during the growing season. We evaluated the use of time series of synthetic aperture radar (SAR) indices to characterize crop development during the growing season. The use of SAR imagery for crop identification addresses issues related to cloudiness. Therefore, it is important to choose the SAR index that is the most stable and has the lowest spatial variability throughout the growing season while being comparable to the normalized difference vegetation index (NDVI). The presented work is devoted to the study of these issues. In this study, the spatial variabilities of different SAR indices time series were compared for a single region for the first time to identify the most stable index for use in precision agriculture, including the in-field heterogeneity of crop sites, crop rotation control, mapping, and other tasks in various agricultural areas. Seventeen Sentinel-1B images of the southern part of the Khabarovsk Territory in the Russian Far East at a spatial resolution of 20 m and temporal resolution of 12 days for the period between 14 April 2021 and 1 November 2021 were obtained and processed to generate vertical–horizontal/vertical–vertical polarization (VH/VV), radar vegetation index (RVI), and dual polarimetric radar vegetation index (DpRVI) time series. NDVI time series were constructed from multispectral Sentinel-2 images using a cloud cover mask. The characteristics of time series maximums were calculated for different types of crops: soybean, oat, buckwheat, and timothy grass. The DpRVI index exhibited the highest stability, with coefficients of variation of the time series that were significantly lower than those for RVI and VH/VV. The main characteristics of the SAR and NDVI time series—the maximum values, the dates of the maximum values, and the variability of these indices—were compared. The variabilities of the maximum values and dates of maximum values for DpRVI were lower than for RVI and VH/VV, whereas the variabilities of the maximum values and the dates of maximum values were comparable for DpRVI and NDVI. On the basis of the DpRVI index, classifications were carried out using seven machine learning methods (fine tree, quadratic discriminant, Gaussian naïve Bayes, fine k nearest neighbors or KNN, random under-sampling boosting or RUSBoost, random forest, and support vector machine) for experimental sites covering a total area of 1009.8 ha. The quadratic discriminant method yielded the best results, with a pixel classification accuracy of approximately 82% and a kappa value of 0.67. Overall, 90% of soybean, 74.1% of oat, 68.9% of buckwheat, and 57.6% of timothy grass pixels were correctly classified. At the field level, 94% of the fields included in the test dataset were correctly classified. The paper results show that the DpRVI can be used in cases where the NDVI is limited, allowing for the monitoring of phenological development and crop mapping. The research results can be used in the south of Khabarovsk Territory and in neighboring territories. [ABSTRACT FROM AUTHOR]

The leaf area index (LAI) is a crucial metric for indicating crop development in the field, essential for both research and the practical implementation of precision agriculture. Unmanned aerial vehicles (UAVs) are widely used for monitoring crop growth due to their rapid, repetitive capture ability and cost-effectiveness. Therefore, we developed a non-destructive monitoring method for peanut LAI, combining UAV vegetation indices (VI) and texture features (TF). Field experiments were conducted to capture multispectral imagery of peanut crops. Based on these data, an optimal regression model was constructed to estimate LAI. The initial computation involves determining the potential spectral and textural characteristics. Subsequently, a comprehensive correlation study between these features and peanut LAI is conducted using Pearson's product component correlation and recursive feature elimination. Six regression models, including univariate linear regression, support vector regression, ridge regression, decision tree regression, partial least squares regression, and random forest regression, are used to determine the optimal LAI estimation. The following results are observed: (1) Vegetation indices exhibit greater correlation with LAI than texture characteristics. (2) The choice of GLCM parameters for texture features impacts estimation accuracy. Generally, smaller moving window sizes and higher grayscale quantization levels yield more accurate peanut LAI estimations. (3) The SVR model using both VI and TF offers the utmost precision, significantly improving accuracy (R2 = 0.867, RMSE = 0.491). Combining VI and TF enhances LAI estimation by 0.055 (VI) and 0.541 (TF), reducing RMSE by 0.093 (VI) and 0.616 (TF). The findings highlight the significant improvement in peanut LAI estimation accuracy achieved by integrating spectral and textural characteristics with appropriate parameters. These insights offer valuable guidance for monitoring peanut growth. [ABSTRACT FROM AUTHOR]

The relative content of chlorophyll, assessed through the soil and plant analyzer development (SPAD), serves as a reliable indicator reflecting crop photosynthesis and the nutritional status during crop growth and development. In this study, we employed machine learning methods utilizing unmanned aerial vehicle (UAV) multi-spectrum remote sensing to predict the SPAD value of litchi fruit. Input features consisted of various vegetation indices and texture features during distinct growth periods, and to streamline the feature set, the full subset regression algorithm was applied for dimensionality reduction. Our findings revealed the superiority of stacking models over individual models. During the litchi fruit development period, the stacking model, incorporating vegetation indices and texture features, demonstrated a validation set coefficient of determination (R2) of 0.94, a root mean square error (RMSE) of 2.4, and a relative percent deviation (RPD) of 3.0. Similarly, in the combined litchi growing period and autumn shoot period, the optimal model for estimating litchi SPAD was the stacking model based on vegetation indices and texture features, yielding a validation set R2, RMSE, and RPD of 0.84, 3.9, and 1.9, respectively. This study furnishes data support for the precise estimation of litchi SPAD across different periods through varied combinations of independent variables. [ABSTRACT FROM AUTHOR]

In high groundwater level mining areas, subsidence resulting from mining can lead to waterlogging in farmland, causing damage to crops and affecting their growth and development, thereby affecting regional food security. Therefore, it is necessary to restore agricultural production in the coal mining subsidence water areas in the densely populated eastern plains. This study focuses on the Yongcheng coal mining subsidence water areas. It utilizes Sentinel-1 and Sentinel-2 data from May to October in the years 2019 to 2022 to monitor the growth and development of crops. The results demonstrated that (1) the accuracy of aquatic crops categorization was improved by adjusting the elevation of the study region with Mining Subsidence Prediction Software (MSPS 1.0). The order of accuracy for classifying aquatic crops using different machine learning techniques is Random Forest (RF) > Classification and Regression Trees (CART) ≥ Support Vector Machine (SVM). Using the RF method, the obtained classification results can be used for subsequent crop growth monitoring. (2) During the early stages of crop growth, when vegetation cover is low, the Radar Vegetation Index (RVI) is sensitive to the volume scattering of crops, making it suitable for tracking the early growth processes of crops. The peak RVI values for crops from May to July are ranked in the following order: rice (2.595), euryale (2.590), corn (2.535), and lotus (2.483). (3) The order of crops showing improved growth conditions during the mid-growth stage is as follows: rice (47.4%), euryale (43.4%), lotus (27.6%), and corn (4.01%). This study demonstrates that in the Yongcheng coal subsidence water areas, the agricultural reclamation results for the grain-focused model with rice as the main crop and the medicinal herb-focused model with euryale as the main crop are significant. This study can serve as a reference for agricultural management and land reclamation efforts in other coal subsidence water areas. [ABSTRACT FROM AUTHOR]

Due to the expanding population and the constantly changing climate, food production is now considered a crucial concern. Although passive satellite remote sensing has already demonstrated its capabilities in accurate crop development monitoring, its limitations related to sunlight and cloud cover significantly restrict real-time temporal monitoring resolution. Considering synthetic aperture radar (SAR) technology, which is independent of the Sun and clouds, SAR remote sensing can be a perfect alternative to passive remote sensing methods. However, a variety of SAR sensors and delivered SAR indices present different performances in such context for different vegetation species. Therefore, this work focuses on comparing various SAR-derived indices from C-band and (Sentinel-1) and X-band (TerraSAR-X) data with the in situ information (phenp; pgy development, vegetation height and soil moisture) in the context of tracking the phenological development of corn, winter wheat, rye, canola, and potato. For this purpose, backscattering coefficients in VV and VH polarizations (σ V V 0 , σ V H 0 ), interferometric coherence, and the dual pol radar vegetation index (DpRVI) were calculated. To reduce noise in time series data and evaluate which filtering method presents a higher usability in SAR phenology tracking, signal filtering, such as Savitzky–Golay and moving average, with different parameters, were employed. The achieved results present that, for various plant species, different sensors (Sentinel-1 or TerraSAR-X) represent different performances. For instance, σ V H 0 of TerraSAR-X offered higher consistency with corn development (r = 0.81), while for canola σ V H 0 of Sentinel-1 offered higher performance (r = 0.88). Generally, σ V V 0 , σ V H 0 performed better than DpRVI or interferometric coherence. Time series filtering makes it possible to increase an agreement between phenology development and SAR-delivered indices; however, the Savitzky–Golay filtering method is more recommended. Besides phenological development, high correspondences can be found between vegetation height and some of SAR indices. Moreover, in some cases, moderate correlation was found between SAR indices and soil moisture. [ABSTRACT FROM AUTHOR]

Ecotourism and agricultural development have been proven to have synergistic effects, although few studies have employed a spatial planning approach to incorporate tourism growth into crop cultivation planning. This study constructed a theoretical framework of environmental suitability—farmland accessibility—tourist's landscape preferences for crop cultivation planning to link regional agriculture and ecotourism development. The spatial planning of rapeseed cultivation in Qinghai Province was chosen as a case study. The main research methods include an environmental suitability analysis based on remote sensing and Maxent modeling, a farmland accessibility analysis based on a GIS platform, and a landscape preference questionnaire survey of tourists. According to the survey's findings, almost 80% of tourists thought rapeseed flowers enhanced the beauty of natural landscapes. This demonstrated the enormous potential of rapeseed fields for fostering ecotourism. Based on environmental factors, the optimum region for rapeseed cultivation covered 5.38% of the study area, or roughly 6327 km2. The comprehensive optimum zone, which encompassed both agricultural accessibility and environmental suitability, was equal to 12.63% of the study area's farming area, or around 929 km2. This study's crop cultivation suitability framework can integrate agricultural and ecotourism development, with substantial implications for achieving coordinated economic, social, and environmental development. [ABSTRACT FROM AUTHOR]

Lodging is one of the typical abiotic adversities during paddy rice growth. In addition to affecting photosynthesis, it can seriously damage crop growth and development, such as reducing rice quality and hindering automated harvesting. It is, therefore, imperative to accurately and in good time acquire crop-lodging areas for yield prediction, agricultural insurance claims, and disaster-management decisions. However, the accuracy requirements for crop-lodging monitoring remain challenging due to complicated impact factors. Aiming at identifying paddy rice lodging on Shazai Island, Guangdong, China, caused by heavy rainfall and strong wind, a decision-tree model was constructed using multiple-parameter information from Sentinel-1 SAR images and the in situ lodging samples. The model innovatively combined the five backscattering coefficients with five polarization decomposition parameters and quantified the importance of each parameter feature. It was found that the decision-tree method coupled with polarization decomposition can be used to obtain an accurate distribution of paddy rice-lodging areas. The results showed that: (1) Radar parameters can capture the changes in lodged paddy rice. The radar parameters that best distinguish paddy rice lodging are VV, VV+VH, VH/VV, and Span. (2) Span is the parameter with the strongest feature importance, which shows the necessity of adding polarization parameters to the classification model. (3) The dual-polarized Sentinel-1 database classification model can effectively extract the area of lodging paddy rice with an overall accuracy of 84.38%, and a total area precision of 93.18%. These observations can guide the future use of SAR-based information for crop-lodging assessment and post-disaster management. [ABSTRACT FROM AUTHOR]

Precision Agriculture (PA) is an approach to maximizing crop productivity in a sustainable manner. PA requires up-to-date, accurate and georeferenced information on crops, which can be collected from different sensors from ground, aerial or satellite platforms. The use of optical and thermal sensors from Unmanned Aerial Vehicle (UAV) platform is an emerging solution for mapping and monitoring in PA, yet many technological challenges are still open. This technical note discusses the choice of UAV type and its scientific payload for surveying a sample area of 5 hectares, as well as the procedures for replicating the study on a larger scale. This case study is an ideal opportunity to test the best practices to combine the requirements of PA surveys with the limitations imposed by local UAV regulations. In the field area, to follow crop development at various stages, nine flights over a period of four months were planned and executed. The usage of ground control points for optimal georeferencing and accurate alignment of maps created by multi-temporal processing is analyzed. Output maps are produced in both visible and thermal bands, after appropriate strip alignment, mosaicking, sensor calibration, and processing with Structure from Motion techniques. The discussion of strategies, checklists, workflow, and processing is backed by data from more than 5000 optical and radiometric thermal images taken during five hours of flight time in nine flights throughout the crop season. The geomatics challenges of a georeferenced survey for PA using UAVs are the key focus of this technical note. Accurate maps derived from these multi-temporal and multi-sensor surveys feed Geographic Information Systems (GIS) and Decision Support Systems (DSS) to benefit PA in a multidisciplinary approach. [ABSTRACT FROM AUTHOR]

Agricultural management decision-making in salinization-prone environments requires efficient soil salinity monitoring methods. This is the case in the B-XII irrigation district in SW Spain, a heavy clay reclaimed marsh area where a shallow saline water table and intensively irrigated agriculture create a fragile balance between salt accumulation and leaching in the root zone, which might be disrupted by the introduction of new crops and increasing climate variability. We evaluated the potential of electromagnetic induction (EMI) tomography for field-scale soil salinity assessment in this hyper-conductive environment, using EMI and limited analytical soil data measured in 2017 and 2020 under a processing tomato–cotton–sugar beet crop rotation. Salinity effects on crop development were assessed by comparing Sentinel 2 NDVI imagery with inverted depth-specific electrical conductivity (EC). Average apparent electrical conductivity (ECa) for the 1-m depth signal was 20% smaller in 2020 than in 2017, although the spatial ECa pattern was similar for both years. Inverted depth-specific EC showed a strong correlation (R ≈ 0.90) with saturated paste extract EC (ECe), [Na+] and sodium absorption ratio (SAR), resulting in linear calibration equations with R2 ≈ 0.8 for both years and leave-one-out cross validation Nash–Sutcliffe Efficiency Coefficient, ranging from 0.57 to 0.74. Overall, the chemical parameter estimation improved with depth and soil wetness (2017), yielding 0.83 < R <0.98 at 0.9 m. The observed spatial EC distributions showed a steadily increasing inverse correlation with NDVI during the growing season, particularly for processing tomato and cotton, reaching R values of −0.71 and −0.85, respectively. These results confirm the potential of EMI tomography for mapping and monitoring soil salinity in the B-XII irrigation district, while it allows, in combination with NDVI imagery, a detailed spatial assessment of soil salinity impacts on crop development throughout the growing season. Contrary to the popular belief among farmers in the area, and despite non-saline topsoil conditions, spatial EC and subsoil salinity patterns were found to affect crop development negatively in the studied field. [ABSTRACT FROM AUTHOR]

Soil moisture content (SMC) is an indispensable basic element for crop growth and development in agricultural production. Obtaining accurate information on SMC in real time over large agricultural areas has important guiding significance for crop yield estimation and production management. In this study, the paper reports on the retrieval of SMC from RADARSAT-2 polarimetric SAR data. The proposed SMC retrieval algorithm includes vegetation correction based on a ratio method and roughness correction based on the optimal roughness method. Three vegetation description parameters (i.e., RVI, LAI, and NDVI) serve as vegetation descriptors in the parameterization of the algorithm. To testify the vegetation correction result of the algorithm, the water cloud model (WCM) was compared with the algorithm. The calibrated integrated equation model (CIEM) was employed to describe the backscattering from the underlying soil. To improve the accuracy of SMC retrieval, the CIEM model was optimized by using the optimal roughness parameter and the normalization method of reference incidence angle. Validation against ground measurements showed a high correlation between the measured and estimated SMC when the NDVI serves as vegetation descriptor (R2 = 0.68, RMSE = 4.15 vol.%, p < 0.01). The overall estimation performance of the proposed SMC retrieval algorithm is better than that of the WCM. It demonstrates that the proposed algorithm has an operational potential to estimate SMC over wheat fields during the growing season. [ABSTRACT FROM AUTHOR]

Drought is a major natural hazard that impacts agriculture, the environment, and socio-economic conditions. In 2018 and 2019, Europe experienced a severe drought due to below average precipitation and high temperatures. Drought stress affects the moisture content and structure of agricultural crops and can result in lower yields. Synthetic Aperture Radar (SAR) observations are sensitive to the dielectric and geometric characteristics of crops and underlying soils. This study uses data from ESA's Sentinel-1 SAR satellite to investigate the influence of drought stress on major arable crops of the Netherlands, its regional variability and the impact of water management decisions on crop development. Sentinel-1 VV, VH and VH/VV backscatter data are used to quantify the variability in the spatio-temporal dynamics of agricultural crop parcels in response to drought. Results show that VV and VH backscatter values are 1 to 2 dB lower for crop parcels during the 2018 drought compared to values in 2017. In addition, the growth season indicated by the cross-ratio (CR, VH/VV) for maize and onion is shorter during the drought year. Differences due to irrigation restrictions are observed in backscatter response from maize parcels. Lower CR values in 2019 indicate the impact of drought on the start of the growing season. Results demonstrate that Sentinel-1 can detect changes in the seasonal cycle of arable crops in response to agricultural drought. [ABSTRACT FROM AUTHOR]

High-frequency Earth observation (EO) data have been shown to be effective in identifying crops and monitoring their development. The purpose of this paper is to derive quantitative indicators of crop productivity using synthetic aperture radar (SAR). This study shows that the field-specific SAR time series can be used to characterise growth and maturation periods and to estimate the performance of cereals. Winter wheat fields on the Rothamsted Research farm in Harpenden (UK) were selected for the analysis during three cropping seasons (2017 to 2019). Average SAR backscatter from Sentinel-1 satellites was extracted for each field and temporal analysis was applied to the backscatter cross-polarisation ratio (VH/VV). The calculation of the different curve parameters during the growing period involves (i) fitting of two logistic curves to the dynamics of the SAR time series, which describe timing and intensity of growth and maturation, respectively; (ii) plotting the associated first and second derivative in order to assist the determination of key stages in the crop development; and (iii) exploring the correlation matrix for the derived indicators and their predictive power for yield. The results show that the day of the year of the maximum VH/VV value was negatively correlated with yield (r = −0.56), and the duration of “full” vegetation was positively correlated with yield (r = 0.61). Significant seasonal variation in the timing of peak vegetation (p = 0.042), the midpoint of growth (p = 0.037), the duration of the growing season (p = 0.039) and yield (p = 0.016) were observed and were consistent with observations of crop phenology. Further research is required to obtain a more detailed picture of the uncertainty of the presented novel methodology, as well as its validity across a wider range of agroecosystems.

Remote sensing observations that vary in response to plant growth and senescence can be used to monitor crop development within and across growing seasons. Identifying when crops reach specific growth stages can improve harvest yield prediction and quantify climate change. Using the Level 2 vegetation optical depth (VOD) product from the European Space Agency’s Soil Moisture and Ocean Salinity (SMOS) satellite, we retrospectively estimate the timing of a key crop development stage in the United States Corn Belt. We employ nonlinear curves nested within a hierarchical modeling framework to extract the timing of the third reproductive development stage of corn (R3) as well as other new agronomic signals from SMOS VOD. We compare our estimates of the timing of R3 to United States Department of Agriculture (USDA) survey data for the years 2011, 2012, and 2013. We find that 87%, 70%, and 37%, respectively, of our model estimates of R3 timing agree with USDA district-level observations. We postulate that since the satellite estimates can be directly linked to a physiological state (the maximum amount of plant water, or water contained within plant tissue per ground area) it is more accurate than the USDA data which is based upon visual observations from roadways. Consequently, SMOS VOD could be used to replace, at a finer resolution than the district-level USDA reports, the R3 data that has not been reported by the USDA since 2013. We hypothesize the other model parameters contain new information about soil and crop management and crop productivity that are not routinely collected by any federal or state agency in the Corn Belt.

Crop phenology is an important parameter for crop growth monitoring, yield prediction, and growth simulation. The dynamic threshold method is widely used to retrieve vegetation phenology from remotely sensed vegetation index time series. However, crop growth is not only driven by natural conditions, but also modified through field management activities. Complicated planting patterns, such as multiple cropping, makes the vegetation index dynamics less symmetrical. These impacts are not considered in current approaches for crop phenology retrieval based on the dynamic threshold method. Thus, this paper aimed to (1) investigate the optimal thresholds for retrieving the start of the season (SOS) and the end of the season (EOS) of different crops, and (2) compare the performances of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) in retrieving crop phenology with a modified version of the dynamic threshold method. The reference data included SOS and EOS ground observations of three major crop types in 2015 and 2016, which includes rice, wheat, and maize. Results show that (1) the modification of the original method ensures a 100% retrieval rate, which was not guaranteed using the original method. The modified dynamic threshold method is more suitable to retrieve crop SOS/EOS because it considers the asymmetry of crop vegetation index time series. (2) It is inappropriate to retrieve SOS and EOS with the same threshold for all crops, and the commonly used 20% or 50% thresholds are not the optimal thresholds for all crops. (3) For single and late rice, the accuracies of the SOS estimations based on EVI are generally higher compared to those based on NDVI. However, for spring maize and summer maize, results based on NDVI give higher accuracies. In terms of EOS, for early rice and summer maize, estimates based on EVI result in higher accuracies, but, for late rice and winter wheat, results based on NDVI are closer to the ground records.

Non-destructive monitoring of crop development is of key interest for agronomy and crop breeding. Crop Surface Models (CSMs) representing the absolute height of the plant canopy are a tool for this. In this study, fresh and dry barley biomass per plot are estimated from CSM-derived plot-wise plant heights. The CSMs are generated in a semi-automated manner using Structure-from-Motion (SfM)/Multi-View-Stereo (MVS) software from oblique stereo RGB images. The images were acquired automatedly from consumer grade smart cameras mounted at an elevated position on a lifting hoist. Fresh and dry biomass were measured destructively at four dates each in 2014 and 2015. We used exponential and simple linear regression based on different calibration/validation splits. Coefficients of determination R² between 0.55 and 0.79 and root mean square errors (RMSE) between 97 and 234 g/m² are reached for the validation of predicted vs. observed dry biomass, while Willmott's refined index of model performance dr ranges between 0.59 and 0.77. For fresh biomass, R² values between 0.34 and 0.61 are reached, with root mean square errors (RMSEs) between 312 and 785 g/m² and dr between 0.39 and 0.66. We therefore established the possibility of using this novel low-cost system to estimate barley dry biomass over time. [ABSTRACT FROM AUTHOR]

Oilseed rape (Brassica napus L.) is one of the three most important oil crops in China, and is regarded as a drought-tolerant oilseed crop. However, it is commonly sensitive to waterlogging, which usually refers to an adverse environment that limits crop development. Moreover, crop growth and soil irrigation can be monitored at a regional level using remote sensing data. High spatial resolution optical satellite sensors are very useful to capture and resist unfavorable field conditions at the sub-field scale. In this study, four different optical sensors, i.e., Pleiades-1A, Worldview-2, Worldview-3, and SPOT-6, were used to estimate the dry above-ground biomass (AGB) of oilseed rape and track the seasonal growth dynamics. In addition, three different soil water content field experiments were carried out at different oilseed rape growth stages from November 2014 to May 2015 in Northern Zhejiang province, China. As a significant indicator of crop productivity, AGB was measured during the seasonal growth stages of the oilseed rape at the experimental plots. Several representative vegetation indices (VIs) obtained from multiple satellite sensors were compared with the simultaneously-collected oilseed rape AGB. Results showed that the estimation model using the normalized difference vegetation index (NDVI) with a power regression model performed best through the seasonal growth dynamics, with the highest coefficient of determination (R2 = 0.77), the smallest root mean square error (RMSE = 104.64 g/m2), and the relative RMSE (rRMSE = 21%). It is concluded that the use of selected VIs and high spatial multiple satellite data can significantly estimate AGB during the winter oilseed rape growth stages, and can be applied to map the variability of winter oilseed rape at the sub-field level under different waterlogging conditions, which is very promising in the application of agricultural irrigation and precision agriculture. [ABSTRACT FROM AUTHOR]

Early warning systems for food security require accurate and up-to-date information on the location of major crops in order to prevent hazards. A recent systematic analysis of existing cropland maps identified priority areas for cropland mapping and highlighted a major need for the Sahelian and Sudanian agrosystems. This paper proposes a knowledge-based approach to map cropland in the Sahelian and Sudanian agrosystems that benefits from the 100-m spatial resolution of the recent PROBA-V sensor. The methodology uses five temporal features characterizing crop development throughout the vegetative season to optimize cropland discrimination. A feature importance analysis validates the efficiency of using a diversity of temporal features. The fully-automated method offers the first cropland map at 100-m using the PROBA-V sensor with an overall accuracy of 84% and an F-score for the cropland class of 74%. The improvements observed compared to existing cropland products are related to the hectometric resolution, to the methodology and to the quality of the labeling layer from which reliable training samples were automatically extracted. Classification errors are mainly explained by data availability and landscape fragmentation. Further improvements are expected with the upcoming enhanced cloud screening of the PROBA-V sensor. [ABSTRACT FROM AUTHOR]

Timely crop yield forecasts at a national level are substantial to support food policies, to assess agricultural production, and to subsidize regions affected by food shortage. This study presents an operational crop yield forecasting system for Poland that employs freely available satellite and agro-meteorological products provided by the Copernicus programme. The crop yield predictors consist of: (1) Vegetation condition indicators provided daily by Sentinel-3 OLCI (optical) and SLSTR (thermal) imagery, (2) a backward extension of Sentinel-3 data (before 2018) derived from cross-calibrated MODIS data, and (3) air temperature, total precipitation, surface radiation, and soil moisture derived from ERA-5 climate reanalysis generated by the European Centre for Medium-Range Weather Forecasts. The crop yield forecasting algorithm is based on thermal time (growing degree days derived from ERA-5 data) to better follow the crop development stage. The recursive feature elimination is used to derive an optimal set of predictors for each administrative unit, which are ultimately employed by the Extreme Gradient Boosting regressor to forecast yields using official yield statistics as a reference. According to intensive leave-one-year-out cross validation for the 2000–2019 period, the relative RMSE for voivodships (NUTS-2) are: 8% for winter wheat, and 13% for winter rapeseed and maize. Respectively, for municipalities (LAU) it equals 14% for winter wheat, 19% for winter rapeseed, and 27% for maize. The system is designed to be easily applicable in other regions and to be easily adaptable to cloud computing environments such as Data and Information Access Services (DIAS) or Amazon AWS, where data sets from the Copernicus programme are directly accessible. [ABSTRACT FROM AUTHOR]

Global crop mapping and monitoring requires high-resolution spatio-temporal information. In this regard, dual polarimetric Synthetic Aperture Radar (SAR) sensors provide high temporal and high spatial resolutions with large swath width. Generally, crop phenological development studies utilized SAR backscatter intensity-based descriptors. However, these descriptors are derived either from the covariance matrix elements or from the eigendecomposition. Therefore, this approach fails to utilize the complete polarization information of the scattered wave. In this study, we propose a target characterization parameter, θ xP that utilizes the 2D Barakat degree of polarization and the elements of the covariance matrix. We also propose an unsupervised clustering scheme using θ xP and the scattering entropy, H xP . We utilize time-series Sentinel-1 data of canola and wheat fields over a Canadian test site to show the sensitivity of θ xP to the development of crop morphology at different phenological stages. During the initial growth stages, θ xP values are low due to the low vegetation density. In contrast, at advanced phenological stages, we observe decreased values of θ xP due to the appearance of complex canopy structure. Similarly, the effectiveness of the unsupervised H xP / θ xP clustering plane is also evident from the temporal clustering plots. This innovative clustering framework is beneficial for the operational use of Sentinel-1 SAR data for agricultural applications. [ABSTRACT FROM AUTHOR]

The net primary productivity (NPP) and aboveground biomass mapping of crops based on remote sensing technology are not only conducive to understanding the growth and development of crops but can also be used to monitor timely agricultural information, thereby providing effective decision making for agricultural production management. To solve the saturation problem of the NDVI in the aboveground biomass mapping of crops, the original CASA model was improved using narrow-band red-edge information, which is sensitive to vegetation chlorophyll variation, and the fraction of photosynthetically active radiation (FPAR), NPP, and aboveground biomass of winter wheat and maize were mapped in the main growing seasons. Moreover, in this study, we deeply analyzed the seasonal change trends of crops' biophysical parameters in terms of the NDVI, FPAR, actual light use efficiency (LUE), and their influence on aboveground biomass. Finally, to analyze the uncertainty of the aboveground biomass mapping of crops, we further discussed the inversion differences of FPAR with different vegetation indices. The results demonstrated that the inversion accuracies of the FPAR of the red-edge normalized vegetation index ( NDVI red - edge ) and red-edge simple ratio vegetation index ( SR red - edge ) were higher than those of the original CASA model. Compared with the reference data, the accuracy of aboveground biomass estimated by the improved CASA model was 0.73 and 0.70, respectively, which was 0.21 and 0.13 higher than that of the original CASA model. In addition, the analysis of the FPAR inversions of different vegetation indices showed that the inversion accuracies of the red-edge vegetation indices NDVI red - edge and SR red - edge were higher than those of the other vegetation indices, which confirmed that the vegetation indices involving red-edge information can more effectively retrieve FPAR and aboveground biomass of crops. [ABSTRACT FROM AUTHOR]

Various drought indices have been used for agricultural drought monitoring, such as Standardized Precipitation Index (SPI), Standardized Precipitation Evapotranspiration Index (SPEI), Palmer Drought Severity Index (PDSI), Soil Water Deficit Index (SWDI), Normalized Difference Vegetation Index (NDVI), Vegetation Health Index (VHI), Vegetation Drought Response Index (VegDRI), and Scaled Drought Condition Index (SDCI). They incorporate such factors as rainfall, land surface temperature (LST), potential evapotranspiration (PET), soil moisture content (SM), and vegetation index to express the meteorological and agricultural aspects of drought. However, these five factors should be combined more comprehensively and reasonably to explain better the dryness/wetness of land surface and the association with crop yield. This study aims to develop the Integrated Crop Drought Index (ICDI) by combining the weather factors (rainfall and LST), hydrological factors (PET and SM), and a vegetation factor (enhanced vegetation index (EVI)) to better express the wet/dry state of land surface and healthy/unhealthy state of vegetation together. The study area was the State of Illinois, a key region of the U.S. Corn Belt, and the quantification and analysis of the droughts were conducted on a county scale for 2004–2019. The performance of the ICDI was evaluated through the comparisons with SDCI and VegDRI, which are the representative drought index in terms of the composite of the dryness and vegetation elements. The ICDI properly expressed both the dry and wet trend of the land surface and described the state of the agricultural drought accompanied by yield damage. The ICDI had higher positive correlations with the corn yields than SDCI and VegDRI during the crucial growth period from June to August for 2004–2019, which means that the ICDI could reflect the agricultural drought well in terms of the dryness/wetness of land surface and the association with crop yield. Future work should examine the other factors for ICDI, such as locality, crop type, and the anthropogenic impacts, on drought. It is expected that the ICDI can be a viable option for agricultural drought monitoring and yield management. [ABSTRACT FROM AUTHOR]

The time series of synthetic aperture radar (SAR) data are commonly and successfully used to monitor the biophysical parameters of agricultural fields. Because, until now, mainly backscatter coefficients have been analysed, this study examines the potentials of entropy, anisotropy, and alpha angle derived from a dual-polarimetric decomposition of Sentinel-1 data to monitor crop development. The temporal profiles of these parameters are analysed for wheat and barley in the vegetation periods 2017 and 2018 for 13 fields in two test sites in Northeast Germany. The relation between polarimetric parameters and biophysical parameters observed in the field is investigated using linear and exponential regression models that are evaluated using the coefficient of determination (R2) and the root mean square error (RMSE). The performance of single regression models is furthermore compared to those of multiple regression models, including backscatter coefficients in VV and VH polarisation as well as polarimetric decomposition parameters entropy and alpha. Characteristic temporal profiles of entropy, anisotropy, and alpha reflecting the main phenological changes in plants as well as the meteorological differences between the two years are observed for both crop types. The regression models perform best for data from the phenological growth stages tillering to booting. The highest R2 values of the single regression models are reached for the plant height of wheat related to entropy and anisotropy with R2 values of 0.64 and 0.61, respectively. The multiple regression models of VH, VV, entropy, and alpha outperform single regression models in most cases. R2 values of multiple regression models of plant height (0.76), wet biomass (0.7), dry biomass (0.7), and vegetation water content (0.69) improve those of single regression models slightly by up to 0.05. Additionally, the RMSE values of the multiple regression models are around 10% lower compared to those of single regression models. The results indicate the capability of dual-polarimetric decomposition parameters in serving as meaningful input parameters for multiple regression models to improve the prediction of biophysical parameters. Additionally, their temporal profiles indicate phenological development dependent on meteorological conditions. Knowledge about biophysical parameter development and phenology is important for farmers to monitor crop growth variability during the vegetation period to adapt and to optimize field management. [ABSTRACT FROM AUTHOR]

Remote sensing data play a crucial role in monitoring crop dynamics in the context of precision agriculture by characterizing the spatial and temporal variability of crop traits. At present there is special interest in assessing the long-term impacts of biochar in agro-ecosystems. Despite the growing body of literature on monitoring the potential biochar effects on harvested crop yield and aboveground productivity, studies focusing on the detailed crop performance as a consequence of long-term biochar enrichment are still lacking. The primary objective of this research was to evaluate crop performance based on high-resolution unmanned aerial vehicle (UAV) imagery considering both crop growth and health through RGB and multispectral analysis, respectively. More specifically, this approach allowed monitoring of century-old biochar impacts on winter wheat crop performance. Seven Red-Green-Blue (RGB) and six multispectral flights were executed over 11 century-old biochar patches of a cultivated field. UAV-based RGB imagery exhibited a significant positive impact of century-old biochar on the evolution of winter wheat canopy cover (p-value = 0.00007). Multispectral optimized soil adjusted vegetation index indicated a better crop development over the century-old biochar plots at the beginning of the season (p-values < 0.01), while there was no impact towards the end of the season. Plant height, derived from the RGB imagery, was slightly higher for century-old biochar plots. Crop health maps were computed based on principal component analysis and k-means clustering. To our knowledge, this is the first attempt to quantify century-old biochar effects on crop performance during the entire growing period using remotely sensed data. Ground-based measurements illustrated a significant positive impact of century-old biochar on crop growth stages (p-value of 0.01265), whereas the harvested crop yield was not affected. Multispectral simplified canopy chlorophyll content index and normalized difference red edge index were found to be good linear estimators of harvested crop yield (p-value(Kendall) of 0.001 and 0.0008, respectively). The present research highlights that other factors (e.g., inherent pedological variations) are of higher importance than the presence of century-old biochar in determining crop health and yield variability. [ABSTRACT FROM AUTHOR]

High-frequency Earth observation (EO) data have been shown to be effective in identifying crops and monitoring their development. The purpose of this paper is to derive quantitative indicators of crop productivity using synthetic aperture radar (SAR). This study shows that the field-specific SAR time series can be used to characterise growth and maturation periods and to estimate the performance of cereals. Winter wheat fields on the Rothamsted Research farm in Harpenden (UK) were selected for the analysis during three cropping seasons (2017 to 2019). Average SAR backscatter from Sentinel-1 satellites was extracted for each field and temporal analysis was applied to the backscatter cross-polarisation ratio (VH/VV). The calculation of the different curve parameters during the growing period involves (i) fitting of two logistic curves to the dynamics of the SAR time series, which describe timing and intensity of growth and maturation, respectively; (ii) plotting the associated first and second derivative in order to assist the determination of key stages in the crop development; and (iii) exploring the correlation matrix for the derived indicators and their predictive power for yield. The results show that the day of the year of the maximum VH/VV value was negatively correlated with yield (r = −0.56), and the duration of "full" vegetation was positively correlated with yield (r = 0.61). Significant seasonal variation in the timing of peak vegetation (p = 0.042), the midpoint of growth (p = 0.037), the duration of the growing season (p = 0.039) and yield (p = 0.016) were observed and were consistent with observations of crop phenology. Further research is required to obtain a more detailed picture of the uncertainty of the presented novel methodology, as well as its validity across a wider range of agroecosystems. [ABSTRACT FROM AUTHOR]

Winter wheat is one of the major cereal crops in China. The spatial distribution of winter wheat planting areas is closely related to food security; however, mapping winter wheat with time-series finer spatial resolution satellite images across large areas is challenging. This paper explores the potential of combining temporally aggregated Landsat-8 OLI and Sentinel-2 MSI data available via the Google Earth Engine (GEE) platform for mapping winter wheat in Shandong Province, China. First, six phenological median composites of Landsat-8 OLI and Sentinel-2 MSI reflectance measures were generated by a temporal aggregation technique according to the winter wheat phenological calendar, which covered seedling, tillering, over-wintering, reviving, jointing-heading and maturing phases, respectively. Then, Random Forest (RF) classifier was used to classify multi-temporal composites but also mono-temporal winter wheat development phases and mono-sensor data. The results showed that winter wheat could be classified with an overall accuracy of 93.4% and F1 measure (the harmonic mean of producer's and user's accuracy) of 0.97 with temporally aggregated Landsat-8 and Sentinel-2 data were combined. As our results also revealed, it was always good to classify multi-temporal images compared to mono-temporal imagery (the overall accuracy dropped from 93.4% to as low as 76.4%). It was also good to classify Landsat-8 OLI and Sentinel-2 MSI imagery combined instead of classifying them individually. The analysis showed among the mono-temporal winter wheat development phases that the maturing phase's and reviving phase's data were more important than the data for other mono-temporal winter wheat development phases. In sum, this study confirmed the importance of using temporally aggregated Landsat-8 OLI and Sentinel-2 MSI data combined and identified key winter wheat development phases for accurate winter wheat classification. These results can be useful to benefit on freely available optical satellite data (Landsat-8 OLI and Sentinel-2 MSI) and prioritize key winter wheat development phases for accurate mapping winter wheat planting areas across China and elsewhere. [ABSTRACT FROM AUTHOR]

Proximal optical sensors (POSs) are effective devices for monitoring the development of crops and the nitrogen (N) status of plants. POSs are both useful and necessary in facilitating the reduction of N losses into the environment and in attaining higher nitrogen use efficiency (NUE). To date, no comparison of these instruments has been made on quinoa. A field experiment conducted in Tuscany, Italy, with different POSs, has assessed the development of quinoa with respect to N status. Three sets of POSs were used (SPAD-502, GreenSeeker, and Canopeo App.) to monitor quinoa development and growth under different types of fertilizers (digestate and urea) and levels of N fertilization (100, 50, and 0 kg N ha−1). The present findings showed that in-season predictions of crop biomass at harvest by SPAD-502 and GreenSeeker optical sensors were successful in terms of the coefficient of determination (R2 = 0.68 and 0.82, respectively) and statistical significance (p < 0.05), while the Canopeo App. was suitable for monitoring the plant´s canopy expansion and senescence. The relative error (RE%) showed a remarkably high performance between observed and predicted values, 5.80% and 4.12% for GreenSeeker and SPAD-502, respectively. Overall, the POSs were effective devices for monitoring quinoa development during the growing season and for predicting dry biomass at harvest. However, abiotic stresses (e.g., heat-stress conditions at flowering) were shown to reduce POSs' accuracy when estimating seed yields at harvest, and this problem will likely be overcome by advancing the sowing date. [ABSTRACT FROM AUTHOR]

The timely monitoring of crop disease development is very important for precision agriculture applications. Remote sensing-based vegetation indices (VIs) can be good indicators of crop disease severity, but current methods are mainly dependent on manual ground survey results. Based on VI normalization, an automated crop disease severity grading method without the use of ground surveys was proposed in this study. This technique was applied to two cotton fields infested with different levels of cotton root rot in south Texas in the United States, where airborne hyperspectral imagery was collected. Six typical VIs were calculated from the hyperspectral imagery and their histograms indicated that VI normalization could eliminate the influences of variable field conditions and the VI value range variations, allowing a potentially broader scope of application. According to the analysis of the obtained results from the spectral dimension, spatial dimension and descriptive statistics, the disease grading results were in general agreement with previous ground survey results, proving the validity of the disease severity grading method. Although satisfactory results could be achieved from different types of VI, there is still room for further improvement through the exploration of more VIs. With the advantages of independence of ground surveys and potential universal applicability, the newly proposed crop disease grading method will be of great significance for crop disease monitoring over large geographical areas. [ABSTRACT FROM AUTHOR]

Remote sensing observations that vary in response to plant growth and senescence can be used to monitor crop development within and across growing seasons. Identifying when crops reach specific growth stages can improve harvest yield prediction and quantify climate change. Using the Level 2 vegetation optical depth (VOD) product from the European Space Agency's Soil Moisture and Ocean Salinity (SMOS) satellite, we retrospectively estimate the timing of a key crop development stage in the United States Corn Belt. We employ nonlinear curves nested within a hierarchical modeling framework to extract the timing of the third reproductive development stage of corn (R3) as well as other new agronomic signals from SMOS VOD. We compare our estimates of the timing of R3 to United States Department of Agriculture (USDA) survey data for the years 2011, 2012, and 2013. We find that 87%, 70%, and 37%, respectively, of our model estimates of R3 timing agree with USDA district-level observations. We postulate that since the satellite estimates can be directly linked to a physiological state (the maximum amount of plant water, or water contained within plant tissue per ground area) it is more accurate than the USDA data which is based upon visual observations from roadways. Consequently, SMOS VOD could be used to replace, at a finer resolution than the district-level USDA reports, the R3 data that has not been reported by the USDA since 2013. We hypothesize the other model parameters contain new information about soil and crop management and crop productivity that are not routinely collected by any federal or state agency in the Corn Belt. [ABSTRACT FROM AUTHOR]

Crop phenology is an important parameter for crop growth monitoring, yield prediction, and growth simulation. The dynamic threshold method is widely used to retrieve vegetation phenology from remotely sensed vegetation index time series. However, crop growth is not only driven by natural conditions, but also modified through field management activities. Complicated planting patterns, such as multiple cropping, makes the vegetation index dynamics less symmetrical. These impacts are not considered in current approaches for crop phenology retrieval based on the dynamic threshold method. Thus, this paper aimed to (1) investigate the optimal thresholds for retrieving the start of the season (SOS) and the end of the season (EOS) of different crops, and (2) compare the performances of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) in retrieving crop phenology with a modified version of the dynamic threshold method. The reference data included SOS and EOS ground observations of three major crop types in 2015 and 2016, which includes rice, wheat, and maize. Results show that (1) the modification of the original method ensures a 100% retrieval rate, which was not guaranteed using the original method. The modified dynamic threshold method is more suitable to retrieve crop SOS/EOS because it considers the asymmetry of crop vegetation index time series. (2) It is inappropriate to retrieve SOS and EOS with the same threshold for all crops, and the commonly used 20% or 50% thresholds are not the optimal thresholds for all crops. (3) For single and late rice, the accuracies of the SOS estimations based on EVI are generally higher compared to those based on NDVI. However, for spring maize and summer maize, results based on NDVI give higher accuracies. In terms of EOS, for early rice and summer maize, estimates based on EVI result in higher accuracies, but, for late rice and winter wheat, results based on NDVI are closer to the ground records. [ABSTRACT FROM AUTHOR]

Agriculture is of huge economic significance in The Netherlands where the provision of real-time, reliable information on crop development is essential to support the transition towards precision agriculture. Optical imagery can provide invaluable insights into crop growth and development but is severely hampered by cloud cover. This case study in the Flevopolder illustrates the potential value of Sentinel-1 for monitoring five key crops in The Netherlands, namely sugar beet, potato, maize, wheat and English rye grass. Time series of radar backscatter from the European Space Agency's Sentinel-1 Mission are analyzed and compared to ground measurements including phenological stage and height. Temporal variations in backscatter data reflect changes in water content and structure associated with phenological development. Emergence and closure dates are estimated from the backscatter time series and validated against a photo archive. Coherence data are compared to Normalized Difference Vegetation Index (NDVI) and ground data, illustrating that the sudden increase in coherence is a useful indicator of harvest. The results presented here demonstrate that Sentinel-1 data have significant potential value to monitor growth and development of key Dutch crops. Furthermore, the guaranteed availability of Sentinel-1 imagery in clouded conditions ensures the reliability of data to meet the monitoring needs of farmers, food producers and regulatory bodies. [ABSTRACT FROM AUTHOR]

Few studies have focused on the potential impacts of topography on regional crop simulation, which might constrain the development of crop models and lead to inaccurate estimations for food security. In this study, we used remote sensing data to calibrate a regional crop model (MCWLA-Rice) for yield simulation in a double-rice crop rotation system in counties of Hunan province dominated by three landforms (plain, hill, and mountain). The calibration scheme with coarse remote sensing data (Global LAnd Surface Satellite, GLASS) greatly improved model accuracy for the double-rice system and is a promising method for yield estimation in large areas. The average improvement in relative root mean square error (RRMSE) was at most 48.00% for early rice and 41.25% for late rice. The average improvement in coefficient of determination (R2) value was at most 0.54 for early rice and 0.19 for late rice. Estimation of yield in counties dominated by different landform types indicated that: (1) MCWLA-Rice tended to be unstable in areas of complex topography and resulted in unbalanced proportions of overestimations and underestimations. (2) Differences in yield simulation between early rice and late rice varied among counties; yield estimates were highest in predominantly hilly counties, followed by counties dominated by plains, and lowest in predominantly mountainous counties. The results indicated that the topography might harm the accuracy of crop model simulations. Integration of topographic factors into crop models may enable yield estimation with enhanced accuracy to promote social development. [ABSTRACT FROM AUTHOR]