Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice sheets”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

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.↗
NASA NTRS1993

Effect of subglacial volcanism on changes in the West Antarctic Ice Sheet

Rapid changes in the West Antarctic Ice Sheet (WAIS) may affect future global sea-level changes. Alley and Whillans note that 'the water responsible for separating the glacier from its bed is produced by frictional dissipation and geothermal heat,' but assume that changes in geothermal flux would ordinarily be expected to have slower effects than glaciological parameters. I suggest that episodic subglacial volcanism and geothermal heating may have significantly greater effects on the WAIS than is generally appreciated. The WAIS flows through the active, largely asiesmic West Antarctic rift system (WS), which defines the sub-sea-level bed of the glacier. Various lines of evidence summarized in Behrendt et al. (1991) indicate high heat flow and shallow asthenosphere beneath the extended, weak lithosphere underlying the WS and the WAIS. Behrendt and Cooper suggest a possible synergistic relation between Cenozoic tectonism, episodic mountain uplift and volcanism in the West Antarctic rift system, and the waxing and waning of the Antarctic ice sheet beginning about earliest Oligocene time. A few active volcanoes and late-Cenozoic volcanic rocks are exposed throughout the WS along both flanks, and geophysical data suggest their presence beneath the WAIS. No part of the rift system can be considered inactive. I propose that subglacial volcanic eruptions and ice flow across areas of locally (episodically?) high heat flow--including volcanically active areas--should be considered possibly to have a forcing effect on the thermal regime resulting in increased melting at the base of the ice streams.

Behrendt, John C.↗
NASA NTRS2021

Surface Melting Over the Greenland Ice Sheet Derived From Enhanced Resolution Passive Microwave Brightness Temperatures (1979–2019)

Surface melting is a major component of the Greenland ice sheet surface mass balance, and it affects sea level rise through direct runoff and the modulation of ice dynamics and hydrological processes, supraglacially, englacially and subglacially. Passive microwave (PMW) brightness temperature observations are of paramount importance in studying the spatial and temporal evolution of surface melting due to their long temporal coverage (1979–present) and high temporal resolution (daily). However, a major limitation of PMW datasets has been the relatively coarse spatial resolution, which has historically been of the order of tens of kilometers. Here, we use a newly released PMW dataset (37 GHz, horizontal polarization) made available through a NASA “Making Earth System Data Records for Use in Research Environments” (MeASUREs) program to study the spatiotemporal evolution of surface melting over the Greenland ice sheet at an enhanced spatial resolution of 3.125 km. We assess the outputs of different detection algorithms using data collected by automatic weather stations (AWSs) and the outputs of the Modèle Atmosphérique Régional (MAR) regional climate model. We found that sporadic melting is well captured using a dynamic algorithm based on the outputs of the Microwave Emission Model of Layered Snowpack (MEMLS), whereas a fixed threshold of 245 K is capable of detecting persistent melt. Our results indicate that, during the reference period from 1979 to 2019 (from 1988 to 2019), surface melting over the ice sheet increased in terms of both duration, up to 4.5 (2.9) d per decade, and extension, up to 6.9 % (3.6 %) of the entire ice sheet surface extent per decade, according to the MEMLS algorithm. Furthermore, the melting season started up to 4.0 (2.5) d earlier and ended 7.0 (3.9) d later per decade. We also explored the information content of the enhanced-resolution dataset with respect to the one at 25 km and MAR outputs using a semi-variogram approach. We found that the enhanced product is more sensitive to local-scale processes, thereby confirming the potential of this new enhanced product for monitoring surface melting over Greenland at a higher spatial resolution than the historical products and for monitoring its impact on sea level rise. This offers the opportunity to improve our understanding of the processes driving melting, to validate modeled melt extent at high resolution and, potentially, to assimilate these data in climate models.

Surface melting↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
NASA NTRS1995

Assessment of climate variability of the Greenland ice sheet: Integration of in situ and satellite data

The proposed research involves the application of multispectral satellite data in combination with ground truth measurements to monitor surface properties of the Greenland ice sheet which are essential for describing the energy and mass of the ice sheet. Several key components of the energy balance are parameterized using satellite data and in situ measurements. The analysis will be done for a ten year time period in order to get statistics on the seasonal and interannual variations of the surface processes and the climatology. Our goal is to investigate to what accuracy and over what geographic areas large scale snow properties and radiative fluxes can be derived based upon a combination of available remote sensing and meterological data sets. Data analysis showed the following results: (1)cloud classification based on longwave sky radiation revealed that overcast sky occurred for 25% of the time in winter, and for 15% in spring and summer respectively (winter and summer both show the same occurrence of clear sky of approximately 26%); (2) comparison of aerodynamic profile method with eddy correlation method to derive sensible and latent heat flux showed good agreement in the diurnal cycle and the turbulent fluxes were underestimated with the aerodynamic method by 10 - 30% as compared to the in situ eddy flux method; (3) the katabatic wind shows a distinct diurnal cycle with a maximum in the morning (7-9 h solar time) and a minimum in the later afternoon (18 h solar time); (4) snow grain size was modeled with a surface energy balance model (SNTHERM) and compared with in situ measurements. Sharp decreases in the modeled snow grain size, caused by accumulation events such as precipitation and deposition, could be verified with observational data; (4) radiative transfer modeling of firn supports our beliefs that the observed trends in 18 and 19 GHz passive microwave brightness temperatures are attributable to accumulation rate changes (modeling also indicates the above relationship is detectable because of the presence of depth hoar; (5) snow melt can be detected by a distinct signal in the passive microwave cross-polarized gradient ratio (19h-37v)/(19h+37v) and has been used for wet/dry snow classification; (6) top of the atmosphere (TOA) broadband albedos were derived from AVHRR visible and near infrared reflectances for the entire ice sheet from May 1990 - June 1991, and the highest albedo values are found along the southeast coast of the ice sheet which is consistent with the summer peak of precipitation due to onshore flow loaded with high water vapor content (TOA albedo values dropped to around 40% along the south-western coast during July and August due to bare ice surfaces); and (7) the net all-wave radiation balance at the top of the atmosphere is negative over the entire ice sheet except for the summer months - June-July-August (in June, the net radiation balance is slightly positive over the dry snow areas--15 W/m2).

Steffen, K.↗
NASA NTRS2001

Determination of Interannual to Decadal Changes in Ice Sheet Mass Balance from Satellite Altimetry

A major uncertainty in predicting sea level rise is the sensitivity of ice sheet mass balance to climate change, as well as the uncertainty in present mass balance. Since the annual water exchange is about 8 mm of global sea level equivalent, the +/- 25% uncertainty in current mass balance corresponds to +/- 2 mm/yr in sea level change. Furthermore, estimates of the sensitivity of the mass balance to temperature change range from perhaps as much as - 10% to + 10% per K. Although the overall ice mass balance and seasonal and inter-annual variations can be derived from time-series of ice surface elevations from satellite altimetry, satellite radar altimeters have been limited in spatial coverage and elevation accuracy. Nevertheless, new data analysis shows mixed patterns of ice elevation increases and decreases that are significant in terms of regional-scale mass balances. In addition, observed seasonal and interannual variations in elevation demonstrate the potential for relating the variability in mass balance to changes in precipitation, temperature, and melting. From 2001, NASA's ICESat laser altimeter mission will provide significantly better elevation accuracy and spatial coverage to 86 deg latitude and to the margins of the ice sheets. During 3 to 5 years of ICESat-1 operation, an estimate of the overall ice sheet mass balance and sea level contribution will be obtained. The importance of continued ice monitoring after the first ICESat is illustrated by the variability in the area of Greenland surface melt observed over 17-years and its correlation with temperature. In addition, measurement of ice sheet changes, along with measurements of sea level change by a series of ocean altimeters, should enable direct detection of ice level and global sea level correlations.

Zwally, H. Jay↗
NASA NTRS1978

An investigation of the astronomical theory of the ice ages using a simple climate-ice sheet model

The astronomical theory of the Quaternary ice ages is incorporated into a simple climate model for global weather; important features of the model include the albedo feedback, topography and dynamics of the ice sheets. For various parameterizations of the orbital elements, the model yields realistic assessments of the northern ice sheet. Lack of a land-sea heat capacity contrast represents one of the chief difficulties of the model.

Pollard, D.↗
NASA NTRS1997

A digital elevation model of the Greenland ice sheet and validation with airborne laser altimeter data

A 2.5 km resolution digital elevation model (DEM) of the Greenland ice sheet was produced from the 336 days of the geodetic phase of ERS-1. During this period the altimeter was operating in ice-mode over land surfaces providing improved tracking around the margins of the ice sheet. Combined with the high density of tracks during the geodetic phase, a unique data set was available for deriving a DEM of the whole ice sheet. The errors present in the altimeter data were investigated via a comparison with airborne laser altimeter data obtained for the southern half of Greenland. Comparison with coincident satellite data showed a correlation with surface slope. An explanation for the behavior of the bias as a function of surface slope is given in terms of the pattern of surface roughness on the ice sheet.

Bamber, Jonathan L.↗
NASA NTRS2010

Development of a Climate-Data Record (CDR) of the Surface Temperature of the Greenland Ice Sheet

Regional "clear sky" surface temperature increases since the early 1980s in the Arctic, measured using Advanced Very High Resolution Radiometer (AVHRR) infrared data, range from 0.57+/-0.02 deg C to 72+/-0.10 deg C per decade. Arctic warming has important implications for ice-sheet mass balance because much of the periphery of the Greenland Ice Sheet is already near 0 deg C during the melt season, and is thus vulnerable to rapid melting if temperatures continue to increase. An increase in melting of the ice sheet would accelerate sea-level rise, an issue affecting potentially billions of people worldwide. To quantify the ice-surface temperature (IST) of the Greenland Ice Sheet, and to provide an IST dataset of Greenland for modelers that provides uncertainties, we are developing a climate-data record (CDR) of daily "clear-sky" IST of the Greenland Ice Sheet, from 1982 to the present using AVHRR (1982 - present) and Moderate-Resolution Imaging Spectroradiometer (MODIS) data (2000 - present) at a resolution of approximately 5 km. Known issues being addressed in the production of the CDR are: time-series bias caused by cloud cover (surface temperatures can be different under clouds vs. clear areas) and cross-calibration in the overlap period between AVHRR instruments, and between AVHRR and MODIS instruments. Because of uncertainties, mainly due to clouds, time-series of satellite IST do not necessarily correspond with actual surface temperatures. The CDR will be validated by comparing results with automatic-weather station data and with satellite-derived surface-temperature products and biases will be calculated.

Hall, Dorthy K.↗
NASA NTRS2016

The Abandoned Ice Sheet Base at Camp Century, Greenland, in a Warming Climate

In 1959 the U.S. Army Corps of Engineers built Camp Century beneath the surface of the northwestern Greenland Ice Sheet. There they studied the feasibility of deploying ballistic missiles within the ice sheet. The base and its wastes were abandoned with minimal decommissioning in 1967, under the assumption they would be preserved for eternity by perpetually accumulating snowfall. Here we show that a transition in ice sheet surface mass balance at Camp Century from net accumulation to net ablation is plausible within the next 75years, under a business-as-usual anthropogenic emissions scenario (Representative Concentration Pathway 8.5). Net ablation would guarantee the eventual remobilization of physical, chemical, biological, and radiological wastes abandoned at the site. While Camp Century and four other contemporaneous ice sheet bases were legally established under a Danish-U.S. treaty, the potential remobilization of their abandoned wastes, previously regarded as sequestered, represents an entirely new pathway of political dispute resulting from climate change.

ice sheet↗
NASA NTRS1993

The diatom record from beneath the West Antarctic Ice Sheet and the global proxy perspective

Recent glaciological evaluation and modeling of the marine-based West Antarctic Ice Sheet (WAIS) support the possibility that the WAIS disintegrated during one or more Pleistocene interglacial period(s). The magnitude of sea level and oxygen isotope variation during certain late-Pleistocene interglacial periods is also consistent with the possibility of major retreat of the WAIS. Although oxygen isotopes from deep-sea sediments provide the best available proxy record for global ice volume (despite the ambiguities in the record), the source of ice volume changes must be hypothesized. Based on the intensity of interglacial isotopic shifts recorded in Southern Ocean marine sedimentary records, stage 11 (400,000 years ago) is the strongest candidate for WAIS collapse, but the records for stages 9, 7, and 5.5 are all consistent with the possibility of multiple late-Pleistocene collapses. Seismic reflection studies through the WAIS have revealed thick successions of strata with seismic characteristics comparable to upper Tertiary marine sediments. Small samples of glacial diamictons from beneath the ice sheet have been collected via hot-water drilled access holes. These sediments include mixed diatom assemblages of varying ages. Late-Miocene diatoms dominate many samples, probably reflecting marine deposition in West Antarctic basins prior to development of a dominantly glacial phase in West Antarctica. In addition to late-Miocene diatoms, samples from Upstream B (1988/89) contain rare post-Miocene diatoms, many of which imply deposition in the West Antarctic interior during one or more Pleistocene deglaciation periods. Age-diagnostic fossils in glacial sediments beneath ice sheets provide relatively coarse chronostratigraphic control, but they do contain direct evidence of regional deglaciation. Thus, sub-glacial till samples provide the evidence regarding the source of ice sheet variability seen in well-dated proxy records. Combined, these independent data sets can provide a more comprehensive and less speculative interpretation of the history of past glacial minima in currently glaciated polar regions.

Scherer, Reed P.↗
NASA NTRS1991

Marine record of late quaternary glacial-interglacial fluctuations in the Ross Sea and evidence for rapid, episodic sea level change due to marine ice sheet collapse

Some of the questions to be addressed by SeaRISE include: (1) what was the configuration of the West Antarctic ice sheet during the last glacial maximum; (2) What is its configuration during a glacial minimum; and (3) has it, or any marine ice sheet, undergone episodic rapid mass wasting. These questions are addressed in terms of what is known about the history of the marine ice sheet, specifically in Ross Sea, and what further studies are required to resolve these problems. A second question concerns the extent to which disintegration of marine ice sheets may result in rises in sea level that are episodic in nature and extremely rapid, as suggested by several glaciologists. Evidence that rapid, episodic sea level changes have occurred during the Holocene is also reviewed.

Anderson, John B.↗
NASA NTRS2008

Recent Ice Sheet and Glacier Elevation Changes in Greenland from Aircraft Laser Altimetry

The Arctic Ice Mapping group (Project AIM) at the NASA Goddard Space Flight Center Wallops Flight Facility has been conducting systematic topographic surveys of the Greenland Ice Sheet (GIS) since 1993, using scanning airborne laser altimeters combined with Global Positioning System (UPS) technology. Earlier surveys showed the ice sheet above 2000-rn elevation to be in balance, but with localized regions of thickening or thinning. Thinning predominates at lower elevations and thinning rates have recently increased, resulting in a negative mass balance for the entire ice sheet. Recently, critical segments of near-coastal flight lines in Greenland were resurveyed. Results from the new data will be presented.

Krabill, William B.↗
NASA NTRS1989

Can Milankovitch orbital variations initiate the growth of ice sheets in a general circulation model?

The possibility of initiating the growth of ice sheets by solar insolation variations is examined. The study is conducted using a climate model with three different orbital configurations corresponding to 116,000 and 106,000 yr before the present and a modified insolation field with greater reductions in summer insolation at high northern latitudes. Despite the reduced summer and fall insolation, the model fails to maintain snow cover through the summer at locations of suspected ice sheet initiation. The results suggest that there is a discrepancy between the model's response to Milankovitch perturbations and the geophysical evidence of ice sheet initiation. If the model results are correct, the growth of ice shown by geophysical evidence would have occurred in an extremely ablative environment, demanding a complicated strategy.

Rind, D.↗
NASA NTRS2013

Insights into Spatial Sensitivities of Ice Mass Response to Environmental Change from the SeaRISE Ice Sheet Modeling Project I: Antarctica

Atmospheric, oceanic, and subglacial forcing scenarios from the Sea-level Response to Ice Sheet Evolution (SeaRISE) project are applied to six three-dimensional thermomechanical ice-sheet models to assess Antarctic ice sheet sensitivity over a 500 year timescale and to inform future modeling and field studies. Results indicate (i) growth with warming, except within low-latitude basins (where inland thickening is outpaced by marginal thinning); (ii) mass loss with enhanced sliding (with basins dominated by high driving stresses affected more than basins with low-surface-slope streaming ice); and (iii) mass loss with enhanced ice shelf melting (with changes in West Antarctica dominating the signal due to its marine setting and extensive ice shelves; cf. minimal impact in the Terre Adelie, George V, Oates, and Victoria Land region of East Antarctica). Ice loss due to dynamic changes associated with enhanced sliding and/or sub-shelf melting exceeds the gain due to increased precipitation. Furthermore, differences in results between and within basins as well as the controlling impact of sub-shelf melting on ice dynamics highlight the need for improved understanding of basal conditions, grounding-zone processes, ocean-ice interactions, and the numerical representation of all three.

SeaRISE↗
NASA NTRS1995

Observations of ice-sheet motion in Greenland using satellite radar interferometry

We present interferometric observations of ice-sheet motion in western Greenland based on pairs of ERS-1 synthetic aperture radar (SAR) images. These observations provide the first detailed regional view of ice motion associated with dynamically supported topography near the margin of an ice sheet. The interferograms of this area are much more complicated than other interferograms of ice sheets presented to date. We devote the largest part of this paper to explaining the source of the complexity in these interferograms. A synthetic interferogram based on a simple model helps to illustrate the effects of different components of the ice velocity field in interferometric data and suggests a method for estimating the large-scale ice velocity field from such data.

Joughin, Ian R.↗
NASA NTRS2021

Effect of Ice Sheet Thickness on Formation of the Hiawatha Impact Crater

The discovery of a large putative impact crater buried beneath Hiawatha Glacier along the margin of the northwestern Greenland Ice Sheet has reinvigorated interest into the nature of large impacts into thick ice masses. This circular structure is relatively shallow and exhibits a small central uplift, whereas a peak-ring morphology is expected. This discrepancy may be due to long-term and ongoing subglacial erosion but may also be explained by a relatively recent impact through the Greenland Ice Sheet, which is expected to alter the final crater morphology. Here we model crater formation using hydrocode simulations, varying pre-impact ice thickness and impactor composition over crystalline target rock. We find that an ice-sheet thickness of 1.5 or 2 km results in a crater morphology that is consistent with the present morphology of this structure. Further, an ice sheet that thick substantially inhibits ejection of rocky material, which might explain the absence of rocky ejecta in most existing Greenland deep ice cores if the impact occurred during the late Pleistocene. From the present morphology of the putative Hiawatha impact crater alone, we cannot distinguish between an older crater formed by a pre-Pleistocene impact into ice-free bedrock or a younger, Pleistocene impact into locally thick ice, but based on our modeling we conclude that latter scenario is possible.

Elizabeth A. Silber↗
NASA NTRS2002

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.↗