radar interferometric mapping of deformation in the year after the landers earthquake although the 1992 landers california earthquake sequence occurred in an area well sampled by geodetic networks1 3 the postseismic deformation in the months following the earthquake has been measured at only 15 geodetic stations4 another shortcoming in the geodetic coverage occurs west of the primary rupture where the existing geodetic observations suggest but cannot resolve sympathetic slip on secondary faults1 such measurements which are needed to place the landers earthquake sequence in the context of a recurring seismic cycle in california can be obtained with the dense spatial coverage provided by satellite radar interferometry5 9 here we present radar maps of the surface deformation field which reveal features that would otherwise have been poorly sampled particularly if the earthquake had occurred in a less accessible area we see triggered slip at the level of several centimetres as far as 100 km from the primary rupture and can resolve the geodetic signal of at least one small magnitude 5 aftershock the amount of surface slip following the main shock is less than a decimetre and is consistent with an exponential decay time of several months for the postseismic deformation

small scale deformations associated with the 1992 landers california earthquake mapped by synthetic aperture radar interferometry phase gradients the landers earthquake mw 7 3 occurred on june 28 1992 and ruptured nearly 100 km of previously mapped and unmapped faults in the mojave desert we use synthetic aperture radar interferometry insar to examine the cumulative surface deformation between april 24 and august 7 1992 in a 100 100 km region surrounding the northern portion of the earthquake rupture also we introduce a technique for manipulating sar interferograms to extract short wavelength displacement information this technique involves computation and subsequent combination of interferometric phase gradient maps the insar results show significant deformation signatures associated with faults fractures dry lake beds and mountainous regions within 75 100 km of the main rupture using the phase gradient method we are able to extract small scale deformation patterns near the main rupture many of the preexisting mapped faults within 50 km of the main rupture experienced triggered slip these include the old woman lenwood johnson valley west calico and calico faults the insar results also indicate right lateral offsets along secondary fractures trending n ne within the left lateral zone of shear between the main rupture and the johnson valley fault additionally there are interesting interferogram fringe signatures surrounding troy dry lake and coyote dry lake that are related to deformation of dry lake beds

monitoring ice sheet behavior from space satellite remote sensing has revolutionized ice sheet research a variety of instruments sensitive to different parts of the electromagnetic spectrum take what the human eye detects as a flat white desert and provide data sets rich in scientific information image based maps of ice sheets are becoming commonplace and have become an integral component of field work more than a pretty picture the digital character of the satellite data from these instruments has become fundamental to the production of elevation motion accumulation and reflectance data sets visible imagery shows the scientist a wealth of features that offer clues to the history and current behavior of the ice sheet radar and microwave imagery provide information from beneath the surface and have been used to estimate snow accumulation rates interferometry principles have recently been applied to measure surface topography and ice motion with unparalleled precision nonimaging instruments also keep a watchful eye monitoring the ice sheet for indications of growth or shrinkage further expansion of the uses of satellite data is anticipated in the future

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

energy density directional spectra of a nearshore wave field measured by interferometric synthetic aperture radar the nearshore wave field within monterey bay california is studied using remote sensing imaging by an interferometric synthetic aperture radar insar and simultaneous ground based measurements it is shown that insar imagery of the ocean surface offers some advantages over conventional synthetic aperture radar because of the direct imaging mechanism of insar quantitative information about the complicated wave field can be obtained the insar derived directional energy density spectra and the corresponding ground based spectra compare well

detection of postseismic fault zone collapse following the landers earthquake stress changes caused by fault movement in an earthquake induce transient aseismic crustal movements in the earthquake source region that continue for months to decades following large events1 4 these motions reflect aseismic adjustments of the fault zone and or bulk deformation of the surroundings in response to applied stresses2 5 7 and supply information regarding the inelastic behaviour of the earth s crust these processes are imperfectly understood because it is difficult to infer what occurs at depth using only surface measurements2 which are in general poorly sampled here we push satellite radar interferometry to near its typical artefact level to obtain a map of the postseismic deformation field in the three years following the 28 june 1992 landers california earthquake from the map we deduce two distinct types of deformation afterslip at depth on the fault that ruptured in the earthquake and shortening normal to the fault zone the latter movement may reflect the closure of dilatant cracks and fluid expulsion from a transiently over pressured fault zone6 8

changes in static stress on southern california faults after the 1992 landers earthquake the magnitude 7 5 landers earthquake of 28 june 1992 was the largest earthquake to strike california in 40 years the slip that occurs in such an earthquake would be expected to induce large changes in the static stress on neighbouring faults these changes in stress should in turn affect the likelihood of future earthquakes stress changes that load faults towards failure have been cited as the cause of small1 5 moderate6 and large7 earthquakes conversely those that relax neighbouring faults have been related to a decrease in seismicity5 here we use an elastic half space model8 to estimate the stress changes produced by the landers earthquake on selected southern california faults including the san andreas we find that the estimated stress changes are consistent with the triggering of four out of the five aftershocks with magnitude greater than 4 5 and that the largest changes 1 10 bar occurring on part of the san bernardino segment of the san andreas fault may have decreased the time to the next magnitude 8 earthquake by about 14 years

post eruptive deformation associated with the 1986 87 and 1989 lava flows of etna detected by radar interferometry we analysed 92 interferograms produced using ers1 sar images taken on etna between may 1992 and october 1993 nineteen show a local subsidence in the eastern flank of the volcano correlated with the location of the 30 october 1986 1 march 1987 and 27 september 9 october 1989 lava flows using fringe unwrapping and data gridding techniques and assuming that over the sampled time window deformation was a linear function of time we derive a map of along range rate of motion the correlation between the deformation field and the area covered by recent lavas suggests that compaction of the lava flows continues several years after the eruptions the area of maximum subsidence 47 mm yr is localised at the narrowing of the 1989 flow between 1500 and 1700 m a s l we observe that subsidence extends outside the lava flows accounting for at least 12 mm yr assuming a relaxation process of the substrate in response to loading produced by recent lavas a simple 1d maxwell visco elastic model predicts a maximum subsidence rate of 25 mm yr and a relaxation time of about 3 5 years the relaxation time agrees with those derived from post eruptive displacements observed by levelling on etna and piton de la fournaise after several eruptions we conclude that at the time of our measurements 25 to 50 of the deformation was related to relaxation of the substrate and the other part due to compaction of the lava flows

sea ice motion from satellite passive microwave imagery assessed with ers sar and buoy motions observing the motion of sea ice from space is analogous to observing wind stress over the wet oceans both provide surface forcing for modeling ocean dynamics ice motion also directly provides the advective component of the equations governing the mass balance of the sea ice cover thus its routine observation from space would be of great value to understanding ice and ocean behavior to demonstrate the feasibility of creating a global multidecadal ice motion record from satellite passive microwave imagery and to quantitatively assess the errors in the estimated ice motions we have tracked ice every 3 days in the arctic ocean and daily in the fram strait and baffin bay during the 8 winter months from october 1992 to may 1993 and daily in the weddell sea during the 8 winter months from march to october 1992 the method which has been well used previously involves finding the spatial offset that maximizes the cross correlation of the brightness temperature fields over 100 km patches in two images separated in time by from 1 to 3 days the resulting ice motions are compared with contemporaneous buoyand sar derived ice motions the uncertainties in the displacement vectors between 5 and 12 km are better than the spatial resolution of the data both 85 ghz data with 12 km spatial resolution and 37 ghz data with 25 km resolution are tracked these trials with the 37 ghz data are new and show quite surprisingly that the error is only about 1 km larger with these data than with the 12 km 85 ghz data errors are typically larger than average in areas of lower ice concentration in the most dynamic regions particularly near the ice edge in the barents and greenland seas and in zones of high shear these passive microwave ice motions show a large increase in spatial detail over motion fields optimally interpolated from buoy and wind observations especially where buoy data are virtually absent such as near coasts and in some passages between the arctic ocean and its peripheral seas the feasibility of obtaining ice motion from the 37 ghz data in addition to the 85 ghz data should allow an important record of ice motion to be established for the duration of the scanning multichannel microwave radiometer smmr special sensor microwave imager ssm i and future microwave sensors that is from 1978 into the next millenium

rapid deformation of the south flank of kilauea volcano hawaii the south flank of kilauea volcano has experienced two large magnitude m 7 2 and m 6 1 earthquakes in the past two decades global positioning system measurements conducted between 1990 and 1993 reveal seaward displacements of kilauea9s central south flank at rates of up to about 10 centimeters per year in contrast the northern side of the volcano and the distal ends of the south flank did not displace significantly the observations can be explained by slip on a low angle fault beneath the south flank combined with dilation deep within kilauea9s rift system both at rates of at least 15 centimeters per year

radar interferometry and its application to changes in the earth s surface geophysical applications of radar interferometry to measure changes in the earth s surface have exploded in the early 1990s this new geodetic technique calculates the interference pattern caused by the difference in phase between two images acquired by a spaceborne synthetic aperture radar at two distinct times the resulting interferogram is a contour map of the change in distance between the ground and the radar instrument these maps provide an unsurpassed spatial sampling density 100 pixels km 2 a competitive precision 1 cm and a useful observation cadence 1 pass month 1 they record movements in the crust perturbations in the atmosphere dielectric modifications in the soil and relief in the topography they are also sensitive to technical effects such as relative variations in the radar s trajectory or variations in its frequency standard we describe how all these phenomena contribute to an interferogram then a practical summary explains the techniques for calculating and manipulating interferograms from various radar instruments including the four satellites currently in orbit ers 1 ers 2 jers 1 and radarsat the next chapter suggests some guidelines for interpreting an interferogram as a geophysical measurement respecting the limits of the technique assessing its uncertainty recognizing artifacts and discriminating different types of signal we then review the geophysical applications published to date most of which study deformation related to earthquakes volcanoes and glaciers using ers 1 data we also show examples of monitoring natural hazards and environmental alterations related to landslides subsidence and agriculture in addition we consider subtler geophysical signals such as postseismic relaxation tidal loading of coastal areas and interseismic strain accumulation we conclude with our perspectives on the future of radar interferometry the objective of the review is for the reader to develop the physical understanding necessary to calculate an interferogram and the geophysical intuition necessary to interpret it

an introduction to the global positioning system and some geological applications receivers equipped to measure dual frequency carrier phase signals from satellites of the global positioning system gps have been capable under special conditions of determining relative horizontal positions among stations separated by one to a few hundred kilometers with a precision of one to several millimeters since the early 1980s the major obstacles to making this capability routine extending it to all parts of the globe and extending it to longer station separations have been equipment cost limitations in the gps satellite constellation arduous data analysis uncertainties in satellite orbits uncertainties in propagation delays associated with variable tropospheric water vapor and difficulties in resolving carrier phase cycle ambiguities recent improvements have occurred in all these areas the increasing ease and reduced cost of gps data acquisition and analysis are having a significant impact on studies of near fault crustal deformation and earthquake processes until recently the province of conventional terrestrial geodetic techniques the enhanced satellite constellation improved models and establishment of global tracking networks have extended several millimeters horizontal positioning capability to station separations of 1000 km or more in virtually all parts of the world this enables study of new classes of tectonic problems that previously were difficult to attack with any geodetic technique examples include a complete kinematic description of ongoing crustal deformation in broad complex continental plate boundary zones and measurement of relative plate motion at convergent boundaries where global models may be poorly constrained

crustal deformation from 1992 to 1995 at the mid atlantic ridge southwest iceland mapped by satellite radar interferometry satellite radar interferometry observations of the reykjanes peninsula oblique rift in southwest iceland show that the reykjanes central volcano subsided at an average rate of up to 13 millimeters per year from 1992 to 1995 in response to use of its geothermal field interferograms spanning up to 3 12 years also include signatures of plate spreading and indicate that the plate boundary is locked at a depth of about 5 kilometers below that depth the plate movements are accommodated by continuous ductile deformation which is not fully balanced by inflow of magma from depth causing subsidence of the plate boundary of about 6 5 millimeters per year

estimates of currents in the nearshore ocean region using interferometric synthetic aperture radar microwave remote sensing of the ocean surface using synthetic aperture radar sar is of great interest due to its high resolution and the potential of nearly instantaneous coverage of large areas the experience accumulated in ocean imaging by sar indicates that there are substantial difficulties in interpretation of the obtained images of ocean wave systems ocean surface currents are almost impossible to obtain from sar imagery the present paper reports on measurements of ocean surface currents based on imaging of the nearshore regions of the monterey bay using an interferometric sar the method is a modification of the conventional sar which employs two spatially separated antennas this modification allows one to construct an interferogram from the two complex images obtained by each antenna which provides direct mapping of the observed surface velocities the present results demonstrate the advantages of this novel technique in studying nearshore ocean surface dynamics in particular spatial distribution of the longshore current velocity is presented

the 1995 grevena northern greece earthquake fault model constrained with tectonic observations and sar interferometry after the 1995 grevena ms 6 6 event in northern greece we mapped the earthquake fault break in detail the surface break is small 8 12 km long 4 cm slip compared to the moment release of the event however the morphologic and tectonic study of the active faults in the field and using the spot satellite imagery suggests that the earthquake ruptured part of a much larger fault system including interconnecting segments we used sar interferometry of the satellite ers 1 imagery to characterize the coseismic displacement field this shows a kidney shaped zone of subsidence reaching 30 cm flanked by an uplift zone reaching 5 cm we reproduce this field using dislocations in an elastic half space and our observations of the fault system this requires 1 m slip from 4 to 15 km depth on a main normal fault segment dipping nnw our preliminary model includes significant ne dipping scissors faulting at the eastern end of the rupture clearly seen in the interferograms

sensing the ups and downs of las vegas insar reveals structural control of land subsidence and aquifer system deformation land subsidence in las vegas nevada united states between april 1992 and december 1997 was measured using spaceborne interferometric synthetic aperture radar the detailed deformation maps clearly show that the spatial extent of subsidence is controlled by geologic structures faults and sediment composition clay thickness the maximum detected subsidence during the 5 75 yr period is 19 cm comparison with leveling data indicates that the subsidence rates declined during the past decade as a result of rising ground water levels brought about by a net reduction in ground water extraction temporal analysis also detects seasonal subsidence and uplift patterns which provide information about the elastic and inelastic properties of the aquifer system and their spatial variability

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

spaceborne and airborne imaging radar observations of sand dunes seasat and aircraft radar imagery of five areas of sand dunes in the southwestern united states and northwestern mexico have been studied and compared to landsat imagery and air photos radar imaging system parameters included two wavelengths 23 5 cm and 3 0 cm and incidence angles ranging from 0 through 70 horizontal polarization was used on both transmit and receive for all radar imagery resolution ranged from 12 m for the best aircraft imagery to about 40 m for seasat sand dunes studied behave as smooth surfaces for both 23 5 cm and 3 0 cm wavelength radars and the reflection of the radar energy from the dunes is specular measured surface irregularities of sand dunes are inadequate to cause either bragg or incoherent scattering of the incident radar beam in the 3 0 to 23 5 cm wavelength region this results in dunes that appear dark when imaged at incidence angles greater than the angle of repose because none of the dune faces are oriented properly normal to the incident radar beam to give a specular return to the antenna a scattering model based on specular return from randomly oriented flat dune faces was used and experimentally verified using aircraft data imaging radar resolution must be adequate to define the changing distribution of dune faces that characterize the dune morphology linear dune features with interdune spacing of less than 3 4 resolution elements are not recognizable distributed dune types such as barchans or stars must be more than about 10 resolution elements in each direction to be defined on the imagery illumination direction of the radar beam is important directional dune features must be oriented within about 60 of perpendicular to the radar illumination direction in order to be imaged if such features have considerable topographic relief the direction of radar illumination is less important because directional trends in unexplored regions cannot be predicted availability of radar imagery from two directions greatly facilitates interpretation of dune morphology and hence conclusions about causative wind regimes

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

mapping small elevation changes over large areas differential radar interferometry a technique that uses synthetic aperture radar sar images to measure very small 1 cm or less surface motions with good resolution 10 m over large swaths 50 km is presented along with experimental results the method could be used for accurate measurements of many geophysical phenomena including swelling and buckling in fault zones residual displacements from seismic events and prevolcanic swelling the method is based on sar interferometry where two images are made of a scene by simultaneously flying two physically separated antennas then the phases of corresponding pixels are differenced and altitude formation is deduced from some simple computation and image rectification it is also possible to use one antenna flown twice over the same scene then if the second flight exactly duplicates the track of the first an interesting possibility occurs there would be no phase changes between the images at all unless there was a physical change in the scene such as ground swelling that would alter the distance from some resolution element to the antenna since the phase changes all occur at the short carrier wavelength the basic limitation on sensitivity is only the phase noise in the system when the two imaging passes are made from flight tracks that are separated which is the case with the seasat images used here it is no longer possible to distinguish surface changes from the parallax caused by topography however with some additional computation a third image made at some other baseline may be used to remove the topography and leave only the surface changes this method was applied using seasat data to an imaging site in imperial valley california where motion effects were observed that were ascribed to the expansion of water absorbing clays phase change images of this area are shown along with associated ground truth about the presence of water problems with the technique are explored along with a discussion of future experimental possibilities on upcoming sar missions like earth observing system eos earth resources satellite ers 1 sir c and the venus imaging radar magellan

effects of subglacial geothermal activity observed by satellite radar interferometry we use one day synthetic aperture radar sar interferograms from data of the earth remote sensing satellites ers 1 and ers 2 to study ice flow and uplift of two surface depressions within the vatnajokull ice cap iceland the ice cauldrons are created by melting at sub glacial geothermal areas meltwater accumulates in a reservoir under the cauldrons over 2 to 3 years until it drains in a jokulhlaup under the ice dam surrounding the reservoir the ice surface in the depressions drops down by several tens of meters during these draining events but rises again as ice flows into the depressions until a jokulhlaup occurs again using sar interferograms we quantify an uplift rate of about 2 to 18 cm day within the jokulhlaup cycle varying with the surface slope of the depressions the uplift rate is high during the first months after a jokulhlaup when the cauldron is relatively deep with steep slopes but the uplift rate decreases as the cauldron is gradually filled a simple axisymmetric model simulating the ice flow into one of the depressions describes quantitatively the filling rate of the cauldron and qualitatively the shape of the ice flow field the best fit model has an ice flow law parameter a0 that is about one order of magnitude lower than typically estimated for temperate glaciers

water rock interaction during meteoric flushing of a limestone implications for porosity development in karstified petroleum reservoirs abstract the lincolnshire limestone comprising a succession of jurassic wackestones packstones and oolitic grainstones forms an important carbonate aquifer in eastern england meteoric waters enter at outcrop and penetrate between confining strata at least 25 km down dip this water dissolves and interacts with the limestone and even water samples collected at or near outcrop are calcite saturated net limestone dissolution is thus a process that is most dominant in the near surface environment water samples taken at increasing distances from outcrop have increasing sr and mg contents and 13c values of dissolved bicarbonate increase from 15 to 8 pdb while 87sr 86 falls from 0 7082 to 0 7077 these data are interpreted to be the result of a coupled carbonate dissolution reprecipitation process modeling of the c and sr isotopic data indicates that primary jurassic carbonate ooids micrite fossils is dissolving preferentially over burial cements isotopically light bulk rock carbon near joint surfaces suggests that reprecipitation of calcite in the form of cement could be concentrated preferentially in and near joints the lincolnshire limestone may be used as an analogue for karstified petroleum reservoirs specifically those which have been buried and lost their unstable carbonate minerals aragonite high mg calcite prior to uplift and karstification the present water ch mical data suggest that in such reservoirs influx of meteoric water at an unconformity creates porosity and enhances permeability through limestone dissolution but this may be concentrated close tens of meters to the unconformity this dissolution also promotes surface erosion which limits the thickness of preserved karstified limestone the dissolution reprecipitation process that occurs at greater distances from outcrop has effects that are far more widespread lens of kilometers leading to an increase in microporosity e g micrite leaching and a possible reduction in joint or fracture porosity such factors should be taken into account when exploring for and appraising karstified petroleum reservoirs

a comparison between short term co seismic and long term one year slip for the landers earthquake measurements from strong motion and sar interferometry the long term relative displacements in the region of the landers earthquake were measured over a time interval of about one year using sar interferometry we have also evaluated the co seismic slip distribution through a strong motion inversion cotton and campillo 1995 and these results are used to compute the static displacement field associated with the co seismic slip we show that the synthetic static and the interferometry measurements are in very good agreement to give an evaluation of the magnitude of the pre and post seismic slip a residual interferogram between data and synthetic co seismic slip has been computed a simple analysis of the residual displacements gives a limiting value of the amplitude of pre and post seismic slip integrated over a period of one year and suggests that it is on the order of a few decimeters

aseismic creep along the san andreas fault northwest of parkfield ca measured by radar interferometry ers 1 radar images acquired 14 months apart studied by differential radar interferometry show the wide area distribution of aseismic creep along the fault segment northwest of parkfield california a sharp discontinuity in the interferometric phase of less than 2 cm equivalent line of sight displacement extends over 80 km in the differential interferogram coincident with the mapped trace of the active san andreas fault and consistent with the expected and measured fault motion although patterns of strain associated with the transition from locked to creeping are not clearly identifiable a decrease in creep displacement from northwest to southeast along the fault is visible the observations are in agreement with a model of elastic deformation constrained by in situ data that supports a maximum expected deformation signature of 10 mm across the image

surface displacement of the 17 may 1993 eureka valley california earthquake observed by sar interferometry satellite synthetic aperture radar sar interferometry shows that the magnitude 6 1 eureka valley earthquake of 17 may 1993 produced an elongated subsidence basin oriented north northwest parallel to the trend defined by the aftershock distribution whereas the source mechanism of the earthquake implies a north northeast striking normal fault the 3 millimeter accuracy of the radar observed displacement map over short spatial scales allowed identification of the main surface rupture associated with the event these observations suggest that the rupture began at depth and propagated diagonally upward and southward on a west dipping north northeast fault plane reactivating the largest escarpment in the saline range

dynamic deformation of etna volcano observed by satellite radar interferometry satellite radar interferometry of mt etna volcano sicily italy reveals a sequence of deformation characterized by deflation during the end of the 1993 eruption inflation from 1993 1995 with an increase in the inflation rate immediately before its resumed eruptive activity in late 1995 this was followed by very low deformation levels during the following year the source of the deformation changed from a depth of 9 km during deflation to more than 11 14 km during the subsequent inflation consistent with a model in which deflation at shallower levels is followed by inflation at greater depth as the volcano system recharges from below before its next eruption this study demonstrates that radar interferometry provides an important contribution towards understanding the dynamic deformation of volcanoes by revealing large scale changes in their pre eruption deformation rates radar interferometry could play an important role in volcano eruption monitoring

on the derivation of coseismic displacement fields using differential radar interferometry the landers earthquake we present a map of the coseismic displacement field resulting from the landers california june 28 1992 earthquake derived using data acquired from an orbiting high resolution radar system we achieve results more accurate than previous space studies and similar in accuracy to those obtained by conventional field survey techniques data from the ers 1 synthetic aperture radar instrument acquired in april july and august 1992 are used to generate a high resolution wide area map of the displacements the data represent the motion in the direction of the radar line of sight to centimeter level precision of each 30 m resolution element in a 113 km by 90 km image our coseismic displacement contour map gives a lobed pattern consistent with theoretical models of the displacement field from the earthquake fine structure observed as displacement tiling in regions several kilometers from the fault appears to be the result of local surface fracturing comparison of these data with global positioning system and electronic distance measurement survey data yield a correlation of 0 96 thus the radar measurements are a means to extend the point measurements acquired by traditional techniques to an area map format the technique we use is 1 more automatic 2 more precise and 3 better validated than previous similar applications of differential radar interferometry since we require only remotely sensed satellite data with no additional requirements for ancillary information the technique is well suited for global seismic monitoring and analysis

synthetic aperture radar interferogram of the 1995 kobe earthquake and its geodetic inversion crustal deformation from the mw 6 9 january 17 jst 1995 kobe japan earthquake was detected by synthetic aperture radar interferometry sar more than four fringes with ellipsoidal shape were observed surrounding the nojima fault on northern awaji island suggesting a shortening in the range of more than 11 75 4 cm on the east side of this fault the sar sensor shows ellipsoidal fringes with an uplift zone in the rokko mountain range and sharp subsidence around the kobe area in addition to sar interferometry leveling and gps surveys were conducted by the geographical survey institute of japan gsi the permanent gps observation sites of the gsi monitored ground displacements from this earthquake with daily updates based on these geodetic data we estimated the faulting mechanism of the earthquake the analytical procedure was an iterative linearized least squares method a multi fault system consisting of six subfaults was constructed the locations of the estimated model faults were close to those of the active faults in this area running from the nojima fault on awaji island to the rokko faults in the kobe area with a change in strike angle near the ohtsuki fault in the kobe area the slip vectors of the nojima fault on awaji island had a large component of reverse dip slip of about 2 m and strike slip of about 2 m in the kobe area the amount of slip ranged from 0 3 to 1 2 m and was small compared with the slip on the nojima fault the results of these slip vectors agreed well with the results of other seismic geodetic and field studies

shallow and peripheral volcanic sources of inflation revealed by modeling two color geodimeter and leveling data from long valley caldera california 1988 1992 we refined the model for inflation of the long valley caldera near mammoth lakes california by combining both geodetic measurements of baseline length and elevation changes baseline length changes measured using a two color geodimeter with submillimeter precision revealed that the resurgent dome started to reinflate in late 1989 measurements between late 1989 and mid 1992 revealed nearly 13 cm of extension across the resurgent dome geodetic leveling surveys with approximately 2 mm precision made in late 1988 and in mid 1992 revealed a maximum of about 8 cm of uplift of the resurgent dome two ellipsoidal sources satisfy both the leveling and two color measurements whereas spherical point sources could not the model s primary inflation source is located 5 5 km beneath the resurgent dome with the two horizontal axes being nearly equal in size and the vertical axis being 4 times the length of the horizontal axes a second ellipsoidal source was added to improve the fit to the two color measurements this secondary source is located at a depth between 10 and 20 km beneath the south moat of the caldera and has the geometry of an elongated ellipsoid or pipe that dips down to the northeast in addition the leveling data suggest dike intrusion beneath mammoth mountain during the 1988 1992 interval which is likely associated with an intense swarm of small earthquakes during the summer of 1989 at that location our analysis shows the dike intrusion to be the shallowest of the three sources with a depth range of 1 3 km below the surface to the top of the intrusion

two examples of the use of sar interferometry on displacement fields of small spatial extent interferometric combination of pairs of sar images acquired by the european ers 1 satellite maps deformation fields associated with two phenomena both of small spatial extension and located in se france the one is rapid terrain deformation caused by a landslide near the city of saint etienne de tinee and the other is slower subsidence caused by underground coal mining near gardanne unlike interferometric measurement of wide field deformation atmospheric propagation heterogeneity is not an accuracy limiting factor although the radar data confirm prior knowledge concerning the landslide such an application of sar interferometry appears difficult under normal conditions of observation using current spaceborne radar systems the study of soil subsidence however can be generalized and improves prior knowledge of the displacement field which has here been modeled assuming elastic deformation in a half space from several sources the two examples help to understand the limits of the interferometric technique

spatial and temporal distribution of slip for the 1992 landers california earthquake we have determined a source rupture model for the 1992 landers earthquake m w 7 2 compatible with multiple data sets spanning a frequency range from zero to 0 5 hz geodetic survey displacements near field and regional strong motions broadband teleseismic waveforms and surface offset measurements have been used explicitly to constrain both the spatial and temporal slip variations along the model fault surface our fault parameterization involves a variable slip multiple segment finite fault model which treats the diverse data sets in a self consistent manner allowing them to be inverted both independently and in unison the high quality data available for the landers earthquake provide an unprecedented opportunity for direct comparison of rupture models determined from independent data sets that sample both a wide frequency range and a diverse spatial station orientation with respect to the earthquake slip and radiation pattern in all models consistent features include the following 1 similar overall dislocation patterns and amplitudes with seismic moments of 7 to 8 10 26 dyne cm seismic potency of 2 3 to 2 7 km 3 2 very heterogeneous unilateral strike slip distributed over a fault length of 65 km and over a width of at least 15 km though slip is limited to shallower regions in some areas 3 a total rupture duration of 24 sec and an average rupture velocity of 2 7 km sec and 4 substantial variations of slip with depth relative to measured surface offsets the extended rupture length and duration of the landers earthquake also allowed imaging of the propagating rupture front with better resolution than for those of prior shorter duration strike slip events our imaging allows visualization of the rupture evolution including local differences in slip durations and variations in rupture velocity rupture velocity decreases markedly at shallow depths as well as near regions of slip transfer from one fault segment to the next as rupture propagates northwestward along the multiply segmented fault length the rupture front slows as it reaches the northern limit of the johnson valley landers faults where slip is transferred to the southern homestead valley fault an abrupt acceleration is apparent following the transfer this process is repeated and is more pronounced as slip is again passed from the northern homestead valley fault to the emerson fault although the largest surface offsets were observed at the northern end of the rupture our modeling indicates that substantial rupture was also relatively shallow less than 10 km in this region

tropospheric corrections of sar interferograms with strong topography application to etna the accuracy of spaceborne geodetic techniques including sar interferometry is limited by the time and spatial variation and altitude dependance of the propagation delay of electomagnetic waves in the lower troposphere particularly in mountainous areas in this paper we use a 1d model developed for tropospheric corrections of gps and doris measurements to correct sar data the differential tropospheric delay is computed at each pixel of the interferogram from ground temperature humidity and pressure using two empirical parameters calibrated from several radio soundings acquired in various latitude and climate conditions it is shown that with such a model given the 3300 meters topography of etna tropospheric variations can generate up to 4 phase rotations between the top and the bottom of the volcano in 16 out of the 20 interferograms processed with images acquired between august 1992 and october 1993 correction of the tropospheric effect reduces the number of fringes associated with the 1991 93 eruption from previous estimates the remaining deformation is consistant with a deforming source located at a depth of 14 1 km during the second half of the eruption the subsidence rate at the top of the volcano is roughly stable at 13 3 mm month these values are in good agreement with tiltmeter data collected on etna during the same period and with the estimated volume of erupted material no significant deformation can be observed during the last month of eruption inflation of the volcano seems to resume immediately after the end of the eruption at a rate of 3 mm month

co seismic displacements of the 1992 landers earthquake sequence we present co seismic displacement vectors derived from global positioning system gps measurements of 92 stations in southern california these gps results are combined with five well determined gps displacement vectors from continuously tracking stations of the permanent gps geodetic array as well as line length changes from usgs geodolite and two color laser trilateration observations to determine a self consistent set of geodetic data for the earthquake these combined displacements are modeled by an elastic dislocation representation of the primary fault rupture planes on average the model residuals are about twice the estimated measurement errors

postseismic deformation following the 1989 m 7 1 loma prieta california earthquake postseismic deformation along a 90 km profile bisecting the projected surface trace of the coseismic rupture of the 1989 loma prieta earthquake has been monitored by frequent gps surveys for 3 3 years following the earthquake in addition to the expected deformation associated with secular strain accumulation on the san andreas and calaveras faults deformation associated with postseismic readjustment has been detected most of that deformation can be attributed to 1 5 m right lateral and 0 9 m reverse postseismic slip on a 5 km wide downdip extension of the loma prieta rupture in addition there seems to be a 0 1 m postseismic collapse of the loma prieta rupture zone in the direction perpendicular to the plane of the rupture the fault normal n48 e surface displacements plotted as a function of time exhibit a curvature suggesting a relaxation time of about 1 4 years similar plots of the fault parallel n42 w displacement components do not exhibit significant curvature presumably the deformation shown in those plots is dominated by secular strain accumulation along the san andreas and calaveras faults rather than postseismic relaxation

remote sensing of ocean currents a method of remotely measuring near surface ocean currents with a synthetic aperture radar sar is described the apparatus consists of a single sar transmitter and two receiving antennas the phase difference between sar image scenes obtained from the antennas forms an interferogram that is directly proportional to the surface current the first field test of this technique against conventional measurements gives estimates of mean currents accurate to order 20 percent that is root mean square errors of 5 to 10 centimeters per second in mean flows of 27 to 56 centimeters per second if the full potential of the method could be realized with spacecraft then it might be possible to routinely monitor the surface currents of the world s oceans

observation and modelling of the saint etienne de tinee landslide using sar interferometry abstract six different interferograms of the la clapiere landslide were derived from ers 1 sar images during the period aug 20 sept 4 1991 the coherence of the associated images is shown to remain significant over most of the surface of the landslide during the two weeks of the survey the interferograms are remarkably similar and indicate steady state displacements over at least 12 days the displacement field derived from the interferograms is shown to be characterized by a non uniform displacement gradient from top to bottom and reveals a significantly faster motion of the western part of the landslide the amplitude of the motion deduced from interferometry is shown to be in good agreement with ground measurements finally sar interferometry is shown to be able to evidence small scale instabilities which may not be observed with discrete ground measurements

changes in state of stress on the southern san andreas fault resulting from the california earthquake sequence of april to june 1992 the april to june 1992 landers earthquake sequence in southern california modified the state of stress along nearby segments of the san andreas fault causing a 50 kilometer segment of the fault to move significantly closer to failure where it passes through a compressional bend near san gorgonio pass the decrease in compressive normal stress may also have reduced fluid pressures along that fault segment as pressures are reequilibrated by diffusion that fault segment should move closer to failure with time that fault segment and another to the southeast probably have not ruptured in a great earthquake in about 300 years

crustal deformation measurements using repeat pass jers 1 synthetic aperture radar interferometry near the izu peninsula japan we have examined the precision of interferometric synthetic aperture radar sar measurements of surface deformation of the earth using 24 cm wavelength data acquired by the japanese earth resources satellite 1 jers 1 spacecraft over the izu peninsula japan radar images acquired from 1993 to 1994 were made into 11 independent interferograms apparent water vapor signatures with equivalent path delay up to 16 cm preclude reliable estimates of small deformation from only one interferogram averaging of interferograms improves the measurement such that residuals between precise leveling data and corresponding radar observations are less than 1 cm rms there is a clear correlation between cloud cover extent on a given day and the spectral power of the water vapor signature in an interferogram formed from the corresponding day s sar image also several water vapor features associated with local topography are consistent with accumulation of moisture on the windward side of an on shore flow analysis shows that decorrelation in mixed conifer and deciduous forest is a weak function of the time between observations imposing a roughly constant level of additional decorrelation relative to scatter stable urban areas furthermore in high relief areas correlation is high when the baseline is tens of meters but drops rapidly with increasing baseline length with strict control of the interferometric baseline length adequate correlation in interferograms over very long time periods is possible even in forested high relief regions observed inflation over ito city is consistent with a localized magma source at depth suggesting a different mechanism from the previous episode in 1989

what drives the north atlantic oscillation s temperature anomaly pattern part i the growth and decay of the surface air temperature anomalies abstractcomposite analysis is used to examine the physical processes that drive the growth and decay of the surface air temperature anomaly pattern associated with the north atlantic oscillation n