mass balance of east antarctic glaciers and ice shelves from satellite data thevelocityandmassdischarge of ninemajoreast antarcticglaciersnot draining into the ross or filchner ronne ice shelves is investigated using interferometric syntheticapertureradar insar datafromtheeuropeanremote sensingsatellite1and2 ers 1 2 andradarsat 1 theglaciersare david ninnis mertz totten scott denman lambert shiraseandstancomb wills insarisusedtolocatetheirgroundinglinewithpre cision ice velocity is measured with either insar or a speckle tracking technique ice thickness is deduced from prior determined ice shelf elevation assuming hydrostatic equi librium mass fluxes are calculated both at the grounding line and at a flux gate located downstream the grounding line flux is compared to a mass input calculated from snow accumulation to deduce the glacier mass balance the calculation is repeated at the flux gate downstream of the grounding line to estimate the average bottom melt rate of the ice shelf under steady state conditions the main results are 1 grounding lines are found several tens of km upstream of prior identified positions not because of a recent ice sheet retreatbutbecause of the inadequacy of prior determined grounding line positions 2 no gross imbalance between outflow and inflow is detected on the nine glaciersbeing investi gated with anuncertainty of 10 20 prior determined largely positive mass imbalances wereduetoanincorrectlocalizationofthegroundingline 3 highratesofbottommelting 24 7 micea 1 are inferred neargrounding zones where ice reaches the deepestdraft a few glaciers exhibit lower bottom melt rates 4 7 micea 1 bottom melting however appearsto be a major source of mass loss on antarctic ice shelves

tide model accuracy in the amundsen sea antarctica from radar interferometry observations of ice shelf motion 1 this study assesses the accuracy of tide model predictions in the amundsen sea sector of west antarctica tide model accuracy in this region is poorly constrained yet tide models contribute to simulations of ocean heat transfer and to the removal of tidal signals from satellite observations of ice shelves we use two satellite based interferometric synthetic aperture radar insar methods to measure the tidal motion of the dotson ice shelf at multiple epochs a single difference technique that measures tidal displacement and a double difference technique that measures changes in tidal displacement we use these observations to evaluate predictions from three tide models tpxo7 1 cats2008a opt and fes2004 all three models perform comparably well exhibiting root mean square deviations from the observations of 9 cm single difference technique and 10 cm double difference technique care should be taken in generalizing these error statistics because 1 the dotson ice shelf experiences relatively small semidiurnal tides and 2 our observations are not sensitive to all tidal constituents an error analysis of our insar based methods indicates measurement errors of 7 and 4 cm for the single and double difference techniques respectively a model based correction for the effect of fluctuations in atmospheric pressure yields an 6 improvement in the agreement between tide model predictions and observations this study suggests that tide model accuracy in the amundsen sea is comparable to other antarctic regions where tide models are better constrained these methods can be used to evaluate tide models in other remote antarctic waters

the response of black rapids glacier alaska to the denali earthquake rock avalanches 1 we describe the impact of three simultaneous earthquake triggered rock avalanches on the dynamics of black rapids glacier alaska by using spaceborne radar imagery and numerical modeling we determined the velocities of the glacier before and after landslide deposition in 2002 by using a combination of ers 1 ers 2 tandem radarsat 1 and alos palsar synthetic aperture radar data ice velocity above the debris covered area of the glacier increased up to 14 after the earthquake but then decreased 20 by 2005 within the area of the debris sheets mean glacier surface velocity increased 44 within 2 years of the landslides at the downglacier end of the lowest landslide where strong differential ablation produced a steep ice cliff velocities increased by 109 over the same period by 2007 ice velocity throughout the debris area had become more uniform consistent with a constant ice flux resulting from drastically reduced ablation at the base of the debris without further analysis we cannot prove that these changes resulted from the landslides because black rapids glacier displays large seasonal and interannual variations in velocity however a full stokes numerical ice flow model of a simplified glacier geometry produced a reversal of the velocity gradient from compressional to extensional flow after 5 years which supports our interpretation that the recent changes in the velocity field of the glacier are related to landslide induced mass balance changes

combining airborne photographs and spaceborne sar data to monitor temperate glaciers potentials and limits monitoring temperate glacier activity has become more and more necessary for economical and security reasons and as an indicator of the local effects of global climate change remote sensing data provide useful information on such complex geophysical objects but they require specific processing techniques to cope with the difficult context of moving and changing features in high relief areas this paper presents the first results of a project involving four laboratories developing and combining specific methods to extract information from optical and synthetic aperture radar sar data two different information sources are processed namely 1 airborne photography and 2 spaceborne c band sar interferometry the difficulties and limitations of their processing in the context of alpine glaciers are discussed and illustrated on two glaciers located in the mont blanc area the results obtained by aerial triangulation techniques provide digital terrain models with an accuracy that is better than 30 cm which is compatible with the computation of volume balance and useful for precise georeferencing and slope measurement updating the results obtained by sar differential interferometry using european remote sensing satellite images show that it is possible to measure temperate glacier surface velocity fields from october to april in one day interferograms with approximately 20 m ground sampling this allows to derive ice surface strain rate fields required to model the glacier flow these different measurements are complementary to results obtained during the summer from satellite optical data and ground measurements that are available only in few accessible points

decorrelation in interferometric radar echoes a radar interferometric technique for topographic mapping of surfaces implemented utilizing a single synthetic aperture radar sar system in a nearly repeating orbit is discussed the authors characterize the various sources contributing to the echo correlation statistics and isolate the term which most closely describes surficial change they then examine the application of this approach to topographic mapping of vegetated surfaces which may be expected to possess varying backscatter over time it is found that there is decorrelation increasing with time but that digital terrain model generation remains feasible the authors present such a map of a forested area in oregon which also includes some nearly unvegetated lava flows such a technique could provide a global digital terrain map

a full coverage high resolution topographic model of greenland computed from a variety of digital elevation data in this paper a new high resolution elevation model of greenland including the ice sheet as well as the ice free regions is presented it is the first published full coverage model computed with an average resolution of 2 km and providing an unprecedented degree of detail the topography is modeled from a wide selection of data sources including satellite radar altimetry from geosat and ers 1 airborne radar altimetry and airborne laser altimetry over the ice sheet and photogrammetric and manual map scannings in the ice free region the ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations the mean accuracy of the ice sheet elevations is estimated to be 12 13 m and it is found that on surfaces of a slope between 0 2 and 0 8 corresponding to approximately 50 of the ice sheet the model presents a 40 improvement over models based on satellite altimetry alone on coastal bedrock the model is compared with stereo triangulated reference points and it is found that the model accuracy is of the order of 25 35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere

three dimensional surface velocities of storstrommen glacier greenland derived from radar interferometry and ice sounding radar measurements non steady state vertical velocities of up to 5 m y 1 exceed the vertical surface parallel flow components over much of the ablation area of storstrommen a large outlet glacier from the east greenland ice sheet neglecting a contribution to the vertical velocity of this magnitude results in substantial errors up to 20 also on the south north component of horizontal velocities derived by satellite synthetic aperture radar interferometry insar measurements in many glacier environments the steady state vertical velocity component required to balance the annual ablation rate is 5 to 10 m y 1 or more this indicates that the surface parallel flow assumption may be problematic also for glaciers in steady state here we derive the three dimensional surface velocity distribution of storstrommen by using the principle of mass conservation to combine insar measurements from ascending and descending satellite tracks with airborne ice sounding radar measurement of ice thickness the results are compared to insar velocities previously derived by using the assumption of surface parallel flow and to velocities obtained by inin situ gps measurements the velocities derived by using the principle of mass conservation are in better agreement with the gps velocities than the previously calculated velocities derived with the assumption of surface parallel flow

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

multi track extraction of two dimensional surface velocity by the combined use of differential and multiple aperture insar in the amery ice shelf east antarctica abstract glaciers and ice streams drain the ice from the ice sheets to the ice shelves and contribute to global sea level changes therefore ice velocity is an essential parameter for the assessment of mass balance in this paper we propose an approach for the multi track extraction of two dimensional ice surface velocity by the combined use of multiple aperture interferometry mai and differential synthetic aperture radar interferometry dinsar the azimuth and ground range ice movements of each individual track are first estimated by the use of mai and dinsar respectively a seamless two dimensional ice velocity map is then achieved by calibration and mosaicking of the estimated velocity of the multi temporal and multi track measurements over a large area finally the dynamic behavior of the amery ice shelf is analyzed by formulating a long term series of ice flow velocities consisting of historical surface velocities 1968 1990 the estimated surface velocity by the proposed method with ers 1 2 data 1996 existing surface velocity products such as the radarsat 1 antarctic mapping project ramp velocities 1997 2000 and rignot s results 2007 2010 as well as the estimated velocities obtained by the use of the offset tracking method from landsat 8 images 2013 2014 the results show that the estimated velocity obtained by the use of the proposed method is superior to the other methods and the long term series of ice velocity measurements shows that the ice velocity in the amery ice shelf has hardly changed over the past 45 years 1968 2014 we conclude that the ice flow velocity of the amery ice shelf has showed a slight deceleration over the past 45 years

application of image cross correlation to the measurement of glacier velocity using satellite image data abstract image to image cross correlation software is applied to pairs of digital satellite images to map the velocity field of moving ice this technique uses small scale glacial surface features such as crevasse scars and snow dunes as markers on the surface of the moving ice displacements of the surface features are mapped by selecting small image areas centered on distinct features or by dividing a large area of densely featured glacial surface into a grid of areas and searching a subsequent image for matching areas using a cross correlation algorithm interpolation of the peak correlation values allows the displacements to be measured to subpixel accuracy resulting in very precise velocity measurements cross correlation is also applied to provide image coregistration in areas devoid of bedrock exposures in such areas subtle large scale topographic undulations in the ice surface related to underlying bedrock structure may be correlated by using large image areas and low pass filtered images both types of applications are demonstrated using ice stream d and ice stream e in west antarctica as test areas a high resolution map of the velocity field of the central portion of ice stream e generated by the displacement measuring technique is presented the use of cross correlation software is a significant improvement over previous manually based photogrammetric methods for velocity measurement and is far more cost effective than in situ methods in remote polar areas

satellite radar interferometry for monitoring ice sheet motion application to an antarctic ice stream satellite radar interferometry sri provides a sensitive means of monitoring the flow velocities and grounding line positions of ice streams which are indicators of response of the ice sheets to climatic change or internal instability the detection limit is about 1 5 millimeters for vertical motions and about 4 millimeters for horizontal motions in the radar beam direction the grounding line detected by tidal motions where the ice goes afloat can be mapped at a resolution of 0 5 kilometer the sri velocities and grounding line of the rutford ice stream antarctica agree fairly well with earlier ground based data the combined use of sri and other satellite methods is expected to provide data that will enhance the understanding of ice stream mechanics and help make possible the prediction of ice sheet behavior

identification and characterization of alpine subglacial lakes using interferometric synthetic aperture radar insar brady glacier alaska usa the temporary storage and subsequent release of water at glacial margins can cause severe flooding in downstream areas and substantially impact glacier dynamics alpine subglacial lakes may not be identified until they become subaerially exposed or release a jokulhlaup we use interferometric synthetic aperture radar insar to identify and characterize three dynamic alpine subglacial lakes of brady glacier alaska usa we quantify changes in vertical displacement of the glacier surface and lake volumes from september 1995 through march 1996 using european remote sensing satellite 1 2 ers 1 2 tandem data in the autumn subsidence ranged from 4 to 26 cm d 1 and the volume of water discharged ranged from 22 000 2000 to 243 000 14 000 m 3 d 1 subsidence and discharge rates declined significantly during the winter and continued at a lesser rate through march application of this technique may allow researchers to locate alpine subglacial lakes years or decades before they begin to release hazardous outburst floods and substantially impact glacier dynamics

estimation of glacier ice extinction using long wavelength airborne pol insar in the recent years there has been increased interest in using synthetic aperture radar sar to study and monitor glaciers and ice sheets for glaciological and climate change research this paper describes the estimation of ice extinctions through the modeling of polarimetric interferometric sar pol insar coherences as a combination of a surface contribution from the snow firn interface and wind induced sastrugi features and a volume response ground to volume scattering ratios derived from a novel polarimetric decomposition are used in conjunction with pol insar coherence magnitudes to invert the extinction of the ice layer the inversion is performed for experimental airborne pol insar data at l band and p band acquired by the german aerospace center s dlr e sar system over the austfonna ice cap in svalbard norway as part of the 2007 icesar campaign extinction dependences on frequency and glacier facie are investigated and validation is performed comparing ground penetrating radar data to sar backscatter and extinction values best extinction results are obtained at shallow incidence angles with small wavenumbers and for low ground to volume scattering ratios for swath areas in mid range to far range accuracies of 25 in extinction are anticipated when averaging over 100 effective looks for a four baseline inversion constraining solutions to vertical wavenumbers of 0 01 kz 0 1 to allow inversion using single baseline pol insar the proposed model is of limited complexity suggested extensions for a more realistic scattering scenario include incorporating multiple englacial ice scattering layers and improving the way multiple baselines are combined

ice velocity measurements of langjokull iceland from interferometric synthetic aperture radar insar we have computed the velocity of langjokull ice cap lic iceland using insar data acquired by the european remote sensing satellite in 1994 and an ice surface elevation model derived from a gps survey conducted in 1997 the velocity map covers 68 of the 925 km 2 of lic eight principal outlet glaciers were identified along which ice velocities increased from zero at the ice divide to a maximum of 75 m a 1 when the insar velocity dataset is compared with surface gps measurements of ice flow determined in summer 2001 there is generally good agreement suggesting interannual and seasonal stability for much of the ice cap during this period however non zero ice velocities at the termini of four of the eight glaciers surveyed suggest that these glaciers were advancing in 1994 comparison of our gps and insar velocities suggests that sudurjokull experienced elevated velocity in late february 1994 which when other studies are considered appears to be consistent with a glacier surge

satellite derived volume loss rates and glacier speeds for the juneau icefield alaska we provide a high resolution map of elevation change rates dh dt at the juneau icefield jif southeastern alaska in order to quantify its contribution to sea level rise between 2000 and 2009 2013 we also produce the first high resolution map of ice speeds at the jif which we use to constrain flux and look for acceleration we calculate dh dt using stacked digital elevation models dems from the advanced spaceborne thermal emission and reflection radiometer aster instrument and the shuttle radar topography mission srtm taking into account srtm c band penetration via comparison with srtm x band elevations overall the jif is losing mass less rapidly 0 13 0 12mw e a than other alaskan icefields 0 79mw e a we determine glacier speeds using pixel tracking on optical image pairs acquired from 2001 to 2010 by aster from radar image pairs acquired between 2007 and 2011 and from radar interferometry in 1995 we detect seasonal speed variations but no interannual acceleration ruling out dynamics as the cause of the observed thinning thinning must therefore be due to the documented warming in the region flux measurements confirm this for mendenhall glacier showing that calving constitutes only 2 5 5 of mass loss there

northern larsen ice shelf antarctica further retreat after collapse changes of larsen ice shelf antarctica and the surrounding glaciers after its collapse in 1995 were investigated using satellite radar imagery with emphasis onchanges in the glaciers which previously nourished the ice shelf north of seal nunataks the large glaciers retreated several kilometres inland of the previousgrounding line the velocity field of drygalski glacier the largest glacier in this area was mapped by means of interferograms derived from pairs of ers syntheticaperture radar images from 1995 and 1999 the main part of the glacier showed significant acceleration of flow over these 4 years with an increase of veocity upto three fold at the terminus similar accelerations were observed by means of interferometry on several other grounded glaciers suggesting that the removal ofice shelves could lead to an effect on eustatic sea level for larsen b the northernmost surviving part of larsen ice shelf the retreat of the ice front to october2000 is documented

estimating the mass change of mountain glacier using a novel method based on insar observations melting glaciers have a direct contribution to sea level rise runoff and glacier lake outburst flood disasters therefore accurate estimations of glacier mass balances with high spatial large region and temporal annual seasonal resolutions are very important however the estimation of glacial mass balances in mountainous regions is usually hampered by remoteness and the lack of high precision topographic data for most mountain ranges this study presents a novel method for estimating mass changes in mountain glaciers using insar data we utilised observations of glacier surface deformation to derive changes in thickness and then calculated the mass change this method can accurately estimate short term mass changes we apply this method to the koxkar glacier in the tien shan mountain range in china and successfully estimate the seasonal mass change using theoretical error and statistical error analysis we determined that the accuracy of this method is much better than other geodetic methods we analyse the spatial characteristics of the mass changes in the ablation zone for the first time the results show considerable spatial and seasonal variability which were mainly from ablation in the summer with the greatest amount of glacial ablation reaching up to 3872 mm water equivalent w e and from accumulation in the winter

greenland ice sheet surface properties and ice dynamics from ers 1 sar imagery c band synthetic aperture radar sar imagery from the european space agency s ers 1 satellite reveals the basic zonation of the surface of the greenland ice sheet the zones have backscatter signatures related to the structure of the snowpack which varies with the balance of accumulation and melt at various elevations the boundaries of zones can be accurately located with the use of this high resolution imagery the images also reveal a large flow feature in northeast greenland that is similar to ice streams in antarctica and may play a major role in the discharge of ice from the ice sheet

positive mass balance during the late 20th century on austfonna svalbard revealed using satellite radar interferometry determining whether increasing temperature or precipitation will dominate the cryo spheric response to climate change is key to forecasting future sea level rise the volume of ice contained in the ice caps and glaciers of the arctic archipelago of svalbard is small compared with that of the greenland or antarctic ice sheets but is likely to be affected much more rapidly in the short term by climate change this study investigates the mass balance of austfonna svalbard s largest ice cap equilibrium line fluxes for the whole ice cap and for individual drainage basins were estimated by combining surface velocities measured using satellite radar interferometry with ice thicknesses derived from radio echo sounding these fluxes were compared with balance fluxes to reveal that during the 1990s the total mass balance of the accumulation zone was 5 6 2 0 10 8 m 3 a 1 three basins in the quiescent phase of their surge cycles contributed 75 of this accumulation the remaining volume may be attributable either to as yet unidentified surge type glaciers or to increased precipitation this result emphasizes the importance of considering the surge dynamics of glaciers when attempting to draw any conclusions on climate change based on snapshot observations of the cryosphere

contribution to the glaciology of northern greenland from satellite radar interferometry interferometric synthetic aperture radar insar data from the ers 1 and ers 2 satellites are used to measure the surface velocity topography and grounding line position of the major outlet glaciers in the northern sector of the greenland ice sheet the mass output of the glaciers at and above the grounding line is determined and compared with the mass input we find that the grounding line output is approximately in balance with the input except for the three largest glaciers for which the mass loss is 4 3 km3 ice year 1 or 11 8 of the mass input along the coast we detect a systematic retreat of the grounding lines between 1992 and 1996 with insar which implies that the outlet glaciers are thinning the inferred coastal thinning is too large to be explained by a few warm summers glacier thinning must be of dynamic origin that is caused by spatial and temporal changes in ice velocity iceberg production from the glaciers is uncharacteristically low it accounts for only 8 of the ice discharge to the ocean about 55 of the ice is lost through basal melting 5 8 m ice year 1 on average from the underside of the floating glacier tongues that are in contact with warm ocean waters mass losses are highest in the first 10 km of floating ice where ice reaches the greatest depths and basal melting is 3 times larger than on average only a small increase in basal melting would suffice to disintegrate the floating glacier tongues copyright 2001 by the american geophysical union

measurement of ice sheet topography using satellite radar interferometry detailed digital elevation models dems do not exist for much of the greenland and antarctic ice sheets radar altimetry is at present the primary in many cases the only source of topographic data over the ice sheets but the horizontal resolution of such data is coarse satellite radar interferometry uses the phase difference between pairs of synthetic aperture radar sar images to measure both ice sheet topography and surface displacement we have applied this technique using ers 1 sar data to make detailed i e 80 m horizontal resolution maps of surface topography in a 100 km by 300 km strip in west greenland extending northward from just above jakobshavns isbrae comparison with a 76 km long line of airborne laser altimeter data shows that we have achieved a relative accuracy of 2 5m along the profile these observations provide a detailed view of dynamically supported topography near the margin of an ice sheet in the final section we compare our estimate of topography with phase contours due to motion and confirm our earlier analysis concerning vertical ice sheet motion and complexity in ers 1 sar interferograms

ice sheet motion and topography from radar interferometry both topography and motion information are present in repeat pass ers 1 interferograms over ice sheets the authors demonstrate that the topography is separable from the surface displacement field when a sequence of radar images are available if the velocity field is constant over the time span of observation the topography can be derived from differential interferograms formed from sequential observations with this measurement a pure displacement field can then be obtained by removal of the topographic contribution to the interferometric phase at each pixel further they discuss how the vertical and horizontal components of displacement affect the interferometrically derived motion field they illustrate their approach with four successive 3 day repeat ers 1 images of a flow feature in northeastern greenland

surface motion of mountain glaciers derived from satellite optical imagery a complete and detailed map of the ice velocity field on mountain glaciers is obtained by cross correlating spot5 optical images this approach offers an alternative to sar interferometry because no present or planned radar satellite mission provides data with a temporal separation short enough to derive the displacements of glaciers the methodology presented in this study does not require ground control points gcps the key step is a precise relative orientation of the two images obtained by adjusting the stereo model of one slave image assuming that the other master image is well georeferenced it is performed with numerous precisely located homologous points extracted automatically the strong ablation occurring during summer time on the glaciers requires a correction to obtain unbiased displacements the accuracy of our measurement is assessed based on a comparison with nearly simultaneous differential gps surveys performed on two glaciers of the mont blanc area alps if the images have similar incidence angles and correlate well the accuracy is on the order of 0 5 m or 1 5 of the pixel size similar results are also obtained without gcps an acceleration event observed in early august for the mer de glace glacier is interpreted in term of an increase in basal sliding our methodology applied to spot5 images can potentially be used to derive the displacements of the earth s surface caused by landslides earthquakes and volcanoes

ice motion of the patagonian icefields of south america 1984 2014 2015 american geophysical union all rights reserved we present the first comprehensive high resolution mosaic of ice velocity of the northern npi and southern patagonian icefields spi from multiple synthetic aperture radar and optical data collected between 1984 and 2014 the results reveal that many of the outlet glaciers extend far into the central ice plateaus which implies that changes in ice dynamics propagate far inside the accumulation area we report pronounced seasonal to interannual variability of ice motion on pio xi and jorge montt a doubling in speed of jorge montt a major slow down of o higgins significant fluctuations of upsala and a deceleration of san rafael which illustrate the need for sustained continuous time series of ice motion to understand the long term evolution of the rapidly thinning icefields the velocity product also resolves major ambiguities in glacier drainage in areas of relatively flat topography illustrating the need to combine topography and flow direction to map drainage basins

changes in west antarctic ice stream velocities observation and analysis 1 we have produced a map of velocity covering much of the siple coast ice streams the map confirms earlier estimates of deceleration on whillans ice stream comparison with bed elevation data indicates that subglacial topography and the location of consolidated sediment play a strong role in determining the location of the tributaries feeding the ice streams force balance estimates based on these data indicate that the tributaries have beds nearly an order of magnitude stronger than those beneath many of the ice streams we have used a theoretical analysis to examine the controls on fast flow this analysis suggests that ice plains very wide ice streams are inherently unstable this instability may be responsible for the current deceleration on the ice plain of whillans ice stream and the shutdown of ice stream c 150 years ago thinning induced reductions in driving stress may also explain some of the observed deceleration particularly in upstream areas the active portions of ice stream c coincide well with the areas where we estimate that melt should be taking place current topography and inferences of large thickening following a shutdown suggest the upstream migration of a stagnation front that initiated at the ice plain uncertainty remains about the basal conditions on ice stream d while the basal resistance on ice stream e is large enough to ensure basal melting

tidal motion ice velocity and melt rate of petermann gletscher greenland measured from radar interferometry over a floating glacier ice tongue or an ice shelf the glacier motion measured by a single repeat pass radar interferogram is difficult to analyze because the long term steady motion of the ice is intermixed with its cyclic downward motion induced by tidal forcing multiple interferograms and a quadruple difference technique are necessary to separate the tidal signal from the long term steady motion of the ice an example of application of this technique is given here using ers 1 radar images of petermann gletscher a major outlet glacier of northern greenland tidal displacements are measured with 5 mm statistical noise the long term ice velocity is measured with a precision of 1 m a 1 the inferred tidal displacements agree well with model predictions from a fixed elastic beam with an elastic damping factor of 0 47 0 01 km 1 the hinge line is mapped with a precision of 20 80 m combining the interferometric ice velocities with ice thickness data the glacier ice discharge is calculated at and below the hinge line at the hinge line the ice flux is 12 1 1 km3 a 1 at the ice front calf ice production is only 0 59 km3 a 1 meaning that 95 of the ice that crosses the grounding line melts before it reaches the calving front assuming steady state conditions the melt rate of the glacier tongue averages 12 1 m a 1 with peak values exceeding 20 m a 1 near the hinge line this high melt rate cannot be accommodated by surface ablation alone only about 2 3 m a 1 and is attributed to pronounced basal melting of the ice tongue basal melting often assumed to be negligible in greenland is the dominant process of mass release from the floating section of petermann gletscher

new velocity map and mass balance estimate of mertz glacier east antarctica derived from landsat sequential imagery automatic feature tracking on two landsat images acquired in january 2000 and december 2001 generates a complete velocity field of mertz glacier east antarctica this velocity field shows two main tributaries to the ice stream between the tributaries a likely obstruction feature in the bedrock results in a slow down of the flow a third landsat image acquired in 1989 and combined with the 2000 image permits the determination of the glacier mean velocity during the 1990s although some parts of the mertz glacier system show evidence of slight speed increase we conclude that the mertz flow speed is constant within our uncertainty 35 m a 1 using this complete velocity field new estimates of the ice discharge flux 17 8 km 3 a 16 4 gt a and of the basal melting of the tongue 11m a of ice are given our results lead to an apparent imbalance of the drainage basin ice discharge 3 5 km 3 a lower than the accumulation considering previous studies in the mertz glacier area we discuss the uncertainty of this imbalance and the problems with accumulation mapping for this region

investigating mountain glacier motion with the method of sar intensity tracking removal of topographic effects and analysis of the dynamic patterns abstract accurate knowledge of glacier motion variations can provide critical information for assessing glacier mass balance changes for fast flowing high mountain glaciers sar intensity tracking may be the most appropriate technique to measure their motions however due to the slant range imaging mode topographic relief can give rise to pixel disalignments in sar image and then bring errors to sar intensity tracking as a result the application of sar intensity tracking has been limited to sar pairs with short baselines because the use of image pairs with long baselines would suffer significant topographic effects over mountainous regions this paper presents an improvement to sar intensity tracking by adding an explicit procedure of topographic correction such correction will not only significantly increase the accuracy of the sar pairs with shorter baselines but also allow the image pairs with longer baselines to be taken advantage of instead of letting them to be dismissed otherwise using the improved sar intensity tracking technique we derive the highly precise motion variation patterns of 11 giant glaciers in tuomuer khan tengri mountain ranges one of the largest glacial centers in tien shan the results manifest that most glaciers there saw flow deceleration during the observation interval 6 january 2007 4 march 2011 except the northern and southern inylchek glacier regarding the glacier motion mechanism ice thermal character and climate features there we deduce that the flow acceleration and deceleration should be the sign of glacier mass gain and loss respectively and the sharp and slight flow decelerations should correspond to remarkable and mild mass loss respectively we confirm these deductions by exploring the close link between the long term motion variations and the dynamic patterns of these glaciers the obvious diversities of the dynamic patterns among the studied glaciers are ascribed to the discrepancies of the glaciers location and scale longer observations are still needed to confirm or redefine these findings in the future though provided that the precipitation keeps increasing whether the glaciers in tuomuer khan tengri mountain ranges will regain mass or not has not been determined yet

geodetically accurate insar data processor we present a new interferometric synthetic aperture radar insar processing approach that capitalizes on the precise orbit tracking that is available with modern radar satellites our method uses an accurate orbit information along with motion compensation techniques to propagate the radar echoes to positions along a noninertial virtual orbit frame in which the location and focusing equations are particularly simple so that images are focused without requiring autofocus techniques and are computed efficiently motion compensation requires two additional focus correction phase terms that are implemented in the frequency domain if the images from an interferometric pair or stack are all computed along the same reference orbit flat earth topographic correction is not needed and image coregistration is simplified obviating many difficulties that are often encountered in insar processing we process several data sets collected by the alos palsar instrument and find that the geodetic accuracy of the radar images is 10 20 m with up to 20 m of additional image distortion needed to align 100 km 100 km scenes with reference digital elevation models we validated the accuracy by using both known radar corner reflector locations and by the registration of the interferograms with digital maps the topography corrected interferograms are free from all geometric phase terms and they clearly show the geophysical observables of crustal deformation atmospheric phase and ionospheric phase

ice flow of the antarctic ice sheet we present a reference comprehensive high resolution digital mosaic of ice motion in antarctica assembled from multiple satellite interferometric synthetic aperture radar data acquired during the international polar year 2007 to 2009 the data reveal widespread patterned enhanced flow with tributary glaciers reaching hundreds to thousands of kilometers inland over the entire continent this view of ice sheet motion emphasizes the importance of basal slip dominated tributary flow over deformation dominated ice sheet flow redefines our understanding of ice sheet dynamics and has far reaching implications for the reconstruction and prediction of ice sheet evolution

evolution of glacier dammed lakes through space and time brady glacier alaska usa abstract glacier dammed lakes and their associated jokulhlaups cause severe flooding in downstream areas and substantially influence glacier dynamics brady glacier in southeast alaska is well suited for a study of these phenomena because it presently dams 10 large 1 km 2 lakes our objectives are to demonstrate how brady glacier and its lakes have co evolved in the past and to apply this knowledge to predict how the glacier and its lakes will likely evolve in the future to accomplish these objectives we georeferenced a variety of maps airphotos and optical satellite imagery to characterize the evolution of the glacier and lakes we also collected bathymetry data and created bathymetric maps of select lakes despite small advances and retreats the main terminus of brady glacier has changed little since 1880 however it downwasted at rates of 2 3 m y between 1948 and 2000 more than the regional average the most dramatic retreat 2 km and downwasting 120 m have occurred adjacent to glacier dammed lakes and are primarily the result of calving brady glacier is a former tidewater glacier with continued downwasting brady glacier may return to a tidewater regime and enter into a phase of catastrophic retreat the situation at brady glacier is not unique and the lessons learned here can be applied elsewhere to identify future glacier dammed lakes jokulhlaups and glacier instability

ice flow dynamics of the greenland ice sheet from sar interferometry synthetic aperture radar sar interferograms produced from esa s ers 1 satellite provide the first synoptic view of ice flow dynamics of the western sector of the greenland ice sheet glacial motion is detected in the radar ranging direction at millimetric scales across a complete sequence of snow accumulation and melting regimes despite significant variations in their radar scattering properties ice flow evolves from a slow regular motion at the higher elevations at lower elevations motion is strongly convoluted by meter scale undulations in surface topography which have a unique interferometric signature that enables a novel approach for retrieving flow direction inferred flow directions combined with surface displacements in the radar ranging direction yield ice velocity estimates that are within 6 of in situ measurements gathered along a 40 km survey line application of repeat pass sar interferometry to the entire greenland ice sheet should enable precise mapping of its ice flow dynamics at an unprecedented level of spatial detail

accuracy of topographic maps derived from ers 1 interferometric radar an interferometric radar technique for topographic mapping of surfaces promises a high resolution approach to the generation of digital elevation models the authors present analyses of data collected by the synthetic aperture radar instrument on board the ers 1 satellite on successive orbits use of a single satellite in a nearly repeating orbit is attractive for reducing cost and spaceborne hardware complexity also it permits inference of changes in the surface from the correlation properties of the radar echoes the data have been reduced to correlation maps and digital elevation models the correlation maps show that temporal correlation decreases significantly with time but not necessarily at a constant well defined rate likely depending on environmental factors when correlation among passes remains high however it is possible to form digital elevation models analyses of noise expected in ers 1 interferometric data collected over alaska and the southwestern united states indicate that maps with relative errors less than 5 m rms are possible in some regions however orbit uncertainties imply that tie points are required in order to reduce absolute height errors to a similar magnitude the authors find that about 6 tie points per 40 spl times 40 km scene with 5 m rms or better height accuracy are needed to keep systematic map height errors below 5 m rms the performance of the ers 1 radar system for topographic applications though useful for a variety of regional and local discipline studies may be improved with respect to temporal decorrelation errors and absolute height acuity by modifying the orbit repeat period and incorporating precise orbit determination techniques the resulting implementation will meet many but not all objectives of a global mapping mission

observation and analysis of ice flow in the largest greenland ice stream we have applied satellite radar interferometry methods to map the velocity field of the recently discovered northeast greenland ice stream we have used these data in conjunction with ice thickness and accumulation data to determine that the ice stream is in balance within the measurement errors we used control methods to invert a finite element model of ice sheet flow constrained by the data to infer the ice stream s basal shear stress distribution our results reveal that flow in a section of the downstream end has much in common with the streaming flow of the ross ice streams of west antarctica e g a weak bed and fast flow in the presence of low driving stresses for several hundred kilometers along the middle of the ice stream the basal shear stress balances the driving stress in the upstream area where the ice stream is first visible in the velocity data the bed appears to be weak which may contribute to the initiation of the ice stream

atmospheric limitations to repeat track radar interferometry in its recent radar imaging mission the shuttle imaging radar satellite sir c devoted three days to repeat track interferometry we have analyzed the data from a test site in the mojave desert of california although good topography 10 m on 21 m postings was obtained most of the error was caused by turbulent water vapor in the lower atmosphere spatial structure of 6 km and all smaller sizes was observed the rms one way time delay was found to be 0 24 cm essentially identical results were obtained at two wavelengths 24 and 5 7 cm

three dimensional glacial flow and surface elevation measured with radar interferometry outlet glaciers which serve to drain ice from ice sheets seem to be dynamically less stable in north greenland than in south greenland1 2 3 storstrommen a large outlet glacier in northeastern greenland which surged between 1978 and 1984 ref 2 has been well studied in general neither glacier surge mechanisms nor the geographical distribution of the surges are well known conventional satellite radar interferometry can provide large scale topography models with high resolution4 and can measure the radar line of sight component of ice flow vectors5 but cannot map full vector flow fields here we present an interferometry method that combines observations from descending and ascending satellite orbits which assuming ice flow parallel to the topographic surface allows us to use the differing view angles to estimate full three dimensional surface flow patterns the accuracy of our technique is confirmed by the good agreement between our radar based flow model and in situ global positioning system gps reference data at storstrommen radar measurements such as these made regularly and at high spatial density have the potential to substantially enhance our understanding of glacier dynamics and ice sheet flow as well as improve the accuracy of glacier mass balance estimates

the integration of jers 1 and ers sar in differential interferometry for measurement of complex glacier motion complex glacier motion in the grove mountains region antarctica is measured using four pass differential synthetic aperture radar sar interferometry insar the components of the motion vector field are resolved using a 44 day separation japanese earth resources satellite 1 jers 1 insar pair and a european remote sensing satellite 1 2 ers 1 2 tandem insar pair the 44 day temporal baseline provides the sensitivity required to observe the range of ice motion around 8 10 m a 1 and the 1 day short baseline provides the best choice for glacier digital elevation model reconstruction it is remarkable that the scattering field of the jers 1 pair remained coherent over the long time interval and the interferometric fringes are clear the overall ice flow is from east to west downslope and towards lambert glacier the regional flow is obstructed by nunataks extending north south with two wide gaps two narrow glaciers flow past the nunataks and rejoin each other at the downstream end regional morphology and the resolved flow in the grove mountains area suggests that ice flow is channeled throughout this eastern flank of lambert glacier

interferometric estimation of three dimensional ice flow using ascending and descending passes satellite radar interferometry sri provides an important new tool for determining ice flow velocity interferometric measurements made from a single track direction are sensitive only to a single component of the three component velocity vector observations from along three different track directions would allow the full velocity vector to be determined a north south looking synthetic aperture radar sar could provide these observations over large portions of the globe but not over large areas of the polar ice sheets the authors develop and demonstrate a technique that allows the three component velocity vector to be estimated from data acquired along two track directions ascending and descending under a surface parallel flow assumption this technique requires that there are accurate estimates of the surface slope which are also determined interferometrically to demonstrate the technique the authors estimate the three component velocity field for the ryder glacier greenland their results are promising although they do not have yet ground truth data with which to determine the accuracy of their estimates

mesozoic tertiary exhumation history of the altai mountains northern xinjiang china new constraints from apatite fission track data abstract this study uses apatite fission track ft analysis to constrain the exhumation history of bedrock samples collected from the altai mountains in northern xinjiang china samples were collected as transects across the main structures related to palaeozoic crustal accretion events ft results and modeling identify three stages in sample cooling history spanning the mesozoic and tertiary stage one records rapid cooling to the low temperature part of the fission track partial annealing zone circa 70 10 c stage two records a period of relative stability with little if any cooling taking place between 75 and 25 20 ma suggesting the altai region had been reduced to an area of low relief support for this can be found in the adjacent junngar basin that received little if any sediment during this interval final stage cooling took place in the miocene at an accelerated rate bringing the sampled rocks to the earth s surface this last stage linked to the far field effects of the himalayan collision most likely generated the surface uplift and relief that define the present day altai mountains

ice sheet velocity mapping a combined interferometric and speckle tracking approach the first and second radarsat antarctic mapping missions amm 1 and 2 have now acquired interferometric synthetic aperture radar sar overmuch of the ice sheet the radarsat 24 day repeat cycle is nearly ideal for measuring slow ice motion e g 1 but application of sar interferometry is limited in faster moving areas with a 1day repeat period ers 1 2 tandem sar data are much better matched to fast motion but are not always available fortunately several authors have demonstrated the ability to measure velocity in fast moving areas by tracking sar speckle from image to image which works well even in the absence of visible features while these estimates have intrinsically lower resolution and poorer accuracy than direct phase measurements they serve well in areas where there are no data suitable for conventional interferometry this paper describes algorithms i have developed for merging interferometric and speckle tracking data from multiple swaths to form a single seamless mosaic of velocity at each point in the mosaic all the available data are combined to produce estimates of the velocity and the associated error this technique is demonstrated using radarsat data collected over lambert glacier antarctica during amm 1 and 2

partitioning of excess argon between alkali feldspars and glass in a young volcanic system abstract argon system behavior in a young volcanic system has been investigated through analysis of alkali feldspar phenocrysts and glass separates in order to constrain the potential effects of excess ar upon 40ar 39ar ages on feldspars and the crystal melt partition coefficient the importance of syn and post emplacement processes was investigated by sampling glass and pumice clasts derived from varying post eruptive conditions including welded sections vapor phase altered sections and heavily weathered portions of the ignimbrite excess 40ar was detected in all glass samples at 1 10 ppb level yet despite this step heating of co existing alkali feldspar revealed a total gas age of 0 644 0 038 ma compared with the generally accepted age of 0 668 0 004 ma for the timing of emplacement of the arico formation this suggests that despite high concentrations of ar dissolved in the glass at quenching the alkali feldspar crystallizing from excess ar rich melt contain very little excess ar no correlation was observed between excess ar in the glass and total gas age obtained from the co existing alkali feldspar and a partition coefficient kd between feldspar and glass was constrained to be less than 10 3 a correlation between the concentrations of excess ar and atmospheric derived 36ar in the glass manifested as a restricted range of 40ar 36ar up to 320 suggests a uniform mechanism of incorporation of ar from these reservoirs measured in all samples prior to eruption in an assessment of post eruptive factors diffusional loss of ar was not significant during the elevated temperatures up to 670 c associated with welding processes at the relevant timescales over which welding occurs since diffusive exchange occurs between crystals and melt over long periods in the magma chamber ages determined on volcanic alkali feldspars reflect the eruption age as a result of the very low kd but in systems containing excess ar the ultimate accuracy of such ages will be limited by the actual value of kd we observe that biotite a mineral with a higher kd yields anomalously old ages in such systems an observation made previously for other volcanic systems diffusive loss of ar from volcanic glass occurs on timescales approximating that of volcanic eruptions suggesting that excess ar measured in the glass separates is representative of a degassed state and indicates even higher pre eruptive ar concentrations present in the magma chamber prior to eruption thus our measured kd is likely to be a maximum value the significance of physical contamination of feldspar by glass is also highlighted and quantified demonstrating that less than 0 01 by weight of adhered phonolite glass was required to elevate a measured feldspar age this is particularly significant for young samples where any such contamination effect would be magnified proportionally to the small amount of in grown radiogenic ar

changes in the west antarctic ice sheet the portion of the west antarctic ice sheet that flows into the ross sea is thinning in some places and thickening in others these changes are not caused by any current climatic change but by the combination of a delayed response to the end of the last global glacial cycle and an internal instability the near future impact of the ice sheet on global sea level is largely due to processes internal to the movement of the ice sheet and not so much to the threat of a possible greenhouse warming thus the near term future of the ice sheet is already determined however too little of the ice sheet has been surveyed to predict its overall future behavior

interferometric radar observations of glaciar san rafael chile 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 after 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 ma 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 slowermoving 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 m d 1 or 6 4 km a 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

the topsar interferometric radar topographic mapping instrument the authors have augmented the nasa dc 8 airsar instrument with a pair of c band antennas displaced across track to form an interferometer sensitive to topographic variations of the earth s surface during the 1991 dc 8 flight campaign data were acquired over several sites in the us and europe and topographic maps were produced from several of these flight lines analysis of the results indicate that statistical errors are in the 2 4 m range while systematic effects due to aircraft motion are in the 10 20 m range the initial results from development of a second generation processor show that aircraft motion compensation algorithms reduce the systematic variations to 2 m while the statistical errors are reduced to 2 3 m

balance velocities of the greenland ice sheet we present a map of balance velocities for the greenland ice sheet the resolution of the underlying dem which was derived primarily from radar altimetery data yields far greater detail than earlier balance velocity estimates for greenland the velocity contours reveal in striking detail the location of an ice stream in northeastern greenland which was only recently discovered using satellite imagery enhanced flow associated with all of the major outlets is clearly visible although small errors in the source data result in less accurate estimates of the absolute flow speeds nevertheless the balance map is useful for ice sheet modelling mass balance studies and field planning

estimation of ice sheet motion using satellite radar interferometry method and error analysis with application to humboldt glacier greenland satellite radar interferometry provides glaciologists with an important new tool for determining the motion and topography of large ice sheets we examine the sources of error in interferometrically derived ice motion measurements including those errors due to inaccurate estimates of the interferometric baseline several simulations are used to assess baseline accuracy in terms of tie point error and the number and distribution of tie points these results give insight into how best to select tie points and also demonstrate the level of accuracy that can be achieved examination of two representative cases likely to occur in mapping ice sheet motion leads to the conclusion that with adequate tie point information ice velocity can be measured accurately to within a few meters per year a method to correct horizontal velocity estimates for the effect of vertical displacement using surface slopes is also developed finally we estimate the single component velocity field for an area on humboldt glacier northern greenland using interferograms formed from ers i sar images we estimate that these velocity measurements are accurate to within 2 3 m year 1

new precipitation and accumulation maps for greenland annual total precipitation and the annual accumulation on the greenland ice sheet are evaluated and presented in two maps the maps are based on accumulation measurements of 25 1 pits and cores obtained from the upper accumulation zone and precipitation measurements made at 35 meteorological stations in the coastal region to construct the accumulation map the annual precipitation was split into solid and liquid precipitation components annual total precipitation exceeding 2500 mm w e occurs on the southeastern tip of greenland while the minimum precipitation is estimated to occur on the northeastern slope of the ice sheet the mean annual precipitation for all of greenland is 340 mm w e the largest annual accumulation of about 1500 mm w e is found on the glaciers in the southeastern corner of greenland while the smallest accumulation is found on the northeastern slope of the ice sheet west of danmarkshavn the mean accumulation on the greenland ice sheet is estimated at 310 mm w e the regional difference in accumulation is examined with respect to the 850 hpa m bar level circulation the present surface topography is found to play an important role in determining regional accumulation on the ice sheet