{"id":969,"date":"2009-08-31T06:10:18","date_gmt":"2009-08-31T11:10:18","guid":{"rendered":"http:\/\/www.realclimate.org\/?p=969"},"modified":"2009-11-08T17:22:41","modified_gmt":"2009-11-08T22:22:41","slug":"ups-and-downs-of-sea-level-projections","status":"publish","type":"post","link":"https:\/\/www.realclimate.org\/index.php\/archives\/2009\/08\/ups-and-downs-of-sea-level-projections\/","title":{"rendered":"Ups and downs of sea level projections"},"content":{"rendered":"<div class=\"kcite-section\" kcite-section-id=\"969\">\n<p><em>By Stefan Rahmstorf and Martin Vermeer<\/em><\/p>\n<p>The scientific sea level discussion has moved a long way since the last IPCC report was published in 2007 (see our <a href=\"http:\/\/www.realclimate.org\/index.php\/archives\/2007\/03\/the-ipcc-sea-level-numbers\/ \">post back then<\/a>). The <a href=\"http:\/\/www.pik-potsdam.de\/news\/press-releases\/files\/synthesis-report-web.pdf\">Copenhagen Synthesis Report <\/a>recently concluded that \u201cThe updated estimates of the future global mean sea level rise are about double the IPCC projections from 2007\u2033. <em>New Scientist<\/em> last month ran a <a href=\"http:\/\/www.newscientist.com\/article\/mg20327151.300-sea-level-rise-its-worse-than-we-thought.html?page=1\">nice article<\/a> on the state of the science, very much in the same vein. But now Mark Siddall, Thomas Stocker and Peter Clark have countered this trend in an article in <em>Nature Geoscience<\/em>, projecting a global rise of only 7 to 82 cm from 2000 to the end of this century.<\/p>\n<p><center><br \/>\n<img decoding=\"async\" data-src=\"http:\/\/www.realclimate.org\/wp-content\/uploads\/atlantis_11.jpg\" alt=\"Coastal erosion: Like the Dominican Republic, many island nations are particularly vulnerable to sea level rise. (c) S.R.\" title=\"atlantis_1\" src=\"data:image\/svg+xml;base64,PHN2ZyB3aWR0aD0iMSIgaGVpZ2h0PSIxIiB4bWxucz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciPjwvc3ZnPg==\" class=\"lazyload\" \/><br \/>\n<small>Coastal erosion: Like the Dominican Republic, many island nations are<br \/>\nparticularly vulnerable to sea level rise. (Photo: S.R.)<\/small><br \/>\n<\/center><\/p>\n<p><!--more--><br \/>\n<strong>Semi-empirical sea level models<\/strong><\/p>\n<p>Siddall et al. use a semi-empirical approach similar to the one Stefan proposed in <em>Science <\/em>in 2007 (let\u2019s call that R07) and to Grinsted et al. (2009), which we <a href=\"http:\/\/www.realclimate.org\/index.php\/archives\/2009\/01\/sea-will-rise-to-levels-of-last-ice-age\/ \">discussed here<\/a>. What are the similarities and where do the differences come from?<\/p>\n<p>For short time scales and small temperature changes everything becomes linear and the two new approaches are mathematically equivalent to R07 (see footnote 1). They can all be described by the simple equation:<\/p>\n<p><center>dS\/dt = a \u0394T(t) + b  &nbsp;&nbsp;&nbsp;      <em>(Eq 1)<\/em><\/center><\/p>\n<p>dS\/dt is the rate of change of sea level S, \u0394T is the warming above some baseline temperature, and <em>a <\/em>and <em>b <\/em>are constants. The baseline temperature can be chosen arbitrarily since any constant temperature offset can be absorbed into <em>b<\/em>. This becomes clear with an example: Assume you want to compute sea level rise from 1900-2000, using as input a temperature time series like the <a href=\"http:\/\/data.giss.nasa.gov\/gistemp\/tabledata\/GLB.Ts+dSST.txt\">global GISS data<\/a>. A clever choice of baseline temperature would then be the temperature around 1900 (averaged over 20 years or so, we\u2019re not interested in weather variability here). Then you can integrate the equation from 1900 to 2000 to get sea level relative to 1900:<\/p>\n<p><center>S(t) = a \u222b\u0394T(t\u2019) dt\u2019 + b t    &nbsp;&nbsp;&nbsp;     <em>(Eq 2)<\/em><\/center><\/p>\n<p>There are two contributions to 20th C sea level rise: one from the warming in the 20th Century (let\u2019s call this the \u201cnew rise\u201d), and a sea level rise that results from any climate changes prior to 1900, at a rate <em>b<\/em> that was already present in 1900 (let\u2019s call this the \u201cold rise\u201d).  This rate is constant for 1900-2000 since the response time scale of sea level is implicitly assumed to be very long in Eq. 1. A simple matlab\/octave code is provided below (2). <\/p>\n<p>If you\u2019re only interested in the total rise for 1900-2000, the temperature integral over the GISS data set is 25 \u00baC years, which is just another way of saying that the mean temperature of the 20th Century was 0.25 \u00baC above the 1900 baseline. The sea level rise over the 20th Century is thus:<\/p>\n<p><center>S(1900-2000) = 25 a + 100 b   &nbsp;&nbsp;&nbsp;      <em>(Eq. 3)<\/em><\/center><\/p>\n<p>Compared to Eq. 1, both new studies introduce an element of non-linearity. In the approach of Grinsted et al, sea level rise may flatten off (as compared to what Eq 1 gives) already on time scales of a century, since they look at a single equilibration time scale \u03c4 for sea level with estimates ranging from 200 years to 1200 years. It is a valid idea that part of sea level rise responds on such time scales, but this is unlikely to be the full story given the long response time of big ice sheets.<\/p>\n<p>Siddall et al. in contrast find a time scale of 2900 years, but introduce a non-linearity in the equilibrium response of sea level to temperature (see their curve in Fig. 1 and footnote 3 below): it flattens off strongly for warm temperatures. The reason for both the long time scale and the shape of their equilibrium curve is that this curve is dominated by ice volume changes. The flattening at the warm end is because sea level has little scope to rise much further once the Earth has run out of ice. However, their model is constructed so that this equilibrium curve determines the rate of sea level rise right from the beginning of melting, when the shortage of ice arising later should not play a role yet. Hence, we consider this nonlinearity, which is partly responsible for the lower future projections compared to R07, physically unrealistic. In contrast, there are some good reasons for the assumption of linearity (see below).<\/p>\n<p><strong>Comparison of model parameters<\/strong><\/p>\n<p>But back to the linear case and Eq. 1: how do the parameter choices compare? <em>a<\/em> is a (more or less) universal constant linking sea level to temperature changes, one could call it the <em>sea level sensitivity<\/em>. <em>b<\/em> is more situation-specific in that it depends both on the chosen temperature baseline and the time history of previous climate changes, so one has to be very careful when comparing <em>b<\/em> between different models.<\/p>\n<p>For R07, and referenced to a baseline temperature for the year 1900, we get<em> a<\/em> = 0.34 cm\/\u00baC\/year and <em>b<\/em> = 0.077 cm\/year. Corresponding values of Grinsted et al. are shown in the table (thanks to Aslak for giving those to us!).<\/p>\n<p>For Siddall et al, <em>a<\/em> = <em>s<\/em>\/\u03c4 where <em>s<\/em> is the slope of their sea level curve, which near present temperatures is 4.8 meters per \u00baC and \u03c4 is the response the time scale. Thus <em>a<\/em> = 0.17 cm\/\u00baC\/year and <em>b<\/em> = 0.04 cm \/year (see table). The latter can be concluded from the fact that their 19th Century sea level rise, with flat temperatures (\u0394T(t) = 0) is 4 cm. Thus, in the model of Siddall et al, sea level (near the present climate) is only half as sensitive to warming as in R07. This is a second reason why their projection is lower than R07.<\/p>\n<table width=\"100%\" border=\"1\" cellspacing=\"0\">\n<tr>\n<td><font size=\"-1\">Model <br \/><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">a [cm\/&ordm;C\/year]<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">b<br \/>\n        [cm \/year]<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">&#8220;new rise&#8221; [cm] (25a)<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">&#8220;old rise&#8221; [cm] (100b)<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">25a+100b<br \/>\n        [cm] <\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">total model rise [cm]<\/div>\n<p><\/font><\/td>\n<\/tr>\n<tr>\n<td><font size=\"-1\">Rahmstorf<\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.34<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.077<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">8.5<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">7.7 <\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">16.2<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">16.2<\/div>\n<p><\/font><\/td>\n<\/tr>\n<tr>\n<td><font size=\"-1\">Grinsted et al &#8220;historical&#8221;<\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.30<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.141<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">7.5<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">14.1<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">21.6<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">21.3<\/div>\n<p><\/font><\/td>\n<\/tr>\n<tr>\n<td><font size=\"-1\">Grinsted et al &#8220;Moberg&#8221;<\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.63<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.085<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">(15.8)<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">(8.5)<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">(24.3) <\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">20.6<\/div>\n<p><\/font><\/td>\n<\/tr>\n<tr>\n<td><font size=\"-1\">Siddall et al<\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.17<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">0.04<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">4.3<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">4<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\">8.3<\/div>\n<p><\/font><\/td>\n<td><font size=\"-1\"><\/p>\n<div align=\"center\"><del>8.3 (?)<\/del> 7.9<\/div>\n<p><\/font><\/td>\n<\/tr>\n<\/table>\n<p><strong><br \/>\nPerformance for 20th Century sea level rise<\/strong><\/p>\n<p>For the 20th Century we can compute the \u201cnew\u201d sea level rise due to 20th Century warming and the \u201cold\u201d rise due to earlier climate changes from Eq. 3. The results are shown in the table. From Grinsted et al, we show two versions fitted to different data sets, one only to \u201chistorical\u201d data using the Jevrejeva et al. (2006) sea level from 1850, and one using the Moberg et al. (2006) temperature reconstruction with the extended Amsterdam sea level record starting in the year 1700.<\/p>\n<p>First note that \u201cold\u201d and \u201cnew\u201d rise are of similar magnitude for the 20th Century because of the small average warming of 0.25 \u00baC. But it is the <em>a<\/em>-term in Eq. (2) that matters for the future, since with future warming the temperature integral becomes many times larger. It is thus important to realise that the total 20th Century rise is not a useful data constraint on <em>a<\/em>, because one can get this right for any value of <em>a<\/em> as long as <em>b<\/em> is chosen accordingly. To constrain the value of <em>a<\/em> \u2013 which dominates the 21st Century projections &#8212; one needs to look at the \u201cnew rise\u201d. How much has sea level rise accelerated over the 20th Century, in response to rising temperatures? That determines how much it will accelerate in future when warming continues.<\/p>\n<p>The Rahmstorf model and the Grinsted \u201chistorical\u201d case are by definition in excellent agreement with 20th Century data (and get similar values of <em>a<\/em>) since they have been tuned to those. The main difference arises from the differences between the two sea level data sets used: Church and White (2006) by Rahmstorf, Jevrejeva et al. (2006) by Grinsted et al. Since the \u201chistorical\u201d case of Grinsted et al. finds a ~1200-year response time scale, these two models are almost fully equivalent on a century time scale (e<sup>-100\/1200<\/sup>=0.92) and give nearly the same results. The total model rise in the last column is just 1.5 percent less than that based on the linear Eq. 3 because of that finite response time scale.<\/p>\n<p>For the Grinsted \u201cMoberg\u201d case the response time scale is only ~210 years, hence our linear approximation becomes bad already on a century time scale (e<sup>-100\/210<\/sup>=0.62, the total rise is 15% less than the linear estimate), which is why we give the linear estimates only in brackets for comparison here.<\/p>\n<p>The rise predicted by Siddall et al is much lower. That is not surprising, since their parameters were fitted to the slow changes of the big ice sheets (time scale \u03c4=2900 years) and don\u2019t \u201csee\u201d the early response caused by thermal expansion and mountain glaciers, which makes up most of the 20th Century sea level rise. What is surprising, though, is that Siddall et al. in their paper claim that their parameter values reproduce 20th Century sea level rise. This appears to be a calculation error (4); this will be resolved in the peer-reviewed literature. Our values in the above table are computed correctly (in our understanding) using the same parameters as used by the authors in generating their Fig.3. Their model with the parameters fitted to glacial-interglacial data thus underestimates 20th Century sea level rise by a factor of two.<\/p>\n<div align=\"center\"><font size=\"-1\"><br \/>\n<figure id=\"attachment_1009\" aria-describedby=\"caption-attachment-1009\" style=\"width: 400px\" class=\"wp-caption aligncenter\"><img decoding=\"async\" data-src=\"http:\/\/www.realclimate.org\/wp-content\/uploads\/icefront4001.jpg\" alt=\"Frosty legacy: We cannot afford to lose even a few percent of the land ice on Earth, which in total would be enough to raise global sea levels by 65 meters. (Calving front in Svalbard, (c) S.R.)\" title=\"icefront400\" width=\"400\" height=\"265\" class=\"size-full wp-image-1009 lazyload\" data-srcset=\"https:\/\/www.realclimate.org\/images\/icefront4001.jpg 400w, https:\/\/www.realclimate.org\/images\/icefront4001-300x198.jpg 300w\" data-sizes=\"(max-width: 400px) 100vw, 400px\" src=\"data:image\/svg+xml;base64,PHN2ZyB3aWR0aD0iMSIgaGVpZ2h0PSIxIiB4bWxucz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciPjwvc3ZnPg==\" style=\"--smush-placeholder-width: 400px; --smush-placeholder-aspect-ratio: 400\/265;\" \/><figcaption id=\"caption-attachment-1009\" class=\"wp-caption-text\">Frosty legacy: We cannot afford to lose even a few percent of the land ice on Earth, which in total would be enough to raise global sea levels by 65 meters. (Calving front in Svalbard, photo by S.R.)<\/figcaption><\/figure><\/font><\/div>\n<p><strong>Future projections<\/strong><\/p>\n<p>It thus looks like R07 and Grinsted et al. both reproduce 20th Century sea level rise and both get similar projections for the 21st Century. Siddall et al. get much lower projections but also strongly under-estimate 20th Century sea level rise. We suspect this will hold more generally: it would seem hard to reproduce the 20th Century evolution (including acceleration) but then get very different results for the 21st Century, with the basic semi-empirical approach common to these three papers. <\/p>\n<p>In fact, the lower part of their 7-82 cm range appears to be rather implausible. At the current rate, 7 cm of sea level rise since 2000 will be reached already in 2020 (see <a href=\"http:\/\/www.pik-potsdam.de\/~stefan\/update_science2007.html\">graph<\/a>). And Eq. 1 guarantees one thing for <em>any <\/em>positive value of <em>a<\/em>: if the 21st Century is warmer than the 20th, then sea level must rise faster. In fact the ratio of new sea level rise in the 21st Century to new sea level rise in the 20th Century according to Eq. 2 is not dependent on <em>a<\/em> or <em>b<\/em> and is simply equal to the ratio of the century-mean temperatures, T<sub>21<\/sub>\/T<sub>20<\/sub> (both measured again relative to the 1900 baseline). For the &#8220;coldest&#8221; IPCC-scenario (1.1 \u00baC warming for 2000-2100) this ratio is 1.3 \u00baC \/ 0.25 \u00baC = 5.2. Thus even in the most optimistic IPCC case, the linear semi-empirical approach predicts about five times the &#8220;new&#8221; sea level rise found for the 20th Century, regardless of parameter uncertainty. In our view, when presenting numbers to the <a href=\"\"http:\/\/www.bristol.ac.uk\/news\/2009\/6484.html\"\">public<\/a> scientists need to be equally cautious about erring on the low as they are on the high side. For society, after all, <em>under<\/em>-estimating global warming is likely the greater danger.<\/p>\n<p><strong>Does the world have to be linear?<\/strong><\/p>\n<p>How do we know that the relationship between temperature rise and sea level rate is linear, also for the several degrees to be expected, when the 20th century has only given us a foretaste of 0.7 degrees? The short answer is: we don&#8217;t. <\/p>\n<p>A slightly longer answer is this. First we need to distinguish two things: linearity in <em>temperature <\/em>(at a given point in time, and all else being equal), and linearity as the system evolves over <em>time<\/em>. The two are conflated in the real world, because temperature is increasing over time.<\/p>\n<p><em>Linearity in temperature<\/em> is a very reasonable assumption often used by glaciologists. It is based on a heat flow argument: the global temperature anomaly represents a heat flow imbalance. Some of the excess heat will go into slowly warming the deep ocean, some will be used to melt land ice, a tiny little bit will hang around in the atmosphere to be picked up by the surface station network. If the anomaly is 2 \u00baC, the heat flow imbalance should be double that caused by a 1 \u00baC anomaly. That idea is supported by the fact that the warming pattern basically stays the same: a 4 \u00baC global warming scenario basically has the same spatial pattern as a 2 \u00baC global warming scenario, only the numbers are twice as big (cf. <a href=\"http:\/\/www.ipcc.ch\/graphics\/ar4-wg1\/jpg\/spm6.jpg\">Figure SMP6<\/a> of the IPCC report). It\u2019s the same for the heating requirement of your house: if the temperature difference to the outside is twice as big, it will lose twice the amount of heat and you need twice the heating power to keep it warm. It\u2019s this \u201clinearity in temperature\u201d assumption that the Siddall et al. approach rejects.<\/p>\n<p><em>Linearity over time<\/em> is quite a different matter. There are many reasons why this cannot hold indefinitely, even though it seems to work well for the past 120 years at least. R07 already discusses this and mentions that glaciers will simply run out of ice after some time. Grinsted et al. took this into account by a finite time scale. We agree with this approach &#8211; we merely have some <a href=\"http:\/\/www.realclimate.org\/index.php\/archives\/2009\/01\/sea-will-rise-to-levels-of-last-ice-age\/ \">reservations <\/a>about whether it can be done with a single time scale, and whether the data they used really allow to constrain this time scale. And there are arguments (<a href=\"http:\/\/www.newscientist.com\/article\/mg19526141.600-huge-sea-level-rises-are-coming--unless-we-act-now.html\">e.g. by Jim Hansen<\/a>) that over time the ice loss may be faster than the linear approach suggests, once the ice gets wet and soft and starts sliding. So ultimately we do not know how much longer the system will behave in an approximately linear fashion, and we do not know yet whether the real sea level rise will then be slower or faster than suggested by the linear approach of Eq. 1.<\/p>\n<p><center><img decoding=\"async\" data-src=\"http:\/\/www.realclimate.org\/wp-content\/uploads\/Greenland600.jpg\" alt=\"Getting soft? Meltwater on the Greenland Ice Sheet. Photo by Ian Joughin.\" title=\"Greenland600\" src=\"data:image\/svg+xml;base64,PHN2ZyB3aWR0aD0iMSIgaGVpZ2h0PSIxIiB4bWxucz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciPjwvc3ZnPg==\" class=\"lazyload\" \/><br \/>\n<small>Getting soft? Meltwater lake and streams on the Greenland Ice Sheet near 68\u00baN at 1000 meters altitude. Photo by Ian Joughin.<\/small><\/center><\/p>\n<p><strong>Can paleoclimatic data help us?<\/strong><\/p>\n<p>Is there hope that, with a modified method, we may successfully constrain sea level rise in the 21st Century from paleoclimatic data? Let us spell out what the question is: How will sea level in the <em>present climate state<\/em> respond on a <em>century time scale<\/em> to a rapid <em>global warming<\/em>? We highlight three aspects here.<\/p>\n<p><em>Present climate state.<\/em> It is likely that a different climate state (e.g. the glacial with its huge northern ice sheets) has a very different sea level sensitivity than the present. Siddall et al. tried to account for that with their equilibrium sea level curve \u2013 but we think the final equilibrium state does not contain the required information about the initial transient sensitivity.<\/p>\n<p><em>Century time scale.<\/em> Sea level responds on various time scales \u2013 years for the ocean mixed layer thermal expansion, decades for mountain glaciers, centuries for deep ocean expansion, and millennia for big ice sheets. Tuning a model to data dominated by a particular time scale \u2013 e.g. the multi-century time scale of Grinsted et al. or the multi-millennia time scale of Siddall et al. \u2013 does not mean the results carry over to a shorter time scale of interest.<\/p>\n<p><em>Global warming.<\/em> We need to know how sea level \u2013 oceans, mountain glaciers, big ice sheets all taken together \u2013 responds to a globally near-uniform forcing (like greenhouse gas or solar activity changes). Glacial-interglacial climate changes are forced by big and highly regional and seasonal orbital insolation changes and do not provide this information. Siddall et al use a local temperature curve from Greenland and assume there is a constant conversion factor to global-mean temperature that applies across the ages and across different mechanisms of climate change. This problem is not discussed much in the paper; it is implicit in their non-dimensional temperature, which is normalised by the glacial-holocene temperature difference. Their best guess for this is 4.2 \u00baC (as an aside, our <a href=\"http:\/\/www.pik-potsdam.de\/~stefan\/Publications\/Journals\/schneider_etal_grl_2006.pdf\">published best guess<\/a> is 5.8 \u00baC, well outside the uncertainty range considered by Siddall et al). But is a 20-degree change in Greenland temperature simply equivalent to a 4.2-degree global change? And how does local temperature translate into a global temperature for <a href=\"http:\/\/www.realclimate.org\/index.php\/archives\/2006\/11\/revealed-secrets-of-abrupt-climate-shifts\/\">Dansgaard-Oeschger events<\/a>, which are generally assumed to be caused by ocean circulation changes and lead to a temperature <a href=\"http:\/\/www.sciencemag.org\/cgi\/content\/short\/282\/5386\/61\">seesaw effect<\/a> between northern and southern hemisphere? What if we used their amplitude to normalise temperature \u2013 given their imprint on global mean temperature is approximately zero?<\/p>\n<p>Overall, we find these problems extremely daunting. For a good constraint for the 21st Century, one would need sufficiently accurate paleoclimatic data that reflect a sea level rise (a drop would not do \u2013 ice melts much faster than it grows) on a century time scale in response to a global forcing, preferably from a climate state similar to ours \u2013 notably with a similar distribution of ice on the planet. If anyone is aware of suitable data, we\u2019d be most interested to hear about them!<\/p>\n<p><strong>Update<\/strong> (8 Sept): We have now received the computer code of Siddall et al (thanks to Mark for sending it). It confirms our analysis above. The code effectively assumes that the warming over each century applies for the whole century. I.e., the time step for the 20<sup>th<\/sup> Century assumes the whole century was 0.74 \u00baC warmer than 1900, rather than just an average of 0.25 \u00baC warmer as discussed above. When this is corrected, the 20<sup>th<\/sup> Century rise reduces from 15 cm to 8 cm in the model (consistent with our linear estimate given above). The 21st Century projections ranging from 32-48 cm in their Table 1 (best estimates) reduce to 24-32 cm.<\/p>\n<p><em>Martin Vermeer is a <a href=\"http:\/\/www.iag-ggos.org\/\">geodesist<\/a> at the Helsinki University of Technology in Finland.<\/em><\/p>\n<p><strong>Footnotes<\/strong><\/p>\n<p><strong>(1)<\/strong> Siddall et al. use two steps. First they determine an equilibrium sea level for each temperature (their Eq 1, and shown in their Fig. 1). Second, they assume an exponential approach of sea level to this equilibrium value in their Eq. 2, which (slightly simplified, for the case of rising sea level) reads:<\/p>\n<p>dS\/dt = (S<sub>e<\/sub>(T) \u2013 S(t)) \/ \u03c4.<\/p>\n<p>Here S is the current sea level (a function of time <em>t<\/em>), S<sub>e<\/sub> the equilibrium sea level (a function of temperature <em>T<\/em>), and \u03c4 the time scale over which this equilibrium is approached (which they find to be 2900 years).<br \/>\nNow imagine the temperature rises. Then S<sub>e<\/sub>(T) increases, causing a rise in sea level dS\/dt. If you only look at short time scales like 100 years (a tiny fraction of those 2900 years response time), S(t) can be considered constant, so the equation simplifies to<\/p>\n<p>dS\/dt = S<sub>e<\/sub>(T)\/ \u03c4 + constant.<\/p>\n<p>Now S<sub>e<\/sub>(T) is a non-linear function, but for small temperature changes (like 1 \u00baC) this can be approximated well by a linear dependence S<sub>e<\/sub>(T) = s * T + constant. Which gives us<\/p>\n<p>dS\/dt = s\/\u03c4 * T + constant, i.e. Eq (1) in the main post above.<\/p>\n<p>R07 on the other hand used:<br \/>\ndS\/dt = <em>a<\/em> * (T \u2013 T0), which is also Eq. (1) above.<br \/>\nNote that <em>a <\/em>= <em>s<\/em>\/\u03c4 and <em>b<\/em> = &#8211;<em>a<\/em>*T0 in our notation.<\/p>\n<p><strong>(2)<\/strong> Here is a very basic matlab\/<a href=\"http:\/\/www.gnu.org\/software\/octave\/\">octave<\/a> script that computes a sea level curve from a given temperature curve according to Eq. 2 above. The full matlab script used in R07, including the data files, is available as <a href=\"http:\/\/www.sciencemag.org\/cgi\/content\/full\/sci;317\/5846\/1866d\/DC1\">supporting online material<\/a> from <em>Science <\/em><br \/>\n<tt><br \/>\n% Semi-empirical sea level model - very basic version<br \/>\nT1900=mean(tempg(11:30)); T=tempg-T1900;<\/p>\n<p>a=0.34;  % sea level sensitivity parameter [cm\/degree\/year]<br \/>\nb=0.077; % note this value depends on a and on the temperature<br \/>\n         % baseline, here the mean 1890-1909<\/p>\n<p>% rate of rise - here you need to put in an annual temperature time series T<br \/>\n% with same baseline as chosen for fitting b!<br \/>\ndSdt = a*T + b;<\/p>\n<p>% integrate this to get sea level over the period covered by the temperature series<br \/>\nS = cumsum(dSdt); plot(S);<\/tt><\/p>\n<p><strong>(3)<\/strong> Here is a matlab\/<a href=\"http:\/\/www.gnu.org\/software\/octave\/\">octave<\/a> script to compute the equilibrium sea level curve of Siddall et al. Note the parameters differ in some cases from those given in the paper \u2013 we obtained the correct ones from Mark Siddall.<br \/>\n<tt><br \/>\n% Siddall et al equilibrium sea level curve, their Fig. 1, NGRIP scenario<br \/>\nA = 15.436083479092469;<br \/>\nb = 0.012630000000000;<br \/>\nc = 0.760400212014386;<br \/>\nd = -73.952809369848552;<\/p>\n<p>Tdash=[-1.5:.05:2];<br \/>\n% Equilibrium sea level curve<br \/>\nSe=A*asinh((Tdash+c)\/b) + d;<br \/>\n% Tangent at current temperature<br \/>\ndSe=A\/sqrt(1+((0+c)\/b)^2)\/b;<br \/>\nSe0= A*asinh((0+c)\/b) + d;<br \/>\nTe=dSe*Tdash + Se0;<br \/>\nplot(Tdash, Se, 'b', Tdash, Te, 'c', Tdash, 0.0*Se, 'k', [0 0], [-150 40], 'k')<br \/>\nxlabel('Dimensionless temperature')<br \/>\nylabel('Equilibrium sea level (m)')<br \/>\nfprintf(1, 'Slope: %f m\/K, Sensitivity: %f cm\/K\/year, zero offset: %f m\\n\\n', dSe\/4.2, 100*dSe\/4.2\/2900, Se0);<br \/>\n<\/tt><\/p>\n<p><strong>(4)<\/strong> We did not yet receive the code at the time of writing, but based on correspondence with the authors conclude that for their values in Fig. 3 and table 1, Siddall et al. integrated sea level with 100-year time steps with a highly inaccurate numerical method, thus greatly overestimating the <em>a<\/em>-term. In their supporting online information they show a different calculation for the 20th Century with annual time steps (their Fig. 5SI). This is numerically correct, giving an <em>a<\/em>-term of about 4 cm, but uses a different value of <em>b<\/em> close to 0.12 cm\/year to obtain the correct total 20th Century rise.<\/p>\n<p><strong>References<\/strong><\/p>\n<p><a href=\"http:\/\/www.pol.ac.uk\/psmsl\/author_archive\/church_white\/GRL_Church_White_2006_024826.pdf\">Church, J. A. &#038; White, N. J. A 20th century acceleration in global sea-level rise. Geophysical Research Letters 33, L01602 (2006).<br \/>\n<\/a><\/p>\n<p><a href=\"http:\/\/www.glaciology.net\/Home\/PDFs\/grinstedclimdyn09sealevel200to2100ad.pdf\">Grinsted, A., Moore, J. C. &#038; Jevrejeva, S. Reconstructing sea level from paleo and projected temperatures 200 to 2100 ad. Climate Dynamics (2009).<\/a><\/p>\n<p><a href=\"http:\/\/www.glaciology.net\/Home\/PDFs\/Jevrejeva_JGR06_-_sea_level_trends_-virtual_station.pdf\">Jevrejeva, S., Grinsted, A., Moore, J. C. &#038; Holgate, S. Nonlinear trends and multiyear cycles in sea level records. Journal of Geophysical Research 111 (2006).<\/a><\/p>\n<p><a href=\"http:\/\/www.nature.com\/nature\/journal\/v433\/n7026\/full\/nature03265.html\">Moberg, A., Sonechkin, D. M., Holmgren, K., Datsenko, N. M. &#038; Karlen, W. Highly variably Northern Hemisphere temperatures reconstructed from low- and high-resolution proxy data. Nature 433, 613-617 (2005).<\/a><\/p>\n<p><a href=\"http:\/\/www.pik-potsdam.de\/~stefan\/Publications\/Nature\/rahmstorf_science_2007.pdf\">Rahmstorf, S. A semi-empirical approach to projecting future sea-level rise. Science 315, 368-370 (2007).<\/a><\/p>\n<p><a href=\"http:\/\/www.pik-potsdam.de\/~stefan\/Publications\/Nature\/rahmstorf_science_response_2007.pdf\">Rahmstorf, S. Response to comments on &#8220;A semi-empirical approach to projecting future sea-level rise&#8221;. Science 317 (2007).<\/a><br \/>\n<a href=\"http:\/\/www.geo.oregonstate.edu\/people\/faculty\/publications\/clarkp\/Siddall-2009-NatureGeo.pdf\"><br \/>\nSiddall, M., Stocker, T. F. &#038; Clark, P. U. Constraints on future sea-level rise from past sea-level change. Nature Geoscience (advance online publication, 26 July 2009).<\/a><\/p>\n<!-- kcite active, but no citations found -->\n<\/div> <!-- kcite-section 969 -->","protected":false},"excerpt":{"rendered":"<p>By Stefan Rahmstorf and Martin Vermeer The scientific sea level discussion has moved a long way since the last IPCC report was published in 2007 (see our post back then). The Copenhagen Synthesis Report recently concluded that \u201cThe updated estimates of the future global mean sea level rise are about double the IPCC projections from [&hellip;]<\/p>\n","protected":false},"author":10,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_exactmetrics_skip_tracking":false,"_exactmetrics_sitenote_active":false,"_exactmetrics_sitenote_note":"","_exactmetrics_sitenote_category":0,"_genesis_hide_title":false,"_genesis_hide_breadcrumbs":false,"_genesis_hide_singular_image":false,"_genesis_hide_footer_widgets":false,"_genesis_custom_body_class":"","_genesis_custom_post_class":"","_genesis_layout":"","footnotes":""},"categories":[19],"tags":[],"class_list":{"0":"post-969","1":"post","2":"type-post","3":"status-publish","4":"format-standard","6":"category-oceans","7":"entry"},"aioseo_notices":[],"post_mailing_queue_ids":[],"_links":{"self":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts\/969","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/users\/10"}],"replies":[{"embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/comments?post=969"}],"version-history":[{"count":93,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts\/969\/revisions"}],"predecessor-version":[{"id":1701,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts\/969\/revisions\/1701"}],"wp:attachment":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/media?parent=969"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/categories?post=969"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/tags?post=969"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}