Ice Sheet Interactions Research Articles

New estimation of critical insolation–CO2 relationship for triggering glacial inception

Abstract. It has been previously proposed that glacial inception represents a bifurcation transition between interglacial and glacial states and is governed by the nonlinear dynamics of the climate–cryosphere system. To trigger glacial inception, the orbital forcing (defined as the maximum of summer insolation at 65° N and determined by Earth’s orbital parameters) must be lower than a critical level, which depends on the atmospheric CO2 concentration. While paleoclimatic data do not provide a strong constraint on the dependence between CO2 and critical insolation, its accurate estimation is of fundamental importance for predicting future glaciations and the effect that anthropogenic CO2 emissions might have on them. In this study, we use the novel Earth system model of intermediate complexity CLIMBER-X with interactive ice sheets to produce a new estimation of the critical insolation–CO2 relationship for triggering glacial inception. We perform a series of experiments in which different combinations of orbital forcing and atmospheric CO2 concentration are maintained constant in time. We analyze for which combinations of orbital forcing and CO2 glacial inception occurs and trace the critical relationship between them, separating conditions under which glacial inception is possible from those where glacial inception is not materialized. We also provide a theoretical foundation for the proposed critical insolation–CO2 relation. We find that the use of the maximum summer insolation at 65° N as a single metric for orbital forcing is adequate for tracing the glacial inception bifurcation. Moreover, we find that the temporal and spatial patterns of ice sheet growth during glacial inception are not always the same but depend on the critical insolation and CO2 level. The experiments evidence the fact that during glacial inception, ice sheets grow mostly in North America, and only under low CO2 conditions are ice sheets also formed over Scandinavia. The latter is associated with a weak Atlantic Meridional Overturning Circulation (AMOC) for low CO2. We find that the strength of AMOC also affects the rate of ice sheet growth during glacial inception.

Rapid Laurentide Ice Sheet growth preceding the Last Glacial Maximum due to summer snowfall

There has been extensive research into the nonlinear responses of the Earth system to astronomical forcing during the last glacial cycle. However, the speed and spatial geometry of ice sheet expansion to its largest extent at the Last Glacial Maximum 21 thousand years ago remains uncertain. Here we use an Earth system model with interactive ice sheets to show that distinct initial North American (Laurentide) ice sheets at 38 thousand years ago converge towards a configuration consistent with the Last Glacial Maximum due to feedbacks between atmospheric circulation and ice sheet geometry. Notably, ice advance speed and spatial pattern in our model are controlled by the amount of summer snowfall, which is dependent on moisture transport pathways from the North Atlantic warm pool linked to ice sheet geometry. The consequence of increased summer snowfall on the surface mass balance of the ice sheet is not only the direct increase in accumulation but the indirect reduction in melt through the snow/ice–albedo feedback. These feedbacks provide an effective mechanism for ice growth for a range of initial ice sheet states and may explain the rapid North American ice volume increase during the last ice age and potentially driving growth during previous glacial periods.

Glacial inception through rapid ice area increase driven by albedo and vegetation feedbacks

Abstract. We present transient simulations of the last glacial inception using the Earth system model CLIMBER-X with dynamic vegetation, interactive ice sheets, and visco-elastic solid Earth responses. The simulations are initialized at the middle of the Eemian interglacial (125 kiloyears before present, ka) and run until 100 ka, driven by prescribed changes in Earth's orbital parameters and greenhouse gas concentrations from ice core data. CLIMBER-X simulates a rapid increase in Northern Hemisphere ice sheet area through MIS5d, with ice sheets expanding over northern North America and Scandinavia, in broad agreement with proxy reconstructions. While most of the increase in ice sheet area occurs over a relatively short period between 119 and 117 ka, the larger part of the increase in ice volume occurs afterwards with an almost constant ice sheet extent. We show that the vegetation feedback plays a fundamental role in controlling the ice sheet expansion during the last glacial inception. In particular, with prescribed present-day vegetation the model simulates a global sea level drop of only ∼ 20 m, compared with the ∼ 35 m decrease in sea level with dynamic vegetation response. The ice sheet and carbon cycle feedbacks play only a minor role during the ice sheet expansion phase prior to ∼ 115 ka but are important in limiting the deglaciation during the following phase characterized by increasing summer insolation. The model results are sensitive to climate model biases and to the parameterization of snow albedo, while they show only a weak dependence on changes in the ice sheet model resolution and the acceleration factor used to speed up the climate component. Overall, our simulations confirm and refine previous results showing that climate–vegetation–cryosphere feedbacks play a fundamental role in the transition from interglacial to glacial states characterizing Quaternary glacial cycles.

A New Eulerian Iceberg Module for Climate Studies

AbstractIcebergs modulate the effective location of freshwater input from ice sheets into the ocean and therefore play an important role for the climate, especially during times of increased ice discharge (e.g., Heinrich events). None of the models participating in the Paleo Modeling Intercomparison Project simulations of the Last Glacial Maximum or the last deglaciation included icebergs. Here, we present a newly developed dynamic/thermodynamic iceberg module that was specifically designed to be incorporated in climate models used for long‐term climate simulations with interactive ice sheets. In contrast to the widely used Lagrangian iceberg models, it is formulated in an Eulerian framework. This simplifies coupling to ocean models and enhances computational efficiency for glacial climates. In a set of sensitivity experiments, where the module was implemented into an Earth System Model, we validate the model for present‐day climate conditions and test its sensitivity to key parameters. Further, we investigate the effect of iceberg hosing on the Atlantic meridional overturning circulation (AMOC) as compared to traditional freshwater hosing. Varying the hosing rate slowly in time yields a good approximation of the hysteresis curve of the AMOC. We find that the sensitivity of the AMOC to iceberg hosing is stronger than to freshwater hosing in the same ocean point, but weaker as compared to a latitude belt forcing in the North Atlantic. This emphasizes the 2e necessity to include interactive icebergs in long‐term coupled climate simulations to realistically represent melt patterns and the response of the AMOC to freshwater input from melting ice sheets.

Surface mass balance and climate of the Last Glacial Maximum Northern Hemisphere ice sheets: simulations with CESM2.1

Abstract. The Last Glacial Maximum (LGM, from ∼26 to 20 ka BP) was the most recent period with large ice sheets in Eurasia and North America. At that time, global temperatures were 5–7 ∘C lower than today, and sea level ∼125 m lower. LGM simulations are useful to understand earth system dynamics, including climate–ice sheet interactions, and to evaluate and improve the models representing those dynamics. Here, we present two simulations of the Northern Hemisphere ice sheet climate and surface mass balance (SMB) with the Community Earth System Model v2.1 (CESM2.1) using the Community Atmosphere Model v5 (CAM5) with prescribed ice sheets for two time periods that bracket the LGM period: 26 and 21 ka BP. CESM2.1 includes an explicit simulation of snow/firn compaction, albedo, refreezing, and direct coupling of the ice sheet surface energy fluxes with the atmosphere. The simulated mean snow accumulation is lowest for the Greenland and Barents–Kara Sea ice sheets (GrIS, BKIS) and highest for British and Irish (BIIS) and Icelandic (IcIS) ice sheets. Melt rates are negligible for the dry BKIS and GrIS, and relatively large for the BIIS, North American ice sheet complex (NAISC; i.e. Laurentide, Cordilleran, and Innuitian), Scandinavian ice sheet (SIS), and IcIS, and are reduced by almost a third in the colder (lower temperature) 26 ka BP climate compared with 21 ka BP. The SMB is positive for the GrIS, BKIS, SIS, and IcIS during the LGM (26 and 21 ka BP) and negative for the NAISC and BIIS. Relatively wide ablation areas are simulated along the southern (terrestrial), Pacific and Atlantic margins of the NAISC, across the majority of the BIIS, and along the terrestrial southern margin of the SIS. The integrated SMB substantially increases for the NAISC and BIIS in the 26 ka BP climate, but it does not reverse the negative sign. Summer incoming surface solar radiation is largest over the high interior of the NAISC and GrIS, and minimum over the BIIS and southern margin of NAISC. Summer net radiation is maximum over the ablation areas and minimum where the albedo is highest, namely in the interior of the GrIS, northern NAISC, and all of the BKIS. Summer sensible and latent heat fluxes are highest over the ablation areas, positively contributing to melt energy. Refreezing is largest along the equilibrium line altitude for all ice sheets and prevents 40 %–50 % of meltwater entering the ocean. The large simulated melt for the NAISC suggests potential biases in the climate simulation, ice sheet reconstruction, and/or highly non-equilibrated climate and ice sheet at the LGM time.

The Subglacial Lake That Wasn't There: Improved Interpretation From Seismic Data Reveals a Sediment Bedform at Isunnguata Sermia

AbstractRadio Echo Sounding (RES) surveys conducted in May 2010 and April 2011 revealed a 2 km2 flat area with increased bed reflectivity at the base of Isunnguata Sermia at the western margin of the Greenland Ice Sheet. This flat reflector was located within a localized subglacial hydraulic potential (hydropotential) minimum, as part of a complex and elongated trough system. By analogy with comparable features in Antarctica, the initial interpretation of such a feature was a potential subglacial lake. In September 2013 a co‐located seismic survey revealed a 1,750 m by 540 and 37 m thick stratified lens‐shaped bedform at the base of a subglacial trough system. Amplitude Versus Angle (AVA) analysis yields a derived reflection coefficient R = 0.09 ± 0.14 indicative of consolidated sediments possibly overlain by dilatant till. The bed and flank on the northern side of the trough consist of unconsolidated, possibly water‐bearing sediments with R = −0.10 ± 0.08, whereas on the southern side it consists of more consolidated material. We interpret the trough as a key component of the wider subglacial drainage network, for which the sediments on its northern side act as a localized water‐storage reservoir. Given the observation of seasonally forming and rapidly draining supraglacial meltwater lakes in this area, we interpret the lens‐shaped bedform as deposited by episodically ponding meltwater within the subglacial trough system. Our results highlight the importance of transient subglacial hydrological and sedimentological processes such as drainage events for the interaction of ice sheets and their substrates, to understand ice dynamics in a warming climate.

The Antarctic contribution to 21st-century sea-level rise predicted by the UK Earth System Model with an interactive ice sheet

Abstract. The Antarctic Ice Sheet will play a crucial role in the evolution of global mean sea level as the climate warms. An interactively coupled climate and ice sheet model is needed to understand the impacts of ice–climate feedbacks during this evolution. Here we use a two-way coupling between the UK Earth System Model and the BISICLES (Berkeley Ice Sheet Initiative for Climate at Extreme Scales) dynamic ice sheet model to investigate Antarctic ice–climate interactions under two climate change scenarios. We perform ensembles of SSP1–1.9 and SSP5–8.5 (Shared Socioeconomic Pathway) scenario simulations to 2100, which we believe are the first such simulations with a climate model that include two-way coupling of atmosphere and ocean models to dynamic models of the Greenland and Antarctic ice sheets. We focus our analysis on the latter. In SSP1–1.9 simulations, ice shelf basal melting and grounded ice mass loss from the Antarctic Ice Sheet are generally lower than present rates during the entire simulation period. In contrast, the responses to SSP5–8.5 forcing are strong. By the end of the 21st century, these simulations feature order-of-magnitude increases in basal melting of the Ross and Filchner–Ronne ice shelves, caused by intrusions of masses of warm ocean water. Due to the slow response of ice sheet drawdown, this strong melting does not cause a substantial increase in ice discharge during the simulations. The surface mass balance in SSP5–8.5 simulations shows a pattern of strong decrease on ice shelves, caused by increased melting, and strong increase on grounded ice, caused by increased snowfall. Despite strong surface and basal melting of the ice shelves, increased snowfall dominates the mass budget of the grounded ice, leading to an ensemble mean Antarctic contribution to global mean sea level of a fall of 22 mm by 2100 in the SSP5–8.5 scenario. We hypothesise that this signal would revert to sea-level rise on longer timescales, caused by the ice sheet dynamic response to ice shelf thinning. These results demonstrate the need for fully coupled ice–climate models in reducing the substantial uncertainty in sea-level rise from the Antarctic Ice Sheet.

Importance of the Antarctic Slope Current in the Southern Ocean Response to Ice Sheet Melt and Wind Stress Change

AbstractWe use two coupled climate models, GFDL‐CM4 and GFDL‐ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre‐industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near‐surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large‐scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean‐state representation and meltwater‐driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.

Distribution Analysis of Local Ice Pressures in the Indentation Test at Various Velocities

When sea ice acts on vertical structures, there are much higher pressures in localized areas known as high-pressure zones (HPZs) than in other areas. The damage failure mode of sea ice varies with the sea ice velocity and affects the distribution of HPZs. In this study, an indentation test that drives ice sheet interaction with a vertical rigid plate (indentor) was designed, and a pressure sensor (consisting of 32 × 32 small pressure units of 100 mm2) was installed on the indentor face to measure the local ice pressure (LIP) at various velocities. The distribution of the LIPs along the ice thickness, the probability distribution of the LIPs and the distribution relationship of the LIPS in space and time were obtained from the measurement. The results show that the LIPs were mainly distributed in the middle of the sea ice, which is consistent with full-scale observations and previous research. Ductile failure of the sea ice results in a larger LIP distribution area than brittle failure at the same threshold kt (kt = σL/σcr). The probability distribution of the LIPs decreases exponentially with increasing pressure and follows a lognormal distribution. The maximum LIP appears at the peak moment of the global force when the sea ice failure mode is mainly ductile failure. However, the maximum LIP may not occur at the peak moment when the sea ice failure is mainly brittle failure and, instead can appear at any moment in the global force time history curve. The HPZ (which is larger than 7/8 times the maximum LIP) area is less than 2% of the nominal contact area at various velocities. The influence of the sea ice velocity on the spatial and temporal distribution of LIP is analyzed, and the results provide a reference for designing structures with local strength in ice regions.

The role of an interactive Greenland ice sheet in the coupled climate-ice sheet model EC-Earth-PISM

Ice sheet processes are often simplified in global climate models as changes in ice sheets have been assumed to occur over long time scales compared to ocean and atmospheric changes. However, numerous observations show an increasing rate of mass loss from the Greenland Ice Sheet and call for comprehensive process-based models to explore its role in climate change. Here, we present a new model system, EC-Earth-PISM, that includes an interactive Greenland Ice Sheet. The model is based on the EC-Earth v2.3 global climate model in which ice sheet surface processes are introduced. This model interacts with the Parallel Ice Sheet Model (PISM) without anomaly or flux corrections. Under pre-industrial climate conditions, the modeled climate and ice sheet are stable while keeping a realistic interannual variability. In model simulations forced into a warmer climate of four times the pre-industrial CO2 concentration, the total surface mass balance decreases and the ice sheet loses mass at a rate of about 500 Gt/year. In the climate warming experiments, the resulting freshwater flux from the Greenland Ice Sheet increases 55% more in the experiments with the interactive ice sheet and the climate response is significantly different: the Arctic near-surface air temperature is lower, substantially more winter sea ice covers the northern hemisphere, and the ocean circulation is weaker. Our results indicate that the melt-albedo feedback plays a key role for the response of the ice sheet and its influence on the changing climate in the Arctic. This emphasizes the importance of including interactive ice sheets in climate change projections.

PARASO, a circum-Antarctic fully coupled ice-sheet–ocean–sea-ice–atmosphere–land model involving f.ETISh1.7, NEMO3.6, LIM3.6, COSMO5.0 and CLM4.5

Abstract. We introduce PARASO, a novel five-component fully coupled regional climate model over an Antarctic circumpolar domain covering the full Southern Ocean. The state-of-the-art models used are the fast Elementary Thermomechanical Ice Sheet model (f.ETISh) v1.7 (ice sheet), the Nucleus for European Modelling of the Ocean (NEMO) v3.6 (ocean), the Louvain-la-Neuve sea-ice model (LIM) v3.6 (sea ice), the COnsortium for Small-scale MOdeling (COSMO) model v5.0 (atmosphere) and its CLimate Mode (CLM) v4.5 (land), which are here run at a horizontal resolution close to 1/4∘. One key feature of this tool resides in a novel two-way coupling interface for representing ocean–ice-sheet interactions, through explicitly resolved ice-shelf cavities. The impact of atmospheric processes on the Antarctic ice sheet is also conveyed through computed COSMO-CLM–f.ETISh surface mass exchange. In this technical paper, we briefly introduce each model's configuration and document the developments that were carried out in order to establish PARASO. The new offline-based NEMO–f.ETISh coupling interface is thoroughly described. Our developments also include a new surface tiling approach to combine open-ocean and sea-ice-covered cells within COSMO, which was required to make this model relevant in the context of coupled simulations in polar regions. We present results from a 2000–2001 coupled 2-year experiment. PARASO is numerically stable and fully operational. The 2-year simulation conducted without fine tuning of the model reproduced the main expected features, although remaining systematic biases provide perspectives for further adjustment and development.

Glacitectonic evidence of ice sheet interaction and retreat across the western part of Dogger Bank (North Sea) during the Last Glaciation

High-resolution 2D seismic data from the western side of Dogger Bank (North Sea) has revealed that the glacigenic sediments of the Dogger Bank Formation record a complex history of sedimentation and penecontemporaneous, large-scale, ice-marginal to proglacial glacitectonism. The resulting complex assemblage of glacial landforms and sediments record the interplay between two separate ice masses revealing that Late Devensian ice sheet dynamics across Dogger Bank were far more complex than previously thought, involving the North Sea lobe of the British and Irish Ice Sheet, advancing from the west, interacting with the Dogger Bank lobe which expanded from the north. The active northward retreat of the Dogger Bank lobe resulted in the development of a complex assemblage of arcuate thrust-block moraines (≤15 km wide, >30 km long) composed of highly folded and thrust sediments, separated by sedimentary basins and meltwater channels filled by outwash. The impact of the North Sea lobe was restricted to the western margin of Dogger Bank and led to deep-seated (100–150 m thick) glacitectonism in response to ice-push from the west. During the earlier expansion of the North Sea lobe, this thrust and fold complex initially occupied a frontal marginal position changing to a more lateral ice-marginal position as the ice sheet continued to expand to the south. The complex structural relationships between the two glacitectonic complexes indicates that these ice masses interacted along the western side of Dogger Bank, with the inundation of this area by ice probably occurring during the last glaciation when the ice sheets attained their maximum extents.

Eighteen-year record of circum-Antarctic landfast-sea-ice distribution allows detailed baseline characterisation and reveals trends and variability

Abstract. Landfast sea ice (fast ice) is an important though poorly understood component of the cryosphere on the Antarctic continental shelf, where it plays a key role in atmosphere–ocean–ice-sheet interaction and coupled ecological and biogeochemical processes. Here, we present a first in-depth baseline analysis of variability and change in circum-Antarctic fast-ice distribution (including its relationship to bathymetry), based on a new high-resolution satellite-derived time series for the period 2000 to 2018. This reveals (a) an overall trend of -882±824 km2 yr−1 (-0.19±0.18 % yr−1) and (b) eight distinct regions in terms of fast-ice coverage and modes of formation. Of these, four exhibit positive trends over the 18-year period and four negative. Positive trends are seen in East Antarctica and in the Bellingshausen Sea, with this region claiming the largest positive trend of +1198±359 km2 yr−1 (+1.10±0.35 % yr−1). The four negative trends predominantly occur in West Antarctica, with the largest negative trend of -1206±277 km2 yr−1 (-1.78±0.41 % yr−1) occurring in the Victoria and Oates Land region in the western Ross Sea. All trends are significant. This new baseline analysis represents a significant advance in our knowledge of the current state of both the global cryosphere and the complex Antarctic coastal system, which are vulnerable to climate variability and change. It will also inform a wide range of other studies.

Climate and ice sheet evolutions from the last glacial maximum to the pre-industrial period with an ice-sheet–climate coupled model

Abstract. The last deglaciation offers an unique opportunity to understand the climate–ice-sheet interactions in a global warming context. In this paper, to tackle this question, we use an Earth system model of intermediate complexity coupled to an ice sheet model covering the Northern Hemisphere to simulate the last deglaciation and the Holocene (26–0 ka). We use a synchronous coupling every year between the ice sheet and the rest of the climate system and we ensure a closed water cycle considering the release of freshwater flux to the ocean due to ice sheet melting. Our reference experiment displays a gradual warming in response to the forcings, with no abrupt changes. In this case, while the amplitude of the freshwater flux to the ocean induced by ice sheet retreat is realistic, it is sufficient to shut down the Atlantic meridional overturning circulation from which the model does not recover within the time period simulated. However, with reduced freshwater flux we are nonetheless able to obtain different oceanic circulation evolutions, including some abrupt transitions between shut-down and active circulation states in the course of the deglaciation. The inclusion of a parameterisation for the sinking of brines around Antarctica also produces an abrupt recovery of the Atlantic meridional overturning circulation, absent in the reference experiment. The fast oceanic circulation recoveries lead to abrupt warming phases in Greenland. Our simulated ice sheet geometry evolution is in overall good agreement with available global reconstructions, even though the abrupt sea level rise at 14.6 ka is underestimated, possibly because the climate model underestimates the millennial-scale temperature variability. In the course of the deglaciation, large-scale grounding line instabilities are simulated both for the Eurasian and North American ice sheets. The first instability occurs in the Barents–Kara seas for the Eurasian ice sheet at 14.5 ka. A second grounding line instability occurs ca. 12 ka in the proglacial lake that formed at the southern margin of the North American ice sheet. With additional asynchronously coupled experiments, we assess the sensitivity of our results to different ice sheet model choices related to surface and sub-shelf mass balance, ice deformation and grounding line representation. While the ice sheet evolutions differ within this ensemble, the global climate trajectory is only weakly affected by these choices. In our experiments, only the abrupt shifts in the oceanic circulation due to freshwater fluxes are able to produce some millennial-scale variability since no self-generating abrupt transitions are simulated without these fluxes.

Review of research on ice loading of offshore structures in level ice fields

The physical and mechanical properties of sea ice and the failure modes of ice sheets significantly affect the ice loading of offshore structures. This paper first introduces the physical and mechanical characteristics of sea ice, then describes the physical processes of the interaction of ice sheets with offshore structures. Finally, the progress of research on the ice loading of offshore structures is briefly analyzed, including the three main approaches of formula estimation, experimental testing and numerical simulation. By discussing the applicability of formula estimation, the research focus and difficulties of experimental testing, and the development trends of numerical simulation, we find that the model testing of ice loading is an effective method. This paper covers the research trends and problems faced by ice loading on offshore structures, and it can provide certain guiding references the related research of scholars in the field.

Antarctic Bottom Water and North Atlantic Deep Water in CMIP6 models

Abstract. Deep and bottom water formation are crucial components of the global ocean circulation, yet they were poorly represented in the previous generation of climate models. We here quantify biases in Antarctic Bottom Water (AABW) and North Atlantic Deep Water (NADW) formation, properties, transport, and global extent in 35 climate models that participated in the latest Climate Model Intercomparison Project (CMIP6). Several CMIP6 models are correctly forming AABW via shelf processes, but 28 models in the Southern Ocean and all 35 models in the North Atlantic form deep and bottom water via open-ocean deep convection too deeply, too often, and/or over too large an area. Models that convect the least form the most accurate AABW but the least accurate NADW. The four CESM2 models with their overflow parameterisation are among the most accurate models. In the Atlantic, the colder the AABW, the stronger the abyssal overturning at 30∘ S, and the further north the AABW layer extends. The saltier the NADW, the stronger the Atlantic Meridional Overturning Circulation (AMOC), and the further south the NADW layer extends. In the Indian and Pacific oceans in contrast, the fresher models are the ones which extend the furthest regardless of the strength of their abyssal overturning, most likely because they are also the models with the weakest fronts in the Antarctic Circumpolar Current. There are clear improvements since CMIP5: several CMIP6 models correctly represent or parameterise Antarctic shelf processes, fewer models exhibit Southern Ocean deep convection, more models convect at the right location in the Labrador Sea, bottom density biases are reduced, and abyssal overturning is more realistic. However, more improvements are required, e.g. by generalising the use of overflow parameterisations or by coupling to interactive ice sheet models, before deep and bottom water formation, and hence heat and carbon storage, are represented accurately.

Accelerated Greenland Ice Sheet Mass Loss Under High Greenhouse Gas Forcing as Simulated by the Coupled CESM2.1‐CISM2.1

AbstractThe Greenland ice sheet (GrIS) is now losing mass at a rate of 0.7 mm of sea level rise (SLR) per year. Here we explore future GrIS evolution and interactions with global and regional climate under high greenhouse gas forcing with the Community Earth System Model version 2.1 (CESM2.1), which includes an interactive ice sheet component (the Community Ice Sheet Model v2.1 [CISM2.1]) and an advanced energy balance‐based calculation of surface melt. We run an idealized 350‐year scenario in which atmospheric CO2 concentration increases by 1% annually until reaching four times pre‐industrial values at year 140, after which it is held fixed. The global mean temperature increases by 5.2 and 8.5 K by years 131–150 and 331–350, respectively. The projected GrIS contribution to global mean SLR is 107 mm by year 150 and 1,140 mm by year 350. The rate of SLR increases from 2 mm yr−1 at year 150 to almost 7 mm yr−1 by year 350. The accelerated mass loss is caused by rapidly increasing surface melt as the ablation area expands, with associated albedo feedback and increased sensible and latent heat fluxes. This acceleration occurs for a global warming of approximately 4.2 K with respect to pre‐industrial and is in part explained by the quasi‐parabolic shape of the ice sheet, which favors rapid expansion of the ablation area as it approaches the interior “plateau.”

An Efficient Ice Sheet/Earth System Model Spin-up Procedure for CESM2-CISM2: Description, Evaluation, and Broader Applicability.

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin‐up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet‐enabled simulations with CESM2‐CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).

Ice-resistant performance analysis of the bucket foundation offshore platform with concrete pillar

The steel bucket foundation platform is an important engineering design for small-scale oil/gas field developments due to its economic and structural characteristics. Severe ice-induced vibrations were observed during an ice sheet interaction with the fixed ice-breaking cone of JZ9-3WHPB platform in the Bohai Sea, China. The steel pillar was replaced by a concrete pillar in order to reduce the ice-induced vibrations for JZ9-3WHPB platform. A good comparison was obtained between the calculated and full-scale observation data; the numerical results show that ice-induced vibrations could be reduced effectively for JZ9-3WHPB platform with concrete pillar.

Equilibrium simulations of Marine Isotope Stage 3 climate

Abstract. An equilibrium simulation of Marine Isotope Stage 3 (MIS3) climate with boundary conditions characteristic of Greenland Interstadial 8 (GI-8; 38 kyr BP) is carried out with the Norwegian Earth System Model (NorESM). A computationally efficient configuration of the model enables long integrations at relatively high resolution, with the simulations reaching a quasi-equilibrium state after 2500 years. We assess the characteristics of the simulated large-scale atmosphere and ocean circulation, precipitation, ocean hydrography, sea ice distribution, and internal variability. The simulated MIS3 interstadial near-surface air temperature is 2.9 ∘C cooler than the pre-industrial (PI). The Atlantic meridional overturning circulation (AMOC) is deeper and intensified by ∼13 %. There is a decrease in the volume of Antarctic Bottom Water (AABW) reaching the Atlantic. At the same time, there is an increase in ventilation of the Southern Ocean, associated with a significant expansion of Antarctic sea ice and concomitant intensified brine rejection, invigorating ocean convection. In the central Arctic, sea ice is ∼2 m thicker, with an expansion of sea ice in the Nordic Seas during winter. Attempts at triggering a non-linear transition to a cold stadial climate state, by varying atmospheric CO2 concentrations and Laurentide Ice Sheet height, suggest that the simulated MIS3 interstadial state in the NorESM is relatively stable, thus underscoring the role of model dependency, and questioning the existence of unforced abrupt transitions in Greenland climate in the absence of interactive ice sheet–meltwater dynamics.