Search NASA⌕ Search

SEARCH · Search NASA

Results for “Icing”

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 595 records · Page 33

Comparison of Ice Shapes on Full-Chord and Truncated Swept Wing Models

A research program was conducted to evaluate the effectiveness of icing tunnel hybrid model design. A hybrid design is where the full-scale leading edge of a wing section is maintained only to a certain percentage of the local chord, while the aft section of the model is redesigned into a shortened or truncated planform. Ice shape geometries on a full-chord length and hybrid versions of the swept CRM65 wing model were obtained in NASA Icing Research Tunnel and compared. For most test conditions, the overall size and shape of the ice shapes compared well. On the main part of the ice shape, the hybrid model had larger chordwise extent, resulting in wider ice shape. For the scallop ice shapes, the ice shape heights at attachment line location were identical. For conditions that produced rime shape, hybrid model ice shapes were larger overall, even at the attachment line location. The only location where there was more ice on the full chord model was on the lower surface, downstream of where the hybrid model no longer had full-scale leading edge. For the test conditions where ice mass measurements were acquired for both models, the full-chord models had 11-22% less ice mass than on the hybrid wing. These differences were larger and opposite of those predicted from LEWICE 3D ice prediction calculations. Further investigation is needed in order determine the cause of the differences between the full-chord and hybrid models.

Icing↗

Comparison of Ice Shapes on Full-Chord and Truncated Swept Wing Models

A research program was conducted to evaluate the effectiveness of icing tunnel hybrid model design. A hybrid design is where the full-scale leading edge of a wing section is maintained only to a certain percentage of the local chord, while the aft section of the model is redesigned into a shortened or truncated planform. Ice shape geometries on a full-chord length and hybrid versions of the swept CRM65 wing model were obtained in NASA Icing Research Tunnel and compared. For most test conditions, the overall size and shape of the ice shapes compared well. On the main part of the ice shape, the hybrid model had larger chordwise extent, resulting in wider ice shape. For the scallop ice shapes, the ice shape heights at attachment line location were identical. For conditions that produced rime shape, hybrid model ice shapes were larger overall, even at the attachment line location. The only location where there was more ice on the full chord model was on the lower surface, downstream of where the hybrid model no longer had full-scale leading edge. For the test conditions where ice mass measurements were acquired for both models, the full-chord models had 11-22% less ice mass than on the hybrid wing. These differences were larger and opposite of those predicted from LEWICE 3D ice prediction calculations. Further investigation is needed in order determine the cause of the differences between the full-chord and hybrid models.

Icing↗

Comparison of Bimodal Ice Shapes on Full-Chord and Hybrid CRM Midspan Models

An ice shape database has been created to document ice accretions on a full-chord length and hybrid versions of a 65 percent scale size of the mid-span section of the Common Research Model resulting from an exposure to a Super-cooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each ice accretion. Drop distribution effects were evaluated by using the same IRT icing conditions except with either monomodal or bimodal drop distributions. Ice shapes resulting from the bimodal distribution as well as from monomodal distributions of identical cloud droplet Median Volume Diameter (MVD) were obtained and compared. Results indicated that the ice shapes resulting from the monomodal and bimodal drop size distributions have similar 2D Maximum Combined Cross Section, or MCCS, ice shape around the model leading edge region and equivalent ice mass measurement per unit span to within the IRT SLD cloud calibration uncertainty, as expected with this approach, but the bimodal distributions have distinctly larger ice feathers that extended further back on the chord of the model. For most test conditions, the overall size and shape of the main ice accretions compared well between the full-chord and hybrid models.

Icing↗

Comparison of Bimodal Ice Shapes on Full-Chord and Hybrid CRM Midspan Models

An ice shape database has been created to document ice accretions on a full-chord length and hybrid versions of a 65 percent scale size of the mid-span section of the Common Research Model resulting from an exposure to a Super-cooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each ice accretion. Drop distribution effects were evaluated by using the same IRT icing conditions except with either monomodal or bimodal drop distributions. Ice shapes resulting from the bimodal distribution as well as from monomodal distributions of identical cloud droplet Median Volume Diameter (MVD) were obtained and compared. Results indicated that the ice shapes resulting from the monomodal and bimodal drop size distributions have similar 2D Maximum Combined Cross Section, or MCCS, ice shape around the model leading edge region and equivalent ice mass measurement per unit span to within the IRT SLD cloud calibration uncertainty, as expected with this approach, but the bimodal distributions have distinctly larger ice feathers that extended further back on the chord of the model. For most test conditions, the overall size and shape of the main ice accretions compared well between the full-chord and hybrid models

Icing , Simulation, Validation↗

Comparison of Bimodal Ice Shapes on Full-Chord and Hybrid CRM Midspan Models

An ice shape database has been created to document ice accretions on a full-chord length and hybrid versions of a 65 percent scale size of the mid-span section of the Common Research Model resulting from an exposure to a Super-cooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each ice accretion. Drop distribution effects were evaluated by using the same IRT icing conditions except with either monomodal or bimodal drop distributions. Ice shapes resulting from the bimodal distribution as well as from monomodal distributions of identical cloud droplet Median Volume Diameter (MVD) were obtained and compared. Results indicated that the ice shapes resulting from the monomodal and bimodal drop size distributions have similar 2D Maximum Combined Cross Section, or MCCS, ice shape around the model leading edge region and equivalent ice mass measurement per unit span to within the IRT SLD cloud calibration uncertainty, as expected with this approach, but the bimodal distributions have distinctly larger ice feathers that extended further back on the chord of the model. For most test conditions, the overall size and shape of the main ice accretions compared well between the full-chord and hybrid models.

Icing↗

Rheology of the Ronne Ice Shelf, Antarctica, Inferred from Satellite Radar Interferometry Data using an Inverse Control Method

The Antarctic Ice Sheet is surrounded by large floating ice shelves that spread under their own weight into the ocean. Ice shelf rigidity depends on ice temperature and fabrics, and is influenced by ice flow and the delicate balance between bottom and surface accumulation. Here, we use an inverse control method to infer the rigidity of the Ronne Ice Shelf that best matches observations of ice velocity from satellite radar interferometry. Ice rigidity, or flow law parameter B, is shown to vary between 300 and 900 kPa a(sup 1/3). Ice is softer along the side margins due to frictional heating, and harder along the outflow of large glaciers, which advect cold continental ice. Melting at the bottom surface of the ice shelf increases its rigidity, while freezing decreases it. Accurate numerical modelling of ice shelf flow must account for this spatial variability in mechanical characteristics.

ice velocity↗

Introduction to Aircraft Icing and NASA's Approach to Understanding It

This presentation provides an introduction to aircraft icing and describes NASA’s approach to understanding it. The presentation covers various topics on icing. A general overview of different types of icing is provided, which includes airframe icing, jet engine icing, and rotorcraft icing. The presentation discusses NASA’s approach to studying the various forms of aircraft icing, which includes conducting flight tests, performing ground-based icing wind tunnel tests, and developing computational tools that simulate the icing physics. An overview on the future of aviation and icing is provided. The presentation is largely high level, but also provides some technical icing physics discussion.

Aircraft Icing↗

Results of the second Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+)

Ocean-driven basal melting of Antarctic ice shelves plays an important role in the mass loss of the Antarctic Ice Sheet. Ice shelf cavity-resolving ocean models are a valuable tool for understanding ice shelf-ocean interactions and for simulating projections of ice shelf and ocean states under future climate. Designed to assess the current state of ice shelf–ocean modelling, the second Ice Shelf–Ocean Model Intercomparison Project, ISOMIP+, consists of 12 ocean model configurations submitted with a common, idealised experimental setup. Here, we focus on the experiments Ocean0–2 (Asay-Davis et al., 2016), which are ocean models with idealised, static ice shelf geometries, but where the ocean reaches a balance with prescribed far-field ocean conditions. Different thermal transfer coefficient values (ranging from 0.011 to 0.2) are used for each model in the melting parameterisation to achieve a common, tuned melt rate since the models cover a range of types of vertical coordinates, ice–ocean boundary layer treatments, and numerical schemes. These model differences lead to spread in the resultant ocean properties, circulation, boundary-layer structure and spatial distribution of melting. We also highlight similarities between models, such as a shared linear relationship across most models between melt rate and overturning and barotropic streamfunctions during the spin-up and spin-down, demonstrating a robust relationship between melt and circulation across models and forcing conditions. The ISOMIP+ results provide a systematic comparison of ice shelf cavity-capable ocean models. However, we also demonstrate the need for realistic ice shelf–ocean model intercomparison projects (some already underway) to assess model biases and inter-model variation against sparse observations. Further research is needed to understand the differences between models and further improve our modelled representations of the ice–ocean boundary layer and ice shelf cavity circulation.

54 ENVIRONMENTAL SCIENCES↗

Sea ice and surface water circulation, Alaskan Continental Shelf

The author has identified the following significant results. The boundaries of land-fast ice, distribution of pack ice, and major polynya were studied in the vicinity of the Bering Strait. Movement of pack ice during 24 hours was determined by plotting the distinctly identifiable ice floes on ERTS-1 imagery obtained from two consecutive passes. Considerably large shallow area along the western Seward Peninsula just north of the Bering Strait is covered by land fast ice. This ice hinders the movement of ice formed in eastern Chukchi Sea southward through the Bering Strait. The movement of ice along the Russian coast is relatively faster. Plotting of some of the ice floes indicated movement of ice in excess of 30 km in and south of the Bering Strait between 6 and 7 March, 1973. North of the Bering Strait the movement approached 18 km. The movement of ice observed during March 6 and 7 considerably altered the distribution and extent of polynya. These features when continually plotted should be of considerable aid in navigation of ice breakers. The movement of ice will also help delineate the migration and distribution of sea mammals.

Wright, F. F.↗

Multilayer sea ice in the Arctic - Model- and satellite-derived

The behavior of sea ice is a strong function of ice strength which, in the central Arctic, depends on the relative fractions of multiyear and first-year ice coverage. A dynamic-thermodynamic sea ice model has been configured to include both multiyear and first-year ice. The model has been run for a 7-year period, 1979-1985, for which a satellite microwave record of multiyear sea ice coverage is available. The simulated concentrations of multiyear ice are slightly higher than the satellite-derived concentrations, although the total coverage of simulated multiyear ice agrees well with the satellite-derived coverage. Both the simulated and the satellite-derived coverages in the early 1980s show large interannual variations that represent potentially large perturbations in the mass balance of Arctic pack ice. These fluctuations are generally consistent in the model and satellite data, and they are generally consistent with year-to-year differences in the large-scale wind forcing. In at least some cases the wintertime distribution of multiyear ice is influenced by the pattern of ice drift as far back as the previous summer. The late December coverage of multiyear ice in the Alaskan region appears to foreshadow summer ice conditions in the same region. The short-term response to an intense synoptic system is considerably larger in the satellite data than in the model results.

Walsh, John E.↗

Equatorial ground ice on Mars: Steady-state stability

Current Martian equatorial surface temperatures are too warm for water ice to exist at the surface for any appreciable length of time before subliming into the atmosphere. Subsurface temperatures are generally warmer still and, despite the presence of a diffusive barrier of porous regolith material, it has been shown by Smoluchowski, Clifford and Hillel, and Fanale et al. that buried ground ice will also sublime and be lost to the atmosphere in a relatively short time. We investigate the behavior of this subliming subsurface ice and show that it is possible for ice to maintain at a steady-state depth, where sublimation and diffusive loss to the atmosphere is balanced by resupply from beneath by diffusion and recondensation of either a deeper buried ice deposits or ground water. We examine the behavior of equatorial ground ice with a numercial time-marching molecular diffusion model. In our model we allow for diffusion of water vapor through a porous regolith, variations in diffusivity and porosity with ice content, and recondensation of sublimed water vapor. A regolith containing considerable amounts of ice can still be very porous, allowing water vapor to diffuse up from deeper within the ice layer where temperatures are warmer due to the geothermal gradient. This vapor can then recondense nearer to the surface where ice had previously sublimed and been lost to the atmosphere. As a result we find that ice deposits migrate to find a steady-state depth, which represents a balance between diffusive loss to the atmosphere through the overlying porous regolith and diffusive resupply through a porous icy regolith below. This depth depends primarily on the long-term mean surface temperature and the nature of the geothermal gradient, and is independent of the ice-free porosity and the regolith diffusivity. Only the rate of loss of ground ice depends on diffusive properties.

Mellon, Michael T.↗

The distribution of ground ice on Mars

A wealth of geologic evidence indicates that subsurface water ice has played an important role in the evolution of Martian landforms. Theoretical models of the stability of ground ice show that in the near-surface regolith ice is currently stable at latitudes poleward of about +/- 40 deg and below a depth of a few centimeters to a meter, with some variations with longitude. If ice is not previously present at a particular location where it is stable, atmospheric water will diffuse into the regolith and condense as ice, driven by the annual subsurface thermal oscillations. The lower boundary of this ice deposit is found to occur at a depth (typically a few meters) where the annual thermal oscillations give way to the geothermal gradient. In the equatorial regions near-surface ice is currently not stable, resulting in the sublimation of any existing ice and subsequent loss to the atmosphere. However, subliming ice might be maintained at a steady-state depth, where diffusion and loss to the atmosphere are balanced by resupply from a possible deeper source of water (either deeper, not yet depleted, ice deposits or ground water). This depth is typically a few tens to hundreds of meters and depends primarily on the surface temperature and the nature of the geothermal gradient, being deeper for a higher surface temperature and a lower geothermal gradient. Such an equatorial deposit is characterized by the regolith ice content being low nearer the surface and increasing with depth in the deposit. Oscillations in the orbit will affect this picture of ground ice in two ways: by causing periodic changes in the pattern of near-surface stability and by producing subsurface thermal waves that may be capable of driving water ice deeper into the regolith.

Mellon, M. T.↗

Role of Wind Tunnels and Computer Codes in the Certification and Qualification of Rotorcraft for Flight in Forecast Icing

The cost and time to certify or qualify a rotorcraft for flight in forecast icing has been a major impediment to the development of ice protection systems for helicopter rotors. Development and flight test programs for those aircraft that have achieved certification or qualification for flight in icing conditions have taken many years, and the costs have been very high. NASA, Sikorsky, and others have been conducting research into alternative means for providing information for the development of ice protection systems, and subsequent flight testing to substantiate the air-worthiness of a rotor ice protection system. Model rotor icing tests conducted in 1989 and 1993 have provided a data base for correlation of codes, and for the validation of wind tunnel icing test techniques. This paper summarizes this research, showing test and correlation trends as functions of cloud liquid water content, rotor lift, flight speed, and ambient temperature. Molds were made of several of the ice formations on the rotor blades. These molds were used to form simulated ice on the rotor blades, and the blades were then tested in a wind tunnel to determine flight performance characteristics. These simulated-ice rotor performance tests are discussed in the paper. The levels of correlation achieved and the role of these tools (codes and wind tunnel tests) in flight test planning, testing, and extension of flight data to the limits of the icing envelope are discussed. The potential application of simulated ice, the NASA LEWICE computer, the Sikorsky Generalized Rotor Performance aerodynamic computer code, and NASA Icing Research Tunnel rotor tests in a rotorcraft certification or qualification program are also discussed. The correlation of these computer codes with tunnel test data is presented, and a procedure or process to use these methods as part of a certification or qualification program is introduced.

Flemming, Robert J.↗

A double-moment multiple-phase four-class bulk ice scheme, Part 1: Description

A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Ferrier, Brad Shoenberg↗

Backscatter from ice growing on shallow tundra lakes near Barrow, Alaska, winter 1991-1992

The timing of freeze-up and break-up of Arctic lake ice is a potentially useful environmental indicator that could be monitored using SAR. In order to do this, it is important to understand how the properties and structure of the ice during its growth and decay affect radar backscatter and thus lake ice SAR signatures. The availability of radiometrically and geometrically calibrated digital SAR data time series from the Alaska SAR Facility has made it possible for the first time to quantify lake ice backscatter intensity (sigma(sup o)) variations. This has been done for ice growing on shallow tundra lakes near Barrow, NW Alaska, from initial growth in September 1991 until thawing and decay in June 1992. Field and laboratory observations and measurements of the lake ice were made in late April 1992. The field investigations of the coastal lakes near Barrow confirmed previous findings that, (1) ice frozen to the lake bottom had a dark signature in SAR images, indicating weak backscatter, while, (2) ice that was floating had a bright signature, indicating strong backscatter. At all sites, regardless of whether the ice was grounded or floating, there was a layer of clear, inclusion-free ice overlaying a layer of ice with dense concentrations of vertically oriented tubular bubbles. At some sites, there was a third layer of porous, snow-ice overlaying the clear ice.

Jeffries, M. O.↗

The Role of Sea Ice in 2 x CO2 Climate Model Sensitivity: Hemispheric Dependencies - Part 2

How sensitive are doubled CO2 simulations to GCM control-run sea ice thickness and extent? This issue is examined in a series of 10 control-run simulations with different sea ice and corresponding doubled CO2 simulations. Results show that with increased control-run sea ice coverage in the Southern Hemisphere, temperature sensitivity with climate change is enhanced, while there is little effect on temperature sensitivity of (reasonable) variations in control-run sea ice thickness. In the Northern Hemisphere the situation is reversed: sea ice thickness is the key parameter, while (reasonable) variations in control-run sea ice coverage are of less importance. In both cases, the quantity of sea ice that can be removed in the warmer climate is the determining factor. Overall, the Southern Hemisphere sea ice coverage change had a larger impact on global temperature, because Northern Hemisphere sea ice was sufficiently thick to limit its response to doubled CO2, and sea ice changes generally occurred at higher latitudes, reducing the sea ice-albedo feedback. In both these experiments and earlier ones in which sea ice was not allowed to change, the model displayed a sensitivity of -0.02 C global warming per percent change in Southern Hemisphere sea ice coverage.

Rind, D.↗

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗

Elevation Change of the Southern Greenland Ice Sheet from Satellite Radar Altimeter Data

Long-term changes in the thickness of the polar ice sheets are important indicators of climate change. Understanding the contributions to the global water mass balance from the accumulation or ablation of grounded ice in Greenland and Antarctica is considered crucial for determining the source of the about 2 mm/yr sea-level rise in the last century. Though the Antarctic ice sheet is much larger than its northern counterpart, the Greenland ice sheet is more likely to undergo dramatic changes in response to a warming trend. This can be attributed to the warmer Greenland climate, as well as a potential for amplification of a global warming trend in the polar regions of the Northern Hemisphere. In collaboration with Drs. Curt Davis and Craig Kluever of the University of Missouri, we are using data from satellite radar altimeters to measure changes in the elevation of the Southern Greenland ice sheet from 1978 to the present. Difficulties with systematic altimeter measurement errors, particularly in intersatellite comparisons, beset earlier studies of the Greenland ice sheet thickness. We use altimeter data collected contemporaneously over the global ocean to establish a reference for correcting ice-sheet data. In addition, the waveform data from the ice-sheet radar returns are reprocessed to better determine the range from the satellite to the ice surface. At JPL, we are focusing our efforts principally on the reduction of orbit errors and range biases in the measurement systems on the various altimeter missions. Our approach emphasizes global characterization and reduction of the long-period orbit errors and range biases using altimeter data from NASA's Ocean Pathfinder program. Along-track sea-height residuals are sequentially filtered and backwards smoothed, and the radial orbit errors are modeled as sinusoids with a wavelength equal to one revolution of the satellite. The amplitudes of the sinusoids are treated as exponentially-correlated noise processes with a time-constant of six days. Measurement errors (e.g., altimeter range bias) are simultaneously recovered as constant parameters. The corrections derived from the global ocean analysis are then applied over the Greenland ice sheet. The orbit error and measurement bias corrections for different missions are developed in a single framework to enable robust linkage of ice-sheet measurements from 1978 to the present. In 1998, we completed our re-evaluation of the 1978 Seasat and 1985-1989 Geosat Exact Repeat Mission data. The estimates of ice thickness over Southern Greenland (south of 72N and above 2000 m) from 1978 to 1988 show large regional variations (+/-18 cm/yr), but yield an overall rate of +1.5 +/- 0.5 cm/yr (one standard error). Accounting for systematic errors, the estimate may not be significantly different from the null growth rate. The average elevation change from 1978 to 1988 is too small to assess whether the Greenland ice sheet is undergoing a long-term change.

Haines, Bruce J.↗