In the global change context, an efficient management of the available resources has become one of the most important topics particularly for sustainable crop development. Many questions concern the evolution of the rice farming systems in Camargue in Southeastern France, which play a crucial role in controlling the soil salinity. Their surface area significantly decreased from 20, 000 ha in 2010 to 14000 ha in 2014. The arrival of the new Sentinel satellites makes it possible to evaluate these crop evolutions. The objectives of this study were to propose operational methodologies to: (1) accurately assess the surface areas of the main crops: rice, wheat and grassland from classifications based on multispectral data, (2) map agricultural practices (sowing and harvest residue burning), and (3) elaborate a farm typology based on variables computed from remote sensing data to better understand the farming strategies. Dense time series of Sentinel images acquired at high spatial resolution (10m) were analyzed for 2016 and 2017. A satisfactory accuracy was obtained for land use classification with 88% of correctly classified fields. The accuracy obtained for the estimation of the sowing date varied according to the studied year from 8 to 12 days, and burned areas were correctly identified (80%). The farm typology allowed to cluster farms at territory level.

The assessment of rice yield at territory level is important for strategic economic decisions. Assessing spatial and temporal yield variability at regional scale is difficult because of the numerous factors involved, including agricultural practices, phenological calendars, and environmental contexts. New remote sensing data acquired at decametric resolution (Sentinel missions) can provide information on this spatial variability. The study objective was thus to evaluate the potential of Sentinel-2 images for monitoring rice cropping systems and yield from farm to region scales. The approach considered both observations and modeling. Indepth farmers surveys were carried out in the Camargue region, Southeastern France. The novelty was to use operational tools (BVNET and PHENOTB) to compute leaf area index, to daily interpolate this biophysical variable from 44 images acquired in 2016 and 2017 for each rice field, and to derive key phenological parameters from the analysis of the temporal profiles. The STICS crop model was spatially used, considering the biophysical variables derived from remote sensing. We tested four simulation strategies, differing in the integration intensity of remote sensing information into the model. Results have shown that (1) Sentinel-2 data allowed distinguishing early and late rice varieties. (2) The phenological stages mapped at the regional level allowed to better understand the agricultural practices of farmers. (3) The assimilation of remote sensing data to the STICS crop model significantly improved yield estimation and provided useful information on the spatial variability observed at regional scale. It was the first time that Sentinel-2 data are used with STICS crop model to assess rice yield at both farm and regional scale in the Camargue area. The proposed method is based on free open data and free access model, easily reproducible in other environmental contexts.

Conventional methods of crop mapping need ground truth information to train the classifier. Thanks to the frequent acquisition allowed by recent satellite missions (Sentinel 2), we can identify temporal patterns that depend on both phenology and crop management. Some of these patterns are specific to a given crop and thus can be used to map it. Thus, we can substitute ground truth information used in conventional methods with agronomic knowledge. This approach was applied to identify irrigated permanent grasslands (IPG) in the Crau area (Southern France), which play a crucial role in groundwater recharge. The grassland is managed by making three mows during the May–October period, which leads to a specific temporal pattern of leaf area index (LAI). The mowing detection algorithm was designed using the temporal LAI signal derived from Sentinel 2 observations. The algorithm includes some filtering to remove noise in the signal that might lead to false mowing detection. A pixel is considered a grassland if the number of detected mows is greater than 1. A data set covering five years (2016–2020) was used. The detection mowing number was conducted at the pixel level, and then the results were aggregated at the plot level. An evaluation data set including 780 plots was used to assess the performances of the classification. We obtained a Kappa index ranging between 0.94 and 0.99 according to the year. These results were better than other supervised classification methods that include training data sets. The analysis of land-use changes shows that misclassified plots concern grasslands managed less intensively with strong intra-parcel heterogeneity due to irrigation defects or year-round grazing. Time series analysis, therefore, allows us to understand different management practices. Real land-use change in use can be observed, but long time series are needed to confirm the change and remove ambiguities with heterogeneous grasslands.