Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ice sheet”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2): Science Requirements, Concept, and Implementation

The Ice, Cloud, and land Elevation Satellite (ICESat) mission used laser altimetry measurements to determine changes in elevations of glaciers and ice sheets, as well as sea ice thickness distribution. These measurements have provided important information on the response of the cryosphere (Earths frozen surfaces) to changes in atmosphere and ocean condition. ICESat operated from 2003-2009 and provided repeat altimetry measurements not only to the cryosphere scientific community but also to the ocean, terrestrial and atmospheric scientific communities. The conclusive assessment of significant ongoing rapid changes in the Earths ice cover, in part supported by ICESat observations, has strengthened the need for sustained, high accuracy, repeat observations similar to what was provided by the ICESat mission. Following recommendations from the National Research Council for an ICESat follow-on mission, the ICESat-2 mission is now under development for planned launch in 2018. The primary scientific aims of the ICESat-2 mission are to continue measurements of sea ice freeboard and ice sheet elevation to determine their changes at scales from outlet glaciers to the entire ice sheet, and from 10s of meters to the entire polar oceans for sea ice freeboard. ICESat carried a single beam profiling laser altimeter that produced approximately 70 m diameter footprints on the surface of the Earth at approximately 150 m along-track intervals. In contrast, ICESat-2 will operate with three pairs of beams, each pair separated by about 3 km across-track with a pair spacing of 90 m. Each of the beams will have a nominal 17 m diameter footprint with an along-track sampling interval of 0.7 m. The differences in the ICESat-2 measurement concept are a result of overcoming some limitations associated with the approach used in the ICESat mission. The beam pair configuration of ICESat-2 allows for the determination of local cross-track slope, a significant factor in measuring elevation change for the outlet glaciers surrounding the Greenland and Antarctica coasts. The multiple beam pairs also provide improved spatial coverage. The dense spatial sampling eliminates along-track measurement gaps, and the small footprint diameter is especially useful for sea surface height measurements in the often narrow leads needed for sea ice freeboard and ice thickness retrievals. The ICESat-2 instrumentation concept uses a low energy 532 nm (green) laser in conjunction with single-photon sensitive detectors to measure range. Combining ICESat-2 data with altimetry data collected since the start of the ICESat mission in 2003, such as Operation IceBridge and ESAs CryoSat-2, will yield a 15+ year record of changes in ice sheet elevation and sea ice thickness. ICESat-2 will also provide information of mountain glacier and ice cap elevations changes, land and vegetation heights, inland water elevations, sea surface heights, and cloud layering and optical thickness.

ICESat-2↗
NASA NTRS2016

Ice Melt, Sea Level Rise and Superstorms: Evidence from Paleoclimate Data, Climate Modeling, and Modern Observations that 2C Global Warming Could Be Dangerous

We use numerical climate simulations, paleoclimate data, and modern observations to study the effect of growing ice melt from Antarctica and Greenland. Meltwater tends to stabilize the ocean column, inducing amplifying feedbacks that increase subsurface ocean warming and ice shelf melting. Cold meltwater and induced dynamical effects cause ocean surface cooling in the Southern Ocean and North Atlantic, thus increasing Earth's energy imbalance and heat flux into most of the global ocean's surface. Southern Ocean surface cooling, while lower latitudes are warming, increases precipitation on the Southern Ocean, increasing ocean stratification, slowing deepwater formation, and increasing ice sheet mass loss. These feedbacks make ice sheets in contact with the ocean vulnerable to accelerating disintegration. We hypothesize that ice mass loss from the most vulnerable ice, sufficient to raise sea level several meters, is better approximated as exponential than by a more linear response. Doubling times of 10, 20 or 40 years yield multi-meter sea level rise in about 50, 100 or 200 years. Recent ice melt doubling times are near the lower end of the 10-40-year range, but the record is too short to confirm the nature of the response. The feedbacks, including subsurface ocean warming, help explain paleoclimate data and point to a dominant Southern Ocean role in controlling atmospheric CO2, which in turn exercised tight control on global temperature and sea level. The millennial (500-2000-year) timescale of deep-ocean ventilation affects the timescale for natural CO2 change and thus the timescale for paleo-global climate, ice sheet, and sea level changes, but this paleo-millennial timescale should not be misinterpreted as the timescale for ice sheet response to a rapid, large, human-made climate forcing. These climate feedbacks aid interpretation of events late in the prior interglacial, when sea level rose to C6-9m with evidence of extreme storms while Earth was less than 1 C warmer than today. Ice melt cooling of the North Atlantic and Southern oceans increases atmospheric temperature gradients, eddy kinetic energy and baroclinicity, thus driving more powerful storms. The modeling, paleoclimate evidence, and ongoing observations together imply that 2 C global warming above the preindustrial level could be dangerous. Continued high fossil fuel emissions this century are predicted to yield (1) cooling of the Southern Ocean, especially in the Western Hemisphere; (2) slowing of the Southern Ocean overturning circulation, warming of the ice shelves, and growing ice sheet mass loss; (3) slowdown and eventual shutdown of the Atlantic overturning circulation with cooling of the North Atlantic region; (4) increasingly powerful storms; and (5) nonlinearly growing sea level rise, reaching several meters over a timescale of 50-150 years. These predictions, especially the cooling in the Southern Ocean and North Atlantic with markedly reduced warming or even cooling in Europe, differ fundamentally from existing climate change assessments. We discuss observations and modeling studies needed to refute or clarify these assertions.

Computer simulation↗
NASA NTRS2001

Balance Mass Flux and Velocity Across the Equilibrium Line in Ice Drainage Systems of Greenland

Estimates of balance mass flux and the depth-averaged ice velocity through the cross-section aligned with the equilibrium line are produced for each of six drainage systems in Greenland. (The equilibrium line, which lies at approximately 1200 m elevation on the ice sheet, is the boundary between the area of net snow accumulation at higher elevations and the areas of net melting at lower elevations around the ice sheet.) Ice drainage divides and six major drainage systems are delineated using surface topography from ERS (European Remote Sensing) radar altimeter data. The net accumulation rate in the accumulation zone bounded by the equilibrium line is 399 Gt/yr and net ablation rate in the remaining area is 231 Gt/yr. (1 GigaTon of ice is 1090 kM(exp 3). The mean balance mass flux and depth-averaged ice velocity at the cross-section aligned with the modeled equilibrium line are 0.1011 Gt kM(exp -2)/yr and 0.111 km/yr, respectively, with little variation in these values from system to system. The ratio of the ice mass above the equilibrium line to the rate of mass output implies an effective exchange time of approximately 6000 years for total mass exchange. The range of exchange times, from a low of 3 ka in the SE drainage system to 14 ka in the NE, suggests a rank as to which regions of the ice sheet may respond more rapidly to climate fluctuations.

Zwally, H. Jay↗
NASA NTRS2023

Evaluating Greenland Surface-Mass-Balance and Firn-Densification Data Using ICESat-2 Altimetry

Surface-mass-balance (SMB) and firn-densification (FD) models are widely used in altimetry studies as a tool to separate atmospheric-driven from ice-dynamics-driven ice-sheet mass changes and to partition observed volume changes into ice-mass changes and firn-air-content changes. Until now, SMB models have been principally validated based on comparison with ice core and weather station data or comparison with widely separated flight radar-survey flight lines. Firn-densification models have been primarily validated based on their ability to match net densification over decades, as recorded in firn cores, and the short-term time-dependent component of densification has rarely been evaluated at all. The advent of systematic ice-sheet-wide repeated ice-surface-height measurements from ICESat-2 (the Ice Cloud, and land Elevation Satellite, 2) allows us to measure the net surface-height change of the Greenland ice sheet at quarterly resolution and compare the measured surface-height differences directly with those predicted by three FD–SMB models: MARv3.5.11 (Modèle Atmosphérique Régional version 3.5.11) and GSFCv1.1 and GSFCv1.2 (the Goddard Space Flight Center FD–SMB models version 1.1 and 1.2). By segregating the data by season and elevation, and based on the timing and magnitude of modelled processes in areas where we expect minimal ice-dynamics-driven height changes, we investigate the models' accuracy in predicting atmospherically driven height changes. We find that while all three models do well in predicting the large seasonal changes in the low-elevation parts of the ice sheet where melt rates are highest, two of the models (MARv3.5.11 and GSFCv1.1) systematically overpredict, by around a factor of 2, the magnitude of height changes in the high-elevation parts of the ice sheet, particularly those associated with melt events. This overprediction seems to be associated with the melt sensitivity of the models in the high-elevation part of the ice sheet. The third model, GSFCv1.2, which has an updated high-elevation melt parameterization, avoids this overprediction.

Benjamin E. Smith↗
NASA NTRS1991

Ice cores and SeaRISE: What we do (and don't) know

Ice core analyses are needed in SeaRISE to learn what the West Antarctic ice sheet and other marine ice sheets were like in the past, what climate changes led to their present states, and how they behave. The major results of interest to SeaRISE from previous ice core analyses in West Antarctic are that the end of the last ice age caused temperature and accumulation rate increases in inland regions, leading to ice sheet thickening followed by thinning to the present.

Alley, Richard B.↗
NASA NTRS1986

A scanning radar altimeter for mapping continental topography

Topographic information constitutes a fundamental data set for the Earth sciences. In the geological and geophysical sciences, topography combined with gravitational information provides an important constraint on the structure and rheologic properties of the crust and lithosphere. Detailed topography data can also be used to map offsets associated with faulting and to reveal the effects of tectonic deformation. In the polar regions, elevation data form a crucial but as yet largely unavailable resource for studying ice sheet mass balance and ice flow dynamics. The vast Antarctic ice sheet is the largest fresh water reservoir on Earth and is an important influence on ocean circulation and global climate. However, our knowledge of its stability is so limited that we cannot even specify whether the Antarctic ice sheet is growing or shrinking. It is clear that there is need for high quality global topography data. A summary of potential applications with their resolution requirements is shown.

Dixon, T. H.↗
NASA NTRS2006

Ice-Borehole Probe

An instrumentation system has been developed for studying interactions between a glacier or ice sheet and the underlying rock and/or soil. Prior borehole imaging systems have been used in well-drilling and mineral-exploration applications and for studying relatively thin valley glaciers, but have not been used for studying thick ice sheets like those of Antarctica. The system includes a cylindrical imaging probe that is lowered into a hole that has been bored through the ice to the ice/bedrock interface by use of an established hot-water-jet technique. The images acquired by the cameras yield information on the movement of the ice relative to the bedrock and on visible features of the lower structure of the ice sheet, including ice layers formed at different times, bubbles, and mineralogical inclusions. At the time of reporting the information for this article, the system was just deployed in two boreholes on the Amery ice shelf in East Antarctica and after successful 2000 2001 deployments in 4 boreholes at Ice Stream C, West Antarctica, and in 2002 at Black Rapids Glacier, Alaska. The probe is designed to operate at temperatures from 40 to +40 C and to withstand the cold, wet, high-pressure [130-atm (13.20-MPa)] environment at the bottom of a water-filled borehole in ice as deep as 1.6 km. A current version is being outfitted to service 2.4-km-deep boreholes at the Rutford Ice Stream in West Antarctica. The probe (see figure) contains a sidelooking charge-coupled-device (CCD) camera that generates both a real-time analog video signal and a sequence of still-image data, and contains a digital videotape recorder. The probe also contains a downward-looking CCD analog video camera, plus halogen lamps to illuminate the fields of view of both cameras. The analog video outputs of the cameras are converted to optical signals that are transmitted to a surface station via optical fibers in a cable. Electric power is supplied to the probe through wires in the cable at a potential of 170 VDC. A DC-to-DC converter steps the supply down to 12 VDC for the lights, cameras, and image-data-transmission circuitry. Heat generated by dissipation of electric power in the probe is removed simply by conduction through the probe housing to the visible features of the lower structure of the ice sheet, including ice layers formed at different times, bubbles, and mineralogical inclusions. At the time of reporting the information for this article, the system was just deployed in two boreholes on the Amery ice shelf in East Antarctica and after successful 2000 2001 deployments in 4 boreholes at Ice Stream C, West Antarctica, and in 2002 at Black Rapids Glacier, Alaska. The probe is designed to operate at temperatures from 40 to +40 C and to withstand the cold, wet, high-pressure [130-atm (13.20-MPa)] environment at the bottom of a water-filled borehole in ice as deep as 1.6 km. A current version is being outfitted to service 2.4-km-deep boreholes at the Rutford Ice Stream in West Antarctica. The probe (see figure) contains a sidelooking charge-coupled-device (CCD) camera that generates both a real-time analog video signal and a sequence of still-image data, and contains a digital videotape recorder. The probe also contains a downward-looking CCD analog video camera, plus halogen lamps to illuminate the fields of view of both cameras. The analog video outputs of the cameras are converted to optical signals that are transmitted to a surface station via optical fibers in a cable. Electric power is supplied to the probe through wires in the cable at a potential of 170 VDC. A DC-to-DC converter steps the supply down to 12 VDC for the lights, cameras, and image-datatransmission circuitry. Heat generated by dissipation of electric power in the probe is removed simply by conduction through the probe housing to the visible features of the lower structure of the ice sheet, including ice layers formed at different times, bubbles, and mineralogical inclusions. At thime of reporting the information for this article, the system was just deployed in two boreholes on the Amery ice shelf in East Antarctica and after successful 2000 2001 deployments in 4 boreholes at Ice Stream C, West Antarctica, and in 2002 at Black Rapids Glacier, Alaska. The probe is designed to operate at temperatures from 40 to +40 C and to withstand the cold, wet, high-pressure [130-atm (13.20-MPa)] environment at the bottom of a water-filled borehole in ice as deep as 1.6 km. A current version is being outfitted to service 2.4-km-deep boreholes at the Rutford Ice Stream in West Antarctica. The probe (see figure) contains a sidelooking charge-coupled-device (CCD) camera that generates both a real-time analog video signal and a sequence of still-image data, and contains a digital videotape recorder. The probe also contains a downward-looking CCD analog video camera, plus halogen lamps to illuminate the fields of view of both cameras. The analog video outputs of the cameras are converted to optical signals that are transmitted to a surface station via optical fibers in a cable. Electric power is supplied to the probe through wires in the cable at a potential of 170 VDC. A DC-to-DC converter steps the supply down to 12 VDC for the lights, cameras, and image-datatransmission circuitry. Heat generated by dissipation of electric power in the probe is removed simply by conduction through the probe housing to the adjacent water and ice.

Behar, Alberto↗
DOE OSTI2024

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↗
NASA NTRS2007

Remote Sensing of Terrestrial Snow and Ice for Global Change Studies

Snow and ice play a significant role in the Earth's water cycle and are sensitive and informative indicators climate change. Significant changes in terrestrial snow and ice water storage are forecast, and while evidence of large-scale changes is emerging, in situ measurements alone are insufficient to help us understand and explain these changes. Imaging remote sensing systems are capable of successfully observing snow and ice in the cryosphere. This chapter examines how those remote sensing sensors, that now have more than 35 years of observation records, are capable of providing information about snow cover, snow water equivalent, snow melt, ice sheet temperature and ice sheet albedo. While significant progress has been made, especially in the last five years, a better understanding is required of the records of satellite observations of these cryospheric variables.

Kelly, Richard↗
NASA NTRS2011

Assessment of Cloud Screening with Apparent Surface Reflectance in Support of the ICESat-2 Mission

The separation of cloud and clear scenes is usually one of the first steps in satellite data analysis. Before deriving a geophysical product, almost every satellite mission requires a cloud mask to label a scene as either clear or cloudy through a cloud detection procedure. For clear scenes, products such as surface properties may be retrieved; for cloudy scenes, scientist can focus on studying the cloud properties. Hence the quality of cloud detection directly affects the quality of most satellite operational and research products. This is certainly true for the Ice, Cloud, and land Elevation Satellite-2 (lCESat-2), which is the successor to the ICESat-l. As a top priority mission, ICESat-2 will continue to provide measurements of ice sheets and sea ice elevation on a global scale. Studies have shown that clouds can significantly affect the accuracy of the retrieved results. For example, some of the photons (a photon is a basic unit of light) in the laser beam will be scattered by cloud particles on its way. So instead of traveling in a straight line, these photons are scattered sideways and have traveled a longer path. This will result in biases in ice sheet elevation measurements. Hence cloud screening must be done and be done accurately before the retrievals.

Yang, Yuekui↗
NASA NTRS2013

Airborne Tomographic Swath Ice Sounding Processing System

Glaciers and ice sheets modulate global sea level by storing water deposited as snow on the surface, and discharging water back into the ocean through melting. Their physical state can be characterized in terms of their mass balance and dynamics. To estimate the current ice mass balance, and to predict future changes in the motion of the Greenland and Antarctic ice sheets, it is necessary to know the ice sheet thickness and the physical conditions of the ice sheet surface and bed. This information is required at fine resolution and over extensive portions of the ice sheets. A tomographic algorithm has been developed to take raw data collected by a multiple-channel synthetic aperture sounding radar system over a polar ice sheet and convert those data into two-dimensional (2D) ice thickness measurements. Prior to this work, conventional processing techniques only provided one-dimensional ice thickness measurements along profiles.

Wu, Xiaoqing↗
DOE OSTI2022

The collapse of the Cordilleran–Laurentide ice saddle and early opening of the Mackenzie Valley, Northwest Territories, Canada, constrained by 10 Be exposure dating

Deglaciation of the northwestern Laurentide Ice Sheet in the central Mackenzie Valley opened the northern portion of the deglacial Ice-Free Corridor between the Laurentide and Cordilleran ice sheets and a drainage route to the Arctic Ocean. In addition, ice sheet saddle collapse in this section of the Laurentide Ice Sheet has been implicated as a mechanism for delivering substantial freshwater influx into the Arctic Ocean on centennial timescales. However, there is little empirical data to constrain the deglaciation chronology in the central Mackenzie Valley where the northern slopes of the ice saddle were located. Here, we present 30 new 10Be cosmogenic nuclide exposure dates across six sites, including two elevation transects, which constrain the timing and rate of thinning and retreat of the Laurentide Ice Sheet in the area. Our new 10Be dates indicate that the initial deglaciation of the eastern summits of the central Mackenzie Mountains began at ~15.8 ka (17.1–14.6 ka), ~1000 years earlier than in previous reconstructions. The main phase of ice saddle collapse occurred between ~14.9 and 13.6 ka, consistent with numerical modelling simulations, placing this event within the Bølling–Allerød interval (14.6–12.9 ka). Our new dates require a revision of ice margin retreat dynamics, with ice retreating more easterly rather than southward along the Mackenzie Valley. In addition, we quantify a total sea level rise contribution from the Cordilleran–Laurentide ice saddle region of ~11.2 m between 16 and 13 ka.

58 GEOSCIENCES↗
NASA NTRS1984

Observing the polar regions from space

The use of active and passive microwave sensors to provide information on the ice cover in the earth's polar regions is described, along with the ice characteristics that can be measured by such techniques. The continental ice sheets of Greenland and Antarctica represent vast reservoirs of fresh water that can significantly impact sea level if the ice sheets are changing in size. Satellite remote sensing gives information on: sea ice extent and physical characteristics; detailed images of ice flows and open-water leads within the ice pack; sea-ice movement; zones of summer melting and snow-accumulation rates on the continental ice sheets; accurate estimates of ice-surface elevation and detection of zones on the ice sheet that are either thickening or thinning; accurate, all-weather mapping of ice coastlines and large crevasses, and estimates of ice discharge rates from the ice sheets.

Thomas, R. H.↗
NASA NTRS2012

Measuring the Surface Temperature of the Cryosphere using Remote Sensing

A general description of the remote sensing of cryosphere surface temperatures from satellites will be provided. This will give historical information on surface-temperature measurements from space. There will also be a detailed description of measuring the surface temperature of the Greenland Ice Sheet using Moderate-Resolution Imaging Spectroradiometer (MODIS) data which will be the focus of the presentation. Enhanced melting of the Greenland Ice Sheet has been documented in recent literature along with surface-temperature increases measured using infrared satellite data since 1981. Using a recently-developed climate data record, trends in the clear-sky ice-surface temperature (IST) of the Greenland Ice Sheet have been studied using the MODIS IST product. Daily and monthly MODIS ISTs of the Greenland Ice Sheet beginning on 1 March 2000 and continuing through 31 December 2010 are now freely available to download at 6.25-km spatial resolution on a polar stereographic grid. Maps showing the maximum extent of melt for the entire ice sheet and for the six major drainage basins have been developed from the MODIS IST dataset. Twelve-year trends of the duration of the melt season on the ice sheet vary in different drainage basins with some basins melting progressively earlier over the course of the study period. Some (but not all) of the basins also show a progressively-longer duration of melt. The consistency of this IST record, with temperature and melt records from other sources will be discussed.

Hall, Dorothy K.↗
DOE OSTI2025

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↗
DOE OSTI2022

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

The Impact of Firn Models on Ultrawideband Brightness Temperatures in the Partially Coherent Model

The density profile of a polar ice sheet is an important parameter for the estimation of ice mass balance. Wave reflections caused by density variations are also a key uncertainty in the retrieval of ice sheet temperature profiles in Ultra-Wide band radiometry. In this paper, we examine different firn density profile models and analyze the subsurface reflections they cause using an analytical partially coherent approach. We also examine firn density profiles obtained from borehole measurements, from past UWBRAD modeling studies, from a community firn model, and from snow radar echo measurements. In previous studies, the ice sheet has been model as a 1D random medium with density variations in depth. However, horizontal density variations also exist, so that the ice sheet is a 3D random medium. Analyses using the partially coherent model show that in the presence of horizontal fluctuations, contributions from short scale variations vanish as the horizontal correlation length decreases due to the diffraction of waves.

Haokui Xu↗
NASA NTRS2019

A Physical Model of Moulin Formation and Evolution

Nearly all proglacial water discharge from the present-day Greenland Ice Sheet is routed englacially via moulins. Identification of these moulins in high-resolution imagery is a frequent topic of study, but the processes controlling how and where moulins form, including on past ice sheets for which remote-sensing data are not available, remain poorly understood. Because moulins may reasonably compose approximately 10-15% of the englacial-subglacial hydrologic system, the evolution and shape of moulins can alter the timing of meltwater inputs to the bed. This evolution can impact both the form of the subglacial hydrologic system and the structure of associated geomorphological structures. Here, we develop a physical model of moulin formation and evolution to constrain the role of englacial processes in controlling the form and structure of the subglacial hydrologic system. Ice deformation within and around a moulin is both viscous and elastic, with the rate of turbulent and heat dissipation from water circulation in the moulin controlling both moulin wall melting and warming of the surrounding ice. We find moulin geometry is responsive to changes in these parameters over hours to days, indicating that diurnal and multi-day variations in surface melt can substantially alter the geometry of a moulin and the pressure-discharge relationship at the bed of the ice sheet. These results should be considered carefully when determining surface water inputs for subglacial hydrologic models. In the future, a parameter space study of these results will be combined with an analytic model to create a predictive, stochastic model of moulin and crevasse locations. This future model will be applicable to constraining the potential for surface-to-bed connections in regions where the exact ice-sheet surface morphology is not known, including ice sheets under future warming atmospheric conditions, and paleo ice sheets, where moulins created modern landforms.

Andrews, Lauren C.↗