Search NASA⌕ Search

SEARCH · Search NASA

Results for “glacier”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Black Carbon Deposited in Hariqin Glacier of the Central Tibetan Plateau Record Changes in the Emission from Eurasia

Black carbon (BC), by the combustion of fossil fuels and biomass, has profound effects on climate change and glacier retreat in industrial eras. In the present study, we report refractory BC (rBC) in an ice core spanning 1850–2014, retrieved from the Hariqin Glacier of the Tanggula Mountains in the central Tibetan Plateau, measured using a single particle soot photometer (SP2). The rBC concentration shows a three-fold increase since the 1950s. The mean rBC concentration was 0.71±0.52 ng mL-1 during 1850s–1940s and 2.11±1.60 ng mL-1 during 1950s–2010s. The substantial increase in rBC since the 1950s is consistent with rBC ice core records from the Tibetan Plateau and Eastern Europe. According to the predominant atmospheric circulation patterns over the glacier and timing of changes in regional emissions, the post-1950 amplification of rBC concentration in the central Tibetan Plateau most likely reflects increases in emissions in Eastern Europe, former USSR, the Middle East, and South Asia. Despite the low-level background rBC concentrations in the ice cores from the Tibetan Plateau, the present study highlights a remarkable increase in anthropogenic BC emissions in recent decades and the consequent influence on glaciers in the Tibetan Plateau.

Black carbon, Ice core, Tibetan Plateau, Industria↗
DOE OSTI2021

Holocene thinning of Darwin and Hatherton glaciers, Antarctica, and implications for grounding-line retreat in the Ross Sea

Chronologies of glacier deposits in the Transantarctic Mountains provide important constraints on grounding-line retreat during the last deglaciation in the Ross Sea. However, between Beardmore Glacier and Ross Island – a distance of some 600 km – the existing chronologies are generally sparse and far from the modern grounding line, leaving the past dynamics of this vast region largely unconstrained. We present exposure ages of glacial deposits at three locations alongside the Darwin–Hatherton Glacier System – including within 10 km of the modern grounding line – that record several hundred meters of Late Pleistocene to Early Holocene thickening relative to present. As the ice sheet grounding line in the Ross Sea retreated, Hatherton Glacier thinned steadily from about 9 until about 3 ka. Our data are equivocal about the maximum thickness and Mid-Holocene to Early Holocene history at the mouth of Darwin Glacier, allowing for two conflicting deglaciation scenarios: (1) ~500 m of thinning from 9 to 3 ka, similar to Hatherton Glacier, or (2) ~950 m of thinning, with a rapid pulse of ~600 m thinning at around 5 ka. We test these two scenarios using a 1.5-dimensional flowband model, forced by ice thickness changes at the mouth of Darwin Glacier and evaluated by fit to the chronology of deposits at Hatherton Glacier. The constraints from Hatherton Glacier are consistent with the interpretation that the mouth of Darwin Glacier thinned steadily by ~500 m from 9 to 3 ka. Rapid pulses of thinning at the mouth of Darwin Glacier are ruled out by the data at Hatherton Glacier. This contrasts with some of the available records from the mouths of other outlet glaciers in the Transantarctic Mountains, many of which thinned by hundreds of meters over roughly a 1000-year period in the Early Holocene. The deglaciation histories of Darwin and Hatherton glaciers are best matched by a steady decrease in catchment area through the Holocene, suggesting that Byrd and/or Mulock glaciers may have captured roughly half of the catchment area of Darwin and Hatherton glaciers during the last deglaciation. An ensemble of three-dimensional ice sheet model simulations suggest that Darwin and Hatherton glaciers are strongly buttressed by convergent flow with ice from neighboring Byrd and Mulock glaciers, and by lateral drag past Minna Bluff, which could have led to a pattern of retreat distinct from other glaciers throughout the Transantarctic Mountains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Streamflow Response to Glacier Mass Loss Varies With Basin Precipitation Across Alaska

Abstract Diminishing glaciers affect streamflow, and given the extent of glaciers in Alaska and adjacent Canada, continued glacier mass loss is likely to have profound effects on ecosystems sensitive to runoff. The effects of glacier mass loss on streamflow are likely to vary across the wide ranges of basin size, glacier cover, and precipitation in this region. In this study, we use U.S. Geological Survey (USGS) streamflow data with satellite‐based glacier volume change estimates to quantify how glacier mass loss subsidized streamflow over the 2000–2019 period for 116 glacierized basins. We examine interannual variability in that subsidy at three USGS‐monitored glaciers to explore the ability of the subsidy to buffer streamflow derived solely from precipitation. We found the relative importance of percent glacier cover on streamflow magnitude increases in drier basins. In the driest basins, glaciers produced 40 times greater percent glacier mass loss subsidies to streamflow for the percent glacier cover compared to the wettest basins. While the subsidy from glacier mass loss buffers interannual variability in streamflow to varying degrees, it can also increase streamflow variability. Smaller amounts of percent glacier cover are needed to produce summer‐melt‐dominated seasonal flow regimes in drier basins than in wetter basins. Decreasing glacier cover will eventually decrease summer streamflow, increasing spring streamflow in drier basins, and attenuating seasonality with increasing spring and autumnal streamflow in wetter basins. Quantifying the downstream effects of continued glacier mass loss without the computational expense of a hydrological model is broadly applicable in this changing climate.

Curran, Janet H. [U.S. Geological Survey Alaska Sc↗
DOE OSTI2023

Late Holocene glacier and climate fluctuations in the Mackenzie and Selwyn mountain ranges, northwestern Canada

Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10 Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10 Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 °C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km 3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10 Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

Abstract. 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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Glacier ice archives nearly 15,000-year-old microbes and phages

Background Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination. Results We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants. Conclusions Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Glaciers in California’s Sierra Nevada are likely disappearing for the first time in the Holocene

Mountain glaciers are sensitive climate indicators. Glaciers in the western United States are projected to disappear by 2100 CE, but whether they were previously absent in the current interglacial (past 11,700 years) remains debated. Here, we present evidence that two of the largest glaciers in California’s Sierra Nevada near Yosemite National Park persisted throughout the Holocene. Cosmogenic in situ carbon-14 and beryllium-10 exposure dating in newly exposed proglacial bedrock indicates continuous Holocene cover, likely by ice. At a nearby smaller glacier, bedrock exposure-burial ages suggest glacier expansion at ~7 thousand years ago, earlier than previously recognized. Moraine exposure ages at these sites and the Sierra Nevada’s largest glacier are distributed across the past several thousand years, suggesting that glaciers were near their preindustrial extent for much of the late Holocene. These findings imply that a glacier-free Sierra Nevada is unprecedented since before the Holocene. Future alpine glacial environments in California are thereby likely a no-analog scenario within the current interglacial.

Jones, Andrew G. [University of Wisconsin, Madison↗
DOE OSTI2021

Local summer insolation and greenhouse gas forcing drove warming and glacier retreat in New Zealand during the Holocene

Geological climate archives from the Holocene Epoch provide baseline information concerning natural climate variability. Temperate mountain glacier extent is sensitive to summer air temperature, thus geological records of past glacier length changes are a useful proxy for this climatic variable. Here we present a new cosmogenic 10 Be chronology of glacier length changes at Dart Glacier in the Southern Alps, New Zealand. Prominent moraines deposited 321 ± 44 yr ago (n = 11) and 7.8 ± 0.3 ka (n = 5) show glaciers during the Little Ice Age were less extensive than during the early Holocene. This pattern of net Holocene glacier retreat is consistent with emerging data from other catchments in New Zealand and across the southern mid-latitudes. Using the physical framework of a transient global climate model simulation, we suggest that cool summers in the early Holocene were promoted by the local summer insolation minimum, together with low atmospheric greenhouse gas concentrations, causing an early Holocene austral glacial maximum. An insolation-driven reduction in seasonality at southern mid-latitudes may reconcile differences between early Holocene temperature reconstructions where climate proxies have different seasonal sensitivities. Furthermore, we suggest that rising greenhouse gas concentrations after 7 ka caused regional-scale glacier retreat and appear to be the dominant driver of multi-millennial summer temperature trends in the southern mid-latitudes during the present interglacial.

58 GEOSCIENCES↗
DOE OSTI2022

Persistent, extensive channelized drainage modeled beneath Thwaites Glacier, West Antarctica

Abstract. Subglacial hydrology is a leading control on basal friction and the dynamics of glaciers and ice sheets. At low discharge, subglacial water flows through high-pressure, sheet-like systems that lead to low effective pressures. However, at high discharge, subglacial water melts the overlying ice into localized channels that efficiently remove water from the bed, thereby increasing effective pressure and basal friction. Recent observations suggest channelized subglacial flow exists beneath Thwaites Glacier, yet it remains unclear if stable channelization is feasible in West Antarctica, where surface melting is nonexistent and water at the bed is limited. Here, we use the MPAS-Albany Land Ice model to run a suite of over 130 subglacial hydrology simulations of Thwaites Glacier across a wide range of physical parameter choices to assess the likelihood of channelization. We then narrow our range of viable simulations by comparing modeled water thicknesses to previously observed radar specularity content, which indicates flat, spatially extensive water bodies at the bed. In all of our data-compatible simulations, stable channels reliably form within 100–200 km of the grounding line and reach individual discharge rates of 35–110 m3 s−1 at the ice–ocean boundary. While only one to two channels typically form across the 200 km width of the glacier in our simulations, their high efficiency drains water across the entire lateral extent of the glacier. We posit the large catchment size of Thwaites Glacier, its funnel-like geometry, and high basal melt rates together accumulate enough water to form stable channels. No simulations resembled observed specularity content when channelization was disabled. Our results suggest channelized subglacial hydrology has two consequences for Thwaites Glacier dynamics: (i) amplifying submarine melting of the terminus and ice shelf while (ii) simultaneously raising effective pressure within 100 km of the grounding line and increasing basal friction. The distribution of effective pressure implied from our modeling differs from parameterizations typically used in large-scale ice sheet models, suggesting the development of more process-based parameterizations may be necessary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗
DOE OSTI2023

Early glacier advance in New Zealand during the Antarctic Cold Reversal

ABSTRACT Glacial landscapes preserve records of past climate change. Investigating the glacier–climate system over the Late Quaternary provides information about past climate change and context for present‐day glacier response to climate warming. Using 28 beryllium‐10 ( 10 Be) surface exposure dates and snowline reconstructions, we present glacier fluctuations and climate changes for the Antarctic Cold Reversal in the Ahuriri River catchment, Southern Alps of New Zealand (44°7′50″S, 169°38′29″E). Prominent terminal and lateral moraine features from the upper right tributary of the Ahuriri River valley have exposure ages of 14.5 ± 0.3, 13.6 ± 0.3 and 12.6 ± 0.2 ka, suggesting retreat of the glacier during the Antarctic Cold Reversal. Maximum elevation of lateral moraines (MELM) and accumulation area ratio (AAR) suggest snowline elevations at these ages were ≤700, ≤630 and ~360 m lower than today, respectively. This equates to air temperatures ≤3.9, ≤3.5 and 2.3 ± 0.7 °C lower than today (1981–2010), assuming no changes in past precipitation. Ice‐sculpted bedrock surfaces bound by a lateral moraine at nearby Canyon Creek have an age of 13.1 ± 0.3 ka, indicating the moraine correlates with those in the Ahuriri upper right tributary. MELM and AAR reconstructions from the Canyon Creek suggest that snowline elevations at 14.5–13.6 ka were ≤500 or ~380 m lower than today, corresponding to air temperatures ≤2.8 or 2.4 ± 0.7 °C lower than the present‐day (1981–2010). Our results provide insight into the structure of the Antarctic Cold Reversal in the Southern Alps, showing that the largest glacier advance occurred at the start of this interval at c . 14.5 ± 0.3 ka and was followed by gradual retreat. We hypothesize that the early cooling and glacier readvance in New Zealand at the onset of the Antarctic Cold Reversal were triggered by a latitudinal shift of the Southern Hemisphere westerly wind belt.

58 GEOSCIENCES↗
DOE OSTI2021

Climate reconstructions for the Last Glacial Maximum from a simple cirque glacier in Fiordland, New Zealand

Mountain glacier records offer important constraints on the timing and magnitude of climate variability during the last glacial cycle. Existing moraine chronologies from the central Southern Alps indicate maximum ice extent was achieved during marine isotope stages 3–4, followed by repeated advances of similar, but gradually declining extent during marine isotope stage 2, until onset of the glacial termination. Questions remain over the precise role of climate in driving these changes, as most existing moraine chronologies come from large, complex former valley glacier systems, where non-climatic influences such as changing bed topography and proglacial lake formation may have influenced glacier length changes. In this work, we address this problem via a new cosmogenic 10Be chronology and equilibrium line altitude reconstruction from a cirque glacier situated in Fiordland, New Zealand. Our chronology shows moraine deposition at 32 ± 11 ka, 18.7 ± 0.2 ka, 18.1 ± 0.1 ka, and c. 17.2 ± 0.3 ka. The simple geometry of the former glacier supports the role of climate in driving a net decline of regional ice volume during marine isotope stages 3–2. Close spacing and good preservation of the 19–17 ka moraines permits 2D glacier reconstruction which suggests the equilibrium line altitude remained depressed by c.1130 m (equivalent to 5.8 ± 0.6 °C colder than present) during this interval. Onset of warming after 17.2 ± 0.2 ka is consistent with climate proxy evidence for a sustained southward shift in the southern westerly winds, which may have promoted deglaciation via shifting ocean currents and promoting increases in atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Long-term shifts in feedbacks among glacier surface change, melt generation and runoff, McMurdo Dry Valleys, Antarctica

Glaciers of the McMurdo dry valleys (MDVs) Antarctica are the main source of streamflow in this polar desert. Because summer air temperatures hover near 0°C small changes in the energy balance strongly affect meltwater generation. In this work, we demonstrate that increased surface roughness, which alters the turbulent transfer of energy between the ice surface and atmosphere, yields a detectable increase in meltwater runoff. At low elevations on the glaciers, basin-like topography became significantly rougher over 13 years between repeat lidar surveys, yielding greater melt. In contrast, the smoother ice at higher elevation exhibited no detectable change in roughness. We pose a conceptual model of the cycle of glacier surface change as a result of climate forcing whereby glacier surfaces transition from being dominated by sublimation to becoming increasingly melt-dominated, which is reversible under prolonged cool periods. This research advances our understanding of warm season effects on polar glaciers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Wintertime Brine Discharge at the Surface of a Cold Polar Glacier and the Unexpected Absence of Associated Seismicity

A subglacial groundwater system beneath Taylor Glacier, Antarctica, discharges hypersaline, iron-rich brine episodically at the glacier surface to create Blood Falls. However, the triggering mechanism for these brine release events is not yet understood. Identifying which fracture processes are observed seismically can help us better characterize the hydrological system at Taylor Glacier, and more generally, provide us with a broader understanding of englacial hydrologic activity in cold glaciers. We document wintertime brine discharge using time-lapse photography. Subfreezing air temperatures during the brine discharge indicate that surface melt-induced hydrofracture is an unlikely trigger for brine release. Further, we analyze local seismic data to test a hypothesis that fracturing generates elevated surface wave energy preceding and/or coinciding with brine release events. Our results show no discernible elevated Rayleigh wave activity prior to or during Blood Falls brine release. Instead, we find a pattern of seismic events dominated by a seasonal signal, with more Rayleigh events occurring in the summer than the winter from the Blood Falls source area. We calculate that the volumetric opening of cracks that would generate Rayleigh waves at our detection limits are of similar size to myriad cracks in glacier ice, lake ice, and frozen sediment in the terminus area. We therefore propose that any fracturing coincident with brine release activity likely consists of a series of smaller opening events that are masked by other seismicity in the local environment.

58 GEOSCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2023

Offshore-onshore record of Last Glacial Maximum–to–present grounding line retreat at Pine Island Glacier, Antarctica

Pine Island Glacier, West Antarctica, is the largest Antarctic contributor to global sea-level rise and is vulnerable to rapid retreat, yet our knowledge of its deglacial history since the Last Glacial Maximum is based largely on marine sediments that record a retreat history ending in the early Holocene. Using a suite of 10 Be exposure ages from onshore glacial deposits directly adjacent to Pine Island Glacier, we show that this major glacier thinned rapidly in the early to mid-Holocene. Our results indicate that Pine Island Glacier was at least 690 m thicker than present prior to ca. 8 ka. We infer that the rapid thinning detected at the site farthest downstream records the arrival and stabilization of the retreating grounding line at that site by 8–6 ka. By combining our exposure ages and the marine record, we extend knowledge of Pine Island Glacier retreat both spatially and temporally: to 50 km from the modern grounding line and to the mid-Holocene, providing a data set that is important for future numerical ice-sheet model validation.

58 GEOSCIENCES↗
DOE OSTI2022

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗