The city of Venice and the surrounding lagoonal ecosystem are highly vulnerable to variations in relative sea level. In the past ∼150 years, this was characterized by an average rate of relative sea-level rise of about 2.5 mm/year resulting from the combined contributions of vertical land movement and sea-level rise. This literature review reassesses and synthesizes the progress achieved in quantification, understanding and prediction of the individual contributions to local relative sea level, with a focus on the most recent studies. Subsidence contributed to about half of the historical relative sea-level rise in Venice. The current best estimate of the average rate of sea-level rise during the observational period from 1872 to 2019 based on tide-gauge data after removal of subsidence effects is 1.23 ± 0.13 mm/year. A higher – but more uncertain – rate of sea-level rise is observed for more recent years. Between 1993 and 2019, an average change of about +2.76 ± 1.75 mm/year is estimated from tide-gauge data after removal of subsidence. Unfortunately, satellite altimetry does not provide reliable sea-level data within the Venice Lagoon. Local sea-level changes in Venice closely depend on sea-level variations in the Adriatic Sea, which in turn are linked to sea-level variations in the Mediterranean Sea. Water mass exchange through the Strait of Gibraltar and its drivers currently constitute a source of substantial uncertainty for estimating future deviations of the Mediterranean mean sea-level trend from the global-mean value. Regional atmospheric and oceanic processes will likely contribute significant interannual and interdecadal future variability in Venetian sea level with a magnitude comparable to that observed in the past. On the basis of regional projections of sea-level rise and an understanding of the local and regional processes affecting relative sea-level trends in Venice, the likely range of atmospherically corrected relative sea-level rise in Venice by 2100 ranges between 32 and 62 cm for the RCP2.6 scenario and between 58 and 110 cm for the RCP8.5 scenario, respectively. A plausible but unlikely high-end scenario linked to strong ice-sheet melting yields about 180 cm of relative sea-level rise in Venice by 2100. Projections of human-induced vertical land motions are currently not available, but historical evidence demonstrates that they have the potential to produce a significant contribution to the relative sea-level rise in Venice, exacerbating the hazard posed by climatically induced sea-level changes.
Introduction
This paper reviews the current knowledge about mean relative sea-level (RSL) changes in the Venice Lagoon on interannual to centennial timescales and the associated contribution from oceanic, land and atmospheric processes. The assessment includes a paleoperspective, considering the Quaternary period. It encompasses an overview of available observed estimates of historical RSL changes in Venice (Sect. ) and quantification of the individual contributions by the major underlying processes, including vertical land motions (Sect. ) and climatic changes (Sect. ). Estimates are supported by a review of downscaling mechanisms of global and large-scale oceanic and atmospheric signals to the Venice Lagoon (Sect. ), with a special focus on processes in the Atlantic and Euro-Mediterranean regions. Estimates of projected long-term future RSL changes based on state-of-the-art models of vertical land motions and of sea-level rise under different scenarios of anthropogenic greenhouse gas emission are discussed, with emphasis on the associated major sources of uncertainty (Sect. ). The review primarily focuses on papers published in the past decade and also aims to define the overarching open research questions and possible approaches for progress (Sect. ).
Given the multidisciplinarity of this review, it is useful to specify the meaning of terms and concepts associated with sea-level changes that are recurrent in this paper and often used inconsistently by different scientific communities (see also , for a broader discussion). Unless otherwise specified in the text, the following definitions apply.
Relative sea-level (RSL) change: change in local sea level relative to the local solid surface , used for tide-gauge data;
Geocentric sea-level (GSL; used for satellite altimetry sea-level data) change: change in local sea level with respect to a geocentric reference, namely a Terrestrial Reference Frame or, equivalently, a reference ellipsoid ;
Subsidence: land surface sinking (; see also );
Vertical land movement (VLM)-corrected RSL: local sea level derived from tide-gauge RSL data corrected for vertical land movements;
Global-mean sea level (GMSL): spatially averaged RSL over the world ocean.
Note that satellite GSL data measure the open sea and therefore do not capture coastal variations. In addition, they are insensitive to the vertical-land-movement component of RSL. Therefore, GSL is not the same as VLM-corrected RSL. A list of abbreviations recurrently used in the paper is provided in Table 1.
List of abbreviations used in the paper.
BPBefore presentOOrder of magnitudeNAONorth Atlantic OscillationRSLRelative sea levelEAWREastern Atlantic–western Russia patternSCAScandinavian patternGIAGlacial isostatic adjustmentGSLGeocentric sea levelCMComune Marino or Comune Alta MareaMTLMean tide levelGMSLGlobal-mean sea levelZMPSZero Mareografico Punta SaluteALESAdaptive leading-edge subwaveformRMSRoot mean squareRADSRadar Altimeter Database SystemESA – SLCCIEuropean Space Agency – Sea Level Climate Change InitiativeCTOHCentre of Topography of the Oceans and the HydrosphereSHYFEMShallow water HYdrodynamic Finite Element ModelAICAkaike information criterionFPEFinal prediction error
The reader is referred to and in this special issue for details about the geographical and historical setting of the Venice Lagoon, the linkage between RSL changes and the phenomenology of surges and extreme water levels affecting the lagoon, about their prediction, and about broader implications for the ecosystems and the historical city.
Monitoring sea-level changes
The monitoring of sea-level changes in Venice relies on both in situ data acquired by tide gauges (Sect. ) and remote sensing observations provided by satellite radar altimetry (Sect. ). Tide gauges record sea-level heights with reference to a permanent benchmark on land. Therefore, they provide measurements of RSL embedding the effects of vertical land motion (Sect. 3). Tide-gauge datasets consist of local, long-term measurements acquired at high frequency and accuracy . Satellite radar altimetry, in contrast, measures GSL changes . These are independent of variations in the local land level changes and hence unaffected by the potentially key component of RSL . These measurements have a lower sampling rate (several days) and a lower accuracy than those provided by tide gauges, but they are representative of wider oceanic areas and have the potential to characterize the evolution of sea-level variability with an almost global coverage.
Tide gauges
Tide gauges have been providing sea-level data in Venice for about 150 years. Historically, the establishment of tide gauges was primarily dictated by navigational needs and tidal measurements, with an operational accuracy of a few centimeters. The first self-recording tide gauge in Venice was installed at Palazzo Loredan, Campo Santo Stefano, in Rio San Vidal, at a distance of about 100 m from the Grand Canal (Fig. ). Systematic measurements began on 27 November 1871. The observations were performed under the responsibility of the Civil Engineering Office (Ufficio del Genio Civile) until 27 July 1896, when the management was taken over by the Italian Military Geographic Institute (Istituto Geografico Militare), which was also in charge of land leveling. Two additional tide gauges became operational in 1888 and 1906. The first one, owned and managed by the Royal Italian Navy, was installed in the Venice Arsenal; the second one was installed in the Grand Canal, near Punta della Salute. The tide gauges at Santo Stefano and in the Arsenal were decommissioned in 1911 and 1917, respectively. In 1923, the gauge on the Grand Canal was moved to the Giudecca Canal side of Punta della Salute. This gauge is still active under the management of the “Istituto Superiore per la Protezione e la Ricerca Ambientale” (ISPRA; Venice branch, https://www.venezia.isprambiente.it/, last access: 5 August 2021). Since 2002, a gauge on the Grand Canal side has again been operational on the site of the previous installation thanks to the Venice municipality (“Centro Previsione e Segnalazione Maree”; https://www.comune.venezia.it/it/content/centro-previsioni-e-segnalazioni-maree, last access: 5 August 2021).
Further details on the tide gauges installed in the Venice Lagoon up to the early 20th century are reported by . reviews the sea-level observations in Venice and summarizes the main development stages of the observational network in the Venice Lagoon, including lists of active and decommissioned tide-gauge stations. provide the most recent review of tide-gauge data for Venice and collect quality-checked published and unpublished records of high and low waters since 1872.
Linking the data from the various tide gauges to provide one continuous dataset of long-term sea-level change requires an accurate knowledge of the corresponding reference levels (or datums) on land. Before the 1910s, the most common vertical reference level in Venice was the so-called “comune alta marea” or “comune marino” (CM). The CM represents the upper edge of the green belt formed by algae on quays and walls and corresponds to mean high water. It was often indicated by an engraved horizontal segment and/or a “C” . According to , Mati established the tide-gauge datum at Santo Stefano at 1.50 m below the CM of 1825. In 1910, a new datum was adopted, namely the mean tide level (MTL) of 1884–1909 (central year 1897), computed from the high and low waters measured at Santo Stefano. According to , it turned out to be 1.2754 m above the tide-gauge datum and 0.2246 m below the CM of the year 1825. This new reference was named the “Zero Mareografico Punta Salute” (ZMPS). Since 1910, the ZMPS has been the standard reference for RSL observations in Venice. The benchmarks of the two tide gauges at Punta della Salute were also connected to the leveling network in 1910 and 1923, respectively. The heights of the various benchmarks and vertical reference levels are shown in the inset of Fig. . The record of high and low waters since 1872 allowed a composite 148-year RSL time series to be developed from 1872 to 2019, with very few gaps (Fig. ). Note that these data consist of MTL which, in principle, differ from the MSL computed from, for example, hourly data because MTL does not account for shallow-water tidal effects. However, the difference is negligible in Venice. In fact, from observations covering the period 1940–2012, obtained MTL-MSL=-0.1±0.1 cm, and the Permanent Service for Mean Sea Level provides a value of 0.0 cm, estimated according to (https://www.psmsl.org/data/obtaining/stations/168.php, last access: 5 August 2021).
Estimates of centennial rates of Venetian RSL and VLM-corrected RSL rise based on tide-gauge data are summarized in Sect. .
Time series of annual-mean RSL measured by tide gauges in Venice and Trieste. Venice data are relative to ZMPS; Trieste data are offset for illustrative purposes. The top left inset defines the reference planes of the tide gauges at Santo Stefano and Punta della Salute (redrawn from ).
Altimetric data
Since the first satellite altimetry missions in the mid-1970s, the accuracy of sea-surface height estimates has increased considerably until high-precision and routinely measured altimetric data were made available in the early 1990s with the launch of the TOPEX/Poseidon mission. Since then, satellite radar altimeters have been providing an operational global monitoring of the GSL . The spatial resolution of these data is controlled by the orbital parameters selected for each mission as the radar altimeters acquire narrow threads of measurements along those portions of the ocean surface that are directly overflown by the satellite (along-track data). Depending on the orbital period, these tracks can be separated by hundreds of kilometers, limiting the actual spatial coverage provided. For example, the Jason-3 mission that continues the climatic sea-level record started in 1993 with TOPEX/Poseidon has only four tracks crossing the Adriatic Sea and does not cover the Venice area. To improve the monitoring, it is possible to take advantage of the data collected by various radar altimeters flying at the same time. However, this requires a characterization of the inter-mission biases and the development of suitable interpolation schemes of the independent ground tracks. Multimission datasets are typically distributed over regular grids.
An additional potential limiting factor for the Venice area is the degradation of the technique performance towards the coast resulting from the contaminating presence of land in the satellite footprint and from the enhanced inhomogeneity of the local ocean surface. Limitations and possible perspectives of coastal altimetry in the Adriatic Sea have been discussed in several studies since the late 1990s . This has motivated further investigations based on the latest coastal altimetry datasets (e.g., CTOH; see ) and/or reprocessing initiatives (e.g., ALES; see ). In the northern Adriatic, these studies analyzed data around Venice and Trieste, including their validation against tide-gauge measurements. The results show that a reasonable increase in quantity and quality of data can be achieved compared to standard products up to a few kilometers from the coastline. Figure illustrates the example of the Gulf of Trieste, where three missions cross the area, and a data gap exists with standard products. In this case, the number of outliers along the Jason-1 and Jason-2 tracks is almost always less than the standard product, and the improvement is clear until 6 km from the coast. The comparative assessments with tide gauges confirm that the correlation of coastal altimetry products is always higher than that of standard products and that the difference in sea-level estimates provided by the two techniques is typically below 10 cm in proximity of the point of closest approach to the tide gauge. Among the most relevant reprocessing efforts of the last years, we should mention the Sea Level Climate Change Initiative (SLCCI) of the European Space Agency, which encompassed nine satellite radar altimetry missions over the period 1993–2015 . The SLCCI product, distributed over a homogenous grid of 0.25∘, contains data close to the coast, e.g., 10 km near Trieste, and was used for the assessment of coastal sea-level trends .
Percentage of outliers along-track for Envisat track 416 (top panel), Jason-1 track 161 (middle panel) and Jason-2 track 196 (bottom panel) from the Sensory Geophysical Data Record (SGDR) (blue line) and ALES (red line) products. Land is shaded in gray. The green line represents the distance from the closest coast (adapted from ).
The sufficient maturity of the algorithms and processing in coastal altimetry offered the opportunity to extend the SLCCI product to the coastal zone. During the bridging phase in 2018 a new product with an along-track spacing of about 350 m for estimating GSL trends was developed in selected regions, including the Adriatic Sea. The experimental dataset only covered the period from July 2002 to June 2016 and the Jason-1 and Jason-2 missions. It combines the post-processing strategy of X-TRACK and the advantage of the ALES re-tracker . The product was tested along track 196 in the Gulf of Trieste. The improvement is particularly good in the entire Gulf of Trieste , confirming what was found by .
Altimetry-based assessments of multidecadal trends of Venetian GSL are summarized in Sect. .
Vertical land movement
Vertical land movement is a critical component of Venetian RSL variability. Therefore, the phenomena that control the vertical land movement are presented in the following together with their relevant timescales and estimated trends. This section includes a paleoperspective on vertical land movements and considers processes whose characteristic timescales extend in some cases largely beyond the observational RSL period. We consider information on such timescales essential in the context of this review to understand ongoing processes and frame them within the correct timescale. The aim is therefore to provide the reader with an overview of the main characteristics of the local vertical land movement, of the methods that allow quantifying it over different time intervals and of the resulting uncertainties. The joint consideration of all these elements determines the constraints on our current ability to make predictions of the future evolution of the local vertical land movement.
The vertical velocity of a given area results from the sum of different velocity components due to tectonics, sediment loading, sediment compaction, glacial isostatic adjustment (GIA) and anthropic activities .
In the Venice area, all the components listed above induce non-negligible displacements even though their magnitude and relative importance have changed over time. The net result is a time-dependent land lowering (subsidence) that enhances RSL. Natural and anthropogenic components are assumed to act on different timescales: millions to thousands of years and hundreds to tens of years, respectively. This assumption allows a separation of the factors controlling sea-level changes if the estimates of vertical land movements over different time spans are available .
Natural land movements
The Venice area is naturally subsiding. This process is characterized by a long-term component controlled by tectonics and geodynamics and sedimentation, active on time spans of about 106–104 years, and a short-term component controlled by glaciation cycles and due to GIA processes acting on periods of 103–104 years .
Depending on the time interval considered, different datasets are available for investigating the rate of vertical land movement. Subsidence rates up to 2 Myr ago can be inferred from the thickness of the different layers of Quaternary sediments. Over this time frame, sedimentation rates are equivalent to subsidence rates since the entire sedimentary sequence was deposited in shallow marine to continental environments . Investigation further in the past is made through seismic lines, which indicate buried interfaces between materials of different acoustic impedances, and drilled cores. Deposition rates can be computed using sedimentological indicators (, and references therein; ; ), nannofossil biostratigraphy, paleomagnetic polarity and magnetic susceptibility . Additional techniques are available for more recent epochs. Radiocarbon dating allows the investigation of organic sediments, mainly peats, up to ∼50 000 years ago , while the depth of archeological remains and historical data provide information for the last few thousand years . Finally, information on the natural component of the contemporary land subsidence is provided by tide-gauge and leveling measurements made before the 1930s, when human activities impacting land subsidence started to develop . The following sections illustrate the evolution of the natural component of subsidence using the marine isotope stage (MIS) 5.5 event (∼130–120 kyr ago) as a reference to separate geologically older and newer RSL changes. Due to its relevance within geophysical studies on sea-level variations, a dedicated section on GIA is also provided.
Pleistocene up to MIS 5.5
The natural subsidence of Venice on timescales from tens of millennia to millions of years is controlled by sedimentary and tectonic and geodynamic processes. Venice is located at the northeastern border of the Po Plain (Figs. and ), which is the foreland basin of two fold-and-thrust belts: the north-northeastern vergent northern Apennines and the southern vergent southern Alps . Figure shows the geometry of the foreland regional monocline related to the subduction of the Adriatic plate (that includes the Po Plain) below the northern Apennines, from the southern Po Plain to the Friuli region, as reconstructed from seismic reflection profiles. The dip of the regional monocline gradually decreases from about 22∘ to close to 0∘. This geometry is consistent with the southward-increasing thickness of Quaternary sediments found in borehole stratigraphies . These data imply that the long-term component of subsidence in the Po Plain and in Venice is almost entirely controlled by the retreat and flexure of the Adriatic plate subducting underneath the Apennines .
Geological setting of the Po Plain area with dominant tectonic features (adapted from ).
derive a more complex evolution for the Venice area from integrated magneto-bio-cyclo-stratigraphy and palynofloral analyses on the VENEZIA-1 borehole, drilled in 1971 by the Consiglio Nazionale delle Ricerche down to a total depth of 950 m. They concluded that the region collapsed about 1.8 Myr ago and was characterized by slow marine sediment accumulation until around 0.8 Myr ago, shoaling rapidly in subsequent times. The initial transition to continental sediments occurred during a glacioeustatic low stand dated at 0.43 or 0.25 Ma. From the VENEZIA-1 record, calculated a total long-term subsidence rate of less than 0.5 mm/year, about half of that proposed earlier on less refined data, and a mean subsidence rate of 0.36 mm/year for the last 600 kyr. The latter value is considerably lower when compared with rates obtained for the Holocene and the upper MIS 5.5. The most reasonable interpretation is that the Mid-Pleistocene rates are unavoidably averaged over many cycles of quiescence and rapid motion; thus they cannot be readily compared to shorter periods, which could experience phases of rapid change induced by both natural and anthropogenic factors.
Concerning natural variations acting on shorter timescales (103–104 years), several transgressive–regressive Pleistocene cycles are recorded in well-core stratigraphies consisting of alternating shallow marine and continental deposits. In the Late Quaternary, the evolution of the Venetian–Friulian Plain was strongly influenced by glacial cycles, and a general regressive trend is apparent . The coastal to shallow marine deposits assigned to MIS 5.5 can be tracked in borehole logs up to 30 km west of the present shoreline. South of the Po Delta, the base of the Tyrrhenian coastal deposits lies at about 125 m below sea level, but its depth rapidly increases toward the south along the Romagna coastal plain . This pattern may reflect the northeasterly retreat of the Adriatic slab .
The MIS 5.5 markers allow the calculation of reliable rates because compaction is negligible, the basal MIS 5.5 unconformity is widely distributed, and the overlying lagoonal paralic sediments in cores are fairly easy to recognize. Several sites related to sea-level position during MIS 5.5 are considered in and . These have a fairly good W–E distribution along the distal sector of the Venetian Plain. The stratigraphic data were obtained from boreholes mainly drilled for the geological map of Italy (Progetto CARG – Cartografia Geologica e Geotematica, Veneto region) and for the mobile-barrier-based protective system (so-called “MOdulo Sperimentale Elettromeccanico”, or MOSE; see ) project by the Venice Water Authority. An error bar of ±2 m was assigned because the sediments are lagoonal. The northwestern Adriatic coast (Friulian–Venetian Plain) shows homogeneous subsidence, with rates ranging between 0.58 and 0.69 mm/year. The MIS 5.5 data from the VENEZIA-1 core provide a rate of 0.69 mm/year .
Late Pleistocene and Holocene
After the Last Glacial Maximum several lagoons developed along the Adriatic Sea, formed by the rivers flowing into the sea. Only two of them, the Grado and Venice lagoons, still exist today, while the rest have been infilled by sedimentation . Recent stratigraphic information about Holocene sea levels (2–6 cal kyr BP) was obtained from lagoonal deposits found in boreholes between the Tagliamento River and the city of Venice. Other data were derived from archeological markers reported in the abundant literature available for the Venice Lagoon and its mainland . The shell base of the lagoon indicates a subsidence rate over the last 7.3 kyr of 1.6 ± 0.3 mm/year (average of ratios of values in columns H and G for sites 26 and 28 in Table 1 of ). The higher Holocene subsidence with respect to the MIS 5.5 is possibly due to sediment compaction, which does not contribute to the long-term rate .
Subsidence rates up to 1.2–1.3 mm/year were calculated by radiocarbon dating on Late Pleistocene and Holocene deposits of the Venice Lagoon e.g., . This estimate is interpreted as the average of a time-varying trend related to periods of excess sedimentation alternating with periods without deposition or even with erosion . Indeed, over relatively short periods, different rates can be observed. For instance, the largest rate, ∼5 mm/year, occurred during the Last Glacial Maximum, which induced the maximum effect of isostatic lowering.
Finally, it is generally assumed that natural subsidence of Venice is continuous in time. However, abrupt catastrophic pulses of subsidence cannot be ruled out as suggested by the sudden disappearance of the island of Malamocco at the beginning of the 12th century. investigated the potential effects of earthquakes on the subsidence of Venice by means of numerical models. The authors concluded that, while the coseismic effects of a single event are unlikely to be detectable, a priori they cannot be considered to be negligible given the number of seismogenic sources within a 100 km distance from the town. These authors, however, suggest that earthquake-induced liquefaction may cause or have caused local acceleration of subsidence in Venice. For example, the destruction and sinking of ancient Malamocco is roughly coincident with a strong earthquake cycle that was associated with phenomena possibly explained by liquefaction of sandy layers.
A summary of the rates of natural subsidence discussed in this and in the previous sections is presented in Table . The values reported in the literature are often presented without indicating the corresponding uncertainty level. In some cases, it is even explicitly stated that the data available to the study did not allow for a quantification of uncertainty e.g., . Uncertainty estimation is further complicated by the fact that subsidence varies not only with time but also in space depending on the local conditions of the subsoil . For what concerns uncertainties in geomorphological and historical markers, proposed a strategy based on archeological metadata and on standard bathymetric corrections for the Holocene and Late Pleistocene . The resulting median uncertainty for the Venice area is 0.2 mm/year (markers 24–30 in Table 1 of ).
Time evolution of the natural component of land subsidence in the Venetian region over geological timescales.
PeriodSubsidence rate [mm/yr]Data sourceReference(s)Last 2 Myr∼0.5Nannofossil biostratigraphy, paleomagnetic polarity, magnetic susceptibility and sedimentological facies of a drilled coreLast 1.43 Myr0.7–1.0Thickness of Pleistocene sediments from seismic lines and boreholesLast 125 kyr0.58–0.69MIS 5.5 paralic deposits in drilled coresLast 40 kyr1.2–1.3Radiocarbon dating on organic remains, mainly peats and shells1Same as previous line4–5 kyr1.1±0.3Geomorphological and archeological markers (average of ratios of values in columns H and G for sites 17, 27 and 30 in their Table 1). Original data from , and Glacial isostatic adjustment
GIA describes the response of the solid Earth to the growth and decay of continental ice sheets because of climate variations (for recent reviews, see, and ). GIA stems from interactions between the cryosphere, the solid Earth and the oceans, involving sluggish deformations of the crust driven by surface mass redistribution, mutual gravitational attraction and rotational variations . The GIA-induced RSL variations are characterized by a strong regional imprint reflecting such interactions. They can be modeled by means of the sea-level equation first introduced by , which is an implicit equation that accounts for variations in the Earth’s topography in response to sea-level change, consistently with changes in the gravity field . Among the processes contributing to present-day RSL change e.g., , GIA is the only one that is sensitive to the solid Earth rheology. Because changes in the Earth system observed by geodetic methods would be unfeasible without taking GIA properly into account (e.g., ), GIA modeling plays an important role in the study of the impacts of contemporary and future climate change.
Due to the widespread evidence of past RSL variations since the Late Holocene across the Mediterranean Sea, much work has been done to reconcile field observations of past RSL variations with GIA modeling predictions (, and references therein). In two recent contributions, attention has been paid to the history of sea level in the northern Adriatic, also providing GIA modeling predictions for the city of Venice. The first one is based on the ice-sheet history “K33_j1b_WS9_6”; it assumes a 65 km thick elastic lithosphere and 1 order of magnitude viscosity increase across the 670 km depth seismic discontinuity. The second one, proposed by and named “ICE-7G_NA(VM7)”, is characterized by a 90 km thick lithosphere and by a comparatively milder viscosity increase (by a factor of ∼3). The two models predict distinctly different histories for the GIA-induced RSL variations during the Holocene: the first shows ∼2.2 m of RSL rise in the last 5000 years, whereas the second indicates essentially unvaried RSL during the same period. Note that in previous work , the GIA predictions for the northern Adriatic had a larger uncertainty, with a range of RSL rise between ∼5 and ∼2 m in the last 5000 years. This shows that GIA models are constantly being updated due to improvements in the constraining RSL datasets and in modeling techniques.
Based on the work quoted above, the rate of long-term RSL change in Venice due to the melting of the Late Pleistocene ice sheets does not appear to be tightly constrained . Further uncertainties arise from the effects of the melting of the Würm Alpine ice sheet, whose chronology remains uncertain regarding several aspects . Nonetheless, the long-term rate of Venetian RSL change due to GIA can be assessed in the range between -0.2 and +0.5 mm/year based on the published works cited above. Estimates by and broadly fall within this range, although these works are more pertinent to the Po Plain scale. Note that since the GIA acts on timescales of millennia, these natural contributions to total RSL change will remain constant over the 21st century.
Anthropogenic subsidence
Anthropogenic land subsidence mainly occurs due to extraction of subsurface fluids causing compaction of unconsolidated sediments. This is a process that is widespread in susceptible areas e.g., . Measurements of piezometric level and of vertical land movements are fundamental to quantitatively constrain these processes. Numerical modeling is often used to link the flow of subsurface fluids to the corresponding geomechanical response of the porous medium, although caution is needed. In fact, the paucity of geological data, the imperfect knowledge of forcing processes, and the geomechanical and hydraulic properties generally require significant modeling assumptions and approximations. These techniques have been used to analyze and control the effects of human activities on subsidence in the Venice area.
Prior to 1930, subsidence rates in the Venice region were similar to Holocene rates, suggesting a limited anthropogenic contribution. This is confirmed by both leveling measurements and differences in RSL trends between Venice and Trieste. The post-1930 period is now considered in more detail.
The 1930–1970 period
In the Po and Veneto plains, anthropogenic activities affecting natural land subsidence mainly began in the 1930s due to the overpumping of groundwater and natural gas to support intense civil and industrial development, as shown by geodetic data and reproduced by numerical models . Between World War II and 1970, anthropogenic subsidence was a problem common to the whole northern Adriatic coastline . However, the nature of the withdrawn fluids varied: artesian water in the Venice area, gas-bearing water in the Po Delta, and both groundwater and gas in the Ravenna region. Anthropogenically driven subsidence rates of 10 to 20 mm/year and even higher occurred in certain locations , dominating the RSL change there.
In the Venice area, large quantities of groundwater were pumped to develop the industrial zone of Marghera after 1930. Groundwater pumping was most intensive after World War II, during the period of greatest industrial growth. The six aquifers found in the upper 350 m were progressively exploited (Fig. a); the most intensively used aquifer was between 200 and 250 m depth due to its productive character .
(a) Comparison of the piezometric level in the Marghera industrial area and in Venice from 1910 to 1980 (redrawn after ); (b) average piezometric level between 1952 and 1975 along a leveling line from the mainland to Venice (redrawn after ). (c) Empirical curve (red line) accounting for subsidence in Venice (updated from ; from 2013 onward, the subsidence trend shown in the figure is the one derived from the GPS data of the station Punta della Salute (PSAL); see Table 3). Black dots represent the annual sea-level difference between Genoa and Venice. Green dots represent the height of various benchmarks; (d) same as (b) but for land subsidence. Height differences were estimated by means of three leveling campaigns performed in 1952, 1969 and 1975, respectively. Note the rebound observed in the early 1970s, in correspondence with the partial recovery of the aquifer.
Recent evolution of land subsidence in the historical city center of Venice as measured using geodetic techniques. point out that the uncertainties associated with their SAR estimates represent the ground motion variability at the city scale and are not related to the measurement accuracy.
PeriodSubsidence rate [mm/year]Data sourceReference(s)At the turn of the 19th and 20th century0.9Long-term interpolation of height benchmarks1931–19702.3Difference in tide-gauge records (with reference to Trieste)2.3Long-term interpolation of height benchmarks1953–19735Leveling-0.02 (uplift)Leveling1992–20020.8 ± 0.7SAR2003–20101.0 ± 0.7SAR2008–20201.7 ± 0.5GPS station VEN1 (Riva dei Sette Martiri)Daily solutions provided by the Nevada Geodetic Laboratory (NGL) ; velocity estimated assuming a white + power-law noise model of a priori unknown spectral index (CATS software; ) (consistent with )2014–20200.9 ± 0.6GPS station PSAL (Punta della Salute)Daily solutions provided by NGL ; velocity estimated assuming a white + power-law noise model of a priori unknown spectral index (CATS software; )
Between 1950 and 1970, human-induced subsidence reached 14 cm at Marghera and averaged 10 cm at Venice (Fig. d). The dramatic effects of this loss of elevation became apparent in the exceptional flood (“acqua alta”) of November 1966. The ground beneath Venice is more sensitive to changes in the hydraulic head because of the occurrence of a larger amount of clay in the subsurface : the ratio between subsidence and piezometric decline was 1/100 at Marghera and 1/50 at Venice .
The post-1970 period
After the 1966 flood, the problem of subsidence in Venice received more attention, and drastic measures were taken after 1970 to reduce both industrial and other groundwater extraction. Groundwater consumption in the Marghera area decreased from 500 L/s in 1969 to 170 L/s in 1975 . A corresponding rapid piezometric rise occurred (Fig. a): in 1978 the hydraulic head rose to ground level, re-establishing the levels existing in 1950 . Land subsidence slowed concurrently and stopped; by 1975, a surface rebound of about 2 cm was recorded (Fig. d), equal to 15 % of the total anthropogenic subsidence experienced. This result is consistent with mathematical model results and was interpreted as the elastic response of cohesive soils after recovery.
In recent years, Global Navigation Satellite System (GNSS) and synthetic aperture radar (SAR) measurements have confirmed that the city of Venice is no longer sinking due to groundwater pumping . However, at the local scale, ground movements are still impacted by anthropogenic activities such as new construction and conservation works dedicated to preserving the Venetian architectural heritage. estimated that about 25 % of the city experienced movements attributable to anthropogenic causes. In most cases (15 %) these displacements induced an increase in local subsidence, but in some areas (10 %) the short-term sinking rate was found to be smaller than the natural one. The measured displacement rates range between -10 and 2 mm/year.
Monitoring land subsidence
Over the 20th century, high-accuracy geodetic techniques became available for monitoring land subsidence with unprecedented temporal and spatial resolution. These new data were key to revealing the increasing impact of human activities on the subsidence rate.
The first direct measurements of changes in land elevation were obtained through leveling campaigns based on both local and national networks . Additional information was derived by comparing the tide-gauge records acquired in Venice with those available in neighboring areas subjected to the same absolute sea-level changes . These techniques allowed an unambiguous detection of the impact of anthropic activities on land lowering. It was possible to identify the increase in subsidence between the 1950s and the 1970s caused by severe groundwater extraction . The maximum rate was observed in Mestre in 1968–1969, when the local subsidence reached 17 mm/year . Leveling lines also provided information on the spatial variability in subsidence over a few tens of kilometers. The cone of land depression was found to spread from Marghera, where most of the pumping occurred, towards the Venice area (Fig. ).
During the following decades, leveling measurements performed in 1973 and 1993 recorded the slowdown in subsidence rates and even a small uplift, which followed the dismissal of artesian wells and the diversification of water supply .
The monitoring capabilities further improved during the 1990s with the development of space techniques such as GNSS and SAR . The latest measurements provided by the integrated use of these techniques confirmed that, in Venice, the anthropogenic subsidence due to activities characterized by large-scale and long-term effects ended a few decades ago. However, relevant trends are still observed locally. Subsidence up to 70 and 20 mm/year is found around the inlets where the MOSE is being constructed and in artificial salt marshes, respectively . In addition, spatial patterns in subsidence have been identified at different scales. The average ranges of subsidence rates observed over the lagoon are 3–4, 1–2 and 2–3 mm/year for the northern, central and southern parts, respectively . This reflects both the increase in the thickness of Holocene deposits from the Venice mainland to the lagoon extremes and residual groundwater pumping in the northeastern sector . In the historic center of Venice, the ancient areas are more stable than those urbanized over the last centuries. This is consistent with the older settlements being developed on well-consolidated sand layers, while recent land claims occurred over areas where consolidation processes are still ongoing. Thanks to the high spatial resolution of SAR images, were able to detect the impact of restoration work and new construction down to the single-building scale. This variability in displacement correlates with the nature of the shallow subsoil, the different phases of growth of the city, and the load and foundation depth of the buildings.
Table presents the evolution of subsidence rates in the historical center of Venice, as measured by geodetic instruments over the last century. The estimates refer to different periods of time as well as different locations, and they result from different techniques, for each of which associated uncertainty is available in the literature. Precise leveling allows the measurement of height differences with a mean error ranging from 0.3 to 1 mm in a line of 1 km . The average uncertainty for the vertical component of the GNSS trends is on the order of 0.3 mm/year for time series spanning a decade or more . This estimate increases to 0.4 mm/year when reference frame uncertainties are considered in the error budget . Finally, propose to present the results of the SAR technique over a selected area in terms of average rate and standard deviation of the spatial variability. By doing so, the technique provides insights on the representativeness of the estimated trend at the investigated spatial scale. Therefore, the two GPS-based trend estimates reported in Table are both consistent with the average value obtained by the SAR technique over the whole city.
Estimation of sea-level changes
This section reviews the estimates of centennial average rates of sea-level rise and recent multidecadal trends as well as interannual to interdecadal variations identified in historical sea-level records for the Venice Lagoon and its surroundings and puts them in the context of observed sea-level changes in the Mediterranean Sea and the global ocean. Average rates of sea-level rise are often calculated using some linear fit to the available data. However, given the variety and nonlinearity of the processes known to contribute to sea-level trends on the considered timescales, such estimates should not be intended as necessarily representing a linear process (see Sect. ).
Average rates of sea-level rise over centennial periods
One of the first estimates of the long-term RSL trend at Venice Punta della Salute was made by , who obtained 2.3 ± 0.2 mm/year performing a least-squares fit of the annual means from 1872 to 1941. Since then, several authors have proposed updated estimates by progressively considering newly acquired data and different approaches. A summary of the long-term RSL trends from tide-gauge data proposed during the last 15 years is presented in Table ; estimates derived from data collected over the considered altimetric period (i.e., since 1993) are pointed out. Analogous VLM-corrected RSL estimates are presented in Table . , , and used a linear fit to analyze the RSL data from about 1910 to 2000, obtaining trends between 2.4 and 2.5 mm/year. also modeled the time series by means of the superposition of a straight line and three empirical mode decomposition components, suggesting a slightly larger trend of 2.78 mm/year. isolated the effect of subsidence on the Venetian RSL time series by deriving an empirical curve from leveling data of benchmarks close to the tide-gauge, GPS and SAR interferometry (InSAR) heights (Fig. c). After removing the estimated subsidence from the tide-gauge data, the centennial average rate of sea-level rise of the corrected time series was 1.23 ± 0.15 mm/year for the period 1872–2012 (see Table ). It should be stressed that the trend analysis has little meaning without the correction of subsidence effects because the linear model is otherwise inadequate to represent the Venice time series. In , the application of the same procedure to the neighboring tide gauge of Marina di Ravenna, also located in a rapidly subsiding area, provided a consistent estimate of 1.22 ± 0.32 mm/year (period 1896–2012).
Linear trends of Venice RSL from tide-gauge measurements estimated by various authors in the last 15 years. Errors are SD (68 % confidence) except where noted in the “Confidence” column. The linear fit of observed sea level is used except where annotated. Estimates are grouped based on the period of analysis: long-term and satellite altimetry period.
PeriodSourceTrend [mm yr-1]ConfidenceNotesLong-term 1909–20002.5 ± 0.12.45 ± 0.091914–20002.43 ± 0.2390 %2.78 ± 0.04Trend derived from fit using straight line plus empirical mode decomposition components1872–2019This study2.53 ± 0.1495 %Deseasoned dataAltimetry period 1993–20156.29 ± 1.5399 %Punta della Salute tide-gauge data5.29 ± 1.2799 %Punta della Salute tide-gauge data with inverse barometer effect (IBE) correction1993–2019This study5.01 ± 1.7595 %Deseasoned data
Linear trends of Venetian sea level from tide-gauge data (VLM-corrected RSL) and from satellite altimetry (GSL) estimated by various authors in the last 15 years. Errors are SD (68 % confidence) except where noted in the “Confidence” column. The linear fit of observed sea level is used.
PeriodSourceTrend [mm yr-1]ConfidenceNotesLong-term 1890–20071.20 ± 0.01Deseasoned data1872–20121.23 ± 0.1595 %Deseasoned data1934–20121.20 ± 0.3595 %Deseasoned data1905–20051.3NAComparison amongAdriatic tide gauges1872–2019This study: VLM-corrected tide-gauge data; subsidence estimated after Zerbini et al. (2017), updated to 20191.23 ± 0.1395 %Deseasoned dataAltimetry period 1993-20085.9 ± 1.495 %Altimetry data with dynamic atmospheric correction5.6 ± 1.695 %Tide-gauge data with dynamic atmospheric correction1993–20154.25 ± 1.2599 %Altimetry point near Venice, IBE removal1993–2019This study2.76 ± 1.7595 %VLM-corrected deseasoned tide-gauge data
These corrections allow for a proper comparison with other tide-gauge records in the Mediterranean Sea that span the whole 20th century and are not affected by significant vertical land motions, namely Trieste in the Adriatic and Marseille and Genoa in the northwestern Mediterranean (; Fig. ). The RSL time series of Trieste, Marseille and Genoa exhibit centennial trends between 1.2 and 1.3 mm/year . The estimates agree with the 20th-century trend of sea-level rise of 1.2 ± 0.1 mm/year reported for the same stations by . The 1σ errors are around 0.1 mm/year according to and , and between 0.10 and 0.24 mm/year (90 % confidence) in . obtained uncertainty values between 0.13 and 0.22 mm/year at the 95 % confidence level considering a reduced number of degrees of freedom due to time autocorrelation. Therefore, the 20th-century average rate of sea-level rise in Venice is consistent within uncertainty with those of Marseille, Genoa and Trieste. Accordingly, EOF analysis on annual means from 1901 to 2012 of the corrected time series of Venice and Marina di Ravenna and those of Marseille, Genoa and Trieste yields a leading mode explaining 62 % of variance and corresponding to coherent sea-level variability in the long time series from Mediterranean tide gauges . estimated a trend of 1.3 mm/year for the period 1905–2005 for the Adriatic Sea using a combination of Adriatic tide gauges and manually removing the land subsidence in Venice.
Tables and include an update of the RSL trend calculation for the period 1872–2019, including a comparison between estimates with and without the subsidence contribution following . The subsidence curve by was updated to 2019 by applying a 1 mm/year trend since 2013 based on the SAR estimate by and on the trend exhibited by the PSAL GPS from 2014 onward (Table ). Our updated estimates agree with previous results concerning the magnitude of the full-period average rates of sea-level rise in RSL (2.53 ± 0.14 mm/year) and VLM-corrected RSL (1.23 ± 0.13 mm/year). Subsidence therefore contributed to about half of the total RSL rise in the period 1872–2019 and explains discrepancies between published Venetian RSL trends.
As illustrated in Sects. and , Mediterranean sea level, and hence Venetian sea-level variations, are tightly connected to sea-level variations in the midlatitude eastern North Atlantic, whose underlying processes differ from those in other oceanic basins. Therefore, any statistical consistency between historical Venetian RSL- or VLM-corrected RSL and GMSL rise should not give the false impression that both variables are interchangeable and that any similarity in the historical period necessarily holds in the future. Still, it is relevant to compare estimates of Venetian RSL- or VLM-corrected RSL and GMSL rise during the 20th century to put local changes in the context of global-mean changes since some of the available projections of future Venetian sea level are directly based upon estimates of the GMSL rise (see, for instance, , and ). Venetian sea-level trends are smaller than GMSL trends reported in the fifth assessment report of the Intergovernmental Panel on Climate Change (IPCC-AR5), quantified as 1.7 (1.5 to 1.9) mm/year (likelihood >90 %, period from 1901 to 2010; see ). They are, however, consistent with revisited estimates of historical GMSL rise that include significantly slower rates than reported by the IPCC-AR5 for the pre-altimetry period, e.g., 1.2 ± 0.2 mm/year (90 % confidence interval; ), 1.1 ± 0.3 mm/year (99 % confidence interval; ) and 1.56 ± 0.33 mm/year (90 % confidence interval; ). Figure revisits Venetian RSL- or VLM-corrected RSL and GMSL trends on timescales ranging from interannual to centennial. Clearly, the significant difference between centennial trends in Venetian RSL and GMSL is strongly damped when the contribution of subsidence is removed, confirming the critical role of vertical land motions in determining local RSL variations. Nonetheless, the Venetian VLM-corrected RSL appears to rise at a lower rate than the GMSL over the second half of the 20th century (Fig. a). The GMSL–Venetian sea-level discrepancy observed in the first portion of the record is resolved when uncertainty in GMSL estimate is considered (not shown).
Comparison between trend estimates of GMSL and Venetian sea levels (RSL and VLM-corrected RSL). (a) Temporal evolution of annual-mean sea levels, including GMSL (from ; dashed lines are upper and lower estimates) and Venetian RSL and VLM-corrected RSL. All anomalies with respect to the 1900–1910 mean. (b–d) Maps illustrating linear sea-level trends, estimated via linear regression, for running windows of variable width along the observation period. In panels (b) and (c), black contours illustrate where the GMSL and Venetian sea-level trend estimates do not overlap within 95 % confidence intervals obtained from the linear regression, while green contours indicate where they differ in the sign.
Rates of sea-level rise during the satellite altimetry era
Sea-level measurements acquired with satellite radar altimetry are available from 1993 on, allowing, together with tide gauges, the estimation of recent multidecadal trends from two independent sources. An overall global-mean GSL trend of about 3 mm/year during the satellite altimetry period is consistently reported by several studies . Regional trends can deviate considerably from the global mean eg.. This is also the case of the Mediterranean Sea, which is subject to pronounced spatial and temporal variability (Fig. a), with the entire area of the Adriatic Sea exhibiting positive GSL trends that peak in the northern part of the basin.
Observed and projected trends of Mediterranean sea-level variations. (a) Sea-level trends in the Mediterranean Sea obtained using the ESA SLCCI V2 product over the period 1993–2015. Dots indicate grid points where the trend is not different from zero within the associated error range, estimated as the 90 % confidence interval from linear least-squares regression analysis. (b) Projected thermosteric basin-average sea-level anomalies for the Mediterranean Sea and associated uncertainties related to the Atlantic hydrographic boundary conditions (blue) and to the socio-economic scenarios based on the Special Report on Emissions Scenarios (red) with the regional ocean model NEMOMED8 for the 2000–2100 period (vs. 1961–1990). Uncertainties are the spread among the simulations with only differing Atlantic boundary conditions (blue) and the spread among the simulations with differing socio-economic scenarios (red). Adapted from . (c) Projected northern Adriatic RSL anomalies and associated uncertainties related to socio-economic scenarios (shading: 5th–95th percentile range; line: median). See main text for details on the calculation of uncertainties. Adapted from .
Several altimetric datasets have been used to estimate sea-level trends in the Venice area. estimated a trend of 5.9 ± 1.4 mm/year over the period 1993–2008 for an along-track point about 80 km away from Venice (see their Table 2). obtained trends of 4.18 ± 0.92 mm/year (period 1993–2014) and 3.40 ± 0.99 mm/year (period 1993–2013) for the closest grid point to the Venice tide gauge in the AVISO and SLCCI V1 products, respectively, with both estimates consistent with each other within errors. A reprocessing of the SLCCI V2 dataset over the period 1993–2015 yielded a trend of 4.25 ± 1.25 mm/year , further reduced to 4.03 ± 1.27 mm/year after removing the seasonal signal . Explanations for the differences between trend estimates in these studies include the different time spans, especially for ; different methodological aspects in the spatial characterization of the study area (e.g., closest point vs. area with a certain radius); and the recurrent reprocessing and continuous improvement of the satellite radar altimetry products.
The altimetric trends derived for Venice are typically consistent with those estimated around Trieste over corresponding time spans . This evidence is supportive of a rather uniform sea-level trend in the northern Adriatic (Fig. a; see also ).
A thorough comparison between tide-gauge and altimetric data in Venice is made possible by the availability of independent information on the evolution of the vertical land motion . For consistency with post-processed altimetry data , the tide-gauge time series need to be preprocessed for removal of the local effect of atmospheric forcing (see Sect. ). The most recent trend estimates by provide values of 6.17 ± 1.50 mm/year from in situ data at the Acqua Alta platform (AAPTF), 14 km offshore the Venice coast, and +5.81 ± 1.47 mm/year at Punta della Salute (inside the city center) during the overlapping altimetry period. After subtracting altimetry and AAPTF tide-gauge time series, the residual time series shows a trend of -2.14 ± 0.65 mm/year. This estimate agrees with the trend of 2.17 mm/year extrapolated from Fig. in that represents a best fit of the benchmarks, GPS and PS-InSAR normalized heights.
Table and Fig. also provide updated RSL trend estimates for the period 1993–2019 based on the Punta della Salute tide-gauge data. Our estimates confirm that for the satellite altimetry period the total RSL trend from the tide gauge (5.01 ± 1.75 mm/year, including subsidence) is consistent with uncertainties with some satellite estimates and the tide-gauge estimate by over similar periods. Our estimate for the VLM-corrected RSL trend for the same period is 2.76 ± 1.75 mm/year, again confirming the results by .
Performance of the linear and quadratic regression models applied on the raw annual time series of Venetian RSL and the corresponding climatic component alone (i.e., VLM-corrected RSL) for the periods 1872–2019 and 1993–2019. R2 is the coefficient of determination, AIC is the Akaike information criterion, and FPE is the final prediction error.
RSL VLM-corrected RSL LinearQuadraticLinearQuadratic1872–2019 R20.890.890.670.69AIC400398570570FPE151548471993–2019 R20.580.590.300.30AIC70697676FPE14141818Interannual to interdecadal variability
In addition to long-term changes, the tide-gauge time series of Venetian RSL is characterized by variability in several interannual to interdecadal periods. Significance of observed fluctuations was tested against a red noise hypothesis, namely a lag-1 autoregressive model characterized by high power spectral density at lower frequencies. It is the most widely used model for geophysical purposes since a large class of geophysical processes produces output statistically compatible with the red noise hypothesis (, and references therein). According to , significance is computed through the distribution of an ensemble of surrogate series describing a red noise process with the lag-1 parameter and variance estimated from the analyzed time series. Hereafter, we indicate detection of a statistically significant spectral component around a period of n years with Onyears, where O means order of magnitude. Based on detrended seasonal Venetian RSL for the period 1872–2003, report spectral peaks in the autumn (OND) time series at O22years and at larger multidecadal timescales, at around O2.4years and at around O3.4years, with secondary peaks at around O5years and O8years. In the winter (JFM) time series, they report significant multidecadal variability at O50years and larger; at O3.4years; O8years; and, less apparent, O5years. also identify an O8years component in Venetian RSL variability. An updated spectral analysis based on the Fourier transform applied on autumn and winter raw detrended (second-order polynomial fit) Venetian RSL indicates that the dominant variability modes contained in the Venetian RSL time series over the time interval 1872–2019 are the interannual components at O2.4years and O5years for both autumn and winter series (95 % confidence level), and they account for about 20 % of the total variance of the records. Moreover, the winter time series features the O8years (∼6 % of the total variance) and O50years (∼9 % of the total variance) spectral peaks as highly significant. A secondary peak at O16years is detected at the 90 % confidence level in the autumn series (∼7 % of the total variance). Removal of subsidence does not change the spectral features of the series, except for the O50years component in the winter series, whose significance then only reaches the 90 % confidence level. , , and consistently identify significant decadal variability in the time series of autumn Venetian surge events for the period 1948–2008 (for an updated analysis see , in this special issue).
Continuous wavelet transform analysis on updated and detrended seasonal time series of the Punta della Salute tide-gauge record confirms the presence of statistically significant interdecadal fluctuations in autumn (O20years, period 1960–2000) and interannual (O5years and O8years, periods 1930–1950 and 1970–1990, respectively) and multidecadal (O50years, since 1950) fluctuations in winter (Fig. ). Such fluctuations, however, appear only over limited periods, typically for a few decades or even less. This intermittent recurrent interdecadal variability can significantly impact sub-centennial trend estimates and contribute to the explanation associated spatial features. For instance, in the period between the mid-1960s and the early 1990s, the RSL time series of Venice and Trieste appear almost stationary (Figs. and ). estimated RSL trends to be zero (within the errors) in the 1960–2000 period at the tide-gauge stations of Trieste, Genoa and Marseille. So, stationary sea level characterized the whole Mediterranean Sea during this period but not the Atlantic Ocean, and proposed explanations include an atmospheric contribution mainly consisting of persistent high pressure and an oceanic contribution due to steric changes in deep water masses . conclude that the mean sea-level rise of 1.3 mm/year in the northern Adriatic in the period 1940–2005 is primarily due to ice cap melting, with no trends in atmospheric drivers, upper ocean temperature and surface layer salinity. The different estimates and attribution of multidecadal trends by and reflect the strong interdecadal variability in Mediterranean sea level and the variety of its driving mechanisms. The former study computes the expansion of the water column missing the contribution of the redistribution of mass and produces negligible values of sea-level rise in shallow-water areas, while the latter accounts for it and explicitly considers the Adriatic Sea. Figure c confirms that often bidecadal trends, but occasionally also longer ones, in annual-mean Venetian VLM-corrected RSL are negative and can differ in sign from the GMSL trend. Accordingly, the linear trend can yield a local sea-level anomaly in Venice from the GMSL of about 10 cm (but up to about 20 cm occasionally) over bidecadal and shorter periods and a rather small anomaly (generally <5 cm) over interdecadal and longer periods (not shown).
Composite analyses of autumn (OND, a, c) and winter (JFM, b, d) reanalyzed mean sea-level pressure (MSLP; shading) and 10 m wind (arrows) around Venetian sea-level (a, b) and NAO anomalies (c, d). Composite years are determined based on detrended seasonal values of Venetian VLM-corrected RSL (second-order polynomial trend) and of the Jones’ NAO index (linear trend) below the 10th percentile and above the 90th percentile. Shown differences are for high minus low sea level (a, b) and for negative minus positive NAO (c, d). Analysis is for the period 1872–2015 and for linearly detrended sea-level pressure and wind data. Only grid points where differences are significant at 90 % confidence according to the Wilcoxon rank test are shown.
Interannual to interdecadal autumn (a, c, e) and winter (b, d, f) Venetian sea-level variability since 1872 and its link with the NAO. (a, b) Detrended (second-order polynomial fit) autumn and winter time series of the VLM-corrected RSL from the Punta della Salute gauge record (VLM estimated from ), (c, d) continuous wavelet spectra. Shading (thick black contour) is the portion of the wavelet spectrum exceeding 90 % (95 %) confidence against red noise (lag-1 autoregressive model) hypothesis (see , for details); black lines: cone of influence where edge effects occur; (e, f) wavelet coherence spectra between Punta della Salute data and the Jones NAO index (arrows indicate the phase, with co-phase pointing to the right; thick black contour: 95 % confidence, in red for Punta della Salute RSL data without removal of VLM and detrended as in the main analysis; black lines: cone of influence).
Climatic drivers of Venetian sea-level variations
Variations in Venetian RSL closely depend on sea-level variations in the Adriatic Sea, which in turn are linked to sea-level variations in the Mediterranean Sea. It is well established that basin-wide sea-level fluctuations in the Mediterranean Sea are associated with sea-level variations in the eastern boundary of the North Atlantic Ocean on a broad range of timescales e.g., . Such linkage is largely determined by the water mass exchange between the Atlantic Ocean and Mediterranean Sea through the Strait of Gibraltar e.g., . Section provides further details on the lateral boundary forcing at the Strait of Gibraltar and its drivers. The water mass signal at Gibraltar then propagates within the Mediterranean Sea almost uniformly, yielding basin-wide sea-level variations (e.g., , and ), and overwhelms steric changes as far as basin-average sea-level trends are concerned .
Spatial heterogeneity in Mediterranean sea-level variations is determined by a combination of local steric and ocean circulation changes as well as local atmospheric mechanical forcing, which are described in Sect. . Another potential contributor to local sea-level changes is freshwater horizontal advection , primarily linked to river runoff. At the basin scale, riverine input to the Mediterranean Sea and other freshwater fluxes resulting in net loss of freshwater are regarded as negligible for sea-level variability as they are quickly compensated by changes in the mass transport through the Strait of Gibraltar . However, the conspicuous and noticeably time-varying inflow of freshwaters from the Italian Po River, whose delta is located only 90 km south of the Venice Lagoon (Fig. ), contributes to the shaping of haline properties in the Adriatic Sea with a detectable imprint on Venetian sea levels .
Venetian RSL variability has been investigated in relation to modes of large-scale climate variability, especially concerning teleconnections during wintertime and with the dominant mode of atmospheric variability over the North Atlantic on interannual timescales known as North Atlantic Oscillation, or NAO. The geostrophic winds caused by the large-scale atmospheric pressure anomalies linked to the NAO drive a barotropic sea-level response in Venice, especially in winter. Given its relevance in the scientific literature, the connection between NAO and Venetian as well as Mediterranean sea level is reviewed in detail in Sect. 5.3. In addition to the NAO, statistical connections identified in the literature between Venetian RSL and climatic modes include the atmospheric patterns known as Scandinavian and eastern Atlantic–western Russia , showing prominent variability at interannual to decadal timescales in the autumn and winter seasons, and the Atlantic Multidecadal Oscillation (AMO) , describing multidecadal fluctuations in North Atlantic sea-surface temperature and influencing atmospheric variability over the Euro-Mediterranean region e.g., .
Lateral boundary forcing at the Strait of Gibraltar
Water mass exchange through the Strait of Gibraltar constitutes the critical lateral boundary forcing to determine how sea-level signal propagates from the midlatitude eastern North Atlantic into the Mediterranean basin. The exchange at Gibraltar consists of a strong surface current of relatively fresh and warm water from the ocean and a deep-water current of salt and cold Mediterranean water, outflowing into the ocean and sinking in the North Atlantic in the form of gravity current. Water mass exchanges across the strait critically depend on the number and location of its hydraulic controls, being sub-maximal if subject to only one control in the western part or maximal if the flow is also controlled in the eastern part, with different implications for the characteristics of the circulation. Local dynamics are strongly influenced by tides, which are responsible for the modulation of the water transport and the hydraulic control as well as for the substantial vertical mixing that has been observed .
Nonseasonal mass exchanges through the strait up to the decadal timescale are modulated by fluctuations in the winds over and near the strait . On longer timescales, oceanic contributions mainly including steric changes and mass addition or removal as well as geological processes such as GIA become predominant. These contributions are specific for each basin of the world ocean and especially so for the subpolar and eastern North Atlantic due to the presence of an active thermohaline circulation and to its proximity to the Greenland ice sheet. There, GIA was predominant over the ocean mass contribution to the determination of the upward sea-level trend over the 20th century; since the late 1950s unbiased estimates of steric changes are also available, indicating a contribution to the sea-level rise comparable to GIA . Steric changes particularly contribute to multidecadal sea-level variability in the different ocean basins , with the Atlantic Ocean accumulating heat at a higher rate than the Pacific and Indian oceans e.g.,. So, even though the overall sea-level trend for the subpolar and eastern North Atlantic over the 20th century is comparable to GMSL trend, primarily driven by changes in the ocean mass, it is caused by a different combination of processes .
Recent progress on the simulation of water mass exchange through the Strait of Gibraltar led to a major improvement in the simulation of Mediterranean circulation. demonstrated that the inclusion of explicit tidal forcing in an eddy-resolving Mediterranean model has important effects on the simulated circulation in addition to the expected intensification of local mixing processes. decompose the Mediterranean sea-level signal into two components: first, variations in the eastern North Atlantic sea level, estimated through global coupled climate models, and second, relative variations in the Mediterranean sea level with respect to the eastern North Atlantic, estimated through regional climate models. More recently, have provided an overview of current methods to implement Atlantic sea-level forcing at the lateral boundary of state-of-the-art regional ocean models for the Mediterranean Sea, concluding that the quality of such forcing is essential for appropriate modeling of Mediterranean sea level.
Air–sea interaction within the Mediterranean basin
Local atmospheric mechanical forcing is primarily exerted through local pressure anomalies, associated with the so-called inverse barometer effect (IBE), and wind anomalies, the latter exerting a dominant effect on Mediterranean coastal sea levels, especially by inducing a barotropic oceanic response e.g., . The IBE is quantified by the hydrostatic equation in about 1 cm of sea-level rise per 1 hPa of sea-level pressure drop. quantify the IBE contribution to decadal winter sea-level variability in Trieste at 25 % for the period 1950–2009. The highest IBE contributions to seasonal Venetian sea-level variability over the period 1872–2003 in autumn (about 32 %) and winter (41.5 %) estimated by can be accounted for by the regression model between local sea-level and local sea-level pressure, which could also embed other contributions than IBE alone e.g., .
As far as wind forcing is concerned, the morphology of the Adriatic basin displaying a NW–SE elongation and a shallow northern portion (Fig. ) enhances the effect of the remarkably strong and frequent meridional sirocco wind: prevailing sirocco-favorable conditions are associated with a strengthened northeastward flow over the central Mediterranean, which then favors the piling of Ionian surface waters toward the northern Adriatic, resulting in an increase in Venetian RSL. This Scirocco situation is illustrated for autumn and summer by maps of seasonal anomalies of near-surface wind and sea-level pressure over the Euro-Mediterranean region obtained by compositing the detrended Venetian VLM-corrected RSL series around years characterized by high and low values (Fig. a, b). The analysis is performed on data from version 3 of the NOAA-CIRES-DO Twentieth Century Reanalysis, which contains objectively analyzed four-dimensional weather maps and their uncertainty from the early 19th century to the 21st century . For each considered season, the years were identified using the 10th and 90th percentiles of the seasonal sea-level time series as thresholds. Results indicate that in both autumn and winter relatively high sea level in Venice is associated with anomalously low sea-level pressure over the Mediterranean region, which intensifies over the northern portion of the basin, and with extensively significant northeastward wind anomalies over the Ionian Sea – corresponding to blowing meridional sirocco wind – as well as with an eastward wind anomaly over the Strait of Gibraltar. The morphology of the Adriatic basin displaying a NW–SE elongation and a shallow northern portion (Fig. ) enhances the effect of the sirocco as prevailing sirocco-favorable conditions favors the piling of Ionian surface waters toward the northern Adriatic, resulting in an increase in Venetian RSL. The anomalous wind patterns agree with those identified by obtained with a different set of data.
Fluctuations in surface heat and freshwater fluxes are important drivers of steric changes, through variations in thermal and haline oceanic properties, as well as of ocean circulation changes by affecting processes such as the intermediate- or deep-water formation and transformation e.g., . An example of ocean circulation effects on local sea-level variability are the negative trends observed in the Ionian Sea and southeast of Crete shown in Fig. a .
Linkage with the NAO and other teleconnection patterns
The NAO is the dominant mode of large-scale interannual to decadal atmospheric variability over the North Atlantic. It is commonly characterized by a teleconnection between tropospheric pressure anomalies over the subpolar Arctic and over the subtropical North Atlantic e.g., . Even though the NAO can be identified throughout the year, its anomalies are more stable in winter. The pressure dipole associated with the NAO is tied to changes in the Euro-Mediterranean weather since it produces a meridional displacement of the primary trajectories of the perturbations originating over the Atlantic and affects the frequency and intensity of blocked regimes in the Euro-Atlantic region. Statistical analysis of atmospheric pressure demonstrates coherent large-scale patterns covering the North Atlantic, Europe and the Mediterranean Sea, which explain significant parts of the atmospheric signal variability at interannual and interdecadal scales, particularly during winter. The large-scale coherency of the atmospheric pressure fields means that several of the local atmospheric parameters as well as the oceanic circulation driven by this forcing become correlated. Such linkage drives coherent sea-level changes within the whole Mediterranean basin and, ultimately, in the Venice Lagoon. Despite being studied so far mainly within the framework of interannual to multidecadal climate variability e.g., , the same connections can be relevant for longer-term trends as well, and we therefore include them in this review.
Numerous studies attribute a significant fraction of interannual to decadal winter sea-level variability in the Mediterranean Sea to the NAO e.g., . The NAO exerts its influence on Mediterranean sea level in multiple ways. First, the larger-scale circulation changes associated with the NAO contribute to setting the wind forcing over and near the Strait of Gibraltar, thereby influencing the Atlantic–Mediterranean water mass exchange on decadal timescales . Then, the NAO modulates sea-level pressure, freshwater and buoyancy fluxes and riverine inputs within the Mediterranean Sea and, in particular, the Adriatic Sea e.g.,.
More specifically for Venetian sea level, estimate that about half of the variability in detrended winter Venetian RSL can be explained linearly by the NAO. An updated analysis between seasonal time series of NAO and detrended (second-order polynomial fit) raw Venetian RSL for the period 1872–2019 confirms large values of the correlation statistics (rJFM=-0.68, p∼0 accounting for autocorrelation in the series, rOND=-0.50, p<0.0001). Results of VLM-corrected RSL confirm the strong connection between NAO and detrended Venetian RSL during the cold semester, particularly in winter.
We update previous results about the statistical connection between NAO and Venetian RSL by also performing a wavelet analysis on the 1872–2019 autumn and winter time series of the NAO index and Venetian sea level. The NAO signal is most stable and yields the clearest teleconnections in winter; the autumn season is considered here due to its importance for the occurrence of surges in Venice . The wavelet coherence spectra in Fig. illustrates the link between NAO and Venetian VLM-corrected RSL fluctuations over a broad range of timescales in the decadal to multidecadal range, when both variables are robustly in rough antiphase in both seasons. In autumn inclusion of subsidence in the RSL series disturbs the multidecadal coherent signal (see the red contours in the bottom panels of Fig. ), confirming the importance of accounting for it in studies on the climatic component of local sea-level variability based on tide-gauge data. Temporal variations in amplitude and phase of wavelet coherence and cross-wavelet spectra have been proposed as evidence of lacking a robust connection between geophysical series e.g., . For instance, both characteristics are evident in the wavelet coherence spectrum of the autumn series at the interannual timescales, posing caveats for the truthfulness of the connection. Similarly, in winter coherent interannual variations in antiphase are only observed intermittently.
We therefore complement the statistical connection between NAO and Venetian sea-level time series with maps of seasonal anomalies of near-surface wind and sea-level pressure over the Euro-Mediterranean region obtained by compositing the data around years characterized by positive and negative values of the NAO index (Fig. c, d). As for compositing around Venetian sea level, for each season, the years were identified using the 10th and 90th percentiles of the seasonal time series of the NAO as thresholds. To ease the comparison with the anomalous patterns around sea level, the maps show differences between values under negative NAO minus values under positive NAO. In both autumn and winter, the patterns of sea-level pressure and wind anomalies under different NAO states superpose well on those that correspond to variations in Venetian sea level, including large-scale low pressures enhanced over the northern Mediterranean Sea and sirocco-like northeastward wind anomalies over the Ionian Sea. However, amplitude and spatial extent of the significant anomalies in the atmospheric forcing fields are larger in winter compared to autumn, confirming the weak imprint of the NAO on Venetian sea level in the latter season.
Overall, the seasonal difference in the strength of the NAO imprint on Venetian sea level confirms previous results and supports the investigation of other large-scale precursors of atmospheric forcing of Venetian sea level in autumn. For instance, identify the SCA and EAWR as being related to the interannual to decadal variability of meridional atmospheric flow over the Adriatic Sea in autumn .
Projections
Estimates of the future RSL change in Venice require that all the different components described in the previous sections are considered over the coming decades and combined . It is useful to distinguish GMSL changes which can be derived from the SROCC (Special Report on Ocean and Cryosphere in a Changing Climate; ), regional deviations from GMSL in the Mediterranean and northern Adriatic, and local vertical-land-movement contributions.
Vertical land movements
Projections of future contribution of vertical land motion are available only for the GIA (see Sect. for a quantification). For the other components of vertical land movements, it is only possible to consider their historical variations to contemplate their potential to affect future RSL changes. Estimates of subsidence at the sub-regional scale can be constrained based on observations of past evolution. The sum of sediment compaction, tectonics and GIA is estimated to be about 1.0 mm/year , with a constant rate on a centennial timescale. At the local scale, subsidence trends of a few centimeters per year are still observed in restricted areas such as at the lagoon inlets interested in the construction of mobile barriers against high tides . Controls on groundwater extraction should prevent a return to the large subsidence rates that occurred between 1930 and 1970, but shallow natural processes, notably consolidation, and anthropogenic activities are observed to continue contributing at rates spatially varying from 10 to -2 mm/year . Hence, human activities can contribute significantly to RSL rise, but detailed scenarios are currently not available in the literature, and this is a generic problem . All these natural and anthropogenic contributions have the potential to increase RSL in Venice, exacerbating the negative impact of local sea-level rise, with a contribution on the order of 10 cm per decade at the sub-regional scale and potentially much larger locally.
Sea-level projections for the northern Adriatic and Venice
Sea-level anomalies linked to GMSL rise propagate from the midlatitudes of the eastern North Atlantic to the Venice Lagoon through the Gibraltar and Otranto straits. Sea-level variations along this path are further determined by the mechanical action of the atmosphere and by steric effects associated with changes in temperature and salinity of the water masses. Therefore, beyond local effects on RSL, the projection of sea-level evolution in the Venice Lagoon depends crucially on processes acting on GMSL, sea-level variations across the world ocean and regional patterns of sea-level change inside the Mediterranean basin.
The recent SROCC report summarizes projected GMSL rise estimates, suggesting a likely range from +29 to +110 cm for the year 2100 with respect the 1986–2005 average depending on the future emission scenario, with the low RCP2.6 and the high RCP8.5 providing the lower and the upper limit, respectively. The contribution of ice sheets, especially Antarctica, and the underlying dynamical processes that have been intensively debated recently e.g., provide the main source of uncertainty in future sea-level change projections . The high uncertainty and strong scientific debate on the contribution of Antarctic ice-sheet melting to the future rate of sea-level rise generates the so-called “deep uncertainty” , i.e., a condition where experts lack sufficient knowledge, or parties to a decision cannot agree upon the system processes and futures see also. While high-end GMSL rises are unlikely, they cannot be excluded e.g., and are particularly important for decision-making applications and adaptation planning, especially in decision contexts with low tolerance to uncertainty . According to , the sea-level rise at the subtropical and midlatitudes of the eastern North Atlantic (and the Mediterranean Sea itself) will only have small deviations (less than 10 %) from the global-mean change.
Analysis of the possible deviations of future mean sea level in the Mediterranean basin from the GMSL has been attempted using dynamical and statistical models. Dynamical models of the circulation inside the Mediterranean Sea allow direct estimates of the mechanical effect of the atmosphere on the circulation (wind stress and IBE). There is consensus that this is a small contribution with changes generally less than 10 cm . The steric effects are computed from temperature and salinity changes using a diagnostic offline computation. This computation obviously depends on the water depth and tends to zero at the coastline. Therefore, for producing a realistic estimate of the actual sea-level change, they should be integrated with the computation of the associated water mass redistribution. In fact, the spatial variations in the resulting sea-level change follow the bathymetry of the Mediterranean Sea and are rather small (<5 cm) over shallow-water areas such as the northern Adriatic Sea e.g.,. The overall steric sea-level change is the consequence of two contrasting effects: thermosteric expansion (associated with warming of water masses) and halosteric contraction (associated with increasing salinity of water masses). Most studies agree that the former is larger, and there is also a freshening effect of the Atlantic inflow across the Strait of Gibraltar whose magnitude is poorly constrained. Several studies contributed to an assessment of the overall steric effects at the Mediterranean scale, with differences depending on periods, models, scenarios, and the representation of water exchange between Mediterranean Sea and Atlantic Ocean . There are various basin-wide Mediterranean steric sea-level projections. For instance, a range of 2 to 7 cm for the mid-century sea-level anomaly (2050 with respect to 2001) under the SRES-A1B scenario was reported . By the end of the 21st century (i.e., 2070–2099), for the steric component alone, found a 13 cm sea-level anomaly (with respect to 1960–1990) under the A2 scenario, while reported sea-level anomalies over the same periods that range between 34 and 49 cm considering a six-member ensemble and the A1, A1B and B1 scenarios (see Fig. b). The magnitude of the steric sea-level rise for the Mediterranean Sea is similar to estimates of the steric GMSL rise, e.g., 29–45 cm (67 % confidence interval) for the RCP8.5 scenario between 2100 and 2000 by .
Estimates of future sea-level rise in the northern Adriatic Sea from a statistical model are provided by . The computation is based on a linear regression linking the deviation of the sub-regional sea level to changes in mean sea-level pressure, water temperature and salinity, meant to represent the mechanical atmospheric forcing, steric effects and the redistribution of mass inside the basin . concluded that regional effects, at the end of the 21st century for the A1B scenario, result in a deviation in the range from 2 to 14 cm from the sea level of the eastern Atlantic outside the Strait of Gibraltar; they also concluded that the main contribution to local sea-level rise is caused by remote effects, such as mass inflow across the Strait of Gibraltar, which was conservatively estimated without considering a likely future acceleration and information from climate projections.
Adopting the “mid-range” and “high-end” climate change scenarios defined by , have evaluated the effect of contemporary and future terrestrial ice melting on future regional RSL evolution across the Mediterranean Sea by solving the sea-level equation including the contribution of glaciers, ice caps, and the Greenland and Antarctic ice sheets. They found that terrestrial ice melting will be responsible for a sea-level rise in the northern Adriatic of ∼8 and ∼17 cm between 2040–2050 relative to 1990–2000 in the two considered scenarios, respectively. Since the sources of terrestrial ice melting are mostly located in the far field of the Mediterranean Sea, these contributions shall be largely uniform across the whole Adriatic Sea.
Northern Adriatic RSL projections informed by climate projections can be obtained by summing up the future regional contributions of sterodynamic effects, which correspond to changes in ocean density and circulation corrected from the IBE; melting of mountain glaciers and ice sheets; land water; and vertical land motions e.g., . provide probabilistic sea-level projections for Venice as part of a global set of local sea‐level projections for three different representative concentration pathways. Their projections build on a combination of expert community assessment (the IPCC-AR5), expert elicitation e.g., and process modeling (e.g., the 5th phase of the Coupled Model Intercomparison Project or CMIP5) for most sea-level contributors. The “background non-climatic local sea-level change” corresponding to GIA, tectonics and other non‐climatic local effects was derived by applying a Gaussian process model to tide-gauge records. This background linear estimate and its uncertainty (0.72 ± 0.33 mm/year for Venice; see Supplement and Table 8 in ) are then included in the projections, together with the other components. The 5th–95th percentile range of the resulting sea-level change projections in the year 2100 is 29–79 cm for RCP2.6 and 41–107 cm for RCP8.5.
However, two critical aspects are needed to obtain reliable sea-level projections for Venice. First, a full characterization of the vertical land motions in Venice requires an explicit local focus beyond what is provided by a mere extrapolation on the site from regional patterns of GIA, tectonics and other non‐climatic contribution alone. In fact, the subsidence estimate by lies in the lower range of the available instrumental estimates, especially as far as the most recent period is considered (see Table 3). Then, the sterodynamic component is derived from the outputs of the coupled climate–model simulations performed within CMIP5. The rather coarse resolution of coupled climate models prevents an accurate representation of small-scale processes (e.g., water exchange at Gibraltar; see ), which in turn affects regional sea-level estimates . Another important caveat for multi-model assessments is their reliability in coastal regions where the contributing models may resolve the coastline and bathymetry peculiarities differently, thus yielding local anomalies in the gridded multi-model output that may reflect a bias originated by heterogeneous spatial resolutions across models rather than a physical process e.g.,.
On this premise, we propose probabilistic projections of northern Adriatic RSL for two climate scenarios (RCP2.6 and RCP8.5) and one high-end scenario following and . The method allows the inflation of the uncertainty in projections of the sterodynamic component by accounting for the low confidence in projections of coastal sea-level rise obtained from the limited number of global circulation models participating in CMIP5 and covering the Mediterranean Sea (see Fig. in ). Specifically, the Mediterranean sterodynamic sea-level projections are estimated by relying on those of the Atlantic area near Gibraltar. Other mass contributions to sea level (i.e., glaciers, ice sheets, land water) have a global effect due to the addition of water mass to the ocean (barystatic sea-level rise) and a regional effect through instantaneous changes in the Earth's gravity, Earth's rotation and solid-Earth deformation (GRD; GRD-induced RSL change). Both contributions are combined into a geographical pattern called the barystatic-GRD fingerprint and are proportional to the land water mass change. For each scenario, uncertainties correspond to the combined uncertainty in each sea-level component calculated as the square root of the sum of the squares of each component uncertainty. Note, however, that contributions that correlate with global-mean air temperature, namely the sterodynamic and ice-sheet surface mass balance components, have correlated uncertainties and are therefore added linearly according to the following equation:
σtot2=(σsterodynamic+σsmb-a +σsmb-g)2+σGlac2+σLW2+σdyn-a2+σdyn-g2,
where smb stands for surface mass balance, Glac for glaciers, LW for land water, and dyn-a and dyn-g for dynamic Antarctic and Greenland (see , for more details). Considering an annual subsidence rate of 1 mm/year, northern Adriatic RSL is projected to rise by 47 cm (likely range 32–62 cm) for the RCP2.6 scenario and by 81 cm (likely range 58–110 cm) for the RCP8.5 scenario by the end of the 21st century with respect to the reference period 1986–2005. These projections start to diverge after 2050 and, excluding the contribution of subsidence of 10.5 cm, yield a total range of local sea-level rise (VLM-corrected RSL) by 2100 between 21 and 100 cm, i.e., ∼10 % lower than GMSL rise (Fig. c). The high-end scenario – obtained by selecting, for each sea-level component, the highest physically based estimate found in the literature (see for details) – shows that, by 2100, northern Adriatic sea level could unlikely but possibly rise by more than 1.8 m. Note that for this high-end scenario, the Antarctic component contributes nearly half of the sea-level change in 2100. This estimate agrees with the analysis by . Our projections largely overlap with those provided by for both scenarios, indicating that they are broadly consistent with each other. Differences are likely due to differences in the methods, models and assumptions employed by both studies, which requires a dedicated investigation to be fully understood.
In conclusion, while our understanding of past RSL change in Venice has improved, the large uncertainty in the magnitude of future GMSL rise remains a major scientific challenge . Hence, public policy needs to recognize this uncertainty and monitor sea-level change as part of the management of this uncertainty, including drawing on relevant experience from elsewhere e.g.,. Regional effects can determine differences on the order of 10 cm between the mean Mediterranean sea level and the GMSL and within different parts of the Mediterranean basin itself as seen in the historical record. Changes over interdecadal periods can also distort the detection of forced trends over rather long periods of time e.g.,. However, scientific literature provides no evidence for a future deviation, on a centennial timescale, of the local sea level at the Venetian coastline from the GMSL that is larger than the abovementioned order of magnitude. RSL can differ by more as land movements and regional atmospheric patterns could provide additional and important contributions.
Gaps of knowledge and opportunities for progress
This literature review identifies several issues where progress is needed, including (1) improving satellite observations of sea-level change and integrating them with tide-gauge measurements, (2) improving monitoring and prediction of vertical land motions, (3) improving the simulation of Mediterranean oceanic circulation, (4) determining the shape of the historical RSL and VLM-corrected RSL trends, and (5) reducing uncertainty in estimates of regional effects on future climate-induced sea-level rise.
Concerning observations, altimeter data are fundamental for providing a regional perspective and reaching robust conclusions on observed local sea-level rise as they provide an independent source of information from the local tide-gauge data. The contribution of coastal altimetry is considered essential to within 0–10 km to link the sea-level changes derived from satellites with those measured at tide-gauge locations . The extension of the satellite-based sea-level record toward the coast with measurement quality comparable to the open ocean allows the detection of differences in the rate of coastal sea-level rise from the open ocean e.g., and their implications for coastal hazards . It has been shown that by improving processing, it is possible to make more accurate sea-level measurements in coastal zones . Progress has been made in fitting the radar signal (the so-called re-tracking) in order to extract a robust estimate of the distance between the satellite and sea surface. The ALES retracker, with a proper threshold on error, recovers significantly more data in the 10 km near the coast . There are several satellite radar altimetry operational products (along-track and gridded) dedicated to the monitoring of open-ocean sea level whose quality is constantly improving. Various experimental coastal altimetry products are also now available and validated in some regions, allowing sea-level research in the coastal zone . An updated table is accessible at https://www.coastalt.eu/community#datasets (last access: 5 August 2021). Importantly, merging altimeter data from different missions requires homogenous reprocessing and minimization of drifts and systematic biases between missions. Overall, progress is expected in terms of improving the altimetry estimates closer to the coast, perhaps with different instruments, and understanding the relationship between the coastal and close-to-the coast sea-level changes in a complex environment .
For the northern Adriatic, an opportunity for progress is provided by the ESA SLCCI extension (CCI+), which will process along-track data from additional satellite missions using re-tracked data, dedicated coastal geophysical corrections and improved editing that will then be combined in a global grid with higher resolution near the coast (Anny Cazenave, personal communication, 2020). Also, the novel GNSS-derived path delay plus (GPD+) correction now provides accurate wet tropospheric delays . The wide-swath interferometry will be deployed for the first time on the Surface Water Ocean Topography mission to provide for the first time sea-level imaging that will solve the limitation of existing satellite sea-level measurements being available only along tracks .
Concerning the monitoring of vertical land motions, integrated systems have been shown to offer the best approach to the study of subsidence : GNSS provides point-wise continuous positioning with respect to a global reference frame; SAR offers spatially dense measurements of surface displacements relative to a ground target selected as a reference point. However, while these techniques can support the investigation of present subsidence patterns with unprecedented detail (i.e., at the scale of individual buildings), future scenarios remain difficult to construct, with the anthropogenic component of vertical land movements being the most difficult to assess . Historical observations showing the potential of anthropogenic subsidence in Venice to be on the order of tens of centimeters per decade demonstrate the need of continued careful regulation of land and groundwater use and monitoring of local subsidence. Together with subsidence process understanding and simulation, this might be used to develop high- and low-subsidence scenarios, respectively.
Concerning the simulation of Mediterranean oceanic circulation, considerable efforts have been invested over recent years into developing and applying regional climate and ocean circulation models approaching the issue of dynamical downscaling from different perspectives e.g., . This included coordinated international regional climate modeling activities (e.g., MedCORDEX) and regional reanalysis for the Mediterranean Sea obtained through numerical simulations with data assimilation . However, despite recent progress in the representation of lateral boundary forcing at the Strait of Gibraltar, there are several aspects that remain poorly understood or are worth deeper investigation. For example, small changes in the salinity difference between Mediterranean and Atlantic waters around a threshold of 2 psu can determine shifts in the simulated hydraulic regime within the Strait of Gibraltar, from sub-maximal to maximal e.g., . Accordingly, a scenario involving a positive trend in the salinity difference can result in a partial isolation of the Mediterranean Sea from the rest of the world ocean . How nonlinear interaction between large-scale ocean variations and local strait phenomena may sustain an abrupt change in the saltwater and freshwater transport between the Mediterranean and the Atlantic and a shift in the Mediterranean mean circulation remains to be investigated in a comprehensive modeling framework. It is also possible that the large-scale ocean circulation in the North Atlantic can impact the Mediterranean Sea level through ocean–atmosphere feedbacks e.g., , but this requires further investigation . Other exchanges are also relevant: for instance, water mass exchanges with the Black Sea through the Turkish Straits remain idealized in current simulations, and their effect is underrepresented in future projections. Water mass exchange between shelf and ocean is performed through cascading processes, which are hardly reproduced by both regional and global solutions e.g., . In this context, the unstructured approach adopted by is promising where the system of inter-connected basins formed by the Mediterranean, the Marmara, the Black and the Azov seas was numerically investigated using a unique computational mesh allowing for a seamless transition between different spatial scales, from narrow straits to open sea (see also , in this special issue).
Poorly simulated internal ocean variability also provides potential weakness to projected circulation changes in the Mediterranean Sea, which calls for a stronger focus on the validation of regional ocean models regarding interior and abyssal dynamics linked to fundamental oceanographic processes. The Adriatic Sea is the only Mediterranean sub-basin in which the evaporation–precipitation–runoff budget is negative: the buoyancy flow at the Otranto Strait is either positive or negative depending on the predominance of production of dense water within the Adriatic or of the inflow of the Levantine Intermediate Water, respectively. Numerical simulations indicate that a nonlinear convection–mixing feedback can favor hysteresis in the Adriatic Sea with multiple equilibria encompassing estuarine and anti-estuarine circulation . Such behavior could have important implications for future sea-level variability in the Venice Lagoon. Overall, even under accurate representation of global steric and mass addition from the Atlantic, projections of Mediterranean sea-level change from current regional ocean models would be reliable only in the basin mean tendencies. Further, comparison of the regional simulations with satellite-derived data highlights local biases in the historical sea-surface height patterns and trends as well as large inter-model heterogeneity in projected changes at the local scale driven by differences in simulated circulation changes. Simulation–data comparison using updated hydrographic datasets, such as EN4 or , could be extremely valuable.
The higher rate of RSL rise observed in recent decades compared to the longer-term estimate (Table 4) raises the practical question of how different statistical models to extract the trend for time series analysis affect the results. Here, trends refer to the long-term movement in a time series, which may be regarded, together with the oscillation and random component, as generating the observed values . We have evaluated the performance of a linear and a quadratic regression model in the RSL times series, including the raw series and the climatic component alone, and for the periods 1872–2019 and 1993–2019. According to several skill metrics including R2, AIC and FPE and references therein (Table 6), the quadratic model only slightly outperforms the linear model for both periods and both series. As far as the higher rates of Venetian RSL rise observed in recent decades are concerned, the simple acceleration expressed statistically in terms of quadratic fitting seems to be insufficient, and further methods could be explored (but as a note of caution see, for example, ). The presence of substantial variations in multidecadal trends of sea level in Venice, together with methodological differences, explains the wide range of estimates of average sea-level rise obtained by different authors considering different periods. Still, certain aspects of the forcing of Venetian sea-level variability require improved quantifications and characterization (for instance regarding the atmospheric forcing of local sea-level variability), and this will better constrain the range of projected future changes.
Concerning understanding the future evolution of Venetian RSL, we identify two critical gaps of knowledge. First, the significant interannual to multidecadal variability observed in Venetian sea level requires new studies to address their implications for near- and mid-term sea-level prediction. A possible focus could be, for instance, the impact on Venetian sea level of a multiannual period of a positive phase of the NAO as observed in the early 1990s. Second, in the long term, there is no consideration of the implications of high-end GMSL rise for RSL in Venice. Research on both aspects would contribute to improved understanding of forcing mechanisms of Venetian sea-level changes, risk assessment and management, and more effective scientific communication.
Finally, this literature review builds on recent studies on local and regional sea-level change projections e.g., to develop explicit local RSL change scenarios for Venice by combining the uncertainty ranges of future projections for the individual contributing processes. Understanding the causes of the differences between the projections obtained in this study and those obtained by can help clarify the possible implications of the choice of methods, models and assumptions to build future scenarios of Venetian sea level and by implication for any site. Further approaches to objectively combining such uncertain estimates are hoped to be tested based on qualitative criteria (e.g., considered process, statistical and numerical framework) or quantitative metrics, such as relative or absolute model skills in representing relevant physical features (e.g., boundary forcing at the Strait of Gibraltar). It will also be important to compare our estimates with updated scenarios of Venetian sea-level future change that are expected from the 6th phase of the Coupled Model Intercomparison Project .
Conclusions
The city of Venice and the surrounding lagoonal ecosystem are critically affected by variations in RSL height driven by a host of diverse processes. These encompass oceanic processes driving sea-level variations from diurnal astronomical oscillations to climatic interannual to multi-centennial fluctuations and vertical land movements causing RSL variations on timescales ranging from a few decades – due to, for example, anthropic activities – to multi-millennial trends due to tectonic activity. This review summarizes and reassesses recent progress in the estimation, understanding and prediction of the individual contributions to RSL by exploiting new observational datasets, improved statistical methods and more realistic numerical simulations of ocean and Earth system components achieved in the past decade.
Our estimate of the historical long-term linear trend of Venetian VLM-corrected RSL is 1.23 ± 0.13 mm/year (from 1872 to 2019). Inclusion of subsidence gives a centennial RSL rise of about 2.5 mm/year. Looking to the future, the effects of both subsidence and climate-induced sea-level rise will have profound implications for Venice. By 2100, natural local subsidence is expected to result in an RSL rise of about 10 cm relative to the late 20th century. Projected climatically induced Venetian sea-level rise from estimates for the GMSL corrected for the regional redistribution of mass contribution components such as glaciers, ice sheets and groundwater is in the range from 21 to 52 cm (from 48 to 100 cm) by 2100 for the RCP2.6 (RCP8.5) socio-economic scenario. Ice-sheet melting provides a highly uncertain contribution. An RSL rise in Venice exceeding 180 cm by 2100 is an unlikely but plausible high-end change under strong melting of Greenland and Antarctica. These estimates neglect the effect of atmospheric forcing of local sea level, which potentially contributes additional uncertainty, estimated here in the range of about 10 cm. The uncertainty range of the RSL rise by 2100 in Venice across all emission scenarios obtained here is thus very large: a minimum low-end scenario is about 10 cm, corresponding to unlikely combinations such as either a low-emission scenario with no subsidence and regional atmospheric interaction producing a negative effect or subsidence continuing at a historical rate and all other effects compensating each other; the upper limit of the likely range considering RCP8.5 projections is about 120 cm. Including unlikely high-end scenarios would raise this limit to 180 cm, which could approach 2 m under adverse regional atmospheric forcing. An additional contribution could be produced by anthropogenically driven subsidence. While, in general, the resulting effect of regional climatic processes could either attenuate or increase regional RSL with respect to GMSL, local subsidence will necessarily exacerbate it. Further, because of the differential rates of subsidence observed across the Venice area, the land movement estimates at the tide gauges may underestimate the risk for other parts of the city.
The review highlights several important advances. Centennial RSL variations are now known to be spatially heterogeneous within the lagoon and the city due to differential vertical ground movements. Hence local trend estimates are not expected to be representative of the city or lagoon as a whole. Due to the nonlinear variations caused by ground subsidence, a linear detrending of Venetian RSL time series is unsuitable unless data are preliminarily corrected for the effect of vertical land motion. Remote climate forcing from the Atlantic sector via atmospheric and oceanic processes critically contributes to interannual to multidecadal Venetian RSL variability. Atlantic hydrographic boundary conditions are a major source of uncertainty for future projection of Mediterranean sea levels: uncertainty in water mass flows at the Strait of Gibraltar yields an ensemble spread between simulations comparable to that determined by uncertainty in greenhouse gas emissions. We confirm the existence of a strong link between interannual and interdecadal variability observed in Venetian sea levels and in the large-scale atmospheric circulation over the North Atlantic during the cold semester, particularly with the North Atlantic Oscillation (about 46 % of shared variability for the detrended winter average time series in the period 1872–2019).
The review highlights several major gaps of knowledge as well. Among these, altimetry data are recorded rather far from Venice and may not represent the lagoon RSL variability; uncertainties in geologic trends remain difficult to assess; a reliable framework is lacking to combine uncertain future estimates of RSL change due to individual contributions, which provides for a major opportunity for progress to better constrain ranges of future projections; and historical evidence demonstrates that subsidence can be temporarily dominated by the anthropic component. This last point shows the importance of sustaining the management regime that brought this anthropogenic component under control across the lagoon and possibly of strengthening it to bring the ongoing small-scale anthropic subsidence under control. Finally, whereas several studies explored scenarios of RSL changes in Venice at the end of the 21st century under global climate change, near-term and transient predictions have not been attempted yet, and high-end scenarios have not been a subject of explicit focus.
Code and data availability
Relevant data used in this study can be obtained at 10.5281/zenodo.5139890.
Author contributions
DZ and SB coordinated the paper with help from FR. Specific contributions to the sections are as follows (LA = leading author, CA = contributing author). Section 1 was written by DZ (LA) with help from SB (CA). Section 2 was written by FR and SZ (LA) with help from MT, SV, GW and SZ (CA). Section 3 was written by SB (LA) with help from FA, EC, GS, GW and SZ (CA). Section 4 was written by DZ and FR (LA) with help from PL (CA). Section 5 was written by DZ and PL (LA) with help from MT, GS, FA, AA, VF, AR and CF (CA). Section 6 was written by DZ and PL (LA) with help from RT, FR, SB and RJN (CA). Section 7 was written by DZ (LA) with help from SB, PL, RJN and CF (CA). Section 8 was written by DZ, PL and RJN (LA) with help from all authors. SR created Figs. 1 and 8. DZ created Figs. 2, 4, 6, 7 and 9. SB created Fig. 5. SV created Fig. 3.
Competing interests
The authors declare that they have no conflict of interest.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Special issue statement
This article is part of the special issue “Venice flooding: understanding, prediction capabilities, and future projections”. It is not associated with a conference.
Acknowledgements
The authors want to thank the European Space Agency (ESA), which funded the Climate Change Initiative to produce a climate quality record of sea level from satellite altimetry. This work is partially funded by the ESA under Phase 2 Bridging Phase to CCI (contract no. 4000109872/13/I-NB – Contract Change Notice 6). The Centro Previsioni e Segnalazioni Maree of the Venice Municipality is acknowledged for providing the tide-gauge data. Thanks are due to Francesco De Biasio (CNR-ISP) for the joint work during the CCI project. Wavelet coherence spectra are calculated using the cross-wavelet and wavelet coherence package by Aslak Grinsted, John Moore and Svetlana Jevrejeva. Linear regressions are performed with the “regress” function of MATLAB. The code to compute regional sea-level projections was formerly developed by Gonéri Le Cozannet. Twentieth Century Reanalysis V3 data were provided by the NOAA/OAR/ESRL PSL, Boulder, Colorado, USA, from their website at https://psl.noaa.gov/ (last access: 5 August 2021). The authors thank the two anonymous reviewers for their comments and insight on the manuscript, which greatly helped to improve the paper.
Financial support
Scientific activity was performed as part of the research program Venezia2021, coordinated by CORILA, with the contribution of the Provveditorato for the Public Works of Veneto, Trentino Alto Adige and Friuli Venezia Giulia. Giogio Spada is funded by a FFABR (Finanziamento delle Attività Base di Ricerca) grant of MIUR (Ministero dell’Istruzione, dell’Università e della Ricerca).
Review statement
This paper was edited by Uwe Ulbrich and reviewed by two anonymous referees.
ReferencesAblain, M., Becker, M., Benveniste, J., Cazenave, A., Champollion, N.,
Ciccarelli, S., Jevrejeva, S., Le Cozannet, G., Nicoletta, L., Loisel, H., Long, N., Maisongrande, P., Mallet, C., Marcos, M., Menendez, M., Meyssignac, B., Plater, A., Raucoules, D., Taramelli, A., Vignudelli, S., Valentini, E., Woodworth, P., and Woppelmann, G.: White paper: Monitoring the evolution of coastal zones under various
forcing factors using space-based observing systems, International Space
Science Institute (ISS), Bern, Switzerland,
available at: https://hal-brgm.archives-ouvertes.fr/hal-01413107 (last access: 5 August 2021), 2016.Adloff, F., Somot, S., Sevault, F., Jordà, G., Aznar, R., Déqué,
M., Herrmann, M., Marcos, M., Dubois, C., Padorno, E., Alvarez-Fanjul, E., and Gomis, D.: Mediterranean
Sea response to climate change in an ensemble of twenty first century
scenarios, Clim. Dynam., 45, 2775–2802, 10.1007/s00382-015-2507-3,
2015.Adloff, F., Jordà, G., Somot, S., Sevault, F., Arsouze, T., Meyssignac, B.,
Li, L., and Planton, S.: Improving sea level simulation in Mediterranean
regional climate models, Clim. Dynam., 51, 1167–1178,
10.1007/s00382-017-3842-3, 2018.
Alessio, S. M.: Digital signal processing and spectral analysis for scientists:
concepts and applications, Springer, Cham, Switzerland, 2016.Allen, M. and Smith, L.: Monte Carlo SSA: Detecting irregular oscillations in
the presence of colored noise, J. Climate, 9, 3373–3404,
10.1175/1520-0442(1996)009<3373:MCSDIO>2.0.CO;2, 1996.Amitai, Y., Ashkenazy, Y., and Gildor, H.: Multiple equilibria and overturning
variability of the Aegean-Adriatic Seas, Global Planet. Change, 151,
49–59, 10.1016/j.gloplacha.2016.05.004, 2017.Amorosi, A., Colalongo, M., Fiorini, F., Fusco, F., Pasini, G., Vaiani, S., and
Sarti, G.: Palaeogeographic and palaeoclimatic evolution of the Po Plain from
150-ky core records, Global Planet. Change, 40, 55–78,
10.1016/S0921-8181(03)00098-5, 2004.Androsov, A., Fofonova, V., Kuznetsov, I., Danilov, S., Rakowsky, N., Harig, S., Brix, H., and Wiltshire, K. H.: FESOM-C v.2: coastal dynamics on hybrid unstructured meshes, Geosci. Model Dev., 12, 1009–1028, 10.5194/gmd-12-1009-2019, 2019.Antonioli, F., Ferranti, L., Fontana, A., Amorosi, A., Bondesan, A.,
Braitenberg, C., Dutton, A., Fontolan, G., Furlani, S., Lambeck, K., Mastronuzzi, G., Monaco, C., Spada, G., and Stocchi, P.:
Holocene relative sea-level changes and vertical movements along the Italian
and Istrian coastlines, Quaternary Int., 206, 102–133,
10.1016/j.quaint.2008.11.008, 2009.Antonioli, F., Anzidei, M., Amorosi, A., Presti, V. L., Mastronuzzi, G.,
Deiana, G., De Falco, G., Fontana, A., Fontolan, G., Lisco, S., Marsico, A., Moretti, M., Orrù, P. E., Sannino, G., Serpelloni, E., and Vecchio, A.:
Sea-level rise and potential drowning of the Italian coastal plains: Flooding
risk scenarios for 2100, Quat. Sci. Rev., 158, 29–43,
10.1016/j.quascirev.2016.12.021, 2017.
Arca, S. and Beretta, G.: Prima sintesi geodetico-geologica sui movimenti
verticali del suolo nell'Italia Settentrionale (1897–1957), Bollettino di
Geodesia e scienze affini, 44, 125–156, 1985.
Armi, L. and Farmer, D. M.: The flow of Mediterranean water through the strait
of Gibraltar. The flow of Atlantic water through the Strait of Gibraltar,
Prog. Oceanogr., 21, 1–105, 1988.Artale, V., Calmanti, S., Malanotte-Rizzoli, P., Pisacane, G., Rupolo, V., and
Tsimplis, M.: The Atlantic and Mediterranean Sea as connected systems, Developments in Earth and Environmental Sciences, 4,
283–323, 10.1016/S1571-9197(06)80008-X, 2006.Artale, V., Calmanti, S., Carillo, A., Dell'Aquila, A., Herrmann, M., Pisacane,
G., Ruti, P. M., Sannino, G., Struglia, M. V., Giorgi, F., Bi, X., Pal, J. S., Rauscher, S., and The PROTHEOUS Group: An
atmosphere–ocean regional climate model for the Mediterranean area:
assessment of a present climate simulation, Clim. Dynam., 35, 721–740,
10.1007/s00382-009-0691-8, 2010.Bakker, A. M., Wong, T. E., Ruckert, K. L., and Keller, K.: Sea-level
projections representing the deeply uncertain contribution of the West
Antarctic ice sheet, Sci. Rep., 7, 1–7,
10.1038/s41598-017-04134-5, 2017.Bamber, J. L. and Aspinall, W. P.: An expert judgement assessment of future sea
level rise from the ice sheets, Nat. Clim. Change, 3, 424–427,
10.1038/nclimate1778, 2013.Barriopedro, D., García-Herrera, R., Lionello, P., and Pino, C.: A
discussion of the links between solar variability and high-storm-surge events
in Venice, J. Geophys. Res.-Atmos., 115, D13101, 10.1029/2009JD013114, 2010.Battistin, D. and Canestrelli, P.: 1872–2004: la serie storica delle maree a
Venezia, Proc. 3rd International Symp. on Land Subsidence,
available at: https://www.comune.venezia.it/content/download-pubblicazioni (last access: 5 August 2021),
2006.Benveniste, J., Cazenave, A., Vignudelli, S., Fenoglio-Marc, L., Shah, R.,
Almar, R., Andersen, O. B., Birol, F., Bonnefond, P., Bouffard, J., Calafat, F., Cardellach, E., Cipollini P., Le Cozannet, G., Dufau, C.,
Fernandes, M. J., Frappart, F., Garrison, J., Gommenginger, C., Han, G., Høyer, J. L., Kourafalou, V., Leuliette, E., Li, Z., Loisel, H., Madsen, K. S., Marcos, M.,
Melet, A., Meyssignac, B., Pascual, A., Passaro, M., Ribó, S., Scharroo, R., Song, Y. T., Speich, S., Wilkin, J.,
Woodworth, P., and Wöppelmann, G.:
Requirements for a coastal hazards observing system,
Front. Mar. Sci., 6, 348, 10.3389/fmars.2019.00348, 2019.Birol, F., Fuller, N., Lyard, F., Cancet, M., Nino, F., Delebecque, C., Fleury,
S., Toublanc, F., Melet, A., Saraceno, M., and Léger, F.: Coastal applications from
nadir altimetry: Example of the X-TRACK regional products, Adv. Space Res., 59, 936–953, 10.1016/j.asr.2016.11.005, 2017.Blewitt, G., Hammond, W. C., and Kreemer, C.: Harnessing the GPS data explosion
for interdisciplinary science, Eos, 99, 1–2, 10.1029/2018EO104623,
2018.Bonaduce, A., Pinardi, N., Oddo, P., Spada, G., and Larnicol, G.: Sea-level
variability in the Mediterranean Sea from altimetry and tide gauges, Clim. Dynam., 47, 2851–2866, 10.1007/s00382-016-3001-2, 2016.
Bondesan, M., Cibin, U., Colalongo, M., Pugliese, N., Stefani, M., Tsakiridis,
E., Vaiani, S., and Vincenzi, S.: Benthic communities and sedimentary facies
recording late Quaternary environmental fluctuations in a Po Delta subsurface
succession (Northern Italy), Proceedings of the 2nd and 3rd Italian Meeting on Environmental Micropaleontology, 11, 21–31, edited by: Coccioni, R., Lirer, F. and Marsili, A., The Grzybowski Foundation Special Publication, 2006.
Bortolami, G., Carbognin, L., and Gatto, P.: The natural subsidence in the
Lagoon of Venice, Italy, in: Proc. 3rd International Symp. on Land Subsidence,
IAHS Publ., 151, 777–784, edited by: Johnson, A.I., Carbognin, L., and
Ubertini, L., Venice, Italy, 1984.Brambati, A., Carbognin, L., Quaia, T., Teatini, P., and Tosi, L.: The Lagoon
of Venice: geological setting, evolution and land subsidence, Episodes, 26,
264–268, 10.18814/epiiugs/2003/v26i3/020, 2003.Brandt, P., Rubino, A., Sein, D. V., Baschek, B., Izquierdo, A., and Backhaus,
J. O.: Sea level variations in the western Mediterranean studied by a
numerical tidal model of the Strait of Gibraltar, J. Phys. Oceanogr., 34, 433–443,
10.1175/1520-0485(2004)034<0433:SLVITW>2.0.CO;2, 2004.Calafat, F. M., Chambers, D. P., and Tsimplis, M. N.: Mechanisms of decadal sea
level variability in the eastern North Atlantic and the Mediterranean Sea,
J. Geophys. Res.-Oceans, 117, C09022, 10.1029/2012JC008285,
2012.Camuffo, D. and Sturaro, G.: Use of proxy-documentary and instrumental data to
assess the risk factors leading to sea flooding in Venice, Glob. Planet Change, 40, 93–103, 10.1016/S0921-8181(03)00100-0, 2004.
Carbognin, L., Gatto, P., Mozzi, G., Gambolati, G., and Ricceri, G.: New trend
in the subsidence of Venice, in: Proceedings of the Second International Symposium on Land Subsidence, held at Anaheim, California, 13–17 December 1976, IAHS
Publication, vol. 121, 65–81, edited by: Rodda, J. C., 1976.
Carbognin, L., Marabini, F., and Tosi, L.: Land subsidence and degradation of
the Venice littoral zone, Italy, in: Fifth International Symposium on Land
Subsidence-FISOLS 1995, The Hague, The Netherlands, October 1995a.
Carbognin, L., Tosi, L., and Teatini, P.: Analysis of actual land subsidence in
Venice and its hinterland (Italy), Land subsidence, 129–137, edited by: Barends, F. B. J., Brouwer, F. J. J., and Schroder, F. H., Balkema, Rotterdam,
1995b.Carbognin, L., Teatini, P., and Tosi, L.: Eustacy and land subsidence in the
Venice Lagoon at the beginning of the new millennium, J. Marine Syst., 51, 345–353, 10.1016/j.jmarsys.2004.05.021, 2004.Carbognin, L., Teatini, P., and Tosi, L.: The impact of relative sea level rise
on the Northern Adriatic Sea coast, Italy, WIT Trans. Ecol. Envir., 127, 137–148, 10.2495/RAV090121, 2009.Carbognin, L., Teatini, P., Tomasin, A., and Tosi, L.: Global change and
relative sea level rise at Venice: what impact in term of flooding, Clim. Dynam., 35, 1039–1047, 10.1007/s00382-009-0617-5, 2010.Carillo, A., Sannino, G., Artale, V., Ruti, P., Calmanti, S., and Dell'Aquila,
A.: Steric sea level rise over the Mediterranean Sea: present climate and
scenario simulations, Clim. Dynam., 39, 2167–2184,
10.1007/s00382-012-1369-1, 2012.Carminati, E. and Di Donato, G.: Separating natural and anthropogenic vertical
movements in fast subsiding areas: the Po plain (N. Italy) case, Geophys. Res. Lett., 26, 2291–2294, 10.1029/1999GL900518, 1999.Carminati, E., Doglioni, C., and Scrocca, D.: Apennines subduction-related
subsidence of Venice (Italy), Geophys. Res. Lett., 30, 1717, 10.1029/2003GL017001, 2003.Carminati, E., Enzi, S., and Camuffo, D.: A study on the effects of seismicity
on subsidence in foreland basins: An application to the Venice area, Global Planet. Change, 55, 237–250, 10.1016/j.gloplacha.2006.03.003,
2007.Carrère, L. and Lyard, F.: Modeling the barotropic response of the global
ocean to atmospheric wind and pressure forcing-comparisons with observations,
Geophys. Res. Lett., 30, 8–8, 10.1029/2002GL016473, 2003.Cazenave, A., Hamlington, B., Horwath, M., Barletta, V. R., Benveniste, J.,
Chambers, D., Döll, P., Hogg, A. E., and Legeais, J. F., Merrifield, M., Meyssignac, B., Mitchum, G., Nerem, S., Pail, R., Palanisamy, H., Paul, F., von Schuckmann, K., and Thompson, P.: Observational requirements for long-term monitoring of the global
mean sea level and its components over the altimetry era, Front. Mar. Sci., 6, 582, 10.3389/fmars.2019.00582, 2019.Chambers, D. P.: Evaluation of empirical mode decomposition for quantifying multi-decadal variations and acceleration in sea level records, Nonlin. Processes Geophys., 22, 157–166, 10.5194/npg-22-157-2015, 2015.Chen, C., Gao, G., Zhang, Y., Beardsley, R. C., Lai, Z., Qi, J., and Lin, H.:
Circulation in the Arctic Ocean: Results from a high-resolution coupled
ice-sea nested Global-FVCOM and Arctic-FVCOM system, Prog. Oceanogr., 141, 60–80, 10.1016/j.pocean.2015.12.002, 2016.Church, J. A., Clark, P. U., Cazenave, A., Gregory, J. M., Jevrejeva, S.,
Levermann, A., Merrifield, M. A., Milne, G. A., Nerem, R. S., Nunn, P. D., Payne, A. J., Pfeffer, W. T., Stammer, D., and Unnikrishnan, A. S.: Sea level change, Tech. rep., Climate Change 2013: The Physical
Science Basis. Contribution of Working Group I to the Fifth Assessment Report
of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United
Kingdom and New York, NY, USA,
available at: http://drs.nio.org/drs/handle/2264/4605 (last access: 5 August 2021), 2013.Cipollini, P., Vignudelli, S., Lyard, F., and Roblou, L.: 15 years of altimetry
at various scales over the Mediterranean, 295–306, in: Remote Sensing of the European Seas, edited by: Barale, V. and Gade, M., Springer, Dordrecht, 10.1007/978-1-4020-6772-3_22, 2008.Cipollini, P., Passaro, M., and Vignudelli, S.: Validation of Total Water Level
Envelope data from the eSurge processor for Envisat in the North Adriatic
Sea, available at: http://www.storm-surge.info/project-docs (last access: 5 August 2021), 2013.Cipollini, P., Benveniste, J., Birol, F., Fernandes, M., Obligis, E., Passaro,
M., Strub, P., Valladeau, G., Vignudelli, S., and Wlikin, J.: Satellite
altimetry in coastal regions, CRC Press, Boca Raton, 10.1201/9781315151779, 2017.
CNR: Relazione sul pozzo Venezia 1-CNR. 1a fase: Operazioni di cantiere e
analisi delle carote. Laboratorio Dinamica Grandi Masse, Consiglio Nazionale
delle Ricerche, Venice, 1971.Criado-Aldeanueva, F., Soto-Navarro, F. J., and García-Lafuente, J.:
Large-scale atmospheric forcing influencing the long-term variability of
Mediterranean heat and freshwater budgets: climatic indices, J. Hydrometeorol., 15, 650–663, 10.1175/JHM-D-13-04.1, 2014.Cuffaro, M., Riguzzi, F., Scrocca, D., Antonioli, F., Carminati, E., Livani,
M., and Doglioni, C.: On the geodynamics of the northern Adriatic plate,
Rend. Lincei, 21, 253–279, 10.1007/s12210-010-0098-9, 2010.Cusinato, E., Zanchettin, D., Sannino, G., and Rubino, A.: Mediterranean
Thermohaline Response to Large-Scale Winter Atmospheric Forcing in a
High-Resolution Ocean Model Simulation, PApGe, 175, 4083–4110,
10.1007/s00024-018-1859-0, 2018.Dangendorf, S., Marcos, M., Wöppelmann, G., Conrad, C. P., Frederikse, T.,
and Riva, R.: Reassessment of 20th century global mean sea level rise,
P. Natl. Acad. Sci. USA, 114, 5946–5951,
10.1073/pnas.1616007114, 2017.De Biasio, F., Baldin, G., and Vignudelli, S.: Revisiting Vertical Land Motion
and Sea Level Trends in the Northeastern Adriatic Sea Using Satellite
Altimetry and Tide Gauge Data, Journal of Marine Science and Engineering, 8,
949, 10.3390/jmse8110949, 2020.DeConto, R. M. and Pollard, D.: Contribution of Antarctica to past and future
sea-level rise, Nature, 531, 591–597, 10.1038/nature17145, 2016.
Dodge, Y., Cox, D., Commenges, D., Davison, A., Solomon, P., and Wilson, S. (Eds.): International Statistical Institute, The Oxford dictionary of
statistical terms, Oxford University Press on demand, New York, ISBN 0198509944, 506 pp., 2003.Dorigo, L.: Le osservazioni mareografiche in Laguna di Venezia, Ist. Veneto
Sc., Lett. Arti, Commissione di studio dei provvedimenti per la conservazione
e difesa della Laguna e della città di Venezia, Rapporti preliminari,
available at: https://www.beic.it (last access: 5 August 2021), 1961.Edwards, T. L., Brandon, M. A., Durand, G., Edwards, N. R., Golledge, N. R.,
Holden, P. B., Nias, I. J., Payne, A. J., Ritz, C., and Wernecke, A.:
Revisiting Antarctic ice loss due to marine ice-cliff instability, Nature,
566, 58–64, 10.1038/s41586-019-0901-4, 2019.Environment Agency: Policy paper: Thames Estuary 2100: 10-Year Review
monitoring key findings, Tech. rep., Environment Agency, UK,
available at: https://www.gov.uk/government/publications/thames-estuary-2100-te2100/thames-estuary-2100-key-findings-from-the-monitoring-review, last access: 23 March 2021.Erkens, G., Bucx, T., Dam, R., De Lange, G., and Lambert, J.: Sinking coastal
cities, P. Int. Ass. Hydrol. Sci., 372, 189–198, 10.5194/piahs-372-189-2015, 2015.ESA CCI: Sea Level Bridging Phase Part 2: Validation. Technical Deliverable v
2.0 2017, available at: https://climate.esa.int/media/documents/D3_SL_cci_BP_Final_Report_Part2_Validation_20190722.pdf (last access: 5 August 2021), 2019.Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958, 10.5194/gmd-9-1937-2016, 2016.Ezer, T., Haigh, I. D., and Woodworth, P. L.: Nonlinear sea-level trends and
long-term variability on western European coasts, J. Coast. Res., 32, 744–755, 10.2112/JCOASTRES-D-15-00165.1, 2016.Farrell, W. and Clark, J. A.: On postglacial sea level, Geophys. J. Int., 46, 647–667,
10.1111/j.1365-246X.1976.tb01252.x, 1976.
Favero, V., Alberotanza, L., and Serandrei Barbero, R.: Aspetti paleoecologici,
sedimentologici e geochimici dei sedimenti attraversati dal pozzo VE 1 bis
CNR Laboratorio per lo Studio della Dinamica delle Grandi Masse, Tech. rep.,
CNR Technical Report 63, Venice, 1973.Fenoglio-Marc, L., Braitenberg, C., and Tunini, L.: Sea level variability and
trends in the Adriatic Sea in 1993–2008 from tide gauges and satellite
altimetry, Phys. Chem Earth PT A/B/C, 40, 47–58,
10.1016/j.pce.2011.05.014, 2012.Fernandes, M. J. and Lázaro, C.: GPD+ wet tropospheric corrections for
CryoSat-2 and GFO altimetry missions, Remote Sens., 8, 851,
10.3390/rs8100851, 2016.Ferranti, L., Antonioli, F., Mauz, B., Amorosi, A., Dai Pra, G., Mastronuzzi,
G., Monaco, C., Orrù, P., Pappalardo, M., Radtke, U., Renda, P., Romano, P., Sansò, P., and Verrubbi, V.: Markers of
the last interglacial sea-level high stand along the coast of Italy: tectonic
implications, Quatern. Int., 145, 30–54,
10.1016/j.quaint.2005.07.009, 2006.Ferrarin, C., Bellafiore, D., Sannino, G., Bajo, M., and Umgiesser, G.: Tidal
dynamics in the inter-connected Mediterranean, Marmara, Black and Azov seas,
Prog. Oceanogr., 161, 102–115, 10.1016/j.pocean.2018.02.006,
2018.
Flemming, N. C.: Predictions of relative coastal sea-level change in the
Mediterranean based on archaeological, historical and tide-gauge data,
chap. 8, 247–281, in: Climate Change and the Mediterrenean, edited by: Jeftic, L., Millman, J. D., and Sestini, G., 247–281, Edward Arnold, London, 1992.Fontana, A., Vinci, G., Tasca, G., Mozzi, P., Vacchi, M., Bivi, G., Salvador,
S., Rossato, S., Antonioli, F., Asioli, A., Bresolin, M., Di Mario, F., and Hajdas, I.: Lagoonal settlements and
relative sea level during Bronze Age in Northern Adriatic: Geoarchaeological
evidence and paleogeographic constraints, Quaternary Int., 439,
17–36, 10.1016/j.quaint.2016.12.038, 2017.Frederikse, T., Landerer, F., Caron, L., Adhikari, S., Parkes, D., Humphrey,
V. W., Dangendorf, S., Hogarth, P., Zanna, L., Cheng, L., and Wu, Y.-H.: The causes
of sea-level rise since 1900, Nature, 584, 393–397,
10.1038/s41586-020-2591-3, 2020.
Fu, L.-L. and Cazenave, A.: Satellite altimetry and earth sciences: a handbook
of techniques and applications, Academic Press, London, 2001.Fukumori, I., Menemenlis, D., and Lee, T.: A near-uniform basin-wide sea level
fluctuation of the Mediterranean Sea, J. Phys. Oceanogr., 37,
338–358, 10.1175/JPO3016.1, 2007.Galassi, G. and Spada, G.: Sea-level rise in the Mediterranean Sea by 2050:
Roles of terrestrial ice melt, steric effects and glacial isostatic
adjustment, Global Planet. Change, 123, 55–66,
10.1016/j.gloplacha.2014.10.007, 2014.
Galloway, D. and Riley, F. S.: San Joaquin Valley, California, Land subsidence in the United States, Circular 1182, US Geological Survey, 1182, 23–34, 1999.Galloway, D. L., Erkens, G., Kuniansky, E. L., and Rowland, J. C.: Preface:
Land subsidence processes, Hydrogeol. J., 24, 547–550,
10.1007/s10040-016-1386-y, 2016.
Gambolati, G. and Gatto, P.: Simulazione della subsidenza di Venezia, Venezia
ei problemi dell'ambiente, Il Mulino, Bologna, 299–360, 1975.Gambolati, G., Gatto, P., and Freeze, R. A.: Mathematical simulation of the
subsidence of Venice: 2. Results, Water Resour. Res., 10, 563–577,
10.1029/WR010i003p00563, 1974.Gambolati, G., Teatini, P., and Ferronato, M.: Anthropogenic land subsidence, 2443–2459, in: Encyclopedia of Hydrological Sciences, edited by: Anderson, M. G., John Wiley & Sons, New Jersey, 10.1002/0470848944.hsa164b,
2006.García Lafuente, J., Bruque Pozas, E., Sánchez Garrido, J. C., Sannino,
G., and Sammartino, S.: The interface mixing layer and the tidal dynamics at
the eastern part of the Strait of Gibraltar, J. Marine Syst.,
117–118, 31–42, 10.1016/j.jmarsys.2013.02.014, 2013.Gatto, P. and Carbognin, L.: The Lagoon of Venice: natural environmental trend
and man-induced modification/La Lagune de Venise: l'évolution naturelle
et les modifications humaines, Hydrolog. Sci. J., 26, 379–391,
10.1080/02626668109490902, 1981.Gómez-Enri, J., González, C., Passaro, M., Vignudelli, S., Álvarez,
O., Cipollini, P., Mañanes, R., Bruno, M., López-Carmona, M., and
Izquierdo, A.: Wind-induced cross-strait sea level variability in the Strait
of Gibraltar from coastal altimetry and in-situ measurements, Remote Sens. Environ., 221, 596–608, 10.1016/j.rse.2018.11.042, 2019.Gomis, D., Ruiz, S., Sotillo, M. G., Álvarez-Fanjul, E., and Terradas, J.:
Low frequency Mediterranean sea level variability: the contribution of
atmospheric pressure and wind, Global Planet. Change, 63, 215–229,
10.1016/j.gloplacha.2008.06.005, 2008.Good, S. A., Martin, M. J., and Rayner, N. A.: EN4: Quality controlled ocean
temperature and salinity profiles and monthly objective analyses with
uncertainty estimates, J. Geophys. Res.-Oceans, 118,
6704–6716, 10.1002/2013JC009067, 2013.Gregory, J. M., Griffies, S. M., Hughes, C. W., Lowe, J. A., Church, J. A.,
Fukimori, I., Gomez, N., Kopp, R. E., Landerer, F., Le Cozannet, G., Ponte, R. M., Stammer, D., Tamisiea, M. E., and van de Wal, R. S. W.:
Concepts and terminology for sea level: mean, variability and change, both
local and global, Surv. Geophys., 40, 1251–1289,
10.1007/s10712-019-09525-z, 2019.Grinsted, A., Moore, J. C., and Jevrejeva, S.: Application of the cross wavelet transform and wavelet coherence to geophysical time series, Nonlin. Processes Geophys., 11, 561–566, 10.5194/npg-11-561-2004, 2004.Haasnoot, M., Kwadijk, J., Van Alphen, J., Le Bars, D., Van Den Hurk, B.,
Diermanse, F., Van Der Spek, A., Essink, G. O., Delsman, J., and Mens, M.:
Adaptation to uncertain sea-level rise; how uncertainty in Antarctic
mass-loss impacts the coastal adaptation strategy of the Netherlands,
Environ. Res. Lett., 15, 034007, 10.1088/1748-9326/ab666c,
2020.Han, W., Stammer, D., Thompson, P., Ezer, T., Palanisamy, H., Zhang, X.,
Domingues, C. M., Zhang, L., and Yuan, D.: Impacts of basin-scale climate
modes on coastal sea level: a review, Surv. Geophys., 40, 1493–1541,
10.1007/s10712-019-09562-8, 2019.Hay, C. C., Morrow, E., Kopp, R. E., and Mitrovica, J. X.: Probabilistic
reanalysis of twentieth-century sea-level rise, Nature, 517, 481–484,
10.1038/nature14093, 2015.Hinkel, J., Church, J. A., Gregory, J. M., Lambert, E., Le Cozannet, G., Lowe,
J., McInnes, K. L., Nicholls, R. J., van der Pol, T. D., and van de Wal, R.:
Meeting user needs for sea level rise information: A decision analysis
perspective, Earth's Future, 7, 320–337, 10.1029/2018EF001071, 2019.Holt, J., Hyder, P., Ashworth, M., Harle, J., Hewitt, H. T., Liu, H., New, A. L., Pickles, S., Porter, A., Popova, E., Allen, J. I., Siddorn, J., and Wood, R.: Prospects for improving the representation of coastal and shelf seas in global ocean models, Geosci. Model Dev., 10, 499–523, 10.5194/gmd-10-499-2017, 2017.Ishii, M., Fukuda, Y., Hirahara, S., Yasui, S., Suzuki, T., and Sato, K.:
Accuracy of global upper ocean heat content estimation expected from present
observational data sets, Sola, 13, 163–167, 10.2151/sola.2017-030,
2017.
ISPRA: Manuale di mareografia e linee guida per i processi di validazione dei
dati mareografici, Manuali e linee guida 77/2012, ISPRA – Istituto Superiore per la Protezione e la Ricerca Ambientale, Roma, ISBN 978-88-448-0532-6,
2012.Jones, P., Jonsson, T., and Wheeler, D.: Extension to the North Atlantic
Oscillation using early instrumental pressure observations from Gibraltar and
South-West Iceland, Int. J. Climatol., 17, 1433–1450,
10.1002/(SICI)1097-0088(19971115)17:13<1433::AID-JOC203>3.0.CO;2-P,
1997.Jordà, G.: Detection time for global and regional sea level trends and
accelerations, J. Geophys. Res.-Oceans, 119, 7164–7174,
10.1002/2014JC010005, 2014.Jordà, G. and Gomis, D.: On the interpretation of the steric and mass
components of sea level variability: The case of the Mediterranean basin,
J. Geophys. Res.-Oceans, 118, 953–963,
10.1002/jgrc.20060, 2013.Jordà, G., Gomis, D., and Álvarez-Fanjul, E.: The VANI2-ERA hindcast of
sea-level residuals: atmospheric forcing of sea-level variability in the
Mediterranean Sea (1958–2008), Sci. Mar., 76, 133–146,
10.3989/scimar.03612.19C, 2012a.Jordà, G., Gomis, D., Álvarez-Fanjul, E., and Somot, S.: Atmospheric
contribution to Mediterranean and nearby Atlantic sea level variability under
different climate change scenarios, Global Planet. Change, 80,
198–214, 10.1016/j.gloplacha.2011.10.013, 2012b.Josey, S. A., Somot, S., and Tsimplis, M.: Impacts of atmospheric modes of
variability on Mediterranean Sea surface heat exchange, J. Geophys. Res.-Oceans, 116, C02032, 10.1029/2010JC006685, 2011.Kent, D. V., Rio, D., Massari, F., Kukla, G., and Lanci, L.: Emergence of
Venice during the Pleistocene, Quat. Sci. Rev., 21, 1719–1727,
10.1016/S0277-3791(01)00153-6, 2002.King, M. A., Altamimi, Z., Boehm, J., Bos, M., Dach, R., Elosegui, P., Fund,
F., Hernández-Pajares, M., Lavallee, D., Cerveira, P. J. M., Penna, N., Riva, R. E. M., Steigenberger, P., van Dam, T., Vittuari, L., Williams, S., and Willis, P.:
Improved constraints on models of glacial isostatic adjustment: a review of
the contribution of ground-based geodetic observations, Surv. Geophys., 31, 465–507, 10.1007/s10712-010-9100-4, 2010.Kopp, R. E., Horton, R. M., Little, C. M., Mitrovica, J. X., Oppenheimer, M.,
Rasmussen, D., Strauss, B. H., and Tebaldi, C.: Probabilistic 21st and 22nd
century sea-level projections at a global network of tide-gauge sites,
Earth's future, 2, 383–406, 10.1002/2014EF000239, 2014.Kopp, R. E., DeConto, R. M., Bader, D. A., Hay, C. C., Horton, R. M., Kulp, S.,
Oppenheimer, M., Pollard, D., and Strauss, B. H.: Evolving understanding of
Antarctic ice-sheet physics and ambiguity in probabilistic sea-level
projections, Earth's Future, 5, 1217–1233, 10.1002/2017EF000663, 2017.Lambeck, K., Antonioli, F., Purcell, A., and Silenzi, S.: Sea-level change
along the Italian coast for the past 10,000 yr, Quat. Sci. Rev.,
23, 1567–1598, 10.1016/j.quascirev.2004.02.009, 2004.Lambeck, K., Antonioli, F., Anzidei, M., Ferranti, L., Leoni, G., Scicchitano,
G., and Silenzi, S.: Sea level change along the Italian coast during the
Holocene and projections for the future, Quaternary Int., 232,
250–257, 10.1016/j.quaint.2010.04.026, 2011.Landerer, F. W. and Volkov, D. L.: The anatomy of recent large sea level
fluctuations in the Mediterranean Sea, Geophys. Res. Lett., 40,
553–557, 10.1002/grl.50140, 2013.Landerer, F. W., Gleckler, P. J., and Lee, T.: Evaluation of CMIP5 dynamic sea
surface height multi-model simulations against satellite observations,
Clim. Dynam., 43, 1271–1283, 10.1007/s00382-013-1939-x, 2014.Legeais, J.-F., Ablain, M., Zawadzki, L., Zuo, H., Johannessen, J. A., Scharffenberg, M. G., Fenoglio-Marc, L., Fernandes, M. J., Andersen, O. B., Rudenko, S., Cipollini, P., Quartly, G. D., Passaro, M., Cazenave, A., and Benveniste, J.: An improved and homogeneous altimeter sea level record from the ESA Climate Change Initiative, Earth Syst. Sci. Data, 10, 281–301, 10.5194/essd-10-281-2018, 2018.
Lempert, R.: Robust Decision Making (RDM) BT–Decision Making under Deep
Uncertainty: From Theory to Practice, edited by:
Marchau, V. A. W. J., Walker, W. E., Bloemen, P. J. T. M., and Popper, S. W., Springer, Berlin, 23–51, 2019.
Lezziero, A.: Indagini paleoambientali nel sottosuolo dell’isola di Burano,
Quaderni di Insula, Venezia, 10, 63–70, 2002.
Lionello, P.: Extreme storm surges in the gulf of venice: Present and future
climate, Flooding and Environmental Challenges for Venice and its Lagoon:
State of Knowledge, in: Flooding and Environmental Challenges for Venice and its Lagoon: State of Knowledge, edited by: Fletcher, C. A. and Spencer, T., 59–69, Cambridge Univ. Press, New York, 2005.Lionello, P., Nicholls, R. J., Umgiesser, G., and Zanchettin, D.: Venice flooding and sea level: past evolution, present issues, and future projections (introduction to the special issue), Nat. Hazards Earth Syst. Sci., 21, 2633–2641, 10.5194/nhess-21-2633-2021, 2021a.Lionello, P., Barriopedro, D., Ferrarin, C., Nicholls, R. J., Orlić, M., Raicich, F., Reale, M., Umgiesser, G., Vousdoukas, M., and Zanchettin, D.: Extreme floods of Venice: characteristics, dynamics, past and future evolution (review article), Nat. Hazards Earth Syst. Sci., 21, 2705–2731, 10.5194/nhess-21-2705-2021, 2021b.Magrini, G. P., De Marchi, L., and Gnesotto, T.: Impianti mareografici eseguiti,
Ricerche Lagunari, available at: https://www.beic.it (last access: 5 August 2021), 1908.Marcos, M. and Tsimplis, M. N.: Coastal sea level trends in Southern Europe,
Geophys. J. Int., 175, 70–82,
10.1111/j.1365-246X.2008.03892.x, 2008.Marcos, M., Jordà, G., and Le Cozannet, G.: Sea level rise and its impacts
on the Mediterranean, The Mediterranean region under climate change. A
scientific update, IRD Éditions, Marseille, 265–275,
10.4000/books.irdeditions.22908, 2016.Mariotti, A. and Dell'Aquila, A.: Decadal climate variability in the
Mediterranean region: roles of large-scale forcings and regional processes,
Clim. Dynam., 38, 1129–1145, 10.1007/s00382-011-1056-7, 2012.Marshall, J., Kushnir, Y., Battisti, D., Chang, P., Czaja, A., Dickson, R.,
Hurrell, J., McCartney, M., Saravanan, R., and Visbeck, M.: North Atlantic
climate variability: phenomena, impacts and mechanisms, Int. J. Climatol., 21,
1863–1898, 10.1002/joc.693, 2001.Marti, F., Cazenave, A., Birol, F., Passaro, M., Léger, F., Niño, F.,
Almar, R., Benveniste, J., and Legeais, J. F.: Altimetry-based sea level
trends along the coasts of western Africa, Adv. Space Res., 68, 504–522, 10.1016/j.asr.2019.05.033, 2019.Martínez-Asensio, A., Marcos, M., Tsimplis, M. N., Gomis, D., Josey, S.,
and Jordà, G.: Impact of the atmospheric climate modes on Mediterranean
sea level variability, Global Planet. Change, 118, 1–15,
10.1016/j.gloplacha.2014.03.007, 2014.Martínez-Asensio, A., Tsimplis, M. N., and Calafat, F. M.: Decadal
variability of European sea level extremes in relation to the solar activity,
Geophys. Res. Lett., 43, 11–744, 10.1002/2016GL071355, 2016.Maslova, V. N., Voskresenskaya, E. N., and Lubkov, A. S.: Multidecadal change
of winter cyclonic activity in the Mediterranean associated with AMO and
PDO, Terr. Atmos. Ocean. Sci., 28, 965–977, 10.3319/TAO.2017.04.23.01, 2017.Massari, F., Rio, D., Barbero, R. S., Asioli, A., Capraro, L., Fornaciari, E.,
and Vergerio, P.: The environment of Venice area in the past two million
years, Palaeogeogr. Palaeocl., 202, 273–308,
10.1016/S0031-0182(03)00640-0, 2004.Melini, D. and Spada, G.: Some remarks on Glacial Isostatic Adjustment
modelling uncertainties, Geophys. J. Int., 218, 401–413,
10.1093/gji/ggz158, 2019.Menemenlis, D., Fukumori, I., and Lee, T.: Atlantic to Mediterranean sea level
difference driven by winds near Gibraltar Strait, J. Phys. Oceanogr., 37, 359–376, 10.1175/JPO3015.1, 2007.Meyssignac, B., Slangen, A. A., Melet, A., Church, J., Fettweis, X., Marzeion,
B., Agosta, C., Ligtenberg, S., Spada, G., Richter, K., Palmer, M. D., Roberts, C. D., and Champollion, N.: Evaluating
model simulations of twentieth-century sea-level rise. Part II: Regional
sea-level changes, J. Climate, 30, 8565–8593,
10.1175/JCLI-D-17-0112.1, 2017.Milne, G. A., Gehrels, W. R., Hughes, C. W., and Tamisiea, M. E.: Identifying
the causes of sea-level change, Nat. Geosci., 2, 471–478,
10.1038/ngeo544, 2009.Mitrovica, J. X., Gomez, N., and Clark, P. U.: The sea-level fingerprint of
West Antarctic collapse, Science, 323, 753–753,
10.1126/science.1166510, 2009.Naranjo, C., Garcia-Lafuente, J., Sannino, G., and Sanchez-Garrido, J.: How
much do tides affect the circulation of the Mediterranean Sea? From local
processes in the Strait of Gibraltar to basin-scale effects, Prog. Oceanogr., 127, 108–116, 10.1016/j.pocean.2014.06.005, 2014.Nicholls, R. J., Hanson, S. E., Lowe, J. A., Warrick, R. A., Lu, X., and Long,
A. J.: Sea-level scenarios for evaluating coastal impacts, WIRES Clim. Change, 5, 129–150,
10.1002/wcc.253, 2014.Nicholls, R. J., Hanson, S. E., Lowe, J. A., Slangen, A. B., Wahl, T., Hinkel,
J., and Long, A. J.: Integrating new sea-level scenarios into coastal risk
and adaptation assessments: An ongoing process, WIRES Clim. Change, 12, e706, 10.1002/wcc.706, 2021.
Oppenheimer, M., Glavovic, B., Hinkel, J., van de Wal, R., Magnan, A. K.,
Abd-Elgawad, A., Cai, R., Cifuentes-Jara, M., Deconto, R. M., Ghosh, T., Hay, J., Isla, F., Marzeion, B., Meyssignac, B., and Sebesvari, Z.: Sea level rise and implications for low lying islands, coasts and
communities, in: IPCC Special Report on the
Ocean and Cryosphere in a Changing Climate, edited by: Roberts, D. C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., Mintenbeck, K., Alegría, A., Nicolai, M., Okem, A., Petzold, J., Rama, B., and Weyer, N. M., Cambridge University Press,
Cambridge, UK, 2019.Palma, M., Iacono, R., Sannino, G., Bargagli, A., Carillo, A., Fekete, B. M.,
Lombardi, E., Napolitano, E., Pisacane, G., and Struglia, M. V.: Short-term,
linear, and non-linear local effects of the tides on the surface dynamics in
a new, high-resolution model of the Mediterranean Sea circulation, Ocean
Dynam., 70, 935–963, 10.1007/s10236-020-01364-6, 2020.Parras-Berrocal, I. M., Vazquez, R., Cabos, W., Sein, D., Mañanes, R., Perez-Sanz, J., and Izquierdo, A.: The climate change signal in the Mediterranean Sea in a regionally coupled atmosphere–ocean model, Ocean Sci., 16, 743–765, 10.5194/os-16-743-2020, 2020.Passaro, M., Cipollini, P., Vignudelli, S., Quartly, G. D., and Snaith, H. M.:
ALES: A multi-mission adaptive subwaveform retracker for coastal and open
ocean altimetry, Remote Sens. Environ., 145, 173–189,
10.1016/j.rse.2014.02.008, 2014.Peltier, W. R.: Global glacial isostasy and the surface of the ice-age Earth:
the ICE-5G (VM2) model and GRACE, Annu. Rev. Earth Planet. Sci., 32,
111–149, 10.1146/annurev.earth.32.082503.144359, 2004.Piecuch, C. G., Dangendorf, S., Gawarkiewicz, G. G., Little, C. M., Ponte,
R. M., and Yang, J.: How is New England coastal sea level related to the
Atlantic meridional overturning circulation at 26∘ N?, Geophys. Res. Lett., 46, 5351–5360, 10.1029/2019GL083073, 2019.
Pirazzoli, P.: Sea level changes. The last 2000 years, John Wiley & Sons, Chichester, 211 pp., 1996.Pisacane, G., Artale, V., Calmanti, S., and Rupolo, V.: Decadal oscillations in
the Mediterranean Sea: a result of the overturning circulation variability in
the eastern basin?, Clim. Res., 31, 257–271, 10.3354/cr031257,
2006.Polli, S.: Il graduale aumento del livello del mare lungo le coste italiane,
Geofisica pura e applicata, 25, 123–129, 10.1007/BF02014061, 1953.Polyakov, I. V., Pnyushkov, A. V., and Timokhov, L. A.: Warming of the
Intermediate Atlantic Water of the Arctic Ocean in the 2000s, J. Climate, 25, 8362–8370, 10.1175/JCLI-D-12-00266.1, 2012.Ponte, R. M., Carson, M., Cirano, M., Domingues, C. M., Jevrejeva, S., Marcos,
M., Mitchum, G., van de Wal, R. S. W., Woodworth, P. L., Ablain, M., Ardhuin,
F., Ballu, V., Becker, M., Benveniste, J., Birol, F., Bradshaw, E., Cazenave,
A., De Mey-Frémaux, P., Durand, F., Ezer, T., Fu, L.-L., Fukumori, I.,
Gordon, K., Gravelle, M., Griffies, S. M., Han, W., Hibbert, A., Hughes,
C. W., Idier, D., Kourafalou, V. H., Little, C. M., Matthews, A., Melet, A.,
Merrifield, M., Meyssignac, B., Minobe, S., Penduff, T., Picot, N., Piecuch,
C., Ray, R. D., Rickards, L., Santamaría-Gómez, A., Stammer, D., Staneva,
J., Testut, L., Thompson, K., Thompson, P., Vignudelli, S., Williams, J.,
Williams, S. D. P., Wöppelmann, G., Zanna, L., and Zhang, X.: Towards
comprehensive observing and modeling systems for monitoring and predicting
regional to coastal sea level, Front. Mar. Sci., 6, 437, 10.3389/fmars.2019.00437, 2019.Quartly, G. D., Legeais, J.-F., Ablain, M., Zawadzki, L., Fernandes, M. J., Rudenko, S., Carrère, L., García, P. N., Cipollini, P., Andersen, O. B., Poisson, J.-C., Mbajon Njiche, S., Cazenave, A., and Benveniste, J.: A new phase in the production of quality-controlled sea level data, Earth Syst. Sci. Data, 9, 557–572, 10.5194/essd-9-557-2017, 2017.Ranger, N., Reeder, T., and Lowe, J.: Addressing “deep” uncertainty over
long-term climate in major infrastructure projects: four innovations of the
Thames Estuary 2100 Project, EURO Journal on Decision Processes, 1, 233–262,
10.1007/s40070-013-0014-5, 2013.Rocco, F. V.: Sea level trends in the Mediterranean from tide gauges and
satellite altimetry, Ph.D. thesis, alma Mater Studiorum, University of
Bologna, Bologna, Italy, available at: https://amslaurea.unibo.it/10172/1/Rocco_Francois_Vincent_tesi.pdf, (last access: 5 August 2021), 2015.Roy, K. and Peltier, W. R.: Relative sea level in the Western Mediterranean
basin: a regional test of the ICE-7G_NA (VM7) model and a constraint on late
Holocene Antarctic deglaciation, Quat. Sci. Rev., 183, 76–87,
10.1016/j.quascirev.2017.12.021, 2018.Rubino, A., Zanchettin, D., Androsov, A., and Voltzinger, N. E.: Tidal Records
as Liquid Climate Archives for Large-Scale Interior Mediterranean
Variability, Sci. Rep., 8, 1–8, 10.1038/s41598-018-30930-8,
2018.
Rusconi, A.: Il Comune Marino a Venezia: ricerche e ipotesi sulle sue
variazioni altimetriche e sui fenomeni naturali che le determinano, Ufficio Idrografico Magistrato Acque, Venice, Italy, Pubbl., 157, 39, 1983.
Salvioni, G.: I movimenti del suolo nell'Italia centro-settentrionale,
Bollettino di geodesia e scienze affini, 16, 325–366, 1957.Sánchez, L., Völksen, C., Sokolov, A., Arenz, H., and Seitz, F.: Present-day surface deformation of the Alpine region inferred from geodetic techniques, Earth Syst. Sci. Data, 10, 1503–1526, 10.5194/essd-10-1503-2018, 2018.Sannino, G., Pratt, L., and Carillo, A.: Hydraulic criticality of the exchange
flow through the Strait of Gibraltar, J. Phys. Oceanogr., 39,
2779–2799, 10.1175/2009JPO4075.1, 2009.Sannino, G., Carillo, A., Pisacane, G., and Naranjo, C.: On the relevance of
tidal forcing in modelling the Mediterranean thermohaline circulation,
Prog. Oceanogr., 134, 304–329, 10.1016/j.pocean.2015.03.002,
2015.Santamaría-Gómez, A., Gravelle, M., Dangendorf, S., Marcos, M., Spada,
G., and Wöppelmann, G.: Uncertainty of the 20th century sea-level rise
due to vertical land motion errors, Earth Planet. Sci. Lett., 473,
24–32, 10.1016/j.epsl.2017.05.038, 2017.Scafetta, N.: Multi-scale dynamical analysis (MSDA) of sea level records versus
PDO, AMO, and NAO indexes, Clim. Dynam., 43, 175–192,
10.1007/s00382-013-1771-3, 2014.Scarascia, L. and Lionello, P.: Global and regional factors contributing to the
past and future sea level rise in the Northern Adriatic Sea, Global Planet. Change, 106, 51–63, 10.1016/j.gloplacha.2013.03.004, 2013.
Scharroo, R., Leuliette, E., Lillibridge, J., Byrne, D., Naeije, M., and
Mitchum, G.: RADS: Consistent Multi-Mission Products, ESASP, in Proceedings
of the Symposium on 20 Years of Progress in Radar Altimetry, 24–29 September
2012, Special Publication SP-710, p. 69, European Space Agency, Venice,
Italy, 2013.Sein, D. V., Mikolajewicz, U., Gröger, M., Fast, I., Cabos, W., Pinto,
J. G., Hagemann, S., Semmler, T., Izquierdo, A., and Jacob, D.: Regionally
coupled atmosphere-ocean-sea ice-marine biogeochemistry model ROM: 1.
Description and validation, J. Adv. Model. Earth Sy., 7,
268–304, 10.1002/2014MS000357, 2015.
Serandrei Barbero, R., Lezziero, A., Albani, A., and Zoppi, U.: Depositi
tardo-pleistocenici ed olocenici nel sottosuolo veneziano: paleoambienti e
cronologia, Il Quaternario, 14, 9–22, 2001.Simoncelli, S., Fratianni, C., Pinardi, N., Grandi, A., Drudi, M., and Oddo,
P.: Mediterranean Sea physical reanalysis (MEDREA 1987–2017) (Version 1),
Copernicus Monitoring Environment Marine Service (CMEMS) [data set],
10.25423/medsea_reanalysis_phys_006_004, 2014.Slangen, A., Katsman, C., Van de Wal, R., Vermeersen, L., and Riva, R.: Towards
regional projections of twenty-first century sea-level change based on IPCC
SRES scenarios, Clim. Dynam., 38, 1191–1209,
10.1007/s00382-011-1057-6, 2012.Slangen, A., Carson, M., Katsman, C., Van de Wal, R., Köhl, A., Vermeersen,
L., and Stammer, D.: Projecting twenty-first century regional sea-level
changes, Clim. Change, 124, 317–332, 10.1007/s10584-014-1080-9,
2014.Slangen, A., Adloff, F., Jevrejeva, S., Leclercq, P., Marzeion, B., Wada, Y.,
and Winkelmann, R.: A review of recent updates of sea-level projections at
global and regional scales, Surv. Geophys., 38, 385–406,
10.1007/s10712-016-9374-2, 2017.Slangen, A. B., Church, J. A., Agosta, C., Fettweis, X., Marzeion, B., and
Richter, K.: Anthropogenic forcing dominates global mean sea-level rise since
1970, Nat. Clim. Change, 6, 701–705, 10.1038/nclimate2991, 2016.Slivinski, L. C., Compo, G. P., Whitaker, J. S., Sardeshmukh, P. D., Giese,
B. S., McColl, C., Allan, R., Yin, X., Vose, R., Titchner, H., Kennedy, J., Spencer, L. J., Ashcroft, L., Brönnimann, S., Brunet, M., Camuffo, D., Cornes, R., Cram, T. A., Crouthamel, R., Domínguez-Castro, F., Freeman, J. E., Gergis, J., Hawkins, E., Jones, P. D., Jourdain, S., Kaplan, A., Kubota, H., Le Blancq, F., Lee, T.-C., Lorrey, A., Luterbacher, J., Maugeri, M., Mock, C. J., Moore, G. W. K., Przybylak, R., Pudmenzky, C., Reason, C., Slonosky, V. C., Smith, C. A., Tinz, B., Trewin, B., Valente, M. A., Wang, X. L., Wilkinson, C., Wood, K., and Wyszyński, P.:
Towards a more reliable historical reanalysis: Improvements for version 3 of
the Twentieth Century Reanalysis system, Q. J. Roy. Meteor. Soc., 145, 2876–2908, 10.1002/qj.3598, 2019.Somot, S., Sevault, F., Déqué, M., and Crépon, M.: 21st century
climate change scenario for the Mediterranean using a coupled
atmosphere–ocean regional climate model, Global Planet. Change, 63,
112–126, 10.1016/j.gloplacha.2007.10.003, 2008.Spada, G.: Glacial isostatic adjustment and contemporary sea level rise: An
overview, Surv. Geophys., 38, 153–185,
10.1007/s10712-016-9379-x, 2017.Spada, G. and Melini, D.: On Some Properties of the Glacial Isostatic
Adjustment Fingerprints, Water, 11, 1844, 10.3390/w11091844, 2019.Spada, G., Stocchi, P., and Colleoni, F.: Glacio–isostatic adjustment in the
po plain and in the northern adriatic region, Pure Appl. Geophys.,
166, 1303–1318, 10.1007/s00024-004-0498-9, 2009.Spada, G., Bamber, J., and Hurkmans, R.: The gravitationally consistent
sea-level fingerprint of future terrestrial ice loss, Geophys. Res. Lett., 40, 482–486, 10.1029/2012GL053000, 2013.Stammer, D. and Cazenave, A.: Satellite altimetry over oceans and land
surfaces, CRC Press, Boca Raton, 10.1201/9781315151779, 2017.Stocchi, P. and Spada, G.: Influence of glacial isostatic adjustment upon
current sea level variations in the Mediterranean, Tectonophysics, 474,
56–68, 10.1016/j.tecto.2009.01.003, 2009.Stocchi, P., Spada, G., and Cianetti, S.: Isostatic rebound following the
Alpine deglaciation: impact on the sea level variations and vertical
movements in the Mediterranean region, Geophys. J. Int.,
162, 137–147, 10.1111/j.1365-246X.2005.02653.x, 2005.Tambroni, N. and Seminara, G.: Are inlets responsible for the morphological
degradation of Venice Lagoon?, J. Geophys. Res.-Earth, 111, F03013, 10.1029/2005JF000334, 2006.Taricco, C., Alessio, S., Rubinetti, S., Zanchettin, D., Cosoli, S.,
Gačić, M., Mancuso, S., and Rubino, A.: Marine Sediments Remotely
Unveil Long-Term Climatic Variability Over Northern Italy, Sci. Rep., 5, 12111, 10.1038/srep12111, 2015.Teatini, P., Ferronato, M., Gambolati, G., Bertoni, W., and Gonella, M.: A
century of land subsidence in Ravenna, Italy, Environ. Geol., 47,
831–846, 10.1007/s00254-004-1215-9, 2005.Teatini, P., Tosi, L., Strozzi, T., Carbognin, L., Cecconi, G., Rosselli, R.,
and Libardo, S.: Resolving land subsidence within the Venice Lagoon by
persistent scatterer SAR interferometry, Phys. Chem. Earth PT A/B/C, 40, 72–79, 10.1016/j.pce.2010.01.002, 2012.Thiéblemont, R., Le Cozannet, G., Toimil, A., Meyssignac, B., and Losada,
I. J.: Likely and High-End Impacts of Regional Sea-Level Rise on the
Shoreline Change of European Sandy Coasts Under a High Greenhouse Gas
Emissions Scenario, Water, 11, 2607, 10.3390/w11122607, 2019.
Torge, W. and Müller, J.: Geodesy, Walter de Gruyter, Berlin, Germany, 1980.Tosi, L., Teatini, P., Carbognin, L., and Brancolini, G.: Using high resolution
data to reveal depth-dependent mechanisms that drive land subsidence: The
Venice coast, Italy, Tectonophysics, 1, 271–284,
10.1016/j.tecto.2009.02.026, 2009.Tosi, L., Teatini, P., Strozzi, T., Carbognin, L., Brancolini, G., and
Rizzetto, F.: Ground surface dynamics in the northern Adriatic coastland over
the last two decades, Rend. Lincei, 21, 115–129,
10.1007/s12210-010-0084-2, 2010.Tosi, L., Teatini, P., and Strozzi, T.: Natural versus anthropogenic subsidence
of Venice, Sci. Rep., 3, 2710, 10.1038/srep02710, 2013.Tosi, L., Lio, C. D., Teatini, P., and Strozzi, T.: Land subsidence in coastal
environments: knowledge advance in the Venice Coastland by TerraSAR-X PSI,
Remote Sens., 10, 1191, 10.3390/rs10081191, 2018.Troccoli, A., Zambon, F., Hodges, K. I., and Marani, M.: Storm surge frequency
reduction in Venice under climate change, Clim. Change, 113, 1065–1079,
10.1007/s10584-011-0093-x, 2012.Tsimplis, M., Álvarez-Fanjul, E., Gomis, D., Fenoglio-Marc, L., and
Pérez, B.: Mediterranean sea level trends: atmospheric pressure and wind
contribution, Geophys. Res. Lett., 32, L20602, 10.1029/2005GL023867,
2005.Tsimplis, M., Shaw, A., Flather, R., and Woolf, D.: The influence of the North
Atlantic Oscillation on the sea-level around the northern European coasts
reconsidered: the thermosteric effects, Philos. T. Roy. Soc. A, 364,
845–856, 10.1098/rsta.2006.1740, 2006.Tsimplis, M., Calafat, F. M., Marcos, M., Jordá, G., Gomis, D.,
Fenoglio-Marc, L., Struglia, M., Josey, S. A., and Chambers, D.: The effect
of the NAO on sea level and on mass changes in the Mediterranean Sea, J. Geophys. Res.-Oceans, 118, 944–952, 10.1002/jgrc.20078,
2013.Tsimplis, M. N. and Baker, T. F.: Sea level drop in the Mediterranean Sea: an
indicator of deep water salinity and temperature changes?, Geophys. Res. Lett., 27, 1731–1734, 10.1029/1999GL007004, 2000.Tsimplis, M. N. and Shaw, A. G.: The forcing of mean sea level variability
around Europe, Global Planet. Change, 63, 196–202,
10.1016/j.gloplacha.2007.08.018, 2008.Tsimplis, M. N., Marcos, M., and Somot, S.: 21st century Mediterranean sea
level rise: steric and atmospheric pressure contributions from a regional
model, Global Planet. Change, 63, 105–111,
10.1016/j.gloplacha.2007.09.006, 2008.Turuncoglu, U. U. and Sannino, G.: Validation of newly designed regional earth
system model (RegESM) for Mediterranean Basin, Clim. Dynam., 48,
2919–2947, 10.1007/s00382-016-3241-1, 2017.Umgiesser, G., Bajo, M., Ferrarin, C., Cucco, A., Lionello, P., Zanchettin, D., Papa, A., Tosoni, A., Ferla, M., Coraci, E., Morucci, S., Crosato, F., Bonometto, A., Valentini, A., Orlić, M., Haigh, I. D., Nielsen, J. W., Bertin, X., Fortunato, A. B., Pérez Gómez, B., Alvarez Fanjul, E., Paradis, D., Jourdan, D., Pasquet, A., Mourre, B., Tintoré, J., and Nicholls, R. J.: The prediction of floods in Venice: methods, models and uncertainty (review article), Nat. Hazards Earth Syst. Sci., 21, 2679–2704,
10.5194/nhess-21-2679-2021, 2021.UNESCO: General Concepts, UNESCO Land Subsidence International Initiative,
available at: https://www.landsubsidence-unesco.org/general-concepts/ (last access: 5 August 2021),
2020.Vecchio, A., Anzidei, M., Serpelloni, E., and 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, 11, 1480,
10.3390/w11071480, 2019.Vignudelli, S.: Analysis of ERS-1 altimeter collinear passes in the
Mediterranean Sea during 1992–1993, Int. J. Remote Sens.,
18, 573–601, 10.1080/014311697218953, 1997.Vignudelli, S., Kostianoy, A. G., Cipollini, P., and Benveniste, J.: Coastal
altimetry, Springer, Berlin, Heidelberg, 10.1007/978-3-642-12796-0, 2011.Vignudelli, S., Birol, F., Benveniste, J., Fu, L.-L., Picot, N., Raynal, M.,
and Roinard, H.: Satellite altimetry measurements of sea level in the coastal
zone, Surv. Geophys., 40, 1319–1349,
10.1007/s10712-019-09569-1, 2019a.Vignudelli, S., De Biasio, F., Scozzari, A., Zecchetto, S., and Papa, A.: Sea
level trends and variability in the Adriatic Sea and around Venice, in:
Fiducial Reference Measurements for Altimetry, 65–74, Springer, Cham, Switzerland, 10.1007/1345_2018_51, 2019b.Volkov, D. L., Baringer, M., Smeed, D., Johns, W., and Landerer, F. W.:
Teleconnection between the Atlantic meridional overturning circulation and
sea level in the Mediterranean Sea, J. Climate, 32, 935–955,
10.1175/JCLI-D-18-0474.1, 2019.Whitehouse, P. L.: Glacial isostatic adjustment modelling: historical perspectives, recent advances, and future directions, Earth Surf. Dynam., 6, 401–429, 10.5194/esurf-6-401-2018, 2018.Williams, S. D.: CATS: GPS coordinate time series analysis software, GPS
solutions, 12, 147–153, 10.1007/s10291-007-0086-4, 2008.Woodworth, P. L., Pouvreau, N., and Wöppelmann, G.: The gyre-scale circulation of the North Atlantic and sea level at Brest, Ocean Sci., 6, 185–190, 10.5194/os-6-185-2010, 2010.Woodworth, P. L.: Differences between mean tide level and mean sea level,
J. Geodesy, 91, 69–90, 10.1007/s00190-016-0938-1, 2017.Wöppelmann, G. and Marcos, M.: Coastal sea level rise in southern Europe
and the nonclimate contribution of vertical land motion, J. Geophys. Res.-Oceans, 117, C01007, 10.1029/2011JC007469, 2012.Wöppelmann, G. and Marcos, M.: Vertical land motion as a key to
understanding sea level change and variability, Rev. Geophys., 54,
64–92, 10.1002/2015RG000502, 2016.Wöppelmann, G., Zerbini, S., and Marcos, M.: Tide gauges and Geodesy: a
secular synergy illustrated by three present-day case studies, C. R. Geosci., 338, 980–991, 10.1016/j.crte.2006.07.006, 2006.Wöppelmann, G., Marcos, M., Coulomb, A., Míguez, B. M., Bonnetain, P.,
Boucher, C., Gravelle, M., Simon, B., and Tiphaneau, P.: Rescue of the
historical sea level record of Marseille (France) from 1885 to 1988 and its
extension back to 1849–1851, J. Geodesy, 88, 869–885,
10.1007/s00190-014-0728-6, 2014.Zanchettin, D.: Aerosol and solar irradiance effects on decadal climate
variability and predictability, Current Climate Change Reports, 3, 150–162,
10.1007/s40641-017-0065-y, 2017.Zanchettin, D., Rubino, A., Traverso, P., and Tomasino, M.: Impact of
variations in solar activity on hydrological decadal patterns in northern
Italy, J. Geophys. Res.-Atmos., 113, D12102, 10.1029/2007JD009157, 2008.Zanchettin, D., Rubino, A., Traverso, P., and Tomasino, M.: Teleconnections
force interannual-to-decadal tidal variability in the Lagoon of Venice
(northern Adriatic), J. Geophys. Res.-Atmos., 114, D07106, 10.1029/2008JD011485, 2009.Zanchettin, D., Bothe, O., Rubino, A., and Jungclaus, J. H.: Multi-model
ensemble analysis of Pacific and Atlantic SST variability in unperturbed
climate simulations, Clim. Dynam., 47, 1073–1090,
10.1007/s00382-015-2889-2, 2016.
Zanchettin, D., Bruni, S., Thiéblemont, R., and Rubinetti, S.: Data from article “Sea-level rise in Venice: historic and future trends (review article)”, Zenodo [data set], 10.5281/zenodo.5139890, 2021.Zanna, L., Khatiwala, S., Gregory, J. M., Ison, J., and Heimbach, P.: Global
reconstruction of historical ocean heat storage and transport, P. Natl. Acad. Sci. USA, 116, 1126–1131,
10.1073/pnas.1808838115, 2019.Zerbini, S., Raicich, F., Prati, C. M., Bruni, S., Del Conte, S., Errico, M.,
and Santi, E.: Sea-level change in the Northern Mediterranean Sea from
long-period tide gauge time series, Earth-Sci. Rev., 167, 72–87,
10.1016/j.earscirev.2017.02.009, 2017.Zezza, F.: The sedimentary structure of Upper Pleistocene–Holocene deposits in
Venice and its effects on the stability of the historic centre, Rend.
Lincei, 21, 211–227, 10.1007/s12210-010-0089-x, 2010.