This study aimed to propose an accurate and cost-effective analytical approach for the delineation of fruit trees in orchards, vineyards, and olive groves in Southern France, considering two locations. A classification based on phenology metrics (PM) derived from the Sentinel-2 time series was developed to perform the classification. The PM were computed by fitting a double logistic model on temporal profiles of vegetation indices to delineate orchard and vineyard classes. The generated PM were introduced into a random forest (RF) algorithm for classification. The method was tested on different vegetation indices, with the best results obtained with the leaf area index. To delineate the olive class, the temporal features of the green chlorophyll vegetation index were found to be the most appropriate. Obtained overall accuracies ranged from 89–96% and a Kappa of 0.86–0.95 (2016–2021), respectively. These accuracies are much better than applying the RF algorithm to the LAI time series, which led to a Kappa ranging between 0.3 and 0.52 and demonstrates the interest in using phenological traits rather than the raw time series of the remote sensing data. The method can be well reproduced from one year to another. This is an interesting feature to reduce the burden of collecting ground-truth information. If the method is generic, it needs to be calibrated in given areas as soon as a phenology shift is expected.

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.

Global radiation is a key climate input in process-based models (PBMs) for forests, as it determines photosynthesis, transpiration and the canopy energy balance. While radiation is highly variable at a fine spatial resolution in complex terrain due to shadowing effects, the data required for PBMs that are currently available over large extents are generally at a spatial resolution coarser than ∼9 km. Downscaling large-scale radiation data to the high resolution available from digital elevation models (DEMs) is therefore of potential importance to refine global radiation estimates and improve PBM estimations. In this study, we introduced a new downscaling model that aims to refine sub-daily global radiation data obtained from climate reanalysis data or projections at large scales to the resolution of a given DEM. First, downscaling involves splitting radiation into a direct and diffuse fraction. The influences of surrounding mountains' shade on direct radiation and the “bowl” (deep valley) effect (or sky-view factor) on diffuse radiation are then considered. The model was evaluated by comparing simulated and observed radiation at the Mont Ventoux study site (southeast of France) using the recent ERA5-Land hourly data available at a 9 km resolution as input and downscaled to different spatial resolutions (from 1 km to 30 m resolution) using a DEM. The downscaling algorithm improved the reliability of radiation at the study site, in particular at scales below 150 m. Finally, by using two different PBMs (CASTANEA, a PBM simulating tree growth, and SurEau, a plant hydraulic model simulating hydraulic failure risk), we showed that accounting for fine-resolution radiation can have a great impact on predictions of forest functions.

By combining a crop model (STICS) and a geochemical model (PHREEQC), a new approach to assess the sustainability of agrosystems is proposed. It is based upon aqueous geochemistry and the stepwise modifications of soil solution during its transfer from the surface till aquifer. Meadows of Crau (SE France), irrigated since the 16th century, were field monitored (2012-2015) and modelled. Except for N, the mineral requirements of hay are largely covered by dissolved elements brought by irrigation water with only slight deficits in K and P, which are compensated by P-K fertilizers and the winter pasture by sheep. N cycle results in a very small nitrate leakage. The main determinants of the chemical composition changes of water are: concentration by evaporation, equilibration with soil pCO2, mineral nutrition of plants, input of fertilizers, sheep grazing, mineral-solution interactions in superficial formations till the aquifer, including ion exchange. Inverse modelling with PHREEQC allows for quantifying these processes. For groundwater, measured composition fit statistically very well with those computed, validating thus this approach. This long-term established agrosystem protects both soil and water resources: soil nutritional status remains constant with even some P and (minor) K fixation in soils; long-term decarbonatation occurs but it is greatly slowed by saturation of irrigation water by carbonate; P fixation in soil protects groundwater from eutrophication.

A range of remote sensing platforms provide high spatial and temporal resolution insights which are useful for monitoring vegetation growth. Very few studies have focused on fruit orchards, largely due to the inherent complexity of their structure. Fruit trees are mixed with inter-rows that can be grassed or non-grassed, and there are no standard protocols for ground measurements suitable for the range of crops. The assessment of biophysical variables (BVs) for fruit orchards from optical satellites remains a significant challenge. The objectives of this study are as follows: (1) to address the challenges of extracting and better interpreting biophysical variables from optical data by proposing new ground measurements protocols tailored to various orchards with differing inter-row management practices, (2) to quantify the impact of the inter-row at the Sentinel pixel scale, and (3) to evaluate the potential of Sentinel 2 data on BVs for orchard development monitoring and the detection of key phenological stages, such as the flowering and fruit set stages. Several orchards in two pedo-climatic zones in southeast France were monitored for three years: four apricot and nectarine orchards under different management systems and nine cherry orchards with differing tree densities and inter-row surfaces. We provide the first comparison of three established ground-based methods of assessing BVs in orchards: (1) hemispherical photographs, (2) a ceptometer, and (3) the Viticanopy smartphone app. The major phenological stages, from budburst to fruit growth, were also determined by in situ annotations on the same fields monitored using Viticanopy. In parallel, Sentinel 2 images from the two study sites were processed using a Biophysical Variable Neural Network (BVNET) model to extract the main BVs, including the leaf area index (LAI), fraction of absorbed photosynthetically active radiation (FAPAR), and fraction of green vegetation cover (FCOVER). The temporal dynamics of the normalised FAPAR were analysed, enabling the detection of the fruit set stage. A new aggregative model was applied to data from hemispherical photographs taken under trees and within inter-rows, enabling us to quantify the impact of the inter-row at the Sentinel 2 pixel scale. The resulting value compared to BVs computed from Sentinel 2 gave statistically significant correlations (0.57 for FCOVER and 0.45 for FAPAR, with respective RMSE values of 0.12 and 0.11). Viticanopy appears promising for assessing the PAI (plant area index) and FCOVER for orchards with grassed inter-rows, showing significant correlations with the Sentinel 2 LAI (R² of 0.72, RMSE 0.41) and FCOVER (R² 0.66 and RMSE 0.08). Overall, our results suggest that Sentinel 2 imagery can support orchard monitoring via indicators of development and inter-row management, offering data that are useful to quantify production and enhance resource management.

Accurate data on crop canopy are among the prerequisites for hydrological modelling, environmental assessment, and irrigation management. In this regard, our study concentrated on an in-depth analysis of optical satellite data of Sentinel-2 (S2) time series of the leaf area index (LAI) to characterise canopy development and inter-row management of grapevine fields. Field visits were conducted in the Ouveze-Ventoux area, South Eastern France, for two years (2021 and 2022) to monitor phenology, canopy development, and inter-row management of eleven selected grapevine fields. Regarding the S2-LAI data, the annual dynamic of a typical grapevine canopy leaf area was similar to a double logistic curve. Therefore, an analytic model was adopted to represent the grapevine canopy contribution to the S2-LAI. Part of the parameters of the analytic model were calibrated from the actual grapevine canopy dynamics timing observation from the field visits, while the others were inferred at the field level from the S2-LAI time series. The background signal was generated by directly subtracting the simulated canopy from the S2 LAI time series. Rainfall data were examined to see the possible explanations behind variations in the inter-row grass development. From the background signals, we could group the inter-row management into three classes: grassed, partially grassed, and tilled, which corroborated our findings on the field. To consider the possibility of avoiding field visits, the model was recalibrated on a grapevine field with a clear canopy signal and applied to two fields with different inter-row management. The result showed slight differences among the inter-row signals, which did not prevent the identification of inter-row management, thus indicating that field visits might not be mandatory.