The Impact of Uncertainties in Ice Sheet Dynamics on Sea-Level Allowances at Tide Gauge Locations

<jats:p>Sea level is projected to rise in the coming centuries as a result of a changing climate. One of the major uncertainties is the projected contribution of the ice sheets in Greenland and Antarctica to sea-level rise (SLR). Here, we study the impact of different shapes of uncertainty distributions of the ice sheets on so-called sea-level allowances. An allowance indicates the height a coastal structure needs to be elevated to keep the same frequency and likelihood of sea-level extremes under a projected amount of mean SLR. Allowances are always larger than the projected SLR. Their magnitude depends on several factors, such as projection uncertainty and the typical variability of the extreme events at a location. Our results show that allowances increase significantly for ice sheet dynamics’ uncertainty distributions that are more skewed (more than twice, compared to Gaussian uncertainty distributions), due to the increased probability of a much larger ice sheet contribution to SLR. The allowances are largest in regions where a relatively small observed variability in the extremes is paired with relatively large magnitude and/or large uncertainty in the projected SLR, typically around the equator. Under the RCP8.5 (Representative Concentration Pathway) projections of SLR, the likelihood of extremes increases more than a factor     10 4     at more than 50–87% of the tide gauges.</jats:p>

section_content: Sea level is projected to rise in the coming century and beyond as a result of a warming climate, with major contributions from warming of the ocean water and increasing land ice mass loss [1]. To protect against or mitigate the effects of sea-level rise (SLR) in highly populated coastal zones, it is important to understand the processes causing SLR, and to reduce the uncertainties for each of the projected contributions. Sea-level rise here is relative sea-level change, i. e., the difference as measured by tide gauges attached to the Earth's surface. Presently, one of the main uncertainties in SLR is the contribution of the Antarctic and Greenland ice sheets [2, 3]. Particularly the timing and magnitude of dynamical changes associated with the ice sheets (e. g., marine ice sheet instability) are still uncertain. This is therefore an active area of research. 

In the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC), the projected ice dynamical contribution to SLR was based on observations and a range of scenarios from models ( [1], Table 13. 5). Around the same time, Bamber and Aspinall [4] presented an estimate based on expert judgement, which had a strongly skewed uncertainty distribution rather than the symmetric distribution presented in IPCC. With the same data but a different way to process the expert opinions, De Vries and Van de Wal [5] presented an estimate with a lower high-end estimate than Bamber and Aspinall [4], though still a strongly skewed uncertainty distribution. Although expert elicitations document the level of consensus on a certain topic, they are sensitive to, for instance, the level of expertise of the interviewees and the way experts are selected. Therefore, process-based estimates are generally a preferable method. Until recently, process-based estimates were only available for specific parts of the ice sheets (e. g., Pine Island Glacier [6] or Thwaites Glacier [7] ). However, Ritz et al. [8] recently presented an estimate for the Antarctic ice dynamical contribution to SLR in the 21st century using a process-based approach. They too found a non-normal distribution of the uncertainties, though less skewed than the De Vries and Van de Wal [5] distribution. 

In this study, we will compare and contrast three different types of uncertainty distributions for the mass loss contribution to SLR due to ice dynamics [1, 5, 8]. We will specifically look at the effect on sea-level allowances rather than focusing on the sea-level projections. A sea-level allowance indicates the height a coastal structure needs to be elevated to maintain the same frequency and likelihood of extreme sea-level events under future SLR [9, 10]. An allowance gives more information than just a 90% or 95% confidence level of projected mean SLR because it also accounts for the local variability in the sea-level extremes. An important assumption for the allowances in general is that the variability in the extremes does not change under the influence of climate change. For the 20th century, this is a reasonable assumption, as Menéndez and Woodworth [11] found that the historical change in extremes is mostly due to changes in mean sea level rather than changes in the storminess contributions to the extremes. For the 21st century, model projections of changes in the wave climate and storm surges are subject to large uncertainties, and it is therefore difficult to say if the assumption will hold in the future [1]. 

We use the regional relative sea-level projections of Slangen et al. [12], which describe the change of the ocean surface relative to the ocean floor. The projections follow the high-end RCP8. 5 (Representative Concentration Pathway [13] ) climate change scenario based on CMIP5 model output (5th phase of the Climate Model Intercomparison Project [14] ). While these projections assume a normal (Gaussian) uncertainty distribution, more recent publications have incorporated skewed ice sheet dynamics distributions into their projections to provide a more complete uncertainty distribution (e. g., [15] [16] [17] [18] ) which can be used for risk-assessment and policy purposes. We follow this line and combine the Slangen et al. [12] sea-level projections with three different ice sheet dynamics scenarios (Section 2. 1) using the framework of De Winter et al. [18]. Whereas De Winter et al. [18] discussed the impact and sensitivity of the high-end regional sea-level projections to different ice sheet uncertainty distributions, here we focus on the effect of the shape of the skewed ice sheet uncertainty distributions on the projected allowances. De Winter et al. [18] showed that the 90th, 95th and 97. 5th percentiles of regional SLR projections are highly sensitive to the uncertainty distribution of ice sheet dynamics. In contrast to De Winter et al. [18], all the ice sheet dynamics scenarios will be shifted to the same median (as described in Section 2. 1) to focus on the effects of the different shapes rather than differences in their magnitude. Previously, regional allowances were mainly computed using normal (Gaussian) uncertainty distributions for the projected SLR [10]. Buchanan et al. [19] used full probability distributions in combination with allowances to compute an "Average Annual Design-Life Level" for the United States coastline. Here, we present allowances for skewed uncertainty distributions for tide gauge stations across the world and compare the different shapes of the distributions (Section 2. 3). 

There are various ways to express sea-level extreme events statistically, depending on the application and the quality of the available data. Methods such as fitting Generalised Extreme Value distributions (of which the Gumbel distribution is a special case) to annual maxima in water levels and fitting Generalised Pareto distributions to peaks over a threshold [20] are often used to represent the statistics of extreme water levels. There is no consensus as to a universal best approach [21], and the method chosen depends on the purpose of the analysis and the data availability and quality. Here, we use a Gumbel parameterisation, which will allow us to compute the allowances following Hunter [9] and Hunter et al. [10]. We use the GESLA-2 tide gauge dataset (Global Extreme Sea Level Analysis Version 2 [22, 23], Section 2. 2), which significantly increases the number of locations with respect to previous studies of allowances [9, 10]. Finally, we show the increase in the frequency of sea level extremes under the RCP8. 5 scenario if the allowances are not applied (Section 3. 3).
section_name: Introduction

section_name: Data and Methodology

section_content: We use the regional sea-level projections for the RCP8. 5 scenario based on CMIP5 model output as presented in Slangen et al. [12], for the period 2010-2100. These projections are similar to the results presented in IPCC AR5 [1] and include the contributions of steric/dynamic change, glacier and ice sheet mass change, groundwater extraction and Glacial Isostatic Adjustment. The global mean projected change for 2010-2100 is 0. 72 ± 0. 28 m, of which the Antarctic dynamics mass loss contribution is 0. 08 ± 0. 06 m and the Greenland dynamics mass loss contribution 0. 05 ± 0. 06 m. All uncertainty distributions in Slangen et al. [12] are assumed to be normal, with means and standard deviations as provided in their Table 1. 

The two dynamical ice sheet contributions will be removed from the total projections of Slangen et al. [12], and replaced by three different scenarios for the ice sheet dynamically-driven mass loss. For clarity, this means that the ice sheet surface mass balance (SMB) contribution to SLR remains unaltered with respect to Slangen et al. [12]. The three scenarios used are:

(1) The IPCC scenario (Figure 1, green, from Church et al. [1] ). We use the cumulative ice sheet dynamical contributions for the Antarctic Ice Sheet (AIS) and Greenland Ice Sheet (GRIS) as presented in IPCC AR5 Table 13. 5 (0. 05 ± 0. 03 m for Greenland, 0. 07 ± 0. 08 for Antarctica), assuming a normal distribution. As Church et al. [1] present the change for 1986-2005 to 2081-2100, we correct the values linearly to estimate the change for 2010 to 2100. The Antarctic contribution is split into 11% for the East Antarctic Ice Sheet (EAIS) and 89% for the West Antarctic Ice Sheet (WAIS), which is the split for the year 2010 as presented in the IMBIE project (Ice sheet Mass Balance Inter-comparison Exercise, [24] ). 

(2) The VW15 scenario (Figure 1, blue) is from De Vries and Van de Wal [5]. We use the skewed uncertainty distribution of rates of SLR in 2100 for WAIS, EAIS and GRIS contributions, based on expert elicitation estimates. VW15 presented a reanalysis of the data in Bamber and Aspinall [4], which more rigorously accounts for the lack of consensus in the expert estimates. To convert mm/year in 2100 into cumulative SLR for 2010-2100, we assume a rate of change in 2010 of 0. 6 mm/year for WAIS, 0. 0 mm/year for EAIS and 0. 3 mm/year for GRIS following VW15 and Bamber and Aspinall [4], and assume a linear change in the rate for the period 2010-2100. 

(3) The R15 scenario for the Antarctic dynamical ice sheet contribution (Figure 1, black) is from Ritz et al. [8]. This is an observation and model-based estimate of the ice dynamical mass loss for 2010-2100. The shape of the Antarctic distributions differs slightly from the original R15 study due to the use of a narrower bandwidth (less smoothing) in the kernel density estimation. For Greenland dynamics, this scenario is the same as the IPCC scenario, as R15 did not calculate the Greenland contribution. 

All the ice dynamics scenarios are shifted to match the median values of the IPCC scenario (Figure 1, dashed lines), in order to make sure that the differences in the allowances are only due to changes in the shape of the distribution rather than the magnitude. We use the median (i. e., the point where the cumulative density function is 0. 5) for this shift, rather than the mean, as the median is less sensitive to extremely large values in the non-normal distributions. Although the published VW15 estimates are for SMB and ice sheet dynamics combined, we here assume that the skewed shape of the distribution is primarily a result of the uncertainty in the dynamical processes and correct for the differences in magnitude by the aforementioned shift in the median. 

The regional SLR for each of the ice sheet dynamics contributions is computed by multiplying each contribution as presented in Figure 1 (Greenland, East Antarctica, West Antarctica) with its respective sea-level "fingerprint" ratio (e. g., [25] ). These fingerprints are computed with a gravitationally consistent sea-level model. The resulting projected ice dynamical contribution to the relative regional SLR was shown in Slangen et al. [12] in their Figure 1e. Typically, sea level falls near the sources of ice mass loss, and rises above average in the far field. The individual fingerprint ratios of each ice sheet may be scaled with the global mean SLR contribution of each ice sheet if the location of the mass change remains the same as the location used to compute the fingerprint. This assumption is most significant in proximity to the ice sheet, where gravitational gradients decrease relatively sharply from the centres of mass loss. We use the SEAWISE model [18] to combine the skewed SLR distributions for Greenland, West-Antarctica and East-Antarctica one by one with the projected non-ice sheet dynamics SLR from Slangen et al. [12] (for which a normal distribution is assumed) at each tide gauge location (Appendix, Figure A1 ). SEAWISE first scales the ice dynamic uncertainty distributions with the individual fingerprints to obtain the ice sheet uncertainty distribution at each location. Then, SEAWISE combines multiple SLR distributions, e. g., for West-Antarctica and Greenland (these can be normal or skewed). The model can combine independent distributions or assign a correlation coefficient to assume a dependency between distributions. Here, all distributions are assumed to be independent. Details of the SEAWISE methodology are presented in the Appendix and in De Winter et al. [18].
section_name: Uncertainty Distributions of Sea-Level Change

section_content: In addition to the regional sea-level projections, we require location-specific information on the statistics of sea-level extremes resulting from tides, storm surges and other (high-frequency) sea-level variability. For this, we use Gumbel scale parameters. The Gumbel scale parameter for a certain location can be derived from historical records of high-frequency sea-level measurements (see Figure 2 for four examples). Here, we use data from the GESLA-2 database (Global Extreme Sea Level Analysis Version 2 [22, 23] ), which contains high-frequency records (at least hourly) from 1300 tide gauge stations as of February 2016. Gumbel distributions were fitted to the annual maxima (expressed relative to the linear trend in mean sea level over the record) of each of the 658 tide gauges, using the Matlab R evfit (Mathworks, Natick, MA, USA) function, as described in Hunter et al. [26]. In a comparison of the Gumbel scale parameters obtained from the evfit method to the ismev method as used in Hunter [9] (both methods are extensively described in Hunter et al. [26] ), we found that the results are very similar to the extent required for the present analysis and we therefore use the values from the Matlab R evfit method only. We select all tide gauge records that contain at least 20 years of data and where each year with data is more than 75% complete (Figure 3 ), which yields a data set of 658 records. Following Hunter et al. [26], we reject four stations that show significant non-Gumbel behaviour, and also leave out two Hudson Bay tide gauges (Canada) as we do not have sea-level projections available there. This leaves 652 records, of which 448 records have ≥30 years of data available, 319 records ≥40 years, 164 records ≥50 years and 94 records ≥60 years. The Gumbel scale parameter describes the variability of the extremes; i. e., a smaller parameter means that the variability in the extremes is small (e. g., Figure 2a, c ) -this does not necessarily mean that the average magnitudes of the extremes are small. A larger scale parameter indicates larger variability in the extremes (e. g., Figure 2b, d ). Sea-level extremes can be a combination of high tide and extreme storm surge, but also a medium tide and an extreme surge, or a high tide and a medium surge. Other factors such as interannual or seasonal variability can have significant effects as well: Barcelona (Figure 2a ), for instance, has a strong seasonal cycle, leading to a more regular occurrence of the annual maximum than for instance in Delfzijl (Figure 2b ). For a detailed discussion of contributing factors to extreme sea levels, we refer the reader to Merrifield et al. [27]. 

Historically, there have been more tide gauge stations in the Northern Hemisphere, resulting in a larger number of estimates of Gumbel scale parameters in the Northern Hemisphere than in the Southern Hemisphere. The Gumbel scale parameters are typically in the range of tenths of meters (0. 12 ± 0. 06 m, mean ± 1σ). Values are generally lower around the equator, increasing towards the mid-latitudes (specifically in the Northern Hemisphere), with the highest values between 40 o and 70 o N (Figure 3 ), where both the tidal range and the storm surge activity tend to be larger [27]. For instance, Gumbel scale parameters are high at the Dutch coast (e. g., Figure 2b, Delfzijl tide gauge), due to a large tidal range and a large variability in storm surge heights. This is the result of the geometry of the semi-enclosed and relatively shallow North Sea basin, resulting in the development of larger storm surges under north western storm wind conditions, while the storm surge levels are considerably lower for other wind directions (e. g., [28] ). This in contrast to the Barcelona tide gauge (Figure 2a ), which has a much smaller tidal range since it is located in the Mediterranean and experiences less extreme storm surge heights.
section_name: The Statistics of Sea-Level Extremes

section_content: We assume a projected sea-level rise by an amount z, where z has an uncertainty distribution function P(z ). If a coastal structure is raised by an allowance a to preserve the frequency of extreme (flooding) events in the future, a is given by

where λ is the scale parameter of the Gumbel distribution (Section 2. 2), which describes the sea level extreme events and which is assumed to remain unchanged under sea-level rise. For the full derivation of the allowance, the reader is referred to Hunter [9] and Hunter et al. [10]. Note that, in these earlier papers, z was defined relative to a projected mean sea level, whereas here we use z as the total projected sea level to simplify the equations. For a normal distribution, the allowance can be reduced to a simple analytical form:

(Equation viii in Supplementary Material of Hunter [9] ), where µ is the mean SLR and σ the standard deviation of SLR when P(z ) follows a normal distribution. However, Equation ( 2 ) is only valid for simple normal distributions of SLR. As we want to test different types of (non-normal) uncertainty distributions, the approach needs to be modified. We first redefine the skewed distributions of the projected SLR (from Section 2. 1) in such a way that they can be used in the allowance computation. This is done by using a nonlinear fitting procedure (Powell's method, Press et al. [29] ), which fits up to four normal distributions (see Figure 4 for an example) to optimally represent the (skewed) sea-level distribution. Then, the equation for the allowances only needs to be modified such that the allowance represents this group of normal distributions rather than just one normal distribution. We start from the description of a single normal distribution with a zero mean and standard deviation σ, where P(z ) is expressed as:

(Eq. ii from Supplementary Material of Hunter [9] ). Then, the distribution P(z ) is modified such that it represents a group of normal distributions rather than one normal distribution:

where each distribution i in the group of normals has a weight w i, a standard deviation σ i, and is centered on µ i, such that P(z ) is described by a group of n normal distributions. The sum of all weights w equals 1, where individual weights may be negative. Now, Equation ( 4 ) is inserted in Equation (1):

This results in the following definition of the allowance for a group of normal distributions:

This equation will be used to compute allowances for non-normal distributions in the remainder of the paper. For one normal distribution (n = 1) with weight 1, Equation ( 6 ) reduces again to Equation (2).
section_name: Allowances Methodology

section_name: Results

section_content: Sea-level projections are computed for each grid point of the RCP8. 5 projections of Slangen et al. [12] in conjunction with the regional patterns of the three scenarios for ice sheet dynamics (Section 2. 1), which leads to three sea-level projection scenarios (Figure 5 ). In the three scenarios, the median values (i. e., the point where the cumulative density function reaches 0. 5) of the uncertainty distributions of WAIS, EAIS and GRIS are set to match the IPCC median values (Figure 1, dashed lines) in order to study only the effect of the change in the shape of the uncertainty distributions, rather than the shifts in magnitude. This is a different approach compared to the work presented in De Winter et al. [18], who did not shift the medians, as the purpose of that study was to look at the sensitivity of the high-end percentiles of regional sea-level projections to different ice sheet scenarios, considering both the shape and the magnitude of the distributions. 

As a result of the differences in the shape of the ice sheet dynamics distributions, the final median SLR is not exactly the same for the three scenarios, even though the medians of the individual contributions are the same. Therefore, the VW15 scenario has a higher median (over the global ocean area) of 0. 81 m compared to 0. 71 m for IPCC and 0. 73 m for R15 (Figure 5, left column). The maps show low or negative SLR around the glaciers and ice sheets. As a result of projected changes in the ocean heat content and circulation, there is a dipolar band associated with a shift in the Antarctic Circumpolar Current stretching south of South Africa and Australia [30] and generally higher values in a wide band around the equator in the Indian and Pacific Oceans. 

The 95th percentile figures of the sea-level projections (Figure 5, right column) show significantly larger values for the VW15 scenario compared to the other two scenarios, reaching projected SLR over 1. 75 m. The global mean 95th percentile of the VW15 scenario is 1. 56 m, which is higher than the other two scenarios (1. 16 m IPCC and 1. 18 m R15). This is a result of the skewness of the VW15 dynamical ice sheet uncertainty distribution, which has a heavier tail that leads to much higher values for the 95th percentile of all three dynamical ice sheet contributions in the VW15 scenario (Figure 1 ). Although the R15 scenario has a higher 95th percentile contribution for the West-Antarctic ice dynamics (+0. 12 m compared to IPCC) and a slightly longer tail to higher values, the shape of the distribution is less skewed than that of the VW15 scenario. The R15 East Antarctic 95th percentile contribution is smaller (-0. 04 m compared to IPCC) due to the distribution being much sharper than the IPCC distribution. Combined, the 95th percentile total SLR of R15 is only marginally larger than that of IPCC (+0. 02 m).
section_name: Sea-Level Projections for 2010-2100

section_content: From the sea-level projections and the Gumbel scale parameters, we compute sea-level allowances for the year 2100 (with respect to 2010) as outlined in Section 2. 3. The allowances can only be computed at tide gauge locations, where Gumbel scale parameters are available. The projected SLR for each location is taken from the regional sea-level maps using a nearest neighbour approach. In some regions with multiple tide gauge stations, this means that the same sea-level projection grid point will be used for different stations. 

Similar to the 95th percentile sea-level projections (Figure 5 ), the allowances show a distinct effect from the skewed distribution of the VW15 scenario, and less so from the more mildly skewed R15 distribution (Figure 6, left column). On average, the VW15 allowances are more than twice as large as the allowances in the other two scenarios (2. 19 m vs. 1. 01 (IPCC) and 1. 03 (R15) m). For the normal and close to normal distributions, it appears that the allowance is slightly less than the 95th percentile SLR at the tide gauge locations (which is 1. 08 m for IPCC and 1. 12 m for R15). For the skewed VW15 distribution, the allowance is much larger than the 95th percentile SLR at the tide gauges (1. 50 m), due to the heavy upper tails of the uncertainty distributions of the dynamical ice sheet contributions. When we compare the allowance to the global mean SLR (Figure 6, middle column) and to the local projected SLR (Figure 6, right column), this shows that an allowance is more than a simple addition of the change in the mean sea level: it accounts for the local distribution of extremes through the Gumbel scale parameter as well. In some regions, the projected local SLR is less than the global mean change in the IPCC and R15 scenarios, for instance in the Baltic Sea and the northwest of the United States. Although the allowance may be less than the global mean SLR in some places, it is always larger than the local SLR (compare middle to right columns, Figure 6 ). The largest differences between the projected SLR (global and regional) and the allowances are in and along the Central Pacific ocean basin, where Gumbel scale parameters (Figure 3 ) tend to be lower (indicating smaller variability in the extremes) and projected SLR and its uncertainty tend to be higher (due to large distance from ice mass loss sources and large steric/dynamic changes). 

A focus on Europe (Figure 7 ) shows that the allowances in this region fall between 0. 5 and 1 m for the IPCC and R15 scenarios, while the allowances are larger for the VW15 scenario, being at least 1 m and as much as 2. 5 m. The allowances pattern generally is opposite to the Gumbel scale parameter (Figure 3 ): regions with higher Gumbel scale parameters tend to have lower allowances and vice versa. Due to Glacial Isostatic Adjustment (ongoing response of the solid Earth following ice sheet unloading after the Last Glacial Maximum [31] ), the projected SLR is lower in the Baltic, leading to lower allowances in all three scenarios. Around the Mediterranean basin, projected SLR is substantial and the present sea-level variability is relatively small (small Gumbel scale parameters), which leads to the highest allowances in Europe (Figure 2a ). 

There is no immediately obvious relation between the Gumbel scale parameter and the projected mean SLR, apart from a tendency for locations with a smaller Gumbel scale parameter (<0. 1 m) to be in locations with an average to above-average SLR projection, for all three scenarios (Figure 8 ). However, if the allowance is included (colours in Figure 8 ), a pattern emerges in the allowance that results from the combination of the projections and the scale parameters. It shows that tide gauges with smaller Gumbel scale parameters and larger projected SLR (uncertainties) require larger allowances to maintain the same likelihood of flooding (following Equation ( 6 )).
section_name: Allowances for Different Uncertainty Distributions

section_content: The frequency of sea-level extremes at any given location will change when mean sea level changes. Looking at Figure 2, the same amount of SLR would have a larger effect on the number of extreme events in places with small variability (both in sense of total range of variability and in the small variability in the extremes, such as Barcelona and Gan, Figure 2a, c ). In these places, even 0. 5 m of SLR would bring the mean sea level outside the range of current sea-level variability. In places such as Delfzijl and Willapa Bay (Figure 2b, d ), which have a larger range already, a mean SLR of 0. 5 m would not be out of the ordinary range of variability. However, it would mean that the extremes, which are already of considerable magnitude, would become even larger. 

The sea-level allowance (Equation ( 6 )) is designed such that a location maintains the same frequency of extremes when the infrastructure is raised by the amount of the allowance. To quantify the effect of not applying the allowance (i. e., the infrastructure is not moved), we compute the change in the frequency of extremes N following

where a is the (skewed or normal) allowance and λ the Gumbel scale parameter. If the frequency of the extremes increases by a factor 500 this means that a 1 in 10,000 year extreme event will by 2100 occur once every 20 years. For a factor 5, the 1 in 10,000 year event would become a 1 in 2000 year event. It is therefore important to consider regional sea-level projections and their uncertainties, and sea-level extreme statistics, which all differ significantly between different locations. The changing frequencies of extreme events for the IPCC and R15 scenarios (Figure 9 ), computed using Equation (7), show larger increases than in (Church et al. [1] their Figure 13. 25), due to the use of the RCP8. 5 scenario here rather than the RCP4. 5 scenario in Figure 13. 25b. Another difference with respect to Church et al. [1] is that there are more stations with updated time series included in Figure 9, which can lead to small variations in the Gumbel scale parameter λ. Since the allowances for the VW15 scenario are on average twice as large as the allowances in the other two scenarios due to the skewed uncertainty distribution, the increase in frequencies is even larger due to the exponent in Equation (7) (Figure 9, Table 1 ). Overall, the largest increases in frequencies are projected in the Equatorial Pacific and along the Southern and Eastern coasts of the U. S., where projected SLR and uncertainties are relatively large and Gumbel scale parameters are small. If Equation (7) results in ratios of more than 10 4, this means that a 1 in 10,000 year event would occur at least once every year by 2100, and a 1 in 100 year event would occur at least 100 times per year by 2100. Factors over 10 4 therefore effectively mean that a coastal structure, designed for any normal range of return periods (i. e., 100-10,000 years), would be lost unless the allowance is applied. For IPCC and R15, we find factors over 10 4 for 50% and 52% of the tide gauge locations, respectively, while, for VW15, this increases to 87% due to the skewed ice sheet contribution (Table 1 ).
section_name: Changes in the Frequency of Extreme Events

section_content: In computing the allowances, several choices and assumptions were made for the projected mean SLR and for the representation of extreme sea levels. 

Regarding the projections, we used the RCP8. 5 scenario as presented in Slangen et al. [12]. The RCP8. 5 scenario is chosen because it is currently the more realistic option among the available scenarios, which is therefore relevant for coastal protection. Compared to the projections presented in [1], there are small differences (for the period 1986-2005 vs. 2080-2100, RCP8. 5 scenario) in the individual contributions: glaciers (+6 cm), Greenland SMB (-1 cm), Antarctic SMB (-4 cm), Greenland Dynamics (+1 cm), Antarctic Dynamics (+2 cm) and groundwater depletion (+4 cm). The main reason is that Slangen et al. [12] used one model or approach for each contribution to sea level, while the IPCC AR5 projections are based on an assessment of multiple estimates for each contribution. As a result, the global mean SLR under the RCP8. 5 scenario is 71 cm in Slangen et al. [12], compared to 63 cm in IPCC AR5 (for 2081-2100 relative to 1986-2005). However, the methodology used for computing regional patterns is the same for the IPCC AR5 and the Slangen et al. [12] projections, and therefore any differences are mainly in the magnitude of the signal rather than the spatial distribution. 

For the contributions of the ice sheet dynamics, we used three different shapes for the uncertainty distributions, from three different sources. The IPCC estimate is a consensus estimate, based on a multitude of models and publications, the VW15 estimate is based on an expert elicitation, and the R15 estimate is based on physical and statistical models. We chose to use these published estimates of the actual ice sheet dynamical contributions rather than completely artificial uncertainty distributions, which probably would have not included a distribution with the shape of the R15 distribution. A systematic sensitivity study could be part of future work. Nevertheless, the current set of distributions shows that a highly skewed distribution leads to much higher estimates of the allowances. It is therefore very important to better understand the shape of the uncertainty distribution using physical understanding and modelling, as the uncertainties are critical to estimating future SLR and to decide on the necessary coastal protection. 

In order to focus on the influence of the shape of the uncertainty distributions rather than the actual magnitude, we shifted all distributions such that they matched the median values of IPCC. If we compare the differences in allowances before and after the shift of the medians, we find that for the VW15 scenario the shift in the median SLR projections is -12 cm, but the impact on the allowances (comparing allowances with and without median shift) is -8 cm. For the R15 scenarios, the medians were shifted by only -3 cm, and the impact on the allowances is -4 cm. 

With regard to the use of the Gumbel scale parameter, which is based on annual maximum sea levels, it is considered to be a method that is relatively straightforward, but perhaps not the ideal representation of sea-level extremes. It is relatively robust for temporal and spatial variations and therefore widely used for global applications [32], but, on the other hand, does not use all the information that is available. Some other potential methods, such as the Generalised Extreme Value distribution or a Generalised Pareto distribution, include more information, as they not only provide a slope, but also a shape of the extreme value distribution [33]. However, these more complicated distributions cannot (yet) be used to compute allowances. In a potential follow-up of this work, it would be useful to examine the effect of different types of extreme value distributions on sea-level allowances, but here we chose to focus on the sensitivity to the shape of the sea-level projection uncertainties alone. 

Generally, the longer the time period used, the better the estimate of the extreme sea-level variability [34]. Here, we use a minimum of 20 years of data in the tide gauge record as a constraint, but find that the average allowance only varies with a maximum of 0. 06 m for the VW15 scenario when selecting longer time series (which results in less tide gauge records): 2. 19 m for ≥20 years vs. 2. 13 for records ≥50 years. For IPCC and R15, the differences are even smaller and the selection of longer time series has a negligible effect on the allowances (respectively, 1. 01 and 1. 03 m for ≥20 years and, respectively, 1. 00 and 1. 02 m for ≥50 years). One caveat for all extreme value statistics is that they depend on the occurrence of extreme events in the observational period, which becomes a more relevant restriction for longer return periods (e. g., 1 in 10,000 years) [35]. As a result, the statistics may change significantly when new extreme events are included in the analysis and it is therefore important to keep them up to date. 

A final assumption regarding the statistics of the extremes is that they do not change over time, i. e., that both the magnitude with respect to the mean sea level and the frequency not change. This assumption is supported by observations and analysis of past extreme sea levels, which show that the observed changes in the frequency of flooding events is mainly caused by the increase in mean sea levels (e. g., [11] ). However, the statistics of storms, and therefore of sea-level extremes, could well change significantly in the future. It is important to keep in mind that sea level can also rise due to local land subsidence, for instance due to groundwater extraction. In some locations, the subsidence even exceeds the climate-driven sea-level rise and therefore amplifies the flooding frequency [36] and changes the location of the coastline due to submergence of the coast. However, information on vertical land motion (VLM, measured with GPS) is not available all over the world and projections are difficult since they directly depend on human decisions. Therefore, VLM and changing coastlines could not yet be taken into account in the sea-level projections.
section_name: Discussion

section_content: We have compared and contrasted three different scenarios for the projected dynamical ice sheet contribution to SLR. For each scenario, we combined information on sea-level extremes at more than 650 tide gauge stations with regional SLR projections in 2100 with respect to 2010, to compute so-called sea-level allowances. An allowance is a guide for the expected height a coastal structure needs to be elevated to keep the same flooding frequency. Allowances depend not only on the mean SLR, but also on the uncertainty in the mean sea-level projection and on the variability and magnitude of sea-level extremes. Each of the three sea-level projection scenarios has a different shape for the uncertainty distributions in the ice sheet dynamics contributions, leading to differences in the allowances. 

We found that the allowances change significantly for our skewed scenario, doubling their values compared to allowances based on a normal distribution of the ice sheet dynamics. This is due to a heavier upper tail associated with higher contributions from ice sheet dynamics (the low-probability/high-impact events). The IPCC-based 95% SLR at the tide gauge stations is 1. 08 m, with an average allowance of 1. 01 m, while the skewed VW15 scenario leads to a projected 95% SLR of 1. 50 m at the tide gauges and a much larger allowance of 2. 19 m. For the R15 scenario, which has a less skewed but still non-normal distribution, the 95% SLR at the tide gauges is 1. 12 m and the average allowance is 1. 03 m. This sensitivity to the skewness of the uncertainty distribution means that it is very important to improve estimates of ice sheet dynamics contributions to SLR, and in particular their uncertainties. As we focused on the shape of the uncertainty distribution here rather than the magnitude, computations of allowances should be repeated with improved ice sheet uncertainty distributions to deliver allowances that are suitable for decision making and coastal protection. 

The sea-level allowances are largest in regions with a relatively small variability in sea-level extremes and SLR projections of large magnitude and/or uncertainty. These are also the regions where frequencies in extremes will increase most if the allowances are not applied (i. e., coastal structure is not elevated), especially for the skewed probabilities in the VW15 scenario. Most of these vulnerable locations are in the tropics, where the largest regional sea-level changes (and uncertainties) are projected due to increasing thermal expansion and above-average regional contributions due to land ice mass loss [12]. This means that these regions would also be particularly vulnerable to larger contributions from the ice sheets. 

Here, the computation of allowances could only be done at tide gauge locations, as this is where high-resolution observational time series are available for reasonable lengths of time. It would, however, be interesting to try and fill in the gaps for the entire coastline, using modelling approaches. Recently, Muis et al. [32] published a global set of extreme sea levels and storm surges (the GTSR data set), including Gumbel scale parameters along the global coastline, which showed good agreement between observations and the model data set. This set could possibly be used in future work to provide sea-level allowances at all coastal locations, which would be helpful for coastal planning purposes.
section_name: Conclusions

section_name: Acknowledgments:

section_content: We gratefully acknowledge the contributors to the GESLA-2 database which are listed on www. gesla. org. A. B. A. S. and R. C. d. W. were funded by the NWO Netherlands Polar Programme, and T. J. R. was funded by the Netherlands Earth Science Center. P. L. W. was funded partly by the Leverhulme Trust.
section_name: Acknowledgments:

section_name: Conflicts of Interest:

section_content: Author Contributions: A. B. A. S. and R. S. W. v. d. W. conceived and designed the experiments; A. B. A. S. performed the experiments, analysed the data and wrote the paper; and all authors contributed to the writing of the paper. R. S. W. v. d. W., R. C. d. W. and T. J. R. designed the SEAWISE model and T. J. R. developed the code. J. R. H., P. L. W. and T. E. contributed with data and analysis tools.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest:

section_content: The following abbreviations are used in this manuscript:
section_name: Abbreviations

section_content: Fifth
section_name: AR5

section_content: SEAWISE is used to combine different uncertainty distributions, e. g., Greenland (P 1 ) and West-Antarctica (P 2 ), into one uncertainty distribution (e. g., P combined1 ). The composed distribution P combined1 (z) (where z is the projected SLR) exists for each z of all contributions of two independent distributions P 1 (z 1 ) and P 2 (z 2 ) for which the summed values z 1 and z 2 add to z. P combined1 (z) is the resulting combined uncertainty distribution of P 1 (z 1 ) and P 2 (z 2 ):

All data within 6σ-boundaries around the mode of each distribution is selected, which usually covers over 99. 9% of the data. The distribution is normalized based on this selected part. In Equation (A1), the interval counter m runs in steps of ∆z between c and d corresponding with z c and z d, respectively the left and right 6σ-boundaries of P 2, while taking z 1 = z -m∆z and z 2 = m∆z. To obtain the entire distribution P combined1 (z), Equation (A1) is calculated for all z between z a at the left 6σ-boundary of P 1 and z b + z d, where z b is the right 6σ-boundary of P 1 and z d the right 6σ-boundary of P 2. In this way, the uncertainty distribution for each SLR contribution can be added, until all uncertainty distributions are combined (Figure A1 ). In addition to providing the combined uncertainty distributions for single points, SEAWISE can produce global maps of the SLR mode, SLR median and a number of SLR nth percentiles or SLR nth σ-intervals [18].
section_name: Appendix A. SEAWISE Methodology