Non-linear climate feedback analysis in an atmospheric general circulation model

 A method is described for evaluating the ‘partial derivatives’ of globally averaged top-of-atmosphere (TOA) radiation changes with respect to basic climate model physical parameters. This method is used to analyse feedbacks in the Australian Bureau of Meteorology Research Centre general circulation model. The parameters considered are surface temperature, water vapour, lapse rate and cloud cover. The climate forcing which produces the changes is a globally uniform sea surface temperature (SST) perturbation. The first and second order differentials of model parameters with respect to the forcing (i.e. SST changes) are estimated from quadratic least square fitting. Except for total cloud cover, variables are found to be strong functions of global SST. Strongly non-linear variations of lapse rate and high cloud amount and height appear to relate to the non-linear response in penetrative convection. Globally averaged TOA radiation differentials with respect to model parameters are also evaluated. With the exception of total cloud contributions, a high correlation is generally found to exist, on the global mean level, between TOA radiation and the respective parameter perturbations. The largest non-linear terms contributing to radiative changes are those due to lapse rate and high cloud. The contributions of linear and non-linear terms to the overall radiative response from a 4 K SST perturbation are assessed. Significant non-linear responses are found to be associated with lapse rate, water vapour and cloud changes. Although the exact magnitude of these responses is likely to be a function of the particular model as well as the imposed SST perturbation pattern, the present experiments flag these as processes which cannot properly be understood from linear theory in the evaluation of climate change sensitivity. Received: 16 January 1997/Accepted: 9 May 1997     

An incentive mechanism for reducing emissions from conversion of intact and non-intact forests

This paper presents a new accounting mechanism in the context of the UNFCCC issue on reducing emissions from deforestation in developing countries, including technical options for determining baselines of forest conversions. This proposal builds on the recent scientific achievements related to the estimation of tropical deforestation rates and to the assessment of ‘intact’ forest areas. The distinction between ‘intact’ and ‘non intact’ forests used here arises from experience with satellite-based deforestation measurements and allows accounting for carbon losses from forest degradation. The proposed accounting system would use forest area conversion rates as input data. An optimal technical solution to set baselines would be to use historical average figures during the time period from 1990 to 2005. The system introduces two different schemes to account for preserved carbon: one for countries with high forest conversion rates where the desired outcome would be a reduction in their rates, and another for countries with low rates. A ‘global’ baseline rate would be used to discriminate between these two country categories (high and low rates). For the hypothetical accounting period 2013–2017 and considering 72% of the total tropical forest domain for which data are available, the scenario of a 10% reduction of the high rates and of the preservation of low rates would result in approximately 1.6 billion tCO₂ of avoided emissions. The resulting benefits of this reduction would be shared between those high-rate countries which reduced deforestation and those low-rate countries which did not increase their deforestation over an agreed threshold (e.g., half of “global” baseline rate).     

Two climatic states and feedbacks on thermohaline circulation in an Earth system model of intermediate complexity

The McGill Paleoclimate Model-2 (MPM-2) is employed to study climate–thermohaline circulation (THC) interactions in a pre -industrial climate, with a special focus on the feedbacks on the THC from other climate system components. The MPM-2, a new version of the MPM, has an extended model domain from 90S to 90N, active winds and no oceanic heat and freshwater flux adjustments. In the MPM-2, there are mainly two stable modes for the Atlantic meridional overturning circulation (MOC) under the ‘present-day’ forcing (present-day solar forcing and the pre-industrial atmospheric CO₂ level of 280 ppm). The ‘on’ mode has an active North Atlantic deep water formation, while the ‘off’ mode has no such deep water formation. By comparing the ‘off’ mode climate state with its ‘on’ mode analogue, we find that there exist many large differences between the two climate states, which originate from large changes in the oceanic meridional heat transports. By suppressing or isolating each process associated with a continental ice sheet over North America, sea ice, the atmospheric hydrological cycle and vegetation, feedbacks from these components on the Atlantic MOC are investigated. Sensitivity studies investigating the role of varying continental ice growth and sea ice meridional transport in the resumption of the Atlantic MOC are also carried out. The results show that a fast ice sheet growth and an enhanced southward sea ice transport significantly favor the resumption of the Atlantic MOC in the MPM-2. In contrast to this, the feedback from the atmospheric hydrological cycle is a weak positive one. The vegetation-albedo feedback could enhance continental ice sheet growth and thus could also favor the resumption of the Atlantic MOC. However, before the shut-down of the Atlantic MOC, feedbacks from these components on the Atlantic MOC are very weak.     

Cause-effect relationships between sea surface temperature, precipitation and sea level along the Bangladesh coast

Summary  The Bangladesh coast, which lies on the confluence of three mighty rivers, the Ganges, the Brahmaputra and the Meghna, with the Himalayas to the north and the Bay of Bengal to the south, is an ideal zone for sea level rise due to enhanced rainfall during the monsoon season from June to September. An attempt has been made here to look into the cause-effect relationships between observed trends in sea surface temperature (SST) over the Bay of Bengal and the trends in monsoon rains and sea level in Bangladesh. The study utilizes the 14-year satellite-derived SSTs over the Bay of Bengal for 1985–1998, the tide gauge stations data along the Bangladesh coast for 1977–1998 and the 31-year monsoon rainfall data for Bangladesh, 1961–1991. Received October 20, 2000     

Local and Global Similarity in Turbulent Transfer of Heat, Water Vapour, And CO 2 in the Dynamic Convective Sublayer Over a Suburban Area

We investigated the ‘local’ and ‘global’ similarity of vertical turbulent transfer of heat, water vapour, and CO₂ within an urban surface layer. The results were derived from field measurements in a residential area of Tokyo, Japan during midday on fair-weather days in July 2001. In this study, correlation coefficients and quadrant analysis were used for the evaluation of ‘global’ similarity and wavelet analysis was employed for investigating ‘local’ similarity. The correlation coefficients indicated that the transfer efficiencies of water vapour and CO₂ were generally smaller than that of heat. Using wavelet analysis, we found that heat is always efficiently transferred by thermal and organized motions. In contrast, water vapour and CO₂, which are passive quantities, were not transferred as efficiently as heat. The quadrant analyses showed that the heat transfer by ejection exceeded that by sweep, and the ratios of ejection to sweep for water vapour and CO₂ transfer were less than that for heat. This indicated that heat is more efficiently transferred by upward motions and supported the findings from wavelet analysis. The differences of turbulent transfer between heat and both CO₂ and water vapour were probably caused both by the active role of temperature and the heterogeneity in the source distribution of scalars     

珠江口近15年海平面变化特点及其与强咸潮发生的关系

通过对珠江口30多年相对海平面和近15年绝对海平面变化的研究,比较1992年12月～2008年12月南海卫星观测和珠江口验潮站观测的海平面变化趋势,认为珠江口的相对海平面（RSL）上升最主要原因是全球气候变暖、海平面上升所致;通过研究29个冬季各月西、北江冬季径流量、海平面、表层盐度的变化趋势,以及强咸潮月份的径流、海平面、盐度的对应关系,得出海平面上升是加大咸潮影响的重要因素。     

Reconstructing sea level from paleo and projected temperatures 200 to 2100 ad

We use a physically plausible four parameter linear response equation to relate 2,000 years of global temperatures and sea level. We estimate likelihood distributions of equation parameters using Monte Carlo inversion, which then allows visualization of past and future sea level scenarios. The model has good predictive power when calibrated on the pre-1990 period and validated against the high rates of sea level rise from the satellite altimetry. Future sea level is projected from intergovernmental panel on climate change (IPCC) temperature scenarios and past sea level from established multi-proxy reconstructions assuming that the established relationship between temperature and sea level holds from 200 to 2100 ad. Over the last 2,000 years minimum sea level (−19 to −26 cm) occurred around 1730 ad, maximum sea level (12–21 cm) around 1150 ad. Sea level 2090–2099 is projected to be 0.9 to 1.3 m for the A1B scenario, with low probability of the rise being within IPCC confidence limits.     

Can local linear stochastic theory explain sea surface temperature and salinity variability?

 Sea surface temperature (SST) and salinity (SSS) time series from four ocean weather stations and data from an integration of the GFDL coupled ocean-atmosphere model are analyzed to test the applicability of local linear stochastic theory to the mixed-layer ocean. According to this theory, mixed-layer variability away from coasts and fronts can be explained as a ‘red noise’ response to the ‘white noise’ forcing by atmospheric disturbances. At one weather station, Papa (northeast Pacific), this stochastic theory can be applied to both salinity and temperature, explaining the relative redness of the SSS spectrum. Similar results hold for a model grid point adjacent to Papa, where the relationships between atmospheric energy and water fluxes and actual changes in SST and SSS are what is expected from local linear stochastic theory. At the other weather stations, this theory cannot adequately explain mixed-layer variability. Two oceanic processes must be taken into account: at Panulirus (near Bermuda), mososcale eddies enhance the observed variability at high frequencies. At Mike and India (North Atlantic), variations in SST and SSS advection, indicated by the coherence and equal persistence of SST and SSS anomalies, contribute to much of the low frequency variability in the model and observations. To achieve a global perspective, TOPEX altimeter data and model results are used to identify regions of the ocean where these mechanisms of variability are important. Where mesoscale eddies are as energetic as at Panulirus, indicated by the TOPEX global distribution of sea level variability, one would expect enhanced variability on short time scales. In regions exhibiting signatures of variability similar to Mike and India, variations in SST and SSS advection should dominate at low frequencies. According to the model, this mode of variability is found in the circumpolar ocean and the northern North Atlantic, where it is associated with the irregular oscillations of the model’s thermohaline circulation. Received: 11 March 1996 / Accepted: 6 September 1996     

The early 1950s regime shift in temperature in Taiwan and East Asia

A summer climate regime shift in temperature in Taiwan and East Asia during the early 1950s was identified in this study. The event was characterized by a cooling land-warming ocean dipole in East Asia and the western North Pacific, marking the decreasing land–sea thermal contrast from the 1940s to the 1950s. The corresponding sea surface temperature (SST) anomaly was characterized by the sign flipping of a north–south dipole in the western North Pacific from north/cool-south/warm to north/warm-south/cool, indicating a weakening north–south SST gradient in the area. The corresponding change in mean sea level pressure was characterized by the rising pressure in continental East Asia and the Philippine Sea, and the falling pressure over the extratropical western North Pacific to the east of Japan. This change was the reflection of a weakening thermal low in the continental East Asia, a weakening monsoon trough in the tropical western North Pacific, a strengthening and southwestward-expanding ridge in the subtropical western North Pacific, and a deepening mid-latitude trough over eastern China and Japan. The phase reversal of the SST anomaly in the western North Pacific exhibited the characteristics of the Pacific Decadal Oscillation (PDO), but lagged behind the phase reversal of the PDO in the extratropical North Pacific by several years. The connection with the PDO is speculated, although the mechanism is not understood. This paper is a contribution to the AMIP-CMIP Diagnostic Sub-project on General Circulation Model Simulation of the East Asian Climate, coordinated by W.-C. Wang.     

Regimes of seasonal air–sea interaction and implications for performance of forced simulations

Sea surface temperature (SST) anomalies can induce anomalous convection through surface evaporation and low-level moisture convergence. This SST forcing of the atmosphere is indicated in a positive local rainfall–SST correlation. Anomalous convection can feedback on SST through cloud-radiation and wind-evaporation effects and wind-induced oceanic mixing and upwelling. These atmospheric feedbacks are reflected in a negative local rainfall–SST tendency correlation. As such, the simultaneous rainfall–SST and rainfall–SST tendency correlations can indicate the nature of local air–sea interactions. Based on the magnitude of simultaneous rainfall–SST and rainfall–SST tendency correlations, the present study identifies three distinct regimes of local air–sea interactions. The relative importance of SST forcing and atmospheric forcing differs in these regimes. In the equatorial central-eastern Pacific and, to a smaller degree, in the western equatorial Indian Ocean, SST forcing dominates throughout the year and the surface heat flux acts mainly as a damping term. In the tropical Indo-western Pacific Ocean regions, SST forcing and atmospheric forcing dominate alternatively in different seasons. Atmospheric forcing dominates in the local warm/rainy season. SST forcing dominates with a positive wind-evaporation feedback during the transition to the cold/dry season. SST forcing also dominates during the transition to the warm/rainy season but with a negative cloud-radiation feedback. The performance of atmospheric general circulation model simulations forced by observed SST is closely linked to the regime of air–sea interaction. The forced simulations have good performance when SST forcing dominates. The performance is low or poor when atmospheric forcing dominates.     

Millennial timescale carbon cycle and climate change in an efficient Earth system model

A new Earth system model, GENIE-1, is presented which comprises a 3-D frictional geostrophic ocean, phosphate-restoring marine biogeochemistry, dynamic and thermodynamic sea-ice, land surface physics and carbon cycling, and a seasonal 2-D energy-moisture balance atmosphere. Three sets of model climate parameters are used to explore the robustness of the results and for traceability to earlier work. The model versions have climate sensitivity of 2.8–3.3°C and predict atmospheric CO₂ close to present observations. Six idealized total fossil fuel CO₂ emissions scenarios are used to explore a range of 1,100–15,000 GtC total emissions and the effect of rate of emissions. Atmospheric CO₂ approaches equilibrium in year 3000 at 420–5,660 ppmv, giving 1.5–12.5°C global warming. The ocean is a robust carbon sink of up to 6.5 GtC year⁻¹. Under ‘business as usual’, the land becomes a carbon source around year 2100 which peaks at up to 2.5 GtC year⁻¹. Soil carbon is lost globally, boreal vegetation generally increases, whilst under extreme forcing, dieback of some tropical and sub-tropical vegetation occurs. Average ocean surface pH drops by up to 1.15 units. A Greenland ice sheet melt threshold of 2.6°C local warming is only briefly exceeded if total emissions are limited to 1,100 GtC, whilst 15,000 GtC emissions cause complete Greenland melt by year 3000, contributing 7 m to sea level rise. Total sea-level rise, including thermal expansion, is 0.4–10 m in year 3000 and ongoing. The Atlantic meridional overturning circulation shuts down in two out of three model versions, but only under extreme emissions including exotic fossil fuel resources.     

Regional modeling of decadal rainfall variability over the Sahel

A regional climate model is used to investigate the mechanism of interdecadal rainfall variability, specifically the drought of the 1970s and 1980s, in the Sahel region of Africa. The model is the National Center for Environmental Prediction’s (NCEPs) Regional Spectral Model (RSM97), with a horizontal resolution of approximately equivalent to a grid spacing of 50 km, nested within the ECHAM4.5 atmospheric general circulation model (AGCM), which in turn was forced by observed sea surface temperature (SST). Simulations for the July–September season of the individual years 1955 and 1986 produced wet conditions in 1955 and dry conditions in 1986 in the Sahel, as observed. Additional July–September simulations were run forced by SSTs averaged for each month over the periods 1950–1959 and the 1978–1987. These simulations yielded wet conditions in the 1950–1959 case and dry conditions in the 1978–1987 case, confirming the role of SST forcing in decadal variability in particular. To test the hypothesis that the SST influences Sahel rainfall via stabilization of the tropospheric sounding, simulations were performed in which the temperature field from the AGCM was artificially modified before it was used to force the regional model. We modified the original 1955 ECHAM4.5 temperature profiles by adding a horizontally uniform, vertically varying temperature increase, taken from the 1986–1955 tropical mean warming in either the AGCM or the NCEP/National Center for Atmospheric Research Reanalysis. When compared to the 1955 simulations without the added tropospheric warming, these simulations show a drying in the Sahel similar to that in the 1986–1955 difference and to the decadal difference between the 1980s and 1950s. This suggests that the tropospheric warming may have been, at least in part, the agent by which the SST increases led to the Sahel drought of the 1970s and 1980s.     

Weakened SST variability in the tropical Atlantic Ocean since 2000

A prominent weakening in equatorial Atlantic sea surface temperature (SST) variability, occurring around the year 2000, is investigated by means of observations, reanalysis products and the linear recharge oscillator (ReOsc) model. Compared to the time period 1982–1999, during 2000–2017 the May–June–July SST variability in the eastern equatorial Atlantic has decreased by more than 30%. Coupled air–sea feedbacks, namely the positive Bjerknes feedback and the negative net heat flux damping are important drivers for the equatorial Atlantic interannual SST variability. We find that the Bjerknes feedback weakened after 2000 while the net heat flux damping increased. The weakening of the Bjerknes feedback does not appear to be fully explainable by changes in the mean state of the tropical Atlantic. The increased net heat flux damping is related to an enhanced response of the latent heat flux to the SST anomalies (SSTa). Strengthened trade winds as well as warmer SSTs are suggested to increase the air–sea specific humidity difference and hence, enhancing the latent heat flux response to SSTa. A combined effect of those two processes is proposed to be responsible for the weakened SST variability in the eastern equatorial Atlantic. The ReOsc model supports the link between reduced SST variability, weaker Bjerknes feedback and stronger net heat flux damping.     

Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers

When density fluctuations of scalars such as CO₂ are measured with open-path gas analyzers, the measured vertical turbulent flux must be adjusted to take into account fluctuations induced by ‘external effects’ such as temperature and water vapour. These adjustments are needed to separate the effects of surface fluxes responsible for ‘natural’ fluctuations in CO₂ concentration from these external effects. Analogous to vertical fluxes, simplified expressions for separating the ‘external effects’ from higher-order scalar density turbulence statistics are derived. The level of complexity in terms of input to these expressions are analogous to that of the Webb–Pearman–Leuning (WPL), and are shown to be consistent with the conservation of dry air. It is demonstrated that both higher-order turbulent moments such as the scalar variances, the mixed velocity-scalar covariances, and the two-scalar covariance require significant adjustments due to ‘external effects’. The impact of these adjustments on the turbulent CO₂ spectra, probability density function, and dimensionless similarity functions derived from flux-variance relationships are also discussed.     

Past and future sea-level rise along the coast of North Carolina,USA

We evaluate relative sea level (RSL) trajectories for North Carolina, USA, in the context of tide-gauge measurements and geological sea-level reconstructions spanning the last ~11,000 years. RSL rise was fastest (~7 mm/yr) during the early Holocene and slowed over time with the end of the deglaciation. During the pre-Industrial Common Era (i.e., 0–1800 CE), RSL rise (~0.7 to 1.1 mm/yr) was driven primarily by glacio-isostatic adjustment, though dampened by tectonic uplift along the Cape Fear Arch. Ocean/atmosphere dynamics caused centennial variability of up to ~0.6 mm/yr around the long-term rate. It is extremely likely (probability P=0.95) that 20th century RSL rise at Sand Point, NC, (2.8 ± 0.5 mm/yr) was faster than during any other century in at least 2,900 years. Projections based on a fusion of process models, statistical models, expert elicitation, and expert assessment indicate that RSL at Wilmington, NC, is very likely (P=0.90) to rise by 42–132 cm between 2000 and 2100 under the high-emissions RCP 8.5 pathway. Under all emission pathways, 21st century RSL rise is very likely (P>0.90) to be faster than during the 20th century. Due to RSL rise, under RCP 8.5, the current ‘1-in-100 year’ flood is expected at Wilmington in ~30 of the 50 years between 2050-2100.     

Calculating Historical Contributions To Climate Change – Discussing The ‘Brazilian Proposal’

This paper discusses methodological issues relevant to the calculation of historical responsibility of countries for climate change (‘The Brazilian Proposal’). Using a simple representation of the climate system, the paper compares contributions to climate change using different indicators: current radiative forcing, current GWP-weighted emissions, radiative forcing from increased concentrations, cumulative GWP-weighted emissions, global-average surface-air temperature increase and two new indicators: weighted concentrations (analogue to GWP-weighted emissions) and integrated temperature increase. Only the last two indicators are at the same time ‘backward looking’ (take into account historical emissions), ‘backward discounting’ (early emissions weigh less, depending on the decay in the atmosphere) and ‘forward looking’ (future effects of the emissions are considered) and are comparable for all gases. Cumulative GWP-weighted emissions are simple to calculate but are not ‘backward discounting’. ‘Radiative forcing’ and ‘temperature increase’ are not ‘forward looking’. ‘Temperature increase’ discounts the emissions of the last decade due to the slow response of the climate system. It therefore gives low weight to regions that have recently significantly increased emissions. Results of the five different indicators are quite similar for large groups (but possibly not for individual countries): industrialized countries contributed around 60% to today’s climate change, developing countries around 40% (using the available data for fossil, industrial and forestry CO₂, CH₄ and N₂O). The paper further argues including non-linearities of the climate system or using a simplified linear system is a political choice. The paper also notes that results of contributions to climate change need to be interpreted with care: Countries that developed early benefited economically, but have high historical emission, and countries developing at a later period can profit from developments in other countries and are therefore likely to have a lower contribution to climate change.     

Changes in Vascular Plant Biodiversity in the Netherlands in the 20th Century Explained by their Climatic and other Environmental Characteristics

In the Netherlands nation-wide databases are available with about 10 million records of occurrences of vascular plant species in the 20th century on a scale of approximately 1 km². These data were analysed with a view to identifying relationships between changes in botanical biodiversity and climatic and other environmental factors. Prior to analysis the data were corrected for several major forms of survey bias. The records were broken down into three periods: 1902–1949, 1975–1984 and 1985–1999. Using multiple regression analysis, differences between successive periods were related to plant functional characteristics as explanatory variables. Between the periods 1902–1949 and 1975–1984 there were small but significant increases in the presence of both thermophilic (‘warm’) and psychrophilic (‘cold’) species. However, in the final decades of the 20th century there was a marked increase in thermophilic species only, coinciding with the marked increase in ambient temperature observed during this period, evidence at least of a rapid response of Dutch flora to climate change. Urbanisation was also examined as an alternative explanation for the increase in thermophilic plant species and was found to explain only 50% of the increased presence of such species in the final decades of the 20th century. Besides temperature-related effects, the most important change during the 20th century was a strong decline in oligotrophic and a marked increase in eutrophic plant species.     

Interannual to decadal variability of the winter Aleutian Low intensity during 1900–2004

A new winter Aleutian Low (AL) intensity index was defined in this paper. A centurial-long time series of this index was constructed using the sea level pressure (SLP) data of nearly 100 years. The features of interannual and decadal variability of the winter AL intensity since 1900 were analyzed by applying the wavelet analysis. The relationship between the winter AL intensity and atmospheric circulation was examined. The cross-wavelet analysis technique was used to further reveal the relationship between the AL intensity and sea surface temperature (SST) in the equatorial eastern Pacific (EEP) and tropical Indian Ocean (TIO) in winter. The results indicate that: 1) On the interannual timescale, the winter AL intensity displays 3–7-yr oscillations, while on the decadal timescale, 8–10-yr and 16–22-yr oscillations are more obvious. 2) Of the linkage to atmospheric circulation, both AO (Arctic Oscillation) and PNA (Pacific North America pattern) are closely associated with winter AL intensity on the interannual timescale, but only PNA contributes to the variation of winter AL intensity on the decadal timescale. 3) As to the ocean impact, winter EEP SST is a major factor affecting the winter AL intensity on the interannual timescale, especially on the 3–7-yr periods. However, on the decadal timescale, though both the TIO and EEP SSTs are associated with the AL intensity in winter, the TIO SST impact is more significant.     

Spatial Variability of Both Turbulent Fluxes and Temperature Profiles in an Urban Roughness Layer

The spatial variability of both turbulent flow statistics in the roughness sublayer (RSL) and temperature profiles within and above the canopy layer (CL) were investigated experimentally in a densely built-up residential area in Tokyo, Japan. Using five towers with measuring devices, each tower isolated from the others by at least 200 m, we collected high-frequency measurements of velocity and temperature at a height z=1.8 z _H, where z _H, the mean building height in the area, is 7.3 m. Also, temperature profiles were measured from z=0.4 to 1.8 z _H. The ‘areal mean’ geometric parameters that were obtained for the areas within 200 m of each tower were fairly homogeneous among the tower sites. The main results are as follows: (1) The spatial variability of all RSL turbulent statistics, except the sensible heat flux, was comparable to that reported in a pine forest. Also, the variability decreased with increasing friction velocity. (2) The spatial variability of the RSL sensible heat flux was larger than that reported in a pine forest. Also, the variability depended on the time of the day and became larger in the morning. The difference among the sites was well related to the areal fraction of vegetation. (3) The spatial variability of the CL temperature profile depended on the time of the day and became larger in the morning. Nevertheless, the spatial standard deviation of CL temperature was always below 0.7 K. (4) It is suggested that the “warming-up” process in the morning when heat storage is dominant increases the spatial variation of RSL sensible heat flux and CL temperature according to the local properties around each tower and the variation decreases once there is further convective mixing in the midday     

Sixty-cm Submersion of Venice Discovered Thanks to Canaletto's Paintings

Relative sea level (RSL) rise is a crucial issue for the safeguard of Venice and its historical buildings. The phenomenon over the last three centuries has been investigated by using a proxy of mean sea level: the height of the algae front on palaces. This indicator was accurately drawn by Canaletto and his pupils in their `photographic' paintings made with an optical camera obscura. The positions of the fronts in the 18th century and the present were compared. The RSL rise is due to a combination of natural and anthropogenic factors, both local and global, which affected the land subsidence. An analysis was performed to establish the long-term trend and distinguish between natural and local man contributions. A prudent scenario for the future would suggest a rate between 1.9 ± 0.4 mm yr^–1 and 2.3 ± 0.4 mm yr^–1.