Engineering Papers⌕ Search

Engineering topics

Neumann, Thomas A.

Publications and source records attributed to Neumann, Thomas A..

Antarctic Ice Shelf Thickness Change from Multimission Lidar Mapping

We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.

Sutterley, Tyler C.↗

Mabel Engineering Flights, 2010-2013: Flight Report

In December 2010, NASA deployed for the first time the Multiple Altimeter Beam Experimental Lidar (MABEL), an airborne simulator for Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) algorithm development. Between 2010 and 2013, engineering flights were conducted in the continental United States, to ready the instrument for deployments to Alaska and Iceland, where flight lines could be designed and flown over sea ice and grounded ice. Ultimately, MABEL engineering missions included: 1) flights based out of NASA Armstrong Flight Research Center (California, formerly Dryden Flight Research Center) in 2010, 2011, and 2012; flights based out of NASA Wallops Flight Facility (Virginia) in 2012; flights based out of NASA Langley Research Center (Virginia) in 2013; and flights based out of the Mojave Air and Space Port (California) in 2013.

Brunt, Kelly M.↗

Assessment of Altimetry Using Ground-Based GPS Data from the 88S Traverse, Antarctica, in Support of ICESat-2

We conducted a 750 km kinematic GPS survey, referred to as the 88S Traverse, based out of South Pole Station, Antarctica, between December 2017 and January 2018. This ground-based survey was designed to validate space-borne altimetry and airborne altimetry developed at NASA. The 88S Traverse intersects 20% of the ICESat-2 satellite orbits on a route that has been flown by two different Operation IceBridge airborne laser altimeters: the Airborne Topographic Mapper (ATM; 26 October 2014) and the University of Alaska Fairbanks (UAF) Lidar (30 November and 3 December 2017). Here we present an overview of the ground-based GPS data quality and a quantitative assessment of the airborne laser altimetry over a flat section of the ice sheet interior. Results indicate that the GPS data are internally consistent (1:1+- 4:1 cm). Relative to the ground-based 88S Traverse data, the elevation biases for ATM and the UAF lidar range from -9:5 to 3.6 cm, while surface measurement precisions are equal to or better than 14.1 cm. These results suggest that the ground-based GPS data and airborne altimetry data are appropriate for the validation of ICESat-2 surface elevation data.

Brunt, Kelly M.↗

MABEL Iceland 2012 Flight Report

In March and April 2012, NASA conducted an airborne lidar campaign based out of Keflavik, Iceland, in support of Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) algorithm development. The survey targeted the Greenland Ice Sheet, Iceland ice caps, and sea ice in the Arctic Ocean during the winter season. Ultimately, the mission, MABEL Iceland 2012, including checkout and transit flights, conducted 14 science flights, for a total of over 80 flight hours over glaciers, icefields, and sea ice.

MABEL↗

Assessment of NASA Airborne Laser Altimetry Data Using Ground-Based GPS Data near Summit Station, Greenland

A series of NASA airborne lidars have been used in support of satellite laser altimetry missions. These airbornelaser altimeters have been deployed for satellite instrument development, for spaceborne data validation, and to bridge the data gap between satellite missions. We used data from ground-based Global Positioning System (GPS) surveys of an 11 km long track near Summit Station, Greenland, to assess the surface elevation bias and measurement precision of three airborne laser altimeters including the Airborne Topographic Mapper (ATM), the Land, Vegetation, and Ice Sensor (LVIS), and the Multiple Altimeter Beam Experimental Lidar (MABEL). Ground-based GPS data from the monthly ground-based traverses, which commenced in 2006, allowed for the assessment of nine airborne lidar surveys associated with ATM and LVIS between 2007 and 2016. Surface elevation biases for these altimeters over the flat, ice-sheet interior are less than 0.12 m, while assessments of measurement precision are 0.09 m or better. Ground-based GPS positions determined both with and without differential post-processing techniques provided internally consistent solutions. Results from the analyses of ground-based and airborne data provide validation strategy guidance for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) elevation and elevation-change data products.

Summit Station↗

Greenland Subglacial Drainage Evolution Regulated by Weakly Connected Regions of the Bed

Penetration of surface meltwater to the bed of the Greenland Ice Sheet each summer causes an initial increase in ice speed due to elevated basal water pressure, followed by slowdown in late summer that continues into fall and winter. While this seasonal pattern is commonly explained by an evolution of the subglacial drainage system from an inefficient distributed to efficient channelized configuration, mounting evidence indicates that subglacial channels are unable to explain important aspects of hydrodynamic coupling in late summer and fall. Here we use numerical models of subglacial drainage and ice flow to show that limited, gradual leakage of water and lowering of water pressure in weakly connected regions of the bed can explain the dominant features in late and post melt season ice dynamics. These results suggest that a third weakly connected drainage component should be included in the conceptual model of subglacial hydrology.

Greenland Ice Sheet↗

Holocene Accumulation and Ice Flow near the West Antarctic Ice Sheet Divide Ice Core Site

The West Antarctic Ice Sheet Divide Core (WDC) provided a high-resolution climate record from near the Ross-Amundsen Divide in Central West Antarctica. In addition, radar-detected internal layers in the vicinity of the WDC site have been dated directly from the ice core to provide spatial variations in the age structure of the region. Using these two data sets together, we first infer a high-resolution Holocene accumulation-rate history from 9.2 thousand years of the ice-core timescale and then confirm that this climate history is consistent with internal layers upstream of the core site. Even though the WDC was drilled only 24 kilometers from the modern ice divide, advection of ice from upstream must be taken into account. We evaluate histories of accumulation rate by using a flowband model to generate internal layers that we compare to observed layers. Results show that the centennially averaged accumulation rate was over 20 percent lower than modern at 9.2 thousand years before present (B.P.), increased by 40 percent from 9.2 to 2.3 thousand years B.P., and decreased by at least 10 percent over the past 2 thousand years B.P. to the modern values; these Holocene accumulation-rate changes in Central West Antarctica are larger than changes inferred from East Antarctic ice-core records. Despite significant changes in accumulation rate, throughout the Holocene the regional accumulation pattern has likely remained similar to today, and the ice-divide position has likely remained on average within 5 kilometers of its modern position. Continent-scale ice-sheet models used for reconstructions of West Antarctic ice volume should incorporate this accumulation history.

ice↗

SIMPL/AVIRIS-NG Greenland 2015: Flight Report

In August 2015, NASA conducted a two-­aircraft, coordinated campaign based out of Thule Air Base, Greenland, in support of Ice, Cloud, and land Elevation Satellite‐2 (ICESat-2) algorithm development. The survey targeted the Greenland Ice Sheet and sea ice in the Arctic Ocean during the summer melt season. The survey was conducted with a photon-counting laser altimeter in one aircraft and an imaging spectrometer in the second aircraft. Ultimately, the mission, SIMPL/AVIRIS-NG Greenland 2015, conducted nine coordinated science flights, for a total of 37 flight hours over the ice sheet and sea ice.

Ice Sheet↗

Determination of Local Slope on the Greenland Ice Sheet Using a Multibeam Photon-Counting Lidar in Preparation for the ICESat-2 Mission

The greatest changes in elevation in Greenland and Antarctica are happening along the margins of the ice sheets where the surface frequently has significant slopes. For this reason, the upcoming Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission utilizes pairs of laser altimeter beams that are perpendicular to the flight direction in order to extract slope information in addition to elevation. The Multiple Altimeter Beam Experimental Lidar (MABEL) is a high-altitude airborne laser altimeter designed as a simulator for ICESat-2. The MABEL design uses multiple beams at fixed angles and allows for local slope determination. Here, we present local slopes as determined by MABEL and compare them to those determined by the Airborne Topographic Mapper (ATM) over the same flight lines in Greenland. We make these comparisons with consideration for the planned ICESat-2 beam geometry. Results indicate that the mean slope residuals between MABEL and ATM remain small (< 0.05◦) through a wide range of localized slopes using ICESat-2 beam geometry. Furthermore, when MABEL data are subsampled by a factor of 4 to mimic the planned ICESat-2 transmit-energy configuration, the results are indistinguishable from the full-data-rate analysis. Results from MABEL suggest that ICESat-2 beam geometry and transmit-energy configuration are appropriate for the determination of slope on ∼90-m spatial scales, a measurement that will be fundamental to deconvolving the effects of surface slope from the ice-sheet surface change derived from ICESat-2.

Greenland↗

Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation for the ICESat-2 Mission

NASA's Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission is a decadal survey mission (2016 launch). The mission objectives are to measure land ice elevation, sea ice freeboard, and changes in these variables, as well as to collect measurements over vegetation to facilitate canopy height determination. Two innovative components will characterize the ICESat-2 lidar: 1) collection of elevation data by a multibeam system and 2) application of micropulse lidar (photon-counting) technology. A photon-counting altimeter yields clouds of discrete points, resulting from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of the returned points to reflectors of interest. The objective of this paper is to derive an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2 data, based on airborne observations with a Sigma Space micropulse lidar. The mathematical algorithm uses spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors, and geostatistical classification parameters and hyperparameters. Validation shows that ground and canopy elevation, and hence canopy height, can be expected to be observable with high accuracy by ICESat-2 for all expected beam energies considered for instrument design (93.01%-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp), and 72.85%-98.68% for 0.48 msp). The algorithm derived here is generally applicable for elevation determination from photoncounting lidar altimeter data collected over forested areas, land ice, sea ice, and land surfaces, as well as for cloud detection.

Algorithms↗

Constraining Landscape History and Glacial Erosivity Using Paired Cosmogenic Nuclides in Upernavik, Northwest Greenland

High-latitude landscape evolution processes have the potential to preserve old, relict surfaces through burial by cold-based, nonerosive glacial ice. To investigate landscape history and age in the high Arctic, we analyzed in situ cosmogenic Be(sup 10) and Al (sup 26) in 33 rocks from Upernavik, northwest Greenland. We sampled adjacent bedrock-boulder pairs along a 100 km transect at elevations up to 1000 m above sea level. Bedrock samples gave significantly older apparent exposure ages than corresponding boulder samples, and minimum limiting ages increased with elevation. Two-isotope calculations Al(sup26)/B(sup 10) on 20 of the 33 samples yielded minimum limiting exposure durations up to 112 k.y., minimum limiting burial durations up to 900 k.y., and minimum limiting total histories up to 990 k.y. The prevalence of BE(sup 10) and Al(sup 26) inherited from previous periods of exposure, especially in bedrock samples at high elevation, indicates that these areas record long and complex surface exposure histories, including significant periods of burial with little subglacial erosion. The long total histories suggest that these high elevation surfaces were largely preserved beneath cold-based, nonerosive ice or snowfields for at least the latter half of the Quaternary. Because of high concentrations of inherited nuclides, only the six youngest boulder samples appear to record the timing of ice retreat. These six samples suggest deglaciation of the Upernavik coast at 11.3 +/- 0.5 ka (average +/- 1 standard deviation). There is no difference in deglaciation age along the 100 km sample transect, indicating that the ice-marginal position retreated rapidly at rates of approx.120 m yr(sup−1).

glacial↗

Changes in Black Carbon Deposition to Antarctica from Two Ice Core Records, A.D. 1850-2000

Continuous flow analysis was based on a steady sample flow and in-line detection of BC and other chemical substances as described in McConnell et al. (2007). In the cold room, previously cut one meter ice core sticks of 3x3cm, are melted continuously on a heated melter head specifically designed to eliminate contamination from the atmosphere or by the external parts of the ice. The melted ice from the most inner part of the ice stick is continuously pumped by a peristaltic pump and carried to a clean lab by Teflon lines. The recorded signal is continuous, integrating a sample volume of about 0.05 mL, for which the temporal resolution depends on the speed of melting, ice density and snow accumulation rate at the ice core drilling site. For annual accumulation derived from the WAIS and Law Dome ice cores, we assumed ~3.1 cm water equivalent uncertainty in each year's accumulation from short scale spatial variability (glaciological noise) which was determined from several measurements of annual accumulation in multiple parallel ice cores notably from the WAIS Divide ice core site (Banta et al., 2008) and from South Pole site (McConnell et al., 1997; McConnell et al., 2000). Refractory black carbon (rBC) concentrations were determined using the same method as in (Bisiaux et al., 2011) and adapted to continuous flow measurements as described by (McConnell et al., 2007). The technique uses a single particle intracavity laser induced incandescence photometer (SP2, Droplet Measurement Technologies, Boulder, Colorado) coupled to an ultrasonic nebulizer/desolvation (CETAC UT5000) Flow Injection Analysis (FIA). All analyses, sample preparation etc, were performed in a class 100 cleanroom using anti contamination "clean techniques". The samples were not acidified.

Bisiaux, Marion M.↗

An Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation of the ICESat-2 Mission

The Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission has been selected by NASA as a Decadal Survey mission, to be launched in 2016. Mission objectives are to measure land ice elevation, sea ice freeboard/ thickness and changes in these variables and to collect measurements over vegetation that will facilitate determination of canopy height, with an accuracy that will allow prediction of future environmental changes and estimation of sea-level rise. The importance of the ICESat-2 project in estimation of biomass and carbon levels has increased substantially, following the recent cancellation of all other planned NASA missions with vegetation-surveying lidars. Two innovative components will characterize the ICESat-2 lidar: (1) Collection of elevation data by a multi-beam system and (2) application of micropulse lidar (photon counting) technology. A micropulse photon-counting altimeter yields clouds of discrete points, which result from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of returned points to reflectors of interest including canopy and ground in forested areas. The objective of this paper is to derive and validate an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2-type data. Data are based on airborne observations with a Sigma Space micropulse lidar and vary with respect to signal strength, noise levels, photon sampling options and other properties. A mathematical algorithm is developed, using spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors and geostatistical classification parameters and hyperparameters. Validation shows that the algorithm works very well and that ground and canopy elevation, and hence canopy height, can be expected to be observable with a high accuracy during the ICESat-2 mission. A result relevant for instrument design is that even the two weaker beam classes considered can be expected to yield useful results for vegetation measurements (93.01-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp9) and 72.85% - 98.68% for 0.48 msp (msp4)). Resampling options affect results more than noise levels. The algorithm derived here is generally applicable for analysis of micropulse lidar altimeter data collected over forested areas as well as other surfaces, including land ice, sea ice and land surfaces.

Herzfeld, Ute C.↗

Greenland Ice Sheet Mass Balance: Distribution of Increased Mass Loss with Climate Warming; 2003-07 Versus 1992-2002

We derive mass changes of the Greenland ice sheet (GIS) for 2003-07 from ICESat laser altimetry and compare them with results for 1992-2002 from ERS radar and airborne laser altimetry. The GIS continued to grow inland and thin at the margins during 2003 07, but surface melting and accelerated flow significantly increased the marginal thinning compared with the 1990s. The net balance changed from a small loss of 7 plus or minus 3 Gt a 1(sup -1) in the 1990s to 171 plus or minus 4 Gt a (sup -1) for 2003-07, contributing 0.5 mm a(sup -1) to recent global sea-level rise. We divide the derived mass changes into two components: (1) from changes in melting and ice dynamics and (2) from changes in precipitation and accumulation rate. We use our firn compaction model to calculate the elevation changes driven by changes in both temperature and accumulation rate and to calculate the appropriate density to convert the accumulation-driven changes to mass changes. Increased losses from melting and ice dynamics (17-206 Gt a(sup-1) are over seven times larger than increased gains from precipitation (10 35 Gt a(sup-1) during a warming period of approximately 2 K (10 a)(sup -1) over the GIS. Above 2000m elevation, the rate of gain decreased from 44 to 28 Gt a(sup-1), while below 2000m the rate of loss increased from 51 to 198 Gt a(sup-1). Enhanced thinning below the equilibrium line on outlet glaciers indicates that increased melting has a significant impact on outlet glaciers, as well as accelerating ice flow. Increased thinning at higher elevations appears to be induced by dynamic coupling to thinning at the margins on decadal timescales.

Zwally, H. Jay↗