COVID-19 lockdown effects on a coastal marine environment: Disentangling perception versus reality

COVID-19 lockdown brought to a drastic reduction of anthropic impacts on the environment worldwide, including the marine-coastal system. Earth-Observation (EO) data have the potential to monitor and diagnose the effects of the lockdown in terms of water quality. Here we connect the dots among some coastal environmental changes that occurred during the Italian COVID-19 lockdown by using EO data, also seeking to assess connectivity between inland and marine systems. We present a holistic analysis of spatial and temporal variability of environmental parameters in the North Adriatic Sea, Mediterranean basin, exploiting the synergy of different satellite sensors, as well as hydrologic data from in situ observations. Our analysis indicates a favourable interplay of environmental variability that resulted in negative anomalies of Chlorophyll-a concentration, with respect to the climatologic values. Peculiar meteo-oceanographic and hydrological conditions made hard to disentangle potential anthropogenic effects. However, a multi-year hierarchical cluster analysis of riverine remote sensing reflectances groups together the optical properties of inland waters during the lockdown. This emergent cluster highlights the possibility of a second-order, anthropogenic effect that, superimposed to the (first-order) environmental natural causes, may have enhanced water quality during the lockdown.


Introduction
Since 23 February 2020, a series of restrictions, including quarantine and social distancing, were enacted in Italy after the identification of two clusters of COVID-19 in its northern plain (Fig. 1). As of March 9th, widespread lockdown extended to the whole country (Lazzerini and Putoto, 2020;Harlan and Pitrelli, 2020;Di Donato et al., 2020), a fact that represented a unique, unplanned experiment of drastic reduction of anthropic impacts on the environment (e.g., traffic, touristic activities, and some non-primary industrial and agricultural activities; see race.esa.int), in particular, over regions of significant polluted inputs (Kalita and Talukdar, 2021).
In such a context, environmental scientists and the Earth-Observation (EO) community wondered about the real effects of the reduction of anthropogenic input in terms of water quality in coastal areas (Braga et al., 2020;Vijay Prakash et al., 2021). Coastal regions are, indeed, characterized by the presence of many industrialized or touristic areas and they were supposed to be one of the most crucial environments that might experience the effect of the 2020 lockdown. However, these areas are also strongly influenced by multiscale meteo-oceanographic and hydrological processes (Vermaat et al., 2008;Delhez and Barth, 2011;Devlin et al., 2015). Connections between natural environmental variability, anthropogenic forcing, and ecological response in these areas can occur across a large range of interacting spatial, temporal, and organizational scales; this makes monitoring activity a rather challenging task (Stock et al., 2011;Ostendorf, 2011;Dickey-Collas, 2014). Disentangling the mutual effect of anthropic pressures and environmental variability requires a multi-scale approach, i.e., the ability of detecting fast responses on a short time scale and monitoring local (i.e., small-scale) processes over large domains, based on consistent datasets that can assess statistically significant anomalies (Halpern et al., 2008;Zingone et al., 2010;Vermeulen et al., 2011). Satellite-based approaches proved to be the most suitable tools for this task (Harvey et al., 2015;Colella et al., 2016;Gohin et al., 2019).
Phytoplankton bloom occurrence is considered a direct proxy for coastal ecosystem health (Conley et al., 2009;Siegel et al., 2011;Blondeau-Patissier et al., 2014). Anthropogenic enrichment of the water column with nutrients can cause an increase in the biomass of phytoplankton, resulting in a range of undesirable disturbances in the marine ecosystem such as an increased turbidity, a decrease of biodiversity, deoxygenation, fish kills and aquatic vegetation loss (Gulati and van Donk, 2002). Satellite Ocean-Colour (OC) observations provide regular Chlorophyll-a (Chl-a) concentration observations at high spatial and temporal frequencies, using a harmonized methodology over large areas. The two replicates of the Copernicus Multi Spectral Instrument (MSI) optical sensor (Drusch et al., 2012), on board the Sentinel-2 satellites, helped monitoring the evolution of the Venice lagoon ( Fig. 1) water transparency as well as of main European rivers largely affected by goods transportation, tourism, and agriculture (Braga et al., 2020;Zambrano-Monserrate et al., 2020;Yunus et al., 2020). On the other hand, the Ocean and Land Colour Instrument (OLCI) sensor (Donlon et al., 2012), on board the Sentinel-3 satellites, offers a time-series of OC observations at medium spatial resolution (300 m) with higher signal-to-noise ratio and higher temporal frequency, allowing for robust statistical analyses on coastal environments, eventually highlighting signals of human pressures on key hot spots.
By using both MSI and OLCI data, we focus on water quality changes occurred in the North Adriatic Sea (NAS; Mediterranean Sea; Fig. 1) and, in particular, off the Venice lagoon and the Po River Delta coastal area, where intense human activities (e.g., industry, aquaculture, tourism) take place. The NAS (Fig. 1) is a semi-enclosed basin, characterized by a large number of river inputs that act in a limited spatial extent (Cozzi and  Giani , 2011;Brando et al., 2015). Along with the freshwater inputs, winddriven forcing can influence the vertical and horizontal current fields, and, consequently, patterns of biogeochemical tracers in this region (Bignami et al., 2007). The Po River is the largest river of Italy and among the main rivers in the Mediterranean Sea in terms of water discharge, debouching in the NAS (Struglia et al., 2004) (Fig. 1). It flows through one of the most urbanized and productive (mainly agriculture and industry) regions of Europe and its discharge exhibits large inter-annual variability (Zanchettin et al., 2008). Because of widespread of pandemic, since 9 March 2020, the whole region that encompasses the Po River basin and the Venice area experienced a major breakdown of its economy. In particular, the unprecedented water transparency in the Venice lagoon was widely covered by international news outlets and social media (www.nytimes.com/2020/07/02/travel/ venice-coronavirus-tourism.html). Similar conditions were observed in the whole coastal region (Depellegrin et al., 2020) and, in general, in many other touristic and industrial hotspots in the whole country (abcnews.go.com/International/wireStory/italys-seas-speak-touristsboats-cleaner-water-70966642).
Our goal is to assess the environmental effects of this unprecedented change of anthropic activities by providing evidence from EO and in situ data. Emergent patterns are investigated in the light of riverine bio-optical and hydrologic data, thus linking river runoffs, meteooceanographic conditions, water quality, and the potential reduction of anthropic inputs and stressors due to the lockdown.

Remote sensing data and processing
Weekly Chl-a concentrations rely on OLCI full resolution observations. We downloaded the full time series (May 2016-present) of OLCI Level-2 full resolution (300 m) data from the EUMETSAT data centre (eumetsat. int, coda.eumetsat.int and codarep.eumetsat.int). Data were then extracted and remapped over a regular equirectangular grid in the NAS. We applied the CMEMS operational algorithms for phytoplankton Chlorophyll retrieval (Volpe et al., 2019), adapted to OLCI bands (the algorithm exploits the 555 nm band, while OLCI has the band centred over the 560 nm channel) and implemented on a daily basis. Chl-a product is, therefore, the result of a weighted average between an algorithm for open waters (Case I) and another for coastal waters (Case II), in which the weighting factor is calculated for every pixel based on the similarity to the Case I and Case II end members. OLCI daily time series was then turned into a weekly time series by averaging on a pixel-by-pixel basis. OLCI revisit-time at the latitude of the Mediterranean Sea is roughly one day so that a weekly product was deemed of appropriate temporal resolution for allowing the target area to be sampled as completely and regularly as possible. For each pixel the average and standard deviation was computed from a data cube of 3 pixels × 3 pixels × 7 days. This averaging reduced random noise and alleviated the lack of data mostly due to clouds. To investigate the impact of the lockdown, all weekly data from 2016 to 2019 concurred to the generation of the weekly climatology, which was then used as reference for the analysis. Finally, we considered the difference between the 2020 weekly observations and the weekly climatologies. This difference was judged against the climatology in order to obtain a relative, dimensionless value.
MSI Level 1 Top-Of the Atmosphere radiances, with a spatial resolution of 10, 20, and 60 m, were downloaded from Copernicus Open Access Hub and atmospherically corrected with the ACOLITE software (Vanhellemont, 2019) in order to calculate spectra of remote sensing reflectances (R rs ) at Pontelagoscuro gauging station before, during, and after the lockdown. Spectral minimum of total absorption determines the peak at the Green band (560 nm) of the spectra is, usually, influenced by a minimum absorption of algal pigments and scattering by inorganic particles and phytoplankton (Stramski et al., 2004;Gitelson et al., 2008). The secondary peak at 704 nm is caused by a small window absorption between the second absorption peak of Chl-a at 665 nm (Gitelson et al., 2008) and the sharply increasing absorption of water molecules. This peak only becomes visible at significant Chl-a concentrations, mostly findable in coastal or inland waters (Gitelson, 1992, Gons et al., 2002, Dall'Olmo and Gitelson, 2005. To detect the characteristic peaks in the spectral reflectance associated with significant Chl-a concentrations in turbid, productive waters, we calculated the Normalized Difference Chlorophyll Index (NDCI) (Mishra and Mishra, 2012) and adapted to MSI bands (Caballero et al., 2020): Remote-sensing reflectances (R rs ) were also converted into Water Turbidity, expressed in formazin nephelometric unit [FNU] (Dogliotti et al., 2015) and into Chl-a concentration (in mg m −3 ) (Pinardi et al., 2018). To assess the accuracy of turbidity products derived from MSI imagery, in situ data collected in the Venice lagoon was used. In particular, turbidity profiles were recorded during field activities from 2017 onwards with a SeaPoint turbidity meter (SeaPoint®). The correlation for 198 match-ups between satellite-derived estimates of turbidity versus in situ measurements was statistically significant with a coefficient of correlation r of 0.962, a coefficient of determination R 2 of 0.925, a RMSE of 2.44 FNU ( Supplementary Fig. S1). These results are consistent with the accuracy assessed for ACOLITE-derived turbidity retrieved from Landsat-8 in the North Adriatic coastal waters (Braga et al., 2017). R rs , Chl-a concentration, and Water Turbidity were extracted from a 5 × 5 pixels region of interest in the Po River, near Pontelagoscuro gauging station.
Satellite-based Sea-Surface Temperature (SST) was obtained from CMEMS (Product SST_MED_SST_L4_NRT_OBSERVATIONS_010_004), by a gap-free, optimally interpolated product and was built by using observations acquired from different platforms, e.g. Envisat, NOAA-18, MetOpA, Aqua, Terra, MetOpB, NPP, Meteosat11, as well as Sentinel 3A and Sentinel 3B satellites (Buongiorno Nardelli et al., 2013;Buongiorno Nardelli et al., 2015). The optimal interpolation procedure consists in the data extraction and quality control, the merging of the multi-platform observations and the interpolation on a regular 0.01°g rid. The gap-free analysis processing is obtained via objective analysis (Le Traon, 1990) collecting SST observations free of diurnal warming effects (from 18:00 to 6:00 UTC) using spatial and temporal decorrelation lengths of 5 km and 1 day, respectively. This SST is thus representative of a foundation SST (i.e., temperature free of diurnal temperature variability) (Minnett et al., 2019).

In situ data
Data of hourly water discharge for the Po River at the Pontelagoscuro gauging station were obtained from the Environmental Protection Agency of the Region of Emilia Romagna (ARPA-ER). Sea temperature and wind direction and speed were measured at AAOT and are provided by the Centro Previsione e Segnalazione Maree of the city of Venice.

Marine currents
The Mediterranean Forecasting System is a hydrodynamic primitive equations model for the Mediterranean Basin (Clementi et al., 2019). The model outputs are distributed through the Copernicus Marine Environment Monitoring Service (CMEMS, Product MEDSEA-ANALYSIS-FORECAST-PHY-006-013) and provide daily to hourly fields of horizontal currents, temperature, salinity and free-surface elevation. The outputs are given on a 1/24°regular grid and 125 unequally spaced vertical levels. The core of the hydrodynamic model is the NEMO (Nucleus for European Modelling of the Ocean) model, the wave component is provided by Wave Watch-III. The numerical simulations assimilate satellite sea-surface temperature, altimeter-derived along-track sea-level anomalies and in situ measured salinity profiles.

Hierarchical cluster analysis
To classify the entire R rs spectra acquired by Sentinel-2 at 9 bands in the range 443-865 nm into cluster trees (i.e., dendrograms), a hierarchical cluster analysis (HCA) was applied for the 2020 R rs spectra (time constrained single-year) and for R rs spectra collected from 2018 to 2020 (multi-year). The use of HCA has been widely stressed for application of remote sensing, using hyperspectral data, especially when coupled with the derivative spectra of R rs (and its derived optical properties) that enhance the signature of absorbing material in water such as phytoplankton Chlorophyll and other pigments. The classification of spectra was performed by using linkage algorithms and similarity criteria (Torrecilla et al., 2011;Organelli et al., 2017). The high similarity we observed in this study was due to the analysis performed on ordinary spectra composed by only 9 bands. Briefly, the unweighted pair-group average linkage algorithm was used to group the spectra according to the average distance between all members. To better characterize the impact of lockdown phases on the Po River optical properties, the single-year HCA was constrained in time in order to compute the similarity between spectra acquired consecutively during 2020. The multi-year HCA was not time constrained, in order to highlight differences of the 2020 R rs spectra with those collected during the other years. The similarity level between spectra was based on the cosine distance which returns values from 0 (i.e., no similarity) to 1 (i.e. the high similarity). Cluster analysis was carried out by the free statistical software PAST version 3.04.
Along with the single-year and multi-year HCA, we performed the analysis of variance (ANOVA) of Chl-a concentration, Water Turbidity, and Po River discharge, calculating the probability values for each single parameter. In particular, for each cluster of the two HCAs we calculated mean values, standard deviation, and probability values (p-value) of Chl-a concentration, Water Turbidity, and Po River discharge.

Anomalies of OLCI-based Chlorophyll-a concentration
We calculated weekly deviations from climatology of Chl-a concentrations by using OLCI data over the NAS. In our study area, strongly influenced by river inputs and human activities, large values of Chl-a concentration result from the discharge of urban sewage, industrial runoffs, and fertilizers from agriculture activities over watersheds (Penna et al., 2004;Colella et al., 2016;Cozzi and Giani, 2011).
Starting from 18 March 2020 (i.e., about a week after the beginning of the Italian 2020 lockdown), the sequence of weekly mean and anomalies highlighted a general, decreasing trend of Chl-a concentration along the Italian coastal zone facing the NAS as well as off the whole Po River Delta (Fig. 2). In particular, around mid-April 2020, the weekly anomalies showed very low values (Fig. 2b, c, d), with negative, relative changes close to 100%.
Although such an intriguing pattern might be due to (or at least enhanced by) natural environmental conditions, such as moderate winds and riverine drought, it is worth investigating the striking coincidence between the consistent Chl-a concentration decrease, occurred form mid-March and the 2020 lockdown, in the light of variability of other natural phenomena. Indeed, while such a decrease of biomass was expected to occur in coastal areas (Braga et al., 2020;Zambrano-Monserrate et al., 2020;Yunus et al., 2020), offshore areas might not be directly affected by inputs of anthropogenic origin.

Meteo-oceanographic conditions
The negative anomalies of Chl-a concentration (Fig. 2) we found off the Venice Lagoon and the Po River Delta likely result from the cumulated  impact of a number of natural phenomena. Low wind magnitudes, recorded at the Acqua Alta Oceanographic Tower (AAOT) station (Fig. 1), from 5 to 20 April 2020, suggested the maintenance of vertical stratification of the water column (Fig. 3a, b), which may have hindered nutrient-rich waters due to mixing processes that would have enhanced the growth of phytoplankton. This hypothesis is confirmed by a significant increase of water temperature at the AAOT station during that period (Fig. 3a, b).
Satellite-derived Sea Surface Temperature (SST) at AAOT (Fig. 3b, validated with in situ measurements at AAOT) shows an increase, in particular, from 5 to 20 April 2020, likely due to the mild wind activity (Fig. 3b, c, d,  e). This condition resulted in a weak alongshore transport of riverine nutrients from the northern sector of the basin (Fig. 3f, g, h). The warm riverplume bulge and the associated SST pattern along the north-west coast of the basin marked a reduction of the well-known cyclonic, cold coastal  current (Bignami et al., 2007;Nof and Pichevin, 2001) that constitutes the natural conveyor belt of the North Adriatic river-plumes, i.e., the primary source of nutrient and sediment loads that flow southward. This likely resulted to a mild riverine nutrient advection along the coastal area, which was further favoured by the absence of wind-induced mixing we recorded during that period (Fig. 3).

River inputs and MSI-derived Water Turbidity
Hydrological data of the Po River at Pontelagoscuro station (~90 km upstream the river mouth) showed that the lowest anomaly of Chl-a concentration off the Po River Delta (Fig. 2) occurred during a minimum stage of water discharge, recorded between 10 and 20 April 2020 (Fig. 4a). The low water discharge we observed (about 40% lower than the 2010-2019 climatological value) confirmed the mild advection of nutrients from the Po River tributaries, as well as the low input of cold riverine water we observed from SST imagery off the Po River delta (Fig. 3c, d, e).
Water Turbidity values, derived from MSI images near Pontelagoscuro gauging station, confirmed a general correlation between nutrient-poor river runoff and synchronous water discharge measurements (Fig. 4b, c, d, e). This correlation seems to be slightly affected by seasonal effects and results to be particularly evident after the lockdown. Before and during the lockdown, indeed, water discharge and turbidity show a weak correlation, likely due to threshold limitations on sediment supply and/or biogeochemical activity (Syvitski and Kettner, 2011;Tesi et al., 2013).

MSI-based spectral classification of the Po River runoff
To investigate the inland water bio-optical properties of the Po River and their eventual relationship with the biogeochemical pattern we (1) (see Methods). Water turbidity was retrieved according to Dogliotti et al. (2015). Colours indicate the pre-lockdown (yellow), the lockdown (blue), and the post-lockdown (red) periods. observed along the coastal areas, we analyzed MSI remote sensing Reflectance (R rs ) spectra, collected before, during, and after the 2020 lockdown, extracted at Pontelagoscuro gauging station ( Supplementary Fig. S1). MSI spectra showed that the R rs peak magnitude in the NIR was often comparable to the green band (560 nm), indicating high Chl-a, while, for some dates, the NIR reflectance values were smaller than the R rs in the red band (665 nm), due to low Chl-a concentration. In particular, four R rs spectra in the pre-lockdown period (January-February 2020; Supplementary Fig. S1) showed similar characteristics in magnitude and shape as those of lockdown period (March-May 2020; Supplementary Fig. S1), except the two spectra with a reflectance peak in the NIR, due to high Chl-a concentration, occurred on 8 and 23 April 2020.
Spectra for the post-lockdown ( Supplementary Fig. S1) exhibited highly variation in both magnitude and shape. A significant R rs peak at 704 nm, respect to neighbouring bands is characteristic of high Chla concentrations, mostly findable in coastal or inland waters (Gitelson, 1992, Gons et al., 2002, Dall'Olmo and Gitelson, 2005. Here, this feature was quantified with the Normalized Difference Chlorophyll Index (NDCI; Fig. 4f). Positive NDCI values were observed at Po River discharges that were lower than 1000 m 3 s −1 . Negative NDCI values were observed for higher discharges. This particular pattern clearly suggests the presence of a Chlorophyll-a absorption occurring during low discharge, while higher energetic river conditions seem to hinder phytoplankton growth (Pinder et al., 1997).
Among the five groups identified by HCA, the first cluster included only spectra corresponding to the pre-lockdown period. The spectra we observed during the lockdown were instead grouped into two different clusters, and the second one showed high similarity with spectra from 23 May to 27 July 2020, which corresponded to the post-lockdown phase. The remaining spectra, acquired during post-lockdown in August, September, and October 2020, were contained in the fourth and fifth clusters ( Fig. 5a; Table 1). Thus, the stop of the human activities due to the COVID-19 pandemic was not sufficient to significantly affect the optical properties of the Po River. This suggests that the spectral variability of R rs was driven, at first order, by natural phenomena.
The five clusters of 2020 MSI R rs spectra identified by HCA resulted from the seasonal variability of Chl-a concentration, water turbidity and discharge, observed for the Po River ( Fig. 5b; Table 1). The spectra from the first, and pre-lockdown cluster were acquired when the river discharge was decreasing and the Water Turbidity was still high (>20 FNU). The interplay between low Chl-a concentrations in the river, low turbidity (<30 FNU), and low discharge (about 1000 m 3 s −1 ) resulted in the second cluster of spectra collected during the lockdown between 19 March and 8 April 2020. Differently, the spectral features of R rs in the third cluster responded to a new increase of the Po River discharge (up to 2500 m 3 s −1 ) and thus of the water turbidity due to snow melting from the Alps and the Apennines (Zanchettin et al., 2008;Grilli et al., 2020). In August 2020, when the impact of water coming from surrounding mountains was low (i.e., riverine discharges < 1000 m 3 s −1 ), high Chl-a concentrations (up to 60 mg m −3 ) determined the R rs spectral optical properties grouped within the fourth cluster ( Fig. 5a, b; Table 1). A new increase in water turbidity and discharge differentiated further the R rs spectral properties from September and Fig. 5. a) Results of the HCA performed on MSI R rs spectra between 443 and 865 nm, extracted from the Po river channel near Pontelagoscuro gauging station for the periods before, during and after the lockdown. b) Chlorophyll concentration, Water Turbidity (derived from MSI), and Po River discharge near the Pontelagoscuro gauging station corresponding to the date of acquisition of the Sentinel2 R rs spectra. Data are shown in chronological order. c) Results of the multi-year HCA performed on MSI R rs spectra between 443 and 865 nm, extracted from the Po River channel at Pontelagoscuro gauging station from 2018 to 2020. d) Means (and standard deviation) of Chl-a concentration, Water Turbidity (derived from MSI), and Po River discharge near the Pontelagoscuro gauging station corresponding to each cluster, identified by the HCA.
October 2020 (i.e. the fifth cluster) as a consequence of the heavy precipitations that characterize the northern part of Italy (Grilli et al., 2020).
From this analysis, we thus observed that the optical characteristics of the Po River in 2020 followed the natural variability of such riverine system, with no significant evidence of the lockdown (Fig. 5a, b). The resulting high-level riverine water quality was primary due to seasonal and hydrological conditions. However, we also notice from the analysis of variance (ANOVA) that, among the 2020 clusters, the Chl-a concentration is the sole parameter that shows statistically significant differences in its variability (Table 1). We therefore hypothesize that changes of Chl-a concentration might not strictly connected to the sole Po River discharge and Water Turbidity changes. This suggests the hypothesis of a second-order, nonseasonal, effect that slightly enhanced the good water quality we observed during the 2020 lockdown.
Accordingly, a multi-year HCA and the related ANOVA, conducted on MSI R rs spectra from 2018 to 2020 at Pontelagoscuro gauging station, seems to highlight a cluster, i.e., Cluster 6 in Fig. 5c, d, that reinforce the hypothesis of an anomalous behaviour of Chl-a concentration during spring 2020 ( Fig. 5c; Table 1). Indeed, Cluster 6 is identified by only R rs spectra 2020, among which three of them are from spring 2020 spectra. Whereas, pre-and post-lockdown spectra from the year 2020 were well grouped with those of the other years and likely mark the seasonal runoff of the Po River. Chl-a concentration of Cluster 6 is particularly low and the ANOVA confirms that difference in Chl-a concentrations among clusters is statistically significant (Table 1). In particular, the probability value of mean Chl-a concentration of Cluster 6, when compared with Cluster 1-to-5 (all grouped together) is~0.06, i.e., a value that mark the anomaly of this 2020 cluster in terms of Chl-concentration, with respect to the other clusters.

Conclusions
During the lockdown, significant negative anomalies of weekly Chl-a concentration were observed over a large part of the coastal, western sector North Adriatic basin. The particular combination of meteo-oceanographic and hydrological conditions was able to explain the negative anomalies of Chl-a concentration offshore the Venice lagoon and off the Po River delta. The mutual effect of environmental parameters, such as near surface winds and river discharge that caused well-stratified thermohaline conditions and a negligible runoff of nutrient inputs, maintained values of Chla concentration fairy low during the lockdown. Moreover, a Hierarchical Cluster Analysis, along with an analysis of variance, suggested that the variability of Chl-a concentration we observe from winter to fall 2020 is in accordance with a seasonal cycle.
Said that, we could not exclude a second order effect due to the decrease of human activities. It was, indeed, reasonable to hypothesize that the favourable environmental conditions we assessed did not hinder the effect of a reduction of anthropogenic inputs. The multi-year Hierarchical Cluster Analysis showed that the most of the spring 2020 R rs spectra (i.e., those that correspond to the 2020 lockdown) tended to group in a specific cluster and, compared with other clusters that range from 2018 to 2020, this cluster showed a statistically significant difference in Chl-a concentration. In other words, such an intriguing, low value of Chl-a concentration we observed in both inland and marine waters during the 2020 lockdown, although mostly explained by seasonal variability of meteo-oceanographic conditions, statistically suggested a secondary, likely anthropogenic effect due to the lockdown restrictions, when analyzed in the frame of multiyear remote sensing and in situ data.
We finally remark that coastal dynamics of river plumes that show a significant reduction of nutrient loads of anthropogenic origin may be at the base of connectivity between inland and offshore water biogeochemical conditions. To test this hypothesis, it will be crucial the ad hoc synergy of additional satellite imagery for inland waters, including hyperspectral products, and in situ measurements of chemical/pollutant components. CRediT authorship contribution statement FF and MHR, delineated and supervised the study. FB, DC, SC, EO, JP, VEB, MB, JC, CG, GMS, and GV collected, processed, and analyzed the Table 1 Mean values and standard deviations (σ) of Chl-a concentration, Water Turbidity, and Po River water discharge related to the cluster we identified in both single-year (2020) and multi-year (2018-2020) HCA. The table shows the probability value (p-value) obtained from the ANOVA; Geen (yellow) shading indicates that the ANOVA p-value is less (greater) than the significant level of 0.05. Both p-values of the single-and multiyear analysis show that Chl-a concentration changes are statistically significant.

Declaration of competing interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.