Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DEM”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

New electric force and charge exchange modules in discrete element model enables particle dynamics simulation in electric field

Discrete element modeling (DEM) is an important technique for particle dynamics simulation. The field of metal additive manufacturing often utilizes DEM to simulate the rheological behaviors of powder. Standard contact and short-range interactions are sufficient in most cases but insufficient to describe the particle dynamics with the influence of an electric field. Modeling such a system requires additional physics to describe the particle–field interactions. The relevant physics has been experimentally understood but is not yet available in DEM. Here, we develop a charge exchange and an electric force module. The electric force module governs particle response to the electric field, while the charge exchange module enables particles to acquire proper charge during contact with charged geometries. We validate the modules against analytical calculations and high-speed videos of electrostatic powder deposition experiments. Notably, the model struggles to capture the initial particle levitation. We later deploy a modified electric field, as supported by static electric field simulation, to better approximate the electric field penetration into the powder layer. This modification improves the model’s capability of simulating realistic particle levitation. The results highlight the challenges of modeling particle behaviors in the electric field while demonstrating the feasibility of obtaining quantitative results, which are difficult to measure experimentally.

charge exchange↗
DOE OSTI2025

Predicting the evolution of biomass bulk density through feedstock preprocessing: Discrete element modeling, regression analysis, and pilot-scale validation

Bulk density is an important material property of biomass feedstocks, influencing handling, storage, transport costs, and conversion efficiency. In this study, predictive regression models for loose and tapped bulk densities of Alamo and Cave-in-Rock switchgrass are developed using a comprehensive dataset generated via calibrated bonded-sphere discrete element method (DEM) simulations. Here, a key contribution of this study is the use of a DEM-based approach, which correlates density with moisture content and particle size distribution parameters and enables analysis across a continuous particle size range, overcoming limitations of purely experimental data. For comparison, regression models are also developed using only experimental data from pilot-scale runs at the Biomass Feedstock National User Facility at Idaho National Laboratory. Validation against pilot-scale data showed reasonable prediction accuracy for both model types, particularly for smaller particle sizes (post-secondary grinding). While the experimental model showed slightly better performance matching the validation data in some cases, the DEM-based model benefits from a much larger dataset, reduced predictor multicollinearity, and continuous parameter coverage, highlighting the utility of validated simulation models for developing robust predictive tools for biomass preprocessing applications.

09 - BIOMASS FUELS↗
DOE OSTI2026

High-resolution mountain topography can inform global snow vulnerability estimates

Snow is changing globally. Computationally intensive snow reanalysis products and downscaled climate model projections allow for the estimation of historical and projected changes in snow over ∼4–10 km resolutions, but these resolutions are coarse relative to the scales needed for water supply and flood planning. Fine-scale digital elevation models (DEMs) are widely available but are underutilized to make first-order assessments of snow vulnerability. Here, we leverage DEMs at a 7.5 arc s (∼250 m) resolution, combining these with historical freezing level height estimates from ERA-5 to derive estimates of changes in the snow-receiving area (SRA) and its variability across global mountain ranges. Results show estimated SRA declines in 29% (1.9 million km2) of the global mountain area from 1982–2020; 66% of the mountainous areas had no change over the historical period. At +1.5 °C of warming relative to the pre-industrial control, global mountain SRA would decline by 9.5% (1.0 million km2) relative to recent conditions. This loss would be approximately doubled with +2 °C of warming. In a +4 °C warming scenario, an additional 34% (3.6 million km2) of SRA would be lost beyond the +2 °C case. Across individual mountain ranges, SRA losses can occur nonlinearly with warming, with some locations that have historically had relatively minor SRA losses at risk of substantially larger losses in warmer climates. Analysis using coarser-resolution DEMs can underestimate or overestimate SRA and its rate of loss, with the largest impacts in relatively warm, low-elevation mountain ranges. Results of this work provide estimates of projected loss in SRA at policy-relevant warming levels; inform the resolutions needed for process-based snow modeling; identify snow vulnerability hotspots; and provide a new integrated approach to snow vulnerability assessment that is achievable at global scales and highlights potential nonlinearities from recent trends to a variety of future warming scenarios.

climate, mountains↗
DOE OSTI2025

XRISM view of a stellar flare: High-resolution Fe K spectra of HR 1099, an RS CVn-type star

A high-resolution X-ray spectroscopic observation was made of the RS CVn-type binary star HR 1099 using the Resolve instrument onboard XRISM for its calibration purposes. During the $\sim$400 ks telescope time covering 1.5 binary orbit, a flare lasting for $\sim$100 ks was observed with a released X-ray radiation energy of ${\sim }10^{34}$ erg, making it the first stellar flare ever observed with an X-ray microcalorimeter spectrometer. The flare peak count rate is 6.4 times higher than that in quiescence and is distinguished clearly in time thanks to the long telescope time. Many emission lines were detected in the 1.7–10 keV range both in the flare and quiescent phases. Using the high spectral resolution of Resolve in the Fe K band (6.5–7.0 keV), we resolved the inner-shell lines of Fe XIX-XXVI as well as the outer-shell lines of Fe XXV-XXVI . These lines have peaks in the contribution functions at different temperatures over a wide range, allowing us to construct the differential emission measure (DEM) distribution over the electron temperature of 1–10 keV (roughly 10–100 MK) based only on Fe lines, thus without an assumption of the elemental abundance. The reconstructed DEM has a bimodal distribution, and only the hotter component increased during the flare. The elemental abundance was derived based on the DEM distribution thus constructed. A significant abundance increase was observed during the flare for Ca and Fe, which are some of the elements with the lowest first ionization potential among those analyzed, but not for Si, S, and Ar. This behavior is seen in some giant solar flares and the present result is a clear example in stellar flares.

Astronomy and AstroPhysics↗
DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗
DOE OSTI2025

TRAILS Output Files

Overview This data repository contains ZIP files that store compressed versions of the output of running the WaterPaths utility planning and management tool in the DU Re-Evaluation mode (to download the tool, please see this GitHub repository). The tool was used to simulate the six-utility North Carolina Research Triangle problem. Details on the contents of each ZIP file can be seen below. Data details Temporal range: Weekly data for 2,344 weeks from 2015 to 2060 (45 years). Spatial range: Six water utilities in the North Carolina Research Triangle region (0: Chapel Hil/OWASA, 1: Durham, 2: Cary, 3: Raleigh, 4: Pittsboro, and 5: Chatham) File types: CSV and OUT Different solutions available The solution numbers correspond to the different pathway strategies (henceforth referred to as "solutions") discussed in paper's main and supporting text (abstract and link to the paper here). They are as follows: Sol92: The Durham-focused pathway strategy Sol132: The Raleigh-focused pathway strategy Sol140: The regionally-robust pathway strategy Objectives files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Objectives_RDMXX_solsXX_to_XX.csv files. Each CSV file will consist of a row representing all the objective values for that specific solution, while every six columns represents the reliability, restriction frequency, infrastructure net present value ($ mil), peak financial cost, worst-case cost, and unit cost ($ per MG; in that order) for each of the six utilities. There will be 1,000 such files, denoting the performance of the six utilities across the 1,000 deeply uncertain states of the world (DU SOWs). Pathway files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Pathways_sXX_RDMXX.out file. Each OUT corresponds to the set of infrastructure being triggered in a specific DU SOW, and each file will have the name file will consist of four tab-delimited columns that are described as follows: Realization: The realization in which an infrastructure options being triggered utility: The utility currently triggering infrastructure week: The week in which a specific infrastructure option is being triggered infra.: The infrastructure option being triggered If the OUT file contains only the header line, no infrastructure was triggered for that specific DU SOW. Policies files These files can be obtained by unzipping Policies.zip. Each of the 1,000 CSV files within the unzipped folder will contain weekly water use restriction policies for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: 0rest_m: restriction multiplier for utility 0 (values between 0 and 1) 1rest_m: restriction multiplier for utility 1 (values between 0 and 1) 2rest_m: restriction multiplier for utility 2 (values between 0 and 1) 3rest_m: restriction multiplier for utility 3 (values between 0 and 1) 4rest_m: restriction multiplier for utility 4 (values between 0 and 1) 5rest_m: restriction multiplier for utility 5 (values between 0 and 1) 0transf: transfer volume for utility 0 (in MGD) 1transf: transfer volume for utility 1 (in MGD) 2transf: transfer volume for utility 2 (in MGD) 3transf: transfer volume for utility 3 (in MGD) 4transf: transfer volume for utility 4 (in MGD) 5transf: transfer volume for utility 5 (in MGD) Water Sources files These files can be obtained by unzipping WaterSources_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each water source for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xvolume: available water volume from source X (in MGD) Xs_area: surface area of source X (in ACF) Xdemand: demand drawn from a water source from source X (in MGD) Xup_spill: upstream spillage from source X (in MGD) Xww_inflow: wastewater inflow from source X (in MGD) Xcatch_inflow: upstream catchment inflow to source X (in MGD) Xevap: evaporation multiplier for source X (values between 0 and 1) Xds_spill: downstream spillage from source X (in MGD) X_Y_alloc_cap: the allocated capacity from source X to utility Y (values between 0 and 1) X_Y_alloc_dem: the allocated demand from source X to utility Y (values between 0 and 1) Xtrmt_alloc_Y: the allocated treatment capacity from source X to utility Y (values between 0 and 1) Utilities files These files can be obtained by unzipping Utilities_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each utility for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xst_vol: total available storage volume of utility X (in MG) Xcapacity: total storage capacity of utility X (in MG) Xnet_inf: : net inflow for all storage infrastructure for utility X (in MGD) Xst_rof: short term ROF for utility X (values between 0 and 1) Xst_stor_rof: short-term storage ROF for utility X (values between 0 and 1) Xst_trmt_rof: short-term treatment ROF for utility X (values between 0 and 1) Xlt_rof: long-term ROF for utility X (values between 0 and 1) Xlt_stor_rof: long-term storage ROF for utility X (values between 0 and 1) Xlt_trmt_rof: long-term treatment ROF for utility X (values between 0 and 1) Xrest_demand: restricted demand for utility X (in MGD) Xunrest_demand: unrestricted demand for utility X (in MGD) Xunfulf_demand: unfulfilled demand for utility X (in MGD) Xwastewater: wastewater return for utility X (in MGD) Xtreat_capacity: total treatment capacity for utility X (in MG) Xcont_fund: reserve (contingency) fund balance for utility X Xins_pout: insurance payout for utility X (% annual volumetric revenue) Xins_price: insurance price for utility X (% annual volumetric revenue) Xinfra_npv: infrastructure net present value for utility ($mil) Xst_vol: total available storage volume of utility X (in MG) Xdebt_serv: debt service for utility X (usually once per year if the infrastructure is triggered; % annual volumetric revenue) Xstor_vol: total stored volume (in MGD) Xobs_ann_dem: observed annual demand for utility X (in MGD) Xproj_dem: projected annual demand for utility X (in MGD) Xpv_debt_serv: present value of debt service payments for utility X (% annual volumetric revenue) Xgross_rev: gross revenue for utility X ($mil) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

Artificial Intelligence↗
NASA NTRS1987

Rectification of terrain induced distortions in radar imagery

This paper describes a technique to generate geocoded synthetic aperture radar (SAR) imagery corrected for terrain induced geometric distortions. This algorithm transforms the raw slant range image, generated by the signal processor, into a map registered product, resampled to either Universal Transverse Mercator (UTM) or Polar Stereographic projections, and corrected for foreshortening. The technique utilizes the space platform trajectory information in conjunction with a digital elevation map (DEM) of the target area to generate an ortho-radar map with near-autonomous operation. The current procedure requires only two to three tie-points to compensate for the platform position uncertainty that results in translational error between the image and the DEM. This approach is unique in that it does not require generation of a simulated radar image from the DEM or a grid of tie-points to characterize the image-to-map distortions. Rather, it models the inherent distortions based on knowledge of the radar data collection characteristics, the signal Doppler parameters, and the local terrain height to automatically predict the registration transformation. This algorithm has been implemented on a minicomputer system equipped with an array processor and a large random-access memory to optimize the throughput.

Kwok, Ronald↗
NASA NTRS1987

Stereocorrelation of Landsat TM images

A digital elevation model (DEM) developed from Landsat TM images of a rugged terrain area in north Georgia by automated stereocorrelation techniques yielded an rms error (z), RMSE(z), value of + or - 42 m. Based on the B/H ratio of 0.18 for the Landsat data, this Z-error corresponds to a planimetric correlation accuracy of about + or - 0.3 pixels, confirming that precise correlation can be achieved with operational satellite data. Contours at a 100-m interval interpolated from the DEM show a deviation of + or - 33 m from reference contours obtained from existing 1:24,000-scale maps. The 28.5-m pixel resolution and the weak B/H ratio impose limitations on the accuracy that can be achieved with Landsat TM data. However, it is anticipated that RMSE(z) values of + or - 10 m or less can be achieved with SPOT-1 panchromatic stereo images of 10-m resolution recorded at B/H ratios of 0.5 to 1.0. DEMs generated by stereocorrelation techniques can be used to create orthoimages, perspective views, and topographic map products.

Ehlers, Manfred↗
NASA NTRS1994

Abundance determination in solar active region

An improved method has been developed to evaluate element abundances from emission line intensities of thin plasmas, depending on the differential emission measure (DEM) of the source. Observations made by the Solar EUV Rocket Telescope and Spectrograph (SERTS) Rocket EUV Spectrograph are used to perform a detailed analysis of the DEM distribution for temperatures larger than 10(exp 5) K in a solar active region. Comparison of the DEM distributions obtained by means of lines from different elements allows the verification of relative abundances for the most common elements of the solar corona, and gives an abundance estimates for the minor components, such as Na, Al, Ar, Cr, Mn and Zn.

Monsignori Fossi, B. C.↗
NASA NTRS1994

The length-scaling properties of topography

The scaling properties of synthetic topographic surfaces and digital elevation models (DEMs) of topography are examined by analyzing their 'structure functions,' i.e., the qth order powers of the absolute elevation differences: delta h(sub q) (l) = E((absolute value of h(x + l) - h(x))(exp q)). We find that the relation delta h(sub 1 l) approximately equal cl(exp H) describes well the scaling behavior of natural topographic surfaces, as represented by DEMs gridded at 3 arc sec. Average values of the scaling exponent H between approximately 0.5 and 0.7 characterize DEMs from Ethiopia, Saudi Arabia, and Somalia over 3 orders of magnitude range in length scale l (approximately 0.1-150 km). Differences in appparent topographic roughness among the three areas most likely reflect differences in the amplitude factor c. Separate determination of scaling properties in the x and y coordinate directions allows us to assess whether scaling exponents are azimuthally dependent (anisotropic) or whether they are isotropic while the surface itself is anisotropic over a restricted range of length scale. We explore ways to determine whether topographic surfaces are characterized by simple or multiscaling properties.

Weissel, Jeffrey K.↗
NASA NTRS1996

Plasma Properties and Magnetic Field Structure of the Solar Corona, Based on Coordinated Max '91 Observations from SERTS, the VLA, and Magnetographs

The plasma properties and magnetic field structure of the solar corona were determined using coordinated observations obtained with NASA/GSFC's Solar EUV Rocket Telescope and Spectrograph (SERTS), the Very Large Array (VLA), and Kitt Peak photospheric longitudinal magnetograms. A problem was identified with the SERTS calibration as determined from laboratory measurements. A revised calibration curve was derived by requiring that the numerous available measured line intensity ratios agreed with their respective theoretical values. Densities were derived from line intensity ratios, and active region densities were found to typically exceed quiet Sun densities by factors of only about 2. The active region density was found to remain constant across the SERTS slit, despite the fact that the emission line intensities vary significantly. This indicates that the product of the path length and the volume filling factor must vary significantly from the active region outskirts to the central core. Filling factors were derived and found to range from much less than one to nearly unity. Wavelength shifts were examined along the SERTS slit in the spatially resolved spectra, but no evidence was found for significant Doppler shifts in active region 7563 or in the quiet Sun. The numerical procedure developed by Monsignori-Fossi and Landini was used to derive the active region and quiet sun differential emission measure (DEM) from the spatially averaged spectra. A DEM was estimated for each spatial pixel in the two dimensional active region images by scaling the averaged active region DEM based upon corresponding pixel intensities of SERTS Mg IX, Fe XV, and Fe XVI images. These results, along with density measurements, were used in an IDL computer code which calculated the temperature dependence of the coronal magnetic field in each spatial pixel by minimizing the difference between the observed and calculated 20 and 6 cm microwave brightness temperatures.

Brosius, Jeffrey W.↗
NASA NTRS1999

The Utility and Validity of Kinematic GPS Positioning for the Geosar Airborne Terrain Mapping Radar System

GeoSAR is an airborne, interferometric Synthetic Aperture Radar (IFSAR) system for terrain mapping, currently under development by a consortium including NASA's Jet Propulsion Laboratory (JPL), Calgis, Inc., a California mapping sciences company, and the California Department of Conservation (CaIDOC), with funding provided by the U.S. Army Corps of Engineers Topographic Engineering Center (TEC) and the U.S. Defense Advanced Research Projects Agency (DARPA). IFSAR data processing requires high-accuracy platform position and attitude knowledge. On 9 GeoSAR, these are provided by one or two Honeywell Embedded GPS Inertial Navigation Units (EGI) and an Ashtech Z12 GPS receiver. The EGIs provide real-time high-accuracy attitude and moderate-accuracy position data, while the Ashtech data, post-processed differentially with data from a nearby ground station using Ashtech PNAV software, provide high-accuracy differential GPS positions. These data are optimally combined using a Kalman filter within the GeoSAR motion measurement software, and the resultant position and orientation information are used to process the dual frequency (X-band and P-band) radar data to generate high-accuracy, high -resolution terrain imagery and digital elevation models (DEMs). GeoSAR requirements specify sub-meter level planimetric and vertical accuracies for the resultant DEMS. To achieve this, platform positioning errors well below one meter are needed. The goal of GeoSAR is to obtain 25 cm or better 3-D positions from the GPS systems on board the aircraft. By imaging a set of known point target corner-cube reflectors, the GeoSAR system can be calibrated. This calibration process yields the true position of the aircraft with an uncertainty of 20- 50 cm. This process thus allows an independent assessment of the accuracy of our GPS-based positioning systems. We will present an overview of the GeoSAR motion measurement system, focusing on the use of GPS and the blending of position data from the various systems. We will present the results of our calibration studies that relate to the accuracy the GPS positioning. We will discuss the effects these positioning, errors have on the resultant DEM products and imagery.

Freedman, Adam↗
NASA NTRS2002

Fractal Analysis of Drainage Basins on Mars

We used statistical properties of drainage networks on Mars as a measure of martian landscape morphology and an indicator of landscape evolution processes. We utilize the Mars Orbiter Laser Altimeter (MOLA) data to construct digital elevation maps (DEMs) of several, mostly ancient, martian terrains. Drainage basins and channel networks are computationally extracted from DEMs and their structures are analyzed and compared to drainage networks extracted from terrestrial and lunar DEMs. We show that martian networks are self-affine statistical fractals with planar properties similar to terrestrial networks, but vertical properties similar to lunar networks. The uniformity of martian drainage density is between those for terrestrial and lunar landscapes. Our results are consistent with the roughening of ancient martian terrains by combination of rainfall-fed erosion and impacts, although roughening by other fluvial processes cannot be excluded. The notion of sustained rainfall in recent Mars history is inconsistent with our findings.

Stepinski, T. F.↗
NASA NTRS1996

Low-Complexity, Digital Encoder/Modulator Developed for High-Data-Rate Satellite B-ISDN Applications

The Space Electronics Division at the NASA Lewis Research Center is developing advanced electronic technologies for the space communications and remote sensing systems of tomorrow. As part of the continuing effort to advance the state-of-the-art in satellite communications and remote sensing systems, Lewis developed a low-cost, modular, programmable, and reconfigurable all-digital encoder-modulator (DEM) for medium- to high-data-rate radiofrequency communication links. The DEM is particularly well suited to high-data-rate downlinks to ground terminals or direct data downlinks from near-Earth science platforms. It can support data rates up to 250 megabits per second (Mbps) and several modulation schemes, including the traditional binary phase-shift keying (BPSK) and quadrature phase-shift keying (QPSK) modes, as well as higher order schemes such as 8 phase-shift keying (8PSK) and 16 quadrature amplitude modulation (16QAM). The DEM architecture also can precompensate for channel disturbances and alleviate amplitude degradations caused by nonlinear transponder characteristics.

Source record↗
NASA NTRS2009

Topography of the Martian Impact Crater Tooting

Tooting crater is approx.29 km in diameter, is located at 23.4degN, 207.5degE, and is classified as a multi-layered ejecta crater [1]. Our mapping last year identified several challenges that can now be addressed with HiRISE and CTX images, but specifically the third dimension of units. To address the distribution of ponded sediments, lobate flows, and volatile-bearing units within the crater cavity, we have focused this year on creating digital elevation models (DEMs) for the crater and ejecta blanket from stereo CTX and HiRISE images. These DEMs have a spatial resolution of approx.50 m for CTX data, and 2 m for HiRISE data. Each DEM is referenced to all of the available individual MOLA data points within an image, which number approx.5,000 and 800 respectively for the two data types

Mouginis-Mark, P. J.↗
NASA NTRS2010

ICESat Lidar and Global Digital Elevation Models: Application to DESDynI

Geodetic control is extremely important in the production and quality control of topographic data sets, enabling elevation results to be referenced to an absolute vertical datum. Global topographic data with improved geodetic accuracy achieved using global Ground Control Point (GCP) databases enable more accurate characterization of land topography and its change related to solid Earth processes, natural hazards and climate change. The multiple-beam lidar instrument that will be part of the NASA Deformation, Ecosystem Structure and Dynamics of Ice (DESDynI) mission will provide a comprehensive, global data set that can be used for geodetic control purposes. Here we illustrate that potential using data acquired by NASA's Ice, Cloud and land Elevation Satellite (ICEsat) that has acquired single-beam, globally distributed laser altimeter profiles (+/-86deg) since February of 2003 [1, 2]. The profiles provide a consistently referenced elevation data set with unprecedented accuracy and quantified measurement errors that can be used to generate GCPs with sub-decimeter vertical accuracy and better than 10 m horizontal accuracy. Like the planned capability for DESDynI, ICESat records a waveform that is the elevation distribution of energy reflected within the laser footprint from vegetation, where present, and the ground where illuminated through gaps in any vegetation cover [3]. The waveform enables assessment of Digital Elevation Models (DEMs) with respect to the highest, centroid, and lowest elevations observed by ICESat and in some cases with respect to the ground identified beneath vegetation cover. Using the ICESat altimetry data we are developing a comprehensive database of consistent, global, geodetic ground control that will enhance the quality of a variety of regional to global DEMs. Here we illustrate the accuracy assessment of the Shuttle Radar Topography Mission (SRTM) DEM produced for Australia, documenting spatially varying elevation biases of several meters in magnitude.

Carabajal, Claudia C.↗
NASA NTRS2009

Elevation Change of Drangajokull, Iceland, from Cloud-Cleared ICESat Repeat Profiles and GPS Ground-Survey Data

Located on the Vestfirdir Northwest Fjords), DrangaJokull is the northernmost ice map in Iceland. Currently, the ice cap exceeds 900 m in elevation and covered an area of approx.l46 sq km in August 2004. It was about 204 sq km in area during 1913-1914 and so has lost mass during the 20th century. Drangajokull's size and accessibility for GPS surveys as well as the availability of repeat satellite altimetry profiles since late 2003 make it a good subject for change-detection analysis. The ice cap was surveyed by four GPS-equipped snowmobiles on 19-20 April 2005 and has been profiled in two places by Ice, Cloud. and land Elevation Satellite (ICESat) 'repeat tracks,' fifteen times from late to early 2009. In addition, traditional mass-balance measurements have been taken seasonally at a number of locations across the ice cap and they show positive net mass balances in 2004/2005 through 2006/2007. Mean elevation differences between the temporally-closest ICESat profiles and the GPS-derived digital-elevation model (DEM)(ICESat - DEM) are about 1.1 m but have standard deviations of 3 to 4 m. Differencing all ICESat repeats from the DEM shows that the overall elevation difference trend since 2003 is negative with losses of as much as 1.5 m/a from same season to same season (and similar elevation) data subsets. However, the mass balance assessments by traditional stake re-measurement methods suggest that the elevation changes where ICESat tracks 0046 and 0307 cross Drangajokull are not representative of the whole ice cap. Specifically, the area has experienced positive mass balance years during the time frame when ICESat data indicates substantial losses. This analysis suggests that ICESat-derived elevations may be used for multi-year change detection relative to other data but suggests that large uncertainties remain. These uncertainties may be due to geolocation uncertainty on steep slopes and continuing cloud cover that limits temporal and spatial coverage across the area.

Shuman, Christopher A.↗