Search

Understanding grounding line dynamics is critical for projecting glacier evolution and sea level rise. Observations from satellite radar interferometry reveal rapid grounding line migration forced by oceanic tides that are several kilometers larger than predicted by hydrostatic equilibrium, indicating the transition from grounded to floating ice is more complex than previously thought. Recent studies suggest seawater intrusion beneath grounded ice may play a role in driving rapid ice loss. Here, we investigate its impact on the evolution of Petermann Glacier, Greenland, using an ice sheet model. We compare model results with observed changes in grounding line position, velocity, and ice elevation between 2010 and 2022. We match the observed retreat, speed up, and thinning using 3‐km‐long seawater intrusion that drive peak ice melt rates of 50 m/yr; but we cannot obtain the same agreement without seawater intrusion. Including seawater intrusion in glacier modeling will increase the sensitivity to ocean warming. Plain Language Summary: Relatively warm seawater melts marine‐terminating glaciers from below. Recent observations suggest that seawater flows below grounded ice at high tide. The presence of seawater at this boundary, referred to as seawater intrusion, has the potential to increase glacier mass loss. We test this hypothesis on Petermann Glacier, Greenland, using an ice sheet flow model. We run the model to reconstruct the glacier's behavior from 2010 to 2022 with and without seawater intrusion. We compare the results with satellite observations of velocity, grounding line position, and ice thinning. When we use enhanced ice melt rates from kilometer‐scale seawater intrusion, we match the observed retreat, speed up, and thinning. When we do not, the model fails to replicate the observations. Seawater intrusion may play a critical role in glacier evolution. Adding this process to ice flow models will increase their sensitivity to ocean warming and projections of ice mass loss and sea level rise. Key Points: Ice melt caused by seawater intrusion in the grounding zone explains the observed grounding line retreat of Petermann GlacierWithout seawater intrusion in the grounding zone, we do not replicate the full extent of observed retreatIncluding seawater intrusion in the grounding zone increases glacier mass loss [ABSTRACT FROM AUTHOR]

On behalf of the journal, AGU, and the scientific community, the editors of Geophysical Research Letters would like to sincerely thank those who reviewed manuscripts for us in 2023. The hours reading and commenting on manuscripts not only improve the manuscripts, but also increase the scientific rigor of future research in the field. With the advent of AGU's data policy, many reviewers have also helped immensely to evaluate the accessibility and availability of data, and many have provided insightful comments that helped to improve the data presentation and quality. We greatly appreciate the assistance of the reviewers in advancing open science, which is a key objective of AGU's data policy. We particularly appreciate the timely reviews in light of the demands imposed by the rapid review process at Geophysical Research Letters. We received 4,512 submissions in 2023 and 5,112 reviewers contributed to their evaluation by providing 8,587 reviews in total. We deeply appreciate their contributions. Plain Language Summary: Individuals in italics provided three or more reviews for GRL in 2023. Key Points: The editors thank the 2023 peer‐reviewers [ABSTRACT FROM AUTHOR]

A realistic initialization of ice flow models is critical for predicting future changes in ice sheet mass balance and their associated contribution to sea level rise. The initial thermal state of an ice sheet is particularly important, as it controls ice viscosity and basal conditions, thereby influencing the overall ice velocity. Englacial and subglacial conditions, however, remain poorly understood due to insufficient direct measurements, which complicate the initialization and validation of thermal models. Here, we investigate the impact of using different geothermal heat flux (GHF) datasets and vertical velocity profiles on the thermal state of the Antarctic ice sheet and compare our modeled temperatures to in situ measurements from 15 boreholes. We find that the temperature profile is more sensitive to vertical velocity than to GHF. The basal temperature of grounded ice and the amount of basal melting are influenced by both selection of GHF and vertical velocity. More importantly, we find that the standard approach, which consists of combining basal sliding speed and incompressibility to derive vertical velocities, provides reasonably good results in fast-flow regions (ice velocity >50 m yr -1) but performs poorly in slow-flow regions (ice velocity <50 m yr -1). Furthermore, the modeled temperature profiles in ice streams, where bed geometry is generated using a mass conservation approach, show better agreement with observed borehole temperatures compared to kriging-based bed geometry. [ABSTRACT FROM AUTHOR]

Ice-sheet flow models capable of accurately projecting their future mass balance constitute tools to improve flood risk assessment and assist sea-level rise mitigation associated with enhanced ice discharge. Some processes that need to be captured, such as grounding-line migration, require high spatial resolution (under the kilometer scale). Conventional ice flow models mainly execute on central processing units (CPUs), which feature limited parallel processing capabilities and peak memory bandwidth. This may hinder model scalability and result in long run times, requiring significant computational resources. As an alternative, graphics processing units (GPUs) are ideally suited for high spatial resolution, as the calculations can be performed concurrently by thousands of threads, processing most of the computational domain simultaneously. In this study, we combine a GPU-based approach with the pseudo-transient (PT) method, an accelerated iterative and matrix-free solution strategy, and investigate its performance for finite elements and unstructured meshes with application to two-dimensional (2-D) models of real glaciers at a regional scale. For both the Jakobshavn and Pine Island glacier models, the number of nonlinear PT iterations required to converge a given number of vertices (N) scales in the order of O(N1.2) or better. We further compare the performance of the PT CUDA C implementation with a standard finite-element CPU-based implementation using the price-to-performance metric. The price of a single Tesla V100 GPU is 1.5 times that of two Intel Xeon Gold 6140 CPUs. We expect a minimum speedup of at least 1.5 times to justify the Tesla V100 GPU price to performance. Our developments result in a GPU-based implementation that achieves this goal with a speedup beyond 1.5 times. This study represents a first step toward leveraging GPU processing power, enabling more accurate polar ice discharge predictions. The insights gained will benefit efforts to diminish spatial resolution constraints at higher computing performance. The higher computing performance will allow for ensembles of ice-sheet flow simulations to be run at the continental scale and higher resolution, a previously challenging task. The advances will further enable the quantification of model sensitivity to changes in upcoming climate forcings. These findings will significantly benefit process-oriented sea-level-projection studies over the coming decades. [ABSTRACT FROM AUTHOR]

Rock avalanche-triggered displacement waves (also termed tsunamis) have recently occurred in Greenland and Alaska, and they illustrate the presence of such hazards in polar regions. To improve understanding of the magnitude of this hazard for these areas, we investigated gigascale subaerial rock avalanches impacting a partially confined water body within the Vaigat strait (western Greenland). We present a new combined subaerial to subaqueous digital elevation model, alongside a new compilation of seismic data, which revealed nine deglacial to Holocene rock avalanche complexes that are between one and two orders of magnitude larger than nearby historical rock avalanches. The three largest complexes have deposit thicknesses up to 300 m, runout distances reaching 19 km, and best-estimate volumes from 1.7 to 8.4 km3. Based on the morphology and the volume-angle of reach relations, it is likely that each complex represents a single or few events, thus making them among the largest displacement wave-generating subaerial to submarine rock avalanches on Earth. We estimated displacement wave magnitude up to 280 m on the opposite shore. The ages of the deposits are poorly constrained but the main rock avalanche activity is referable to early Holocene times. With significant climatic changes predicted in the Arctic, we recommend that hazard assessments account for events not only from the historical record but also those from the recent geological past. [ABSTRACT FROM AUTHOR]