The impact of meltwater discharge from the Greenland ice sheet on the Atlantic nutrient supply to the northwest European shelf
- 1Ocean Department, Max-Planck-Institute for Meteorology, Bundesstr. 53, 20146 Hamburg, Germany
- 2Institute of Coastal Research, Helmholtz-Zentrum Geesthacht, Max-Planck-Str. 1, 21502 Geesthacht, Germany
Correspondence: Moritz Mathis (email@example.com)
Projected future shoaling of the wintertime mixed layer in the northeast (NE) Atlantic has been shown to induce a regime shift in the main nutrient supply pathway from the Atlantic to the northwest European shelf (NWES) near the end of the 21st century. While reduced winter convection leads to a substantial decrease in the vertical nutrient supply and biological productivity in the open ocean, vertical mixing processes at the shelf break maintain a connection to the subpycnocline nutrient pool and thus productivity on the shelf. Here, we investigate how meltwater discharge from the Greenland ice sheet (GIS), not yet taken into account, impacts the mixed layer shoaling and the regime shift in terms of spatial distribution and temporal variability. To this end, we have downscaled sensitivity experiments by a global Earth system model for various GIS melting rates with a regionally coupled ocean–atmosphere climate system model. The model results indicate that increasing GIS meltwater discharge leads to a general intensification of the regime shift. Atlantic subpycnocline water masses mixed up at the shelf break become richer in nutrients and thus further limit the projected nutrient decline on the shelf. Moreover, the stronger vertical nutrient gradient through the pycnocline results in an enhanced interannual variability of on-shelf nutrient fluxes which, however, do not significantly increase variations in nutrient concentrations and primary production on the shelf. Due to the impact of the GIS meltwater discharge on the NE Atlantic mixed layer depth, the regime shift becomes initiated earlier in the century. The effect on the onset timing, though, is found to be strongly damped by the weakening of the Atlantic meridional overturning circulation. A GIS melting rate that is even 10 times higher than expected for emission scenario Representative Concentration Pathway (RCP) 8.5 would not lead to an onset of the regime shift until the 2070s.
Previous climate change impact studies identified a general weakening of the biological productivity of the outer northwest European shelf (NWES) as a regional response to a globally warming climate (e.g., Holt et al., 2012, 2016; Gröger et al., 2013; Wakelin et al., 2015; Schrum et al., 2016; Mathis et al., 2018). The main driver has been attributed to a reduction in the oceanic nutrient import from the adjacent northeast (NE) Atlantic. Under present-day conditions, up to 80 %–90 % of the nutrient inventory on the shelf is advected from the open ocean via cross-shelf break transport of nutrient-rich Atlantic water masses (e.g., Vermaat et al., 2008; Liu et al., 2010; Thomas et al., 2010). The nutrient concentrations of the water masses flushing the shelf are primarily controlled by the maximum depth of the wintertime mixed layer (MLD) west of the shelf break (e.g., Holliday, 2003; Williams et al., 2011; Holt et al., 2012; Gröger et al., 2013; Mathis et al., 2018). The warming and freshening of the upper North Atlantic projected by global climate models induces a weakening of the buoyancy-driven convection and thus a reduction of the wintertime MLD and upper-ocean nutrient concentration (e.g., Yool et al., 2015; Fu et al., 2016; Alexander et al., 2018). Accordingly, previous downscaling simulations mainly based on the Special Report on Emissions Scenarios (SRES) emission scenario A1B indicated that a nutrient decline in the upper NE Atlantic would lead to a similar nutrient decline on the NWES, limiting net primary production and atmospheric CO2 uptake in open shelf areas (e.g., Holt et al., 2012; Gröger et al., 2013; Wakelin et al., 2015; Schrum et al., 2016). The parent global models used to provide the forcing data, however, did not include interactive ice sheet models and hence miss an important factor in their projections of future changes in the North Atlantic circulation and water mass characteristics: the melting of the Greenland ice sheet (GIS). In particular, the effect of increasing meltwater discharge to the subpolar North Atlantic on the ocean mixed layer shoaling is not taken into account and therefore is not included in the downscaling simulations for the NWES.
Observational data have provided evidence for a significant increase in GIS mass loss during the recent decades (Luthcke et al., 2006; Vaughan et al., 2013; van den Broeke et al., 2017; Chen, 2019). Modeling studies considered by the Fifth Assessment Report of the United Nations Intergovernmental Panel on Climate Change (IPCC AR5) agree that the GIS will continue to decrease in area and volume in a warmer climate as a consequence of increased melting rates not compensated by increased snowfall and amplified by positive feedbacks (Mikolajewicz et al., 2007; Collins et al., 2013; Church et al., 2013; Vizcaino et al., 2015). Enhanced GIS meltwater discharge has been shown to decrease wintertime MLDs in the North Atlantic due to the haline-induced reduction of surface buoyancy (Böning et al., 2016; Oliver et al., 2018). To that effect, numerous other studies suggest that substantial meltwater input to the North Atlantic will weaken the Atlantic meridional overturning circulation (AMOC) or even lead to a collapse (Schiller et al., 1997; Stouffer et al., 2006; Weijer et al., 2012; Sévellec et al., 2017; Liu et al., 2018). The resulting decrease in northward heat transport might further affect the surface buoyancy and thus the permanent stratification and MLD (Hu et al., 2009; Woollings et al., 2012; Wouters et al., 2012; Liu et al., 2017).
In a recent study by Mathis et al. (2019), it has been shown that for a critical ML shoaling according to emission scenario Representative Concentration Pathway (RCP) 8.5, upper North Atlantic water masses lose their dominant influence on shaping biogeochemical conditions on the NWES. This is due to a regime shift in the main supply pathway of Atlantic nutrients to the shelf. In the open ocean, the ML shoaling leads to a reduction of the vertical nutrient supply from intermediate depths to the euphotic zone. At the shelf break, however, various mixing processes maintain a connection to the Atlantic subpycnocline nutrient pool, e.g., due to the interaction of tidal currents with the abrupt change in the topography, internal waves and instabilities of the slope current. Nutrient-enriched water masses mixed up near the shelf break spread to the NWES and lead to a weaker nutrient decline in open shelf areas, compared to the adjacent upper NE Atlantic. As a consequence, the projected weakening of the biological productivity on the shelf is not as strong as would have resulted from a sole influence of upper NE Atlantic water masses. Finally, nutrient concentrations in the upper NE Atlantic reach lower levels than on the shelf, leading to a reversal of the ocean–shelf nutrient gradient and to the development of a nutrient front along the continental margin.
Furthermore, during the shallow-ML regime, the nutrient transport to the NWES is subject to a rapid increase in interannual to multidecadal variability. As near the end of the 21st century the ML in the NE Atlantic becomes as shallow as the shelf edge, upper-ocean conditions are more sensitive to variations in the atmospheric forcing. In particular, a positive MLD anomaly (deepening) leads to an erosion of warm and saline subpycnocline water masses and initiates a positive feedback on the surface heat flux, the upper-ocean buoyancy and the MLD. The resulting enhanced variability of the on-shelf nutrient transport affects the variability in pre-bloom nutrient concentrations and annual primary production on the NWES.
Similar to previous downscaling studies, though, the freshwater cycle in the underlying simulations by Mathis et al. (2019) did not account for a future increase in GIS meltwater discharge, ignoring its impact on the density distribution in the upper North Atlantic. Hence, it is not yet clear how and to what extent GIS meltwater discharge impacts the proposed regime shift in the Atlantic nutrient supply to the NWES. In particular, it has been shown by Mathis et al. (2019) that during the shallow-ML regime the transfer of anomalies in the upper NE Atlantic circulation and water mass characteristics to the NWES becomes undermined by the associated shift in the dominant ocean–shelf exchange processes from horizontal advection to vertical mixing. In the present study, our general goal is to explore the influence of GIS meltwater discharge on the signal transfer to the NWES and thus to advance our understanding of the coupling mechanisms in the transition zone between the open ocean and the shallow shelf under a strong radiative forcing scenario. Specifically, we focus on the following questions: Does the regime shift become initiated earlier in the 21st century due to the impact on the MLD? Does the ML become shallower than the depth of the shelf edge, providing a direct connection between the shelf and the Atlantic subpycnocline nutrient pool? Does the enhanced variability of the nutrient concentrations at the shelf break become even stronger due to the impact on the permanent stratification? If such changes are found, what are the consequences for the nutrient budget and primary production in open shelf areas such as the northern North Sea, which are predominantly influenced by Atlantic inflow? In the present study, we address these questions by an investigation of the impact of GIS meltwater discharge on both the mean Atlantic nutrient supply to the NWES and its temporal variability on interannual to multidecadal scales.
The original study by Mathis et al. (2019), henceforth referred to as M19, was based on a dynamical downscaling of global climate projections for emission scenario RCP8.5 by the Max-Planck-Institute Earth System Model (MPI-ESM-LR). It was shown that the global model was not able to capture the drastic changes in the on-shelf nutrient transport near the end of the 21st century. By contrast, pre-bloom nutrient concentrations on the NWES simply reflected the continuous upper-ocean nutrient decline of the NE Atlantic due to an underrepresentation of local processes relevant to adequately simulating a comparatively small region such as the NWES. Therefore, in the present study, we utilize the same high-resolution regionally coupled climate system model as in M19 and conducted a set of downscaling simulations for emission scenario RCP8.5 with additional freshwater runoff along the coast of Greenland. To account for uncertainties due to natural variability, we performed ensembles of three realizations for experiments with and without GIS meltwater discharge as well as a pre-industrial control run. The downscaling model system is briefly described in Sect. 2.1, and the experiment design is introduced in Sect. 2.2. Results are presented and discussed in Sects. 3 and 4, respectively. Concluding remarks are given in Sect. 5.
2.1 Model description
The downscaling simulations presented here are performed with a regionally coupled ocean–atmosphere climate system model, consisting of the global Max-Planck-Institute Ocean Model (MPIOM) with a zoom on the NWES, the Hamburg Ocean Carbon Cycle Model (HAMOCC), the Regional Atmosphere Model (REMO) and the global Hydrological Discharge Model (HD). A schematic of the ocean grid and coupling domain is shown, e.g., in M19. Simulated physical and biogeochemical conditions on the NWES have been evaluated in various studies employing MPIOM/HAMOCC in coupled and uncoupled modes (Gröger et al., 2013; Mathis et al., 2015, 2018, 2019; Pätsch et al., 2017; Hátún et al., 2017).
MPIOM (Maier-Reimer, 1997; Marsland et al., 2003; Jungclaus et al., 2013) is the ocean–sea ice component of the global Max-Planck-Institute Earth System Model (MPI-ESM). The primitive equations of oceanic motion are discretized on an Arakawa C grid with z coordinates and free surface. Hydrostatic and Boussinesq approximations are applied. A second-order total variation diminishing scheme (Sweby, 1984) is used to simulate tracer and momentum advection. Vertical mixing is implemented after Pacanowski and Philander (1981), with an additional parameterization for wind-induced stirring (Marsland et al., 2003). The full lunisolar tidal potential is calculated according to Thomas et al. (2001).
Using a global ocean component in the downscaling model avoids problematic influences of open lateral boundary conditions necessary to run conventional regional ocean models (Mathis et al., 2018). Higher grid resolutions in the subpolar North Atlantic and NWES regions are achieved by using a stretched grid configuration with non-diametrical poles located over central Europe and North America (Mathis et al., 2018, 2019). In the present study, MPIOM has been run at a nominal horizontal grid resolution of 0.6∘, yielding mesh sizes of about 5 to 12 km on the NWES. The vertical is resolved by 30 levels with eight levels in the upper 100 m.
Biogeochemical processes in the ocean are simulated by HAMOCC (Maier-Reimer et al., 2005; Ilyina et al., 2013). Marine biology dynamics is represented by nutrients, phytoplankton, zooplankton, detritus and dissolved organic matter (Six and Maier-Reimer, 1996). Biogeochemical cycles and trophic levels are connected by nutrient uptake and remineralization of organic matter (Six and Maier-Reimer, 1996; Kloster et al., 2006). A sediment module accounts for deposition and dissolution processes as well as the biogeochemistry of pore water in the upper bioturbated layers (Heinze et al., 1999). Sediment resuspension is based on the concept of incipient motion (Wilcock, 2004; Mathis et al., 2019). Tracer transport is simulated by MPIOM.
Over the EURO-CORDEX domain (e.g., Kotlarski et al., 2014), MPIOM is interactively coupled with REMO (Jacob and Podzun, 1997; Jacob et al., 2001) at a coupling time step of 1 h. REMO has been run at a uniform horizontal resolution of about 25 km and 27 hybrid levels. Sea surface fluxes of heat, momentum and freshwater are transferred to MPIOM by the OASIS3 coupler (Valcke, 2013). In turn, sea surface temperature and sea ice cover and thickness are passed from MPIOM to REMO. At the land–atmosphere interface, a simple bucket scheme is applied for the treatment of soil hydrology. Freshwater drainage and runoff are calculated by the HD model (Hagemann and Dümenil-Gates, 2001).
As forcing data for the downscaling model, we use 6-hourly atmospheric output from the global Earth system model MPI-ESM in its low-resolution version (LR; Giorgetta et al., 2013). The data are used as lateral boundary conditions for REMO and as input for the calculation of atmosphere–ocean fluxes by MPIOM outside the domain covered by REMO.
2.2 Experiment design
The future regime shift in Atlantic nutrient supply to the NWES (M19) has been identified for emission scenario RCP8.5. To investigate the impact of GIS meltwater discharge in the present study, we therefore conduct RCP8.5 simulations with various amounts of additional freshwater runoff along the coast of Greenland (Fig. 1 and Table 1). The spatial distribution of the runoff (Fig. 2) follows the climatology by Bamber et al. (2012), based on satellite observations and regional climate modeling. The seasonal cycle (Fig. A1) has been derived by Martin et al. (2019). The annual mean GIS freshwater flux according to this study corresponds to 0.05 Sv and is comparable to the estimate of about 0.04 Sv starting from the year 2010 by Yang et al. (2016) and Bamber et al. (2018). In our simulations, the prescribed freshwater fluxes enter the surface layer of the ocean model, thus ignoring that many marine-terminating outlet glaciers have a grounding line depth several hundred meters below sea level (An et al., 2017; Morlighem et al., 2017).
All experiments are performed with both the global MPI-ESM-LR model to create forcing data and the MPIOM/HAMOCC/REMO downscaling model. The simulations cover the period 1920–2150 in consistency with M19. The MPI-ESM-LR simulations are initialized by 1920 conditions of the MPI-ESM-LR consortium simulations which contributed to CMIP5 (Coupled Model Intercomparison Project, phase 5). The original CMIP5 simulations could not be used here because they were run on a former high-performance computer (HPC) and hence are inconsistent with our GIS discharge experiments. A bitwise reproduction of a simulation is not possible on different HPCs, leading to independent realizations even when started from identical initial conditions. To minimize the effect of long-term model drift, the MPIOM/HAMOCC/REMO simulations are initialized by the end of the pre-industrial control simulation (year 2150) used in M19.
Experiment C0 is the control run at constant pre-industrial atmospheric greenhouse gas concentrations without additional GIS meltwater discharge. This experiment allows us to account for remaining model drift as well as to distinguish anthropogenic climate change signals from natural variability.
Experiment E0 is a combination of the historical period 1920–2005 and the RCP8.5 scenario period 2006–2150. Also, here, no additional GIS meltwater discharge is taken into account, and hence this experiment is comparable to the original simulation by M19. However, in the present study, we have performed an ensemble of three realizations (Table 1) differing in the initial conditions used for the parent MPI-ESM-LR simulations. This allows us to assess the influence of internal variability of the forcing global model.
In experiment E010, a linear increase of GIS meltwater discharge from 0 Sv in the year 2006 to 0.1 Sv in 2100 is prescribed and held constant thereafter. The melting rate of 0.1 Sv is adopted from simulations by interactive ice sheet models for a quadrupling of pre-industrial atmospheric CO2 concentrations (e.g., Mikolajewicz et al., 2007; Vizcaíno et al., 2010) and used in several other studies to investigate GIS melting impacts, e.g., on the strength of the AMOC (Stouffer et al., 2006; Yu et al., 2016; Swingedouw et al., 2015). The assumption of a linear increase is an idealized approach to deal with the uncertainty in the construction of a hydrological sensitivity parameter, often defined as a constant freshwater discharge per degree atmospheric warming (e.g., Zickfeld et al., 2008; Kuhlbrodt et al., 2009), and has likewise been applied, e.g., in Jungclaus et al. (2006). For this experiment, we have also performed an ensemble of three realizations, branching from ensemble E0 in the year 2006. Thus, E0 and E010 represent the core experiments to investigate the impact of GIS meltwater discharge on the regime shift in the nutrient supply to the NWES.
In experiment E025, we increase the GIS meltwater discharge to 0.25 Sv during the same period as in E010 (2006–2100). This experiment is designed to further stress results from experiment E010 and to gain more insight into the involved physical processes and driving mechanisms. Besides, the higher discharge value accounts for uncertainties in GIS melting rates projected by ice sheet models (Church et al., 2013; Agarwal et al., 2015). The range covered by experiments E010 and E025 is similar to that implicated in past AMOC shutdowns in climate history (e.g., Hemming, 2004; Jongma et al., 2013; Ziemen et al., 2019). E025 can therefore be considered as a high-end GIS melting estimate and sensitivity experiment. A similarly high discharge rate of 0.2 Sv has been considered as upper limit, e.g., by Kwiatkowski et al. (2019).
It shall be shown in Sect. 3.3 that the impact of GIS meltwater discharge on the onset of the regime shift is comparatively weak. In experiment E100, we therefore provoke a more distinct signal by further increasing the meltwater discharge to an extreme value of 1.0 Sv in the year 2100. This high discharge rate is purely motivated by process understanding and exceeds any present estimate of GIS runoff during the 21st century. In fact, given a present-day GIS volume of about 2.9 × 1015 m3, it would lead to a complete disintegration of the ice sheet in the first half of the 22nd century, depending on the surface mass balance.
2.3 Model evaluation
Numerous other studies have investigated the impact of an additional freshwater hosing to the northern North Atlantic on the ocean circulation and biological productivity in the context of potential future climate change. To evaluate our downscaling simulations, we briefly compare here projected change signals most relevant for our study, i.e., the shoaling of the wintertime MLD, the weakening of the AMOC and the decline of upper-ocean nutrient concentrations, as well as some variables of general interest, such as sea surface temperature (SST) and net primary production.
Maximum (March) MLD in the NE Atlantic simulated by our downscaling model is about 800 m for 1971–2000 (Fig. 3a) and is generally overestimated by about 200 m (also in the main convection sites) compared to the observational climatology by de Boyer Montégut et al. (2004). An overestimation of wintertime MLD and deep water formation has also been found as a common feature in the majority of global Earth system models of CMIP5 generation (Heuzé, 2017). Partly, this is due to the often applied density difference criterion, which leads to comparatively deep MLDs (e.g., Costoya et al., 2014; Courtois et al., 2017). In the present study, we have used a density threshold of Δσθ = 0.125 kg m−3 (e.g., Suga et al., 2004).
During the 21st century, NE Atlantic MLD decreases by about 400 m under RCP8.5 conditions without additional GIS meltwater release (E0; Fig. 3b) and by about 500 m with an increasing GIS melting rate of up to 0.1 Sv (E010; Fig. 3c). The MLD changes projected by other studies reflect the wide range of simulated ocean stability across the applied circulation models. In a stronger warming scenario by Kuhlbrodt et al. (2009), North Atlantic MLD decreases but deep water formation does not stop completely (in the Nordic Seas). With a maximum hosing of 0.5 Sv, there is no MLD deeper than 200 m anywhere north of 45∘ N. These results are similar to our simulations with both significantly weaker radiative forcing and freshwater input.
The present-day strength of the AMOC has been estimated from observations to be around 16–20 Sv (Talley et al., 2003; Kanzow et al., 2010; Fu et al., 2018) and is well captured by our downscaling model, simulating 17.7 ±1.1 Sv at 30∘ N during 1971–2000. Associated with the range of MLD sensitivity, Earth system models show a large range of AMOC sensitivity to additional freshwater input to the North Atlantic. Accordingly, previous studies can be grouped by their simulated AMOC responses into low freshwater impact (< 30 % AMOC weakening within 100 years due to ∼0.1 Sv freshwater hosing), moderate impact (30 %–60 % AMOC weakening) and high impact (> 60 % AMOC weakening). Low impacts have been projected, e.g., by Mikolajewicz et al. (2007); Driesschaert et al. (2007) and Hu et al. (2009), low to moderate impacts (depending on the model), e.g., by Stouffer et al. (2006) and Swingedouw et al. (2013), moderate impacts, e.g., by Jungclaus et al. (2006) and Swingedouw et al. (2015), and high impacts, e.g., by Zickfeld et al. (2008) and Kuhlbrodt et al. (2009). The additional AMOC weakening of about 10 % in our hosing experiment (E010; see Sect. 3.1) classifies our downscaling model as a low impact model with respect to AMOC sensitivity.
Furthermore, in our simulations, the hosing causes a freshening of the entire subpolar and arctic North Atlantic compared to E0. The freshening of the subpolar region is also consistent among global Earth system models (Swingedouw et al., 2013). Some models, however, show a positive salinity anomaly in the Nordic Seas because more Atlantic waters reach the Arctic. This is also associated with a warming of the Nordic Seas but a cooling everywhere else, which is also indicated in our simulations. With additional freshwater input, the cooling becomes stronger and the warming of the Nordic Seas weaker. Experiments by Kuhlbrodt et al. (2009), by contrast, show amplified cooling in particular in the Nordic Seas.
Similarly, the impact of climate change on marine primary production and upper-ocean nutrient concentrations has considerable spread among global Earth system models (Bopp et al., 2013; Fu et al., 2016). Nevertheless, in the subpolar North Atlantic, the majority of models show a negative response under emission scenario RCP8.5 due to a weakening of the vertical nutrient supply from intermediate depths to the upper ocean. In agreement with these trends, the biogeochemical response in our simulations to a MLD decrease and AMOC weakening is a general reduction in surface nutrient concentrations and net primary production in the NE Atlantic by about 60 % and 40 %, respectively (see Sect. 3.1). Comparable changes have been found by Zickfeld et al. (2008) and Kwiatkowski et al. (2019). A rather strong signal of up to 90 % reduction in primary production was simulated by Kuhlbrodt et al. (2009) due to their stronger warming scenario and higher freshwater input. Nevertheless, similar to our study, the climate-induced changes showed a distinctly stronger impact on primary production than the relative changes induced by freshwater hosing.
It has to be noted, however, that the cited model experiments differ in the precise location, timing and amount of the imposed freshwater hosing. Unless an interactive ice sheet model is involved, the additional freshwater is typically released either uniformly over a wide region in the North Atlantic or restricted to the coast of Greenland with a sudden introduction of a constant discharge rate or a transient increase over several decades or centuries under past, present-day or future climatic conditions. The differences in the experiment design impede direct comparison, often requiring scaling assumptions, and may contribute to the wide spread in the simulated responses of the ocean circulation. After all, with respect to MLD and AMOC perturbations, our downscaling model is comparatively sensitive to changes in the atmospheric forcing and comparatively insensitive to additional freshwater input to the northern North Atlantic.
In M19, the transition to the future shallow-ML regime was identified to happen near the end of the 21st century. In the present study, climate change signals are therefore referred to as differences between the distant future period (2101–2150) and the recent past (1971–2000). Time series are shown for 1960–2150, excluding the period 1920–1959 considered as spinup. For our analysis, oceanic nutrients in general are represented by dissolved phosphate concentrations, consistent with M19. Nitrate concentrations, for example, are subject to nitrogen fixation and denitrification processes, making phosphate a more appropriate candidate for the investigation of the impact of climatic changes on oceanic tracer advection.
The bathymetry of the NWES region is shown in Fig. 4. Transect S1 is used to illustrate changes in the vertical nutrient distribution near the shelf break. Budgets for the northern North Sea (Sect. 3.4) are calculated over an area bounded in the north by the solid line (Orkney–Shetland and Shetland–Norway at 60.4∘ N) and in the south by the 50 m isobath.
3.1 Changes in the mean ocean–shelf nutrient transport
The general freshening and warming of the upper ocean in the subpolar North Atlantic during emission scenario RCP8.5 leads to a decrease in sea surface salinity, a strengthening of the permanent stratification (Fig. A2a, b) and a reduction of the wintertime MLD in the NE Atlantic (Fig. 5a, b). The additional freshwater hosing due to GIS meltwater discharge further intensifies these changes (Figs. 5c, d and A2c, d).
In experiment E0, the shoaling of the wintertime ML results in a limited vertical nutrient supply to the euphotic zone (Fig. 5) and thus weakens the biological productivity during the following summer season (Gröger et al., 2013; Mathis et al., 2018, 2019). As near the end of the 21st century the ML becomes as shallow as the NWES, mixing processes at the shelf break maintain local vertical nutrient fluxes. Nutrient-enriched water masses spread from the shelf break to the interior of the shelf and yield weaker reductions in nutrient concentrations and productivity there (M19). When GIS meltwater is taken into account (experiments E010 and E025), the MLD in the NE Atlantic decreases further and nutrient concentrations in the upper ocean decline more (Fig. 5c, d). Surprisingly though, nutrient concentrations on the shelf decrease even less than in the experiment without GIS meltwater discharge (E0). The regime shift in E0 is characterized by a disruption of the ocean–shelf signal transfer by mixing processes at the shelf break. The contrasting impact on the nutrient distribution seen in the meltwater discharge experiments (E010 and E025) indicates that also here the change signal in the NE Atlantic is not coherently transferred to the shelf but there are other mechanisms involved influencing the on-shelf nutrient transport.
These opposite effects in the open ocean and shallow shelf are a consequence of the effect of the GIS meltwater on the North Atlantic circulation at a basin-wide scale. In E0, the climate change signal from the atmosphere (surface warming and increased moisture transport) drives a general slowdown of the AMOC (Fig. 6a). The resulting weakening of the northward heat transport enhances the meridional temperature gradient in the ocean and atmosphere and induces a northward shift of the atmospheric wind field, as indicated in most CMIP3 and CMIP5 models (e.g., Yin, 2005; Woollings et al., 2012; Hu et al., 2013; Fischer et al., 2017). The associated strengthening of the westerlies at midlatitudes by up to 30 % and the increase in wind stress curl over the NE Atlantic by up to 15 % drive an expansion of the subpolar gyre (SPG) towards the east (Fig. 7). Moreover, in the western part of the gyre system, the boundary between the SPG and subtropical gyre moves northward.
In the upper ocean, the widening of the SPG leads to a reduction in the northward intrusion of saline subtropical water to the NE Atlantic and to a stronger eastward propagation of low-salinity subpolar water. The related decrease in the density strengthens the stratification and reduces the MLD in addition to the climate change signal from the atmosphere. Accordingly, GIS meltwater discharge causes a stronger salinity drop in the SPG region and hence amplifies the impact on the stratification and MLD (Fig. 5c, d). The resulting weakening of the vertical nutrient supply from the ocean's interior leads to lower nutrient concentrations in the NE Atlantic wintertime ML when more GIS meltwater is added.
The additional source for the nutrients on the NWES (Fig. 5c, d), by contrast, is found in Atlantic water masses below the pycnocline. The weakening of the AMOC (Fig. 6a) is associated with a general weakening of the North Atlantic deep circulation. Deep water formation in the Labrador Sea, Irminger Sea and Nordic Seas breaks down around 2080 in all experiments. Subpycnocline water masses thus become older with respect to the time since they were last in contact with the atmosphere. Continuous remineralization of sinking particulate organic matter is at the expense of dissolved oxygen (e.g., Keeling et al., 2010) but releases nutrients (Fig. 6b, c). When GIS meltwater is added, the density of the upper ocean decreases more, leading to a slightly stronger slowdown of the AMOC (Fig. 6a). Thus, more organic material is remineralized, enhancing subpycnocline oxygen utilization and increasing nutrient concentrations. During the future shallow-ML regime, finally, the nutrient-enriched subpycnocline water masses are mixed to the upper water column near the shelf break and spread over the shelf, causing elevated nutrient concentrations there.
Ocean–shelf exchange is mainly governed by wind-driven Ekman transport, internal tidal waves and cross-shelf break transport induced by the northward-flowing slope current (e.g., Huthnance et al., 2009; Simpson and Sharples, 2012; Ruiz-Castillo et al., 2019). The simulated vertical profile of cross-shelf break transport during winter (DJFM) integrated along the 200 m isobath between 47 and 61.5∘ N (Fig. 8a) shows the typical structure of a strong wind-driven on-shelf surface flow (0–40 m), an on-shelf flow in the interior of the water column (40–130 m) and an off-shelf Ekman drain at the bottom (130–200 m) due to basal stress of the slope current (Graham et al., 2018). Independent of the GIS melting rates, the net on-shelf winter transport through this section weakens under RCP8.5 by about 23 % (from 1.11 in 1971–2000 to 0.86 Sv in 2101–2150). Most prominent changes are indicated in the on-shelf interior and off-shelf bottom transports (Fig. 8a), induced by a weakening of the slope current (Fig. 9a).
The corresponding phosphate (PO4) profiles at the shelf edge (Fig. 8b) illustrate the projected general nutrient decline in the upper water column as well as the influence of the increasing nutrient concentrations at deeper levels. When GIS meltwater is added, the nutrient enrichment of subpycnocline Atlantic water masses results in a relative increase of future nutrient concentrations near the shelf break, mitigating the projected nutrient decline on the shelf. Net on-shelf PO4 fluxes decrease from 839 (1971–2000) to 311 mol s−1 (2101–2150) in E0, 265 in E010 and 211 mol s−1 in E025. This increasing reduction among the experiments, however, is dominated by the decreasing volume transports (−0.25 Sv in E0, −0.30 in E010, −0.32 in E025) and does not reflect the changes in the nutrient concentrations. As future mean pre-bloom PO4 concentrations on the shelf range between 0.55 and 0.7 mmol m−3 in E0, with maximum values of 0.8 mmol m−3 in E025 (Fig. 5b–d), the on-shelf transport in the interior of the water column is an essential feature to establish these elevated concentrations. It is therefore important to understand the weakening of the slope current and its impact on the cross-shelf break exchange.
Primarily, the northward-flowing slope current along the continental margin of the NWES is driven by the meridional density gradient at intermediate depths to the west of the shelf break (Marsh et al., 2017). The region of influence has been identified to span from about 45 to 62∘ N and zonally as far as 28∘ W. The meridional density gradient in this region induces a geostrophic flow towards the east which turns northward as it approaches the shelf break, forming the core of the slope current. Under RCP8.5 conditions (E0), the meridional density difference (between north and south) decreases from about 0.51 (1971–2000) to 0.35 kg m−3 (2101–2150; Fig. 9b and Table 2). The slowdown of the AMOC and the expansion of the SPG entail an eastward migration of old Labrador Sea water and water of Arctic and Nordic Seas origin flowing through the Denmark Strait (Fig. A3) as well as a reduction of the northward migration of warm and saline subtropical water, leading to a freshening and cooling of NE Atlantic intermediate waters (Fig. 10). A stronger impact on the density though is indicated in the northern region of influence (−0.35 kg m−3; Table 2), where the cooling effect is weakest and the water temperature thus still increases due to the increasing radiative forcing. In the south, the remaining influence of warm subtropical water and Mediterranean overflow water induces a more moderate change (−0.19 kg m−3). Because of the stronger density drop in the north than in the south, the meridional density gradient decreases.
As GIS meltwater is added, both the freshening and cooling of the NE Atlantic intermediate water become stronger but their effects on the water density roughly cancel (Table 2). The decrease of the meridional density gradient is therefore very similar for all experiments with and without meltwater discharge. Moreover, most significant changes happen during the period when the SPG widens (Fig. 11). Around the year 2080, the SPG reaches its maximum size, limited by the topographic boundary of the European continent (Fig. 7). After 2080, the meridional density gradient as well as the core of the slope current do not change further (Fig. 9), resulting in similar subsurface transport profiles among the experiments during the new shallow-ML regime (Fig. 8a).
Regarding the nutrient fluxes to the shelf, though, the weakening of the interior on-shelf transport from historical to future conditions is partly compensated by the influence of nutrient-rich subpycnocline water masses mixed up at the shelf break. When GIS meltwater is taken into account, subpycnocline nutrient concentrations become even higher and thus lead to an intensification of the developed ocean–shelf nutrient front (Fig. 5c, d).
3.2 Changes in the variability of ocean–shelf nutrient transport
The interannual to multidecadal variability of the advective nutrient supply to the NWES has been shown to increase rapidly during the shallow-ML regime (M19). The main mechanism behind this is enhanced MLD variations near the shelf break which control the nutrient concentration of the on-shelf transport along the continental margin. The strongest signal though was found at the Celtic shelf break. This feature is confirmed by our additional RCP8.5 realizations (Figs. 12a and A5a), showing an increase in the standard deviation of surface PO4 concentrations at the Celtic shelf break (Fig. 4, S1) by a factor of 4.4, from 0.027 (1971–2000) to 0.120 mmol m−3 (2101–2150, E0).
We find corresponding anomalies in the future sea level pressure (SLP) field to be related to variations in the strength and orientation of the SLP gradient between the Iceland Low and Azores High (Fig. 13). The drivers for the atmospheric anomalies are sporadic southeastward displacements of the Iceland Low, not represented by any of the leading variability modes of North Atlantic SLP (not shown). The resulting SLP anomalies over the NE Atlantic cause near-surface variations in the meridional wind component, the wind speed and air temperature, affecting the surface heat flux between the ocean and atmosphere and thus the ocean MLD (Fig. A4). In E0, the variability of the atmospheric forcing does not significantly increase through the scenario period. During the shallow-ML regime, however, upper-ocean conditions are more sensitive to variations in the atmospheric forcing. A positive SLP anomaly is associated with southward advection of cold Arctic air masses and enhanced wind speeds. The enhanced surface cooling promotes convection in the upper ocean and deepens the ML. In particular, near the shelf break, the deepened ML entrains warm and saline subpycnocline water to the upper ocean and leads to positive SST and sea surface salinity (SSS) anomalies. Part of the anomalous heat is released to the atmosphere while the salt remains in the water, which increases the upper-ocean density. This mechanism has been shown to feed back to the vertical convection and MLD (M19).
When GIS meltwater is taken into account (E010 and E025), there is a weak increase in the future variability of the atmospheric forcing over the NE Atlantic (Table 3). Moreover, the mean wintertime MLD in the NE Atlantic is more shallow than in E0 (Fig. 5c, d), enhancing the influence of the atmospheric forcing on the MLD (Table 3). However, the impact on the variability of upper-ocean nutrient concentrations is very weak. The main impact is due to the stronger vertical nutrient gradient, resulting from the shallower MLD and higher subpycnocline nutrient concentrations (Sect. 3.1). The stronger nutrient gradient translates MLD variations to enhanced variations in the vertical nutrient supply. Accordingly, for the shallow-ML regime, the standard deviation of surface PO4 concentrations at the Celtic shelf break increases further to 0.148 mmol m−3 for E010 (Fig. 12b) and 0.153 mmol m−3 for E025 (Fig. 12c). Comparing the standard deviations given here and in Table 3, the increase in the variability of upper-ocean nutrient concentrations by 23 % for E010 (compared to E0) decomposes into contributions of only 2 % from the atmospheric forcing and 21 % from the vertical nutrient gradient.
The described mechanisms similarly affect the vertical distribution of salinity in the open ocean as well as the variability of the on-shelf salt flux (Fig. A7). In the NE Atlantic, the freshening of the upper ocean under RCP8.5 leads to the development of a vertical salinity gradient through the pycnocline of 1.20 psu (2101–2150) in E0. In E010 and E025, this gradient strengthens further to 1.65 and 1.94 psu, respectively, due to the additional freshwater hosing at the sea surface. Moreover, in all experiments, the mixing of saline subpycnocline water masses to the upper ocean near the shelf break establishes an ocean–shelf salinity front with higher SSS in the Celtic Sea than in the NE Atlantic (Fig. A2). Finally, the variability of SSS at the Celtic shelf break increases substantially from 0.03 psu (standard deviation 1971–2000) to 0.25 psu (2101–2150) in E0 and further to 0.33 in E010 and 0.36 psu in E025. Unlike nutrient concentrations, however, the stronger increase in the variability of salinity in E010 and E025 is induced by the stronger change (freshening) of the upper ocean. Subpycnocline salinities rather show a weak decrease the more GIS meltwater is added, which is discussed in the next section.
3.3 Changes in the timing of the ocean–shelf exchange regime shift
The shallow-ML regime is fully established when the ML near the shelf break becomes as shallow as the depth of the shelf edge, i.e., about 150–200 m. As can be expected, the impact of GIS meltwater discharge on the stratification strengthening and ML shoaling (Sect. 3.1) implies that the MLD in the NE Atlantic decreases faster when GIS meltwater is added, passing this critical depth earlier in the 21st century (Fig. A6). Nevertheless, the changes in the regime shift timing are surprisingly weak, given the melting rates vary considerably between the experiments. For a melting rate of 0.1 Sv (E010), the impact is only about 10 years (see also PO4 time series shown in Figs. 12 and 16).
The freshwater hosing to the upper ocean leads to a reduction in upper-ocean salinity (Fig. A2c, d). Besides the freshening of the upper ocean, however, the salinity of subpycnocline water masses shows a slight decrease too (Fig. A7). The weakening of the AMOC reduces the northward migration of saline subtropical water and the SPG expansion fosters eastward migration of fresh Labrador Sea and Arctic waters at intermediate depths. This subpycnocline freshening weakens the GIS meltwater impact on the stratification strengthening and ML shoaling, imposing a damping effect on the regime shift timing. Even in a hypothetical extreme scenario of 1.0 Sv meltwater discharge (E100), the shallow-ML regime does not become initiated earlier than the 2070s (Fig. A5d).
3.4 Impact on the biological productivity of the northern North Sea
In M19, the projected changes in the Atlantic nutrient supply to the NWES have been shown to impact the biological productivity in the northern North Sea, as being the largest shelf area primarily influenced by Atlantic inflow. The import of nutrient-rich subpycnocline water masses mixed up at the shelf break mitigates the general weakening in the northern North Sea primary production to be expected from the reduction of upper-ocean nutrient concentrations in the NE Atlantic. Nevertheless, the enhanced variability in the on-shelf nutrient transport also leads to enhanced variations in primary production, as the local influence of the strength of the seasonal stratification decreases to a subordinate factor. In this subsection, we therefore investigate the impact of GIS meltwater discharge on the Atlantic nutrient import and net primary production in the northern North Sea.
A more detailed analysis of the driving mechanisms of the slope current reveals that subtle changes in the northern part of the slope current affect the inflow of Atlantic water masses to the northern North Sea (Fig. 14). Until about the year 2080, the warming and freshening of the upper ocean due to the climate change signal in the atmospheric forcing are stronger on the shallow NWES than in the open NE Atlantic. The ocean–shelf gradient in density and sea surface height (SSH) increases and induces a geostrophic flow along the slope (Simpson and Sharples, 2012; Marsh et al., 2017). Moreover, the westerly winds become stronger over the northern NWES region (Sect. 3.1). The wind component parallel to the shelf break increases by 20 %–30 % and thus contributes further to a temporary strengthening of the slope current. During this period, the enhancing effects of the SSH gradient and wind forcing outweigh the weakening effect of the meridional density gradient (Sect. 3.1), leading to an increase in the Fair Isle inflow to the North Sea and no change in the East Shetland inflow (Fig. 14a, b). These inflows are further maintained by the enhanced wind-driven Ekman transport.
Around 2080, the expansion of the SPG reaches the continental margin and causes sudden decreases in the ocean–shelf gradients of salinity (by 0.22 psu), density (by 0.17 kg m−3) and SSH (by 0.08 m) within one decade. After this period, the weak meridional density gradient west of the shelf break (Fig. 9b) dominates the slope current also in the northern region, and the Fair Isle and East Shetland inflows start to weaken as well.
The inflow along the western side of the Norwegian Trench is mainly governed by topographic steering as the slope current follows the sharp topographic turn to the right at about 62∘ N. Our simulations indicate a substantial weakening of this inflow (Fig. 14c) and are in line with the findings by Holt et al. (2018). The proposed driving mechanism is an increase in the deformation radius due to the strengthening of the stratification. In addition to the general weakening of the slope current, a smaller fraction of the slope current is then able to follow the topographic turn. In our simulations, the permanent stratification around the entrance of the Norwegian Trench strengthens by about 5 g m−4 in experiment E0 and 7 g m−4 in E010 (Fig. A2), and the baroclinic deformation radius increases by about 4–5 km in E0 (similar to 3–4 km in Holt et al., 2018) and 6–7 km in E010. The maximum deformation radii at the end of the simulations are about 10 and 12 km, respectively. During the shallow-ML regime (2101–2150), the correlation between the Norwegian Trench inflow and the deformation radius is about −0.35 (detrended time series) for both experiments. After around 2120, the remaining inflow to the Norwegian Trench does not penetrate further south than about 60∘ N, implying a reduction down to 0 in Fig. 14c.
In Holt et al. (2018), the changes in the stratification and deformation radius are amplified by a strengthening of the circulation in the Nordic Seas, in particular the East Iceland Current, leading to a southward intrusion of cold and fresh Arctic water masses into the northern inflow region of the North Sea. Similarly, in our simulations, the Nordic Seas circulation strengthens and a larger fraction of the East Iceland Current passes the Faroe Islands southeastward, joining the slope current to the north of the Shetland Islands (Fig. A8). Northeast of the Faroe Islands, the East Iceland Current intensifies by about 30 % (in E0 and E010). We thus confirm the enhanced influence of Arctic water masses on the Norwegian Trench inflow projected by Holt et al. (2018). A sudden shutdown due to positive feedback with the accumulation of coastal North Sea water of low salinity, however, does not occur in our simulations. The changes in the circulation are rather smooth over several decades (see Fig. 14c for the Norwegian Trench inflow; not shown for the East Iceland Current) and simulated to happen later than in Holt et al. (2018) by about 30 years. Moreover, a reversal of the Norwegian Trench inflow to a persistent outflow, as also projected by Tinker et al. (2016), is not indicated but may happen after 2150.
Since the future changes in the NE Atlantic circulation and density structure are rather similar in all experiments with and without GIS meltwater discharge, the impact of the freshwater hosing on the slope current and transport rates to the NWES is very low. The main mechanism for the increase in nutrient concentrations on the shelf when GIS meltwater is added thus is the increase in nutrient concentrations of subpycnocline Atlantic water masses (Sect. 3.1). It has been shown by M19 that during the shallow-ML regime the sustained connection to the Atlantic subpycnocline nutrient pool promotes primary production on the NWES. In our new simulations, accordingly, the further increase in subpycnocline nutrient concentrations as a response to added GIS meltwater discharge leads to an increase in pre-bloom nutrient concentrations in the northern North Sea by about 4 % for E010 and 6 % for E025 (relative to E0 2101–2150) and hence to higher net primary production (NPP; Fig. 15a, c). The projected decrease in northern North Sea productivity during the scenario period by 32 % for E0 (2101–2150 minus 1971–2000) therefore weakens to 29 % for E010 and to 28 % for E025. Similarly, the increase in the variability of Atlantic nutrient supply tied to added GIS meltwater leads to higher variability in both pre-bloom nutrient concentrations and net primary production by up to 16 % under the shallow-ML regime (relative to E0 2101–2150) with considerable uncertainty due to internal variability (Fig. 15b, d). These results demonstrate how under decreasing wintertime MLDs the signal transfer from the open ocean to the shelf becomes more sensitive to climate-induced large-scale changes in the North Atlantic mean circulation and interannual variability.
Our simulations indicate a moderate reduction of net on-shelf transport from the NE Atlantic to the NWES under future RCP8.5 conditions. Considering the net winter transport across the continental margin from 44 to 61.5∘ N (excluding the Norwegian Trench inflow), the reduction of the on-shelf component from about 3.12 (1971–2000) to 2.47 Sv (2101–2150) reduces the flushing of the open shelf areas during one winter season from about 45 % to 35 %, fairly independent of the considered amount of GIS meltwater. Nevertheless, even though future on-shelf transport weakens, nutrient concentrations on the NWES do not decline as strongly as in the upper NE Atlantic, owing to mixing processes at the shelf break (in agreement with M19). The maintained connection to nutrient-rich Atlantic subpycnocline water masses gains particular importance as the southward turn of the slope current into the Norwegian Trench is found to weaken substantially during the first half of the 22nd century, thus closing the only direct inflow of deeper Atlantic water to the NWES (in qualitative agreement with Holt et al., 2018). As an implication of the weaker circulation, coastal areas in the southern and eastern North Sea become more influenced by river runoff and nutrient loads. In addition to the decrease in the nutrient transport across the shelf break, this further intensifies the mean nutrient gradient between the inner (dominated by river loads) and outer (dominated by Atlantic inflow) shelf areas as well as its interannual variability. For instance, in our experiments, both the mean pre-bloom nutrient gradient between the northern North Sea and the German Bight and its standard deviation increase by about 20 % during the scenario period (2101–2150 minus 1971–2000).
When GIS meltwater is added, our simulations suggest that the projected nutrient decline in the open shelf areas becomes even weaker due to an increase in Atlantic subpycnocline nutrient concentrations. This increase is induced by the stronger slowdown of the North Atlantic deep circulation as a consequence of a stronger AMOC weakening. Moreover, the subpycnocline nutrient enrichment may be underestimated because our model system does not account for the effect of biologically relevant substances transported into the ocean by meltwater and iceberg calving due to microbial activity and hydrolysis reactions at the interface between land ice and the bedrock (such as dissolved iron, silicate and nitrogen; Bhatia et al., 2013; Duprat et al., 2016; Wadham et al., 2016; Hatton et al., 2019). Part of this additional nutrient input to the upper ocean would be consumed by primary producers and exported to deeper levels. A direct connection of NWES water to Atlantic subpycnocline water is established if MLDs are less than about 130 m, which is the maximum depth of the mean on-shelf transport (Fig. 8a). The MLD passes this threshold only in the most extreme experiment of 1.0 Sv meltwater discharge (E100; Fig. A6d). Accordingly, NWES nutrient concentrations in this experiment show the strongest increase by up to 0.13 mmol P m−3 compared to the reference RCP8.5 simulation without GIS meltwater discharge (E0).
Furthermore, during the shallow-ML regime, the high variability in on-shelf wintertime nutrient fluxes enhances further when GIS meltwater is added. This results mainly from an increase in the vertical nutrient gradient through the pycnocline and is reinforced by a stronger impact of variations in the atmospheric forcing on the shallower ML. In the upper 30 m of the water column, variations in nutrient concentrations at the shelf break are negatively correlated with the on-shelf transport (about −0.5 for E0 during 2101–2150) mediated by the wind-driven advection of nutrient-poor Atlantic surface water. In the deeper water column (> 30 m), by contrast, this correlation is positive (about +0.4) and increases when GIS meltwater is added (to about +0.6 for E010). Thus, the increasing variability in the nutrient concentrations implies an increasing mean nutrient flux to the shelf via the on-shelf transport in the interior water column (Sect. 3.1). In this way, the enhanced variability in the nutrient concentrations at the shelf break contributes to the elevated mean concentrations on the shelf, in addition to the nutrient enrichment of Atlantic subpycnocline water masses. The contribution to the mean nutrient flux (Sect. 3.1) though is only about 3 % for E010.
The SLP anomaly pattern associated with nutrient anomalies at the shelf break is best represented by a non-stationary combination of the North Atlantic Oscillation (NAO) and Scandinavian (SCAN) teleconnection patterns (Chafik et al., 2017). It also resembles the European–North Atlantic sector of the East Atlantic–west Russian pattern, a higher-order atmospheric mode of the Northern Hemisphere which occasionally mixes with the NAO (Barnston and Livezey, 1987; Fagherazzi et al., 2005). The other two characteristic pressure centers over west Siberia and China, however, are not shown in the larger SLP anomaly structure corresponding to Fig. 13 (simulated by MPI-ESM; not shown).
While the strength of the NAO was a major factor in the single realization analyzed by M19, we find a similar meridionally oriented anomaly pattern only in one of the seven simulations presented here (E0, E010, E025). Otherwise, the SLP gradient shows a distinct southwest to northeast orientation, thus dominating the ensemble mean shown in Fig. 13. This means that pronounced anomalies in on-shelf nutrient fluxes are mostly induced by anomalies in the meridional wind component over the NE Atlantic (NAO/SCAN) but can also be triggered by anomalies in the wind speed of the predominant westerlies (NAO).
The impact of GIS meltwater discharge leading to both higher mean pre-bloom nutrient concentrations and annual net primary production in the northern North Sea is significant with respect to the ensemble spreads. The ranges between the ensemble minima and maxima do not overlap among the simulations with and without GIS meltwater discharge (Fig. 15). The increase in the variability, however, is not significant as the ensemble spreads clearly overlap for both variables, nutrient concentration and primary production. This means that while we have gained understanding about the combined effect of increasing radiative forcing and GIS meltwater discharge on physical processes governing the variability of ocean–shelf exchange, whether or not the related changes would have a significant impact on the shelf conditions is strongly influenced by natural variability.
These effects on the ocean–shelf nutrient exchange have been investigated by ensemble simulations with 0.1 and 0.25 Sv meltwater discharge along the coast of Greenland. Future changes in the Greenland surface mass balance and peripheral glaciers projected by ice sheet models show a similarly substantial spread, ranging from no change at all to twice the ensemble mean (Church et al., 2013; Agarwal et al., 2015; Vizcaino et al., 2015). The higher discharge experiment (E025) is additionally motivated by indications that ice sheet models might underestimate future GIS melting rates (Burkett et al., 2014). Moreover, our downscaling model system as well as the parent global model use a fixed GIS grounding line and topography and are therefore not able to account for all positive feedbacks between runoff-induced changes in the ocean and atmosphere and the melting of the GIS (Driesschaert et al., 2007; Castro de la Guardia et al., 2015). They also do not account for additional meltwater discharge due to icebergs, which brings more freshwater to the Labrador Sea and the interior of the SPG (Marson et al., 2018). Besides, in global circulation models (GCMs), the simulated AMOC weakening varies between 10 % and > 60 % after 100 years for a typical freshwater hosing experiment of 0.1 Sv to the northern North Atlantic (e.g., Stouffer et al., 2006; Swingedouw et al., 2015; Yu et al., 2016, see also Sect. 2.3). Reasons for the uncertainty in simulated AMOC responses have been linked to the grid resolution of the ocean model (Swingedouw et al., 2013; Weijer et al., 2012; Böning et al., 2016) and the implemented mixing schemes and convection parameterizations (Nilsson et al., 2003; Yu et al., 2008; Marzeion and Levermann, 2009), affecting the simulated export of the added freshwater from the SPG region to adjacent ocean basins. These uncertainties related to GIS meltwater discharge have to be assessed in addition to the uncertainties in the atmospherically driven AMOC weakening and stratification strengthening deduced from multi-model ensembles. Under scenario RCP8.5, the AMOC weakening by the year 2100 ranges between 12 % and 54 % of the individual CMIP5 model's historical mean (Weaver et al., 2012; Cheng et al., 2013), and the upper-ocean stratification increases by 6 %–30 % (Fu et al., 2016) with a similar uncertainty in MLD shoaling (Bopp et al., 2001; Capotondi et al., 2012; Yool et al., 2015; Alexander et al., 2018).
Though in the present study we have not utilized a GCM of typical CMIP or IPCC configurations, we expect comparable uncertainties in the projected AMOC weakening, stratification strengthening and MLD decrease in our downscaling simulations due to the influences of the involved particular model components and the particular forcing global model. These uncertainties are mainly associated with the magnitudes of the simulated changes, since the directions of the changes are consistent within the cited multi-model ensembles. Low uncertainty is therefore associated with the qualitative impacts of GIS meltwater discharge and the involved driving mechanisms, i.e., the general intensification of the regime shift due to stronger stratification, shallower MLD and higher subpycnocline nutrient concentrations. Nevertheless, relative to the change signals, the biases and uncertainties in our model results are partly of opposite directions. The simulated North Atlantic circulation is more sensitive to the atmospheric forcing but less sensitive to additional freshwater input, compared to other climate system models (Sect. 2.3). Present-day wintertime MLD seems to be overestimated but there is no consistent relationship between MLD biases and other variables such as water density or salinity (Heuzé, 2017). It is therefore hard to speculate about different regime shift characteristics as simulated by other climate system models, calling for the need to conduct further downscaling experiments and intermodel comparisons.
In the parent global model (MPI-ESM-LR), the persistent dominance of Atlantic upper-ocean conditions on shaping shelf conditions shown by M19 is also seen in the respective freshwater hosing experiments carried out in preparation of the downscaling simulations presented here (Fig. 16). As the future stratification strengthens and the MLD decreases when GIS meltwater is added, surface nutrient concentrations in the NE Atlantic decrease to lower and lower levels, driving a stronger nutrient decline on the NWES as well. The meltwater impact though is generally very weak, as has also been found, e.g., by Zickfeld et al. (2008) and Kuhlbrodt et al. (2009). The regime shift in the on-shelf nutrient transport cannot be captured by the global model because of its coarse grid resolution (50–70 km at the shelf break) and the neglect of intertidal currents.
We have made recent RCP8.5 downscaling simulations for the physical and biogeochemical state of the NWES (Mathis et al., 2018, 2019) more consistent by taking into account the effects of GIS meltwater discharge. The considered GIS melting rates cover a wide range projected by ice sheet models (Burkett et al., 2014). Furthermore, we gained knowledge about the influence of the parent global model's internal variability on the characteristics of the proposed regime shift in Atlantic nutrient supply by utilizing an ensemble of three realizations for experiments with 0 and 0.1 Sv GIS melting rates.
We found that for the high-emission scenario RCP8.5, GIS meltwater discharge influences the nutrient supply from the NE Atlantic to the NWES. It potentially intensifies the regime shift near the end of the 21st century in terms of spatial distribution and temporal variability and leads to an earlier onset, depending on GIS melting rates. Nevertheless, most striking features associated with the regime shift are induced by changes in upper-ocean conditions of the NE Atlantic tied to a critical reduction of the wintertime MLD, a slowdown of the AMOC and a widening of the SPG. These are pronounced changes already in the RCP8.5 simulations without additional GIS meltwater discharge, primarily caused by a thermally induced buoyancy reduction (Liu et al., 2017). The impact of GIS meltwater discharge on the regime shift characteristics is therefore not as substantial as the general difference between the present-day deep-ML and future shallow-ML regimes. The simulated intensity and timing of the regime shift, however, are subject to uncertainties in the projected GIS melting rate, AMOC weakening and NE Atlantic stratification strengthening and MLD shoaling. In particular, it has been demonstrated that under sufficiently shallow wintertime MLDs, changes in the upper-ocean conditions of the NE Atlantic are no longer transmitted to the NWES with an intimate correspondence. Rather, the conditions of on-shelf transport across the shelf break become heavily modified by subpycnocline water masses and hence controlled by the interplay between the deeper ocean, the ocean mixed layer and the atmospheric forcing.
In open shelf areas, the intensification of the regime shift due to GIS meltwater discharge leads to higher mean pre-bloom concentrations and fuels higher productivity during the following plankton growing season (relative to experiments without GIS meltwater discharge). The impact on the variability in nutrient concentrations and productivity, however, is low. The main reasons are a strong influence of less variable Atlantic subpycnocline water masses on the net on-shelf nutrient transport and the influence of local physical conditions of the upper water column on the shelf (SST and depth of seasonal stratification) on the plankton spring bloom and summer growth (M19). The ensemble approach allowed us to render the impact on the variability insignificant with respect to the ensemble spread. The higher mean pre-bloom concentrations and annual primary production, by contrast, were revealed to be robust signals.
The increasing influence of Atlantic subpycnocline water masses on the NWES physical and biogeochemical state affects not only the supply of oceanic nutrients to the shelf. Any conservative constituent of the on-shelf transport potentially has an impact on the NWES marine ecosystem, such as salinity, temperature, oxygen and alkalinity. Further simulations with regional ecosystem models driven by boundary conditions from high-resolution future projections as presented here are therefore highly appreciated to more comprehensively assess the impact of the shallow-ML regime on the NWES biological environment and the efficacy of the shelf carbon pump.
Primary data and scripts used in the analysis that may be useful in reproducing the work are archived by the Max-Planck-Institute for Meteorology and can be obtained by contacting firstname.lastname@example.org.
MM and UM jointly conceived and designed the model experiments. Simulations and analyses were carried out by MM. Results and conclusions were developed and discussed by both authors. The manuscript was written by MM with contributions from UM.
The authors declare that they have no conflict of interest.
This work was funded by the BMBF-funded joint research projects RACE – Regional Atlantic Circulation and Global Change and RACE – Synthesis. Computational resources were made available by the German Climate Computing Center (DKRZ) through support from the German Federal Ministry of Education and Research (BMBF). The used GIS meltwater climatology was provided by Torge Martin from the Helmholtz Centre for Ocean Research Kiel (GEOMAR).
This research has been supported by the German Federal Ministry of Education and Research (BMBF) (grant nos. 03F0729D and 03F0824D).
The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.
This paper was edited by Markus Meier and reviewed by two anonymous referees.
Agarwal, N., Jungclaus, J. H., Köhl, A., Mechoso, C. R., and Stammer, D.: Additional contributions to CMIP5 regional sea level projections resulting from Greenland and Antarctic ice mass loss, Environ. Res. Lett., 10, 074008, https://doi.org/10.1088/1748-9326/10/7/074008, 2015. a, b
Alexander, M. A., Scott, J. D., Friedland, K. D., Mills, K. E., Nye, J. A., Pershing, A. J., and Thomas, A. C.: Projected sea surface temperatures over the 21st century: Changes in the mean, variability and extremes for large marine ecosystem regions of Northern Oceans, Elem. Sci. Anth., 6, https://doi.org/10.1525/elementa.191, 2018. a, b
An, L., Rignot, E., Elieff, S., Morlighem, M., Millan, R., Mouginot, J., Holland, D. M., Holland, D., and Paden, J.: Bed elevation of Jakobshavn Isbræ, West Greenland, from high‐resolution airborne gravity and other data, Geophys. Res. Lett., 44, 3728–3736, 2017. a
Bamber, J., van den Broeke, M., Ettema, J., Lenaerts, J., and Rignot, E.: Recent large increases in freshwater fluxes from Greenland into the North Atlantic, Geophys. Res. Lett., 39, L19501, https://doi.org/10.1029/2012GL052552, 2012. a, b, c
Bamber, J. L., Tedstone, A. J., King, M. D., Howat, I. M., Enderlin, E., van den Broeke, M. R., and Noel, B.: Land ice freshwater budget of the Arctic and North Atlantic Oceans: 1. Data, methods, and results, J. Geophys. Res.-Oceans, 123, 1827–1837, 2018. a
Barnston, A. G. and Livezey, R. E.: Classification, seasonality and persistence of low‐frequency atmospheric circulation patterns, Mon. Weather Rev., 115, 1083–1126, 1987. a
Bhatia, M., Kujawinski, E. B., Das, S. B., Breier, C. F., Henderson, P. B., and Charette, M. A.: Greenland meltwater as a significant and potentially bioavailable source of iron to the ocean, Nat. Geosci., 6, 274–278, 2013. a
Böning, C. W., Behrens, E., Biastoch, A., Getzlaff, K., and Bamber, J. L.: Emerging impact of Greenland meltwater on deepwater formation in the North Atlantic Ocean, Nat. Geosci., 9, 523–527, 2016. a, b
Bopp, L., Monfray, P., Aumont, O., Dufresne, J. L., Treut, H. L., Madec, G., Terray, L., and Orr, J. C.: Potential impact of climate change on marine export production, Global Biogeochem. Cy., 15, 81–99, 2001. a
Bopp, L., Resplandy, L., Orr, J. C., Doney, S. C., Dunne, J. P., Gehlen, M., Halloran, P., Heinze, C., Ilyina, T., Séférian, R., Tjiputra, J., and Vichi, M.: Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models, Biogeosciences, 10, 6225–6245, https://doi.org/10.5194/bg-10-6225-2013, 2013. a
Burkett, V. R., Suarez, A. G., Bindi, M., Conde, C., Mukerji, R., Prather, M. J., Clair, A. L. S., and Yohe, G. W.: Point of departure, in: Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Field, C. B., Barros, V. R., Dokken, D. J., Mach, K. J., Mastrandrea, M. D., Bilir, T. E., Chatterjee, M., Ebi, K. L., Estrada, Y. O., Genova, R. C., Girma, B., Kissel, E. S., Levy, A. N., MacCracken, S., Mastrandrea, P. R., and White, L. L., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 169–194, 2014. a, b
Capotondi, A., Alexander, M. A., Bond, N. A., Curchitser, E. N., and Scott, J. D.: Enhanced upper ocean stratification with climate change in the CMIP3 models, J. Geophys. Res., 117, C04031, https://doi.org/10.1029/2011JC007409, 2012. a
Castro de la Guardia, L., Hu, X., and Myers, P. G.: Potential positive feedback between Greenland Ice Sheet melt and Baffin Bay heat content on the west Greenland shelf, Geophys. Res. Lett., 42, 4922–4930, 2015. a
Chafik, L., Nilsen, J. E. O., and Dangendorf, S.: Impact of North Atlantic Teleconnection Patterns on Northern European Sea Level, Journal of Marine Science and Engineering, 5, https://doi.org/10.3390/jmse5030043, 2017. a
Cheng, W., Chiang, J. C., and Zhang, D.: Atlantic Meridional Overturning Circulation (AMOC) in CMIP5 Models: RCP and Historical Simulations, J. Climate, 26, 7187–7197, 2013. a
Church, J. A., Clark, P. U., Cazenave, A., Gregory, J. M., Jevrejeva, S., Levermann, A., Merrifield, M. A., Milne, G. A., Nerem, R. S., Nunn, P. D., Payne, A. J., Pfeffer, W. T., Stammer, D., and Unnikrishnan, A. S.: Sea Level Change, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G. K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1137–1216, 2013. a, b, c
Collins, M., Knutti, R., Arblaster, J., Dufresne, J. L., Fichefet, T., Friedlingstein, P., Gao, X., Gutowski, W. J., Johns, T., Krinner, G., Shongwe, M., Tebaldi, C., Weaver, A. J., and Wehner, M.: Long-term Climate Change: Projections, Commitments and Irreversibility, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G. K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1029–1136, 2013. a
Costoya, X., deCastro, M., Gómez-Gesteira, M., and Santos, F.: Mixed layer depth trends in the Bay of Biscay over the period 1975–2010, PLoS ONE, 9, e99321, https://doi.org/10.1371/journal.pone.0099321, 2014. a
Courtois, P., Hu, X., Pennelly, C., Spence, P., and Myers, P. G.: Mixed layer depth calculation in deep convection regions in ocean numerical models, Ocean Model., 120, 60–78, 2017. a
de Boyer Montégut, C., Madec, G., Fischer, A. S., Lazar, A., and Iudicone, D.: Mixed layer depth over the global ocean: An examination of profile data and a profile‐based climatology, J. Geophys. Res., 109, C12003, https://doi.org/10.1029/2004JC002378, 2004. a
Driesschaert, E., Fichefet, T., Goosse, H., Huybrechts, P., Janssens, I., Mouchet, A., Munhoven, G., Brovkin, V., and Weber, S. L.: Modeling the influence of Greenland ice sheet melting on the Atlantic meridional overturning circulation during the next millennia, Geophys. Res. Lett., 34, L10707, https://doi.org/10.1029/2007GL029516, 2007. a, b
Duprat, L. P. A. M., Bigg, G. R., and Wilton, D. J.: Enhanced Southern Ocean marine productivity due to fertilization by giant icebergs, Nat. Geosci., 9, 219–221, 2016. a
Fischer, M., Domeisen, D. I. V., Müller, W. A., and Baehr, J.: Changes in the seasonal cycle of the Atlantic meridional heat transport in a RCP 8.5 climate projection in MPI-ESM, Earth Syst. Dynam., 8, 129–146, https://doi.org/10.5194/esd-8-129-2017, 2017. a
Fu, W., Randerson, J. T., and Moore, J. K.: Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in the CMIP5 models, Biogeosciences, 13, 5151–5170, https://doi.org/10.5194/bg-13-5151-2016, 2016. a, b, c
Fu, Y., Karstensen, J., and Brandt, P.: Atlantic Meridional Overturning Circulation at 14.5∘ N in 1989 and 2013 and 24.5∘ N in 1992 and 2015: volume, heat, and freshwater transports, Ocean Sci., 14, 589–616, https://doi.org/10.5194/os-14-589-2018, 2018. a
Giorgetta, M. A., Jungclaus, J., Reick, C. H., Legutke, S., Bader, J., Böttinger, M., Brovkin, V., Crueger, T., Esch, M., Fieg, K., Glushak, K., Gayler, V., Haak, H., Hollweg, H. D., Ilyina, T., Kinne, S., Kornblueh, L., Matei, D., Mauritsen, T., Mikolajewicz, U., Mueller, W., Notz, D., Pithan, F., Raddatz, T., Rast, S., Redler, R., Roeckner, E., Schmidt, H., Schnur, R., Segschneider, J., Six, K. D., Stockhause, M., Timmreck, C., Wegner, J., Widmann, H., Wieners, K. H., Claussen, M., Marotzke, J., and Stevens, B.: Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5, J. Adv. Model. Earth Sy., 5, 572–597, 2013. a
Graham, J. A., Rosser, J. P., O'Dea, E., and Hewitt, H. T.: Resolving Shelf Break Exchange Around the European Northwest Shelf, Geophys. Res. Lett., 45, 12386–12395, 2018. a
Gröger, M., Maier-Reimer, E., Mikolajewicz, U., Moll, A., and Sein, D.: NW European shelf under climate warming: implications for open ocean – shelf exchange, primary production, and carbon absorption, Biogeosciences, 10, 3767–3792, https://doi.org/10.5194/bg-10-3767-2013, 2013. a, b, c, d, e
Hagemann, S. and Dümenil-Gates, L.: Validation of the hydrological cycle of ECMWF and NCEP reanalyses using the MPI hydrological discharge model, J. Geophys. Res.-Atmos., 106, 1503–1510, 2001. a
Hatton, J. E., Hendry, K. R., Hawkings, J. R., Wadham, J. L., Opfergelt, S., Kohler, T. J., Yde, J. C., Stibal, M., and Žárský, J. D.: Silicon isotopes in Arctic and sub-Arctic glacial meltwaters: the role of subglacial weathering in the silicon cycle, P. Roy. Soc. A, 475, https://doi.org/10.1098/rspa.2019.0098, 2019. a
Hátún, H., Azetsu-Scott, K., Somavilla, R., Rey, F., Johnson, C., Mathis, M., Mikolajewicz, U., Coupel, P., Tremblay, J. É., Hartman, S., Pacariz, S. V., Salter, I., and Ólafsson, J.: The subpolar gyre regulates silicate concentrations in the North Atlantic, Sci. Rep., 7, https://doi.org/10.1038/s41598-017-14837-4, 2017. a
Heinze, C., Maier-Reimer, E., Winguth, A. M. E., and Archer, D.: A global oceanic sediment model for long-term climate studies, Global Biogeochem. Cy., 13, 221–250, 1999. a
Hemming, S. R.: Heinrich events: Massive late Pleistocene detritus layers of the North Atlantic and their global climate imprint, Rev. Geophys., 42, RG1005, https://doi.org/10.1029/2003RG000128, 2004. a
Holt, J., Butenschön, M., Wakelin, S. L., Artioli, Y., and Allen, J. I.: Oceanic controls on the primary production of the northwest European continental shelf: model experiments under recent past conditions and a potential future scenario, Biogeosciences, 9, 97–117, https://doi.org/10.5194/bg-9-97-2012, 2012. a, b, c
Holt, J., Schrum, C., Cannaby, H., Daewel, U., Allen, I., Artioli, Y., Bopp, L., Butenschon, M., Fach, B. A., Harle, J., Pushpadas, D., Salihoglu, B., and Wakelin, S.: Potential impacts of climate change on the primary production of regional seas: A comparative analysis of five European seas, Prog. Oceanogr., 140, 91–115, 2016. a
Holt, J., Polton, J., Huthnance, J., Wakelin, S., O'Dea, E., Harle, J., Yool, A., Artioli, Y., Blackford, J., Siddorn, J., and Inall, M.: Climate-driven change in the North Atlantic and Arctic Oceans can greatly reduce the circulation of the North Sea, Geophys. Res. Lett., 45, 11827–11836, 2018. a, b, c, d, e, f
Hu, A., Meehl, G., Han, W., and Yin, J.: Transient response of the MOC and climate to potential melting of the Greenland ice sheet in the 21st century, Geophys. Res. Lett., 36, L10707, https://doi.org/10.1029/2009GL037998, 2009. a, b
Hu, Y., Tao, L., and Liu, J.: Poleward expansion of the Hadley circulation in CMIP5 simulations, Adv. Atmos. Sci., 30, 790–795, 2013. a
Ilyina, T., Six, K. D., Segschneider, J., Maier-Reimer, E., Li, H., and Núñez-Riboni, I.: The global ocean biogeochemistry model HAMOCC: Model architecture and performance as component of the MPI-Earth System Model in different CMIP5 experimental realizations, J. Adv. Model. Earth Sy., 5, 287–315, 2013. a
Jacob, D. and Podzun, R.: Sensitivity studies with the regional climate model REMO, Meteorol. Atmos. Phys., 63, 119–129, 1997. a
Jacob, D., van den Hurk, B. J. J. M., Andræ, U., Elgered, G., Fortelius, C., Graham, L. P., Jackson, S. D., Karstens, U., Köpken, C., Lindau, R., Podzun, R., Rockel, B., Rubel, F., Sass, B. H., Smith, R. N. B., and Yang, X.: A comprehensive model inter-comparison study investigating the water budget during the BALTEX-PIDCAP period, Meteorol. Atmos. Phys., 77, 19–43, 2001. a
Jongma, J. I., Renssen, H., and Roche, D. M.: Simulating Heinrichevent 1 with interactive icebergs, Clim. Dynam., 40, 1373–1385, 2013. a
Jungclaus, J., Fischer, N., Haak, H., Lohmann, K., Marotzke, J., Matei, D., Mikolajewicz, U., Notz, D., and von Storch, J. S.: Characteristics of the ocean simulations in MPIOM, the ocean component of the MPI Earth System Model, J. Adv. Model. Earth Sy., 5, 422–446, 2013. a
Jungclaus, J. H., Haak, H., Esch, M., Roeckner, E., and Marotzke, J.: Will Greenland melting halt the thermohaline circulation?, Geophys. Res. Lett., 33, L17708, https://doi.org/10.1029/2006GL026815, 2006. a, b
Kanzow, T., Cunningham, S. A., Johns, W. E., Hirschi, J. J., Marotzke, J., Baringer, M. O., Meinen, C. S., Chidichimo, M. P., Atkinson, C., Beal, L. M., Bryden, H. L., and Collins, J.: Seasonal Variability of the Atlantic Meridional Overturning Circulation at 26.5∘ N, J. Climate, 23, 5678–5698, 2010. a
Keeling, R. F., Körtzinger, A., and Gruber, N.: Ocean deoxygenation in a warming world, Annu. Rev. Mar. Sci., 2, 199–229, 2010. a
Kloster, S., Feichter, J., Maier-Reimer, E., Six, K. D., Stier, P., and Wetzel, P.: DMS cycle in the marine ocean-atmosphere system – a global model study, Biogeosciences, 3, 29–51, https://doi.org/10.5194/bg-3-29-2006, 2006. a
Kotlarski, S., Keuler, K., Christensen, O. B., Colette, A., Déqué, M., Gobiet, A., Goergen, K., Jacob, D., Lüthi, D., van Meijgaard, E., Nikulin, G., Schär, C., Teichmann, C., Vautard, R., Warrach-Sagi, K., and Wulfmeyer, V.: Regional climate modeling on European scales: a joint standard evaluation of the EURO-CORDEX RCM ensemble, Geosci. Model Dev., 7, 1297–1333, https://doi.org/10.5194/gmd-7-1297-2014, 2014. a
Kuhlbrodt, T., Rahmstorf, S., Zickfeld, K., Vikebø, F. B., Sundby, S., Hofmann, M., Link, P. M., Bondeau, A., Cramer, W., and Jaeger, C.: An integrated assessment of changes in the thermohaline circulation, Climatic Change, 96, 489–537, 2009. a, b, c, d, e, f
Kwiatkowski, L., Naar, J., Bopp, L., Aumont, O., Defrance, D., and Couespel, D.: Decline in Atlantic Primary Production Accelerated by Greenland Ice Sheet Melt, Geophys. Res. Lett., 46, https://doi.org/10.1029/2019GL085267, 2019. a, b
Liu, K. K., Atkinson, L., Quiñones, R., and Talaue-McManus, L.: Biogeochemistry of continental margins in a global context, in: Carbon and nutrient fluxes in continental margins: A global synthesis, edited by: Liu, K. K., Atkinson, L., Quiñones, R., and Talaue-McManus, L., Springer, Berlin, 3–24, 2010. a
Liu, W., Xie, S. P., Liu, Z., and Zhu, J.: Overlooked possibility of a collapsed Atlantic meridional overturning circulation in warming climate, Sci. Adv., 3, e1601666, https://doi.org/10.1126/sciadv.1601666, 2017. a, b
Liu, Y., Hallberg, R., Sergienko, O., Samuels, B. L., Harrison, M., and Oppenheimer, M.: Climate response to the meltwater runoff from Greenland ice sheet: evolving sensitivity to discharging locations, Clim. Dynam., 51, 1733–1751, 2018. a
Luthcke, S. B., Zwally, H. J., Abdalati, W., Rowlands, D. D., Ray, R. D., Nerem, R. S., Lemoine, F. G., McCarthy, J. J., and Chinn, D. S.: Recent Greenland ice mass loss by drainage system from satellite gravity observations, Science, 314, 1286–1289, 2006. a
Maier-Reimer, E.: Design of a closed boundary regional model of the Arctic Ocean, Workshop on polar processes in global climate, Amer. Meteor. Soc., Boston, 72–73, 1997. a
Maier-Reimer, E., Kriest, I., Segschneider, J., and Wetzel, P.: The Hamburg Ocean Carbon Cycle Model HAMOCC5.1 – Technical Description Release 1.1, Berichte zur Erdsystemforschung, 14, 50 pp., Max Planck Institute for Meteorology, Hamburg, Germany, 2005. a
Marsh, R., Haigh, I. D., Cunningham, S. A., Inall, M. E., Porter, M., and Moat, B. I.: Large-scale forcing of the European Slope Current and associated inflows to the North Sea, Ocean Sci., 13, 315–335, https://doi.org/10.5194/os-13-315-2017, 2017. a, b
Marsland, S. J., Haak, H., Jungclaus, J. H., Latif, M., and Roeske, F.: The Max-Planck-Institute global ocean/sea ice model with orthogonal curvilinear coordinates, Ocean Model., 5, 91–127, 2003. a, b
Marson, J. M., Myers, P. G., Hu, X., and Sommer, J. L.: Using vertically integrated ocean fields to characterize Greenland icebergs' distribution and lifetime, Geophys. Res. Lett., 45, 4208–4217, 2018. a
Martin, T., Biastoch, A., Harlaß, J., Lohmann, G., and U. Mikolajewicz, X. W.: Weak response to enhanced Greenland melting by complex climate models with non-eddying ocean, in preparation, 2019. a, b, c, d
Marzeion, B. and Levermann, A.: Stratification-dependent mixing may increase sensitivity of a wind-driven Atlantic overturning to surface freshwater flux, Geophys. Res. Lett., 36, L20602, https://doi.org/10.1029/2009GL039947, 2009. a
Mathis, M., Elizalde, A., Mikolajewicz, U., and Pohlmann, T.: Variability patterns of the general circulation and sea water temperature in the North Sea, Prog. Oceanogr., 135, 91–112, 2015. a
Mathis, M., Elizalde, A., and Mikolajewicz, U.: Which complexity of regional climate system models is essential for downscaling anthropogenic climate change in the Northwest European Shelf?, Clim. Dynam., 50, 2637–2659, 2018. a, b, c, d, e, f, g
Mikolajewicz, U., Vizcaíno, M., Jungclaus, J., and Schurgers, G.: Effect of ice sheet interactions in anthropogenic climate change simulations, Geophys. Res. Lett., 34, L18706, https://doi.org/10.1029/2007GL031173, 2007. a, b, c
Morlighem, M., Williams, C. N., Rignot, E., An, L., Arndt, J. E., Bamber, J. L., Catania, G., Chauché, N., Dowdeswell, J. A., Dorschel, B., Fenty, I., Hogan, K., Howat, I., Hubbard, A., Jakobsson, M., Jordan, T. M., Kjeldsen, K. K., Millan, R., Mayer, L., Mouginot, J., Noël, B. P. Y., O'Cofaigh, C., Palmer, S., Rysgaard, S., Seroussi, H., Siegert, M. J., Slabon, P., Straneo, F., van den Broeke, M. R., Weinrebe, W., Wood, M., and Zinglersen, K. B.: BedMachine v3: Complete bed topography and ocean bathymetry mapping of Greenland from multi‐beam echo sounding combined with mass conservation, Geophys. Res. Lett., 44, 418–422, 2017. a
Nilsson, J., Broström, G., and Walin, G.: The thermohaline circulation and vertical mixing: does weaker density stratification give stronger overturning?, J. Phys. Oceanogr., 33, 2781–2795, 2003. a
Oliver, H., Luo, H., Castelao, R. M., van Dijken, G. L., Mattingly, K. S., Rosen, J. J., Mote, T. L., Arrigo, K. R., Rennermalm, A. K., Tedesco, M., and Yager, P. L.: Exploring the potential impact of Greenland meltwater on stratification, photosynthetically active radiation, and primary production in the Labrador Sea, J. Geophys. Res.-Oceans, 123, 2570–2591, 2018. a
Pacanowski, R. C. and Philander, S. G. H.: Parameterization of vertical mixing in numerical models of tropical oceans, J. Phys. Ocean., 11, 1443–1451, 1981. a
Pätsch, J., Burchard, H., Dieterich, C., Gräwe, U., Gröger, M., Mathis, M., Kapitza, H., Bersch, M., Moll, A., Pohlmann, T., Su, J., Ho-Hagemann, H. T. M., Schulz, A., Elizalde, A., and Eden, C.: An evaluation of the North Sea circulation in global and regional models relevant for ecosystem simulations, Ocean Model., 116, 70–95, 2017. a
Ruiz-Castillo, E., Sharples, J., Hopkins, J., and Woodward, M.: Seasonality in the cross-shelf physical structure of a temperate shelf sea and the implications for nitrate supply, Prog. Oceanogr., 177, 101985, https://doi.org/10.1016/j.pocean.2018.07.006, 2019. a
Schiller, A., Mikolajewicz, U., and Voss, R.: The stability of the North Atlantic thermohaline circulation in a coupled ocean-atmosphere general circulation model, Clim. Dynam., 13, 325–347, 1997. a
Schrum, C., Lowe, J., Meier, H. E. M., Grabemann, I., Holt, J., Mathis, M., Pohlmann, T., Skogen, M. D., Sterl, A., and Wakelin, S.: Projected Change – North Sea, in: North Sea Region Climate Change Assessment, edited by: Quante, M. and Colijn, F., Springer, Berlin, Heidelberg, 175–217, 2016. a, b
Sévellec, F., Fedorov, A. V., and Liu, W.: Arctic sea-ice decline weakens the Atlantic meridional overturning circulation, Nat. Clim. Change, 7, 604–610, 2017. a
Stouffer, R. J., Yin, J., Gregory, J. M., Dixon, K. W., Spelman, M. J., Hurlin, W., Weaver, A. J., Eby, M., Flato, G. M., Hasumi, H., Hu, A., Jungclaus, J. H., Kamenkovich, I. V., Levermann, A., Montoya, M., Murakami, S., Nawrath, S., Oka, A., Peltier, W. R., Robitaille, D. Y., Sokolov, A., Vettoretti, G., and Weber, S. L.: Investigating the Causes of the Response of the Thermohaline Circulation to Past and Future Climate Changes, J. Climate, 19, 1365–1387, 2006. a, b, c, d
Suga, T., Motoki, K., Aoki, Y., and Macdonald, A. M.: The North Pacific climatology of winter mixed layer and mode waters, J. Phys. Oceanogr., 34, 3–22, 2004. a
Sweby, P. K.: High resolution schemes using flux limiters for hyperbolic conservation laws, Siam J. Numer. Anal., 21, 995–1011, 1984. a
Swingedouw, D., Rodehacke, C. B., Behrens, E., Menary, M., Olsen, S. M., Gao, Y., Mikolajewicz, U., Mignot, J., and Biastoch, A.: Decadal fingerprints of freshwater discharge around Greenland in a multi-model ensemble, Clim. Dynam., 41, 695–720, 2013. a, b, c
Swingedouw, D., Rodehacke, C. B., Olsen, S. M., Menary, M., Gao, Y., Mikolajewicz, U., and Mignot, J.: On the reduced sensitivity of the Atlantic overturning to Greenland ice sheet melting in projections: a multi-model assessment, Clim. Dynam., 44, 3261–3279, 2015. a, b, c
Talley, L. D., Reid, J. L., and Robbins, P. E.: Data-based meridional overturning streamfunctions for the global ocean, J. Climate, 16, 3213–3226, 2003. a
Thomas, H., Bozec, Y., de Baar, H., Elkalay, K., Frankignoulle, M., Kühn, W., Lenhart, H., Moll, A., Pätsch, J., Radach, G., Schiettecatte, L. S., and Borges, A. V.: Carbon and nutrient budgets of the North Sea, in: Carbon and nutrient fluxes in continental margins: A global synthesis, edited by: Liu, K. K., Atkinson, L., Quiñones, R., and Talaue-McManus, L., Springer, Berlin, 346–355, 2010. a
Thomas, M., Sündermann, J., and Maier-Reimer, E.: Consideration of ocean tides in an OGCM and impacts on subseasonal to decadal polar motion excitation, Geophys. Res. Lett., 28, 2457–2460, 2001. a
Tinker, J., Lowe, J., Pardaens, A., Holt, J., and Barciela, R.: Uncertainty in climate projections for the 21st century northwest European shelf seas, Prog. Oceanogr., 148, 56–73, 2016. a
van den Broeke, M., Box, J., Fettweis, X., Hanna, E., Noël, B., Tedesco, M., van As, D., van de Berg, W. J., and van Kampenhout, L.: Greenland Ice Sheet Surface Mass Loss: Recent Developments in Observation and Modeling, Current Climate Change Reports, 3, 345–356, 2017. a
Vaughan, D., Comiso, J., Allison, I., Carrasco, J., Kaser, G., Kwok, R., Mote, P., Murray, T., Paul, F., Ren, J., Rignot, E., Solomina, O., Steffen, K., and Zhang, T.: Observations: Cryosphere, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G. K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, UK and New York, NY, USA, 317–382, 2013. a
Vermaat, J. E., McQuatters-Gollop, A., Eleveld, M. A., and Gilbert, A. J.: Past, present and future nutrient loads of the North Sea: Causes and consequences, Estuar. Coast. Shelf Sci., 80, 53–59, 2008. a
Vizcaíno, M., Mikolajewicz, U., Jungclaus, J., and Schurgers, G.: Climate modification by future ice sheet changes and consequences for ice sheet mass balance, Clim. Dynam., 34, 301–324, 2010. a
Vizcaino, M., Mikolajewicz, U., Ziemen, F., Rodehacke, C. B., Greve, R., and van den Broeke, M. R.: Coupled simulations of Greenland Ice Sheet and climate change up to A.D. 2300, Geophys. Res. Lett., 42, 3927–3935, 2015. a, b
Wadham, J. L., Hawkings, J., Telling, J., Chandler, D., Alcock, J., O'Donnell, E., Kaur, P., Bagshaw, E., Tranter, M., Tedstone, A., and Nienow, P.: Sources, cycling and export of nitrogen on the Greenland Ice Sheet, Biogeosciences, 13, 6339–6352, https://doi.org/10.5194/bg-13-6339-2016, 2016. a
Wakelin, S. L., Artioli, Y., Butenschön, M., Allen, J. I., and Holt, J. T.: Modelling the combined impacts of climate change and direct anthropogenic drivers on the ecosystem of the northwest European continental shelf, J. Marine Syst., 152, 51–63, 2015. a, b
Weaver, A. J., Sedláček, J., Eby, M., Alexander, K., Crespin, E., Fichefet, T., Philippon-Berthier, G., Joos, F., Kawamiya, M., Matsumoto, K., Steinacher, M., Tachiiri, K., Tokos, K., Yoshimori, M., and Zickfeld, K.: Stability of the Atlantic meridional overturning circulation: A model intercomparison, Geophys. Res. Lett., 39, L20709, https://doi.org/10.1029/2012GL053763, 2012. a
Weijer, W., Maltrud, M. E., Hecht, M. W., Dijkstra, H. A., and Kliphuis, M. A.: Response of the Atlantic Ocean circulation to Greenland Ice Sheet melting in a strongly-eddying ocean model, Geophys. Res. Lett., 39, L09606, https://doi.org/10.1029/2012GL051611, 2012. a, b
Wilcock, P.: Sediment Transport Seminar, University of California at Berkeley, Lectures 1 and 3, available at: http://calm.geo.berkeley.edu/geomorph/wilcock/wilcock.html (last access: 17 January 2020), 2004. a
Williams, R. G., McDonagh, E., Roussenov, V. M., Torres-Valdes, S., King, B., Sanders, R., and Hansell, D. A.: Nutrient streams in the North Atlantic: Advective pathways of inorganic and dissolved organic nutrients, Global Biogeochem. Cy., 25, 16, https://doi.org/10.1029/2010GB003853, 2011. a
Woollings, T., Gregory, J. M., Pinto, J. G., Reyers, M., and Brayshaw, D. J.: Response of the North Atlantic storm track to climate change shaped by ocean-atmosphere coupling, Nat. Geosci., 5, 313–317, 2012. a, b
Wouters, B., Drijfhout, S. S., and Hazeleger, W.: Interdecadal North Atlantic meridional overturning circulation variability in EC-EARTH, Clim. Dynam., 39, 2695–2712, 2012. a
Yang, Q., Dixon, T. H., Myers, P. G., Bonin, J., Chambers, D., van den Broeke, M. R., Ribergaard, M. H., and Mortensen, J.: Recent increases in Arctic freshwater flux affects Labrador Sea convection and Atlantic overturning circulation, Nat. Commun., 7, 10525, https://doi.org/10.1038/ncomms10525, 2016. a
Yu, L., Gao, Y. Q., Wang, H. J., and Drange, H.: Revisiting effect of ocean diapycnal mixing on atlantic meridional overturning circulation recovery in a freshwater perturbation simulation, Adv. Atmos. Sci., 25, 597–609, 2008. a
Yu, L., Gao, Y., and Otterå, O. H.: The sensitivity of the Atlantic meridional overturning circulation to enhanced freshwater discharge along the entire, eastern and western coast of Greenland, Clim. Dynam., 46, 1351–1369, 2016. a, b
Zickfeld, K., Eby, M., and Weaver, A. J.: Carbon‐cycle feedbacks of changes in the Atlantic meridional overturning circulation under future atmospheric CO2, Global Biogeochem. Cy., 22, GB3024, https://doi.org/10.1029/2007GB003118, 2008. a, b, c, d
Ziemen, F. A., Kapsch, M.-L., Klockmann, M., and Mikolajewicz, U.: Heinrich events show two-stage climate response in transient glacial simulations, Clim. Past, 15, 153–168, https://doi.org/10.5194/cp-15-153-2019, 2019. a