Variability of Sea-Surface Height (SSH) from ocean dynamic processes is an important component of sea-level change. In this study we dynamically downscale a present-day control simulation of a climate model to replicate sea-level variability in the Northwest European shelf seas. The simulation can reproduce many characteristics of sea-level variability exhibited in tide gauge and satellite altimeter observations. We examine the roles of lateral ocean boundary conditions and surface atmospheric forcings in determining the sea-level variability in the model interior using sensitivity experiments. Variability in the oceanic boundary conditions leads to uniform sea-level variations across the shelf. Atmospheric variability leads to spatial SSH variability with a greater mean amplitude. We separate the SSH variability into a uniform loading term (change in shelf volume with no change in distribution), and a spatial redistribution term (with no volume change). The shelf loading variance accounted for 80% of the shelf mean total variance, but this drops to ~ 60% around Scotland and in the southeast North Sea. We analyse our modelled variability to provide a useful context to coastal planners and managers. Our 200-year simulation allows the distribution of the unforced trends (over 4–21 year) of sea-level changes to be quantified. We found that the 95th percentile change over a 4-year period can lead to coastal sea-level changes of ~ 58 mm, which must be considered when using smooth sea level projections. We also found that simulated coastal SSH variations have long correlation length-scales, suggesting that observations of interannual sea-level variability from tide gauges are typically representative of > 200 km of the adjacent coast. This helps guide the use of tide gauge variability estimates.

Sea level on the northwestern European shelf (NWES) varies substantially from year to year. Removing explained parts of interannual sea level variability from observations helps to improve estimates of long-term sea level trends. To this end, the contributions of different drivers to interannual sea level variability need to be understood and quantified. We quantified these contributions for the entire NWES by performing sensitivity experiments with a high-resolution configuration of the Regional Ocean Modeling System (ROMS). The lateral and atmospheric boundary conditions were derived from reanalyses. We compared our model results with satellite altimetry data and used our sensitivity experiments to show that nonlinear feedbacks cause only minor interannual sea level variability on the shelf. This indicates that our experiments can be used to separate the effects of different drivers. We find that wind dominates the variability of annual mean sea level in the southern and eastern North Sea (up to 4.7-cm standard deviation), whereas the inverse barometer effect dominates elsewhere on the NWES (up to 1.7-cm standard deviation). In contrast, forcing at the lateral ocean boundaries results in small and coherent variability on the shelf (0.5-cm standard deviation). Variability driven by buoyancy fluxes ranges from 0.5- to 1.3-cm standard deviation. The results of our sensitivity experiments explain the (anti)correlation between interannual sea level variability at different locations on the NWES and can be used to estimate sea level rise from observations in this region with higher accuracy.

Changes in ocean properties and circulation lead to a spatially non-uniform pattern of ocean dynamic sea-level change (DSLC). The projections of ocean dynamic sea level presented in the IPCC AR5 were constructed with global climate models (GCMs) from the Coupled Model Intercomparison Project 5 (CMIP5). Since CMIP5 GCMs have a relatively coarse resolution and exclude tides and surges it is unclear whether they are suitable for providing DSLC projections in shallow coastal regions such as the Northwestern European Shelf (NWES). One approach to addressing these shortcomings is dynamical downscaling – i.e. using a high-resolution regional model forced with output from GCMs. Here we use the regional shelf seas model AMM7 to show that, depending on the driving CMIP5 GCM, dynamical downscaling can have a large impact on DSLC simulations in the NWES region. For a business-as-usual greenhouse gas concentration scenario, we find that downscaled simulations of twenty-first century DSLC can be up to 15.5 cm smaller than DSLC in the GCM simulations along the North Sea coastline owing to unresolved processes in the GCM. Furthermore, dynamical downscaling affects the simulated time of emergence of sea-level change (SLC) above sea-level variability, and can result in differences in the projected change of the amplitude of the seasonal cycle of sea level of over 0.3 mm/yr. We find that the difference between GCM and downscaled results is of similar magnitude to the uncertainty of CMIP5 ensembles used for previous DSLC projections. Our results support a role for dynamical downscaling in future regional sea-level projections to aid coastal decision makers.