Search NASA⌕ Search

SEARCH · Search NASA

Results for “MASS BALANCE”

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 145 records · Page 8

Spatial and Temporal Antarctic Ice Sheet Mass Trends, Glacio-Isostatic Adjustment, and Surface Processes from a Joint Inversion of Satellite Altimeter, Gravity, and GPS Data

We present spatiotemporal mass balance trends for the Antarctic Ice Sheet from a statistical inversion of satellite altimetry, gravimetry, and elastic-corrected GPS data for the period 2003-2013. Our method simultaneously determines annual trends in ice dynamics, surface mass balance anomalies, and a time-invariant solution for glacio-isostatic adjustment while remaining largely independent of forward models. We establish that over the period 2003-2013, Antarctica has been losing mass at a rateof -84 +/- 22 Gt per yr, with a sustained negative mean trend of dynamic imbalance of -111 +/- 13 Gt per yr. West Antarctica is the largest contributor with -112 +/- 10 Gt per yr, mainly triggered by high thinning rates of glaciers draining into the Amundsen Sea Embayment. The Antarctic Peninsula has experienced a dramatic increase in mass loss in the last decade, with a mean rate of -28 +/- 7 Gt per yr and significantly higher values for the most recent years following the destabilization of the Southern Antarctic Peninsula around 2010. The total mass loss is partly compensated by a significant mass gain of 56 +/- 18 Gt per yr in East Antarctica due to a positive trend of surface mass balance anomalies.

Martin-Espanol, Alba↗

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↗

Unprecedented Atmospheric Conditions (1948–2019) Drive the 2019 Exceptional Melting Season Over the Greenland Ice Sheet

Understanding the role of atmospheric circulation anomalies on the surface mass balance of the Greenland ice sheet (GrIS) is fundamental for improving estimates of its current and future contributions to sea level rise. Here, we show, using a combination of remote sensing observations, regional climate model outputs, reanalysis data, and artificial neural networks, that unprecedented atmospheric conditions (1948–2019) occurring in the summer of 2019 over Greenland promoted new record or close-to-record values of surface mass balance (SMB), runoff, and snowfall. Specifically, runoff in 2019 ranked second within the 1948–2019 period (after 2012) and first in terms of surface mass balance negative anomaly for the hydrological year 1 September 2018–31 August 2019. The summer of 2019 was characterized by an exceptional persistence of anticyclonic conditions that, in conjunction with low albedo associated with reduced snowfall in summer, enhanced the melt–albedo feedback by promoting the absorption of solar radiation and favored advection of warm, moist air along the western portion of the ice sheet towards the north, where the surface melt has been the highest since 1948. The analysis of the frequency of daily 500 hPa geopotential heights obtained from artificial neural networks shows that the total number of days with the five most frequent atmospheric patterns that characterized the summer of 2019 was 5 standard deviations above the 1981–2010 mean, confirming the exceptional nature of the 2019 season over Greenland.

Greenland ice sheet↗

Influence of nonconservative differencing on transonic streamline shapes

A computer program recently developed by South and Brandt (1976) which contained the Murman (1973) conservative finite-difference scheme is easily modified to use the Garabedian and Korn (1971) nonconservative finite difference scheme. This program solves the transonic small disturbance equation for only symmetric flow, but incorporates several iterative solution techniques. Results are presented for the case where the equally spaced computational grid extended to infinity in both the streamwise and normal directions. Streamline shapes are obtained along several grid lines by a streamwise integration of the normal component of the perturbation velocity. Comparison cases are run for a 10% thick parabolic arc airfoil at zero incidence for freestream Mach numbers of 0, 0.70, 0.84, and 0.95. It is shown that the use of a nonconservative finite-difference scheme in transonic flow calculations destroys the global mass balance when shocks are present. This lack of mass balance may prove to be more crucial in the case of an unconfined external flow.

Newman, P. A.↗

Determining Greenland Ice Sheet Accumulation Rates from Radar Remote Sensing

An important component of NASA's Program for Arctic Regional Climate Assessment (PARCA) is a mass balance investigation of the Greenland Ice Sheet. The mass balance is calculated by taking the difference between the snow accumulation and the ice discharge of the ice sheet. Uncertainties in this calculation include the snow accumulation rate, which has traditionally been determined by interpolating data from ice core samples taken throughout the ice sheet. The sparse data associated with ice cores, coupled with the high spatial and temporal resolution provided by remote sensing, have motivated scientists to investigate relationships between accumulation rate and microwave observations.

Jezek, Kenneth C.↗

Determining Greenland Ice Sheet Accumulation Rates from Radar Remote Sensing

An important component of NASA's Program for Arctic Regional Climate Assessment (PARCA) is a mass balance investigation of the Greenland Ice Sheet. The mass balance is calculated by taking the difference between the areally Integrated snow accumulation and the net ice discharge of the ice sheet. Uncertainties in this calculation Include the snow accumulation rate, which has traditionally been determined by interpolating data from ice core samples taken from isolated spots across the ice sheet. The sparse data associated with ice cores juxtaposed against the high spatial and temporal resolution provided by remote sensing , has motivated scientists to investigate relationships between accumulation rate and microwave observations as an option for obtaining spatially contiguous estimates. The objective of this PARCA continuation proposal was to complete an estimate of surface accumulation rate on the Greenland Ice Sheet derived from C-band radar backscatter data compiled in the ERS-1 SAR mosaic of data acquired during, September-November, 1992. An empirical equation, based on elevation and latitude, is used to determine the mean annual temperature. We examine the influence of accumulation rate, and mean annual temperature on C-band radar backscatter using a forward model, which incorporates snow metamorphosis and radar backscatter components. Our model is run over a range of accumulation and temperature conditions. Based on the model results, we generate a look-up table, which uniquely maps the measured radar backscatter, and mean annual temperature to accumulation rate. Our results compare favorably with in situ accumulation rate measurements falling within our study area.

Jezek, Kenneth C.↗

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↗

Ice Surface Elevation Changes in East Antarctica from Satellite Altimetry

Estimates of the overall mass balance and seasonal and inter-annual variations in the surface mass balance are obtainable from time-series of ice surface elevations measured by satellite altimetry. Beginning in 2001, NASA's ICESat laser altimeter and lidar mission will significantly improve the range accuracy, the orbit accuracy, and the spatial coverage for measurement of ice sheet elevations (to 86 S) , as compared to previous radar altimeters designed for ocean measurements The radar altimeters on Seasat and Geosat provided ice sheet measurements to 72 S, and on ERS-1 and ERS-2 to 81 S. Although radar altimetry has significant limitations in coverage (due to loss of tracking) and accuracy over sloping surfaces, information on ice-sheet surface-elevation changes has been derived for parts of Antarctica. Recently, the accuracy of the ice measurements by Seasat (3 months of 1978) and Geosat (1985 to 1989) have been improved by new calculations of the satellite orbit heights and other altimeter corrections. Residual orbit errors and inter-satellite biases are evaluated by crossover analysis and by global adjustments to an ocean surface derived from altimeter data. The standard deviation of the orbit error is less than 9 cm, and the long-term trend in the error appears to be less than 1 cm/yr. Orbit errors can be further reduced by adjustment to the ocean surface, but false signals of several cm/yr may be also introduced by the adjustments. These false signals are caused mainly by residual errors in the altimeter corrections over the ocean, and secondary by real changes in the ocean surface elevation. Maps of ice sheet elevation changes north of 72 S are derived from Seasat-Geosat crossovers and from 4.5 years of Geosat crossovers. A notable ice thinning rate of about 50 cm/yr is found at elevations below 2200 meters between 70 and 72 S to the East of the Amery ice shelf, in both the Seasat-Geosat and Geosat-Geosat time intervals Above 2200 meters, to the ridge around 95 E, the thinning is about 10 to 20 cm/yr. For the Seasat/Geosat interval, a thickening is noted on the area East of the ridge between about 95 E and 115 E. Preliminary results from analysis of ERS-1 and ERS-2 radar altimeter data are also described.

Zwally, H. Jay↗

Southeast Greenland Winter Precipitation Strongly Linked to the Icelandic Low Position

Greenland's largest precipitation flux occurs in its southeast (SE) region during the winter, controlled primarily by easterly winds and frequent cyclogenesis in the North Atlantic. Several studies have attempted to link SE Greenland precipitation to the North Atlantic Oscillation (NAO) but results are inconsistent. This work uses reanalysis, automatic weather station data, and regional climate model output to show that the east-west position of the Icelandic low is a better predictor of SE Greenland precipitation (average correlation of r equals -0.48 in DJF (December January February)) than climate indices such as the NAO (r equals -0.06 in DJF). In years when the Icelandic low is positioned extremely west, moisture transport increases up to approximately 40 percent (or up to 40 kilograms per meter per second) off the SE Greenland coast compared to when the low is in an extreme east position. Furthermore, in years when the Icelandic low is positioned extremely west, storm track density and intensity increase just off the SE coast of Greenland. Thus, the Icelandic low's longitudinal position dominates SE Greenland ice sheet's wintertime precipitation, a positive term in the ice sheet mass balance. Given SE Greenland's importance in the overall ice sheet mass balance, the position of the Icelandic low is therefore important for making projections of future sea level.

Arctic↗

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗

VRA Modeling, phase 1

The destruction of organic contaminants in waste water for closed systems, such as that of Space Station, is crucial due to the need for recycling the waste water. A co-current upflow bubble column using oxygen as the gas phase oxidant and packed with catalyst particles consisting of a noble metal on an alumina substrate is being developed for this process. The objective of this study is to develop a plug-flow model that will predict the performance of this three phase reactor system in destroying a multicomponent mixture of organic contaminants in water. Mass balances on a series of contaminants and oxygen in both the liquid and gas phases are used to develop this model. These mass balances incorporate the gas-to-liquid and liquid-to-particle mass transfer coefficients, the catalyst effectiveness factor, and intrinsic reaction rate. To validate this model, a bench scale reactor has been tested at Michigan Technological University at elevated pressures (50-83 psig,) and a temperature range of 200 to 290 F. Feeds consisting of five dilute solutions of ethanol (approx. 10 ppm), chlorobenzene (approx. 20 ppb), formaldehyde (approx. 100 ppb), dimethyl sulfoxide (DMSO approx. 300 ppb), and urea (approx. 20 ppm) in water were tested individually with an oxygen mass flow rate of 0.009 lb/h. The results from these individual tests were used to develop the kinetic parameter inputs necessary for the computer model. The computer simulated results are compared to the experimental data obtained for all 5 components run in a mixture on the differential test column for a range of reactor contact times.

Kindt, Louis M.↗

Estimation of Thermodynamic and Dynamic Contributions to Sea Ice Growth in the Central Arctic Using ICESat-2 and MOSAiC SIMBA Buoy Data

The fine spatial resolution of the ICESat-2 (IS2) satellite altimeter allows monitoring the evolution of sea ice thickness with detailed dynamic information (e.g. ridges and leads). In this study, we first assess the ability of IS2 to estimate thermodynamic ice growth and dynamic thickening during the ice-growing season in the central Arctic Ocean. As an indicator of the thermodynamic ice growth, we use 10 thermistor string-based sea ice mass balance array (SIMBA) buoys deployed at a scale of ~ 50 km from the Icebreaker Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. We collect IS2 data within 20 km buffer distance from the individual buoys, and calculate the mode, median, and mean of the IS2-derived ice thickness. The IS2 modal thickness shows the least bias (−0.169 m) with the buoy ice thickness, representing level ice thickness. In addition, the increasing rate of the IS2 modal thickness is close to the thermodynamic ice growth with a small bias of −0.054 cm/day. However, the increasing rates of the IS2 median and mean thickness are greater than the thermodynamic ice growth by about 0.114 cm/day and 0.198 cm/day, respectively, because they also include ice growth caused by thickness redistribution during dynamic deformation. The dynamic contributions may account for 26.1 ± 10.3% and 34.4 ± 10.1% of the total increase of the IS2 median and mean thickness, respectively. Within a ~ 50 km radius area from the MOSAiC Central Observatory, IS2 measurements exhibit that the ridge fraction increased from <2% in November to ~ 4% in March (~0.029%/day of average increasing rate) and ridge height increased about 0.047 cm/day during the same period. However, lead formation does not show significant contributions to the dynamic ice thickening because leads are temporary features lasting only 2–3 days. Although there are considerable uncertainties in IS2 ice thickness estimation and IS2-buoy thickness comparison, this study emphasizes the importance of combining measurements by IS2 and SIMBA buoys to explain the regional sea ice mass balance with separating the thermodynamic and dynamic contributions.

satellite altimeter↗

Resource Tracking Model Updates and Trade Studies

The Resource tracking model has been updated to capture system manager and project manager inputs. Both the Trick/GUNNS RTM simulator and the RTM mass balance spreadsheet have been revised to address inputs from system managers and to refine the way mass balance is illustrated. The revisions to the RTM included addition of a Plasma Pyrolysis Assembly (PPA) to recover hydrogen from Sabatier reactor methane which was vented in the prior version of the RTM. The effect of the PPA on the overall balance of resources in an exploration vehicle is illustrated in the increased recycle of vehicle oxygen. Additionally simulation of EVAs conducted from the exploration module was added. Since the focus of the exploration module is to provide a habitat during deep space operations the EVA simulation approach to EVA is based on ISS EVA protocol and processes. Case studies have been run to show the relative effect of performance changes on vehicle resources.

Chambliss, Joe↗

Resource Tracking Model Updates and Trade Studies

The Resource Tracking Model has been updated to capture system manager and project manager inputs. Both the Trick/General Use Nodal Network Solver Resource Tracking Model (RTM) simulator and the RTM mass balance spreadsheet have been revised to address inputs from system managers and to refine the way mass balance is illustrated. The revisions to the RTM included the addition of a Plasma Pyrolysis Assembly (PPA) to recover hydrogen from Sabatier Reactor methane, which was vented in the prior version of the RTM. The effect of the PPA on the overall balance of resources in an exploration vehicle is illustrated in the increased recycle of vehicle oxygen. Case studies have been run to show the relative effect of performance changes on vehicle resources.

Chambliss, Joe↗

Glacier and Ice Shelves Studies Using Satellite SAR Interferometry

Satellite radar interferometry is a powerful technique to measure the surface velocity and topography of glacier ice. On ice shelves, a quadruple difference technique separates tidal motion from the steady creep flow deformation of ice. The results provide a wealth of information about glacier grounding lines , mass fluxes, stability, elastic properties of ice, and tidal regime. The grounding line, which is where the glacier detaches from its bed and becomes afloat, is detected with a precision of a few tens of meters. Combining this information with satellite radar altimetry makes it possible to measure glacier discharge into the ocean and state of mass balance with greater precision than ever before, and in turn provide a significant revision of past estimates of mass balance of the Greenland and Antarctic Ice Sheets. Analysis of creep rates on floating ice permits an estimation of basal melting at the ice shelf underside. The results reveal that the action of ocean water in sub-ice-shelf cavities has been largely underestimated by oceanographic models and is the dominant mode of mass release to the ocean from an ice shelf. Precise mapping of grounding line positions also permits the detection of grounding line migration, which is a fine indicator of glacier change, independent of our knowledge of snow accumulation and ice melting. This technique has been successfully used to detect the rapid retreat of Pine Island Glacier, the largest ice stream in West Antarctica. Finally, tidal motion of ice shelves measured interferometrically provides a modern, synoptic view of the physical processes which govern the formation of tabular icebergs in the Antarctic.

Rignot, Eric↗

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗

Flutter calculations in three degrees of freedom

The present paper is a continuation of the general study of flutter published in NACA reports nos. 496 and 685. The paper is mainly devoted to flutter in three degrees of freedom (bending, torsion, and aileron) for which a number of selected cases have been calculated and presented in graphical form. The results are analyzed and discussed with regard to the effects of structural damping, of fractional-span ailerons, and of mass-balancing. The analysis shows that more emphasis should be put on the effect of structural damping and less on mass-balancing. The conclusion is drawn that a definite minimum amount of structural damping, which is usually found to be present, is essential in the calculations for an adequate description of the flutter case. Theoretical flutter predictions are thus brought into closer agreement with the facts of experience. A brief discussion is included of a particular biplane that had experienced flutter at about 200 miles per hour. Some simplifications have been achieved in the method of calculation. (author)

Theodorsen, Theodore↗

A Preliminary Model for Spacecraft Propulsion Performance Analysis Based on Nuclear Gain and Subsystem Mass-Power Balances

Rapid transportation of human crews to destinations throughout the solar system will require propulsion systems having not only very high exhaust velocities (i.e., I(sub sp) >= 10(exp 4) to 10(exp 5) sec) but also extremely low mass-power ratios (i.e., alpha <= 10(exp -2) kg/kW). These criteria are difficult to meet with electric propulsion and other power-limited systems, but may be achievable with propulsion concepts that use onboard power to produce a net gain in energy via fusion or some other nuclear process. This paper compares the fundamental performance of these gain-limited systems with that of power-limited systems, and determines from a generic power balance the gains required for ambitious planetary missions ranging up to 100 AU. Results show that energy gain reduces the required effective mass-power ratio of the system, thus enabling shorter trip times than those of power-limited concepts.

Chakrabarti, S.↗