Crop NDVI Monitoring Based on Sentinel 1

Rice yield prediction through integration of biophysical parameters with SAR and optical remote sensing data using machine learning models

In an era marked by growing global population and climate variability, ensuring food security has become a paramount concern. Rice, being a staple crop for billions of people, requires accurate and timely yield prediction to ensure global food security. This study was undertaken across two rice crop seasons in the Udham Singh Nagar district of Uttarakhand state to predict rice yield at 45, 60 and 90 days after transplanting (DAT) through machine learning (ML) models, utilizing a combination of optical and Synthetic Aperture Radar (SAR) data in conjunction with crop biophysical parameters. Results revealed that the ML models were able to provide relatively accurate early yield estimates. For summer rice, eXtreme gradient boosting (XGB) was the best-performing model at all three stages (45, 60, and 90 DAT), while for kharif rice, the best-performing models at 45, 60, and 90 DAT were XGB, Neural network (NNET), and Cubist, respectively. The combined ranking of ML models showed that prediction accuracy improved as the prediction date approaches harvest, and the best prediction of yield was observed at 90 DAT for both summer and kharif rice. Overall rankings indicate that for summer rice, the top three models were XGB, NNET, and Support vector regression, while for kharif rice, these were Cubist, NNET, and Random Forest, respectively. The findings of this study offer valuable insights into the potential of the combined use of remote sensing and biophysical parameters using ML models, which enhances food security planning and resource management by enabling more informed decision-making by stakeholders such as farmers, policy planners as well as researchers.

Assessing impact of agroecological interventions in Niger through remotely sensed changes in vegetation

Water scarcity is a major challenge in the Sahel region of West Africa. Water scarcity in combination with prevalent soil degradation has severely reduced the land productivity in the region. The decrease in resiliency of food security systems of marginalized population has huge societal implications which often leads to mass migrations and conflicts. The U.S. Agency for International Development (USAID) and development organizations have made major investments in the Sahel to improve resilience through land rehabilitation activities in recent years. To help restore degraded lands at the farm level, the World Food Programme (WFP) with assistance from USAID's Bureau for Humanitarian Assistance supported the construction of water and soil retention structures called half-moons. The vegetation growing in the half-moons is vitally important to increase agricultural productivity and feed animals, a critical element of sustainable food security in the region. This paper investigates the effectiveness of interventions at 18 WFP sites in southern Niger using vegetative greenness observations from the Landsat 7 satellite. The pre - and post-intervention analysis shows that vegetation greenness after the half-moon intervention was nearly 50% higher than in the pre-intervention years. The vegetation in the intervened area was more than 25% greener than the nearby control area. Together, the results indicate that the half-moons are effective adaptations to the traditional land management systems to increase agricultural production in arid ecosystems, which is evident through improved vegetation conditions in southern Niger. The analysis shows that the improvement brought by the interventions continue to provide the benefits. Continued application of these adaptation techniques on a larger scale will increase agricultural production and build resilience to drought for subsistence farmers in West Africa. Quantifiable increase in efficacy of local-scale land and water management techniques, and the resulting jump in large-scale investments to scale similar efforts will help farmers enhance their resiliency in a sustainable manner will lead to a reduction in food security shortages.

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.

Crop classification by using dual-pol SAR vegetation indices derived from Sentinel-1 SAR-C data

In the following study, an attempt is made for crop classification of rainy season through analyzing time-series Sentinel-1 SAR data of May 2020 to September 2020. The SVI_DP index derived from dual-pol (VV and VH) bands consisting of NRPB ([Formula: see text]), DPDD [Formula: see text]), IDPDD ([Formula: see text]), and VDDPI [Formula: see text] ratios are utilized for discriminating inter-vegetative boundaries of crop pixels. This study was conducted near Karnal city region, Karnal district, Haryana, India. The Sentinel-1 data has the capability to penetrate thick cloud cover and provide high revisit frequency data for rain-fed crops. Obtained classification achieved higher accuracy in both RF (93.77%) and SVM (93.50%) classifiers. Obtained linear regression statistics of mean raster imagery reveals that IDPDD index is much sensitive to other crop which has highest standard deviations in σ_vh° and σ_vv° bands throughout the period, and high R² with σ_vh° (0.70), VV (0.58), NRPB (0.693), and DPDD (0.697) indices. In contrast to this, IDPDD index has least correlation (< 0.289) with σ_vh°, σ_vv°, EVI 2, NRPB, and DPDD indices for water body which has smooth surface and lowest SAR backscattering with minimum standard deviations in the same period.

Sentinel-1 to NDVI for Agricultural Fields Using Hyperlocal Dynamic Machine Learning Approach

The normalized difference vegetation index (NDVI) is a key parameter in precision agriculture. It has been used globally since the 1970s as a proxy to monitor crop growth and correlates to the crop coefficient (Kc), leaf area index (LAI), crop cover, and more. Yet, it is susceptible to clouds and other atmospheric conditions that might alter the crop’s real NDVI value. Synthetic Aperture Radar (SAR), on the other hand, can penetrate clouds and is hardly affected by atmospheric conditions, but it is sensitive to the physical structure of the crop and therefore does not give a direct indication of the NDVI. Several SAR indices and methods have been suggested to estimate NDVIs via SAR; however, they tend to work for local spatial and temporal conditions and do not work well globally. This is because they are not flexible enough to capture the changing NDVI–SAR relationship throughout the crop-growing season. This study suggests a new method for converting Sentinel-1 to NDVIs for Agricultural Fields (SNAF) by utilizing a hyperlocal machine learning approach. This method generates multiple on-the-fly disposal field- and time-specific models for every available Sentinel-1 image across 2021. Each model learns the field-specific NDVI (from Sentinel-2 and Landsat-8) –SAR (Sentinel-1) relationship based on recent NDVI and SAR time series and consequently estimates the optimal NDVI value from the current SAR image. The SNAF was tested on 548 commercial fields from 18 countries with 28 crop types and, based on 6880 paired NDVI–SAR images, achieved an RMSE, bias, and R2 of 0.06, 0.00, and 0.92, respectively. The outcome of this study aspires to a persistent seamless stream of NDVI values, regardless of the atmospheric conditions, illumination, or local conditions, which can assist in agricultural decision making.

Biophysical Impact of Sunflower Crop Rotation on Agricultural Fields

Crop rotation is an important determining factor of crop productivity. Sustainable agriculture requires correct rules of crop rotation. Failure to comply with these rules can lead to deterioration of soil biochemical characteristics and land degradation. In Ukraine as well as in many other countries, sunflower monocropping is common practice and the impact of this fact should be studied to find the most precise rules to save the economic potential of land and minimize land degradation factors. This research provides an evaluation of the sunflower monocropping effect on the vegetation indices obtained from MODIS vegetation indices datasets for Ukraine as one of the countries with the biggest sunflower export in Europe. The crop rotation schemes are represented by their area proportions at the village level calculated based on the crop classification maps for 2016 to 2020. This representation gives the possibility to use regression models and f-test feature importance analysis to measure the impact of 3-year and 5-year crop rotation sequences. For these purposes, we use several models: a four-year binary representation model (model A1) and a model with all possible three-year crop rotation scheme representations (model B). These models gave the possibility to evaluate crop rotation schemes based on their biophysical impact on the next sunflower plantings and found that sunflower planting with an interval of three or more years is optimal in terms of the sustainability of soil fertility.

A Flexible Multi-Temporal and Multi-Modal Framework for Sentinel-1 and Sentinel-2 Analysis Ready Data

The rich, complementary data provided by Sentinel-1 and Sentinel-2 satellite constellations host considerable potential to transform Earth observation (EO) applications. However, a substantial amount of effort and infrastructure is still required for the generation of analysis-ready data (ARD) from the low-level products provided by the European Space Agency (ESA). Here, a flexible Python framework able to generate a range of consistent ARD aligned with the ESA-recommended processing pipeline is detailed. Sentinel-1 Synthetic Aperture Radar (SAR) data are radiometrically calibrated, speckle-filtered and terrain-corrected, and Sentinel-2 multi-spectral data resampled in order to harmonise the spatial resolution between the two streams and to allow stacking with multiple scene classification masks. The global coverage and flexibility of the framework allows users to define a specific region of interest (ROI) and time window to create geo-referenced Sentinel-1 and Sentinel-2 images, or a combination of both with closest temporal alignment. The framework can be applied to any location and is user-centric and versatile in generating multi-modal and multi-temporal ARD. Finally, the framework handles automatically the inherent challenges in processing Sentinel data, such as boundary regions with missing values within Sentinel-1 and the filtering of Sentinel-2 scenes based on ROI cloud coverage.

Land Cover Changes after the Massive Rohingya Refugee Influx in Bangladesh: Neo-Classic Unsupervised Approach

The Rohingya refugee influx to Bangladesh in 2017 was a historical incident; the number of refugees was so massive that significant impacts to local communities was inevitable. The Bangladesh government provided land in a preserved area for constructing makeshift camps for the refugees. Previous studies have revealed the land cover changes and impacts of the refugee influx around campsites, especially with regard to local forest resources. Our aim is to establish a convenient approach of providing up-to-date information to monitor holistic local situations. We employed a classic unsupervised technique—a combination of k-means clustering and maximum likelihood estimation—with the latest rich time-series satellite images of Sentinal-1 and Sentinal-2. A combination of VV and normalized difference water index (NDWI) images was successful in identifying built-up/disturbed areas, and a combination of VH and NDWI images was successful in differentiating wetland/saltpan, agriculture /open field, degraded forest/bush, and forest areas. By doing this, we provided annual land cover classification maps for the entire Teknaf peninsula for the pre- and post-influx periods with both fair quality and without prior training data. Our analyses revealed that on-going impacts were still observed by May 2021. As a simple estimation of the intervention consequence, the built-up/disturbed areas increased 6825 ha (compared with the 2015–17 period). However, while the impacts on the original forest were not found to be significant, the degraded forest/bush areas were largely degraded by 4606 ha. These cultivated lands would be used for agricultural activities. This is in line with the reported farmers’ increased income, despite local people with other occupations that are all equally facing the decreases in income. The convenience of our unsupervised classification approach would help keep accumulating a time-series land cover classification, which is important in monitoring impacts on local communities.

A Convolutional Neural Network Algorithm for Soil Moisture Prediction from Sentinel-1 SAR Images

Achieving the rational, optimal, and sustainable use of resources (water and soil) is vital to drink and feed 9.725 billion by 2050. Agriculture is the first source of food production, the biggest consumer of freshwater, and the natural filter of air purification. Hence, smart agriculture is a “ray of hope” in regard to food, water, and environmental security. Satellites and artificial intelligence have the potential to help agriculture flourish. This research is an essential step towards achieving smart agriculture. Prediction of soil moisture is important for determining when to irrigate and how much water to apply, to avoid problems associated with over- and under-watering. This also contributes to an increase in the number of areas being cultivated and, hence, agricultural productivity and air purification. Soil moisture measurement techniques, in situ, are point measurements, tedious, time-consuming, expensive, and labor-intensive. Therefore, we aim to provide a new approach to detect moisture content in soil without actually being in contact with it. In this paper, we propose a convolutional neural network (CNN) architecture that can predict soil moisture content over agricultural areas from Sentinel-1 images. The dual-pol (VV–VH) Sentinel-1 SAR data have being utilized (V = vertical, H = horizontal). The CNN model is composed of six convolutional layers, one max-pooling layer, one flatten layer, and one fully connected layer. The total number of Sentinel-1 images used for running CNN is 17,325 images. The best values of the performance metrics (coefficient of determination (R2=0.8664), mean absolute error (MAE=0.0144), and root mean square error (RMSE=0.0274)) have been achieved due to the use of Sigma naught VH and Sigma naught VV as input data to the CNN architecture (C). Results show that VV polarization is better than VH polarization for soil moisture retrieval, and that Sigma naught, Gamma naught, and Beta naught have the same influence on soil moisture estimation.

Improving the Accuracy of Multiple Algorithms for Crop Classification by Integrating Sentinel-1 Observations with Sentinel-2 Data

The availability of an unprecedented amount of open remote sensing data, such as Sentinel-1 and -2 data within the Copernicus program, has boosted the idea of combining the use of optical and radar data to improve the accuracy of agricultural applications such as crop classification. Sentinel-1’s Synthetic Aperture Radar (SAR) provides co- and cross-polarized backscatter, which offers the opportunity to monitor agricultural crops using radar at high spatial and temporal resolution. In this study, we assessed the potential of integrating Sentinel-1 information (VV and VH backscatter and their ratio VH/VV with Sentinel-2A data (NDVI) to perform crop classification and to define which are the most important input data that provide the most accurate classification results. Further, we examined the temporal dynamics of remote sensing data for cereal, horticultural, and industrial crops, perennials, deciduous trees, and legumes. To select the best SAR input feature, we tried two approaches, one based on classification with only SAR features and one based on integrating SAR with optical data. In total, nine scenarios were tested. Furthermore, we evaluated the performance of 22 nonparametric classifiers on which most of these algorithms had not been tested before with SAR data. The results revealed that the best performing scenario was the one integrating VH and VV with normalized difference vegetation index (NDVI) and cubic support vector machine (SVM) (the kernel function of the classifier is cubic) as the classifier with the highest accuracy among all those tested.

Kc and LAI Estimations Using Optical and SAR Remote Sensing Imagery for Vineyards Plots

Daily or weekly irrigation monitoring conducted per sub-field or management zone is an important factor in vine irrigation decision-making. The objective is to determine the crop coefficient (Kc) and the leaf area index (LAI). Since the 1990s, optic satellite imagery has been utilized for this purpose, yet cloud-cover, as well as the desire to increase the temporal resolution, raise the need to integrate more imagery sources. The Sentinel-1 (a C-band synthetic aperture radar—SAR) can solve both issues, but its accuracy for LAI and Kc mapping needs to be determined. The goals of this study were as follows: (1) to test different methods for integrating SAR and optic sensors for increasing temporal resolution and creating seamless time-series of LAI and Kc estimations; and (2) to evaluate the ability of Sentinel-1 to estimate LAI and Kc in comparison to Sentinel-2 and Landsat-8. LAI values were collected at two vineyards, over three (north plot) and four (south plot) growing seasons. These values were converted to Kc, and both parameters were tested against optic and SAR indices. The results present the two Sentinel-1 indices that achieved the best accuracy in estimating the crop parameters and the best method for fusing the optic and the SAR data. Utilizing these achievements, the accuracy of the Kc and LAI estimations from Sentinel-1 were slightly better than the Sentinel-2′s and the Landsat-8′s accuracy. The integration of all three sensors into one seamless time-series not only increases the temporal resolution but also improves the overall accuracy.

Comparing Sentinel-1 and -2 Data and Indices for Agricultural Land Use Monitoring

Agricultural vegetation development and harvest date monitoring over large areas requires frequent remote sensing observations. In regions with persistent cloud coverage during the vegetation season this is only feasible with active systems, such as SAR, and is limited for optical data. To date, optical remote sensing vegetation indices are more frequently used to monitor agricultural vegetation status because they are easily processed, and the characteristics are widely known. This study evaluated the correlations of three Sentinel-2 optical indices with Sentinel-1 SAR indices over agricultural areas to gain knowledge about their relationship. We compared Sentinel-2 Normalized Difference Vegetation Index, Normalized Difference Water Index, and Plant Senescence Radiation Index with Sentinel-1 SAR VV and VH backscatter, VH/VV ratio, and Sentinel-1 Radar Vegetation Index. The study was conducted on 22 test sites covering approximately 35,000 ha of four different main European agricultural land use types, namely grassland, maize, spring barley, and winter wheat, in Lower Saxony, Germany, in 2018. We investigated the relationship between Sentinel-1 and Sentinel-2 indices for each land use type considering three phenophases (growing, green, senescence). The strength of the correlations of optical and SAR indices differed among land use type and phenophase. There was no generic correlation between optical and SAR indices in our study. However, when the data were split by land use types and phenophases, the correlations increased remarkably. Overall, the highest correlations were found for the Radar Vegetation Index and VH backscatter. Correlations for grassland were lower than for the other land use types. Adding auxiliary data to a multiple linear regression analysis revealed that, in addition to land use type and phenophase information, the lower quartile and median SAR values per field, and a spatial variable, improved the models. Other auxiliary data retrieved from a digital elevation model, Sentinel-1 orbit direction, soil type information, and other SAR values had minor impacts on the model performance. In conclusion, despite the different nature of the signal generation, there were distinct relationships between optical and SAR indices which were independent of environmental variables but could be stratified by land use type and phenophase. These relationships showed similar patterns across different test sites. However, a regional clustering of landscapes would significantly improve the relationships.

New JAXA High-Resolution Land Use/Land Cover Map for Vietnam Aiming for Natural Forest and Plantation Forest Monitoring

Highly detailed and accurate forest maps are important for various applications including forest monitoring, forestry policy, climate change, and biodiversity loss. This study demonstrates a comprehensive and geographically transferable approach to produce a 12 category high-resolution land use/land cover (LULC) map over mainland Vietnam in 2016 by remote sensing data. The map included several natural forest categories (evergreen broadleaf, deciduous (mostly deciduous broadleaf), and coniferous (mostly evergreen coniferous)) and one category representing all popular plantation forests in Vietnam such as acacia (Acacia mangium, Acacia auriculiformis, Acacia hybrid), eucalyptus (Eucalyptus globulus), rubber (Hevea brasiliensis), and others. The approach combined the advantages of various sensor data by integrating their posterior probabilities resulting from applying a probabilistic classifier (comprised of kernel density estimation and Bayesian inference) to each datum individually. By using different synthetic aperture radar (SAR) images (PALSAR-2/ScanSAR, PALSAR-2 mosaic, Sentinel-1), optical images (Sentinel-2, Landsat-8) and topography data (AW3D30), the resultant map achieved 85.6% for the overall accuracy. The major forest classes including evergreen broadleaf forests and plantation forests had a user’s accuracy and producer’s accuracy ranging from 86.0% to 95.3%. Our map identified 9.55 × 106 ha (±0.16 × 106 ha) of natural forests and 3.89 × 106 ha (±0.11 × 106 ha) of plantation forests over mainland Vietnam, which were close to the Vietnamese government’s statistics (with differences of less than 8%). This study’s result provides a reliable input/reference to support forestry policy and land sciences in Vietnam.

Integration of Time Series Sentinel-1 and Sentinel-2 Imagery for Crop Type Mapping over Oasis Agricultural Areas

Timely and accurate crop type mapping is a critical prerequisite for the estimation of water availability and environmental carrying capacity. This research proposed a method to integrate time series Sentinel-1 (S1) and Sentinel-2 (S2) data for crop type mapping over oasis agricultural areas through a case study in Northwest China. Previous studies using synthetic aperture radar (SAR) data alone often yield quite limited accuracy in crop type identification due to speckles. To improve the quality of SAR features, we adopted a statistically homogeneous pixel (SHP) distributed scatterer interferometry (DSI) algorithm, originally proposed in the interferometric SAR (InSAR) community for distributed scatters (DSs) extraction, to identify statistically homogeneous pixel subsets (SHPs). On the basis of this algorithm, the SAR backscatter intensity was de-speckled, and the bias of coherence was mitigated. In addition to backscatter intensity, several InSAR products were extracted for crop type classification, including the interferometric coherence, master versus slave intensity ratio, and amplitude dispersion derived from SAR data. To explore the role of red-edge wavelengths in oasis crop type discrimination, we derived 11 red-edge indices and three red-edge bands from Sentinel-2 images, together with the conventional optical features, to serve as input features for classification. To deal with the high dimension of combined SAR and optical features, an automated feature selection method, i.e., recursive feature increment, was developed to obtain the optimal combination of S1 and S2 features to achieve the highest mapping accuracy. Using a random forest classifier, a distribution map of five major crop types was produced with an overall accuracy of 83.22% and kappa coefficient of 0.77. The contribution of SAR and optical features were investigated. SAR intensity in VH polarization was proved to be most important for crop type identification among all the microwave and optical features employed in this study. Some of the InSAR products, i.e., the amplitude dispersion, master versus slave intensity ratio, and coherence, were found to be beneficial for oasis crop type mapping. It was proved the inclusion of red-edge wavelengths improved the overall accuracy (OA) of crop type mapping by 1.84% compared with only using conventional optical features. In comparison, it was demonstrated that the synergistic use of time series Sentinel-1 and Sentinel-2 data achieved the best performance in the oasis crop type discrimination.