AG_57.04.14_SPADA_finalonline_Layout 6 ANNALS OF GEOPHYSICS, 57, 4, 2014, S0432; doi:10.4401/ag-6548 S0432 Anomalous secular sea-level acceleration in the Baltic Sea caused by isostatic adjustment Giorgio Spada1,*, Marco Olivieri2, Gaia Galassi1 1 Università di Urbino, Dipartimento di Scienze di Base e Fondamenti (DiSBeF), Urbino, Italy 2 Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Italy ABSTRACT Observations from the global array of tide gauges show that global sea- level has been rising at an average rate of 1.5-2 mm/yr during the last ~150 years [Douglas 1991, Spada and Galassi 2012]. Although a global sea-level acceleration was initially ruled out [Douglas 1992], subse- quent studies [Douglas 1997, Church and White 2006, Jevrejeva et al. 2008, Church and White 2011] have coherently proposed values of ~1 mm/year/century [Olivieri and Spada 2013]. More complex non-linear trends and abrupt sea-level variations have now also been recognized. Globally, these could manifest a regime shift between the late Holocene and the current rhythms of sea-level rise [Gehrels and Woodworth 2013], while locally they result from ocean circulation anomalies, steric effects and wind stress [Bromirski et al. 2011, Merrifield 2011]. Although iso- static readjustment affects the local rates of secular sea-level change [Milne and Mitrovica 1998, Peltier 2004], a possible impact on regional acceleration has been so far discounted [Douglas 1992, Jevrejeva et al. 2008, Woodworth et al. 2009] since the process evolves on a millennium time scale [Turcotte and Schubert 2002]. Here we report a previously un- noticed anomaly in the long-term sea-level acceleration of the Baltic Sea tide gauge records, and we explain it by the classical post-glacial rebound theory and numerical modeling of glacial isostasy. Contrary to previous assumptions, our findings demonstrate that isostatic compensation plays a role in the regional secular sea-level acceleration. While studying the regional patterns of the global sea-level acceleration (GSLA) [Olivieri and Spada 2013], we noted that accelerations in secular tide gauge time series from the Baltic Sea tend to exceed, albeit slightly, those from records of comparable length and quality from elsewhere. Tide gauges from this region are gen- erally not exploited in GSLA studies [Church and White 2006, Woodworth et al. 2009] because these are believed to be too strongly influenced by local circula- tion [Jevrejeva et al. 2008] and possibly by GIA. To con- firm the existence of an anomalous acceleration, we have now compared the sea-level trends obtained from 60 validated Baltic Sea records (Figure 1a, BS) with those shown by the global set of tide gauge time-series previously employed by Douglas [1991, 1997] and obey- ing rigorous quality constraints in terms of record length and completeness. Contrary to the former, the latter are not from regions that were deeply covered by ice at the Last Glacial Maximum (~21,000 years ago) where isostatic disequilibrium is still causing crustal up- lift, and consequently a relative sea-level fall of several millimeters per year [Steffen and Wu 2011]. We have augmented the global set with ten records of compa- rable length, recently shown [Spada and Galassi 2012] to minimize the contamination from Glacial Isostatic Adjustment (GIA) and the sensitivity to the model adopted to perform the GIA correction. The resulting global set consists of 32 records (Figure 1b, GL); it in- cludes sites largely insensitive to GIA and thus more suitable to constrain the long term sea-level changes driven by current climate variations. Following Douglas [1992], sea-level accelerations in Figure 1 have been computed as twice the quadratic term of a second-degree best-fitting polynomial. The distribution of the accelerations obtained for sets BS and GL are shown by the histograms of Figure 2a (here, we have found convenient to express accelerations in units of mmyc = mm/yr/century). By a Student t-test with Welch correction [Welch 1947] we have verified that the weighed averages of the two distributions, shown Figure 2b, are significantly different at the 95% confidence level. This supports the hypothesis that the sea-level accelerations for the BS tide gauges are, on the average, larger than those belonging to the GL set. Ac- cording to our computations, the anomaly (i.e. the dif- ference in the average sea-level accelerations for the two sets), can be quantified as Na = (0.2 ± 0.2) mmyc. The anomaly is small, but it is not negligible compared to Article history Received April 17, 2014; accepted June 11, 2014. Subject classification: Sea-level change, Glacial isostatic adjustment, Rheology. GSLA estimates so far obtained in the literature, which are close to 1 mmyc [Olivieri and Spada 2013]. The positive sign for Na can be explained, in a straightforward way, invoking the physics of post- glacial rebound [Turcotte and Schubert 2002]. In for- merly glaciated regions like the Baltic Sea, relative sea-level has been gradually falling since the end of deglaciation ~8,000 years ago, as a direct consequence of the post-glacial uplift of the Earth’s crust [Steffen and Wu 2011]. According to the classical post-glacial re- bound theory for a Newtonian viscous Earth [Turcotte and Schubert 2002], after unloading, the free surface re- laxes toward isostatic equilibrium following an expo- nential law. Hence the rate of relative sea-level fall was larger immediately after the end of ice melting and it is now progressively increasing toward less and less neg- ative values. At the BS tide gauges in Figure 1a, this today implies an excess positive sea-level acceleration compared to distant sites, which are less (or not at all) influenced by the post-glacial rebound process; this ex- plains why the observed anomaly Na is positive. By more quantitative arguments, it is possible to predict the order of magnitude of the anomalous ac- celeration. For a simplified half-space model, the post- glacial displacement u(t) of the Earth’s free surface decays with time exponentially [Turcotte and Schubert 2002]. Hence, the rate of vertical uplift v(t) scales with v(t) = −u(t)/xr , where t is time and xr is the Maxwell re- laxation time of the Earth’s mantle [Turcotte and Schu- bert 2002]. Since vertical uplift dominates the relative sea-level variation s(t) in formerly glaciated regions [Steffen and Wu 2011], neglecting gravitational effects SPADA ET AL. 2 Figure 1. Location of tide gauge instruments. (a) Location of Baltic Sea (BS) tide gauges (number of records: 60, minimum record length: 100 years, average completeness: 92%). In (b) the global (GL) set is considered (32, 56 years, 58%). Shades of blue in (b) depict the thickness of the ice sheets at the Last Glacial Maximum according to the GIA model ICE-5G(VM2) [Peltier 2004]. Colors of squares and diamonds de- note the local sea-level acceleration (1 mmyc = 1 millimeter per year per century). All time series and tide gauge station information are ex- tracted from the PSMSL database of annual RLR (revised local reference) records [Woodworth and Player 2003]. Figure 2. Distribution of observed sea-level accelerations. (a) Sea-level accelerations observed for the Baltic Sea (BS, blue) and the global (GS, red) sets of tide gauges. (b) Average sea-level accelerations for the two sets and their standard deviations. For BS and GL, the accelerations have average values of (0.8 ± 0.1) and (0.6 ± 0.1) mmyc, respectively. Their difference defines the anomalous sea-level acceleration Na = (0.2 ± 0.2) mmyc. 3 [Peltier 2004] we can write s(t) = −u(t), which taking the second derivative gives a(t) = v(t)/xr , where a(t) is the local sea-level acceleration. We note that the rate of present day sea-level change and its acceleration are negatively correlated, namely r(t) = −a(t)xr. The relax- ation time xr is obtained analytically in terms of the bulk properties of the Earth and the spatial wavelength of the ice load [Turcotte and Schubert 2002]. Assum- ing a characteristic uplift rate v ~10 mm/yr, represen- tative of the Baltic region [Steffen and Wu 2011], for xr= 4,400 years [Turcotte and Schubert 2002] (suitable for the size of Fennoscandia) and a bulk mantle viscos- ity [Haskell 1936] of 1021 Pa · s, we obtain a ~0.2 mmyc at present time, consistent with the amplitude of the observed anomaly Na in Figure 2b. Since a(t) scales with v(t), the anomalous acceleration is doomed to disappear as soon as the surface of the Earth will reach, in a few millennia, a full isostatic compensation (v(t) = 0). To substantiate the order-of-magnitude estimates above, we have employed a realistic GIA model to com- pute and visualize the “fingerprints” [Milne and Mitro- vica 1998, Mitrovica et al. 2001] describing patterns of present-day sea-level rates and acceleration associated with GIA. Here we have implemented model ICE- 5G(VM2) [Peltier 2004] for the melting history of the late-Pleistocene ice sheets in the open source code SELEN [Spada and Stocchi 2007, Spada et al. 2012]. Our computations account for the horizontal migration of shorelines and for the impact of Earth rotation varia- tions on sea-level change, and assume a three-layer vis- coelastic mantle with Maxwell rheology with viscosi- ties obtained by volume-averaging profile VM2 [Peltier 2004]. The two fingerprints are shown in Figures 3a and 3b, respectively; the second is obtained from the first by numerical time differentiation. We have verified that their global shape would not change perturbing the vis- cosity profile of the mantle within reasonable bounds, nor adopting alternative GIA models, characterized by different assumptions about the Earth’s viscosity pro- file and including a distinct time chronology for the melting of the late-Pleistocene ice sheets [Lambeck and Johnston 1989, Lambeck et al. 1989]. The GIA fingerprint for the rate of sea-level change in Figure 3a clearly reflects the direct effect of unload- ing, which dominates across the previously glaciated re- gions causing an enhanced sea-level fall. Finer details of the map have been quantitatively explained in terms of “continental levering” (tilting of coastal regions due to ocean water loading) and “ocean syphoning” (migration of water from the far field of former ice sheets to fill the space left by the collapsing isostatic forebulges) [Mitro- vica and Milne 2002]. This pattern is indeed well known since it is commonly employed to decontaminate the rates of sea-level variations observed at tide gauges from the effects of GIA [Douglas 1991, Spada and Galassi 2012]. However, the fingerprint of sea-level acceleration in Figure 3b is shown here for the first time, since pre- vious studies have assumed that GIA has negligible ef- fects on the GSLA [Douglas 1992, 1997, Church and White 2006, Jevrejeva et al. 2008, Woodworth et al. SECULAR SEA-LEVEL ACCELERATION AND GIA Figure 3. GIA sea-level fingerprints. GIA fingerprints for the rate of sea-level rise (a) and its time-derivative (sea-level acceleration) (b), computed at present time for model ICE-5G(VM2) [Peltier 2004] and obtained by code SELEN [Spada and Stocchi 2007, Spada et al. 2012]. In (a) the dots show the PSMSL tide gauges, and the actual rates in polar regions far exceed the range of the color table. Since relative sea-level is defined by the offset between the geoid height and the Earth’s solid surface [Milne and Mitrovica 1998], these fingerprints are also defined across the continents. 2009, Olivieri and Spada 2013] and its importance has been overlooked. Similarly to the correction for the rate of sea-level rise (Figure 3a), that for the acceleration (3b) has a van- ishing average across the oceans, and the mean value at the locations of the global set of PSMSL tide gauges (see Figure 3a) is effectively negligible (<10−2 mmyc), which justifies the assumptions in previous GSLA stud- ies. However, the regional impact of GIA-induced sea- level acceleration is apparent, especially in areas which are still subject to significant isostatic disequilibrium like the Baltic Sea, the Hudsons Bay region, northern Greenland and West Antarctica. Across the Baltic Sea, the GIA acceleration amounts to ~0.1-0.2 mmyc, in agreement with the order of magnitude observed in this study and with the classical post-glacial rebound theory. The anomaly attains its largest amplitude (>0.3 mmyc) in West Antarctica and in central Canada. Since century-long tide gauge observations from the Hud- sons Bay region are not available, it is impossible to compare predictions for this region with instrumental data. The negative correlation between the two finger- prints in Figures 3a and 3b is apparent and consistent with the extremely simplified post-glacial rebound the- ory discussed above. However, their ratio, which can be interpreted as the GIA “relaxation time”, is not constant since the complex rheological layering of the GIA model employed here implies a multi-exponential re- laxation [Peltier 2004, Spada and Stocchi 2007]. Using directly the results from Figure 3, for the Baltic Sea re- gion we obtain xr ~5,000 years, in agreement with the classical Haskell theory [Haskell 1936]. In response to GIA, tide gauge records located along the coasts of the Baltic Sea exhibit a small - but significant - long-term sea-level acceleration in excess to those in the far field of previously glaciated regions. The sign and the amplitude of the anomaly is consis- tent with the post-glacial rebound theory and with re- alistic numerical predictions of GIA models routinely employed to decontaminate the tide gauges observa- tions from the GIA effects [Peltier 2004]. Model com- putations predict the existence of anomalies of similar amplitude in other regions of the globe where GIA is still particularly vigorous at present, but no long-term instrumental observations are available to support their existence. We confirm that a GIA correction for secular sea-level acceleration is not required in GSLA assess- ments because its average value is vanishingly small at the locations of the PSMSL tide gauges [Douglas 1992]. Nevertheless, here we have shown that GIA is con- tributing significantly on a regional scale, and therefore it should be recognized as one of the processes respon- sible for local, long-term sea-level acceleration. Acknowledgements. Party supported by a research grant of Dipartimento di Scienze di Base e Fondamenti (DiSBeF) of the Uni- versity of Urbino. All figures have been drawn using the GMT pack- age [Wessel and Smith 1998]. Florence Colleoni is acknowledged for advice. The tide gauge data are obtained from the Permanent Serv- ice for Mean Sea Level (PSMSL) database (http://www.psmsl.org/). References Bromirski, P.D., A.J. Miller, R.E. Flick and G. Auad (2011). Dynamical suppression of sea level rise along the Pacific coast of North America: Indications for imminent acceleration, J. Geophys. Res., 116, C07005. Church, J.A., and N.J. White (2006). A 20th century ac- celeration in global sea-level, Geophys. Res. Lett., 33, L01602. Church, J.A., and N.J. White (2011). Sea-level rise from the late 19th to the early 21st century, Surv. Geo- phys., 32 (4), 585-602. Douglas, B. (1991). Global sea sevel rise, J. Geophys. Res., 96, 6981-6992. Douglas, B. (1992). Global sea level acceleration, J. Geo- phys. Res., 97 (C8), 12699-12706. Douglas, B. (1997). Global sea-level rise: a redetermi- nation, Surv. Geophys., 18, 279-292. Gehrels, W.R., and P.L. Woodworth (2013). When did modern rates of sea-level rise start?, Global Planet. Change, 100, 263-277. Haskell, N.A. (1936). The motion of a viscous fluid under a surface load, Part II, Physics, 7, 56-61. Jevrejeva, S., J.C. Moore, A. Grinsted and P.L. Wood- worth (2008). Recent global sea level acceleration started over 200 years ago?, Geophys. Res. Lett., 35, L08715. Lambeck, K., and P. Johnston (1989). The viscosity of the mantle: Evidence from analyses of glacial re- bound phenonena, In: I. Jackson (ed.), The Earth’s Mantle: Composition, Structure, and Evolution, Cambridge Univ. Press, New York, 461-502. Lambeck, K., C. Smither and P. Johnston (1989). Sea- level change, glacial rebound and mantle viscosity for northern Europe, Geophys. J. Int., 134, 102-144. Merrifield, M.A. (2011). A shift in western tropical Pa- cific sea level trends during the 1990s, J. Clim., 24 (15), 4126-4138. Milne, G., and J. Mitrovica (1998). Postglacial sea-level change on a rotating Earth, Geophys. J. Int., 133, 1-10. Mitrovica, J., M.E. Tamisiea, J.L. Davis and G. Milne (2001). Recent mass balance of polar ice sheets in- ferred from patterns of global sea-level change, Na- ture, 409, 1026-1029. Mitrovica, J.X., and G. Milne (2002). On the origin of late Holocene sea-level highstands within equatorial ocean basins, Quat. Sci. Rev., 21, 2179-2190. Olivieri, M., and G. Spada (2013). Intermittent sea-level SPADA ET AL. 4 5 acceleration, Global Planet. Change. 109, 64-72. Peltier, W. (2004). Global glacial isostasy and the sur- face of the Ice-Age Earth: the ICE-5G(VM2) model and GRACE, Annu. Rev. Earth Pl. Sc., 32, 111-149. Spada, G., and P. Stocchi (2007). SELEN: a Fortran 90 program for solving the “Sea Level Equation”, Comput. and Geosci., 33, 538-562. Spada, G., and G. Galassi (2012). New estimates of sec- ular sea-level rise from tide gauge data and GIA modeling, Geophys. J. Int., 191 (3), 1067-1094. Spada, G., D. Melini, G. Galassi and F. Colleoni (2012). Modeling sea level changes and geodetic variations by glacial isostasy: the improved SELEN code; http://arxiv.org/abs/1212.5061. Steffen, H., and P. Wu (2011). Glacial isostatic adjust- ment in Fennoscandia - A review of data and mod- eling, Journal of Geodynamics, 52, 169-204. Turcotte, D.L., and G. Schubert (2002). Geodynamics, Cambridge University Press, Cambridge. Welch, B.L. (1947). The generalization of “Student’s” problem when several different population vari- ances are involved, Biometrika, 34, 28-35. Wessel, P., and W.H.F. Smith (1998). New, improved version of generic mapping tools released, EOS, 79, 579. Woodworth, P.L., and R. Player (2003). The Permanent Service for Mean Sea Level: an update to the 21st century, J. Coastal Res., 19, 287-295. Woodworth, P., et al. (2009). Evidence for the acceler- ations of sea level on multi-decade and century timescales, Int. J. Climatol., 29, 777-789. *Corresponding author: Giorgio Spada, Università di Urbino, Dipartimento di Scienze di Base e Fondamenti (DiSBeF), Urbino, Italy; email: giorgio.spada@gmail.com. © 2014 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. SECULAR SEA-LEVEL ACCELERATION AND GIA << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice