Search

This research investigates the use of Interferometric Synthetic Aperture Radar (InSAR) to generate a time-series of snow water equivalent (SWE) for dry snow within the Kuparuk watershed, North Slope, Alaska, during the winter of 1993/1994. Maps depicting relative change in phase and the theoretical relative change in SWE between satellite acquisitions are created for 3-day periods at the end of March 1994 using both ascending and descending ERS-1 overpasses. The theoretical coefficient relating relative change in phase and relative change in SWE for C-band is found to be at least twice as large as what is expected when using a simple single-layer snow model for this study area and time period. Without any direct measurements of SWE on the ground, station measurements of snow depth and hourly wind are linked to each 3-day relative change in phase map. Along with a qualitative assessment, quantitative measures of the rate and magnitude of phase change around these stations are directly compared to the hourly wind data for a given 3-day period. InSAR-derived maps acquired around a measured precipitation event show a considerable relationship to the predominant direction of strong winds over each 3-day period while maps acquired around no measureable precipitation depict much less correlation between phase change and predominant direction of strong winds. Despite limited ground measurements to infer snowpack conditions, these results show continued promise for the InSAR technique to measure changes in snowpack conditions (e.g. SWE) at much higher resolutions than manual sampling methods or passive microwave remote sensing. The extension of this technique to current L-band InSAR satellite platforms is also discussed.

The spatiotemporal variability of subsidence around the East Mesa Geothermal Field (EMGF) near the All American Canal (AAC) has been measured using 30 temporally averaged interferograms from 1992 to 2000. Deformation rate maps from two shorter time periods indicated the maximum subsidence rate of the EMGF was reduced from -43 mm year-1 (1992-1997) to -34 mm year-1 (1996-2000) corresponding to decreasing net geothermal water production. The maximum subsidence on the East Highline Canal was -9.5 ± 0.5 and -2.4 ± 0.8 cm for each shorter time period. Interferometric synthetic aperture radar (InSAR) stacking demonstrated its utility in monitoring subsidence of the canal caused by the nearby geothermal plant at regional coverage superior to ground levelling networks. Such data on the subsidence of surface and subsurface hydrodynamics along the US-Mexico border are scarce, and are particularly significant in a zone of present and likely future acute water resource management sensitivity.

Investigations into glacial changes, including understanding variations in the rates of glacial volume and surface-elevation changes, have increased over the past decade. This study uses historical glacier elevation data in the form of topographic maps from 1970 and a digital elevation model from the mid-1980s to calculate surface-elevation and volume changes for Ha-Iltzuk Icefield, southwest British Columbia, Canada. Results indicate that the icefield thinned at an average rate of 0.76±0.25 ma–1 during this period. A previous study of Ha-Iltzuk Icefield also using the geodetic method found a thinning rate of 1.0±0.20ma–1 between the mid-1980s and 1999, indicating a slight increase in the amount of icefield thinning. Within the ablation zone, thinning increased with decreasing elevation at a rate of 1.9±0.68 ma–1 km–1 between these two periods (1970 to mid-1980s versus mid-1980s to 1999). Analysis of meteorological data suggests that increases in both temperature and rainfall, as well as decreases in snowfall, likely contributed to the increased thinning rate.

Increased mass losses from the Greenland ice sheet and inferred contributions to sea-level rise have heightened the need for hydrologic observations of meltwater exiting the ice sheet. We explore whether temporal variations in ice-sheet surface hydrology can be linked to the development of a downstream sediment plume in Kangerlussuaq Fjord by comparing: (1) plume area and suspended sediment concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery and field data; (2) ice-sheet melt extent from Special Sensor Microwave/Imager (SSM/I) passive microwave data; and (3) supraglacial lake drainage events from MODIS. Results confirm that the origin of the sediment plume is meltwater release from the ice sheet. Interannual variations in plume area reflect interannual variations in surface melting. Plumes appear almost immediately with seasonal surface-melt onset, provided the estuary is free of landfast sea ice. A seasonal hysteresis between melt extent and plume area suggests late-season exhaustion in sediment supply. Analysis of plume sensitivity to supraglacial events is less conclusive, with 69% of melt pulses and 38% of lake drainage events triggering an increase in plume area. We conclude that remote sensing of sediment plume behavior offers a novel tool for detecting the presence, timing and interannual variability of meltwater release from the ice sheet.

This paper describes satellite data collected as part of the 2002/03 Cold Land Processes Experiment (CLPX). These data include multispectral and hyperspectral optical imaging, and passive and active microwave observations of the test areas. The CLPX multispectral optical data include the Advanced Very High Resolution Radiometer (AVHRR), the Landsat Thematic Mapper/Enhanced Thematic Mapper Plus (TM/ETM+), the Moderate Resolution Imaging Spectroradiometer (MODIS), and the Multi-angle Imaging Spectroradiometer (MISR). The spaceborne hyperspectral optical data consist of measurements acquired with the NASA Earth Observing-1 (EO-1) Hyperion imaging spectrometer. The passive microwave data include observations from the Special Sensor Microwave Imager (SSM/I) and the Advanced Microwave Scanning Radiometer (AMSR) for Earth Observing System (EOS; AMSR-E). Observations from the Radarsat synthetic aperture radar and the SeaWinds scatterometer flown on QuikSCAT make up the active microwave data.

This study adjusts and compares digital elevation models (DEMs) created from photogrammetric and interferometric synthetic aperture radar techniques to determine volume and surface elevation changes of five icefields in a remote region of southwest British Columbia, Canada, between the mid-1980s and 1999. Preliminary differences between the DEMs in ice-free and vegetation-free areas indicated variable elevation offsets with increasing altitude (11 m km−1) and with increasing slope (2.7 m (10°)−1). Results indicate a surface elevation change of −6.0 ± 2.7 m (−0.5 ± 0.2 m a−1) and a total volume loss of −19.4 ± 8.8 km3 (−1.5 ± 0.7 km3 a−1), which represents a potential sea-level rise contribution of 0.004 ± 0.002 mm a−1. Temperature and snowfall data from four nearby meteorological stations indicate that increased temperatures and decreased snowfall throughout the late 1980s and 1990s are a likely cause of the thinning. Glacier terminus positions were compared between a historical map (1927) and satellite images (1974, 1990/91 and 2000/01). All observed glaciers were in retreat between 1927 and 1974, as well as between 1990/91 and 2000/01, but many glaciers advanced or significantly slowed in their retreat between 1974 and 1990/91.

Thermally incised meltwater channels that flow each summer across melt-prone surfaces of the Greenland ice sheet have received little direct study. We use high-resolution WorldView-1/2 satellite mapping and in situ measurements to characterize supraglacial water storage, drainage pattern, and discharge across 6,812 km(2) of southwest Greenland in July 2012, after a record melt event. Efficient surface drainage was routed through 523 high-order stream/river channel networks, all of which terminated in moulins before reaching the ice edge. Low surface water storage (3.6 ± 0.9 cm), negligible impoundment by supraglacial lakes or topographic depressions, and high discharge to moulins (2.54-2.81 cm⋅d(-1)) indicate that the surface drainage system conveyed its own storage volume every

Large variations in both snow water equivalent (SWE) and snow slope stability are known to exist in the alpine snowpack, caused by wind, topographic and microclimatic effects. This variability makes extrapolation of point measurements of snowpack properties difficult and prone to error, but these types of measurements are used to estimate SWE and stability across entire mountain ranges. Radar technology provides a promising alternative to point measurements, because large areas can be covered quickly and non-intrusively. There is great potential for obtaining information on a large spatial scale from airborne applications. Frequency-modulated continuous wave (FMCW) radar measurements were made from the ground in several different alpine snowpacks, along with manual and in situ electrical measurements. The surface and ground reflections from the radar data, combined with an average density estimate, can provide a useful estimate of SWE. In addition, the locations of internal reflections are highly correlated with both visually identified layers and measured changes in in situ dielectric properties.

This research investigates the utility of passive microwave remote sensing instruments to accurately determine snow water equivalent (SWE) over large spatial extents. Three existing Special Sensor Microwave Imager (SSM/I) snow water equivalent algorithms produced by Chang, Tait and Goodison were evaluated for their ability to determine snow water equivalent in a snowpack containing substantial depth hoar, large faceted snow crystals. The Kuparuk River Watershed (8140 km2) test site on the North Slope of Alaska was chosen for its snowpack containing a think depth hoar layer and long history of ground truth data. A new regional snow water equivalent algorithm was developed to determine if it could produce better results than the existing algorithms in an area known to contain significant depth hoar. The four algorithms were tested to see how well they could determine snow water equivalent: (1) on a per pixel basis, (2) across swath-averaged spatial bands of approximately 850 km2, and (3) on a watershed scale. The algorithms were evaluated to see if they captured the annual spatial distribution in snow water equivalent over the watershed. Results show that the algorithms developed by Chang and from this research are generally within 3 cm of the spatially averaged snow water equivalents over the entire watershed. The algorithms produced by Chang, Tait, and in this research were able to predict the basin-wide ground measured snow water equivalent value within a percent error range from −32.4% to 24.4% in the years with a typical snowpack. None of the algorithms produce accurate results on a pixel-by-pixel scale, with errors ranging from −26% to 308%.

Ground-based frequency-modulated continuous wave (FMCW) radar measurements made throughout the 2002-03 NASA Cold Lands Processes (CLPX) mission in Colorado, USA, were designed to study the major electromagnetic transitions within the alpine snowpack over a wide range of conditions, and their effect on measurements made with different active radar measurement parameters. Measurements during 2002 determined that measurements at C-, X-, and Ku-band frequencies were necessary to retrieve the most snowpack information in a wide range of conditions. Measurements at four different incidence angles indicated that the snowpack layering was still visible at 15° but that the rough surface scattering of the snow-ground interface dominated the signal above 30° incidence. Measurements during 2003 were focused on characterizing the microwave response at C-, X-, and Ku-band frequencies at four different sites with different snowpack conditions, indicating that the optimal measurement parameters vary with snowpack conditions. Measurements at different bandwidths illustrate the effect of bandwidth on vertical resolution. This ground-truth data will help interpretation of air- and space-borne active and passive microwave radar measurements that were made coincident with this study. In addition, this work may help guide future researchers when choosing FMCW radar measurement parameters, depending on the type of snowpack information their study requires.

The Intensity–Hue–Saturation (IHS) transformation is used to integrate the high spectral resolution, provided by hyperspectral data (Airborne Visible Infrared Imaging Spectrometer, AVIRIS), and the surface texture information, derived from radar data (Topographic Synthetic Aperture Radar, TOPSAR), into a single image of an urban area. This transformed image is superimposed on the Digital Elevation Model (DEM) data derived from TOPSAR data to create a 3D perspective view. The ambiguity of several urban land cover types is resolved to a larger degree using the higher spectral and spatial resolutions and the synergistic visual content provided by the fused image in a 3D perspective. For urban areas that are at risk from geological hazards (e.g., avalanches, mudflows, and debris flows), the fused image provides a cost-effective product, rich in the information necessary for assessment and mitigation of these hazards.

Seward Glacier, on the Alaskan/Yukon border along the Gulf of Alaska, sits atop an important structural and morphological junction in the Saint Elias orogen. It is situated at the intersection between the Fairweather and Bagley strike–slip faults, and in the hanging wall of the Malaspina and Chugach–Saint Elias thrust faults. An ice surface velocity map of Seward Glacier derived from interferometric synthetic aperture (InSAR) reveals a complex flow pattern, which implies there is a previously unmapped northwest-trending supra-/subsurface ridge crossing the Seward. Analysis of additional remote-sensing images, ASTER, ERS SAR and the InSAR coherence map, confirms this observation. The presence of this ridge leads to a set of tectonic models describing the possible interaction of the underlying faults.

Daily acquisitions from satellite microwave sensors can be used to observe the spatial and temporal characteristics of the Arctic sea-ice snowmelt onset because the initial presence of liquid water in a dry snowpack causes a dramatic change in the active-and passive-microwave response. A daily sequence of backscatter coefficient images from the NASA scatterometer (NSCAT) clearly shows the spatially continuous progression of decreasing backscatter associated with snowmelt onset across the Arctic Ocean during spring 1997. A time series of the active NSCAT backscatter and a scattering index from the passive Special Sensor Microwave/Imager (SSM/I) show similar trends during the time of the melt onset. An NSCATsnowmelt-onset detection algorithm is developed using the derivative of the backscatter with respect to time to select a melt-onset date for each pixel, generating a melt map for the Arctic sea ice. Comparison between this melt map and one previously generated from an SSM/I scattering index shows the NSCAT algorithm predicts the onset occurs 1−10 days earlier than the SSM/I-based algorithm for most portions of multi-year ice.

Repeat-pass L-band interferometric synthetic aperture radar (InSAR) data for part of Hielo Patagónico Sur, Chile, were collected by the space-shuttle-based Spaceborne Imaging Radar C (SIR-C) over a 4 day span in October 1994. Three co-registered complex SAR images are used to generate phase-coherence maps, a digital elevation model (DEM) and an ice-velocity map. The phase-coherence maps indicate low coherence in the 5–15 km approaching the termini due to large displacements, ice deformation and melting. However, the coherence is high over nearly all of the remaining imaged icefield. Ice-velocity precision is greater than 2 cm d−1, while the DEM is good to about 25 m. A flow divide between two of the glaciers is mapped by locating a narrow band of near-zero ice velocity. Horizontal ice-surface velocity profiles calculated along flowlines show there is a high degree of spatial variability reaching a peak value of 5.5 m d−1 located 3.5 km from the terminus of Glaciar Europa. Longitudinal strain rates along the center lines calculated from these velocities at the locations of the initiation of crevassing are used to compute the tensile strength of ice (169–224 kPa).

C-band radar backscatter generally seems to be inversely related to accumulation rate across Greenland's dry snow zone. A coupled snow metamorphosis-backscatter model was constructed in order to investigate the observation. Accumulation rate and mean annual temperature were the coupled-model inputs. In a local area where accumulation rate is known very accurately the model predicts the backscatter level along several transects within the calibration uncertainties of the ERS-1 SAR data. The analysis was applied to the remainder of the dry snow zone. The model showed that both accumulation rate and temperature can modulate surface backscatter across Greenland's dry snow zone with accumulation rate being the primary influence.

The evolution of four dynamic radar glacier zones at the surface of an alpine icefield in British Columbia is monitored using a time series of 35 First European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) images acquired from 1992 to 1994. These zones result from changing wetness and textural properties, and appear to represent: (1) cold snow with no liquid water present; (2) an initial melt front with an upper boundary near the elevation of the 0° isotherm; (3) metamorphosed, rapidly melting first-year snow with a rough or pitted surface; and (4) bare ice. This interpretation is aided by temperature and runoff data, air photographs and field measurements of snowpack properties acquired with two ERS-1 SAR scenes, ice-surface elevations derived from 1:50 000 topographic maps and simulations of radar backscatter from a geometric optics model of surface scattering. Meltwater production is affected by the development of zones (2), (3) and (4), which form, migrate up-elevation and disappear each year between April and September.

The evolution of four dynamic radar glacier zones at the surface of an alpine icefield in British Columbia is monitored using a time series of 35 First European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) images acquired from 1992 to 1994. These zones result from changing wetness and textural properties, and appear to represent: (1) cold snow with no liquid water present; (2) an initial melt front with an upper boundary near the elevation of the 0° isotherm; (3) metamorphosed, rapidly melting first-year snow with a rough or pitted surface; and (4) bare ice. This interpretation is aided by temperature and runoff data, air photographs and field measurements of snowpack properties acquired with two ERS-1 SAR scenes, ice-surface elevations derived from 1:50 000 topographic maps and simulations of radar backscatter from a geometric optics model of surface scattering. Meltwater production is affected by the development of zones (2), (3) and (4), which form, migrate up-elevation and disappear each year between April and September.

The space-shuttle-based SIR-C/X-SAR synthetic aperture radar (SAR) imaged part of Hielo Patagónico Sur (HPS; southern Patagonia ice field, South America) for five successive days during missions in April and October 1994. A significant meteorological event occurred during each mission, including a major storm in April and a sharp temperature decrease in October. Changes in backscatter are observed for both episodes in X-SAR returns from the mid-portions of one of the two large outlet glaciers in the study area. Ground-station and satellite meteorological, and hydrological data are combined with the daily X-SAR images to interpret changes in glacier surface conditions caused by meteorological events. Effects interpreted from the April storm are (1) wind- and precipitation-influenced surface roughening of a wet snowpack, and (2) the deposition of new wet snow at lower elevation and its subsequent retreat up-glacier. An abrupt decrease in regional temperature during October is thought to reduce the snow wetness and increase grain-size. The changes in the radar-defined glacier zones due to the April precipitation event are subtle, while the October temperature drop causes significant backscatter increases. Our results suggest that trends in HPS glacier surface and near-surface conditions observable from spaceborne SARs are not significantly masked by precipitation events.

The first topographic and ice-motion maps of the northwestern flank of Hielo Patagónico Norte (HPN, northern Patagonia Icefield), in Chile, were produced using satellite synthetic-aperture interferometric radar data acquired by NASA’s Spaceborne Imaging Radar C instrument in October 1994. The topographic map has a 10 m vertical precision with a 30 m horizontal spacing, which should be sufficient to serve as a reference for monitoring future mass changes of the icefield. The ice-motion map is accurate to within 4 mm d−1 (or 1 m a−1). The radar-derived surface topography and ice velocity are used to estimate the ice discharge from the accumulation area of four outlet glaciers, and the calving flux and mass balance of Glaciar San Rafael. The results demonstrate the use of SAR interferometry for monitoring glaciological parameters on a spatial and temporal scale unattainable by any other means.

Interferometric radar observations of Glaciar San Rafael, Chile, were collected in October 1994 by NASA’s Spaceborne Imaging Radar C (SIR-C) at both L- (24 cm) and C-band frequency (5.6 cm), with vertical transmit and receive polarization. The C-band data did not yield good geophysical products, because the temporal coherence of the signal was significantly reduced alter 24 h. The L-band data were, however, successfully employed to map the surface topography of the icefield with a 10 m uncertainty in height, and measure ice velocity with a precision of 4 mm d−1 or 1.4 m a−1. The corresponding error in strain rates is 0.05 a −1 at a 30 m horizontal spacing. The one-dimensional interferometric velocities were subsequently converted to horizontal displacements by assuming a flow direction and complemented by feature-tracking results near the calving front. The results provide a comprehensive view of the ice-flow dynamics of Glaciar San Rafael. The glacier has a core of rapid flow, 4.5 km in width and 3.5° in average slope, surrounded by slower-moving ice, not by rock. Ice velocity is 2.6 m d−1 or 0.95 km a−1 near the equilibrium-line altitude (1200 m), increasing rapidly before the glacier enters the narrower terminal valley, to reach 17.5 md−1 or 6.4 k ma−1 at the calving front. Strain rates are dominated by lateral shearing at the glacier margins (0.4–0.7 a−1), except for the terminal-valley section, where longitudinal strain rates average close to 1 a−1. This spectacular longitudinal increase in ice velocity in the last few kilometers may be a fundamental feature of tidewater glaciers.

Using multitemporal ERS 1 synthetic aperture radar (SAR) satellite imagery and simultaneous ground measurements of streamflow, a strong correlation (R2 = 0.89) was found between water surface area and discharge for a braided glacial river in British Columbia, Canada. Satellite-derived effective width (We) was found to vary with discharge (Q) as We = 27.5Q0.42, where We is defined as the total water surface area within a 10 km × 3 km control section, divided by the section length. This “area/discharge rating curve” yields instantaneous discharge estimates with a mean error of ±275 m3/s for ground-measured flows that ranged from 242 to 6350 m3/s.

Our poor understanding of tidewater glacier dynamics remains the primary source of uncertainty in sea level rise projections. On the ice sheets, mass lost from tidewater calving exceeds the amount lost from surface melting. In Alaska, the magnitude of calving mass loss remains unconstrained, yet immense calving losses have been observed. With 20% of the global new-water sea level rise coming from Alaska, partitioning of mass loss sources in Alaska is needed to improve sea level rise projections. Here we present the first regionally comprehensive map of glacier flow velocities in Central Alaska. These data reveal that the majority of the regional downstream flux is constrained to only a few coastal glaciers. We find regional calving losses are 17.1 Gt a−1, which is equivalent to 36% of the total annual mass change throughout Central Alaska. Alaskan mountain glaciers are losing ice and contribute to sea level rise, but contributions from specific ice-loss mechanisms are not known. Here, calving losses in Central Alaska are found to equal 36% of the net regional mass change each year and regional flux is dictated largely by snow accumulation rates.

[1] Past estimates of Greenland Ice Sheet accumulation rates have been multiyear climatologies based on ice/firn cores and coastal precipitation records. Existing annually resolved estimates have incompletely quantified uncertainty, owing primarily to incomplete spatial coverage. This study improves upon these shortcomings by calibrating annual (1958–2007) solid precipitation output from the Fifth Generation Mesoscale Model modified for polar climates (Polar MM5) using firn core and meteorological station data. The calibration employs spatial interpolation of regionally derived linear correction functions. Residual uncertainties exhibit coherent spatial patterns, which are modeled via spatial interpolation of root mean squared errors. Mean 1958–2007 Greenland Ice Sheet annual accumulation rate is 337 ± 48 mm/yr water equivalent (w.e.) or 591 ± 83 Gt/yr. Annual estimates contain one standard deviation uncertainties of 74 mm/yr w.e., 22%, or 129 Gt/yr. Accumulation rates in southeast Greenland are found to exceed 2000 mm/yr w.e. and to dominate interannual variability in Greenland Ice Sheet total accumulated mass, representing 31% of the whole. Accumulation rates in the southeast are of sufficient magnitude to affect the sign of Greenland mass balance during some years. The only statistically significant temporal change in total ice sheet accumulation in the 1958–2007 period occurred between 1960 and 1972, when a simultaneous accumulation increase and decrease occurred in west and east Greenland, respectively. No statistically significant uniform change in ice sheet‐wide accumulation is evident after 1972. However, regional changesdooccur,includingan accumulation increaseonthewestcoast post‐1992.Thehigh accumulation rates of 2002–2003 appear to be confined to the southeast.

Global totals of ice area and volume for mountain glaciers are small compared to the Antarctic and Greenland ice sheets, however, their present contribution to sea level rise (SLR) is on par with the ice sheets. Alaska Glaciers provide the largest contribution of all mountain glaciers to SLR. We discuss the use of spaceborne remote sensing techniques to study Alaska glacier contributions to sea level and suspected changes in ice velocity. Direct calculation of ice volume changes and rates of change on Kenai Peninsula glaciers are based on differences in USGS and SRTM DEMs and airborne Lidar profiles. A preliminary study of ice velocities on the Seward Glacier using L-band InSAR from PALSAR is presented.

[1] Temperate mountain glaciers are thinning at high rates and significantly contributing to sea level rise. Due to these glaciers' remote locations, remote sensing and digital elevation models (DEMs) are an effective way to calculate their thinning rates and contribution to sea level rise. Comparisons of 1950s United States Geological Survey (USGS) and 2000 Shuttle Radar Topographic Mission (SRTM) DEMs with DEMs produced from Lidar profiles of glacier center-lines indicate thinning rates from the mid-1990s to 1999 (−0.72 ± 0.13 m y−1) accelerated by a factor of 1.5 as compared with 1950 to mid-1990s (−0.47 ± 0.01 m y−1) for glaciers on the Harding Icefield. Overall, comparison of USGS and SRTM DEMs indicate the Harding Icefield and Grewingk-Yalik Glacier Complex, Alaska, are thinning −0.61 ± 0.12 m y−1 from 1950 to 1999.

[1] An important component of ice sheet mass balance calculations is the snow accumulation rate, traditionally determined across the ice sheet from interpolated ice core data. We derive an accumulation rate map within the Greenland dry-snow facies, using C-band radar backscatter and the mean annual temperature. Values for the backscatter were obtained from the ERS-1 synthetic aperture radar mosaic 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 (5.3 GHz) 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, consistent with the Greenland dry-snow facies. On the basis of our modeled results, a look-up table is generated, which uniquely maps the measured radar backscatter and mean annual temperature to accumulation rate. Our results compare favorably (15.4% average difference) with 124 in situ accumulation rate measurements falling within our study area. We attribute some of the differences to (1) unreliable in situ data points, which were based on a single year's accumulation and (2) a systematic underestimate of accumulation rate by our forward model in the low-accumulation regions associated with the northern portion of our map. With questionable data points and those above 77°N removed, the average differences were less than 10.3%, as determined from 73 ice cores.

Ice sheet backscatter around South Pole Station was measured from several different azimuthal angles using RADARSAT-1 SAR data. The authors observed strong azimuthal anisotropy (3-4 dB) between radar look directions across and along the prevailing wind direction. They speculate that the azimuthal variation of backscatter around the South Pole area may be largely attributed to systematically oriented sastrugi and buried sastrugi. Radar backscatter is well correlated with surface topography over a spatial wavelength band of 15 to 20 km. They propose that grain size variations control backscatter patterns associated with undulating surface topography.

Snow facies, or melt zones, on the Greenland Ice Sheet are correlated with brightness patterns on spaceborne radar imagery. To better relate radar scattering signatures to individual zones, the authors investigate the relationship between radar scattering and firn physical properties. They operated surface-based radar systems at five different sites on the Greenland Ice Sheet. In the dry snow zone the angular variation of the backscattering coefficient indicates that at high incidence angles the return energy is dominated by volume scattering. They also observe a systematically higher backscattering coefficient when the physical temperature is low and all the other parameters are identical. In the wet snow zone the backscatter coefficient is temporally dependant as the temperature variations affect the moisture content and modifies the spectral backscattered energy. The snow and ice surface roughnesses determine the angular response. In the percolation zone, the backscatter coefficient is independent of the frequency and the angular response is also controlled by the ice structures (surface roughness, lenses, and pipes).

The authors obtained in-situ radar backscatter data at several sites on the Greenland Ice Sheet to better relate scattering signatures to firn physical properties. Backscatter measurements were made using an ultra-wide band radar system (.5-18 GHz) for incidence angles between 0 and 50 degrees, with respect to the snow surface. For large angles of incidence (>25/spl deg/) the measured backscatter is dominated by volumetric scattering, as opposed to surface and specular backscatter occurring at near normal incidence. Using the data the authors present a method for estimating the volumetric backscatter coefficient as a function of depth in firn. Results show trends consistent with both seasonal variations and long-term grain growth, which is primarily influenced by the accumulation rate.

An X-band FM-CW radar was used to determine the feasibility of observing annual snow-accumulation layers in Antarctica with a high-resolution inexpensive radar system. The formation of layering boundaries, their resultant electromagnetic discontinuity and their detection by reflected energy are presented. Large returns from depths corresponding to reasonable positions for annual layers were found. The average accumulation rates calculated from the radar returns agree with those measured in a previous pit study done in the same area. The detection of the annual accumulation layers with this system implies a simple, inexpensive mobile radar could be used to profile large areas allowing the distorting effects of local topography to be removed.This type of system with a concurrent pit study could provide insight into the effect of sub-surface strata on spaceborne or airborne microwave remote sensing.

An interferometric method is used to derive ice motion from RADARSAT data collected during the Antarctic Mapping Mission. Although one cannot solve for both topography and ice motion using one interferometric pair, it is possible to use a coarsely sampled digital terrain model to estimate ice motion using an image registration method. Less accurate than the usual fringe counting method for estimation of radial displacement, the image registration method allows useful motion estimation in both range and azimuth. The method is described and some results shown for a large area (/spl sim/17000 km/sup 2/) including ice flow into the Filchner Ice Shelf.

The entire continent of Antarctica was mapped at a 25 meter resolution with synthetic aperture radar (SAR) during the Radarsat Antarctic Mapping Project (RAMP) over a 30 day period in the fall of 1997 providing a static "snapshot" of the ice sheet. Since Radarsat-l has a 24-day orbit cycle, repeat-pass interferometric SAR (InSAR) data was also acquired. The extensive InSAR data will provide a view of ice sheet kinematics. The authors present the results of analysis of an InSAR pair for the Recovery Glacier, East Antarctica. Analysis of the Recovery Glacier data set has demonstrated the excellent quality of InSAR data available from the Radarsat Antarctic Mapping Mission and its application to glacier dynamics. Phase coherence is high over most this highly dynamic area even though the temporal baseline is 24 days. The coherence image enhances subtle ice sheet topography and indicates the location of shear margins and the grounding line. A radar line-of sight ice velocity map identified a large region of slow moving ice that drains the Shackelton Range. Longitudinal stain rates along a flow line within the main ice stream were calculated directly from the wrapped phase.

Improved understanding of Greenland ice sheet hydrology is critically important for assessing its impact on current and future ice sheet dynamics and global sea level rise. This has motivated the collection and integration of in situ observations, model development, and remote sensing efforts to quantify meltwater production, as well as its phase changes, transport, and export. Particularly urgent is a better understanding of albedo feedbacks leading to enhanced surface melt, potential positive feedbacks between ice sheet hydrology and dynamics, and meltwater retention in firn. These processes are not isolated, but must be understood as part of a continuum of processes within an integrated system. This letter describes a systems approach to the study of Greenland ice sheet hydrology, emphasizing component interconnections and feedbacks, and highlighting research and observational needs.

The northwest directed motion of the Pacific plate is accompanied by migration and collision of the Yakutat terrane into the cusp of southern Alaska. The nature and magnitude of accretion and translation on upper crustal faults and folds is poorly constrained, however, due to pervasive glaciation. In this study we used high-resolution topography, geodetic imaging, seismic, and geologic data to advance understanding of the transition from strike-slip motion on the Fairweather fault to plate margin deformation on the Bagley fault, which cuts through the upper plate of the collisional suture above the subduction megathrust. The Fairweather fault terminates by oblique-extensional splay faulting within a structural syntaxis, allowing rapid tectonic upwelling of rocks driven by thrust faulting and crustal contraction. Plate motion is partly transferred from the Fairweather to the Bagley fault, which extends 125 km farther west as a dextral shear zone that is partly reactivated by reverse faulting. The Bagley fault dips steeply through the upper plate to intersect the subduction megathrust at depth, forming a narrow fault-bounded crustal sliver in the obliquely convergent plate margin. Since . 20 Ma the Bagley fault has accommodated more than 50 km of dextral shearing and several kilometers of reverse motion along its southern flank during terrane accretion. The fault is considered capable of generating earthquakes because it is linked to faults that generated large historic earthquakes, suitably oriented for reactivation in the contemporary stress field, and locally marked by seismicity. The fault may generate earthquakes of Mw

The extent to which northern peatlands respond to or influence climate change is an unresolved question in Arctic science. Recent studies in Alaska, Canada, and Fennoscandia have raised concerns that northern peatlands, while currently a net sink or minor source of atmospheric CO2 , may become a significant CO2 source under a warming climate. Expanding peatlands emit methane but sequester atmospheric carbon through long-term accumulation of undecomposed plant matter. Drier conditions may reverse this process by increasing temperatures and lowering the peatland water table, causing anaerobic decomposition of stored peat and subsequent outgassing of CO2.While this process would likely reduce methane emissions and possibly enhance C uptake from increased soil nutrient mineralization rates [Oechel and Vourlitis, 1994], many scientists now believe that warming and drying of northern peatlands will liberate stored C for uptake by the atmosphere and biosphere.

The Luneburg lens, which is commonly used to calibrate radar scatterometers, may not be appropriate for calibration of dual-antenna systems. Theory and results from measurements of the Luneburg lens reradiation pattern indicate that gross errors may be incurred when using the lens to calibrate a dual-antenna system. Measurements of the reradiation pattern indicate that special attention must be paid to the angular separation of the antennas when using the lens to calibrate a dual-antenna radar system. The angular separation of the antennas must be less than the half-power beamwidth of the lens to ensure accurate repeatable measurements. A formula for estimating the half-power beamwidth is given. >