Engineering topics

Brooke Medley

Publications and source records attributed to Brooke Medley.

At least 19 records

NASA NTRS2024

A Physics-Based Antarctic Melt Detection Technique: Combining Advanced Microwave Scanning Radiometer 2, Radiative-Transfer Modeling, and Firn Modeling

Surface melt on ice shelves has been linked to hydrofracture and subsequent ice shelf breakup. Since the 1990s, scientists have been using microwave radiometers to detect melt on ice shelves and ice sheets by applying various statistical thresholding techniques to identify significant increases in brightness temperature that are associated with melt. In this study, instead of using a fixed threshold, we force the Snow Radiative Transfer model (SMRT) with outputs from the Community Firn Model (CFM) to create a dynamic, physics-based threshold for melt. In the process, we also combine our method with statistical thresholding techniques and produce microwave grain size information in the process. We run this “Hybrid Method” across the Larsen C ice shelf as well as thirteen sites on the Antarctic Ice Sheet. Melt and non-melt days from the Hybrid Method and three statistical thresholding techniques match with surface energy balance within 94±1%; the effect of melt on the passive microwave is mostly binary and thus largely detectable by statistical thresholding techniques as well as physics-based techniques. Rather than always replacing statistical thresholding techniques with the Hybrid Method, we recommend using the Hybrid Method in studies where melt volume or grain size is of interest. In this study, we show that the Hybrid Method can be used to (a) model dry snow brightness temperatures of Antarctic snow and (b) derive a measure of grain size; therefore, it is an important step forwards towards using firn and radiative transfer modeling to quantify melt rather than to simply detect melt days.

Marissa E Dattler↗

NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗

NASA NTRS2023

A Wind-Driven Snow Redistribution Module for Alpine3d V3.3.0: Adaptations Designed for Downscaling Ice Sheet Surface Mass Balance

Ice sheets gain mass via snow accumulation at the ice sheet surface, which is the primary component of surface mass balance. On the Antarctic ice sheet, winds redistribute snow resulting in surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 km. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to MERRA-2, a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 mm w.e.yr−1 (13%) and increases variance explained by 24%. Despite these improvements, Alpine3D still underestimates observed accumulation variability, thus providing an opportunity for future model improvement.

Eric Keenan↗

NASA NTRS2023

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗

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 NTRS2022

Simulations of Firn Processes Over the Greenland and Antarctic Ice Sheets: 1980-2021

Conversion of altimetry-derived ice-sheet volume change to mass requires an understanding of the evolution of the combined ice and air content within the firn column. In the absence of suitable techniques to observe the changes to the firn column across the entirety of an ice sheet, the firn column processes are typically modeled. Here, we present new simulations of firn processes over the Greenland and Antarctic ice sheets (GrIS and AIS) using the Community Firn Model and atmospheric reanalysis variables for more than four decades. A data set of more than 250 measured depth–density profiles from both ice sheets provides the basis of the calibration of the dry-snow densification scheme. The resulting scheme results in a reduction in the rate of densification, relative to a commonly used semi-empirical model, through a decreased dependence on the accumulation rate, a proxy for overburden stress. The 1980–2020 modeled firn column runoff, when combined with atmospheric variables from MERRA-2, generates realistic mean integrated surface mass balance values for the Greenland (+390 Gt yr(exp −1)) and Antarctic (+2612 Gt yr(exp −1)) ice sheets when compared to published model-ensemble means. We find that seasonal volume changes associated with firn air content are on average approximately 2.5 times larger than those associated with mass fluxes from surface processes for the AIS and 1.5 times larger for the GrIS; however, when averaged over multiple years, ice and air-volume fluctuations within the firn column are of comparable magnitudes. Between 1996 and 2019, the Greenland Ice Sheet lost nearly 5% of its firn air content, indicating a reduction in the total meltwater retention capability. Nearly all (94 %) of the meltwater produced over the Antarctic Ice Sheet is retained within the firn column through infiltration and refreezing.

Greenland↗

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 NTRS2022

Response of the East Antarctic Ice Sheet to Past and Future Climate Change

The East Antarctic Ice Sheet (EAIS) contains the vast majority of Earth’s glacier ice (~52 metres sea-level equivalent), but is often viewed as less vulnerable to global warming than the West Antarctic or Greenland ice sheets. However, some regions of the EAIS have lost mass over recent decades, prompting the need to re-evaluate its sensitivity to climate change. Here we review the EAIS’s response to past warm periods, synthesise current observations of change, and evaluate future projections. Some marine-based catchments that underwent significant mass loss during past warm periods are currently losing mass, but most projections indicate increased accumulation across the EAIS over the 21st Century, keeping the ice sheet broadly in balance. Beyond 2100, high emissions scenarios generate increased ice discharge and potentially several metres of sea-level rise within just a few centuries, but substantial mass loss could be averted if the Paris Agreement to limit warming below 2°C is satisfied.

Antarctic↗

NASA NTRS2022

Response of the East Antarctic Ice Sheet to Past and Future Climate Change

East Antarctic Ice Sheet↗

NASA NTRS2021

High-Spatial-Resolution Mass Rates From GRACE and GRACE-FO: Global and Ice Sheet Analyses

GRACE and GRACE-FO gravimetry data are unique in its ability to observe global water mass variability and has frequently been applied for the determination of global mass rates to quantify the impacts of climate change and human activity on global freshwater availability and land ice evolution. However, its relatively coarse spatial resolution of ∼300-500 km at monthly intervals precludes its direct use to measure mass changes in smaller drainage basins and individual glacier systems. We present a new global high resolution mascon product that improves spatial resolution and signal recovery by estimating trends directly from more than 16 years of Level 1B GRACE and GRACE-FO data. We compare this product to independent gravimetry spherical harmonic trend estimates and to ice sheet mass rates derived from ICESat and ICESat-2 laser altimetry. We demonstrate enhanced global signal recovery and improved agreement with ICESat/ ICESat-2 and discuss likely causes for the remaining discrepancies between the gravimetry and laser altimetry solutions.

Bryant D. Loomis↗

NASA NTRS2021

The Scientific Legacy of NASA’s Operation IceBridge

The National Aeronautics and Space Administration (NASA)’s Operation IceBridge (OIB) was a 13-year (2009–2021) airborne mission to survey land and sea ice across the Arctic, Antarctic, and Alaska. Here, we review OIB’s goals, instruments, campaigns, key scientific results, and implications for future investigations of the cryosphere. OIB’s primary goal was to use airborne laser altimetry to bridge the gap in fine-resolution elevation measurements of ice from space between the conclusion of NASA’s Ice, Cloud, and land Elevation Satellite (ICESat; 2003–2009) and its follow-on, ICESat-2 (launched 2018). Additional scientific requirements were intended to contextualize observed elevation changes using a multisensor suite of radar sounders, gravimeters, magnetometers, and cameras. Using 15 different aircraft, OIB conducted 968 science flights, of which 42% were repeat surveys of land ice, 42% were surveys of previously unmapped terrain across the Greenland and Antarctic ice sheets, Arctic ice caps, and Alaskan glaciers, and 16% were surveys of sea ice. The combination of an expansive instrument suite and breadth of surveys enabled numerous fundamental advances in our understanding of the Earth’s cryosphere. For land ice, OIB dramatically improved knowledge of interannual outlet-glacier variability, ice-sheet, and outlet-glacier thicknesses, snowfall rates on ice sheets, fjord and sub-ice-shelf bathymetry, and ice-sheet hydrology. Unanticipated discoveries included a reliable method for constraining the thickness within difficult-to-sound incised troughs beneath ice sheets, the extent of the firn aquifer within the Greenland Ice Sheet, the vulnerability of many Greenland and Antarctic outlet glaciers to ocean-driven melting at their grounding zones, and the dominance of surface-melt-driven mass loss of Alaskan glaciers. For sea ice, OIB significantly advanced our understanding of spatiotemporal variability in sea ice freeboard and its snow cover, especially through combined analysis of fine-resolution altimetry, visible imagery, and snow radar measurements of the overlying snow thickness. Such analyses led to the unanticipated discovery of an interdecadal decrease in snow thickness on Arctic sea ice and numerous opportunities to validate sea ice freeboards from satellite radar altimetry. While many of its data sets have yet to be fully explored, OIB’s scientific legacy has already demonstrated the value of sustained investment in reliable airborne platforms, airborne instrument development, interagency and international collaboration, and open and rapid data access to advance our understanding of Earth’s remote polar regions and their role in the Earth system.

Joseph A. MacGregor↗

NASA NTRS2021

Microwave Radiometry at Frequencies From 500 to 1400 MHz: An Emerging Technology for Earth Observations

Microwave radiometry has provided valuable spaceborne observations of Earth's geophysical properties for decades. The recent SMOS, Aquarius, and SMAP satellites have demonstrated the value of measurements at 1400 MHz for observing surface soil moisture, sea surface salinity, sea ice thickness, soil freeze/thaw state, and other geophysical variables. However, the information obtained is limited by penetration through the subsurface at 1400 MHz and by a reduced sensitivity to surface salinity in cold or wind-roughened waters. Recent airborne experiments have shown the potential of brightness temperature measurements from 500–1400 MHz to address these limitations by enabling sensing of soil moisture and sea ice thickness to greater depths, sensing of temperature deep within ice sheets, improved sensing of sea salinity in cold waters, and enhanced sensitivity to soil moisture under vegetation canopies. However, the absence of significant spectrum reserved for passive microwave measurements in the 500–1400 MHz band requires both an opportunistic sensing strategy and systems for reducing the impact of radio-frequency interference. Here, we summarize the potential advantages and applications of 500–1400 MHz microwave radiometry for Earth observation and review recent experiments and demonstrations of these concepts. We also describe the remaining questions and challenges to be addressed in advancing to future spaceborne operation of this technology along with recommendations for future research activities.

microwave radiometry↗

NASA NTRS2021

Physics-Based SNOWPACK Model Improves Representation of Near-Surface Antarctic Snow and Firn Density

Estimates of snow and firn density are required for satellite-altimetry-based retrievals of ice sheet mass balance that rely on volume-to-mass conversions. Therefore, biases and errors in presently used density models confound assessments of ice sheet mass balance and by extension ice sheet contribution to sea level rise. Despite this importance, most contemporary firn densification models rely on simplified semi-empirical methods, which are partially reflected by significant modeled density errors when compared to observations. In this study, we present a new drifting-snow compaction scheme that we have implemented into SNOWPACK, a physics-based land surface snow model. We show that our new scheme improves existing versions of SNOWPACK by increasing simulated near-surface (defined as the top 10 m) density to be more in line with observations (near-surface bias reduction from −44.9 to −5.4 kg m−3). Furthermore, we demonstrate high-quality simulation of near-surface Antarctic snow and firn density at 122 observed density profiles across the Antarctic ice sheet, as indicated by reduced model biases throughout most of the near-surface firn column when compared to two semi-empirical firn densification models (SNOWPACK mean bias=−9.7 kg m−3, IMAU-FDM mean bias=−32.5 kg m−3, GSFC-FDM mean bias=15.5 kg m−3). Notably, our analysis is restricted to the near surface where firn density is most variable due to accumulation and compaction variability driven by synoptic weather and seasonal climate variability. Additionally, the GSFC-FDM exhibits lower mean density bias from 7–10 m (SNOWPACK bias=−22.5 kg m−3, GSFC-FDM bias=10.6 kg m−3) and throughout the entire near surface at high-accumulation sites (SNOWPACK bias=−31.4 kg m−3, GSFC-FDM bias=−4.7 kg m−3). However, we found that the performance of SNOWPACK did not degrade when applied to sites that were not included in the calibration of semi-empirical models. This suggests that SNOWPACK may possibly better represent firn properties in locations without extensive observations and under future climate scenarios, when firn properties are expected to diverge from their present state.

Eric Keenan↗

NASA NTRS2021

Atmospheric River Precipitation Contributed to Rapid Increases in Surface Height of the West Antarctic Ice Sheet in 2019

Estimating the relative contributions of the atmospheric and dynamic components of ice-sheet mass balance is critical for improving projections of future sea level rise. Existing estimates of changes in Antarctic ice-sheet height, which can be used to infer changes in mass, are only accurate at multiyear time scales. However, NASA's Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) laser altimetry mission now allows us to accurately measure changes in ice-sheet height at subannual time scales. Here, we use ICESat-2 data to estimate height changes over Antarctica between April 2019 and June 2020. These data show widespread increases in surface height over West Antarctica during the 2019 austral winter. Using climate reanalysis data, we show that 41% of increases in height during winter were from snow accumulation via extreme precipitation events—63% of these events were associated with landfalling atmospheric rivers (ARs) which occurred only 5.1% of the time.

Susheel Adusumilli↗

NASA NTRS2020

Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves

Ocean-driven basal melting of Antarctica’s floating ice shelves accounts for about half of their mass loss in steady state, where gains in ice-shelf mass are balanced by losses. Ice-shelf thickness changes driven by varying basal melt rates modulate mass loss from the grounded ice sheet and its contribution to sea level, and the changing meltwater fluxes influence climate processes in the Southern Ocean. Existing continent-wide melt-rate datasets have no temporal variability, introducing uncertainties in sea level and climate projections. Here, we combine surface height data from satellite radar altimeters with satellite-derived ice velocities and a new model of firn-layer evolution to generate a high-resolution map of time-averaged (2010–2018) basal melt rates and time series (1994–2018) of meltwater fluxes for most ice shelves. Total basal meltwater flux in 1994 (1,090 ± 150 Gt/yr) was similar to the steady-state value (1,100 ± 60 Gt/yr), but increased to 1,570 ± 140 Gt/yr in 2009, followed by a decline to 1,160 ± 150 Gt/yr in 2018. For the four largest ‘cold-water’ ice shelves, we partition meltwater fluxes into deep and shallow sources to reveal distinct signatures of temporal variability, providing insights into climate forcing of basal melting and the impact of this melting on the Southern Ocean.

Susheel Adusumilli↗

NASA NTRS2020

Pervasive ice sheet mass loss reflects competing ocean and atmosphere processes

Quantifying changes in Earth’s ice sheets and identifying the climate drivers are central to improving sea level projections. We provide unified estimates of grounded and floating ice mass change from 2003 to 2019 using NASA’s Ice, Cloud and land Elevation Satellite (ICESat) and ICESat-2 satellite laser altimetry. Our data reveal patterns likely linked to competing climate processes: Ice loss from coastal Greenland (increased surface melt), Antarctic ice shelves (increased ocean melting), and Greenland and Antarctic outlet glaciers (dynamic response to ocean melting) was partially compensated by mass gains over ice sheet interiors (increased snow accumulation). Losses outpaced gains, with grounded-ice loss from Greenland (200 billion tonnes per year) and Antarctica (118 billion tonnes per year) contributing 14 millimeters to sea level. Mass lost from West Antarctica’s ice shelves accounted for more than 30% of that region’s total.

Ben Smith↗

NASA NTRS2019

A New Regional Climate Model for POLAR‐CORDEX: Evaluation of a 30‐Year Hindcast with COSMO‐CLM^2 Over Antarctica

Continent‐wide climate information over the Antarctic Ice Sheet (AIS) is important to obtain accurate information of present climate and reduce uncertainties of the ice sheet mass balance response and resulting global sea level rise to future climate change. In this study, the COSMO‐CLM^2 Regional Climate Model is applied over the AIS and adapted for the specific meteorological and climatological conditions of the region. A 30‐year hindcast was performed and evaluated against observational records consisting of long‐term ground‐based meteorological observations, automatic weather stations, radiosoundings, satellite records, stake measurements and ice cores. Reasonable agreement regarding the surface and upper‐air climate is achieved by the COSMO‐CLM^2 model, comparable to the performance of other state‐of‐the‐art climate models over the AIS. Meteorological variability of the surface climate is adequately simulated, and biases in the radiation and surface mass balance are small. The presented model therefore contributes as a new member to the COordinated Regional Downscaling EXperiment project over the AIS (POLAR‐CORDEX) and the CORDEX‐CORE initiative.

Antarctica↗

NASA NTRSDate not supplied

PolarMERRA: A collaborative effort for NASA’s polar reanalysis

The polar atmosphere is often poorly represented in global models, adversely impacting analyses in these regions and forecasts on much larger scales. PolarMERRA is an initiative that brings together the modeling expertise from NASA’s global modeling and assimilation office (GMAO) and observational expertise from NASA’s Cryospheric sciences lab to assess and improve the representation of polar processes in the NASA’s modelling system in preparation for upcoming and next generation reanalyses. Model improvements in resolution, physics, and data assimilation have been guided by assessment of model skill in the polar regions. Preliminary results from a novel validation framework indicate a reduction in near surface air temperature, turbulent fluxes, and precipitation biases.

Chelsea Parker↗