Subglacial hydrology plays a key role in many glaciological processes, including ice dynamics via the modulation of basal sliding. Owing to the lack of an overarching theory, however, a variety of model approximations exist to represent the subglacial drainage system. The Subglacial Hydrology Model Intercomparison Project (SHMIP) provides a set of synthetic experiments to compare existing and future models. We present the results from 13 participating models with a focus on effective pressure and discharge. For many applications (e.g. steady states and annual variations, low input scenarios) a simple model, such as an inefficient-system-only model, a flowline or lumped model, or a porous-layer model provides results comparable to those of more complex models. However, when studying short term (e.g. diurnal) variations of the water pressure, the use of a two-dimensional model incorporating physical representations of both efficient and inefficient drainage systems yields results that are significantly different from those of simpler models and should be preferentially applied. The results also emphasise the role of water storage in the response of water pressure to transient recharge. Finally, we find that the localisation of moulins has a limited impact except in regions of sparse moulin density.

Abstract. Subglacial hydrology has a strong influence on glacier and ice sheet dynamics, particularly through the dependence of sliding velocity on subglacial water pressure. Significant challenges are involved in modeling subglacial hydrology, as the drainage geometry and flow mechanics are constantly changing, with complex feedbacks that play out between water and ice. A clear tradition has been established in the subglacial hydrology modeling literature of distinguishing between channelized (efficient) and sheetlike (inefficient or distributed) drainage systems or components and using slightly different forms of the governing equations in each subsystem to represent the dominant physics. Specifically, many previous subglacial hydrology models disregard opening by melt in the sheetlike system or redistribute it to adjacent channel elements in order to avoid runaway growth that occurs when it is included in the sheetlike system. We present a new subglacial hydrology model, SHAKTI (Subglacial Hydrology and Kinetic, Transient Interactions), in which a single set of governing equations is used everywhere, including opening by melt in the entire domain. SHAKTI employs a generalized relationship between the subglacial water flux and the hydraulic gradient that allows for the representation of laminar, turbulent, and transitional regimes depending on the local Reynolds number. This formulation allows for the coexistence of these flow regimes in different regions, and the configuration and geometry of the subglacial system evolves naturally to represent sheetlike drainage as well as systematic channelized drainage under appropriate conditions. We present steady and transient example simulations to illustrate the features and capabilities of the model and to examine sensitivity to mesh size and time step size. The model is implemented as part of the Ice Sheet System Model (ISSM).

We explore potential changes in Greenland ice sheet form and flow associated with increasing ice temperatures and relaxing effective ice viscosities. We define “thermal‐viscous collapse” as a transition from the polythermal ice sheet temperature distribution characteristic of the Holocene to temperate ice at the pressure melting point and associated lower viscosities. The conceptual model of thermal‐viscous collapse we present is dependent on: (1) sufficient energy available in future meltwater runoff, (2) routing of meltwater to the bed of the ice sheet interior, and (3) efficient energy transfer from meltwater to the ice. Although we do not attempt to constrain the probability of thermal‐viscous collapse, it appears thermodynamically plausible to warm the deepest 15% of the ice sheet, where the majority of deformational shear occurs, to the pressure melting point within four centuries. First‐order numerical modeling of an end‐member scenario, in which prescribed ice temperatures are warmed at an imposed rate of 0.05 K/a, infers a decrease in ice sheet volume of 5 ± 2% within five centuries of initiating collapse. This is equivalent to a cumulative sea‐level rise contribution of 33 ± 18 cm. The vast majority of the sea‐level rise contribution associated with thermal‐viscous collapse, however, would likely be realized over subsequent millennia.

No abstract

Results are presented and compared for the Community Earth System Model version 2 (CESM2) simulations of the middle Holocene (MH, 6 ka) and Last Interglacial (LIG, 127 ka). These simulations are designated as Tier 1 experiments (midHolocene and lig127k) for the Coupled Model Intercomparison Project phase 6 (CMIP6) and the Paleoclimate Modeling Intercomparison Project phase 4 (PMIP4). They use the low-top, standard 1°version of CESM2 contributing to CMIP6 DECK, historical, and future projection simulations, and to other modeling intercomparison projects. The midHolocene and lig127k provide the opportunity to examine the responses in CESM2 to the orbitally induced changes in the seasonal and latitudinal distribution of insolation. The insolation anomalies result in summer warming over the Northern Hemisphere continents, reduced Arctic summer minimum sea ice, and increased areal extent of the North African monsoon. The Arctic remains warm throughout the year. These changes are greater in the lig127k than midHolocene simulation. Other notable changes are reduction of the Niño3.4 variability and Drake Passage transport and a small increase in the Atlantic Meridional Overturning Circulation from the piControl to midHolocene to lig127k simulation. Comparisons to paleo-data and to simulations from previous model versions are discussed. Possible reasons for mismatches with the paleo-observations are proposed, including missing processes in CESM2, simplifications in the CMIP6 protocols for these experiments, and dating and calibration uncertainties in the data reconstructions. Plain Language Summary The modeling of past climates, using physically based tools and comparing to paleo-observations, is increasingly seen as a strong test of the models that are used for the projection of future climate changes. We report on simulations by one of the latest large-scale climate models, the Community Earth System Model version 2 (CESM2), for the two most recent warm epochs: the current interglacial-the Holocene, which began 11,650 years agoand the previous interglacial-the Last Interglacial, which extended from 129,000 to 116,000 years ago. The simulations find summer warming over the Northern Hemisphere continents, reduced Arctic summer sea ice, and increased areal extent of the North African monsoon. Cryospheric and oceanic feedbacks drive year-round Arctic warmth. The CESM2 middle Holocene and Last Interglacial simulations are an important resource for the community participating in the Coupled Model Intercomparison Project phase 6 (CMIP6).