Next Article in Journal
Effective Removal of Pb2+ from Aqueous Solution Using Magnetic Mesoporous Silica Prepared by Rubidium-Containing Biotite Leaching Residues and Wastewater
Next Article in Special Issue
Spatial-Temporal Evolution Characteristics and Driving Force Analysis of NDVI in the Minjiang River Basin, China, from 2001 to 2020
Previous Article in Journal
Evaluating Monthly Flow Prediction Based on SWAT and Support Vector Regression Coupled with Discrete Wavelet Transform
Previous Article in Special Issue
Realizing Sustainable Development of Yellow River Basin by Horizontal Eco-Compensation Based on Integrated Water Rights (IWRs) Transactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evidence and Implications of Hydrological and Climatic Change in the Reno and Lamone River Basins and Related Coastal Areas (Emilia-Romagna, Northern Italy) over the Last Century

Department of Biological, Geological and Environmental Sciences, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Water 2022, 14(17), 2650; https://doi.org/10.3390/w14172650
Submission received: 21 July 2022 / Revised: 24 August 2022 / Accepted: 25 August 2022 / Published: 28 August 2022

Abstract

:
Climate change and human activities have consequences on coastal areas as they affect hydrological processes in the related river basins. The riverine sediment supply to the beaches of the Emilia-Romagna coast, a highly urbanized area with high economic and naturalistic value, has been heavily impacted by human activities throughout the catchment, reducing solid transport to the coast and increasing the threat of coastal erosion and flooding. Despite the introduction of safeguard policies in the early 1980s and the consequent stoppage of such activities, the expected return in solid transport has not yet been reflected at the coast. To better understand the various processes acting at the river basin scale, we utilized empirical mode decomposition to analyze the variability in different parameters (river discharge, rainfall, air temperature, and sea level) from the headwaters to the coast of the Reno and Lamone rivers over the last century. The anthropogenic footprint, linked to the large-scale dimming/brightening phenomenon, is visible in the long-term trends. Moreover, natural signals with variable periodicity are evident and partially correlated with two major climate modes (North Atlantic Oscillation and Atlantic Multidecadal Oscillation). The coupled interactions among these processes, combined with the changes in land use and evapotranspiration during the last century, have resulted in the prolonged scarcity of river sediment supply and a long-term trend of erosion of the coastal area.

1. Introduction

Riverine sediment discharge, the foremost feeding source for low-elevation coastal zones (LECZ) [1,2], has declined drastically worldwide over the last century (e.g., [3,4,5,6,7,8]) and will likely continue to decline [9,10]. The amount and role of sediment supplies at the coast varies with the local climate, geological setting, and degree of human impact within each catchment area [11,12,13]. LECZs worldwide have reached a critical tipping point, shifting from a natural pseudo-stable state since the Holocene to a state controlled by human dynamics [14] owing to land use modifications throughout their catchment areas [15,16]. Furthermore, rising sea levels and other climate-driven effects are expected to exacerbate such impacts on LECZs [17]. Despite climate-driven modifications remaining influential, they may locally be of secondary importance [18] since human effects may more heavily interact with hydrological processes. Anthropogenic modifications over catchment areas, together with the hydrological and climatic changes, and their related connection with atmospheric dynamics, should be thus taken into account in order to better understand which processes are most likely to drive the evolution of coastal areas.
The Emilia-Romagna (ER) coast is a 130 km long LECZ south of the Po River Delta, Northern Italy, facing the northern Adriatic Sea (Figure 1a), which has undergone severe anthropization since the 1950s [19]. Increased subsidence due to underground fluid exploitation [20,21,22,23], land reclamation, and the widespread use of coastal defense structures [24,25,26] have affected the ER coastal plain, making it sensitive to floods [27]. Furthermore, the heavy, anthropogenically driven reduction in sediment delivery by rivers (mainly due to sediment digging and river regulation), mostly between the 1940s and 1980s, has led to generalized erosion and shoreline retreat [28,29,30]. Although safeguard policies for riverbed excavation, introduced by the ER regional government during the early 1980s, were expected to gradually reverse this, only the amount of suspended material has slowly increased, whereas fluvial bedload has decreased [31,32]. Granulometric analysis on samples collected since the 1970s in the submerged beach of the ER coast, at or below 3 m depth, has shown a generalized overall decrease in grain size until the present day [33,34,35,36]. Indeed, a recovery in bedload sediment supply has not yet occurred since the introduction of the safeguard policies; thus, most of the ER shoreline has become threatened by erosion and flooding [34]. Previous analyses have been conducted over the regional scale in order to better elucidate the hydrological, environmental, and climatic evolution of the area, finding a significant decrease in river discharge [37] and indicating that sediments are currently not being transferred downstream as they were previously. While a clear increase in air temperature has also been observed for the last 70 years (e.g., [38]), trends in rainfall have been found to be statistically significant only at certain locations [39,40,41], highlighting a great spatial variability for this parameter. This effect has primarily been linked to the frequencies of diverse weather regimes which, locally, have led to differentiation in both total annual rainfall amount and its variability through time [42]. Generally speaking, significant/non-significant decreases in rainfall are a matter of debate on both the Italian [43,44,45,46] and Mediterranean scales [47,48], for which the contrast has been mainly attributed to the length of the time interval considered in the analysis [45].
In this study, the Reno and Lamone river basins (Figure 1a), which feed an approximately 20 km stretch of the ER coast between the Porto Garibaldi and Porto Corsini harbors (Figure 1b), have been considered. The Reno River, with a length of 210 km and a basin area of approximately 4630 km2, is the main river of the ER region after the Po River, and the sixth nationwide in terms of basin size. It drains both mountains and plains and is joined by many tributaries from the Apennines and artificial channels in the plains (Figure 1c). The basin of the Lamone River (approximately 520 km2, Figure 1a) is much smaller than the Reno River basin; nevertheless, its sediment supply also provides an input for the coastal budget of this area. Coastal prograding at the mouth of the Reno River during the early 20th century was previously associated with a sustained sediment supply [33], whereas from the late 1930s onward, progressive erosion has occurred [34]. The coastal areas of these river basins and adjacent wetlands (the Comacchio Valleys and the Bellocchio marsh; Figure 1b), represent a significant economic and naturalistic resource (Ramsar Convention site); however, these areas are among the most threatened by shoreline retreat and flooding along the Mediterranean and European coasts [49]. Furthermore, these coasts are significantly impacted by economic activities that cause land subsidence, such as gas exploitation by ENI (Ente Nazionale Idrocarburi) in the “Dosso degli Angeli” reservoir (Figure 1b) [50].
In this paper, the linear and nonlinear evolution of environmental-climatic parameters for both river basins have been analyzed by considering river discharge, rainfall, and air temperature datasets, representative of the river basins over the last century. The main goal of the study is to better understand the implications of the hydrological and climatic variations occurring over the river basins and at the coast. According to Preciso et al. [37], the Reno River catchment no longer represents the primary sediment source for the coast, primarily because of reforestation within the catchment since the 1950s. Moreover, in several Mediterranean locations, changes in land use, together with increased evapotranspiration and soil erosion control measures, have reduced runoff, and thus, sediment supply [51,52,53,54]. Due to this, land use changes and evapotranspiration datasets were also taken into account in our analyses in order to obtain a more holistic view of the changes occurring throughout the catchment areas. Moreover, the local sea level, here affected by both natural and anthropogenic components (see [55,56] for details), was considered, as it is a critical factor in coastal evolution.
In detail, this study addresses the following questions: (1) How have anthropogenic and natural signals acted on hydrological processes at the river-basin scale, from the headwater to the coast, over the last century? (2) Why, despite safeguards for river basin management, has the solid river supply to coastal sites not yet recovered? (3) How have climate change and anthropogenic footprints affected sea-level changes in the studied LECZ?

2. Materials and Methods

2.1. Data

The average monthly river discharge, rainfall, and air temperature data were obtained from the Dext3r-SIMC platform [57] and the annals of the National Hydrological Service [58]. Only long-term data within both the river basins were obtained, in order to consider a wider, historical view of the surface processes. Furthermore, liquid transport time series were reconstructed from three selected river gauges (Figure 1c); however, these datasets contained gaps (Figure 2). To rebuild a representative time series, thereby increasing the signal-to-noise ratio, data from the twelve longest and more complete individual rain gauges (out of a total of 92 stations available throughout the two basins) were considered and stacked. Although these datasets were affected by a common gap in the year 1924, minor gaps in each time series were filled using records from the closest station (excluded from the analysis) within 10 km after careful comparison of the overlapping data. To further examine rainfall variability, the number of wet days (≥1 mm/day) were also considered. The minimum and maximum temperatures were used to produce a monthly climatological assessment of the river basins. Additionally, a stacked curve was computed to achieve a basin-representative climate time series that considered data from nine stations (Figure 1c).
The relative sea level data (RSL, which measures the sea level with respect to a benchmark on land) were obtained from the Marina di Ravenna tide gauge (RTG; Figure 1b), which is currently managed by ISPRA (Istituto Superiore per la Protezione e la Ricerca Ambientale). The gauge has recorded data since 1873, with a gap from 1922 to 1933 and discontinuities since 2016 (Figure 2); hence, data after 2015 were discarded. The RTG time series was assembled using data from various sources (1873–1979 and 2001–2016 from the Permanent Service for Mean Sea Level (PSMSL) data bank [59,60]; 1971–2001 from Romagnoli et al. [61], and homogenized considering data overlaps and following Bruni et al. [56]).
River discharge, rainfall, temperature, and RTG time series were compared with three climate indices, namely the North Atlantic Oscillation (NAO), Atlantic Multidecadal Oscillation (AMO), and Western Mediterranean Oscillation (WeMO). The NAO is among the most significant modes of climate variability in the Northern Hemisphere, with seasonally fluctuating manifestations, especially during the boreal winter [62,63]. The winter (December–March) station-based NAO dataset was used; this represents the normalized sea-level pressure between Lisbon (Portugal) and Stykkishólmur (Iceland) since 1864, by the National Centers for Environmental Prediction/National Center for Atmospheric Research [64], based on Hurrell et al. [63]. The AMO [65] is another prominent climate variability mode in the Northern Hemisphere, defined by variability in sea surface temperature (SST) in the North Atlantic Ocean and characterized by a dominant periodicity of approximately 60–70 years. The AMO is the area-weighted average SST (between 0° and 70° N) with the linear trend removed. The annual AMO index, available since 1856, was provided by the National Oceanic and Atmospheric Administration/Physical Sciences Division [66]. The WeMO index [67] is defined as a dipole structure formed (in its positive phase) by an anticyclone over the Azores and a depression over north-western Italy; its index is the result of the difference between the standardized values of surface atmospheric pressure in San Fernando (Spain) and Padua (Italy). In this study, we also accounted for the WeMO since it is considered by some authors more significant than the NAO to explain rainfall anomalies over the north-western Mediterranean (e.g., [68] and references within).
The Standardized Precipitation Evapotranspiration Index (SPEI) [69] was used for long-term drought analysis (1921–2019). The SPEI allows for the identification of drought severity using its length and intensity at various timescales, considering the effects of both temperature and precipitation. This index combines the strength of the multi-temporal Standardized Precipitation Index [70] and the sensitivity of the Palmer Drought Severity Index [71]. The Global SPEI database [72] used in this study covers the globe at a spatial resolution of 0.5 degrees and offers timescales from 1 to 48 months. The potential evapotranspiration in the SPEI database is based on the FAO-56 Penman-Monteith method [73]. Despite its coarser spatial resolution with respect to local datasets, the SPEI has been adopted as it has been thoroughly validated by several authors (e.g., [74,75,76,77,78,79,80,81]) to reveal moisture anomalies for environmental and agricultural applications. Since our study primarily considers the long-term (decadal to centennial) alternation of drought and wet conditions rather than short-term (monthly) events, the 48 month-average time series have been used. The Bologna grid cell was considered as representative of the Lamone and Reno river basins because the variability between adjacent cells at this spatial scale was negligible.
Snow cover percentage over the period 1950–2022 was obtained from the ERA5-Land reanalysis dataset [82,83,84] and limited to the area under analysis. This parameter represents the fraction (0–100%) of the grid cell occupied by snow and is provided with a spatial resolution of 0.1 × 0.1 degrees.
Land use was evaluated using land use maps for 1976, 1994, 2003, 2008, 2014, and 2017, produced by the GIS and Statistic Service of the ER Region [85]. These maps were developed utilizing aerial photography and are characterized by their minimal detectable area: 0.375 (1976), 1.56 (1994 and 2003), and 0.16 (2008, 2014, and 2017) hectares. The land use for these products was mapped using the Corine Land Cover legend until level three.

2.2. Methods

All datasets in this study were homogenized into monthly means. The time series were plotted using a 1-year moving average low-pass filter to highlight the interannual variability. Modified non-parametric Mann-Kendall (MK) tests [86,87,88] were performed to assess temporal trends (statistically significant at the 95% confidence interval). MK analysis was implemented using the Theil-Sen estimator [89,90] to evaluate linear trends. This method is insensitive to outliers and more accurate than simple linear regression [91]. To account for serial autocorrelation and properly estimate trend-related errors, the procedure described in Zervas [92] has been followed.
Considering the stacked rainfall time series and the monthly sum of wet days, the Simple Daily Intensity Index (SDII), which is the ratio of total rainfall to the number of wet days, was computed based on Frich et al. [93]. Only the SDII values above the 95th and 99th percentiles were considered.
For RSL analysis, several authors have addressed the removal or assessment of vertical crustal movements (VLM) from tide gauge data [55,94,95,96,97]. Optimally, VLM contributions can be distinguished by Global Navigation Satellite System monitoring of tide gauge data [98]; however, this began in the late 1990s [99,100]. To remove the VLM components, the RTG time series was subtracted from the Trieste Tide Gauge (TTG; location in Figure 1a) dataset (data 1875–2021 from the PSMSL archive); from this, a locally weighted polynomial regression was created and subsequently removed from the original RTG datasets. This procedure relies on the assumption that TTG location is considered a relatively stable area in terms of vertical land movements, with a weak uplift in the range of a few tenths of mm·year1 or negligible (e.g., [101,102,103,104,105,106]).
To assess nonlinear and non-stationary signals, empirical mode decomposition (EMD) [107] was applied to decompose monthly mean river discharge, rainfall, air temperature, and sea level into finite empirical orthogonal intrinsic mode functions (IMFs). Each IMF describes cyclic variations, though not necessarily constant phases and amplitudes, representative of the oscillating mode from the highest to the lowest frequency. The lasting, non-oscillatory mode is the residual (RES), usually associated with the long-term trend of the signal, which is typically monotonic and provides information about the reliability of the linear model for nonperiodic components. Previous studies have employed EMD to correlate nonlinear variations among different phenomena [108,109,110,111,112,113,114]; moreover, it is well suited for analyzing non-stationary time series.
The non-parametric Kendall rank correlation was employed to examine statistical relationships for NAO, AMO, and WeMO with our time series. This method has significantly smaller gross error sensitivity and asymptotic variance than Spearman rho or Pearson r and provides significantly more accurate p-values [115]. To evaluate autocorrelation, the time series were linearly detrended, and their effective sample size, based on the lag-1 autocorrelation coefficient, was considered [116,117]. The statistical significance of the correlation results was set at 95% CI.
Periodic signals were characterized by computing Lomb-Scargle periodograms [118,119] for each time series and both the NAO and AMO indices, providing a description of the overall frequencies. The periodograms were confined to the standard normalized power interval 0–1, according to Baluev [120], and limited to >4 years to focus on low-frequency signals.
In the land use analysis, only polygons within the Reno and Lamone river basins were considered. A further subdivision was made between the polygons within the alluvial plain and Apennine portions of the maps to highlight the various changes between the two sectors.

3. Results

3.1. Long-Term Linear Analysis

MK tests on all three river gauge stations indicated a statistically significant decrease in river liquid discharge (Figure 3a and Figure A1) over the last century. The Casalecchio di Reno time series, which is the most complete, showed a decrease of −0.18 ± 0.08 m3/s·year−1 during 1921–2021 (Figure 3a). This station was more reliable than that at Ponte Bastia (Figure 1c) since, although the latter could depict the behavior of the lower Reno River basin, it was affected by massive anthropogenic contamination, mainly due to the Cavo Napoleonico (Figure 1c), one of the most significant hydraulic structures in the Po Plain, whose function since the late 1960s has been that of diverting the Reno River water into the Po River during floods and vice versa during drought periods. The river gauge at Pracchia also provides a nearly complete time series of the Reno River runoff (Figure 2). Despite its location in the high Apennine and the relatively limited drainage area upstream, a significant, slight decrease in river discharge occurred. In terms of seasonality (by considering the sole Casalecchio di Reno time series) the monthly mean data show clear and significant decreasing trends, with the lowest discharge values during summer (namely June, July, and August), marked decrease trends for both winter (December, January, and February) and autumn (September, October, and November), and a stronger decrease during spring (March, April, and May) (Figure A2a).
The stacked monthly rainfall time series (Figure 3b) showed a negative but not significant trend (−0.07 ± 0.08 mm·year−1); the same was true for most of the individual rain gauges (Figure A3). Seasonality within the rainfall stacked curve shows the same behavior as that described by Pavan et al. [40], with the minimum values observed during summer and the maximum in autumn (Figure A2b); however, only the other two seasons (winter and spring), comparable in terms of both the amount of precipitation and trends, denote a significant decreasing trend.
The monthly mean temperature trends over the Reno and Lamone river basins during the last century almost all showed statistically significant increases (Figure A4). The stacked temperature curve (Figure 3c) shows a statistically significant increase of 0.016 ± 0.007 °C·year−1, which corresponds to a warming of over 1.5 °C during the last century. The variability of temperature highlights an increasing and significant trend for all seasons (Figure A2c). During autumn, however, the trend is very slight, while stronger warming for the last century is observed in winter. Monthly mean temperatures ranged from 9 to over 14 °C at the highest station in the Apennines (Maresca) and urban areas in the plains (Bologna, Ferrara, and Imola), respectively, showing thermal variability due to both elevation and urbanization. Moreover, the latter has caused an urban heat island effect, which affects local and regional air temperatures near densely settled areas [121]; its main impact is a large increase in nighttime temperatures [122]. In addition, the local topography and land use have led to large variability in temperature, with the minimum (5.5 °C at Maresca) and maximum (18.6 °C at Bologna) averages differing by over 13 °C. However, grouping datasets from similar geographic contexts allows for a gross assessment of the variability induced by land use and local topography. Long-term minimum and maximum averages from the Apennine datasets were 0.2 and 2.4 °C lower than those from the rural plains, respectively. In contrast, maximum averages in the rural plains were not substantially different from those in urban areas (0.16 °C), whereas minimum averages were approximately 1.6 °C higher in urban areas.
The four-year low-pass filtered time series of wet days in Figure 3d shows a nonlinear, cyclic pattern, with minima during the 1920s, 1940s, and the early 1980s–2010. However, the >95th and >99th percentile values of the SDII do not give a clear pattern throughout the century; rather, they suggest an intensification of their occurrence, beginning in the 1970s.
Figure 3e shows the evolution of the SPEI during 1921–2019 in 48-month intervals. Paulo et al. [123] categorized the SPEI drought classes as non-drought (greater than −0.5), mild (−0.5 to −1), moderate (−1.5 to −1), severe (−1.5 to 2), and extreme (less than −2) conditions. In our data, moist and dry periods fluctuated throughout the century, in agreement with the wet days (Figure 3d); however, the wet periods have declined in recent decades. Since the 1980s, drought conditions have been increasingly persistent. Furthermore, extreme drought conditions only occurred during the 1940s and recently.
Trends in ground snow cover percentage (Figure A2d) over the area under study denote a marked and significant decrease in winter (−0.58 ± 0.18 %·year−1) over the period 1950–2022, corresponding to a decrease in the winter snow surface of about 63% in the last 72 years; no significant trends emerged from both autumn and spring.
RSL data from the RTG time series (Figure 4a) clearly indicates a higher rate of sea-level rise during 1940–1980, compared to the prior and subsequent periods. This behavior witnesses the well-known subsidence phenomena due to local human activities (exploitation of gas and water from the underground) which overlapped with the natural land subsidence of the area [124,125]. However, subsidence has slowed significantly since the 1980s, mainly due to safeguarding policies introduced by the regional administration. Accounting for the TTG time series (Figure 4b) as a stable reference at the century scale, Zerbini et al. [126] evaluated the natural subsidence at the RTG site at around 1.88 mm·year−1, by considering the period prior to the start of the mining activity (1873–1922); this value is in agreement with previous geological investigations [127,128,129] and with the first portion of the TTG-RTG time series of Figure 4c. The estimated trend of the RTG time series (1875–2016), detrended for VLMs (Figure 4d) is 1.3 ± 0.2 mm·year−1, which is comparable with the TTG estimate of 1.33 ± 0.17 mm·year−1 (1875–2021).

3.2. Nonlinear Signal

The EMD-decomposed time series of river discharge, rainfall, and temperature are shown in Figure 5a–c respectively, along with the sea level time series from the corrected RTG (Figure 5d). EMD decomposition produced seven IMFs for each of these four time series; however, to focus on long-term variability, only the lowest frequency IMFs (IMF4, IMF5, and IMF6) and RES were considered.
The RES component in river discharge was not perfectly linear but decreased over time; this decrease was 90% and 190% greater during 1993–2021 compared to that during 1950–1980 and 1921–1950, respectively. The RES in both rainfall and temperature show non-monotonic behavior, suggesting the presence of cycles longer than the time series themselves. In detail, rainfall decreased during 1950–2000 but otherwise increased. In contrast, temperature has increased since the 1920s; although it dropped slightly during the 1960s and 1970s, it has strongly increased since the 1980s. Conversely, the RES in sea level consistently increased over the last 140 years. RES trends (Table 1) were statistically significant in all cases, even for rainfall, in contrast to the original time series (see Section 4.1).
The time series compared with the NAO and AMO climate indices (Figure 5e,f) show significant (95% CI) negative correlations of NAO phases with river discharge, rainfall, and sea level of −0.31, −0.16, and −0.34, respectively. In contrast, air temperature was positively correlated with both NAO (0.26) and AMO (0.21) fluctuations, and the sea level correlated with the AMO (0.14) index. However, only the air temperature-NAO correlation was statistically significant at the 95% CI, while the correlations with the AMO index become significant only if the 90% CI is considered. No correlations were found between the AMO with either river discharge or rainfall, nor was the WeMO index with any variable considered. Since no correlation was found on any annual time series representative of this area, the latter is not further discussed in this paper.
Several periodicities “contaminate” the time series, as detected by the LSP analysis (Table 2). Most are common to all the datasets, albeit with different values. Periodicities of 4–5 and 8–9 years were detected in all the time series and are dominant within the IMF4 of rainfall and temperature and the IMF4 of sea level, respectively (Figure 5). Additionally, an approximately 6-year periodicity, which is dominant in the IMF4 of river discharge, is not present in the temperature and sea level. At lower frequencies, only the approximately 12-year period is common among all the time series, while periods of 10–11, approximately 14, and 16–17 years are dominant within IMF5 of rainfall, river discharge and temperature, and sea level, respectively. A periodicity of 22–23 years is common to all time series (dominant in the IMF6 of river discharge and temperature) when considering >20 years, whereas in rainfall and sea level, a powerful periodicity of approximately 40 years characterizes IMF6. Furthermore, two periodicities prevail at the longest wavelength: 53–56 years for river discharge and temperature, and approximately 67 years for sea level. All these periodicities occur in the NAO and AMO time series (Table 2), often with values common to both indices, and fluctuate in phase or anti-phase (as described above) with the IMFs of all the decomposed time series.

3.3. Land Use Changes

Recent land use changes within the Reno and Lamone River basins were categorized using the first level of the Corine Land Cover code as artificial areas, agricultural areas, forests and semi-natural areas, wetlands, and water bodies (see also Heymann et al. [130]). The land use changes in the Apennines for both basins (Figure 6) highlight that agricultural areas have been replaced since the 1970s, mainly by forests and semi-natural areas (approximately 70%) and some artificial surfaces (30%). Recent reforestation in the upper ER river basins has been driven by the migration of rural populations to urban areas and consequent cropland abandonment [32,34].
Trends in the plain sector have been nearly constant since the second postwar period, when the landscape began changing, following a rapid shift from a largely peasant economy to an industrial and touristic economy [131,132]. This has caused a gradual reduction in agricultural areas and intensive urban sprawl. This “uncontrolled spread of towns and villages into undeveloped areas” [133] has increased soil sealing due to impermeable artificial materials, which are considered the largest issue for pervasive soil degradation [134]. Moreover, the area of sealed soils (excluding areas such as gardens, urban wastelands, and parks) has increased by approximately 123% since 1976. Thus, sealed surfaces accounted for approximately 83% of the total artificial areas in all six analyzed maps. Figure 6 shows that forests and wetlands represent very little of the land use in the plains, which is a consequence of human activities beginning with progressive deforestation by the Romans during the first century B.C. and land reclamation of swamps and salt marshes, begun in the 1800s and completed in 1964. Finally, the slight increase in water bodies in the plains was primarily due to progressive increases in artificial basins (over 120%) and, to a lesser extent, in waterways and irrigation canals (approximately 45% increase).

4. Discussion

4.1. Climatic and Hydrological Variability Induced by the Anthropogenic Dimming/Brightening Phenomenon

Surface water runoff has decreased in several areas of the ER region in recent decades [135]. Local water scarcity has been driven by the onset of persistent drought conditions, which began in the 1980s and have worsened, as also suggested by the SPEI (Figure 3e); conditions are projected to further deteriorate over the next decades [136,137].
In our analysis, air temperature (Figure 3c) showed an increasing trend over the last century, as previously observed throughout the ER region by Antolini et al. [38]. However, the RES curve of temperature (Figure 5c) showed an increase from the 1920s to 1950s, a generally stable (or slight cooling) period until the late 1970s, and rapid warming thereafter. This general behavior in temperature, as previously observed throughout Italy [138,139,140], Europe [141,142], and worldwide [143], has been linked to a dimming/brightening phenomenon [144,145] from surface solar radiation (SSR) changes over the last century [146,147] driven by anthropogenic air pollution, which altered the transparency of the atmosphere [148,149,150]. The marked change around 1980 led to a “global warming-type drought” [151,152]; its effects on water evaporation and local temperature have been further accentuated by the presence of artificial surfaces.
Despite the lack of statistical significance (Figure A3), the decrease in rainfall over the last century over large portions of the ER was mainly attributed to a decrease in spring and winter rainfall [38,39,40,153,154], as confirmed by our analysis (Figure A2b). However, when the oscillations are removed, a significant decrease and non-monotonic behavior are visible (Figure 5b). The latter is likely due to the anthropogenic dimming/brightening process described above, as any change in SSR leads to large modifications in the water cycle [155,156]. Wild [157] have found evidence of concurrent decreases in rainfall and SSR from the 1950s to 1980s across the Northern Hemisphere, followed by rainfall intensification from the 1980s onward because of subsequent brightening and increased evapotranspiration [158]. This is consistent with our results, although rainfall resumed long after the 1980s (Figure 5b). This may be due to climatological moisture divergence over the Mediterranean, which transports water vapor outside the region; therefore, moisture does not fully translate to regional increases in rainfall [159,160,161]. Moreover, this phenomenon is likely connected to a possible increase in rainfall extremes [161,162], as also suggested by the increase in the >95% SDII values from the 1970s (Figure 3d). This could be strictly linked to the concordant reduction in rainy days, and in this case, not necessarily associated to a substantial decrease in rainfall totals. It can be argued that, following this decline in the 1970s, the amount of total monthly rainfall has not changed considerably (Figure 3b and Figure 5b), and therefore falls in a shorter time.

4.2. Anthropogenic Influences at the River-Basin Scale

The river discharge time series in this study shows a significant decreasing trend since the 1920s (Figure 3a and Figure A1), in agreement with Preciso et al. [37]. Unfortunately, data on bedload materials are lacking, and the relationship between them and river flow discharge cannot be verified [32]. As stated in Section 1, the Reno catchment no longer represents the primary sediment source for the coast, primarily because of reforestation in the catchment since the 1950s [37]. This was confirmed by our analysis (Figure 6); increased forest cover represents the prevailing land use change in the Apennines portion of both basins (see [163]), which likely influenced the reduction in discharge and sediment production. The decrease in discharge showed a negative acceleration in the RES curve (Figure 5a) during 1950–1980, compared to the previous 30 years, which is due to the influence of both river regulation and reforestation, as pastures, shrubs, and crops, consume much less water than forests [164,165,166]. However, the huge decrease during the following years (1980–2021; Table 1 and Section 4.2) cannot be attributed to river regulation. First, the increase in reforested areas has remained constant, although it has slowed since the mid-1990s. Second, the notable increase in temperature since the 1980s (Figure 5c) has significantly increased the atmospheric evaporative demand and led to persistent local drought conditions (Figure 3e), which, coupled with evapotranspiration, may have driven further declines in discharge. The non-monotonic variation in rainfall (Figure 5b) also contributed to the decreased river discharge since the 1950s; however, its significant but slight changes played a secondary role to the above processes. Furthermore, another consequence of the increasing temperatures may be the decreasing extent of snow cover and number of days with snow present on the Italian peninsula [41] and throughout the Northern Hemisphere [167] over the last century, especially since the late 1970s. Declining snow cover extent throughout the upper catchment areas during winter (Figure A2d) could be a further driver of the reduced river discharge by weakening or even depriving the fluvial systems of water supplied by snowmelt [168,169]. The decline in snow cover extent throughout the upper portion of the studied catchments area can be ascribed as the main driver for the deterioration in spring trends of river discharge (Figure A2a), by weakening, or even depriving, the fluvial systems of the water supply due to the melting of snow.
Once eroded, soil particles can undergo transportation-deposition cycles throughout their travel downstream and either remain within the basin or outflow from it. Thus, processes throughout the plains, usually of anthropogenic origin, are crucial for sediment supply. Our analysis noted that agricultural areas are substantially decreasing in the plains of both basins (Figure 6) and are being replaced by artificial surfaces and, to a lesser extent, water bodies. The increase in the latter, despite the small area, may have significantly influenced discharge by withholding water from the system. Artificial basins (+120% during 1976–2017) and irrigation channels (+45%) have been continuously implemented throughout the plains to satisfy increasing water demand. Intensive groundwater withdrawals in these areas, mainly for industrial and domestic needs [170], also contributed significantly to the decreased discharge, as over-exploitation of aquifers weakens (and often prevents) natural river recharge from the aquifers. The increased temperatures and UHI effect also directly influence the evaporative processes and water demand for both agriculture and domestic needs in the plains. These processes have been magnified during recent decades owing to urban sprawl, the increase in artificial surfaces (Figure 6), and soil sealing (+123% during 1976–2017). The latter reduces and sometimes prevents water infiltration into soils and, therefore, aquifer recharge, increasing both water scarcity and dependence on irrigation and artificial reservoirs [171].

4.3. Natural Contributions from Coupled Ocean/Atmosphere Processes

The NAO influences the weather and hydrology by controlling the wind intensity and direction and the interactions among air masses over much of Europe, and the AMO influences several hydroclimatic variables, including Mediterranean air temperature (e.g., [172,173,174,175]). Additionally, the sea level in the Mediterranean is influenced by the NAO, which drives atmospheric sea-level pressure changes [63,176] and the net water flux at Gibraltar by directly altering wind and oceanic circulation near the strait [177,178] and indirectly by impacting regional river runoff, evaporation, and rainfall [179]. Moreover, the AMO influences the SST and sea evaporation [174,179].
Previous studies have found significant correlations between the NAO and Mediterranean surface air temperatures [62,174,180,181,182,183]. Our findings of a positive correlation between temperature and AMO cyclicity (Table 2) agree; however, they cannot be confirmed because of the statistical non-significance of the correlation at the 95% CI. Mariotti and Dell’Aquila [174] have also found no correlation between the AMO variability and rainfall over the Mediterranean; our study confirms this apparent lack of influence.
Several authors have found significant negative correlations between rainfall and NAO over Southern Europe and the Mediterranean, with the positive phases of the NAO linked to drier conditions and vice versa [62,63,139,159,181,184,185]. Generally, positive and negative NAO phases decrease and increase rainfall, respectively, over Southern Europe through anomalous behavior of the Atlantic westerlies [184,186], thereby severely impacting regional hydrological cycles [159]. This is reflected by river discharge variability, as indicated by our analysis, which shows a significant negative correlation between NAO phases and river runoff (Figure 5a,e). Negative correlations between river discharge and NAO have been found throughout Europe [187,188,189,190]. However, López-Moreno et al. [191] note that, unlike rainfall, river discharge is characterized by a nonlinear response of short, but intense, discharge anomalies during negative NAO phases and persistent, but weaker, anomalies during positive phases. The same authors highlight that during positive NAO years, surface runoff and aquifer recharge are severely reduced, whereas during negative NAO years, river basins receive more rainfall, leading to soil saturation and rapid rainfall transfer to runoff. Positive anomalies, such as the strong peak observed in river discharge in 2010 (Figure 5a), may be a consequence of the negative NAO phase that occurred around 2010, which might have left a record in the river basins and coastal systems. Granulometric analysis performed by ARPAE (Agenzia Regionale per la Prevenzione, l’Ambiente e l’Energia—Emilia Romagna) on samples collected in 2012 in the submerged beach of the studied coastal tract [25] showed larger grain sizes than those from later [36] or prior [33,34,35] regional sampling. Excluding the possible effects of local interventions (i.e., beach replenishment) [36], these variations in grain size may be related to the large, positive anomaly in river discharge that occurred in 2010. This natural anomaly could have momentarily strengthened the solids’ contribution to the coast by enhancing the ability of rivers to carry coarser sediment across the submerged beach through the river mouth bypassing and subsequent redistribution by longshore currents. Thus, the generalized increase in grain size observed in the studied coastal stretch and elsewhere in the nearshore ER coast in 2012 could have been associated with this short-term event, which temporarily disrupted the tendency towards gradually finer and weaker sediment supply observed in recent decades.
At the ER coast, persistent sea-level rise (Figure 5d), periodically altered by various cyclicities, inevitably impacted shoreline evolution. The large positive sea-level anomaly shown in the IMF4 in Figure 5d, has also been observed in previous studies throughout the Mediterranean [178,192,193] and has been attributed to a strongly negative NAO [194]. The dominant periodicity observed in sea-level fluctuations was approximately 8–9 years, which is comparable to the approximately 10-year periodicity found by Bonaduce et al. [112]. Additionally, Galassi and Spada [111] associated the 2010 sea-level anomalies with the coincidence of both negative NAO and positive AMO phases, which has a periodicity of approximately 20 years.
Since the 1990s, sea level has risen rapidly throughout the Mediterranean, which has been linked to increased SST [195], mostly controlled by increasing air temperatures. This also happened during the early 20th century, the other major period of constant rise in global sea and air temperatures [196]. The sea level oscillations at the RTG during these periods are described by IMF6 (Figure 5d), which has a dominant periodicity of approximately 40 years (Table 2). Indeed, major sea-level changes at the RTG over the last 140 years (excluding the effects of local VLMs) were described by multiple phases of acceleration, during the early 20th century, the 1970s–1980s, and 1990s to mid-2000s. Moreover, this behavior is highly phase-correlated with the rainfall IMF6 (Figure 5b), which is also characterized by a similarly dominant periodicity. The NAO may be the driver of this process, because it has a similar periodicity (approximately 37–40 years) and a good negative correlation with both sea level and rainfall. For the same study area, Meli et al. [114] noticed a slight negative acceleration in sea level change since the mid-2000s and attributed it to the possible presence of periodic signals at periods longer than the 26-year period they considered. Therefore, this recent change in the local sea level trend could be associated with the above periodicities, as the final portion of the sea-level IMF6 denotes a temporally concordant, progressive flattening of the rising trend that began in the early 1990s.

5. Conclusions

The environmental and climatic data considered in this study provide a holistic view of the hydrological and climatic processes that have affected the Reno and Lamone river basins of the ER coastal areas in Northern Italy over the last century. The analysis provided responses to the questions posed in Section 1.
  • The anthropogenic footprint, attributed to effects such as land use changes and the large-scale anthropogenic dimming/brightening phenomenon, acting at both the river basin and regional-to-global scales, profoundly impacts the catchment dynamics by driving long-term, nonlinear trends, upon which natural oscillations are superimposed. Interactions with major climate modes, in fact, affect the signals over various periodicities. Both positive and negative correlations among some of the studied parameters and the main climatic indexes (NAO, AMO, WeMo) are evidenced;
  • The marked negative acceleration in river discharge provides an indication of the lack of recovery in terms of sediment supply to the coast, despite safeguard policies introduced by the ER regional administration in the early 1980s. This decline resulted from river regulation and land use changes during 1950–1980 and related implications. Moreover, since the 1980s, local air temperatures have increased significantly, leading to persistent drought conditions contributing to the drastic reduction in river discharge. However, periodic natural signals can significantly restore the river discharge, as observed during the strongly negative NAO event in 2010;
  • Persistent sea-level rise has affected the coastal site under study over the last 140 years. Apart from the extremely high rates during the 1950s–1980s due to anthropogenic-induced land subsidence, sea level is periodically amplified or reduced by natural fluctuations, such as that observed in 2010, and lower frequency fluctuations in the local sea level, as during the early 20th century, 1970s–1980s, and 1990s to mid-2000s. This nonlinear behavior, locally enhanced by subsidence that overwhelms the river sediment inputs, notably impacted the ER coast and may continue to do so.
Overall, the interconnection found between various climatic, hydrological, and anthropogenic processes with individual nonlinear behaviors illustrates the need for a holistic approach when considering complex environments such as LECZs and heavily anthropized areas. Expanded observation and monitoring at the scale of entire river basins is necessary to better understand the natural and anthropogenic influences on coastal behavior. This approach should be thoroughly considered in coastal management and in adaptation strategies, especially when facing a changing climate and increased vulnerability in the most threatened environments, such as LECZs.

Author Contributions

Conceptualization, M.M. and C.R.; methodology, M.M.; formal analysis, M.M.; investigation, M.M. and C.R.; data curation, M.M.; writing—original draft preparation, M.M. and C.R; writing—review and editing, M.M. and C.R.; visualization, M.M.; supervision, C.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. Ph.D. fellowship of M.M. is financed through the Strategic Development Project of the BiGeA Department, University of Bologna.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Nunzio De Nigris, Maurizio Morelli, and Margherita Aguzzi from ARPAE (Agenzia Regionale per la Prevenzione, l’Ambiente e l’Energia—Emilia Romagna) for making coastal data available. We thank Giovanni Liguori for helpful feedback on technical aspects, Erica and Anna Visibelli for valuable help on Python codes, Marco Olivieri for his help on the EMD analysis and Victor Malagon Santos for his valuable suggestions on an earlier version of this paper. Three anonymous reviewers are acknowledged for their fruitful comments and suggestions that helped improving the paper. This work was carried out in the framework of the Ph.D. course in Future Earth, Climate Change and Societal Challenges (University of Bologna).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Historical time series of Reno River discharge, and related linear trends, at Pracchia and Ponte Bastia (for locations see Figure 1).
Figure A1. Historical time series of Reno River discharge, and related linear trends, at Pracchia and Ponte Bastia (for locations see Figure 1).
Water 14 02650 g0a1
Figure A2. The seasonal trends computed for river discharge (a), rainfall (b), air temperature (c), and snow cover (d).
Figure A2. The seasonal trends computed for river discharge (a), rainfall (b), air temperature (c), and snow cover (d).
Water 14 02650 g0a2
Figure A3. Rainfall time series, and related linear trends, for the twelve rain gauge stations considered in this work (see Figure 1 for locations).
Figure A3. Rainfall time series, and related linear trends, for the twelve rain gauge stations considered in this work (see Figure 1 for locations).
Water 14 02650 g0a3
Figure A4. Monthly maximum (red), minimum (blue), and mean (yellow) temperature time series, with estimated linear trends (black dashed line computed over the monthly mean) for the nine stations considered in this work.
Figure A4. Monthly maximum (red), minimum (blue), and mean (yellow) temperature time series, with estimated linear trends (black dashed line computed over the monthly mean) for the nine stations considered in this work.
Water 14 02650 g0a4

References

  1. Nicholls, R.J.; Wong, P.P.; Burkett, V.; Codignotto, J.; Hay, J.; McLean, R.; Ragoonaden, S.; Woodroffe, C.D. Coastal systems and low-lying areas. In Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Parry, M.L., Canziani, O.F., Palutikof, J.P., Van Der Linden, P., Hanson, C.E., Eds.; Cambridge University Press: Cambridge, UK, 2007; pp. 315–357. [Google Scholar]
  2. Syvitski, J.; Kettner, A.; Overeem, I.; Hutton, E.; Hannon, M.; Brakenridge, G.R.; Day, J.; Vörösmarty, C.; Saito, Y.; Giosan, L.; et al. Sinking deltas due to human activities. Nat. Geosci. 2009, 2, 681–686. [Google Scholar] [CrossRef]
  3. Llorens, P.; Queralt, I.; Plana, F.; Gallart, F. Studying solute and particulate sediment transfer in a small Mediterranean mountainous catchment subject to land abandonment. Earth Surf. Process. Landf. 1997, 22, 1027–1035. [Google Scholar] [CrossRef]
  4. Syvitski, J.P.M.; Harvey, N.; Wolanski, E.; Burnett, W.C.; Perillo, G.M.E.; Gornitz, V.; Arthurton, R.K.; Bokuniewicz, H.; Campbell, J.W.; Cooper, L.; et al. Dynamics of the coastal zone. In Coastal Fluxes in the Anthropocene; Crossland, C.J., Kremer, H.H., Lindeboom, H., Marshall Crossland, J.I., Le Tissier, M.D.A., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; pp. 39–94. [Google Scholar]
  5. Walling, D.E. Human impact on land-ocean sediment transfer by the world’s rivers. Geomorphology 2006, 79, 192–216. [Google Scholar] [CrossRef]
  6. Syvitski, J.P.M.; Kettner, A. Sediment flux and the Anthropocene. Phil. Trans. R. Soc. A 2011, 369, 957–975. [Google Scholar] [CrossRef]
  7. Weston, N.B. Declining Sediments and Rising Seas: An Unfortunate Convergence for Tidal Wetlands. Estuaries Coasts 2014, 37, 1–23. [Google Scholar] [CrossRef]
  8. Buendia, C.; Bussi, G.; Tuset, J.; Vericat, D.; Sabater, S.; Palau, A.; Batalla, R.J. Effects of afforestation on runoff and sediment load in an upland Mediterranean catchment. Sci. Total Environ. 2016, 540, 144–157. [Google Scholar] [CrossRef]
  9. Zarfl, C.; Lumsdon, A.E.; Berlekamp, J.; Tydecks, L.; Tockner, K. A global boom in hydropower dam construction. Aquat. Sci. 2015, 77, 161–170. [Google Scholar] [CrossRef]
  10. Dunn, F.E.; Darby, S.E.; Nicholls, R.J.; Cohen, S.; Zarfl, C.; Fekete, B.M. Projections of declining fluvial sediment delivery to major deltas worldwide in response to climate change and anthropogenic stress. Environ. Res. Lett. 2019, 14, 084034. [Google Scholar] [CrossRef]
  11. Milliman, J.D. Delivery and fate of fluvial water and sediment to the sea: A marine geologist’s view of European rivers. Sci. Mar. 2001, 65, 121–132. [Google Scholar] [CrossRef]
  12. Vörösmarty, C.J.; Meybeck, M.; Fekete, B.; Sharma, K.; Green, P.; Syvitski, J.P.M. Anthropogenic sediment retention: Major global impact from registered river impoundments. Glob. Planet. Chang. 2003, 39, 169–190. [Google Scholar] [CrossRef]
  13. Syvitski, J.P.M.; Vörösmarty, C.J.; Kettner, A.J.; Green, P. Impact of humans on the flux of terrestrial sediment to the Global Coastal Ocean. Science 2005, 308, 376–380. [Google Scholar] [CrossRef] [PubMed]
  14. Tessler, Z.D.; Vörösmarty, C.J.; Grossberg, M.; Gladkova, I.; Aizenman, H.; Syvitski, J.P.M.; Foufoula-Georgiou, E. Profiling risk and sustainability in coastal deltas of the world. Science 2015, 349, 638–643. [Google Scholar] [CrossRef] [PubMed]
  15. Milliman, J.D.; Syvitski, J.P.M. Geomorphic/tectonic control of sediment discharge to the ocean: The importance of small mountainous rivers. J. Geol. 1992, 100, 525–544. [Google Scholar] [CrossRef]
  16. Darby, S.E.; Hackney, C.R.; Leyland, J.; Kummu, M.; Lauri, H.; Parson, D.R.; Best, J.L.; Nicholas, A.P.; Aalto, R. Fluvial sediment supply to a mega-delta reduced by shifting tropical-cyclone activity. Nature 2016, 539, 276–279. [Google Scholar] [CrossRef]
  17. Fox-Kemper, B.; Hewitt, H.T.; Xiao, C.; Aðalgeirsdóttir, G.; Drijfhout, S.S.; Edwards, T.L.; Golledge, N.R.; Hemer, M.; Kopp, R.E.; Krinner, G.; et al. Ocean, Cryosphere and Sea Level Change. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M.I., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2021; pp. 1211–1362. [Google Scholar]
  18. Darby, S.E.; Dunn, F.E.; Nicholls, R.J.; Rahman, M.; Riddy, L. A first look at the influence of anthropogenic climate change on the future delivery of fluvial sediment to the Ganges-Brahmaputra-Meghna delta. Environ. Sci. Process. Impact 2015, 17, 1587–1600. [Google Scholar] [CrossRef]
  19. Lorito, S.; Calabrese, L.; Perini, L.; Cibin, U. Uso del suolo della costa. In Il Sistema Mare-Costa dell’Emilia-Romagna; Perini, L., Calabrese, L., Eds.; Edizioni Pendragon: Bologna, Italy, 2010; pp. 109–118. [Google Scholar]
  20. Carbognin, L.; Gatto, P.; Mozzi, G. Case history no.9.15: Ravenna, Italy. In Guidebook to Studies of Land Subsidence due to Ground-Water Withdrawal; Poland, J.F., Ed.; UNESCO: Paris, France, 1984; pp. 291–305. [Google Scholar]
  21. Bertoni, W.; Brighenti, G.; Gambolati, G.; Ricceri, G.; Vullermin, F. Land subsidence due to gas production in the on- and off-shore natural gas fields of the Ravenna area, Italy. In Proceedings of the 5th International Symposium on Land Subsidence, The Hague, The Netherlands, 16–20 October 1995. [Google Scholar]
  22. Gambolati, G.; Teatini, P.; Tomasi, L.; Gonella, M. Coastline regression of the Romagna Region, Italy, due to sea level rise and natural and anthropogenic land subsidence. Water Resour. Res. 1999, 35, 163–184. [Google Scholar] [CrossRef]
  23. Teatini, P.; Ferronato, M.; Gambolati, G.; Bertoni, W.; Gonella, M. A century of land subsidence in Ravenna, Italy. Environ. Geol. 2005, 47, 831–846. [Google Scholar] [CrossRef]
  24. Elfrink, B.; Christensen, E.D.; Broker, I. Coastal morphodynamics in subsiding areas. In CENAS: Coastline Evolution of the Upper Adriatic Sea Due to Sea Level Rise and Natural and Anthropogenic Land Subsidence; Gambolati, G., Ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1998; pp. 235–270. [Google Scholar]
  25. Aguzzi, M.; Bonsignore, F.; De Nigris, N.; Morelli, M.; Paccagnella, T.; Romagnoli, V.; Unguendoli, S. Stato del Litorale Emiliano-Romagnolo al 2012. Erosione e Interventi di Difesa; Arpae Emilia-Romagna: Bologna, Italy, 2016; p. 230. [Google Scholar]
  26. Billi, P.; Salemi, E.; Preciso, E.; Ciavola, P.; Armaroli, C. Field measurement of bedload in a sand-bed river supplying a sediment starving beach. Z. Für Geomorphol. 2017, 61, 207–223. [Google Scholar] [CrossRef]
  27. Antonioli, F.; Anzidei, M.; Amorosi, A.; Lo Presti, V.; Mastronuzzi, G.; Deiana, G.; De Falco, G.; Fontana, A.; Fontolan, G.; Lisco, S.; et al. Sea-level rise and potential drowning of the Italian coastal plains: Flooding risk scenarios for 2100. Quat. Sci. Rev. 2017, 158, 29–43. [Google Scholar] [CrossRef]
  28. Dal Cin, R. I litorali del delta del Po e alle foci dell’Adige e del Brenta: Caratteri tessiturali e dispersione dei sedimenti, cause dell’arretramento e previsioni sull’evoluzione futura. Boll. Soc. Geol. Ital. 1983, 102, 9–56. [Google Scholar]
  29. Simeoni, U.; Bondesan, M. The role and responsibility of man in the evolution of the Adriatic alluvial coasts of Italy. In Transformations and Evolution of the Mediterranean Coastline; Briand, F., Maldonado, A., Eds.; Musée Océanographique: Monaco-Ville, Monaco, 1997; pp. 111–132. [Google Scholar]
  30. Preti, M.; De Nigris, N.; Morelli, M.; Monti, M.; Bonsignore, F.; Aguzzi, M. Stato del Litorale Emiliano-Romagnolo All’anno 2007 e Piano Decennale di Gestione; Arpa Emilia-Romagna: Bologna, Italy, 2008. [Google Scholar]
  31. Ente Regionale per le Politiche Ambientali (IDROSER). Il Trasporto Solido Fluviale Nei Bacini Tributari Dell’Adriatico. Regione Emilia-Romagna, Piano Progettuale per la Difesa della Costa Emiliano-Romagnola; Regione Emilia-Romagna: Bologna, Italy, 1983; p. 429. [Google Scholar]
  32. Agenzia Regionale per la Protezione Ambientale (ARPA). Stato Del Litorale Emiliano-Romagnolo All’anno 2000; Arpa Emilia-Romagna: Bologna, Italy, 2002. [Google Scholar]
  33. Bondesan, M.; Calderoni, G.; Dal Cin, R. Il litorale delle province di Ferrara e Ravenna (Alto Adriatico); evoluzione morfologica e distribuzione dei sedimenti. Bol. Soc. Geol. Ital. 1978, 97, 247–287. [Google Scholar]
  34. Ente Regionale per le Politiche Ambientali (IDROSER). Progetto di Piano per la Difesa dal Mare e la Riqualificazione Ambientale del Litorale della Regione Emilia-Romagna, Relazione Generale; Regione Emilia-Romagna: Bologna, Italy, 1996; p. 365. [Google Scholar]
  35. Simeoni, U.; Atzeni, P.; Bonora, N.; Borasio, E.; Del Grande, C.; Gabbianelli, G.; Gonella, M.; Tessari, U.; Valpreda, E.; Zamariolo, A. Integrated Management Study of Comacchio Coast (Italy). J. Coast. Res. 2002, 36, 686–693. [Google Scholar] [CrossRef]
  36. Aguzzi, M.; Costantino, R.; De Nigris, N.; Morelli, M.; Romagnoli, C.; Unguendoli, S.; Vecchi, E. Stato del Litorale Emiliano-Romagnolo al 2018. Erosione e Interventi di Difesa; Arpae Emilia-Romagna: Bologna, Italy, 2020; p. 222. [Google Scholar]
  37. Preciso, E.; Salemi, E.; Billi, P. Land use changes, torrent control works and sediment mining: Effect on channel morphology and sediment flux, case study of the Reno River (Northern Italy). Hydrol. Processes 2012, 26, 1134–1148. [Google Scholar] [CrossRef]
  38. Antolini, G.; Auteri, L.; Pavan, V.; Tomei, F.; Tomozeiu, R.; Marletto, V. A daily high-resolution gridded climatic data set for Emilia-Romagna, Italy, during 1961–2010. Int. J. Clim. 2016, 36, 1970–1986. [Google Scholar] [CrossRef]
  39. Pavan, V.; Tomozeiu, R.; Cacciamani, C.; Di Lorenzo, M. Daily precipitation observations over Emilia-Romagna: Mean values and extremes. Int. J. Climatol. 2008, 28, 2065–2079. [Google Scholar] [CrossRef]
  40. Pavan, V.; Antolini, G.; Barbiero, R.; Berni, N.; Brunier, F.; Cacciamani, C.; Cagnati, A.; Cazzuli, O.; Cicogna, A.; De Luigi, C.; et al. High resolution climate precipitation analysis for north-central Italy, 1961–2015. Clim. Dyn. 2019, 52, 3435–3453. [Google Scholar] [CrossRef]
  41. Fioravanti, G.; Piervitali, E.; Desiato, F.; Perconti, W.; Fraschetti, P. Variazioni e Tendenze Degli Estremi di Temperature e Precipitazione in Italia; ISPRA: Roma, Italy, 2013; p. 60. [Google Scholar]
  42. Cacciamani, C.; Nanni, S.; Tibaldi, S. Mesoclimatology of winter temperature and precipitation in the Po Valley of Northern Italy. Int. J. Climatol. 1994, 14, 777–814. [Google Scholar] [CrossRef]
  43. Brunetti, M.; Maugeri, M.; Monti, F.; Nanni, T. Temperature and precipitation variability in Italy in the last two centuries from homogenized instrumental time series. Int. J. Climatol. 2006, 26, 345–381. [Google Scholar] [CrossRef]
  44. Desiato, F.; Fioravanti, G.; Fraschetti, P.; Piervitali, E.; Pavan, V. Gli Indicatori del Clima in Italia nel 2014; ISPRA: Roma, Italy, 2015. [Google Scholar]
  45. Billi, P.; Fazzini, M. Global change and river flow in Italy. Glob. Planet. Chang. 2017, 155, 234–246. [Google Scholar] [CrossRef]
  46. Toreti, A.; Fioravanti, G.; Perconti, W.; Desiato, F. Annual and seasonal precipitation over Italy from 1961 to 2006. Int. J. Climatol. 2009, 29, 1976–1987. [Google Scholar] [CrossRef]
  47. Trenberth, K.E.; Jones, P.D.; Ambenje, P.; Bojariu, R.; Easterling, D.; Klein Tank, A.; Parker, D.; Rahimzadeh, F.; Renwick, J.A.; Rusticucci, M.; et al. Observations: Surface and atmospheric climate change. In Climate Change 2007: The Physical Science Basis, Contribution of Working Group I to the Fourth Assessment Report of the IPCC; Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K.B., Tignor, M., Miller, H.L., Eds.; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  48. Norrant, C.; Douguedroit, A. Monthly and daily precipitation trends in the Mediterranean (1950–2000). Theor. Appl. Climatol. 2006, 83, 89–106. [Google Scholar] [CrossRef]
  49. Eurosion. Living with Coastal Erosion in Europe: A Guide to Coastal Erosion Management Practices in Europe: Lessons Learned. A Case Study of Marina di Ravenna and Lido Adriano; Directorate General Environment European Commission: Brussels, Belgium, 2004. [Google Scholar]
  50. Comerci, V.; Vittori, E. The Need for a Standardized Methodology for Quantitative Assessment of Natural and Anthropogenic Land Subsidence: The Agosta (Italy) Gas Field Case. Remote Sens. 2019, 11, 1178. [Google Scholar] [CrossRef] [Green Version]
  51. Lasanta, T.; Beguería, S.; García-Ruiz, J.M. Geomorphic and hydrological effects of traditional shifting agriculture in a Mediterranean mountain, Central Spanish Pyrenees. Mt. Res. Dev. 2006, 26, 146–152. [Google Scholar] [CrossRef]
  52. Falcucci, A.; Maiorano, L.; Boitani, L. Changes in land-use/land-cover patterns in Italy and their implications for biodiversity conservation. Landsc. Ecol. 2007, 22, 617–631. [Google Scholar] [CrossRef]
  53. Béjar, M.; Vericat, D.; Batalla, R.J.; Gibbins, C.N. Variation in flow and suspended sediment transport in a montane river affected by hydropeaking and instream mining. Geomorphology 2018, 310, 69–83. [Google Scholar] [CrossRef]
  54. Spalevic, V.; Barovic, G.; Vujacic, D.; Curovic, M.; Behzadfar, M.; Djurovic, N.; Dudic, B.; Billi, P. The Impact of Land Use Changes on Soil Erosion in the River Basin of Miocki Potok, Montenegro. Water 2020, 12, 2973. [Google Scholar] [CrossRef]
  55. Zerbini, S.; Raicich, F.; Prati, C.M.; Bruni, S.; Del Conte, S.; Errico, M.; Santi, E. Sea-level change in the Northern Mediterranean Sea from long-period tide gauge time series. Earth Sci. Rev. 2017, 167, 72–87. [Google Scholar] [CrossRef]
  56. Bruni, S.; Zerbini, S.; Raicich, F.; Errico, M. Rescue of the 1873-1922 high and low waters of the Porto Corsini/Marina di Ravenna (northern Adriatic, Italy) tide gauge. J. Geod. 2019, 93, 1227–1244. [Google Scholar] [CrossRef]
  57. Dext3r-SIMC Platform. Available online: https://simc.arpae.it/dext3r/ (accessed on 19 April 2021).
  58. Arpae Website. Available online: www.arpae.it (accessed on 12 February 2021).
  59. Holgate, S.J.; Matthews, A.; Woodworth, P.L.; Rickards, L.J.; Tamisiea, M.E.; Bradshaw, E.; Foden, P.R.; Gordon, K.M.; Jevrejeva, S.; Pugh, J. New Data System and Products at the Permanent Service for Mean Sea Level. J. Coast. Res. 2013, 29, 493–504. [Google Scholar] [CrossRef]
  60. PSMSL Website. Available online: http://www.psmsl.org/data/obtaining/ (accessed on 3 January 2022).
  61. Romagnoli, C.; Zerbini, S.; Raicich, F.; Richter, B.; Lago, L.; Domenichini, D.; Thomsen, I. Sea level variations from tide gauges and GPS in the Northern Adriatic. In Proceedings of the EGS XXVII General Assembly, Nice, France, 21–26 April 2002. [Google Scholar]
  62. Hurrell, J.W.; van Loon, H. Decadal variations in climate associated with the north Atlantic oscillation. Clim. Chang. 1997, 36, 301–326. [Google Scholar] [CrossRef]
  63. Hurrell, J.W.; Kushnir, Y.; Ottersen, G.; Visbeck, M. An Overview of the North Atlantic Oscillation. In The North Atlantic Oscillation: Climatic Significance and Environmental Impact; Hurrell, J.W., Kushnir, Y., Ottersen, G., Visbeck, M., Eds.; American Geophysical Union: Washington, DC, USA, 2003. [Google Scholar] [CrossRef]
  64. NCAR Website. Available online: https://climatedataguide.ucar.edu/climate-data/hurrell-north-atlantic-oscillation-nao-index-station-based (accessed on 14 March 2021).
  65. Kerr, R.A. A North Atlantic Climate Pacemaker for the Centuries. Science 2000, 288, 1984–1985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. NOAA/PSD Website. Available online: https://psl.noaa.gov/data/timeseries/AMO/ (accessed on 14 March 2021).
  67. Universitat de Barcelona Website. Available online: http://www.ub.edu/gc/wemo/ (accessed on 29 June 2022).
  68. Lopez-Bustins, J.A.; Lemus-Canovas, M. The influence of the Western Mediterranean Oscillation upon the spatio-temporal variability of precipitation over Catalonia (northeastern of the Iberian Peninsula). Atmos. Res. 2020, 236, 104819. [Google Scholar] [CrossRef]
  69. Vicente-Serrano, S.M.; Begueria, S.; López-Moreno, J.I. A multiscalar drought index sensitive to global warming: The standardized precipitation evapotranspiration index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef]
  70. McKee, T.B.; Doesken, N.J.; Kleist, J. The relationship of drought frequency and duration to time scales. In Proceedings of the 8th Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993. [Google Scholar]
  71. Palmer, W.C. Meteorological Drought; U.S. Government Printing Office: Washington, DC, USA, 1965.
  72. CSIC Website. Available online: https://spei.csic.es/database.html (accessed on 18 April 2021).
  73. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  74. Wang, K.; Li, Z.; Cribb, M. Estimation of evaporative fraction from a combination of day and night land surface temperatures and NDVI: A new method to determine the Priestley–Taylor parameter. Remote Sens. Environ. 2006, 102, 293–305. [Google Scholar] [CrossRef]
  75. Vicente-Serrano, S.M.; Beguería, S.; Lorenzo-Lacruz, J.; Camarero, J.J.; López-Moreno, J.I.; Azorin-Molina, C.; Revuelto, J.; Morán-Tejeda, E.; Sanchez-Lorenzo, A. Performance of Drought Indices for Ecological, Agricultural, and Hydrological Applications. Earth Interact. 2012, 16, 1–27. [Google Scholar] [CrossRef]
  76. McEvoy, D.J.; Huntington, J.L.; Abatzoglou, J.T.; Edwards, L.M. An Evaluation of Multiscalar Drought Indices in Nevada and Eastern California. Earth Interact. 2012, 16, 1–18. [Google Scholar] [CrossRef]
  77. Zipper, S.C.; Qiu, J.; Kucharik, C.J. Drought effects on US maize and soybean production: Spatiotemporal patterns and historical changes. Environ. Res. Lett. 2016, 11, 094021. [Google Scholar] [CrossRef]
  78. Chen, T.; Xia, G.; Liu, T.; Chen, W.; Chi, D. Assessment of drought impact on main cereal crops using a standardized precipitation evapotranspiration index in Liaoning Province, China. Sustainability 2016, 8, 1069. [Google Scholar] [CrossRef]
  79. Zhang, Q.; Kong, D.; Singh, V.P.; Shi, P. Response of vegetation to different time-scales drought across China: Spatiotemporal patterns, causes and implications. Glob. Planet. Chang. 2017, 152, 1–11. [Google Scholar] [CrossRef]
  80. Páscoa, P.; Gouveia, C.M.; Russo, A.; Trigo, R.M. The role of drought on wheat yield interannual variability in the Iberian Peninsula from 1929 to 2012. Int. J. Biometeorol. 2017, 61, 439–451. [Google Scholar] [CrossRef]
  81. Labudová, L.; Labuda, M.; Takáč, A. Comparison of SPI and SPEI applicability for drought impact assessment on crop production in the Danubian lowland and the east Slovakian lowland. Theor. Appl. Climatol. 2017, 128, 491–506. [Google Scholar] [CrossRef]
  82. Muñoz Sabater, J. ERA5-Land Monthly Averaged Data from 1981 to Present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/10.24381/cds.68d2bb30?tab=overview (accessed on 10 July 2022).
  83. Muñoz Sabater, J. ERA5-Land Monthly Averaged Data from 1950 to 1980. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/10.24381/cds.68d2bb30?tab=overview (accessed on 10 July 2022).
  84. C3S Website. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land-monthly-means?tab=overview (accessed on 10 July 2022).
  85. Emilia-Romagna Region Website. Available online: https://geoportale.regione.emilia-romagna.it/ (accessed on 4 March 2021).
  86. Mann, H.B. Non-parametric tests against trend. Econometrica 1945, 13, 163–171. [Google Scholar] [CrossRef]
  87. Kendall, M.G. Rank Correlation Methods, 4th ed.; Charles Griffin: London, UK, 1975. [Google Scholar]
  88. Hamed, K.H.; Rao, A.R. A modified Mann-Kendall trend test for autocorrelated data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  89. Theil, H. A Rank-Invariant Method of Linear and Polynomial Regression Analysis; Nederl. Akad. Wetensch., Proc.; Springer: Berlin/Heidelberg, Germany, 1950. [Google Scholar]
  90. Sen, P.K. Estimates of the Regression Coefficient Based on Kendall’s Tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  91. Wilcox, R.R. Theil–Sen Estimator, Fundamentals of Modern Statistical Methods: Substantially Improving Power and Accuracy; Springer: New York, NY, USA, 2001. [Google Scholar]
  92. Zervas, C.E. Sea level variations of the United States, 1854–1999. In NOAA Technical Report NOS CO-OPS 36; U.S. Department of Commerce; National Oceanic and Atmospheric Administration; National Ocean Service: Silver Spring, MD, USA, 2011. [Google Scholar]
  93. Frich, P.; Alexander, L.V.; Della-Marta, P.; Gleason, B.; Haylock, M.; Klein Tank, A.M.G.; Peterson, T. Observed coherent changes in climatic extremes during the second half of the twentieth century. Clim. Res. 2002, 19, 193–212. [Google Scholar] [CrossRef]
  94. Tsimplis, M.N.; Raicich, F.; Fenoglio-Marc, L.; Shaw, A.G.P.; Marcos, M.; Somot, S.; Bergamasco, F. Recent developments in understanding sea level rise in the Adriatic coasts. Phys. Chem. Earth Parts A/B/C 2012, 40–41, 59–71. [Google Scholar] [CrossRef]
  95. Wöppelmann, G.; Marcos, M. Vertical land motion as a key to understanding sea level change and variability. Rev. Geophys. 2016, 54, 64–92. [Google Scholar] [CrossRef]
  96. Watson, P.J. An Assessment of the Utility of Satellite Altimetry and Tide Gauge Data (ALT-TG) as a Proxy for Estimating Vertical Land Motion. J. Coast. Res. 2019, 35, 1131–1144. [Google Scholar] [CrossRef]
  97. Oelsmann, J.; Passaro, M.; Dettmering, D.; Schwatke, C.; Sánchez, L.; Seitz, F. The zone of influence: Matching sea level variability from coastal altimetry and tide gauges for vertical land motion estimation. Ocean Sci. 2021, 17, 35–57. [Google Scholar] [CrossRef]
  98. Santamaría-Gómez, A.; Gravelle, M.; Collilieux, X.; Guichard, M.; Martín Míguez, B.; Tipahaneau, P.; Wöppelmann, G. Mitigating the effects of vertical land motion in tide gauge records using a state-of-the-art GPS velocity field. Glob. Planet. Chang. 2012, 98–99, 6–17. [Google Scholar] [CrossRef]
  99. Zerbini, S.; Plag, H.-P.; Baker, T.; Becker, M.; Billiris, H.; Bürki, B.; Kahle, H.-G.; Marson, I.; Pezzoli, L.; Richter, B.; et al. Sea level in the Mediterranean: A first step towards separating crustal movements and absolute sea-level variations. Glob. Planet. Chang. 1996, 14, 1–48. [Google Scholar] [CrossRef]
  100. Bouin, M.-N.; Wöppelmann, G. Land motion estimates from GPS at tide gauges: A geophysical evaluation. Geophys. J. Int. 2010, 180, 193–209. [Google Scholar] [CrossRef]
  101. Carbognin, L.; Teatini, P.; Tosi, L. The impact of relative sea level rise on the Northern Adriatic Sea coast, Italy. In Management of Natural Resources, Sustainable Development and Ecological Hazards II; Brebbia, C.A., Jovanovic, N., Tiezzi, E., Eds.; WIT Press: Southampton, UK, 2009. [Google Scholar]
  102. Furlani, S.; Biolchi, S.; Cucchi, F.; Antonioli, F.; Busetti, M.; Melis, R. Tectonic effects on Late Holocene sea level changes in the Gulf of Trieste (NE Adriatic Sea, Italy). Quat. Int. 2011, 232, 144–157. [Google Scholar] [CrossRef]
  103. Farolfi, G.; Del Ventisette, C. Contemporary crustal velocity field in Alpine Mediterranean area of Italy from new geodetic data. GPS Solut. 2016, 20, 715–722. [Google Scholar] [CrossRef]
  104. Vacchi, M.; Marriner, N.; Morhange, C.; Spada, G.; Fontana, A.; Rovere, A. Multiproxy assessment of Holocene relative sea-level changes in the western Mediterranean: Sea-level variability and improvements in the definition of the isostatic signal. Earth Sci. Rev. 2016, 155, 172–197. [Google Scholar] [CrossRef]
  105. Vilibić, I.; Šepić, J.; Pasarić, M.; Orlić, M. The Adriatic Sea: A long-standing laboratory for sea level studies. Pure Appl. Geophys. 2017, 174, 3765–3811. [Google Scholar] [CrossRef]
  106. Sánchez, L.; Völksen, C.; Sokolov, A.; Arenz, H.; Seitz, F. Presentday surface deformation of the Alpine region inferred from geodetic techniques, Earth Syst. Sci. Data 2018, 10, 1503–1526. [Google Scholar]
  107. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. A 1998, 45, 903–995. [Google Scholar] [CrossRef]
  108. Ezer, T.; Corlett, W.B. Is sea level rise accelerating in the Chesapeake Bay? A demonstration of a novel new approach for analyzing sea level data. Geophys. Res. Lett. 2012, 39, L19605. [Google Scholar] [CrossRef]
  109. Ezer, T.; Atkinson, L.P.; Corlett, W.B.; Blanco, J.L. Gulf Stream’s induced sea level rise and variability along the US mid-Atlantic coast. J. Geophys. Res. Oceans 2013, 118, 685–697. [Google Scholar] [CrossRef]
  110. Spada, G.; Galassi, G.; Olivieri, M. A study of the longest tide gauge sea-level record in Greenland (Nuuk/Godthab, 1958–2002). Glob. Planet. Chang. 2014, 118, 42–51. [Google Scholar] [CrossRef]
  111. Galassi, G.; Spada, G. Linear and non-linear sea-level variations in the Adriatic Sea from tide gauge records (1872–2012). Ann. Geophys. 2014, 57, P0658. [Google Scholar] [CrossRef]
  112. Bonaduce, A.; Pinardi, N.; Oddo, P.; Spada, G.; Larnicol, G. Sea-level variability in the Mediterranean Sea from altimetry and tide gauges. Clim. Dyn. 2016, 47, 2851–2866. [Google Scholar] [CrossRef]
  113. Vecchio, A.; Anzidei, M.; Serpelloni, E.; Florindo, F. Natural Variability and Vertical Land Motion Contributions in the Mediterranean Sea-Level Records over the Last Two Centuries and Projections for 2100. Water 2019, 11, 1480. [Google Scholar] [CrossRef]
  114. Meli, M.; Olivieri, M.; Romagnoli, C. Sea-Level Change along the Emilia-Romagna Coast from Tide Gauge and Satellite Altimetry. Remote Sens. 2021, 13, 97. [Google Scholar] [CrossRef]
  115. Croux, C.; Dehon, C. Influence function of the Spearman and Kendall correlation measures. Stat. Methods Appl. 2010, 19, 497–515. [Google Scholar] [CrossRef]
  116. Bartlett, M.S. Some aspects of the time-correlation problem in regard to tests of significance. J. R. Stat. Soc. 1935, 98, 536–543. [Google Scholar] [CrossRef]
  117. Santer, B.D.; Wigley, T.M.L.; Boyle, J.S.; Gaffen, D.J.; Hnilo, J.J.; Nychka, D.; Parker, D.E.; Taylor, K.E. Statistical significance of trends and trend differences in layer-average atmospheric temperature time series. J. Geophys. Res. 2000, 105, 7337–7356. [Google Scholar] [CrossRef]
  118. Lomb, N.R. Least-squares frequency analysis of unequally spaced data. Astrophys. Space Sci. 1976, 39, 447–462. [Google Scholar] [CrossRef]
  119. Scargle, J.D. Studies in astronomical time series analysis. II—Statistical aspects of spectral analysis of unevenly spaced data. Astrophys. J. 1982, 263, 835–853. [Google Scholar] [CrossRef]
  120. Baluev, R.V. Assessing the statistical significance of periodogram peaks. Mon. Not. R. Astron. Soc. 2008, 385, 1279–1285. [Google Scholar] [CrossRef]
  121. Oke, T.R. The energetic basis of the urban heat island. Q. J. R. Meteorol. Soc. 1982, 108, 1–24. [Google Scholar] [CrossRef]
  122. Swaid, H. Nocturnal variation of air-surface temperature gradients for typical urban and rural surfaces. Atmos. Environ. Part B Urban Atmos. 1991, 25, 333–341. [Google Scholar] [CrossRef]
  123. Paulo, A.A.; Rosa, R.D.; Pereira, L.S. Climate trends and behaviour of drought indices based on precipitation and evapotranspiration in Portugal. Nat. Haz. Earth Syst. Sci. 2012, 12, 1481–1491. [Google Scholar] [CrossRef]
  124. Gatto, P.; Carbognin, L. The Lagoon of Venice: Natural environmental trend and man-induced modification. Hydrol. Sci. J. 1981, 26, 379–391. [Google Scholar] [CrossRef]
  125. Tosi, L.; Teatini, P.; Strozzi, T.; Carbognin, L.; Brancolini, G.; Rizzetto, F. Ground surface dynamics in the northern Adriatic coastland over the last two decades. Rend. Fis. Acc. Lincei 2010, 21, 115–129. [Google Scholar] [CrossRef]
  126. Zerbini, S.; Bruni, S.; Raicich, F. Tide gauge data archaeology provides natural subsidence rates along the coasts of the Po Plain and of the Veneto–Friuli Plain, Italy. Geophys. J. Int. 2021, 225, 253–260. [Google Scholar] [CrossRef]
  127. Gambolati, G.; Teatini, P. Numerical analysis of land subsidence due to natural compaction of the upper Adriatic Sea Basin. In CENAS: Coastline Evolution of the Upper Adriatic Sea Due to Sea Level Rise and Natural and Anthropogenic Land Subsidence; Gambolati, G., Ed.; Springer: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  128. Carminati, E.; Doglioni, C.; Scrocca, D. Apennines subduction-related subsidence of Venice (Italy). Geophys. Res. Lett. 2003, 30, 1717. [Google Scholar] [CrossRef]
  129. Antonioli, F.; Ferranti, L.; Fontana, A.; Amorosi, A.; Bondesan, A.; Braitenberg, C.; Dutton, A.; Fontolan, G.; Furlani, S.; Lambeck, K.; et al. Holocene relative sea-level changes and vertical movements along the Italian and Istrian coastlines. Quat. Int. 2009, 206, 102–133. [Google Scholar] [CrossRef]
  130. Heymann, Y.; Steenmans, C.; Croisille, G.; Bossard, M. CORINE Land Cover. Technical Guide; Office for Official Publications of European Communities: Luxembourg, 1994. [Google Scholar]
  131. Foley, J.A.; Defries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef]
  132. Blondel, J. The ‘design’ of Mediterranean landscapes: A millennial story of humans and ecological systems during the historic period. Hum. Ecol. 2006, 34, 713–729. [Google Scholar] [CrossRef]
  133. Federal Office for the Environment (FOEN). Urban Sprawl in Europe; Publications Office of the European Union: Luxembourg, 2016. [Google Scholar]
  134. Scalenghe, R.; Marsan, F.A. The anthropogenic sealing of soils in urban areas. Landsc. Urban Plan. 2009, 90, 1–10. [Google Scholar] [CrossRef]
  135. Agenzia Regionale per la Protezione Ambientale (ARPA). Annuario Regionale dei Dati Ambientali; ARPA: Bologna, Italy, 2010. [Google Scholar]
  136. Todaro, V.; D’Oria, M.; Secci, D.; Zanini, A.; Tanda, M.G. Climate Change over the Mediterranean Region: Local Temperature and Precipitation Variations at Five Pilot Sites. Water 2022, 14, 2499. [Google Scholar] [CrossRef]
  137. Mariotti, A.; Zeng, N.; Yoon, J.; Artale, V.; Navarra, A.; Alpert, P.; Li, L.Z.X. Mediterranean water cycle changes: Transition to drier 21st century conditions in observations and CMIP3 simulations. Environ. Res. Lett. 2008, 3, 044001. [Google Scholar] [CrossRef]
  138. Maugeri, M.; Nanni, T. Surface Air Temperature Variations in Italy: Recent Trends and an Update to 1993. Theor. Appl. Climatol. 1998, 61, 191–196. [Google Scholar] [CrossRef]
  139. Brunetti, M.; Maugeri, M.; Nanni, T. Variations of Temperature and Precipitation in Italy from 1866 to 1995. Theor. Appl. Climatol. 2000, 65, 165–174. [Google Scholar] [CrossRef]
  140. Toreti, A.; Desiato, F. Temperature trend over Italy from 1961 to 2004. Theor. Appl. Climatol. 2008, 91, 51–58. [Google Scholar] [CrossRef]
  141. Balling, R.C., Jr.; Vose, R.S.; Weber, G.-R. Analysis of long-term European temperature records: 1751–1995. Clim. Res. 1998, 10, 193–200. [Google Scholar] [CrossRef]
  142. Klein Tank, A.M.G.; Wijngaard, J.B.; Koennen, G.P.; Boehm, R.; Demaree, G.; Gocheva, A.; Mileta, M.; Pashiardis, S.; Hejkrlik, L.; Kern-Hansen, C.; et al. Daily dataset of 20th-century surface air temperature and precipitation series for the European Climate Assessment. Int. J. Climatol. 2002, 22, 1441–1453. [Google Scholar] [CrossRef]
  143. Jones, P.D.; Parker, D.E.; Osborn, T.J.; Briffa, K.R. Global and Hemisphere Temperature Anomalies Land and Marine Instrumental Records; ESS-DIVE Repository; Carbon Dioxide Information Analysis Center (CDIAC); Oak Ridge National Laboratory (ORNL): Oak Ridge, TN, USA, 2005.
  144. Stanhill, G.; Cohen, S. Global dimming: A review of the evidence for a widespread and significant reduction in global radiation with discussion of its probable causes and possible agricultural consequences. Agric. For. Meteor. 2001, 107, 255–278. [Google Scholar] [CrossRef]
  145. Wild, M.; Gilgen, H.; Roesch, A.; Ohmura, A.; Long, C.N.; Dutton, E.G.; Forgan, B.; Kallis, A.; Russak, V.; Tsvetkov, A. From Dimming to Brightening: Decadal Changes in Solar Radiation at Earth’s Surface. Science 2005, 308, 847–850. [Google Scholar] [CrossRef] [PubMed]
  146. Kvalevåg, M.M.; Myrhe, G. Human impact on direct and diffuse solar radiation during the industrial era. J. Clim. 2007, 20, 4874–4883. [Google Scholar] [CrossRef] [Green Version]
  147. Kim, D.Y.; Ramanathan, V. Solar radiation budget and radiative forcing due to aerosols and clouds. J. Geophys. Res. 2008, 113, D02203. [Google Scholar] [CrossRef]
  148. Ruckstuhl, C.; Philipona, R.; Behrens, K.; Coen, M.C.; Dürr, B.; Heimo, A.; Mätzler, C.; Nyeki, S.; Ohmura, A.; Vuilleumier, L.; et al. Aerosol and cloud effects on solar brightening and the recent rapid warming. Geophys. Res. Lett. 2008, 35, L12708. [Google Scholar] [CrossRef]
  149. Streets, D.G.; Yan, F.; Chin, M.; Diehl, T.; Mahowald, N.; Schultz, M.; Wild, M.; Wu, Y.; Yu, C. Anthropogenic and natural contributions to regional trends in aerosol optical depth, 1980-2006. J. Geophys. Res. Atmos. 2009, 114, D00D18. [Google Scholar] [CrossRef]
  150. Wang, K.C.; Dickinson, R.E.; Liang, S.L. Clear sky visibility has decreased over land globally from 1973 to 2007. Science 2009, 323, 1468–1470. [Google Scholar] [CrossRef]
  151. Breshears, D.D.; Cobb, N.S.; Rich, P.M.; Price, K.P.; Allen, C.D.; Balice, R.G.; Romme, W.H.; Kastens, J.H.; Floyd, M.L.; Belnap, J.; et al. Regional vegetation die-off in response to global-change-type drought. Proc. Nat. Acad. Sci. USA 2005, 102, 15144–15148. [Google Scholar] [CrossRef]
  152. Adams, H.D.; Guardiola-Claramonte, M.; Barron-Gafford, G.A.; Villegas, J.C.; Breshears, D.D.; Zou, C.B.; Troch, P.A.; Huxman, E. Temperature sensitivity of drought-induced tree mortality portends increased regional die-off under global-change-type drought. Proc. Nat. Acad. Sci. USA 2009, 106, 7063–7066. [Google Scholar] [CrossRef]
  153. Tomozeiu. R.; Busuioc, A.; Marletto, V.; Zinoni, F.; Cacciamani, C. Detection of changes in the summer precipitation time series of the region Emilia-Romagna, Italy. Theor. Appl. Climatol. 2000, 67, 193–200. [Google Scholar] [CrossRef]
  154. Crespi, A.; Brunetti, M.; Lentini, G.; Maugeri, M. 1961–1990 high resolution monthly precipitation climatologies for Italy. Int. J. Climatol. 2018, 38, 878–895. [Google Scholar] [CrossRef]
  155. Ramanathan, V.; Crutzen, P.J.; Jiehl, J.T.; Rosenfeld, D. Aerosols, climate, and the hydrological cycle. Science 2001, 294, 2119–2124. [Google Scholar] [CrossRef] [PubMed]
  156. Wild, M.; Liepert, B. The Earth radiation balance as driver of the global hydrological cycle. Environ. Res. Lett. 2010, 5, 025003. [Google Scholar] [CrossRef]
  157. Wild, M. Enlightening Global Dimming and Brightening. Bull. Am. Meteor. Soc. 2012, 93, 27–37. [Google Scholar] [CrossRef]
  158. Wild, M. Global dimming and brightening: A review. J. Geophys. Res. 2009, 114, D00D16. [Google Scholar] [CrossRef]
  159. Mariotti, A.; Struglia, M.V.; Zeng, N.; Lau, K.-M. The Hydrological Cycle in the Mediterranean Region and Implications for the Water Budget of the Mediterranean Sea. J. Clim. 2002, 15, 1674–1690. [Google Scholar] [CrossRef]
  160. Seager, R.; Liu, H.; Henderson, N.; Simpson, I.; Kelley, C.; Shaw, T.; Kushnir, Y.; Ting, M. Causes of increasing aridification of the Mediterranean region in response to rising greenhouse gases. J. Clim. 2014, 27, 4655–4676. [Google Scholar] [CrossRef]
  161. Mariotti, A.; Pan, Y.; Zeng, N.; Alessandri, A. Long-term climate change in the Mediterranean region in the midst of decadal variability. Clim. Dyn. 2015, 44, 1437–1456. [Google Scholar] [CrossRef]
  162. Toreti, A.; Naveau, P.; Zampieri, M.; Schindler, A.; Scoccimarro, E.; Xoplaki, E.; Dijkstra, H.A.; Gualdi, S.; Luterbacher, J. Projections of global changes in precipitation extremes from Coupled Model Intercomparison Project Phase 5 Models. Geophys. Res. Lett. 2013, 40, 4887–4892. [Google Scholar] [CrossRef]
  163. Benini, L.; Bandini, V.; Marazza, D.; Contin, A. Assessment of land use changes through an indicator-based approach: A case study from the Lamone river basin in Northern Italy. Ecol. Indic. 2010, 10, 4–14. [Google Scholar] [CrossRef]
  164. Scanlon, B.R.; Jolly, I.; Sophocleous, M.; Zhang, L. Global impacts of conversion from natural to agricultural ecosystems on water resources: Quantity versus quality. Water Resour. Res. 2007, 43, W03437. [Google Scholar] [CrossRef]
  165. Filoso, S.; Bezerra, M.O.; Weiss, K.C.B.; Palmer, M.A. Impacts of forest restoration on water yield: A systematic review. PLoS ONE 2017, 12, e0183210. [Google Scholar] [CrossRef] [PubMed]
  166. Teuling, A.J.; de Badts, E.; Jansen, F.A.; Fuchs, R.; Buitink, J.; van Dijke, A.J.; Sterling, S. Climate change, re-/afforestation, and urbanisation impacts on evapotranspiration and streamflow in Europe. Hydrol. Earth Syst. Sci. 2019, 23, 3631–3652. [Google Scholar] [CrossRef] [Green Version]
  167. Milner, A.M.; Khamis, K.; Battin, T.J.; Brittain, J.E.; Barrand, N.E.; Füreder, L.; Cauvy-Fraunié, S.; Gíslason, G.M.; Jacobsen, D.; Hannah, D.M.; et al. Glacier shrinkage driving global changes in downstream systems. Proc. Nat. Acad. Sci. USA 2017, 114, 9770–9778. [Google Scholar] [CrossRef]
  168. López-Moreno, J.I.; García-Ruiz, J.M.; Beniston, M. Environmental change and water management in the Pyrenees. Facts and future perspectives for Mediterranean mountains. Glob. Planet. Chang. 2008, 66, 300–312. [Google Scholar] [CrossRef]
  169. García-Ruiz, J.M.; López-Moreno, J.I.; Vicente-Serrano, S.M.; Lasanta-Martínez, T.; Beguería, S. Mediterranean water resources in a global change scenario. Earth Sci. Rev. 2011, 105, 121–139. [Google Scholar] [CrossRef]
  170. Istituto Nazionale di Statistica (ISTAT), 2010. Censimento Agricoltura 2010. Database. Available online: http://dati-censimentoagricoltura.istat.it/ (accessed on 7 September 2021).
  171. European Commission. Guidelines on Best Practice to Limit, Mitigate or Compensate Soil Sealing; Publications Office of the European Union: Luxembourg, 2012. [Google Scholar]
  172. Kingston, D.G.; Lawler, D.M.; McGregor, G.R. Linkages between atmospheric circulation, climate and streamflow in the northern North Atlantic: Research prospects. Prog. Phys. Geogr. 2006, 30, 143–174. [Google Scholar] [CrossRef]
  173. Marullo, S.; Artale, V.; Santoleri, R. The SST multi-decadal variability in the Atlantic- Mediterranean region and its relation to AMO. J. Clim. 2011, 24, 4385–4401. [Google Scholar] [CrossRef]
  174. Mariotti, A.; Dell’Aquila, A. Decadal climate variability in the Mediterranean region: Roles of large-scale forcings and regional processes. Clim. Dyn. 2012, 38, 1129–1145. [Google Scholar] [CrossRef]
  175. Börgel, F.; Frauen, C.; Neumann, T.; Markus Meier, H.E. The Atlantic Multidecadal Oscillation controls the impact of the North Atlantic Oscillation on North European climate. Environ. Res. Lett. 2020, 15, 104025. [Google Scholar] [CrossRef]
  176. Tsimplis, M.N.; Josey, S.A. Forcing of the Mediterranean Sea by atmospheric oscillations over the North Atlantic. Geophys. Res. Lett. 2001, 28, 803–806. [Google Scholar] [CrossRef]
  177. Menemenlis, D.; Fukumori, J.; Lee, T. Atlantic to Mediterranean Sea Level Difference Driven by Winds near Gibraltar Strait. J. Phys. Oceanogr. 2007, 37, 359–376. [Google Scholar] [CrossRef]
  178. Tsimplis, M.N.; Calafat, F.M.; Marcos, M.; Jordà, G.; Gomis, D.; Fenoglio-Marc, L.; Struglia, M.V.; Josey, S.A.; Chambers, D.P. The effect of the NAO on sea level and on mass changes in the Mediterranean Sea. J. Geophys. Res. 2013, 118, 944–952. [Google Scholar] [CrossRef] [Green Version]
  179. Fenoglio-Marc, L.; Mariotti, A.; Sannino, G.; Meyssignac, B.; Carillo, A.; Struglia, M.V.; Rixen, M. Decadal variability of net water flux at the Mediterranean Sea Gibraltar Strait. Glob. Planet. Chang. 2013, 100, 1–10. [Google Scholar] [CrossRef]
  180. Della-Marta, P.M.; Luterbacher, J.; von Weissenfluh, H.; Xoplaki, E.; Brunet, M.; Wanner, H. Summer heat waves over western Europe 1880-2003, their relationship to large-scale forcings and predictability. Clim. Dyn. 2007, 29, 251–275. [Google Scholar] [CrossRef]
  181. López-Moreno, J.I.; Vicente-Serrano, S.M.; Morán-Tejeda, E.; Lorenzo-Lacruz, J.; Kenawy, A.; Beniston, M. Effects of the North Atlantic Oscillation (NAO) on combined temperature and precipitation winter modes in the Mediterranean mountains: Observed relationship and projections for the 21st century. Glob. Planet. Chang. 2011, 77, 62–76. [Google Scholar] [CrossRef]
  182. Rust, H.W.; Richling, A.; Bissolli, P.; Ulbrich, U. Linking teleconnection patterns to European temperature—a multiple linear regression model. Meteorol. Z. 2015, 24, 411–423. [Google Scholar] [CrossRef]
  183. Di Bacco, M.; Scorzini, A.R. Recent changes in temperature extremes across the north-eastern region of Italy and their relationship with large-scale circulation. Clim. Res. 2020, 81, 167–185. [Google Scholar] [CrossRef]
  184. Hurrell, J.W. Decadal trends in the North-Atlantic Oscillation: Regional temperatures and precipitation. Science 1995, 269, 676–679. [Google Scholar] [CrossRef]
  185. Xoplaki, E.; González-Rouco, J.F.; Luterbacher, J.; Wanner, H. Wet season Mediterranean precipitation variability: Influence of large-scale dynamics and trends. Clim. Dyn. 2004, 23, 63–78. [Google Scholar] [CrossRef]
  186. Serreze, M.C.; Carse, F.; Barry, R.G.; Rogers, J.C. Icelandic low cyclone activity: Climatological features, linkages with the NAO, and relationships with recent changes in the Northern Hemisphere circulation. J. Clim. 1997, 10, 453–464. [Google Scholar] [CrossRef]
  187. Cullen, H.M.; Kaplan, A.; Arkin, P.; Demenocal, P.B. Impact of the North Atlantic Oscillation on Middle Eastern climate and streamflow. Clim. Chang. 2002, 55, 315–338. [Google Scholar] [CrossRef]
  188. Struglia, M.V.; Mariotti, A.; Filograsso, A. River discharge into the Mediterranean Sea: Climatology and aspects of the observed variability. J. Clim. 2004, 17, 4740–4751. [Google Scholar] [CrossRef]
  189. Trigo, R.M.; Pozo-Vasquez, D.; Osborn, T.J.; Castro-Diez, Y.; Gamiz-Fortis, S.; Esteban-Parra, M.J. North Atlantic Oscillation influence on precipitation, river flow and water resources in the Iberian Peninsula. Int. J. Climatol. 2004, 24, 925–944. [Google Scholar] [CrossRef]
  190. Rimbu, N.; Boroneant, C.; Buta, C.; Dima, M. Decadal variability of the Danube River flow in the lower basin and its relation with the North Atlantic Oscillation. Int. J. Climatol. 2006, 22, 1169–1179. [Google Scholar] [CrossRef]
  191. López-Moreno, J.I.; Beguería, S.; Vicente-Serrano, S.M.; García-Ruiz, J.M. Influence of the North Atlantic Oscillation on water resources in central Iberia: Precipitation, streamflow anomalies, and reservoir management strategies. Water Resour. Res. 2007, 43, W09411. [Google Scholar] [CrossRef]
  192. Calafat, F.; Chambers, D.; Tsimplis, M. Mechanism of decadal sea level variability in the eastern North Atlantic and the Mediterranean Sea. J. Geophys. Res. Oceans 2012, 117, C9. [Google Scholar] [CrossRef]
  193. Mohamed, B.; Abdallah, A.M.; El-Din, K.A.; Nagy, H.; Shaltout, M. Inter-Annual Variability and Trends of Sea Level and Sea Surface Temperature in the Mediterranean Sea over the Last 25 years. Pure Appl. Geophys. 2019, 176, 3787–3810. [Google Scholar] [CrossRef]
  194. Landerer, F.W.; Volkov, D.L. The anatomy of recent large sea level fluctuations in the Mediterranean Sea. Geophys. Re. Lett. 2013, 40, 553–557. [Google Scholar] [CrossRef]
  195. Cazenave, A.; Cabanes, C.; Dominh, K.; Mangiarotti, S. Recent sea level changes in the Mediterranean Sea revealed by TOPEX/POSEIDON satellite altimetry. Geophys. Res. Lett. 2001, 28, 1607–1610. [Google Scholar] [CrossRef]
  196. Jones, P.D.; Osborn, T.J.; Briffa, K.R.; Folland, C.K.; Horton, E.B.; Alexander, L.V.; Parker, D.E.; Rayner, N.A. Adjusting for sampling density in grid box land and ocean surface temperature time series. J. Geophys. Res. 2001, 106, 3371–3380. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Location of the Reno and Lamone river basins in the Emilia-Romagna region, Northern Italy. TTG: Trieste tide gauge; (b) the coastal portion of the river basins. RTG: Marina di Ravenna tide gauge; (c) detailed map of the Reno and Lamone river basins and tributaries, with position of gauges and stations used in this study.
Figure 1. (a) Location of the Reno and Lamone river basins in the Emilia-Romagna region, Northern Italy. TTG: Trieste tide gauge; (b) the coastal portion of the river basins. RTG: Marina di Ravenna tide gauge; (c) detailed map of the Reno and Lamone river basins and tributaries, with position of gauges and stations used in this study.
Water 14 02650 g001
Figure 2. Duration of the considered records.
Figure 2. Duration of the considered records.
Water 14 02650 g002
Figure 3. (a) Reno River discharge at the Casalecchio river gauge and (b) the stacked rainfall time series, both with 1-year low-pass filter time series (solid black and blue lines, respectively) and trend estimates (dashed red lines); (c) the stacked temperature time series (black), its trend estimate (dashed black line), and the 1-year moving average time series of max (red), mean (yellow) and min (blue) temperatures observed at the nine stations in this study; (d) the number of wet days filtered with a 4-year moving window (solid dark-blue line) and the SDII, showing only instances above the 95th (light blue) and 99th (red) percentiles; (e) the 48-month SPEI time series for the study site.
Figure 3. (a) Reno River discharge at the Casalecchio river gauge and (b) the stacked rainfall time series, both with 1-year low-pass filter time series (solid black and blue lines, respectively) and trend estimates (dashed red lines); (c) the stacked temperature time series (black), its trend estimate (dashed black line), and the 1-year moving average time series of max (red), mean (yellow) and min (blue) temperatures observed at the nine stations in this study; (d) the number of wet days filtered with a 4-year moving window (solid dark-blue line) and the SDII, showing only instances above the 95th (light blue) and 99th (red) percentiles; (e) the 48-month SPEI time series for the study site.
Water 14 02650 g003
Figure 4. (a) RSL recorded at the Marina di Ravenna tide gauge (RTG, 1873–2016) and related 1-year moving window time series (yellow line); (b) the TTG, 1875–2021), related trend estimate (dashed red line), and 1-year moving average time series (yellow line); (c) the difference between the TTG and RTG time series with the best fit (solid red) line; (d) the RTG time series corrected for vertical land movements, long-term trend estimate (red dashed line), and 1-year moving average curve (yellow line).
Figure 4. (a) RSL recorded at the Marina di Ravenna tide gauge (RTG, 1873–2016) and related 1-year moving window time series (yellow line); (b) the TTG, 1875–2021), related trend estimate (dashed red line), and 1-year moving average time series (yellow line); (c) the difference between the TTG and RTG time series with the best fit (solid red) line; (d) the RTG time series corrected for vertical land movements, long-term trend estimate (red dashed line), and 1-year moving average curve (yellow line).
Water 14 02650 g004
Figure 5. The IMFs and RES curves for river discharge (a), rainfall (b), temperature (c) and sea level (d), from the EMD analysis; the winter NAO (e) and AMO (f) indices time series, with a 10-year filter.
Figure 5. The IMFs and RES curves for river discharge (a), rainfall (b), temperature (c) and sea level (d), from the EMD analysis; the winter NAO (e) and AMO (f) indices time series, with a 10-year filter.
Water 14 02650 g005
Figure 6. Land use changes (in hectares) within the Reno and Lamone river basins since 1976, for the Apennines (left panel) and alluvial plain (right panel).
Figure 6. Land use changes (in hectares) within the Reno and Lamone river basins since 1976, for the Apennines (left panel) and alluvial plain (right panel).
Water 14 02650 g006
Table 1. Datasets considered in this study and related residual trends.
Table 1. Datasets considered in this study and related residual trends.
StationVariableResidual
Rate (*·Year−1)
Period
Casalecchio di Reno
(Reno River)
River
discharge
(m3/s)
−0.202 ± 0.0021921–2021
−0.101 ± 0.0041921–1950
−0.190 ± 0.0011950–1980
−0.268 ± 0.0021993–2021
Stacked rainfallRainfall(mm)−0.101 ± 0.0051921–2021
0.157 ± 0.0041921–1950
−0.263 ± 0.0051950–2000
0.200 ± 0.0062000–2021
Stacked temperatureTemperature
(°C)
0.032 ± 0.0021925–2021
0.040 ± 0.0011925–1950
−-0.007 ± 0.0011950–1980
0.071 ± 0.0011980–2021
Marina di Ravenna
(VLM corrected)
Sea level
(mm)
0.786 ± 0.0081875–2016
Note: * symbol: unit of measurement of the specific variable.
Table 2. Most evident periodicities that affect the variables considered in this study, and their related Standard Normalized Power (Snp, according to Baluev [120]). Indications in parentheses for some periodicities highlight their dominance within a specific IMF.
Table 2. Most evident periodicities that affect the variables considered in this study, and their related Standard Normalized Power (Snp, according to Baluev [120]). Indications in parentheses for some periodicities highlight their dominance within a specific IMF.
VariablePeriodSnp
NAO4–50.041
~60.042
7–80.087
8–90.023
~120.012
~140.022
16–170.011
22–230.04
37–400.04
AMO4–50.015
~60.018
7–80.025
8–90.05
10–110.047
16–170.018
20–220.012
~670.5
River discharge4–50.056
~6 (IMF 4)0.073
8–90.037
10–110.072
~120.068
~14 (IMF 5)0.058
16–170.028
22–23 (IMF 6)0.056
33–350.052
53–560.097
Rainfall4–5 (IMF 4)0.094
~60.038
8–90.035
10–11 (IMF 5)0.052
~120.027
16–170.019
23–260.021
38–41 (IMF 6)0.061
Temperature4–5 (IMF 4)0.04
8–90.03
~120.042
~14 (IMF 5)0.06
16–170.02
22–23 (IMF 6)0.07
30–330.022
53–560.19
Sea level4–50.039
8–9 (IMF 4)0.043
10–110.02
~120.038
~140.076
16–17 (IMF 5)0.079
22–230.06
28–300.13
40–42 (IMF 6)0.29
~670.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Meli, M.; Romagnoli, C. Evidence and Implications of Hydrological and Climatic Change in the Reno and Lamone River Basins and Related Coastal Areas (Emilia-Romagna, Northern Italy) over the Last Century. Water 2022, 14, 2650. https://doi.org/10.3390/w14172650

AMA Style

Meli M, Romagnoli C. Evidence and Implications of Hydrological and Climatic Change in the Reno and Lamone River Basins and Related Coastal Areas (Emilia-Romagna, Northern Italy) over the Last Century. Water. 2022; 14(17):2650. https://doi.org/10.3390/w14172650

Chicago/Turabian Style

Meli, Matteo, and Claudia Romagnoli. 2022. "Evidence and Implications of Hydrological and Climatic Change in the Reno and Lamone River Basins and Related Coastal Areas (Emilia-Romagna, Northern Italy) over the Last Century" Water 14, no. 17: 2650. https://doi.org/10.3390/w14172650

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop