Articles | Volume 21, issue 8
Review article
 | Highlight paper
01 Sep 2021
Review article | Highlight paper |  | 01 Sep 2021

Sea-level rise in Venice: historic and future trends (review article)

Davide Zanchettin, Sara Bruni, Fabio Raicich, Piero Lionello, Fanny Adloff, Alexey Androsov, Fabrizio Antonioli, Vincenzo Artale, Eugenio Carminati, Christian Ferrarin, Vera Fofonova, Robert J. Nicholls, Sara Rubinetti, Angelo Rubino, Gianmaria Sannino, Giorgio Spada, Rémi Thiéblemont, Michael Tsimplis, Georg Umgiesser, Stefano Vignudelli, Guy Wöppelmann, and Susanna Zerbini

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.

1 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. 2) and quantification of the individual contributions by the major underlying processes, including vertical land motions (Sect. 3) and climatic changes (Sect. 4). Estimates are supported by a review of downscaling mechanisms of global and large-scale oceanic and atmospheric signals to the Venice Lagoon (Sect. 5), 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. 6). 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. 7).

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 Gregory et al.2019, 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 (Gregory et al.2019), 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 (Gregory et al.2019);

  • Subsidence: land surface sinking (UNESCO2020; see also Gregory et al.2019);

  • 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.

Table 1List of abbreviations used in the paper.

Download Print Version | Download XLSX

The reader is referred to Lionello et al. (2020b, a) and Umgiesser et al. (2020) 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.

2 Monitoring sea-level changes

The monitoring of sea-level changes in Venice relies on both in situ data acquired by tide gauges (Sect. 2.1) and remote sensing observations provided by satellite radar altimetry (Sect. 2.2). 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 (Zerbini et al.2017). Satellite radar altimetry, in contrast, measures GSL changes (Fu and Cazenave2001; Stammer and Cazenave2017). These are independent of variations in the local land level changes (Gregory et al.2019) and hence unaffected by the potentially key component of RSL (Wöppelmann and Marcos2016; Gregory et al.2019). 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.

2.1 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. 1). 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,, 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”;, last access: 5 August 2021).

Figure 1Maps of the study area and major locations and geographical features mentioned in the paper. Top: the Mediterranean Sea (main panel) and satellite altimetry tracks over the northern Adriatic Sea (blue: Envisat 416; red: Jason1-151; green: Jason2-196) (inset); bottom left: the Venice Lagoon; bottom right: the historical city of Venice. Tide-gauge stations are indicated with yellow dots (1: Trieste; 2: Marina di Ravenna; 3: Genoa: 4: Marseille). The map for the bottom right panel is extracted from Google Earth; © 2018 Landsat/Copernicus, © 2018 CNES/Airbus, © 2018 Digital Globe, © 2018 TerraMetric.

Further details on the tide gauges installed in the Venice Lagoon up to the early 20th century are reported by Magrini and Gnesotto (1908). Dorigo (1961) 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. Battistin and Canestrelli (2006) 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” (Rusconi1983; Camuffo and Sturaro2004). According to Dorigo (1961), 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 Dorigo (1961), 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. 2. 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. 2). 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, Zerbini et al. (2017) 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 Woodworth (2017) (, 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. 4.1.

Figure 2Time 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 Battistin and Canestrelli2006).

2.2 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 (Cazenave et al.2019). 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 (Cipollini et al.2008; Fenoglio-Marc et al.2012; Vignudelli1997; Vignudelli et al.2011, 2019a). This has motivated further investigations (Cipollini et al.2013; Passaro et al.2014) based on the latest coastal altimetry datasets (e.g., CTOH; see Birol et al.2017) and/or reprocessing initiatives (e.g., ALES; see Passaro et al.2014). 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 3 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 (Legeais et al.2018). 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 (Rocco2015; Vignudelli et al.2019b).

Figure 3Percentage 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 Passaro et al.2014).

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 (Birol et al.2017) and the advantage of the ALES re-tracker (Passaro et al.2014). The product was tested along track 196 in the Gulf of Trieste. The improvement is particularly good in the entire Gulf of Trieste (ESA CCI2019), confirming what was found by Passaro et al. (2014).

Altimetry-based assessments of multidecadal trends of Venetian GSL are summarized in Sect. 4.2.

3 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 (Carminati and Di Donato1999; Pirazzoli1996).

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 (Carminati and Di Donato1999).

3.1 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 106104 years, and a short-term component controlled by glaciation cycles and due to GIA processes acting on periods of 103104 years (Antonioli et al.2017; Cuffaro et al.2010; Stocchi et al.2005).

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 (Massari et al.2004). 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 (Antonioli et al.2009, and references therein; Carminati and Di Donato1999; Favero et al.1973), nannofossil biostratigraphy, paleomagnetic polarity and magnetic susceptibility (Kent et al.2002). Additional techniques are available for more recent epochs. Radiocarbon dating allows the investigation of organic sediments, mainly peats, up to 50 000 years ago (Bortolami et al.1984), while the depth of archeological remains and historical data provide information for the last few thousand years (Flemming1992). 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 (Gatto and Carbognin1981). 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.

3.1.1 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. 1 and 4), which is the foreland basin of two fold-and-thrust belts: the north-northeastern vergent northern Apennines and the southern vergent southern Alps (Carminati et al.2003). Figure 4 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 (Carminati and Di Donato1999). 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 (Carminati et al.2003; Cuffaro et al.2010).

Figure 4Geological setting of the Po Plain area with dominant tectonic features (adapted from Cuffaro et al.2010).

Kent et al. (2002) 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 (CNR1971) 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, Kent et al. (2002) 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 (103104 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 (Massari et al.2004). 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 (Amorosi et al.2004; Bondesan et al.2006). This pattern may reflect the northeasterly retreat of the Adriatic slab (Cuffaro et al.2010; Ferranti et al.2006).

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 Antonioli et al. (2009) and Lambeck et al. (2011). 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 Lionello et al.2020a) 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 (Ferranti et al.2006).

3.1.2 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 (Tambroni and Seminara2006). 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 (Antonioli et al.2009; Fontana et al.2017; Lambeck et al.2011). 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 Antonioli et al.2009). 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 (Gatto and Carbognin1981; Tosi et al.2009).

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.,  Bortolami et al.1984; Gatto and Carbognin1981). 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 (Bortolami et al.1984). 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. Carminati et al. (2007) 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 2. 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.,  Carminati and Di Donato1999). 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 (Brambati et al.2003). For what concerns uncertainties in geomorphological and historical markers, Antonioli et al. (2009, 2017) proposed a strategy based on archeological metadata and on standard bathymetric corrections for the Holocene and Late Pleistocene (Ferranti et al.2006; Lambeck et al.2004). The resulting median uncertainty for the Venice area is 0.2 mm/year (markers 24–30 in Table 1 of Antonioli et al.2009).

Kent et al. (2002)Carminati et al. (2003)Antonioli et al. (2009)Bortolami et al. (1984)Bortolami et al. (1984)Antonioli et al. (2009)Antonioli et al. (2009)Lezziero (2002)Serandrei Barbero et al. (2001)

Table 2Time evolution of the natural component of land subsidence in the Venetian region over geological timescales.

Download Print Version | Download XLSX

3.1.3 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 Spada2017, and Whitehouse2018). 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 (Melini and Spada2019; Spada and Melini2019). 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 Farrell and Clark (1976), 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 (Peltier2004). Among the processes contributing to present-day RSL change (e.g.,  Milne et al.2009), 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., King et al.2010), 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 (Antonioli et al.2009, 2017, 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 (Lambeck et al.2011) 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 Roy and Peltier (2018) 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 (Lambeck et al.2004), 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 (Tosi et al.2013). Further uncertainties arise from the effects of the melting of the Würm Alpine ice sheet, whose chronology remains uncertain regarding several aspects (Spada et al.2009). 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 Carminati and Di Donato (1999) and Stocchi and Spada (2009) 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.

3.2 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.,  Gambolati et al.2006; Galloway and Riley1999; Erkens et al.2015; Galloway et al.2016). 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 (Dorigo1961; ISPRA2012; Wöppelmann et al.2006) and differences in RSL trends between Venice and Trieste. The post-1930 period is now considered in more detail.

3.2.1 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 (Gambolati et al.1974; Gambolati and Gatto1975; Carbognin et al.1976). Between World War II and 1970, anthropogenic subsidence was a problem common to the whole northern Adriatic coastline (Tosi et al.2010). 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 (Carminati and Di Donato1999; Teatini et al.2005), 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. 5a); the most intensively used aquifer was between 200 and 250 m depth due to its productive character (Carbognin et al.1976).

Figure 5(a) Comparison of the piezometric level in the Marghera industrial area and in Venice from 1910 to 1980 (redrawn after Gatto and Carbognin1981); (b) average piezometric level between 1952 and 1975 along a leveling line from the mainland to Venice (redrawn after Carbognin et al.1976). (c) Empirical curve (red line) accounting for subsidence in Venice (updated from Zerbini et al.2017; 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.

Zerbini et al. (2017)Carbognin et al. (2004)Zerbini et al. (2017)Gatto and Carbognin (1981)Carbognin et al. (1995a, b)Tosi et al. (2013)Tosi et al. (2013)(Blewitt et al.2018)Williams2008Santamaría-Gómez et al.2017(Blewitt et al.2018)Williams2008

Table 3Recent evolution of land subsidence in the historical city center of Venice as measured using geodetic techniques. Tosi et al. (2013) 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.

Download Print Version | Download XLSX

Between 1950 and 1970, human-induced subsidence reached 14 cm at Marghera and averaged 10 cm at Venice (Fig. 5d). 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 (Zezza2010): the ratio between subsidence and piezometric decline was 1/100 at Marghera and 1/50 at Venice (Gatto and Carbognin1981).

3.2.2 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 (Gatto and Carbognin1981). A corresponding rapid piezometric rise occurred (Fig. 5a): in 1978 the hydraulic head rose to ground level, re-establishing the levels existing in 1950 (Gatto and Carbognin1981). Land subsidence slowed concurrently and stopped; by 1975, a surface rebound of about 2 cm was recorded (Fig. 5d), 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 (Tosi et al.2013). 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. Tosi et al. (2018) 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.

3.3 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 (Salvioni1957; Gambolati et al.1974; Gatto and Carbognin1981; Arca and Beretta1985; Carbognin et al.1995a, b). 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 (Carbognin et al.2004; Zerbini et al.2017). 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 (Gambolati et al.1974). The maximum rate was observed in Mestre in 1968–1969, when the local subsidence reached 17 mm/year (Brambati et al.2003). 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. 5). 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 (Carbognin et al.1995a, b).

The monitoring capabilities further improved during the 1990s with the development of space techniques such as GNSS and SAR (Teatini et al.2012; Tosi et al.2013, 2018; Zerbini et al.2017). 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 (Tosi et al.2018). 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 (Tosi et al.2018). 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 (Tosi et al.2013). 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, Tosi et al. (2018) 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 3 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 (Torge and Müller1980). 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 (Santamaría-Gómez et al.2017). This estimate increases to 0.4 mm/year when reference frame uncertainties are considered in the error budget (Santamaría-Gómez et al.2017). Finally, Tosi et al. (2013) 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 3 are both consistent with the average value obtained by the SAR technique over the whole city.

4 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. 6).

4.1 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 Polli (1953), 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 4; 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 5. Marcos and Tsimplis (2008), Wöppelmann and Marcos (2012), and Vecchio et al. (2019) used a linear fit to analyze the RSL data from about 1910 to 2000, obtaining trends between 2.4 and 2.5 mm/year. Vecchio et al. (2019) 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. Zerbini et al. (2017) 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. 5c). 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 5). 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 Zerbini et al. (2017), 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).

Marcos and Tsimplis (2008)Wöppelmann and Marcos (2012)Vecchio et al. (2019)Vignudelli et al. (2019b)

Table 4Linear 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.

Download Print Version | Download XLSX

Carbognin et al. (2010)Zerbini et al. (2017)Scarascia and Lionello (2013)Fenoglio-Marc et al. (2012)Vignudelli et al. (2019b)

Table 5Linear 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.

Download Print Version | Download XLSX

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 (Carbognin et al.2009; Wöppelmann et al.2014; Zerbini et al.2017; Sánchez et al.2018; Fig. 1). The RSL time series of Trieste, Marseille and Genoa exhibit centennial trends between 1.2 and 1.3 mm/year (Marcos and Tsimplis2008; Wöppelmann and Marcos2012; Zerbini et al.2017). 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 Marcos et al. (2016). The 1σ errors are around 0.1 mm/year according to Marcos and Tsimplis (2008) and Wöppelmann and Marcos (2012), and between 0.10 and 0.24 mm/year (90 % confidence) in Vecchio et al. (2019). Zerbini et al. (2017) 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 (Zerbini et al.2017). Scarascia and Lionello (2013) 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 4 and 5 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 Zerbini et al. (2017). The subsidence curve by Zerbini et al. (2017) was updated to 2019 by applying a 1 mm/year trend since 2013 based on the SAR estimate by Tosi et al. (2013) and on the trend exhibited by the PSAL GPS from 2014 onward (Table 3). 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. 5 and 6, 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, Troccoli et al.2012, and Carbognin et al.2010). 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 Church et al.2013). 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; Hay et al.2015), 1.1 ± 0.3 mm/year (99 % confidence interval; Dangendorf et al.2017) and 1.56 ± 0.33 mm/year (90 % confidence interval; Frederikse et al.2020). Figure 6 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. 6a). 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).

Figure 6Comparison 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 Frederikse et al.2020; 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.

4.2 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 (Hay et al.2015; Chen et al.2016; Dangendorf et al.2017; Quartly et al.2017). Regional trends can deviate considerably from the global mean (eg. Scharroo et al.2013; Legeais et al.2018; Cazenave et al.2019). This is also the case of the Mediterranean Sea, which is subject to pronounced spatial and temporal variability (Fig. 7a), with the entire area of the Adriatic Sea exhibiting positive GSL trends that peak in the northern part of the basin.

Figure 7Observed 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 Slangen et al. (2016). (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 Thiéblemont et al. (2019).

Several altimetric datasets have been used to estimate sea-level trends in the Venice area. Fenoglio-Marc et al. (2012) 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). Rocco (2015) 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 (Vignudelli et al.2019b), further reduced to 4.03 ± 1.27 mm/year after removing the seasonal signal (ESA CCI2019). Explanations for the differences between trend estimates in these studies include the different time spans, especially for Fenoglio-Marc et al. (2012); 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 (Fenoglio-Marc et al.2012). This evidence is supportive of a rather uniform sea-level trend in the northern Adriatic (Fig. 7a; see also Bonaduce et al.2016).

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 (Fenoglio-Marc et al.2012; Wöppelmann and Marcos2016). For consistency with post-processed altimetry data (Carrère and Lyard2003), the tide-gauge time series need to be preprocessed for removal of the local effect of atmospheric forcing (see Sect. 5.2). The most recent trend estimates by Vignudelli et al. (2019b) 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. 3 in Zerbini et al. (2017) that represents a best fit of the benchmarks, GPS and PS-InSAR normalized heights.

Table 4 and Fig. 6 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 Vignudelli et al. (2019b) 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 Vignudelli et al. (2019b).

Table 6Performance 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.

Download Print Version | Download XLSX

4.3 Interannual 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 (Allen and Smith1996, and references therein). According to Grinsted et al. (2004), 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 On years, where O means order of magnitude. Based on detrended seasonal Venetian RSL for the period 1872–2003, Zanchettin et al. (2009) report spectral peaks in the autumn (OND) time series at O22 years and at larger multidecadal timescales, at around O2.4 years and at around O3.4 years, with secondary peaks at around O5 years and O8 years. In the winter (JFM) time series, they report significant multidecadal variability at O50 years and larger; at O3.4 years; O8 years; and, less apparent, O5 years. Carbognin et al. (2010) also identify an O8 years 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.4 years and O5 years 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 O8 years (6 % of the total variance) and O50 years (9 % of the total variance) spectral peaks as highly significant. A secondary peak at O16 years 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 O50 years component in the winter series, whose significance then only reaches the 90 % confidence level. Lionello (2005), Barriopedro et al. (2010), Troccoli et al. (2012) and Martínez-Asensio et al. (2016) consistently identify significant decadal variability in the time series of autumn Venetian surge events for the period 1948–2008 (for an updated analysis see Lionello et al.2020b, 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 (O20 years, period 1960–2000) and interannual (O5 years and O8 years, periods 1930–1950 and 1970–1990, respectively) and multidecadal (O50 years, since 1950) fluctuations in winter (Fig. 9). 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. 2 and 6). Marcos and Tsimplis (2008) 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 (Tsimplis and Baker2000; Tsimplis et al.2005; Gomis et al.2008). Scarascia and Lionello (2013) 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 Marcos and Tsimplis (2008) and Scarascia and Lionello (2013) 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 6c 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).

Figure 8Composite 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.

Figure 9Interannual 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 Zerbini et al.2017), (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 Grinsted et al.2004, 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).


5 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.,  Fukumori et al.2007; Calafat et al.2012; Landerer and Volkov2013; Marcos et al.2016; Volkov et al.2019). Such linkage is largely determined by the water mass exchange between the Atlantic Ocean and Mediterranean Sea through the Strait of Gibraltar (e.g.,  Brandt et al.2004; Fukumori et al.2007; Menemenlis et al.2007; Gomis et al.2008; Calafat et al.2012; Landerer and Volkov2013; Tsimplis et al.2013; Adloff et al.2018; Volkov et al.2019). Section 5.1 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., Fukumori et al.2007, and Calafat et al.2012), and overwhelms steric changes as far as basin-average sea-level trends are concerned (Jordà and Gomis2013). 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. 5.2. Another potential contributor to local sea-level changes is freshwater horizontal advection (Jordà and Gomis2013), 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 (Adloff et al.2018). 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. 1), contributes to the shaping of haline properties in the Adriatic Sea (Taricco et al.2015) with a detectable imprint on Venetian sea levels (Zanchettin et al.2009).

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 (Zanchettin et al.2009), showing prominent variability at interannual to decadal timescales in the autumn and winter seasons, and the Atlantic Multidecadal Oscillation (AMO) (Scafetta2014), describing multidecadal fluctuations in North Atlantic sea-surface temperature and influencing atmospheric variability over the Euro-Mediterranean region (e.g.,  Mariotti and Dell'Aquila2012; Maslova et al.2017).

5.1 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 (Armi and Farmer1988) as well as for the substantial vertical mixing that has been observed (García Lafuente et al.2013).

Nonseasonal mass exchanges through the strait up to the decadal timescale are modulated by fluctuations in the winds over and near the strait (Fukumori et al.2007; Landerer and Volkov2013; Volkov et al.2019). 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 (Frederikse et al.2020). Steric changes particularly contribute to multidecadal sea-level variability in the different ocean basins (Frederikse et al.2020), with the Atlantic Ocean accumulating heat at a higher rate than the Pacific and Indian oceans (e.g., Zanna et al.2019). 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 (Frederikse et al.2020).

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. Sannino et al. (2015) 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. Marcos et al. (2016) 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, Adloff et al. (2018) 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.

5.2 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.,  Calafat et al.2012; Jordà et al.2012a, b). The IBE is quantified by the hydrostatic equation in about 1 cm of sea-level rise per 1 hPa of sea-level pressure drop. Calafat et al. (2012) 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 Zanchettin et al. (2009) 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.,  Woodworth et al.2010).

As far as wind forcing is concerned, the morphology of the Adriatic basin displaying a NW–SE elongation and a shallow northern portion (Fig. 1) 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. 8a, 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 (Slivinski et al.2019). 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. 1) 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 Zanchettin et al. (2009) 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.,  Calafat et al.2012; Jordà and Gomis2013; Cusinato et al.2018). 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. 7a (Bonaduce et al.2016).

5.3 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.,  Marshall et al.2001). 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.,  Zanchettin2017; Han et al.2019), 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.,  Tsimplis et al.2006, 2013; Gomis et al.2008; Tsimplis and Shaw2008; Calafat et al.2012; Martínez-Asensio et al.2014; Ezer et al.2016; Rubino et al.2018). 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 (Fukumori et al.2007; Landerer and Volkov2013; Volkov et al.2019). 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., Zanchettin et al.2008; Josey et al.2011; Tsimplis et al.2013; Criado-Aldeanueva et al.2014; Cusinato et al.2018). More specifically for Venetian sea level, Zanchettin et al. (2009) 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 (Jones et al.1997) 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 (Lionello et al.2020b, a). The wavelet coherence spectra (Grinsted et al.2004) in Fig. 9 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. 9), 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.,  Zanchettin et al.2016; Piecuch et al.2019). 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. 8c, 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 (Zanchettin et al.2009) and supports the investigation of other large-scale precursors of atmospheric forcing of Venetian sea level in autumn. For instance, Zanchettin et al. (2009) identify the SCA and EAWR as being related to the interannual to decadal variability of meridional atmospheric flow over the Adriatic Sea in autumn (Zanchettin et al.2009).

6 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 (Nicholls et al.2021). It is useful to distinguish GMSL changes which can be derived from the SROCC (Special Report on Ocean and Cryosphere in a Changing Climate; Oppenheimer et al.2019), regional deviations from GMSL in the Mediterranean and northern Adriatic, and local vertical-land-movement contributions.

6.1 Vertical land movements

Projections of future contribution of vertical land motion are available only for the GIA (see Sect. 6.2 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 (Antonioli et al.2017; Tosi et al.2013), 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 (Tosi et al.2018). 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 (Tosi et al.2013). 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 (Nicholls et al.2014). 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.

6.2 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.,  Kopp et al.2014; DeConto and Pollard2016; Edwards et al.2019) provide the main source of uncertainty in future sea-level change projections (Bakker et al.2017; Oppenheimer et al.2019). 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” (Lempert2019), i.e., a condition where experts lack sufficient knowledge, or parties to a decision cannot agree upon the system processes and futures (see also Haasnoot et al.2020). While high-end GMSL rises are unlikely, they cannot be excluded (e.g., Nicholls et al.2014; Kopp et al.2017) and are particularly important for decision-making applications and adaptation planning, especially in decision contexts with low tolerance to uncertainty (Hinkel et al.2019). According to Slangen et al. (2017), 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 (Tsimplis and Shaw2008; Tsimplis et al.2008; Jordà et al.2012b; Jordà2014; Adloff et al.2018). 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., Tsimplis et al.2008). 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 (Adloff et al.2018). 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 (Carillo et al.2012). By the end of the 21st century (i.e., 2070–2099), for the steric component alone, Tsimplis et al. (2008) found a 13 cm sea-level anomaly (with respect to 1960–1990) under the A2 scenario, while Adloff et al. (2015) 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. 7b). 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 Kopp et al. (2014).

Estimates of future sea-level rise in the northern Adriatic Sea from a statistical model are provided by Scarascia and Lionello (2013). 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 (Jordà and Gomis2013). Scarascia and Lionello (2013) 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 Spada et al. (2013), Galassi and Spada (2014) 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.,  Slangen et al.2012, 2014; Kopp et al.2014; Gregory et al.2019). Kopp et al. (2014) 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., Bamber and Aspinall2013) 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 Kopp et al.2014) 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 Kopp et al. (2014) 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 Parras-Berrocal et al.2020), which in turn affects regional sea-level estimates (Marcos and Tsimplis2008; Slangen et al.2017). 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., Landerer et al.2014).

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 Meyssignac et al. (2017) and Thiéblemont et al. (2019). 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. 2 in Thiéblemont et al.2019). 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 (Mitrovica et al.2009; Gregory et al.2019) 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:

(1) σ tot 2 = ( σ sterodynamic + σ smb-a  + σ smb-g ) 2 + σ Glac 2 + σ LW 2 + σ dyn-a 2 + σ dyn-g 2 ,

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 Church et al.2013, 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. 7c). The high-end scenario – obtained by selecting, for each sea-level component, the highest physically based estimate found in the literature (see Thiéblemont et al.2019 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 Scarascia and Lionello (2013). Our projections largely overlap with those provided by Kopp et al. (2014) 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 (Oppenheimer et al.2019). 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., Ranger et al.2013; Environment Agency2021). 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., Jordà2014). 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.

7 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 (Ponte et al.2019). 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., Marti et al.2019) and their implications for coastal hazards (Ablain et al.2016; Benveniste et al.2019). It has been shown that by improving processing, it is possible to make more accurate sea-level measurements in coastal zones (Cipollini et al.2017). 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 (Passaro et al.2014). 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 (Gómez-Enri et al.2019). An updated table is accessible at (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 (De Biasio et al.2020).

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 (Fernandes and Lázaro2016). 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 (Vignudelli et al.2019a).

Concerning the monitoring of vertical land motions, integrated systems have been shown to offer the best approach to the study of subsidence (Wöppelmann and Marcos2016; Zerbini et al.2017; Tosi et al.2009): 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 (Nicholls et al.2021). 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.,  Somot et al.2008; Sannino et al.2009; Artale et al.2010; Naranjo et al.2014; Sein et al.2015; Turuncoglu and Sannino2017; Androsov et al.2019; Palma et al.2020). 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 (Simoncelli et al.2014). 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.,  Artale et al.2006). 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 (Tsimplis and Baker2000). 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.,  Marshall et al.2001; Volkov et al.2019), but this requires further investigation (Piecuch et al.2019). 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.,  Polyakov et al.2012; Holt et al.2017). In this context, the unstructured approach adopted by Ferrarin et al. (2018) 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 Umgiesser et al.2020, 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 (Pisacane et al.2006; Amitai et al.2017). 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 (Good et al.2013) or Ishii et al. (2017), 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 (Dodge et al.2003). 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 (Alessio2016, 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, Chambers2015). 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.,  Slangen et al.2012; Kopp et al.2014; Oppenheimer et al.2019; Thiéblemont et al.2019) 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 Kopp et al. (2014) 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 (Eyring et al.2016).

8 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 (Zanchettin et al.2021).

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.


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.


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 (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.


Ablain, 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: (last access: 5 August 2021), 2016. a

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,, 2015. a

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,, 2018. a, b, c, d, e

Alessio, S. M.: Digital signal processing and spectral analysis for scientists: concepts and applications, Springer, Cham, Switzerland, 2016. a

Allen, M. and Smith, L.: Monte Carlo SSA: Detecting irregular oscillations in the presence of colored noise, J. Climate, 9, 3373–3404,<3373:MCSDIO>2.0.CO;2, 1996. a

Amitai, Y., Ashkenazy, Y., and Gildor, H.: Multiple equilibria and overturning variability of the Aegean-Adriatic Seas, Global Planet. Change, 151, 49–59,, 2017. a

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,, 2004. a

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,, 2019. a

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,, 2009. a, b, c, d, e, f, g, h, i, j

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,, 2017. a, b, c, d

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. a

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. a

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,, 2006. a

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,, 2010. a

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,, 2017. a

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,, 2013. a

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,, 2010. a

Battistin, D. and Canestrelli, P.: 1872–2004: la serie storica delle maree a Venezia, Proc. 3rd International Symp. on Land Subsidence, available at: (last access: 5 August 2021), 2006. a, b

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,, 2019. a

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,, 2017. a, b

Blewitt, G., Hammond, W. C., and Kreemer, C.: Harnessing the GPS data explosion for interdisciplinary science, Eos, 99, 1–2,, 2018. a, b

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,, 2016. a, b

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. a

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. a, b, c, d, e

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,, 2003. a, b

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,<0433:SLVITW>2.0.CO;2, 2004. a

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,, 2012. a, b, c, d, e, f, g

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,, 2004. a

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. a, b, c

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. a, b, c

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. a, b, c

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,, 2004. a, b

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,, 2009. a

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,, 2010. a, b, c

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,, 2012. a

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,, 1999. a, b, c, d, e, f, g

Carminati, E., Doglioni, C., and Scrocca, D.: Apennines subduction-related subsidence of Venice (Italy), Geophys. Res. Lett., 30, 1717,, 2003. a, b, c

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,, 2007. a

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,, 2003. a

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,, 2019. a, b

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,, 2015. a

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,, 2016. a

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: (last access: 5 August 2021), 2013. a, b

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,, 2008. a

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: (last access: 5 August 2021), 2013. a

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,, 2017. a

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. a

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,, 2014. a

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,, 2010. a, b, c, d

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,, 2018. a, b

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,, 2017. a, b

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,, 2020. a

DeConto, R. M. and Pollard, D.: Contribution of Antarctica to past and future sea-level rise, Nature, 531, 591–597,, 2016. a

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. a

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: (last access: 5 August 2021), 1961. a, b, c, d

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,, 2019. a

Environment Agency: Policy paper: Thames Estuary 2100: 10-Year Review monitoring key findings, Tech. rep., Environment Agency, UK, available at:, last access: 23 March 2021. a

Erkens, G., Bucx, T., Dam, R., De Lange, G., and Lambert, J.: Sinking coastal cities, P. Int. Ass. Hydrol. Sci., 372, 189–198,, 2015. a

ESA CCI: Sea Level Bridging Phase Part 2: Validation. Technical Deliverable v 2.0 2017, available at: (last access: 5 August 2021), 2019. a, b

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,, 2016. a

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,, 2016. a

Farrell, W. and Clark, J. A.: On postglacial sea level, Geophys. J. Int., 46, 647–667,, 1976. a

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. a

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,, 2012. a, b, c, d, e, f

Fernandes, M. J. and Lázaro, C.: GPD+ wet tropospheric corrections for CryoSat-2 and GFO altimetry missions, Remote Sens., 8, 851,, 2016. a

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,, 2006. a, b, c

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,, 2018. a

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. a

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,, 2017. a

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,, 2020. a, b, c, d, e

Fu, L.-L. and Cazenave, A.: Satellite altimetry and earth sciences: a handbook of techniques and applications, Academic Press, London, 2001. a

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,, 2007. a, b, c, d, e

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,, 2014. a

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. a

Galloway, D. L., Erkens, G., Kuniansky, E. L., and Rowland, J. C.: Preface: Land subsidence processes, Hydrogeol. J., 24, 547–550,, 2016. a

Gambolati, G. and Gatto, P.: Simulazione della subsidenza di Venezia, Venezia ei problemi dell'ambiente, Il Mulino, Bologna, 299–360, 1975. a

Gambolati, G., Gatto, P., and Freeze, R. A.: Mathematical simulation of the subsidence of Venice: 2. Results, Water Resour. Res., 10, 563–577,, 1974. a, b, c

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,, 2006. a

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,, 2013. a

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,, 1981. a, b, c, d, e, f, g, h, i

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,, 2019. a

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,, 2008. a, b, c

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,, 2013. a

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,, 2019. a, b, c, d, e, f, g, h

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,, 2004. a, b, c

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,, 2020. a

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,, 2019. a

Hay, C. C., Morrow, E., Kopp, R. E., and Mitrovica, J. X.: Probabilistic reanalysis of twentieth-century sea-level rise, Nature, 517, 481–484,, 2015. a, b

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,, 2019. a

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,, 2017. a

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,, 2017. a

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. a

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,<1433::AID-JOC203>3.0.CO;2-P, 1997. a

Jordà, G.: Detection time for global and regional sea level trends and accelerations, J. Geophys. Res.-Oceans, 119, 7164–7174,, 2014. a, b

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,, 2013. a, b, c, d

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,, 2012a. a

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,, 2012b. a, b

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,, 2011. a

Kent, D. V., Rio, D., Massari, F., Kukla, G., and Lanci, L.: Emergence of Venice during the Pleistocene, Quat. Sci. Rev., 21, 1719–1727,, 2002. a, b, c, d

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,, 2010. a

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,, 2014. a, b, c, d, e, f, g, h, i

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,, 2017. a

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,, 2004. a, b

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,, 2011. a, b, c

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,, 2013. a, b, c, d

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,, 2014. a

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,, 2018. a, b

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. a

Lezziero, A.: Indagini paleoambientali nel sottosuolo dell’isola di Burano, Quaderni di Insula, Venezia, 10, 63–70, 2002. a

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. a

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,, 2021a. a, b, c

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,, 2021b. a, b, c

Magrini, G. P., De Marchi, L., and Gnesotto, T.: Impianti mareografici eseguiti, Ricerche Lagunari, available at: (last access: 5 August 2021), 1908. a

Marcos, M. and Tsimplis, M. N.: Coastal sea level trends in Southern Europe, Geophys. J. Int., 175, 70–82,, 2008. a, b, c, d, e, f, g

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,, 2016. a, b, c

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,, 2012. a

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,, 2001. a, b

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,, 2019. a

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,, 2014. a

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,, 2016. a

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,, 2017. a

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,, 2004. a, b

Melini, D. and Spada, G.: Some remarks on Glacial Isostatic Adjustment modelling uncertainties, Geophys. J. Int., 218, 401–413,, 2019. a

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,, 2007. a

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,, 2017. a

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,, 2009. a

Mitrovica, J. X., Gomez, N., and Clark, P. U.: The sea-level fingerprint of West Antarctic collapse, Science, 323, 753–753,, 2009. a

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,, 2014. a

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,, 2014. a, b

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,, 2021. a, b

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. a, b, c, d

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,, 2020. a

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,, 2020. a

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,, 2014. a, b, c, d, e, f

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,, 2004. a

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,, 2019. a, b

Pirazzoli, P.: Sea level changes. The last 2000 years, John Wiley & Sons, Chichester, 211 pp., 1996. a

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,, 2006. a

Polli, S.: Il graduale aumento del livello del mare lungo le coste italiane, Geofisica pura e applicata, 25, 123–129,, 1953. a

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,, 2012. a

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,, 2019. a

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,, 2017. a

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,, 2013. a

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:, (last access: 5 August 2021), 2015. a, b

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,, 2018. a

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,, 2018. a

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. a

Salvioni, G.: I movimenti del suolo nell'Italia centro-settentrionale, Bollettino di geodesia e scienze affini, 16, 325–366, 1957. a

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,, 2018. a

Sannino, G., Pratt, L., and Carillo, A.: Hydraulic criticality of the exchange flow through the Strait of Gibraltar, J. Phys. Oceanogr., 39, 2779–2799,, 2009. a

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,, 2015. a

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,, 2017. a, b, c

Scafetta, N.: Multi-scale dynamical analysis (MSDA) of sea level records versus PDO, AMO, and NAO indexes, Clim. Dynam., 43, 175–192,, 2014. a

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,, 2013. a, b, c, d, e, f, g

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. a

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,, 2015. a

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. a

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],, 2014. a

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,, 2012. a, b

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,, 2014. a

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,, 2017. a, b

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,, 2016. a

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,, 2019. a

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,, 2008. a

Spada, G.: Glacial isostatic adjustment and contemporary sea level rise: An overview, Surv. Geophys., 38, 153–185,, 2017. a

Spada, G. and Melini, D.: On Some Properties of the Glacial Isostatic Adjustment Fingerprints, Water, 11, 1844,, 2019. a

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,, 2009. a

Spada, G., Bamber, J., and Hurkmans, R.: The gravitationally consistent sea-level fingerprint of future terrestrial ice loss, Geophys. Res. Lett., 40, 482–486,, 2013. a

Stammer, D. and Cazenave, A.: Satellite altimetry over oceans and land surfaces, CRC Press, Boca Raton,, 2017. a

Stocchi, P. and Spada, G.: Influence of glacial isostatic adjustment upon current sea level variations in the Mediterranean, Tectonophysics, 474, 56–68,, 2009. a

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,, 2005. a

Tambroni, N. and Seminara, G.: Are inlets responsible for the morphological degradation of Venice Lagoon?, J. Geophys. Res.-Earth, 111, F03013,, 2006. a

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,, 2015. a

Teatini, P., Ferronato, M., Gambolati, G., Bertoni, W., and Gonella, M.: A century of land subsidence in Ravenna, Italy, Environ. Geol., 47, 831–846,, 2005. a

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,, 2012. a

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,, 2019. a, b, c, d, e

Torge, W. and Müller, J.: Geodesy, Walter de Gruyter, Berlin, Germany, 1980. a

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,, 2009. a, b

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,, 2010. a

Tosi, L., Teatini, P., and Strozzi, T.: Natural versus anthropogenic subsidence of Venice, Sci. Rep., 3, 2710,, 2013. a, b, c, d, e, f, g, h, i, j, k

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,, 2018. a, b, c, d, e, f

Troccoli, A., Zambon, F., Hodges, K. I., and Marani, M.: Storm surge frequency reduction in Venice under climate change, Clim. Change, 113, 1065–1079,, 2012. a, b

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,, 2005. a

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,, 2006. a

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,, 2013. a, b, c

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,, 2000. a, b

Tsimplis, M. N. and Shaw, A. G.: The forcing of mean sea level variability around Europe, Global Planet. Change, 63, 196–202,, 2008. a, b

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,, 2008. a, b, c

Turuncoglu, U. U. and Sannino, G.: Validation of newly designed regional earth system model (RegESM) for Mediterranean Basin, Clim. Dynam., 48, 2919–2947,, 2017. a

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,, 2021. a, b

UNESCO: General Concepts, UNESCO Land Subsidence International Initiative, available at: (last access: 5 August 2021), 2020. a

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,, 2019. a, b, c, d

Vignudelli, S.: Analysis of ERS-1 altimeter collinear passes in the Mediterranean Sea during 1992–1993, Int. J. Remote Sens., 18, 573–601,, 1997. a

Vignudelli, S., Kostianoy, A. G., Cipollini, P., and Benveniste, J.: Coastal altimetry, Springer, Berlin, Heidelberg,, 2011. a

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,, 2019a. a, b

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,, 2019b. a, b, c, d, e, f, g

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,, 2019. a, b, c, d, e

Whitehouse, P. L.: Glacial isostatic adjustment modelling: historical perspectives, recent advances, and future directions, Earth Surf. Dynam., 6, 401–429,, 2018. a

Williams, S. D.: CATS: GPS coordinate time series analysis software, GPS solutions, 12, 147–153,, 2008. a, b

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,, 2010. a

Woodworth, P. L.: Differences between mean tide level and mean sea level, J. Geodesy, 91, 69–90,, 2017. a

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,, 2012. a, b, c, d

Wöppelmann, G. and Marcos, M.: Vertical land motion as a key to understanding sea level change and variability, Rev. Geophys., 54, 64–92,, 2016. a, b, c

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,, 2006. a

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,, 2014. a

Zanchettin, D.: Aerosol and solar irradiance effects on decadal climate variability and predictability, Current Climate Change Reports, 3, 150–162,, 2017. a

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,, 2008. a

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,, 2009. a, b, c, d, e, f, g, h, i

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,, 2016.  a

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],, 2021. a

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,, 2019. a

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,, 2017. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s

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,, 2010. a

Short summary
Relative sea level in Venice rose by about 2.5 mm/year in the past 150 years due to the combined effect of subsidence and mean sea-level rise. We estimate the likely range of mean sea-level rise in Venice by 2100 due to climate changes to be between about 10 and 110 cm, with an improbable yet possible high-end scenario of about 170 cm. Projections of subsidence are not available, but historical evidence demonstrates that they can increase the hazard posed by climatically induced sea-level rise.
Final-revised paper