The objective of this study is to address the capabilities of multi-temporal optical images to estimate the fine-scale yield variability of wheat, over a study site located in southwestern France. The methodology is based on the Landsat-8 and Sentinel-2 satellite images acquired after the sowing and before the harvest of the crop throughout four successive agricultural seasons, the reflectance constituting the input variables of a statistical algorithm (random forest). The best performances are obtained when the Normalized Difference Vegetation Index (NDVI) is combined with the yield maps collected during the crop rotation, the agricultural season 2014 showing the lower level of performances with a coefficient of determination (R²) of 0.44 and a root mean square error (RMSE) of 8.13 quintals by hectare (q.h⁻¹) (corresponding to a relative error of 12.9%), the three other years being associated with values of R² close or upper to 0.60 and RMSE lower than 7 q.h⁻¹ (corresponding to a relative error inferior to 11.3%). Moreover, the proposed approach allows estimating the crop yield throughout the agricultural season, by using the successive images acquired from the sowing to the harvest. In such cases, early and accurate yield estimates are obtained three months before the end of the crop cycle. At this phenological stage, only a slight decrease in performance is observed compared to the statistic obtained just before the harvest.

Sentinel-2 (S2) earth observation satellite mission, launched in 2015, is foreseen to promote within-field decisions in Precision Agriculture (PA) for both: (1) optimizing crop production; and (2) regulating environmental impacts. In this second scope, a set of Leaf Area Index (LAI) derived from S2 type time-series (2006–2010, using Formosat-2 satellite) is used to spatially constrain the within-field crop growth and the related nitrogen contamination of surface water simulated at a small experimental catchment scale with the distributed agro-hydrological model Topography Nitrogen Transfer and Transformation (TNT2). The Soil Water Holding Capacity (SWHC), represented by two parameters, soil depth and retention porosity, is used to fit the yearly maximum of LAI (LAX) at each pixel of the satellite image. Possible combinations of soil parameters, defining 154 realistic SWHC found on the study site are used to force spatially homogeneous SWHC. LAX simulated at the pixel level for the 154 SWHC, for each of the five years of the study period, are recorded and hereafter referred to as synthetic LAX. Optimal SWHC_{year_I,pixel_j}, corresponding to minimal difference between observed and synthetic LAX_{year_I,pixel_j}, is selected for each pixel, independent of the value at neighboring pixels. Each re-estimated soil maps are used to re-simulate LAX_{year_I}. Results show that simulated and synthetic LAX_{year_I,allpixels} obtained from SWHC_{year_I,allpixels} are close and accurately fit the observed LAX_{year_I,allpixels} (RMSE = 0.05 m²/m² to 0.2 and R² = 0.99 to 0.94), except for the year 2008 (RMSE = 0.8 m²/m² and R² = 0.8). These results show that optimal SWHC can be derived from remote sensing series for one year. Unique SWHC solutions for each pixel that limit the LAX error for the five years to less than 0.2 m²/m² are found for only 10% of the pixels. Selection of unique soil parameters using multi-year LAX and neighborhood solution is expected to deliver more robust soil parameters solutions and need to be assessed further. The use of optical remote sensing series is then a promising calibration step to represent crop growth within crop field at catchment level. Nevertheless, this study discusses the model and data improvements that are needed to get realistic spatial representation of agro-hydrological processes simulated within catchments.

Unmanned aerial vehicles (UAVs) provide images at decametric spatial resolutions. Their flexibility, efficiency, and low cost make it possible to apply UAV remote sensing to multisensor data acquisition. In this frame, the present study aims at employing RGB UAV images (at a 3 cm resolution) and multispectral images (at a 16 cm resolution) with related vegetation indices (VIs) for mapping surfaces according to their illumination. The aim is to map land cover in order to access temperature distribution and compare NDVI and MTVI2 dynamics as a function of their illuminance. The method, which is based on a linear discriminant analysis, is validated at different periods during the phenological cycle of the crops in place. A model based on a given date is evaluated, as well as the use of a generic model. The method provides a good capacity of separation between four classes: vegetation, no-vegetation, shade, and sun (average kappa of 0.93). The effects of agricultural practices on two adjacent plots of maize respectively submitted to conventional and conservation farming are assessed. The transition from shade to sun increases the brightness temperature by 2.4 °C and reduces the NDVI by 26% for non-vegetated surfaces. The conservation farming plot is found to be 1.9 °C warmer on the 11th of July 2019, with no significant difference between vegetation in the sun or shade. The results also indicate that the NDVI of non-vegetated areas is increased by the presence of crop residues on the conservation agriculture plot and by the effect of shade on the conventional plot which is different for MTVI2.

Recent advances in sensors onboard harvesting machines allow accessing the intra-plot variability of yields, spatial scale fully compatible with numerous on-going satellite missions. The aim of this study is to estimate the sunflower yield at the intra-plot spatial scale using the multi-temporal images provided by the Landsat-8 and Sentinel-2 missions. The proposed approach is based on a statistical algorithm, testing different sampling strategies to partition the dataset into independent training and testing sets: A random selection (testing different ratio), a systematic selection (focusing on different plots) and a forecast procedure (using an increasing number of images). Emphasis is put on the use of high spatial and temporal resolution satellite data acquired throughout two agricultural seasons, on a study site located in southwestern France. Ground measurements consist in intra-plot yields collected by a surveying harvesting machine with GPS system on track mode. The forecast of yield throughout the agricultural season provides early accurate estimation two months before the harvest, with R² equal to 0.59 or 0.66 and root mean square error (RMSE) of 4.7 or 3.4 q ha⁻¹, for the agricultural seasons 2016 and 2017 respectively. Results obtained with the random selection or the systematic selection will be developed later, in a longer paper.