Neural operators for accelerating UQ: application to ice-sheet modeling
Explore the source record for details and available documents.
Perego, Mauro [Sandia National Laboratories (SNL-N↗
SEARCH · Search NASA
Results for “ice sheets”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Results of the second Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+)
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 (Asay-Davis et al., 2016), 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.
54 ENVIRONMENTAL SCIENCES↗
Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt
Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.
54 ENVIRONMENTAL SCIENCES↗
Glacial geology of the Hudson Mountains, Amundsen Sea sector, West Antarctica
The Hudson Mountains are situated in the eastern Amundsen Sea sector of the West Antarctic Ice Sheet, adjacent to Pine Island Glacier. They form a volcanic field of 17 stratovolcanoes and parasitic vents, preserved as nunataks. Two former tributaries of Pine Island Glacier (Larter and Lucchitta glaciers) flow through the mountains. Here we present a detailed study of the glacial geology of the area. We describe field observations and measurements of geomorphological features from 15 of the nunataks, meltwater ponds found on the surface of three nunataks and supraglacial features (ice dolines) from two sites near the present grounding line. Together these provide constraints on the past ice sheet extent, flow pathways and thermal regime, and enhance our understanding of the present hydrological regime – all of which are important as context for the observed modern ice sheet behaviour. We find evidence suggesting that all nunataks in the Hudson Mountains were covered by ice during the Last Glacial Maximum (defined here as 26.5-19 ka) and have since deglaciated. Faceted and polished erratic cobbles and boulders of exotic lithologies (syenites, alkali granites, granites, granodiorites, tonalites and gabbros) are numerous and perched on nunatak surfaces. A marked difference between the dominant erratic lithologies on nunataks adjacent to Pine Island Glacier (granite) and Lucchitta Glacier (granodiorite-tonalite) indicates that the ice sheet was transporting clasts from at least two distinct upstream source regions. The similarity in degree of weathering suggests, however, that all the erratics were transported by one phase of (warm-based) glaciation; their presence on or close to the summits of all except one nunatak indicates that the ice sheet during that time was at least 700 m thicker than present. These results are consistent with ice sheet model simulations which suggest that all nunataks in the Hudson Mountains were completely submerged by the Last Glacial Maximum ice sheet.
58 GEOSCIENCES↗
Direct Geologic Constraints on the Timing of Late Holocene Ice Thickening in the Amundsen Sea Embayment, Antarctica
Abstract Constraining past West Antarctic Ice Sheet (WAIS) change helps validate numerical models simulating future ice sheet dynamics. Following rapid deglaciation during the mid‐Holocene, ice near Thwaites Glacier was ∼35 m thinner than present; however, the timing of ice regrowth to its present configuration remains unknown. To fill this knowledge gap, we present cosmogenic nuclide exposure ages of cobbles from the surface of a moraine situated between Thwaites and Pope glaciers. We infer that the moraine formed and stabilized in the Late Holocene (∼1.4 ka) when a small glacier thickened. We also present a novel reconstruction of WAIS volume constrained by sea‐level data, which demonstrates that moraine formation coincided with a large‐scale WAIS readvance. Our new geologic constraints will help inform models of the solid Earth response to surface mass loading, improving our understanding of ice sheet dynamics in a vulnerable part of WAIS.
Nichols, Keir A. [Imperial College London London U↗
Assessing the suitability of sites near Pine Island Glacier for subglacial bedrock drilling aimed at detecting Holocene retreat–readvance
Abstract. Unambiguous identification of past episodes of ice sheet thinning below the modern surface and grounding line retreat inboard of present requires recovery and exposure dating of subglacial bedrock. Such efforts are needed to understand the significance and potential future reversibility of ongoing and projected change in Antarctica. Here we evaluate the suitability for subglacial bedrock drilling of sites in the Hudson Mountains, which are located in the Amundsen Sea sector of West Antarctica. We use an ice sheet model and field data – geological observations, glaciological observations and bedrock samples from nunataks, and ground-penetrating radar from subglacial ridges – to rate each site against four key criteria: (i) presence of ridges extending below the ice sheet, (ii) likelihood of increased exposure of those ridges if the grounding line was inboard of present, (iii) suitability of bedrock for drilling and geochemical analysis, and (iv) accessibility for aircraft and drilling operations. Our results demonstrate that although no site in the Hudson Mountains is perfect for this study when assessed against all criteria, the accessibility, N–S orientation and basaltic bedrock lithology of Winkie Nunatak's southernmost ridge (74.86° S, 99.77° W) make it a feasible site both for drilling and subsequent cosmogenic nuclide analysis. Furthermore, the ridge is strewn with glacial erratics at all elevations, providing valuable constraints on its early Holocene deglacial history. Based on our experiences during this study, we conclude with a series of recommendations for assessing site suitability for future bedrock drilling campaigns. We emphasise the importance of consulting a range of expertise prior to drilling and ensuring that sufficient field reconnaissance is undertaken (including obtaining detailed grids of radar survey data and bedrock samples).
58 GEOSCIENCES↗
Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model
Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.
58 GEOSCIENCES↗
Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model
Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.
54 ENVIRONMENTAL SCIENCES↗
The Association Between Cloud Droplet Number over the Summer Southern Ocean and Air Mass History
The cloud properties and governing processes in Southern Ocean marine boundary layer clouds have emerged as a central issue in understanding the Earth's climate sensitivity. While our understanding of Southern Ocean cloud feedbacks have evolved in the most recent climate model intercomparison, the background properties of simulated summertime clouds in the Southern Ocean are not consistent with measurements due to known biases in simulating cloud condensation nuclei concentrations. This paper presents several case studies collected during the Capricorn 2 and Marcus campaigns held aboard Australian research vessels in the Austral Summer of 2018. Combining the surface–observed cases with MODIS data along forward and backward air mass trajectories, we demonstrate the evolution of cloud properties with time. These cases are consistent with multi–year statistics showing that long trajectories of air masses over the Antarctic ice sheet are critical to creating high droplet number clouds in the high latitude summer Southern Ocean. We speculate that secondary aerosol production via the oxidation of biogenically derived aerosol precursor gasses over the high actinic flux region of the high latitude ice sheets is fundamental to maintaining relatively high droplet numbers in Southern Ocean clouds during Summer.
54 ENVIRONMENTAL SCIENCES↗
Weather and climate extremes in a changing Arctic
Weather and climate extremes are increasingly occurring in the Arctic. Here, in this Review, we evaluate historical and projected changes in rare Arctic extremes across the atmosphere, cryosphere and ocean and elucidate their driving mechanisms. Clear shifts occur in mean and extreme distributions after ~2000. For instance, pre-2000 to post-2000 observational probabilities of 1.5 standard deviation events increase by 20% for atmospheric heat waves, 76.7% for Atlantic layer warm events, 83.5% for Arctic sea ice loss and 62.9% for Greenland Ice Sheet melt extent — in many cases, low probability, rare extreme events in the early period become the norm in the latter period. These observed changes can be explained using a ‘pushing and triggering’ concept, representing interplay between external forcing and internal variability: long-term warming destabilizes the climate system and ‘pushes’ it to a new state, allowing subsequent variability associated with large-scale atmosphere–ocean–ice interactions and synoptic systems to ‘trigger’ extreme events over different timescales. Ongoing anthropogenic warming is expected to further increase the frequency and magnitude of extremes, such that simulated probabilities of 1.5 standard deviation events increase by 72.6% for atmospheric heat waves, 68.7% for Atlantic layer warm events and 93.3% for Greenland Ice Sheet melt rate between historic (1984–2014) and future (2069–2099) periods under a very high emission scenario. Future research should prioritize the development of physically based metrics, enhance high-resolution observation and modelling capabilities and improve understanding of multiscale Arctic climate drivers.
atmospheric dynamics↗
Rates of Sea‐Level Rise Are Highly Sensitive to Ice Viscosity Parameters in Model Benchmarks
Glacier flow plays a major role in current and future rates of globally averaged sea-level rise. The viscosity of glacial ice, controlling the rate of flow, decreases as stress increases and is highly sensitive to the value of the stress exponent, $n$, in the constitutive equation for viscous flow. Glaciologists and climate modelers almost exclusively assume $n=3$ when modeling ice flow and projecting sea-level rise through forward modeling. However, recent work suggests that $n\approx 4$ better fits observations, prompting the question: How sensitive are projections of sea-level rise to the value of $n$? We use an established community ice flow model and standard benchmark experiments designed as an idealized representation of Pine Island Glacier, West Antarctica. While initializing an $n=3$ model to match observations of an $n=4$ ice sheet is possible, we find that incorrectly assuming $n=3$ when in fact $n=4$ dramatically underestimates rates of sea-level rise. The scale of this error grows nonlinearly with the magnitude of the climate forcing, acting to increase projection uncertainties. Additionally, we find that models often account for this stress-dependent rheology mismatch during model initialization in a way that masks this rheological effect in the short term while leaving model outputs vulnerable to larger biases in longer-term projections. Initializations to observations of Pine Island Glacier display similar rheology-mismatch fingerprints to our idealized example.
climate sensitivity↗
A high-resolution pseudo-polygon discrete element model for regional sea ice
Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.
54 ENVIRONMENTAL SCIENCES↗
Basal Melting and Oceanic Observations Beneath Central Fimbulisen, East Antarctica
Abstract Basal melting of ice shelves is fundamental to Antarctic ice sheet mass loss, yet direct observations remain sparse. We present the first year‐round melt record (2017–2021) from a phase‐sensitive radar on Fimbulisen, one of the fastest flowing ice shelves in Dronning Maud Land, East Antarctica. The observed long‐term mean ablation rate at 350 m depth below the central ice shelf was 1.0 ± 0.5 m yr −1 , marked by substantial sub‐weekly variability ranging from 0.4 to 3.5 m yr −1 . 36‐h filtered basal melt rate fluctuations closely align with ocean velocity. On seasonal time scales, melt rates peak during austral spring to autumn (September–March), driven by both elevated ocean velocities and thermal driving near the base. The combined effect of thermal driving and current speed explains the majority of the melt rate variability ( r = 0.84), highlighting the dominant role of shear‐driven turbulence. This relationship enables parameterization of melt rates for the decade‐long ocean record (2010–2021), although deviations appear under low and high forcing conditions. Both observed and parameterized melt rates show similar yearly mean magnitudes compared to satellite‐derived melt rates but with a tenfold lower seasonal amplitude and a 3‐month delay in seasonality. These detailed concurrent ice–ocean observations provide essential validation data for remote sensing and numerical models that aim to quantify and project ice‐shelf response to a change in ocean forcing. In situ measurements and continued monitoring are crucial for accurately assessing and modeling future basal melt rates, and for understanding the complex dynamics driving ice‐shelf stability and sea‐level change.
54 ENVIRONMENTAL SCIENCES↗
Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0
Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.
58 GEOSCIENCES↗
Preliminary Screening of Features, Events, and Processes for an Arctic-Focused Climate Intervention Performance Assessment
Geoengineering, the deliberate large-scale intervention in Earth's climate system, holds significant potential in the rapidly warming Arctic, where temperatures currently rise at more than twice the global average, accelerating ice sheet and permafrost melt. This contributes to global sea-level rise and releases methane, a potent greenhouse gas. Strategies like solar radiation management (SRM) and carbon dioxide removal (CDR) could mitigate these effects; for instance, SRM techniques aim to reflect a portion of the sun's energy back into space, potentially slowing ice melt and stabilizing permafrost. However, geoengineering in the Arctic faces challenges, including potential unintended consequences on the fragile ecosystem, disruption of local weather patterns, and impacts on indigenous communities. Effective governance requires robust international cooperation, environmental impact assessments, and regulatory frameworks. Despite these challenges, geoengineering's potential benefits make it a critical research area. This report explores application of the Performance Assessment (PA) methodology to Arctic Climate Intervention, providing an initial screening of relevant features, events, and processes (FEPs). At the core of the PA approach is the identification and evaluation of FEPs that could impact the performance of the intervention scheme. Here we provide an initial screening of FEPs to consider in the application of PA to Arctic Climate Intervention.
54 ENVIRONMENTAL SCIENCES↗
Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM
The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.
58 GEOSCIENCES↗
Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers
Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.
atmospheric rivers↗
Subglacial Discharge Effects on Antarctic Ice‐Shelf Basal Melt and the Southern Ocean in a Global, Coupled Ocean—Sea‐Ice Model
Subglacial freshwater from beneath Antarctica enters the ocean at depth, enhancing ice-shelf melting and affecting Southern Ocean properties. To study these effects in an Antarctic-wide context, we use a continental-scale subglacial hydrology model that calculates grounding line freshwater flux for a global, coupled ocean—sea-ice model. We find that subglacial discharge impacts melt rates primarily through continental shelf temperature modification, contrasting with findings from regional studies that do not permit large-scale adjustments. The consequence is that Antarctic melt rates scale with subglacial discharge more strongly than inferred from regional studies. We also find that the addition of buoyancy at depth facilitates heat upwelling to the surface, resulting in higher sea ice volume downstream of cold ice shelves and lower sea ice volume downstream of warm ice shelves. This highlights the drawbacks of simplifications in previous global studies that deposit Antarctic meltwater at the ocean surface and find uniform ocean surface cooling and sea-ice growth. While the patterns we find are robust, we conclude that the addition of subglacial discharge at present-day rates has a small effect on basal melt rates, hydrography, and sea ice. However, stronger discharge can have significant effects and can even accelerate a shift from low to high melting for ice shelves close to such a tipping point. Finally, the importance of feedbacks between enhanced cavity overturning and continental shelf conditions poses a complication for parameterizing subglacial discharge effects on melting for ice-sheet projections that do not include a coupled ocean component.
54 ENVIRONMENTAL SCIENCES↗