Ocean-driven basal melting of Antarctic ice shelves plays an important role in the mass loss of the Antarctic Ice Sheet. Ice shelf cavity-resolving ocean models are a valuable tool for understanding ice shelf-ocean interactions and for simulating projections of ice shelf and ocean states under future climate. Designed to assess the current state of ice shelf-ocean modelling, the second Ice Shelf-Ocean Model Intercomparison Project, ISOMIP+, consists of 12 ocean model configurations submitted with a common, idealised experimental setup. Here, we focus on the experiments Ocean0-2 , which are ocean models with idealised, static ice shelf geometries, but where the ocean reaches a balance with prescribed far-field ocean conditions. Different thermal transfer coefficient values (ranging from 0.011 to 0.2) are used for each model in the melting parameterisation to achieve a common, tuned melt rate since the models cover a range of types of vertical coordinates, ice-ocean boundary layer treatments, and numerical schemes. These model differences lead to spread in the resultant ocean properties, circulation, boundary-layer structure and spatial distribution of melting. We also highlight similarities between models, such as a shared linear relationship across most models between melt rate and overturning and barotropic streamfunctions during the spin-up and spin-down, demonstrating a robust relationship between melt and circulation across models and forcing conditions. The ISOMIP+ results provide a systematic comparison of ice shelf cavity-capable ocean models. However, we also demonstrate the need for realistic ice shelf-ocean model intercomparison projects (some already underway) to assess model biases and inter-model variation against sparse observations. Further research is needed to understand the differences between models and further improve our modelled representations of the ice-ocean boundary layer and ice shelf cavity circulation.

We have conducted an ensemble of idealised climate overshoot simulations in a state-of-the-art Earth System Model in which global mean temperature is increased at a constant rate to global warming levels (GWLs) from 1.5 to 5 degrees C by prescribing constant CO2 emissions, followed by a 50 year period of zero CO2 emissions started at different GWLs, and then a period of cooling in which CO2 is removed from the atmosphere. We give an overview of the ice sheet, sea-ice and sea level responses in these simulations highlighting long term responses at different GWLs and discuss the degree to which those responses can be simply reversed as global surface temperature cools. We show a broad divide between the two hemispheres, in which Northern Hemisphere polar processes have larger direct responses to warming which are more simply reversible than those in the Southern Hemisphere.Cessation of CO2 emissions at most GWLs stabilises surface temperatures at high northern latitudes, although a slow warming trend continues at high southern latitudes. Northern Hemisphere sea-ice extent and Greenland surface mass balance both stabilise under zero CO2 emissions and appear to return toward preindustrial levels along with surface temperatures when CO2 is sequestered, but Southern Hemisphere sea-ice continues to decline under zero emissions and does not simply increase again as global temperatures cool. Likewise, the Antarctic circumpolar westerlies exhibit strengthening and shift poleward under global warming, but do not simply return to their preindustrial state when the climate is cooled, with implications for ocean circulation and Antarctic surface mass balance. The thermosteric contributions to global mean sea level from ocean heat uptake and the Greenland ice sheet continue when CO2 emissions are ceased or reversed, at rates which slowly decline on centennial timescales. The net mass balance of Antarctica and its contribution to sea level do not simply scale with global temperature; they result from a complex interaction between the basal and surface mass balances, the ice-dynamic response to those forcings and significant trends inherent to the initialisation of our ice sheet model state.

The configuration of the Northern Hemisphere ice sheets during the Last and Penultimate Glacial Maxima (LGM; PGM) influenced millennial-scale climate changes, as well as solid Earth and sea level changes that occurred during the subsequent deglaciations, due to their effects on the atmosphere, ocean circulation, and the solid Earth. Thus, realistic simulations of these ice sheets are crucial for the initialisation of deglaciation experiments that can help improve our understanding of interactions between the climate, ice sheets and sea levels. Here, we produce the first large ensembles of complex coupled atmosphere-ice sheet model (FAMOUS-BISICLES) simulations of the PGM and LGM, varying 12 uncertain parameters that control the ice sheet albedo, ice dynamics and climate. We quantify the sensitivity to input parameters using Gaussian Process emulators to perform a Sobol sensitivity analysis. Albedo parameters have the largest influence on ice volumes for both ice sheets and time periods. Parameters controlling precipitation and sliding have a larger effect on Eurasian than North American ice sheet size due to the differences in geographical and climatic settings. Out of 120 parameter combinations, we find 4 that produce LGM and PGM ice volumes and extents compatible with palaeo-evidence. The resulting ice sheet configurations provide new and improved reconstructions of PGM Northern Hemisphere ice sheets for use as inputs in climate, ice sheet and sea level models.

The Land and Land Ice Theme in the Coupled Model Intercomparison Project Phase 7 (CMIP7) represents the current understanding of physical processes in land surface ecosystems, hydrology, cryosphere, and their physical interactions with other Earth system components. Simulations from Earth system models (ESMs) could provide crucial information for assessing planetary safety, such as critical tipping elements, and be used to inform climate risks for improving climate impact assessments and policy decisions. This paper presents a collaborative effort to identify scientific opportunities in the Land and Land Ice Theme of the CMIP7 Data Request. The proposed opportunities build upon advances in ESMs, including new freshwater system and land ice processes being included in CMIP7, as well as the scientific community's demand for high-frequency and sub-grid-scale land surface outputs. In total, 25 variable groups that contain 716 variables have been identified to be potentially available to the broad scientific audience for performing analysis in land-atmosphere coupling, hydrological processes and freshwater systems, glacier and ice sheet mass balance and their influence on the sea levels, land use, and plant phenology. Key reflections from this data request effort include advocacy for closer engagement between the user community and modeling groups, reduction in the technical barriers to tracking existing parameters and defining new variables, and more streamlined variable management. These will be essential to enhance the usability and reliability of CMIP7 outputs for climate and Earth system research and applications to a broad audience that relies on the CMIP7 endeavor.

There is a strong concern about how fast and how much the global sea level will rise in the next few decades due to the current global warming. However, the projection range is large due to uncertainties about the future evolution of the Antarctic ice sheet, particularly the possibility of the ice shelves' collapse. With the base submerged in seawater, these ice shelves are strongly influenced by the surrounding oceanic conditions, which can be split into two regimes: cold or warm cavity. When ice shelves are exposed to warm water, basal melt increases sharply, leading to a loss of buttressing of the grounded ice upstream and potentially the collapse of the shelf. Results from TIPMIP (Tipping Points Modelling Intercomparison Project) idealised experiments carried out by the UK Earth System Model (UKESM) with an interactive Antarctic ice sheet component suggest that tipping points for several ice shelves will be reached in the future at relatively high global warming levels (GWLs). However, it is questionable whether the warming thresholds for tipping that we find are realistic due to the model biases and other uncertainties. This study focuses on exploring the uncertainty in the climate simulated by UKESM and assessing the consequences of ice shelf tipping on the wider Earth System. To do this, we induce the cavity regime shift at a lower GWL than the reported threshold by mimicking the key climate change forcing identified from the higher GWL experiments via an artificial freshwater around the Antarctic ice sheet margin. By doing so, we obtain a pair of low GWL experiments with and without ice shelf tipping, which allows us to isolate the impact of ice shelf tipping on the Earth System. In addition, the experiment setup also allows us to explore the consequences of different scenarios of various freshwater hosing values. The preliminary results indicate that the excessive freshwater induces an expansion of Southern Ocean sea ice, leading to a cooling trend in global mean temperature.

After having lived as the dominant human species in Europe for over 200 kyr, Homo neanderthalensis (the Neanderthals) disappeared around 40 kyr BP (Before Present) Higham T (2014). Competition with Homo sapiens, who arrived in Europe around the same time, is often invoked to explain this extinction. Others have argued that climate change may have reduced the living space of this population making its disappearance more rapid. In order to test the climate change hypothesis we modelled the Neanderthals' ecological niches in Europe between 90 and 50 kyr BP through paleoenvironmental reconstructions and Eco-Cultural Niche Modelling. We selected five environmental variables (orographic height, mean annual precipitation, mean temperature of the coldest month, carrying capacity and friction, see below) from climate model simulations of 5 periods between 90 and 50 kyr BP in Europe. We used Structural Similarity (SSIM) index to compare the probability maps of suitable niches to Neanderthals performed by Maxent. After a strong initial environmental change between the first (P1 = 90 to 83 kyr BP) and second (P2 = 83 to 69 kyr BP) periods, our results show that large areas highly suitable for Neanderthal occupation persisted across Europe. As our results show an increase/stability of the areas suitable to Neanderthals, the question of the cause of the decrease or displacement of the Neanderthal population towards southern Europe after this climatic change remains open.

Abstract. Simulations of past periods are useful for testing the ability of numerical models to simulate ice sheet changes under significantly different climate conditions to present day. This can help improve projections of future sea level rise made by these same models and avoid over-tuning to particular (e.g. modern) stationary climate conditions. The Last Glacial Maximum (LGM; ~21 thousand years ago (ka)) has been extensively used for this purpose since it is relatively well constrained by empirical evidence. However, less is known about the Penultimate Glacial Maximum (PGM; ~140 ka) and why the vast ice sheets covering much of the Northern Hemisphere (NH), differed to the LGM. The answer likely lies, at least in part, in the different orbital configurations between the two periods, and the resulting impact on climate-ice sheet interactions. Here, we perform and compare the first large ensembles of coupled climate-ice sheet (FAMOUS-BISICLES) simulations of the LGM and PGM to better understand how NH ice sheets interact with the climate and quantify how sensitive the simulations are to the choice of uncertain model inputs, including physical parameter values. Specifically, we vary 12 uncertain parameters that control the model representations of ice sheet albedo, ice dynamics and climate. The ensembles are evaluated against palaeo-evidence of global mean temperature, ice volume and extent to calibrate the model and find combinations of parameters that simultaneously yield plausible ice sheets and climates for both periods. The sensitivity of the North American ice sheet and the Eurasian ice sheet during the LGM and PGM, to each of the 12 parameter values, is explored using Gaussian Process emulators to perform a Sobol sensitivity analysis. From the whole ensemble, we find two simulations that meet our evaluation constraints for the LGM ice sheets. The parameter values that influence the albedo of the ice sheet have the largest influence on the resulting ice sheet volumes, but several other parameters display different sensitivity indices depending on the ice sheet (North American versus Eurasian) and time period (PGM versus LGM). This includes parameters that affect the cloud liquid water, lapse rate, basal sliding and downscaling elevation heights.

The configuration of the Northern Hemisphere ice sheets during the Penultimate Glacial Maximum differed to the Last Glacial Maximum. However, the reasons for this are not yet fully understood. These differences likely contributed to the varied deglaciation pathways experienced following the glacial maxima and may have had consequences for the interglacial sea level rise. To understand the differences between the North American Ice Sheet at the Last and Penultimate glacial maxima (21 and 140 ka), we perform two perturbed-physics ensembles of 62 simulations using a coupled atmosphere-ice sheet model, FAMOUS-ice, with prescribed surface ocean conditions, in which the North American and Greenland ice sheets are dynamically simulated with the Glimmer ice sheet model. We apply an implausibility metric to find ensemble members that match reconstructed ice extent and volumes at the Last and Penultimate glacial maxima. We use a resulting set of "plausible" parameters to perform sensitivity experiments to decompose the role of climate forcings (orbit, greenhouse gases) and initial conditions on the final ice sheet configurations. This confirms that the initial ice sheet conditions used in the model are extremely important in determining the difference in final ice volumes between both periods due to the large effect of the ice-albedo feedback. In contrast to evidence of a smaller Penultimate North American Ice Sheet, our results show that the climate boundary conditions at these glacial maxima, if considered in isolation, imply a larger Penultimate Glacial Maximum North American Ice Sheet than at the Last Glacial Maximum by around 6 m sea level equivalent. This supports the notion that the growth of the ice sheet prior to the glacial maxima is key in explaining the differences in North American ice volume.

The Last Glacial Maximum (LGM) was characterised by huge ice sheets covering the Northern Hemisphere, especially over North America, and by its cold climate. Previous authors have performed numerical simulations of the LGM to better understand coupled climate–ice sheet systems. However, the results of such simulations are sensitive to many model parameters. Here, we perform a 200-member ensemble of simulations of the North American and Greenland ice sheets and climate of the LGM with a coupled ice sheet–atmosphere–slab ocean model (FAMOUS-BISICLES) to explore sensitivities of the coupled climate–ice system to 16 uncertain parameters. In the ensemble of simulations, the global mean surface temperature is primarily controlled by the combination of parameters in the large-scale condensation scheme and the cumulus convection scheme. In simulations with plausible LGM global mean surface temperatures, we find that the albedo parameters have only a small impact on the Greenland ice volume due to the limited area of surface ablation associated with the cold climate. Instead, the basal sliding law controls the ice volume by affecting ice transport from the interior to the margin. On the other hand, like the Greenland ice sheet in future climate change, the LGM North American ice sheet volume is controlled by parameters in the snow and ice albedo scheme. Few of our simulations produce an extensive North American ice sheet when the global temperature is above 12 °C. Based on constraints on the LGM global mean surface temperature, the ice volume and the southern extent of the North American ice sheet, we select 16 acceptable simulations. These simulations lack the southern extent of ice compared to reconstructions, but they show reasonable performance on the ice sheet configuration and ice streams facing Baffin Bay and the Arctic Ocean. The strong sensitivities of the North American ice sheet to albedo at the LGM may imply a potential constraint on the future Greenland ice sheet by constraining the albedo schemes.

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from - 5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice-climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

Abstract Objectives Despite the success of immune checkpoint blockade, most metastatic melanoma patients fail to respond to therapy or experience severe toxicity. Assessment of biomarkers and immunophenotypes before or early into treatment will help to understand favourable responses and improve therapeutic outcomes. Methods We present a high‐dimensional approach for blood T‐cell profiling using three multi‐parameter cytometry panels: (1) a TruCount panel for absolute cell counts, (2) a 27‐colour spectral panel assessing T‐cell markers and (3) a 20‐colour spectral panel evaluating intracellular cytokine expression. Pre‐treatment blood mononuclear cells from patients and healthy controls were cryopreserved before staining across 11 batches. Batch effects were tracked using a single‐donor control and the suitability of normalisation was assessed. The data were analysed using manual gating and high‐dimensional strategies. Results Batch‐to‐batch variation was minimal, as demonstrated by the dimensionality reduction of batch‐control samples, and normalisation did not improve manual or high‐dimensional analysis. Application of the workflow demonstrated the capacity of the panels and showed that patients had fewer lymphocytes than controls (P = 0.0027), due to lower naive CD4+ (P = 0.015) and CD8+ (P = 0.011) T cells and follicular helper T cells (P = 0.00076). Patients showed trends for higher proportions of Ki67 and IL‐2‐expressing cells within CD4+ and CD8+ memory subsets, and increased CD57 and EOMES expression within TCRγδ+ T cells. Conclusion Our optimised high‐parameter spectral cytometry approach provided in‐depth profiling of blood T cells and found differences in patient immunophenotype at baseline. The robustness of our workflow, as demonstrated by minimal batch effects, makes this approach highly suitable for the longitudinal evaluation of immunotherapy effects.

One of the largest sources of uncertainty when projecting the Antarctic contribution to sea-level rise is the ocean-induced melt at the base of Antarctic ice shelves. This is because resolving the ocean circulation and the ice-ocean interactions occurring in the cavity below the ice shelves is computationally expensive.Instead, for large ensembles and long-term projections of the ice-sheet evolution, ice-sheet models currently rely on parameterisations to link the ocean temperature and salinity in front of ice shelves to the melt at their base. However, current physics-based parameterisations struggle to accurately simulate basal melt patterns.As an alternative approach, we explore the potential use of a deep feedforward neural network as a basal melt parameterisation. To do so, we train a neural network to emulate basal melt rates simulated by highly-resolved circum-Antarctic ocean simulations. We explore the influence of different input variables and show that the neural network struggles to generalise to ice-shelf geometries unseen during training, while it generalises better on timesteps unseen during training. We also test the parameterisation on separate coupled ocean-ice simulations to assess the neural network’s performance on independent data.

Testing the ability of climate-ice sheet coupled models to simulate past ice sheets and climates can provide a way to evaluate the models. One example is the Last Glacial Maximum (LGM), when huge ice sheets covered the Northern Hemisphere, especially over the North America. Here, we perform 200 ensemble member simulations of the North American and Greenland ice sheets and climate of the LGM with an ice sheet-atmosphere-slab ocean coupled model Famous-BISICLES. 16 parameters associated with climate and ice dynamics are varied. The simulated results are evaluated against the LGM global temperature, the total ice volume and the ice extent at the southern margin of the North American ice sheet. In the ensemble simulations, the global temperature is controlled by the combination of precipitation efficiency in the large-scale condensation and entrainment rate in the cumulus convection. Under reasonable LGM global temperature, we find that the surface albedo and Weertman coefficient in the basal sliding law control the North American ice volume. In contrast, the ice volume of Greenland is found to be controlled by the Weertman coefficient. Based on the constraints, the model produces 6 good simulations with reasonable global temperature and North American ice sheet. We also find that warm summer surface temperature biases at the ice sheet interior as well as downscaling of surface mass balance based on altitude can cause strong local ice melting. This implies the need of better representing the atmospheric conditions and surface mass balance in the ice sheet interior.

Coupled climate-ice sheet models are crucial to evaluating climate-ice feedbacks' role in future ice sheet evolution. Such models are calibrated to reproduce modern-day ice sheets, but current observations alone are insufficient to constrain the strength of climate-ice feedbacks. The extent of the Northern Hemisphere ice sheets during the last glacial maximum, ~20,000 years ago, is well known and could provide a benchmark for calibrating coupled climate-ice sheet models. We test this with the FAMOUS-ice coupled Climate-Ice Sheet model (Smith et al., 2020), a fast GCM coupled to the Glimmer ice sheet model. We ran Last Glacial Maximum simulations using FAMOUS-ice with interactive North American Ice Sheet, following the PMIP4 protocol (Kageyama et al., 2018). We find that the standard model setup, calibrated to produce a good present-day Greenland (Smith et al., 2020), produced a collapsed North American ice sheet at the Last Glacial Maximum. We ran ensembles of hundreds of simulations to explore the influence of uncertain ice sheet, albedo, atmospheric, and oceanic parameters on the ice sheet extent. The North American continent deglaciated rapidly in most of our simulations, leaving only a handful of useful simulations out of 280. We thus developed a method to efficiently identify regions of the parameter space that can produce a reasonable ice-sheet extent. This involved emulating the equilibrium ice volume and area as a function of the surface mass balance at the start of our simulations. We then ran three waves of short simulations for 20-50 years to identify parameter values and surface mass balance conditions potentially suitable to grow a realistic ice sheet. This enabled us to find ~160 simulations with good ice extent. Through analysis of these simulations, we find that albedo parameters determine the majority of uncertainty when simulating the Last Glacial Maximum North American Ice Sheets. The differences in cloud cover over the ablation zones of the North American and Greenland ice sheet explains why the ice sheets have different sensitivities to surface mass balance parameters. Based on our work, we propose that the Last Glacial Maximum can provide an “out-of-sample” target to avoid over calibrating coupled climate-ice sheet models to the present day. References: Kageyama, M. et al. The PMIP4 contribution to CMIP6 – Part 4: Scientific objectives and experimental design of the PMIP4-CMIP6 Last Glacial Maximum experiments and PMIP4 sensitivity experiments. Geosci. Model Dev. 10, 4035–4055 (2017). Smith, R. S., George, S., and Gregory, J. M.: FAMOUS version xotzt (FAMOUS-ice): a general circulation model (GCM) capable of energy- and water-conserving coupling to an ice sheet model, Geosci. Model Dev., 14, 5769–5787, https://doi.org/10.5194/gmd-14-5769-2021, 2021.

BACKGROUND:The addition of stereotactic ablative radiotherapy (SABR) to immune checkpoint inhibitors (ICIs) has the potential to significantly improve outcomes in the treatment of metastatic melanoma. We analysed peripheral blood immune cells of patients receiving combination SABR and ICI to detect the effect of treatment and identify potential biomarkers that predict outcome.METHODS:24 polymetastatic melanoma patients participated in the SABR IMPACT trial, receiving standard dose immunotherapy with anti-PD-1 and/or anti-CTLA-4 and stereotactic ablative radiotherapy to one site. Comprehensive immunophenotyping of T-cells was performed with flow cytometry on blood samples from 13 patients at baseline and following the first 4 cycles of treatment.RESULTS:Following four cycles of immunotherapy and SABR, the proportion of naïve subsets were reduced within both the CD4 and CD8 T-cell lineages. Independently, SABR resulted in increased expression of PD-1 (p = 0.019) and ICOS (p = 0.046) on the CD8+ T-cells, accompanied by a reduction in regulatory T-cell frequencies (p = 0.048). A multivariate discriminant analysis revealed a baseline signature of lower levels of CD8+ naive T-cells and higher expression of TIM-3 on regulatory T-cells and memory T-cells better predicted response.CONCLUSION:The combination of immunotherapy and SABR changed the immunophenotype of blood T cells, with some shifts attributable to SABR. Importantly, we identified a T-cell signature at baseline that best predicted response. Validation of these findings in an independent cohort could confirm these as biomarkers at baseline or early during treatment, and whether these can be utilised to stratify patients for high or low intensity treatment to reduce toxicity.

Here we present the experimental design and results from a new mid-Pliocene simulation using the latest version of the UK’s physical climate model, HadGEM3-GC31-LL, conducted under the auspices of CMIP6/PMIP4/PlioMIP2. Although two other paleoclimate simulations have been recently run using this model, they both focused on more recent periods within the Quaternary and therefore this is the first time this version of the UK model has been run this far back in time. The mid-Pliocene Warm Period, ~3 Ma, is of particular interest because it represents a time period when the Earth was in equilibrium with CO2 concentrations roughly equivalent to those of today, providing a possible analogue for current and future climate change. The implementation of the Pliocene boundary conditions is firstly described in detail, based on the PRISM4 dataset, including CO2, ozone, orography, ice mask, lakes, vegetation fractions and vegetation functional types. These were incrementally added into the model, to change from a preindustrial setup to a Pliocene setup. The results of the simulation are then presented, which are firstly compared with the model’s pre-industrial simulation, secondly with previous versions of the same model and with available proxy data, and thirdly with all other models included in PlioMIP2. Firstly, the comparison with preindustrial suggests that the Pliocene simulation is consistent with current understanding and existing work, showing warmer and wetter conditions, and with the greatest warming occurring over high latitude and polar regionsthe. The global mean surface air temperature anomaly at the end of the Pliocene simulation is 5.1°C, which is the 2nd highest of all models included in PlioMIP2 and is consistent with the fact that HadGEM3-GC31-LL has one of the highest Effective Climate Sensitivities of all CMIP6 models. Secondly, the comparison with previous generation models and with proxy data suggests a clear increase in global sea surface temperatures as the model has undergone development. Up to a certain level of warming, this results in a better agreement with available proxy data, and the “sweet spot” appears to be the previous CMIP5 generation of the model, HadGEM2-AO. The most recent simulation presented here, however, appears to show poorer agreement with the proxy data compared with HadGEM2, and may be overly sensitive to the Pliocene boundary conditions resulting in a climate that is too warm. Thirdly, the comparison with other models from PlioMIP2 further supports this conclusion, with HadGEM3-GC31-LL being one of the warmest and wettest models in all of PlioMIP2 and, if all the models are ordered according to agreement with proxy data, HadGEM3-GC31-LL ranks approximately halfway among them. A caveat to these results is the relatively short run length of the simulation, meaning the model is not in full equilibrium. Given the computational cost of the model it was not possible to run for longer; a Gregory plot analysis indicates that had it been allowed to come to full equilibrium, the final global mean surface temperature could have been approximately 1.5°C higher.

As part of a project working to improve coupled climate-ice sheet modelling of the response of ice sheets to changes in climate across different periods since the Last Glacial Maximum, we present simulations of the modern Greenland climate and ice sheet using the FAMOUS-BISICLES model. FAMOUS-BISICLES, a variant of FAMOUS-ice (Smith et al., 2021a), is a low resolution (7.5°X5°) global climate model that is two-way coupled to a higher resolution (minimum grid spacing of 1.2 km) adaptive mesh ice sheet model, BISICLES. It uses a system of elevation classes to downscale the lower resolution atmospheric variables onto the ice sheet grid and calculates surface mass balance using a multilayer snow model. FAMOUS-ice is computationally affordable enough to simulate the millennial evolution of the coupled climate-ice sheet system, and has been shown to simulate Greenland well in previous work using the Glimmer shallow ice model (Gregory et al., 2020). The ice sheet volume and area are sensitive to a number of parametrisations related to atmospheric and snow surface processes and ice sheet dynamics. Based on that, we designed a perturbed parameters ensemble using a Latin Hypercube sampling technique and ran simulations with climate forcings appropriate for the late 20th century. The ice sheet area and volume are most correlated to parameters that set the snow/firn albedo while the relationship is less simple for parameters related to clouds and precipitation. We compare FAMOUS-ice SMB and coupled behaviour against the more sophisticated, higher resolution, CMIP6-class UKESM-ice coupled climate ice sheet model for a late 20th century simulation as well as an abrupt 4XCO2 experiment. Our simulations produce a large range of climate and ice sheet behaviours, including a stable control state for the modern Greenland, and we have been able to highlight the sensitivity of the system to other sets of parameters and future changes in climate.

Understanding the response of ice sheets to global temperature changes is a critical issue for the climate community. To accurately simulate future ice sheet evolution, we need to know the strength of feedbacks between the climate and ice sheets. Testing the ability of coupled climate-ice sheet models to simulate past ice sheet extent can provide a way to evaluate the models and ground truth projections. One example is the Last Glacial Maximum (LGM), when huge ice sheets covered the Northern Hemisphere, especially over the North America. Here, we performed simulations of the North American ice sheet and climate of the LGM with a recently updated ice sheet-atmosphere coupled model Famous-Ice (Smith et al. 2021, Gregory et al. 2020). The model consists of a low-resolution atmospheric general circulation model Famous (Smith et al. 2008) and an ice sheet model BISICLES (Cornford et al. 2013). It calculates the surface mass balance over ice sheets based on an energy budget scheme and incorporates an updated albedo scheme, which accounts for albedo changes associated with modifications in surface air temperature, grain size and density of the snow. The atmospheric model reproduces the surface mass balance of the modern Greenland ice sheet reasonably well (Smith et al. 2021). Simulations of projections of future sea-level rise (Gregory et al. 2020) and the LGM (Gandy et al. in prep) have also been performed with Famous-Ice using a different ice sheet model GLIMMER. We present simulations of the LGM with interactive ice sheets in North America and Greenland using FAMOUS-BISICLES. Uncertain input parameters controlling the surface temperatures and ice albedo are varied in our simulations. The global temperature is specified by applying fixed sea surface temperature in the atmospheric model producing a global cooling that ranges from -3K to -6.5K in the simulations. The bare ice minimum albedo is varied from 0.2 to 0.7, which corresponds to the range in PMIP3 models. Our results show a better representation of North American ice sheet when forced with a colder LGM (-6.5K) and high bare ice albedo. We will further discuss potential roles of model biases and compare our results with simulations performed with FAMOUS-GLIMMER (Gandy et al. in prep).

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.

To better understand the processes contributing to future climate change, palaeoclimate model simulations are an important tool because they allow testing of the models’ ability to simulate very different climates than that of today. As part of CMIP6/PMIP4, the latest version of the UK’s physical climate model, HadGEM3-GC31-LL (hereafter, for brevity, HadGEM3), was recently used to simulate the mid-Holocene (~6 ka) and Last Interglacial (~127 ka) simulations and the results were compared to the preindustrial era, previous versions of the same model and proxy data (see Williams et al. 2020, Climate of the Past). Here, we use the same model to go further back in time, presenting the results from the mid-Pliocene Warm Period (~3.3 to 3 ma, hereafter the “Pliocene” for brevity). This period is of particular interest when it comes to projections of future climate change under various scenarios of CO2 emissions, because it is the most recent time in Earth’s history when CO2 levels were roughly equivalent to today. In response, albeit due to slower mechanisms than today’s anthropogenic fossil fuel driven-change, during the Pliocene global mean temperatures were 2-3°C higher than today, more so at the poles. Here, we present results from the HadGEM3 Pliocene simulation. The model is responding to the Pliocene boundary conditions in a manner consistent with current understanding and existing literature. When compared to the preindustrial era, global mean temperatures are currently ~5°C higher, with the majority of warming coming from high latitudes due to polar amplification from a lack of sea ice. Relative to other models within the Pliocene Modelling Intercomparison Project (PlioMIP), this is the 2nd warmest model, with the majority of others only showing up to a 4.5°C increase and many a lot less. This is consistent with the relatively high sensitivity of HadGEM3, relative to other CMIP6-class models. When compared to a previous generation of the same UK model, HadCM3, similar patterns of both surface temperature and precipitation changes are shown (relative to preindustrial). Moreover, when the simulations are compared to proxy data, the results suggest that the HadGEM3 Pliocene simulation is closer to the reconstructions than its predecessor.