Recent literature reflects the substantial progress in combining spatial, temporal and spectral capacities for remote sensing applications. As a result, new issues are arising, such as the need for methodologies that can process simultaneously the different dimensions of satellite information. This paper presents PLS regression extended to three-way data in order to integrate multiwavelengths as variables measured at several dates (time-series) and locations with Sentinel-2 at a regional scale. Considering that the multi-collinearity problem is present in remote sensing time-series to estimate one response variable and that the dataset is multidimensional, a multiway partial least squares (N-PLS) regression approach may be relevant to relate image information to ground variables of interest. N-PLS is an extension of the ordinary PLS regression algorithm where the bilinear model of predictors is replaced by a multilinear model. This paper presents a case study within the context of agriculture, conducted on a time-series of Sentinel-2 images covering regional scale scenes of southern France impacted by the heat wave episode that occurred on 28 June 2019. The model has been developed based on available heat wave impact data for 107 vineyard blocks in the Languedoc-Roussillon region and multispectral time-series predictor data for the period May to August 2019. The results validated the effectiveness of the proposed N-PLS method in estimating yield loss from spectral and temporal attributes. The performance of the model was evaluated by the R2 obtained on the prediction set (0.661), and the root mean square of error (RMSE), which was 10.7%. Limitations of the approach when dealing with time-series of large-scale images which represent a source of challenges are discussed; however, the N–PLS regression seems to be a suitable choice for analysing complex multispectral imagery data with different spectral domains and with a clear temporal evolution, such as an extreme weather event.

Unexpected climatic conditions or extreme climatic events in vineyards are a worldwide problem that requires accurate spatial and temporal monitoring. Satellite-based remote sensing is an important source of data to assess this challenge in a climate-change context. This paper provides a first insight into the capacity of a multiway analysis method applied to Sentinel-2 time series to assess heatwave impacts in vineyards at a regional scale. Multi-way partial least squares (N-PLS) regression was used as a supervised technique to predict the intensity of damage caused to vineyards by the heatwave phenomenon that impacted the vineyards in the south of France in 2019. The model was developed based on available ground truth data of yield losses for 107 vineyard blocks in the Languedoc-Roussillon region and multispectral time-series predictor data for the period May to August 2019. The model showed a performance accuracy (R2) of 0.56 in the calibration set and of 0.66 in the validation set, with a standard error of cross-validation in the calibration set of 12.4% and a standard error of the prediction of yield losses in the validation set of 10.7. The model was applied at a regional scale on 4978 vineyard blocks to predict yield losses using spectral and temporal attributes. The prediction of the yield loss due to heat stress at a regional scale was related to the spatial pattern of maximum temperatures recorded during the extreme weather event. This relation was confirmed by a chi-square test (p < 5%). The introduction of N-PLS insights into the analysis enables the characterisation of heat stress responses in vineyards and the identification of spectro-temporal profiles relevant for understanding the effects of heatwaves on vine blocks at a regional scale.

Monitoring wine-growing regions and maximising the value of production based on their region/local specificities requires accurate spatial and temporal monitoring. The increasing amount and variability of information from remote sensing data is a potential tool to assess this challenge for the grape and wine industry. This article provides a first insight into the capacity of a multiway analysis method applied to Sentinel-2 time series to assess the value of simultaneously considering spectral and temporal information to highlight site-specific canopy evolution in relation to environmental factors and management practices, which present a large diversity at this regional scale. Parallel Factor Analysis (PARAFAC) was used as an unsupervised technique to recover pure spectra and temporal signatures from multi-way spectral imagery of vineyards in the Languedoc-Roussillon region in the south of France. The model was developed using a time series of Sentinel-2 satellite imagery collected over 4978 vineyard blocks between May 2019 and August 2020. From the Sentinel-2 (spectral and temporal) signal, the PARAFAC analysis allowed the identification of spectral and temporal profiles in the form of pure components, which corresponded to vegetation and soil. The PARAFAC analysis also identified that two of the pure spectra were strongly related to characteristics and dynamics of vineyard cultivation at a regional scale. A conceptual framework was proposed in order to simultaneously consider both vegetation and soil profiles and to summarise the mass of data accordingly. This methodology allowed the computation of a concentration index that characterised how close a field was to a vegetation or a soil profile over the season. The concentration indices were validated for the vegetation and the soil over two growing seasons (2019 and 2020) with geostatistical analysis. A non-random distribution of the concentration index at the regional scale was assumed to highlight a strongly spatially organised phenomenon related to spatially organised environmental factors (soil, climate, training system, etc.). In a second step, spatial patterns of indices were subjected to the expertise of a panel of advisors of the wine industry in order to validate them in relation to vine-growing conditions. Results showed that the introduction of the PARAFAC method opened up the possibility to identify relevant spectro-temporal profiles for vine monitoring purposes.

Multispectral image time-series have been promising for some years; yet, the substantial advance of the technology involved, with unprecedented combinations of spatial, temporal, and spectral capabilities for remote sensing applications, raises new challenges, in particular, the need for methodologies that can process the different dimensions of satellite information. Considering that the multi-collinearity problem is present in remote sensing time-series, regression models are widespread tools to model multi-way data. This paper presents the results of the analysis of a high order data of Sentinel-2-time series, conducted in the framework of extreme weather event. A feature extraction method for multi-way data, N-CovSel was used to identify the most relevant features explaining the loss of yield in Mediterranean vineyards during the 2019 heatwave. Different regression models (uni-way and multi-way) from features extracted from the N-CovSel algorithm were calibrated based on available heat wave impact data for 107 vineyard blocks in the Languedoc-Roussillon region and multispectral time-series predictor data for the period May to August. The performance of the models was evaluated by the r2 and the root mean square of error (RMSE) as follows: for the temporal N-PLS model (r2 = 0.62—RMSE = 11%), for the spatial N-PLS model (r2 = 0.61—RMSE = 12%) and the temporal-spectral PLS model (r2 = 0.63—RMSE = 11%). The results validated the effectiveness of the proposed N-CovSel algorithm in order to reduce the number of total variables and restricting it to the most significant ones. The N-CovSel algorithm seems to be a suitable choice to interpret complex multispectral imagery by temporally discriminating the most appropriate spectral information.