combination of precise leveling and insar data to constrain source parameters of the mw 6 5 26 december 2003 bam earthquake we used new precise leveling data acquired 40 days after the bam earthquake in combination with radar interferometry observations from both ascending and descending orbits to investigate static deformation associated with the 2003 bam earthquake we invert this geodetic data set to gain insight into the fault geometry and slip distribution of the rupture the best fitting dislocation model is a steeply east dipping right lateral strike slip fault that has a size of 11 by 8 km and strikes n2 w we find that such smooth geometry fits available geodetic data better than previously proposed models for this earthquake our distributed slip model indicates a maximum strike slip of 3 m occurring about 3 to 5 km deep the slip magnitude and depth of faulting taper to the north where the fault approaches the bam city inclusion of crustal layering increases the amount of maximum slip inferred at depth by about 4

monitoring landslide displacements by using ground based synthetic aperture radar interferometry application to the ruinon landslide in the italian alps 1 synthetic aperture radar sar data are collected by a ground based radar system forming the synthetic aperture by the sliding of the antennas on a linear rail coherent sar processing converts the raw data into a complex image the phase of each image pixel contains information on the target sensor distance and can be exploited as a ranging tool the interferometric technique based on the comparison between paired and coherent sar images taken at different times permits the quantitative extraction of this information thus allowing the monitoring of the morphological changes the portable device used in this application was developed by the joint research center ispra italy specifically for measurements in the field it is known as linear sar and it is able to provide 17 ghz measurements with a 2 8 m synthetic aperture a measurement campaign lasting about 1 week was performed between july and august 2000 for monitoring superficial displacements at the ruinon landslide a 30 million m3 rockslide in the italian alps two sequences of interferograms are presented and discussed the interpretation of the sequences has allowed us to derive multitemporal deformation maps of the test area thus showing the entire displacement field of those landslide sectors characterized by higher radar reflectivity and coherence displacement rates up to 1 2 mm h 1 have been measured with a pixel resolution of 5 m and a measurement precision of 0 75 mm the results have been validated by using ground truth data obtained through automatic extensometers and topographic measurements discrepancies are limited to a few millimeters

geohazards monitoring in roma from insar and in situ data outcomes of the pangeo project within the pangeo project financed by the european commission under the 7th framework program the geological survey of italy ispra and the urban planning department of the city of roma developed a geodatabase and map of the geological hazards for the territory of roma integrating remotely sensed data psinsar permanent scatterer interferometry synthetic aperture radar and in situ geological information numerous thematic layers maps and inventories of hazards e g landslides sinkholes cavities geological and hydrogeological data added to historical and recent urbanization information were compared to the permanent scatterer ps data from the european remote sensing satellites ers 1 2 1992 2000 and envisat 2002 2005 descending scenes in order to produce a ground stability layer gsl based on the ps data most of the territory appears stable almost 70 of ps velocities are within 1 mm year about 14 of the pss show positive line of sight los velocities measured along the los of the satellite between 1 and 3 mm year and more than 2 exceed 3 mm year more than 11 of pss show negative los velocities between 1 and 3 mm year while about 3 exceed 3 mm year with tens of the pss showing velocities over 20 mm year the gsl is comprised of polygons or multi polygons multipart polygons grouping individual polygons under a single identifier geohazard enclosing areas where geohazards have been pointed out by ps data and or in situ surveys observed instabilities and by polygons enclosing areas potentially affected by geohazards based on the available knowledge of the territory potential instabilities in roma s gsl 18 multi polygons covering ca 600 km2 related to observed instabilities have been outlined where ground movements could be detected through insar data or where landslides and sinkholes are known to have occurred other 13 multi polygons covering nearly 900 km2 concern areas where the potential occurrence of geohazards was inferred by combining geological and or geothematic data potential instabilities the geohazards mapped in roma have been landslides collapsible grounds compressible grounds groundwater abstraction mining man made ground tectonic movements and volcanic inflation deflation the lattermost is the likely cause of the significant uplift observed in the alban hills area however this paper focuses on two more currently impending hazards subsidence and sinkholes in general sinkhole prone areas areas of dense underground cavities are hard to discern from satellite data but can be revealed by ruling out other potential causes of observed ground movement based on in situ data subsiding zones are effectively detected by the available psinsar dataset over a total extent of about 60 km2 mostly overlapping the recent alluvial areas of the tiber and its tributaries pss show a very different behaviour inside and outside the historical centre inside loading by anthropogenic construction and man made ground since ancient times has led to an almost complete consolidation of the recent river deposits marked by modest to absent subsidence in contrast outside subsidence clearly stands out with negative los velocities that although generally within several mm year can locally exceed 25 mm year ps data have provided motion information at both a regional and local scale up to the scale of a single building closer to the sea in the tiber delta area velocities increase especially above recently reclaimed marsh areas rich in peat and organic clays velocities can change significantly over short distances as in the international airport area reflecting the local stratigraphic setting the same occurs in the two subsiding areas located within the alban hills volcanic complex which is generally affected by ground uplift as a whole psinsar ground motion velocities offer a significant contribution to susceptibility and hazard recognition studies in particular such a method provides a fast and effective tool available to local authorities to monitor ground and building behaviour possibly allowing for timely prevention activities especially when coupled with appropriate in situ knowledge

geodetic imaging reservoir monitoring using satellite interferometry summary fluid fluxes within subsurface reservoirs give rise to surface displacements particularly over periods of a year or more observations of such deformation provide a powerful tool for mapping fluid migration within the earth providing new insights into reservoir dynamics in this paper we use interferometric synthetic aperture radar insar range changes to infer subsurface fluid volume strain at the coso geothermal field furthermore we conduct a complete model assessment using an iterative approach to compute model parameter resolution and covariance matrices the method is a generalization of a lanczos based technique which allows us to include fairly general regularization such as roughness penalties we find that we can resolve quite detailed lateral variations in volume strain both within the reservoir depth range 0 4 2 5 km and below the geothermal production zone 2 5 5 0 km the fractional volume change in all three layers of the model exceeds the estimated model parameter uncertainty by a factor of two or more in the reservoir depth interval 0 4 2 5 km the predominant volume change is associated with northerly and westerly oriented faults and their intersections however below the geothermal production zone proper the depth range 2 5 5 0 km there is the suggestion that both north and northeast trending faults may act as conduits for fluid flow

coupled magma chamber inflation and sector collapse slip observed with synthetic aperture radar interferometry on mt etna volcano 1 volcanoes deform dynamically due to changes in both their magmatic system and instability of their edifice mt etna features vigorous and almost continuous eruptive activity from its summit craters and periodic flank eruptions even though its shape is that of a large stratovolcano its structure features two rift systems and a flank collapse structure similar to hawaiian shield volcanoes we analyze european remote sensing ers satellite differential interferometric synthetic aperture radar insar data 1993 1996 for mt etna spanning its quiescence from 1993 through the initiation of renewed eruptive activity in late 1995 we use synthetic aperture radar sar data from both ascending and descending ers satellite tracks comparison of independent interferograms covering the first 2 years of the inflationary period shows a pattern consistent with inflation of the volcano calculation of the tropospheric path delay based on meteorological data does not change this interpretation interferograms from late summer 1995 1996 show no significant deformation joint inversion of interferograms from ascending and descending satellite tracks require both inflation from a spheroidal magmatic source located beneath the summit at 5 km below sea level and displacement of the east flank of etna along a basal decollement both sources of deformation were contemporaneous within the resolution of our data and suggest that inflation of the central magma chamber acted to trigger slip of etna s eastern flank these results demonstrate that flank instability and recharge of a volcano s magma system must both be considered toward understanding how volcanoes work and in their hazard evaluation

estimation of seismic and aseismic deformation in mexicali valley baja california mexico in the 2006 2009 period using precise leveling dinsar geotechnical instruments data and modeling ground deformation and seismicity in mexicali valley baja california mexico the southern part of the mexicali imperial valley are influenced by active tectonics and human activity in this study data from two successive leveling surveys in 2006 and 2009 2010 are used to estimate the total deformation occurred in mexicali valley during 2006 2009 the leveling data span more than 3 5 years and include deformation from several natural and anthropogenic sources that acted at different temporal and spatial scales during the analyzed period because of its large magnitude the aseismic anthropogenic deformation caused by fluid extraction in the cerro prieto geothermal field obscures the deformation caused by other mechanisms and sources the method of differential interferograms stacking was used to estimate the aseismic interseismic tectonic and anthropogenic components of the observed displacement using sar images taken in 2007 during a period when no significant seismicity occurred in the study area after removing the estimated aseismic signal from the leveling data residual vertical displacement remained and to identify possible sources and mechanisms of this displacement a detailed analysis of records from tiltmeters and creepmeters was performed the results of this analysis suggest that the residual displacement is mainly caused by moderate sized seismicity in the area of study modeling of the vertical ground deformation caused by the coseismic slip on source fault primary mechanism of the two most important earthquakes may 24 2006 mw 5 4 and december 30 2009 mw 5 8 was performed the modeling results together with the analysis of geotechnical instruments data suggests that this moderate sized seismicity influences the deformation in the study area by coseismic slip on the source fault triggered slip on secondary faults and soft sediments deformation

pine island glacier antarctica a study of the catchment using interferometric synthetic aperture radar measurements and radar altimetry we have mapped the catchment region of pine island glacier antarctica using ers 1 and ers 2 radar altimetry and synthetic aperture radar sar interferometry the radar altimeter data were converted to a digital elevation model dem using an algorithm we developed that corrects for slopes and grids the surface elevation estimates in a single step we used the dem and precise postprocessed orbital ephemerides to generate terrain corrected and geolocated complex radar imagery interferometric pairs on both ascending and descending orbits corrected by means of the dem yielded components of ice surface motion in the radar look direction where possible we combined velocity components from ascending and descending tracks to produce vector velocities elsewhere the velocity was calculated by assuming the direction of motion is down the surface slope the velocity map reveals a system of tributaries that channel ice from the basin like catchment area into the fast flowing downstream trunk of pine island glacier none of the tributaries shows a clear onset region instead the ice speed gradually increases along with a gradual decrease in driving stress along each tributary we combined the velocity data with sparse information on ice thicknesses and snow accumulation rates to calculate an approximate mass balance we find that the region is in balance within an uncertainty of 30 of the total mass throughput

insar assessment of surface deformations in urban coastal terrains associated with groundwater dynamics monitoring ground deformations arising from groundwater dynamics in dense urban coastal terrains is crucial for the sustainable development of infrastructures in these highly populated areas the city of montreal which is located in the saint laurent plain in eastern canada with its fast growing populations is a unique case study for other similar cities in coastal terrains the city undergoes high level house foundation damages with densities reaching up to 89 repairs km2 resulting from time dependent ground deformations that are correlated to groundwater dynamics and evapotranspiration using radarsat 2 c band synthetic aperture radar interferometry we observe 3 to 5 mm ground line of sight displacement variations temporally outphased by few months relative to the 2 m subartesian aquifer hydraulic head variations the deformations are observed over a 60 km2 area located in the central part of the montreal island in canada from 2008 to 2010 we observe displacements of 1 mm year uplift in the areas covered by 15 m thick clay layer these displacements are well correlated to the number of house repairs we also observe 2 mm year subsidence on elevated terrains associated with evapotranspiration the amplitudes of the displacements observed during this two year study are significant when integrated over the average lifetime of urban structures we conclude that the observed ground deformations are related to the seasonal variation of hydraulic head in most of the areas of montreal moreover wetter climate forecasts over upcoming decades for this area will accentuate groundwater level fluctuations thus more ground deformations are foreseen and have to be considered in future infrastructure design standards

early landslide detection based on d insar technique at the wudongde hydropower reservoir differential interferometric synthetic aperture radar d insar is an effective technique for geological disaster prevention in this study an early landslide detection method based on the d insar technique is proposed first sar data are used to monitor landslide deformation in the wudongde hydropower reservoir then the ground surface deformation is derived by using the d insar technique which includes both soil and non soil deformations to remove the non soil deformation the vertical deformation filter method and slope deformation filter method for extracting soil deformations are established and the error of d insar results is analyzed by using two methods permanent scatterers insar and corner reflectors insar the errors caused by the atmospheric delay are estimated quantitatively additionally the relationship between the landslide deformation and the temporal spatial change of the sar data is analyzed and used to identify landslides finally through the proposed early landslide detection method the jinpingzi landslide in the wudongde area is investigated to validate the proposed method and the other unknown landslides in this area are also detected which demonstrates the feasibility of the proposed method the outcome of this study provides a workable technique for early detection of landslides in valley reservoir areas

quantification of everest region glacier velocities between 1992 and 2002 using satellite radar interferometry and feature tracking quincey d j luckman a benn d 2009 quantification of everest region glacier velocities between 1992 and 2002 using satellite radar interferometry and feature tracking journal of glaciology 55 192 596 606 sponsorship knowledge transfer project no 3742

strain accumulation in southern california 1973 1980 frequent surveys of seven trilateration networks in southern california over the interval 1973 1980 suggest that a regional increment in strain may have occurred in 1978 1979 prior to 1978 and after late 1979 the strain accumulation has been predominantly a uniaxial north south compression this secular trend was interrupted sometime in 1978 1979 by an increment in both north south and east west extension in five of the seven networks the onset of this change appears to have occurred first in the networks farthest south the changes occurred without any unusual seismicity within the networks but the overall seismicity in southern california was unusually low prior to and has been unusually high since the occurrence the average principal strain rates for the seven networks in the 1973 1980 interval are 0 17 mu strain yr north south contraction and 0 08 mu strain yr east west extension although the observed increment in strain could be related to unidentifed systematic error in the measuring system a careful review of the measurements and comparison with three other measuring systems reveal no appreciable cumulative systematic error

permanent scatterers analysis for atmospheric correction in ground based sar interferometry ground based synthetic aperture radar gb sar interferometry has already been recognized as a powerful tool complementary or alternative to spaceborne sar interferometry for terrain monitoring and for detecting structural changes in buildings it has been noted that in spite of the very short range compared with the satellite configuration in gb sar measurement the disturbances due to atmospheric effects cannot be neglected either the analysis of the interferometric phases of very coherent points called permanent scatterers pss allows the evaluation of the atmospheric disturbance and the possibility of removing it in this paper the ps analysis is carried out both on a test site facility and on a real campaign citrin valley italy that provided data with a temporal baseline of about ten months

early detection and in depth analysis of deformation phenomena by radar interferometry the spaceborne differential interferometric sar dinsar technique may play an important role in the measurement of land deformation phenomena especially in urban suburban and industrial areas in order to fully exploit its capabilities and increase its operational use this paper proposes two complementary levels of analysis the first one is a comprehensive low cost screening in order to detect unknown subsidence phenomena over large areas this analysis which provides a first estimation of deformations has to be performed using a limited set of sar images the second level includes an in depth quantitative analysis based on large image stacks which requires more image acquisition and data processing resources this type of analysis can be typically used to study deformation phenomena of special relevance where an effective support to the decision makers requires a fully quantitative estimation of deformations the two levels of analysis were proved over a test area of about 2100 km2 located in catalonia spain without any a priori information seven deformation phenomena were detected using a reduced set of interferograms furthermore the quantitative analysis of a subsidence of small spatial extent which was based on ascending and descending datasets confirmed the capability of dinsar to quantitatively assess deformation phenomena

ground deformation monitoring by using the permanent scatterers technique the example of the oltrepo pavese lombardia italy abstract the applicability of the permanent scatterers synthetic aperture radar interferometry psinsar technique for detecting and monitoring ground displacements was tested in the oltrepo pavese territory northern italy southern lombardia which could be representative of similar geological contexts in the italian apennines the study area which extends for almost 1100 km 2 is characterized by a complex geological and structural setting and the presence of clay rich sedimentary formations these characteristics make the oltrepo pavese particularly prone to several geological hazards shallow and deep landslides subsidence and swelling shrinkage of the clayey soils the psinsar technique used in this study overcomes most of the limitations of conventional interferometric approaches by identifying within the area of interest a set of radar benchmarks ps where very precise displacement measurements can be carried out more than 90 000 ps were identified by processing synthetic aperture radar sar images acquired from 1992 to 2001 by the european remote sensing satellites ers the psinsar application at a sub regional scale detected slow ground deformations ranging from 5 to 16 mm year and resulting from various processes landslides swelling shrinkage of clay soils and water pumping the ps displacements were analysed by collecting data obtained through geological geomorphologic field surveys geotechnical analysis of the soils and the information was integrated within a landslide inventory and the damaged building inventory despite the limited number of landslide bodies with ps 7 of the inventoried landslides the ps data helped to revise the state of activity of several landslides furthermore some previously unknown unstable slopes were detected two areas of uplift and two areas of subsidence were identified

atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps interferogram images derived from repeat pass spaceborne synthetic aperture radar systems exhibit artifacts due to the time and space variations of atmospheric water vapor other tropospheric variations such as pressure and temperature also induce distortions but the effects are smaller in magnitude and more evenly distributed throughout the interferogram than the wet troposphere term spatial and temporal changes of 20 in relative humidity lead to 10 cm errors in deformation products and perhaps 100 m of error in derived topographic maps for those pass pairs with unfavorable baseline geometries in wet regions such as hawaii these are by far the dominant errors in the spaceborne imaging radar c and x band synthetic aperature radar sir c x sar interferometric products the unknown time delay from tropospheric distortion is independent of frequency and thus multiwavelength measurements such as those commonly used to correct radar altimeter and global positioning system gps ionospheric biases cannot be used to rectify the error in the topographic case the errors may be mitigated by choosing interferometric pairs with relatively long baselines as the error amplitude is inversely proportional to the perpendicular component of the interferometer baseline for the sir c x sar hawaii data we found that the best longest baseline pair produced a map supporting 100 m contouring whereas the poorest baseline choice yielded an extremely noisy topographic map even at this coarse contour interval in the case of deformation map errors the result is either independent of baseline parameters or else very nearly so here the only solution is averaging of independent interferograms so in order to create accurate deformation products in wet regions many multiple passes may be required rules for designing optimal data acquisition and processing sequences for interferometric analyses in nondesert parts of the world are 1 to use the longest radar wavelengths possible within ionospheric scintillation and faraday rotation limits 2 for topography maximize interferometer baseline within decorrelation limits and 3 for surface deformation use multiple observations and average the derived products following the above recipe yields accuracies of 10 m for digital elevation models and 1 cm for deformation maps even in very wet regions such as hawaii

crustal deformation and source models of the yellowstone volcanic field from geodetic data 1 geodetic observations comprising interferometric synthetic aperture radar insar global positioning system gps and precision leveling measurements are used to infer volume change in the subsurface associated with the dynamics of the yellowstone volcanic system we focus primarily on the yellowstone caldera and its related magmatic hydrothermal and fault systems it appears that known faults play a significant role in controlling crustal volume increases and decreases due to the migration of volcanic and hydrothermal fluids for example over 5 cm of subsidence from 1992 to 1995 is associated with source volume changes 6 10 km beneath the nw trending elephant back fault zone and a north trending fault cutting across the caldera furthermore we are able to image an episode of fluid intrusion near the northern edge of the caldera the intrusion is elongated in the north south direction and is parallel to the north trending volume decrease the primary intrusion and related hydrothermal activity occurred between 1996 and 2000 though the volume changes appear to have continued shallowed and changed shape between 2000 and 2002 there is evidence that the intrusive activity influenced extensional faults to the north of the caldera

landslide activity as a geoindicator in italy significance and new perspectives from remote sensing landslide activity responds to rapid environmental changes and represents a relevant geoindicator in mountainous or hilly areas this paper discusses the socio economic relevance of landslide hazard in italy and the problems encountered in establishing relationships between landslide frequency climate and vegetation changes at different time scales landslides blocking a river channel have been carefully taken into account because they are usually characterized by high intensities with regard to the involved masses and movement velocities and their occurrences are often datable via radiocarbon dating this is due to the recovery of organic matter in the landslide dammed lakes for these reasons they can be considered important geoindicators in the wider category of slope failures the marked effects of the anthropogenic activity on slope instability processes in the last 50 years are discussed with reference to two case histories the chianti hills in tuscany and the cinque terre national park in liguria finally two novel techniques of remote sensing are proposed as tools for a systematic monitoring of slope instability at different time and spatial scales both techniques are based on the interferometric synthetic aperture radar sar technology and differ on the type of platform satellite and ground based used to acquire data

phase gradient approach to stacking interferograms the phase gradient approach is used to construct averages and differences of interferograms without phase unwrapping our objectives for change detection are to increase fringe clarity and decrease errors due to tropospheric and ionospheric delay by averaging many interferograms the standard approach requires phase unwrapping scaling the phase according to the ratio of the perpendicular baseline and finally forming the average or difference however unique phase unwrapping is usually not possible since the phase gradient due to topography is proportional to the perpendicular baseline phase unwrapping is unnecessary prior to averaging or differencing phase unwrapping may be needed to interpret the results but it is delayed until all of the largest topographic signals are removed we demonstrate the method by averaging and differencing six interferograms having a suite of perpendicular baselines ranging from 18 to 406 m cross spectral analysis of the difference between two tandem interferograms provides estimates of spatial resolution which are used to design prestack filters a wide range of perpendicular baselines provides the best topographic recovery in terms of accuracy and coverage outside of mountainous areas the topography has a relative accuracy of better than 2 m residual interferograms single interferogram minus stack have tilts across the unwrapped phase that are typically 50 mm in both range and azimuth reflecting both orbit error and atmospheric delay smaller scale waves with amplitudes of 15 mm are interpreted as atmospheric lee waves a few global positioning system gps control points within a frame could increase the precision to 20 mm for a single interferogram further improvements may be achieved by stacking residual interferograms

ice velocity determined using conventional and multiple aperture insar abstract we combine conventional interferometric synthetic aperture radar insar and multiple aperture insar mai to determine the ice surface velocity on the langjokull and hofsjokull ice caps in iceland in 1994 this approach allows the velocity of the 20 principal ice cap outlet glaciers to be fully resolved we show that mai leads to displacement estimates of finer resolution 15 versus 150 m and superior precision 5 versus 15 cm to those afforded by the alternative technique of speckle tracking using sar data acquired in ascending and descending orbits we show that ice flows within 15 of the direction of maximum surface slope across 66 of the ice caps it is therefore possible to determine ice displacement over the majority of the ice caps using a single sar image pair thereby reducing errors associated with temporal fluctuations in ice flow

form and flow of the devon island ice cap canadian arctic 1 in this study 3370 km of 100 mhz ice penetrating radar data were acquired from devon ice cap arctic canada bed returns were obtained from 90 of flight tracks mean crossing point errors in ice surface elevation and ice thickness were 7 8 m digital elevation models of ice cap surface and bed elevation and ice thickness were produced and can be used as boundary conditions in numerical modeling devon ice cap including 1960 km2 of contiguous stagnant ice to its west is 14 010 km2 the ice cap proper is 12 050 km2 its largest drainage basin is 2630 km2 the ice cap crest has a maximum measured elevation of 1921 m maximum recorded ice thickness is 880 m ice cap volume is 3980 km3 about 10 mm sea level equivalent the bed 8 of which lies below sea level is an upland plateau dissected by steep sided valleys that control the locations of the major outlet glaciers which dominate ice cap drainage about 73 km 4 of the ice cap margin ends in tidewater the margin is not floating icebergs of 100 m were observed offshore only a single outlet glacier showed signs of past surge activity all major outlet glaciers along the eastern ice cap margin have retreated 1 3 km since 1960 synthetic aperture radar interferometic velocity structure shows slow undifferentiated flow predominating in the west and center with fast flowing outlet glaciers and intervening slow flowing ridges typical elsewhere outlet glacier velocities are 7 10 times higher than in areas of undifferentiated flow this velocity structure of fast flowing units within slower flowing ice appears typical of many large 103 km2 arctic ice caps

insar mapping and modelling of an active iranian salt extrusion abstract interferometric synthetic aperture radar insar imaging is a powerful technique that is increasingly used to extract detailed 3d information on earth surface structures including exposed diapiric surfaces we have used insar to map for the first time the cumulative surface deformation in a 6 km 6 km region surrounding an active salt diapir syahoo and its associated surface salt flow or namakier exposed in the zagros mountains of southern iran images provided by the european space agency were acquired in 12 increments ranging in length between 35 and 1248 days over a 14 year interval between july 1992 and may 2006 the deformation of the salt surface is non steady with extrapolated rates of displacement varying between surficial uplifts of 1 4 mm day 1 511 mm a 1 and subsidence of 2 2 mm day 1 803 mm a 1 growth of a central topographic dome occurs following short wet intervals to create a salt fountain morphology which then slowly decays during the intervening long dry periods salt associated with dynamic bulging9 of the central dome during wet intervals may flow laterally via gravity spreading into the surrounding salt sheet resulting in deflation and subsidence of the dome which is counteracted by growth and inflation of the adjacent namakiers salt bulges9 that migrate down the namakier resulting in local inflationary and deflationary cycles of the surface may be regarded as episodic pulses of gravity spreading areas of inflation and deflation are also observed to commonly reverse during dry to wet periods as the overall salt system effectively self regulates as it continually strives for dynamic equilibrium as long as the source of salt remains undepleted gravity spreading of the dome ultimately results in more buoyancy driven salt flowing up the diapiric neck to replenish and feed the extrusion and maintain the gross fountain morphology

integrated satellite interferometry tropospheric noise gps estimates and implications for interferometric synthetic aperture radar products interferometric synthetic aperture radar insar like other astronomic and space geodetic techniques is limited by the spatially and temporally variable delay of electromagnetic waves propagating through the neutral atmosphere statistical analysis of these variations from a wide variety of instruments reveals a power law dependence on frequency that is characteristic of elementary kolmogorov turbulence a statistical model for a major component of the delay fluctuations the wet component has previously been developed by treuhaft and lanyi 1987 for very long baseline interferometry a continuous global positioning system gps network is now in place in southern california that allows estimation of along with geodetic parameters the total delay due to the atmosphere above each site on a subhourly basis these measurements are shown to conform to the treuhaft and lanyi tl statistical model both temporally and spatially the tl statistical model is applied to the problem of insar and used to produce the covariance between two points separated in time and or space the error due to the atmospheric variations for sar products such as topography and surface deformation is calculated via propagation of errors there are two methods commonly cited to reduce the effect of atmospheric distortion in products from sar interferometry stacking and calibration stacking involves averaging independent interferograms to reduce the noise calibration involves removing part or all of the delay using data from an independent source such as total zenith delay estimates from continuous gps networks despite the relatively poor spatial density of surface measurements calibration can be used to reduce noise if the measurements are sufficiently accurate reduction in tropospheric noise increases with increasing number of measurement points and increasing accuracy up to a maximum of n where n is the number of points stacking and calibration are shown to be complementary and can be used simultaneously to reduce the noise to below that achievable by either method alone

dem error correction in insar time series we present a mathematical formulation for the phase due to the errors in digital elevation models dems in synthetic aperture radar sar interferometry insar time series obtained by the small baseline sb or the small baseline subset method we show that the effect of the dem error in the estimated displacement is proportional to the perpendicular baseline history of the set of sar acquisitions this effect at a given epoch is proportional to the perpendicular baseline between the sar acquisition at that epoch and the reference acquisition therefore the dem error can significantly affect the time series results even if sb interferograms are used we propose a new method for dem error correction of insar time series which operates in the time domain after inversion of the network of interferograms for the displacement time series this is in contrast to the method of berardino 2002 in which the dem error is estimated in the interferogram domain we show the effectiveness of this method using simulated insar data we apply the new method to fernandina volcano in the galapagos islands and show that the proposed dem error correction improves the estimated displacement significantly

mapping landslide surface displacements with time series sar interferometry by combining persistent and distributed scatterers a case study of jiaju landslide in danba china abstract insar technology provides a powerful tool to detect potentially unstable slopes across wide areas and to monitor surface displacements of a single landslide however conventional time series insar methods such as persistent scatterer interferometry psi and small baseline subset sbas can rarely identify sufficient measurement points mps in mountainous areas due to decorrelations caused by steep terrain and vegetation coverage in this study we developed a new insar approach coherent scatterer insar csi to map landslide surface displacements in the radar line of sight los direction by combining persistent scatterers ps and distributed scatterers ds the key ideas of csi include the employment of the generalized likelihood ratio glr test for the identification of statistically homogeneous pixels shps and the use of the phase linking algorithm to estimate optimal phase for each ds pixel the joint exploitation of ps and ds targets dramatically increases the spatial density of mps which makes the phase unwrapping more reliable to demonstrate the effectiveness of the csi approach we applied it to retrieve the historical displacements of the jiaju landslide in danba county of southwest china using 19 l band alos palsar images 2006 2011 and nine c band envisat asar images 2007 2008 multiple comparisons clearly illustrated the big advantages of csi over psi and sbas in mapping landslide displacements with more details owing to much higher 10 times mp density furthermore the superiority of l band sar data over c band for landslide investigation in rural environments was confirmed quantitative validation of the csi results for palsar data against in situ gps measurements suggested an accuracy of about 10 5 millimeters per year mm year in terms of root mean square error rmse afterwards the spatial temporal characteristics of the jiaju landslide surface displacements were summarized with a new upper boundary for the active northern part delineated particularly the northern part of the landslide moved faster than the southern part exhibiting a maximum los displacement rate of around 120 mm year subsequently the fluvial erosion by the dajinchuan river was identified as the predominant impact factor for the instability of the jiaju landslide finally the major problems and challenges for the application of csi method were discussed and the conclusions were given

analysis of ground deformation detected using the sbas dinsar technique in umbria central italy ground deformation affecting the umbria region central italy in the 9 year period from 1992 to 2000 was investigated through multi temporal differential synthetic aperture radar interferometry dinsar for the purpose the small baseline subset sbas technique was adopted which allows studying the temporal evolution of the detected deformation at two spatial scales a low resolution regional scale and a full resolution local scale for the analysis sar data acquired by the european remote sensing ers 1 2 satellites along ascending and descending orbits were used the detected deformation was analysed to investigate its relevance to geophysical geomorphologic and human induced processes that may result in hazardous conditions to the population of umbria low resolution deformation data were used to i determine the amount of displacement caused by the umbria marche earthquake sequence from september 1997 to april 1998 in the foligno area ii determine the number and percentage of the known landslides that can be monitored by the dinsar technology in the investigated area and iii identify and measure subsidence induced by exploitation of a confined aquifer in the valle umbra results indicate that earthquakes moved through the foligno area westwards up to 3 9 cm and with an uplift reaching 1 7 cm intersection in a gis of the low resolution deformation maps with a detailed landslide inventory map allowed the determination that the portion of landslides that can be monitored by the sbas dinsar technique in umbria ranges from 2 7 to 3 4 and the percentage of the total landslide area ranges from 10 4 to 12 8 in the valle umbra a dependency was found between the time and the amount of detected ground deformation and the record of water withdrawal the full resolution deformation data were used to investigate the movement of the ivancich landslide in the assisi municipality joint analysis of the spatial and the temporal characteristics of the ground displacement allowed the formulation of a hypothesis on the landslide geometry and deformation pattern

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

interferometric radar measurement of ocean surface currents we describe a new method of measuring surface currents using an interferometric synthetic aperture radar an airborne implementation has been tested over the san francisco bay near the time of maximum tidal flow resulting in a map of the east west component of the current only the line of sight component of velocity is measured by this technique where the signal to noise ratio was strongest statistical fluctuations of less than 4 cm s 1 were observed for ocean patches of 60 60 m

land subsidence in the cerro prieto geothermal field baja california mexico from 1994 to 2005 an integrated analysis of dinsar leveling and geological data abstract cerro prieto is the oldest and largest mexican geothermal field in operation and has been producing electricity since 1973 the large amount of geothermal fluids extracted to supply steam to the power plants has resulted in considerable deformation in and around the field the deformation includes land subsidence and related ground fissuring and faulting these phenomena have produced severe damages to the local infrastructure such as roads irrigation canals and other facilities in this paper the technique of differential synthetic aperture radar interferometry dinsar is applied using c band envisar asar data acquired between 2003 and 2006 to determine the extent and amount of land subsidence in the mexicali valley near cerro prieto geothermal field the dinsar results were compared with published data from precise leveling surveys 1994 1997 and 1997 2006 and detailed geological information in order to improve understanding of the temporal and spatial distributions of anthropogenic subsidence in the mexicali valley the leveling and dinsar data were modeled to characterize the observed deformation in terms of fluid extraction the results confirm that the tectonic faults control the spatial extent of the observed subsidence these faults likely act as groundwater flow barriers for aquifers and reservoirs the shape of the subsiding area coincides with the cerro prieto pull apart basin in addition the changes in spatial pattern and rate of the subsidence are correlated with the development of the cerro prieto geothermal field

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

a statistical study of ionospheric profile parameters derived from millstone hill incoherent scatter radar measurements diurnal seasonal and solar activity variations of the bottomside electron density profile parameters b0 and b1 representing the f2 layer thickness and shape are studied using a large incoherent scatter radar dataset for millstone hill covering the period 1976 2002 these results are compared with the latest iri model our statistical study is characterized by morning and afternoon falls in the diurnal variation of b0 for seasons other than summer and a similar to 15 change in b1 over a solar cycle features not fully well represented by the standard iri model the standard iri b1 however is very close to observations in terms of the diurnal variation

the lowest place on earth is subsiding an insar interferometric synthetic aperture radar perspective since the early 1990s sinkholes and wide shallow subsidence features wssfs have become major problems along the dead sea shores in israel and jordan sinkholes are readily observed in the field but their locations and timing are unpredictable wssfs are often difficult to observe in the field however once identified they delineate zones of instability and increasing hazard in this study we identify characterize and measure rates of subsidence along the dead sea shores by the interferometric synthetic aperture radar insar technique we analyze 16 sar scenes acquired during the years 1992 to 1999 by the european remote sensing ers 1 and ers 2 satellites the interferograms span periods of between 2 and 71 months wssfs are observed in the lisan peninsula and along the dead sea shores in a variety of appearances including circular and elongate coastal depressions a few hundred meters to a few kilometers in length depressions in ancient alluvial fans and depressions along salt diapir margins phase differences measured in our interferograms correspond to subsidence rates generally in the range of 0 20 mm yr within the studied period with exceptional high rates that exceed 60 mm yr in two specific regions during the study period the level of the dead sea and of the associated ground water has dropped by 6 m this water level drop within an aquifer overlying fine grained marly layers would be expected to have caused aquifer system consolidation resulting in gradual subsidence comparison of our insar observations with calculations of the expected consolidation shows that in areas where marl layers are known to compose part of the upper 30 m of the profile estimated consolidation settlements are of the order of the measured subsidence our observations also show that in certain locations subsidence appears to be structurally controlled by faults seaward landslides and salt domes gradual subsidence is unlikely to be directly related to the sinkholes excluding the use of the wssfs features as predictable precursors to sinkhole formation

flow dynamics and iceberg calving rates of devon ice cap nunavut canada the surface velocity field of devon ice cap nunavut canada was mapped using interferometric synthetic aperture radar insar ascending european remote sensing satellite 1 and 2 ers 1 2 tandem mode data were used for the western and southeast sectors and 3 day repeat pass ers 1 imagery for the northeast sector speckle tracking procedures were used with radarsat 1 imagery to obtain surface velocities over the terminus of belcher glacier a major calving front where decorrelation between ers data occurred the insar data highlight a significant contrast in ice flow dynamics between the east and west sides of the ice cap ice movement west of the main north south divide is dominated by relatively uniform sheet flow but three fast flowing outlet glaciers that extend 14 23 km beyond the ice cap margin also drain this region several outlet glaciers that extend up to 60km inland from the eastern margin drain the eastern side of the ice cap the dominant ice flow regimes were classified based on the relationship between the driving stress averaged over a length scale often ice thicknesses and the ratio of surface velocity to ice thickness the mapped distribution of flow regimes appears to depict the spatial extent of basal sliding across the ice cap this is supported by a close relationship between the occurrence of flow stripes on the ice surface and flow regimes where basal sliding was found to be an important component of the glacier motion iceberg calving rates were computed using measured surface velocities and ice thicknesses derived from airborne radio echo sounding the volume of ice calved between 1960 and 1999 was estimated to be 20 5 4 7 km 3 or 0 57km 3 a 1 approximately 89 of this loss occurred along the eastern margin the largest single source is belcher glacier which accounts for 50 of the total amount of ice calved

geodetic airborne laser altimetry of breidamerkurjokull and skeidararjokull iceland and jakobshavns isbrae west greenland

combined gps and insar models of postseismic deformation from the northridge earthquake models of combined global positioning system gps and interferometric synthetic aperture radar insar data collected in the region of the northridge earthquake indicate that significant afterslip on the main fault occurred following the earthquake additional shallow deformation occurred to the west of the main rupture plane both data sets are consistent with logarithmic time dependent behavior following the earthquake indicative of afterslip rather than postseismic relaxation aftershocks account for only about 10 of the postseismic motion the two data sets are complimentary in determining the postseismic processes fault afterslip and shallow deformation dominate the deformation field in the two years following the earthquake lower crustal deformation may play an important role later in the earthquake cycle

poroelastic rebound along the landers 1992 earthquake surface rupture maps of surface displacement following the 1992 landers california earthquake generated by interferometric processing of ers 1 synthetic aperture radar sar images reveal effects of various postseismic deformation processes along the 1992 surface rupture the large scale pattern of the postseismic displacement field includes large lobes mostly visible on the west side of the fault comparable in shape with the lobes observed in the coseismic displacement field this pattern and the steep displacement gradient observed near the emerson camp rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture models show that horizontal slip occurring on a buried dislocation in a poisson s material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10 20 km from the fault yet this pattern is not observed in the postseismic interferograms as previously proposed to explain local strain in the fault step overs peltzer et al 1996b we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault compensating the vertical ground shift produced by fault afterslip such a rebound is explained by the gradual change of the crustal rocks poisson s ratio value from undrained coseismic to drained postseismic conditions as pore pressure gradients produced by the earthquake dissipate using the poisson s ratio values of 0 27 and 0 31 for the drained and undrained crustal rocks respectively elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the sar data and the horizontal displacement measured at global positioning system gps sites along a 60 km long transect across the emerson fault savage and svarc 1997 using a detailed surface slip distribution on the homestead valley kickapoo and johnson valley faults we modeled the poroelastic rebound in the homestead valley pull apart a poisson s ratio value of 0 35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3 5 years following the earthquake this large value implies a seismic velocity ratio vp vs of 2 1 consistent with the observed low vs values of fault zone guided waves at shallow depth li et al 1997 the sar data also reveal postseismic creep along shallow patches of the eureka peak and burnt mountain faults with a characteristic decay time of 0 8 years coseismic dilatant hardening locking process followed by post seismic pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault

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

estimation of rock glacier surface deformation using sar interferometry data the detection and quantification of surface deformation of an active rock glacier using the differential synthetic aperture radar sar interferometry d insar technique is presented an average deformation rate of 6 mm 35 days in the radar line of ight was estimated for the summer of 1992 the maximum deformation rate 18 mm 35 days was identified at the upper part of the rock glacier whereas the deformation rate at the snout of the rock glacier was about 10 mm 35 days the spatial distribution of the surface deformation in the d insar displacement map is smooth and supports the idea that ice is the stress transferring medium in rock glaciers

spatial analysis of surface deformation distribution detected by persistent scatterer interferometry in lanzhou region china persistent scatterer synthetic aperture radar interferometry ps insar is a remote sensing method that can be used to detect surface deformation which is an indicator of potential geohazards by capturing such deformations over time it is possible to obtain valuable information regarding geohazards such as landslides this study focused on the use of ps insar to investigate the distribution and causes of surface deformation in the lanzhou region of gansu province in china between 2003 and 2010 41 advanced synthetic aperture radar images were captured by the envisat satellite and analyzed using ps insar and the correlation between the observed surface deformation and topographic geologic and anthropogenic factors was derived based on a geographic information system platform it was found that the largest number and highest density of surface deformations occurred at elevations of 1486 1686 m it was also established that slope ranges of 25 30 and 35 40 are threshold values at which surface deformation changes abruptly and that slopes with north and northwest aspects are most prone to surface deformation the lithologies most susceptible to surface deformation are clay sandy soil and loess the normalized difference vegetation index indicated that surface deformation occurred most often in areas with sparse vegetation anthropogenic activities e g construction and wastewater discharge could be inferred as causal mechanisms of surface deformation comparison of the distributions of geohazards and surface deformation showed considerable consistency which proves surface deformation can induce geohazards these results could help governments improve urban planning and geohazard mitigation

ice shelf dynamics near the front of the filchner ronne ice shelf antarctica revealed by sar interferometry fifteen synthetic aperture radar sar images of the ronne ice shelf also referred to as the filchner ronne ice shelf antarctica obtained by the european remote sensing satellites ers 1 and 2 are used to study ice shelf dynamics near two ends of the iceberg calving front interferograms constructed from these sar images are used to resolve the ice shelf displacement along several directions in response to both ocean tide and long term creep flow tidal motion is separated from creep flow using differential interferometry i e two or more interferograms in which fringe patterns common to all are predominantly associated with creep flow creep flow velocities thus determined compare well with prior ice shelf velocity surveys using these data we studied the influence of large scale rifts ice rises and coastal separation on the ice shelf flow many of the large rifts that appear to form the boundaries where tabular icebergs may eventually detach from the ice shelf are filled with a melange of sea ice ice shelf debris and wind blown snow the interferograms show that this melange tends to deform coherently in response to the ice shelf flow and has sufficient strength to trap large tabular ice shelf fragments for several decades before the fragments eventually become icebergs in many instances the motion of the tabular fragments is a rigid body rotation about a vertical axis that is driven by velocity shear within the melange the mechanical role of the rift filling melange may be to bind tabular ice shelf fragments to the main ice shelf before they calve this suggests two possible mechanisms by which climate could influence tabular iceberg calving first spatial gradients in oceanic and atmospheric temperature may determine where the melange melts and thus the location of the iceberg calving margin second melting or weakening of ice melange as a consequence of climate change could trigger a sudden or widespread release of tabular icebergs and lead to rapid ice shelf disintegration

deformation time series generation in areas characterized by large displacement dynamics the sar amplitude pixel offset sbas technique we exploit the amplitude information of a sequence of synthetic aperture radar sar images acquired at different times in order to generate displacement time series in areas characterized by large and or rapid deformation the size of which is on the order of the image s pixel dimensions we follow the same rationale of the small baseline subset sbas differential sar interferometry dinsar approach by coupling the available sar images into pairs characterized by a small separation between the acquisition orbits we exploit the amplitudes of the selected image pairs in order to calculate the relative across track range and along track azimuth pixel offsets po finally we apply the sbas inversion strategy to retrieve the range and azimuth displacement time series this approach referred to as pixel offset po sbas technique has been applied to a set of 25 envisat sar observations of the sierra negra caldera galapagos islands spanning the 2003 2007 time interval the retrieved deformation time series show the capability of the technique to detect and measure the large displacements affecting the inner part of the caldera that in correspondence to the october 2005 eruption reached several meters moreover by comparing the po sbas results to continuous gps measurements we estimate that the accuracy of the po sbas time series is on the order of 1 30th of a pixel for both range and azimuth directions

topographic mapping from interferometric synthetic aperture radar observations interferometric synthetic aperture radar observations provide a means for obtaining high resolution topographic terrain maps from data acquired simultaneously at two slightly displaced antennas calculation of the three dimensional coordinates of all the points in a radar image can be made from the combination of along track slant range and interferometer fringe measurements thus the result compensates for the layover exhibited by conventional radar maps and removes the false targets generated by multiple signal paths to a given object in the scene we have derived a topographic map of a portion of the san francisco bay area utilizing data that were recorded by a radar system mounted on a nasa cv990 aircraft and processed by a general purpose digital computer this map displays the height above sea level of a region approximately 11 km by 10 km in size sampled on an 11 m pixel grid uncertainties in the estimated height result from imprecise knowledge of the observational geometry radar bandwidth limitations and finite signal to noise ratios for our system each factor contributed approximately 3 4 m to a relative rms error of about 6 m additional global errors of the order of 10 m can result from inaccurate monitoring of aircraft attitude the statistical variation of the height measurements from the ocean portion of the map which we can presume to be quite flat is 2 10 m rms which is consistent with the theoretical estimate comparison of our results with u s geological survey contour maps indicates a high degree of correlation between the two sets of altitude data

actively growing anticlines beneath catania from the distal motion of mount etna s decollement measured by sar interferometry and gps we used sar interferometry and gps to measure the rise of an anticline beneath the urban area of catania italy which originates from outward thrusting above the basal decollement of etna volcano the anticline grows at a rate of over 0 015 m a 1 of se thrusting and about 0 007 m a 1 of relative uplift by relating this growth to the simultaneous extension measured at the summit of the volcano we demonstrate the occurrence of active volcanic spreading this process may exert a control on eruptions earthquakes sector collapses and landslides

estimation of small surface displacements in the upper rhine graben area from a combined analysis of ps insar levelling and gnss data the intra plate deformation of the upper rhine graben urg located in central europe is investigated using geodetic measurement techniques we present a new approach to calculate a combined velocity field from insar levelling and gnss measurements as the expected tectonic movements in the urg area are small less than 1mma 1 the best possible solutions for linear velocity rates from single technique analyses are estimated in a first step second we combine the velocity rates obtained from insar line of sight velocity rates in ascending and descending image geometries levelling vertical velocity rates and gnss horizontal velocity rates using least squares adjustment lsa focusing on the northern urg area we analyse sar data on four different image stacks ers ascending ers descending envisat ascending envisat descending using the persistent scatterer ps approach the linear velocity rates in ascending and descending image geometries respectively are estimated in an lsa from joint time series analysis of ers and envisat data vertical velocity rates from levelling are obtained from a consistent adjustment of more than 40 000 measured height differences using a kinematic displacement model horizontal velocity rates in east and north direction are calculated from a time series analysis of daily coordinate estimates at 76 permanently operating gnss sites in the urg region as the locations at which the measurement data of ps insar levelling and gnss reside do not coincide spatial interpolation is needed during several steps of the rigorous processing we use ordinary kriging to interpolate from a given set of data points to the locations of interest with a special focus on the modeling and propagation of errors the final 3 d velocity field is calculated at a 200 m grid which carries values only close to the location of ps points resulting in a mean horizontal and vertical precision of 0 30 and 0 13 mm a 1 respectively the vertical component of the combined velocity field shows a significant subsidence of about 0 5 mm a 1 in the northern part of the graben coinciding with a well known quaternary basin structure horizontal displacement rates of up to 0 8 mm a 1 in southeast direction are observed outside the graben in reasonable alignment with the average direction of maximum horizontal stress within the graben the velocity directions rotate toward east in the non subsiding part while an opposite trend is observed in the subsiding part of the graben the complexities of the observed velocity field are compatible to the geomechanical situation in our investigation area which is characterized by a transition from a restraining to a releasing bend setting glacial isostatic adjustment is another potential source influencing the observed velocity field as well as anthropogenic signals due to mining oil exploration and groundwater usage that have been identified in some places

land subsidence caused by the east mesa geothermal field california observed using sar interferometry interferometric combination of pairs of synthetic aperture radar sar images acquired by the ers 1 satellite maps the deformation field associated with the activity of the east mesa geothermal plant located in southern california sar interferometry is applied to this flat area without the need of a digital terrain model several combinations are used to ascertain the nature of the phenomenon short term interferograms reveal surface phase changes on agricultural fields similar to what had been observed previously with seasat radar data long term 2 years interferograms allow the study of land subsidence and improve prior knowledge of the displacement field and agree with existing sparse levelling data this example illustrates the power of the interferometric technique for deriving accurate industrial intelligence as well as its potential for legal action in cases involving environmental damages

monitoring very slow slope movements by means of sar interferometry a case study from a mass waste above a reservoir in the otztal alps austria the application of radar interferometry to detect slope movements on the order of millimeters to centimeters per year is demonstrated the deformation field of a slope above a hydropower reservoir in the austrian alps was derived from sar interferometric pairs of the satellites ers 1 and ers 2 acquired between july 1992 and august 1998 above the treeline it was possible to map the motion from interferograms covering time spans up to three years whereas ground based geodetic measurements focussed on the movements of the lower slope section the interferometric analysis shows that the mass wasting processes affect the entire slope of 1000 m vertical extent significant interannual differences of the displacement rates became apparent which were found to be related to the pattern of summer rainfall

radar interferogram filtering for geophysical applications the use of sar interferometry is often impeded by decorrelation from thermal noise temporal change and baseline geometry power spectra of interferograms are typically the sum of a narrow band component combined with broad band noise we describe a new adaptive filtering algorithm that dramatically lowers phase noise improving both measurement accuracy and phase unwrapping while demonstrating graceful degradation in regions of pure noise the performance of the filter is demonstrated with sar data from the ers satellites over the jakobshavns glacier of greenland

deformation and seismicity in the coso geothermal area inyo county california observations and modeling using satellite radar interferometry interferometric synthetic aperture radar insar data collected in the coso geothermal area eastern california during 1993 1999 indicate ground subsidence over a 50 km 2 region that approximately coincides with the production area of the coso geothermal plant the maximum subsidence rate in the peak of the anomaly is 3 5 cm yr 1 and the average volumetric rate of subsidence is of the order of 10 6 m 3 yr 1 the radar interferograms reveal a complex deformation pattern with at least two irregular subsidence peaks in the northern part of the anomaly and a region of relative uplift on the south we invert the insar displacement data for the positions geometry and relative strengths of the deformation sources at depth using a nonlinear least squares minimization algorithm we use elastic solutions for a prolate uniformly pressurized spheroidal cavity in a semi infinite body as basis functions for our inversions source depths inferred from our simulations range from 1 to 3 km which corresponds to the production depths of the coso geothermal plant underpressures in the geothermal reservoir inferred from the inversion are of the order of 0 1 1 mpa except a few abnormally high underpressures that are apparently biased toward the small source dimensions analysis of the insar data covering consecutive time intervals indicates that the depths and or horizontal extent of the deformation sources may increase with time this increase presumably reflects increasing volumes of the subsurface reservoir affected by the geothermal exploitation we show that clusters of microearthquakes associated with the geothermal power operation may result from perturbations in the pore fluid pressure as well as normal and shear stresses caused by the deflation of the geothermal reservoir

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

coseismic and early post seismic slip associated with the 1999 izmit earthquake turkey from sar interferometry and tectonic field observations summary we use combined tectonic field observations and sar data to determine an improved model of the slip associated with the 1999 izmit earthquake which ruptured the north anatolian fault at the eastern end of the sea of marmara the leading goal is to understand the main features of the coseismic and post seismic deformation which are captured together in the sar data to achieve this we make a critical analysis of the ers1 2 sar data which allows atmospheric effects to be identified and removed we also use detailed field mapping and measurements of the earthquake surface rupture dislocations in elastic half space and a forward modelling strategy allow us to obtain a slip model by steps a trial and error approach is combined with conventional inversion techniques to determine the slip in the different regions of the fault the sar data are well explained with three main zones of high slip along the fault releasing a total moment of 2 3 10 20 nm mw 7 6 which is higher than the seismological estimates 1 7 2 0 10 20 n m the inhomogeneous slip distribution correlates with fault segments identified at the surface the izmit rupture appears to have extended 30 km west of the hersek peninsula into the sea of marmara with slip tapering fro m2mt ozero the western end of the rupture is located 40 km sse from istanbul we show that some features seen near to mudurnu and gevye and previously interpreted as slip on secondary faults are explained mostly as atmospheric effects correlated with the topography using our approach and the available gps data we obtain a slip model that represents the coseismic slip alone which suggests that the moment release during the main shock was 1 9 10 20 nm mw 7 5 consistent with the seismological estimates we conclude that the sar data include the effects of 2 m of fast after slip during the month following the main shock within a zone of the fault located 12 24 km below the epicentral region near the hypocentre at a depth of 18 km the fault appears to have experienced dynamic slip of 1 m associated with the main shock followed by 2 m of rapidly decelerating post seismic shear during the following month we suggest that the distribution of heterogeneous slip and loading along the different fault segments may be important factors controlling the propagation of large earthquake ruptures along the north anatolian fault

reaction of lake environment on the climatic cooling transition from the eemian interglacial to early vistulian on the basis of solniki palaeolake sediments ne poland abstract results of cladocera and organic matter analyses from the eemian lake deposits in solniki ne poland are presented we focused on the reaction of the lake s environment to progressive cooling reflected in the changes of the organic matter and the cladocera community approximately 3 m of sediments accumulated from the optimum of the eemian interglacial e5 to the early vistulian ev2 were analyzed results of analyses of organic matter correlated with cladocera data enabled a reconstruction of varying environmental conditions and geneses of organic matter the organic carbon content is closely related to climatic conditions the lowest values are recorded in the coldest periods during the ending of the eemian interglacial and the herning stadial of the early vistulian the highest ones occur unexpectedly in the eemian post optimum at the end of the spruce phase and at the beginning of the pine phase during the eemian optimum e5 biodiversity of cladocera species was maximal 17 but the frequency of specimens was low at the end of the optimum e5 an improvement of living conditions more optimal for various ecological groups of cladocera took place algae remained a significant component of organic matter the early eemian post optimum e6 is characterized as a period of the high number of zooplankton species the highest frequency as well as biodiversity at that time probably a higher water level and an optimal temperature formed very favorable conditions for the fauna development during this period maximal values of toc tn and the slight increase of c n ratio suggest a higher primary productivity with a small supply of organic matter of terrestrial origin coexistence of two species of the family bosminidae and maximum values of 15 n suggest a high trophy mesotrophy high resolution and detailed analyses of the lake sediments accumulated during the pollen zone e7 allowed to identity a few minor environmental fluctuations recorded in varied cladocera frequency toc and c n ratio in the final period of the eemian interglacial and in the early vistulian herning stadial a drastic reduction in frequency of cladocera and the features of om the lowest amounts of toc and tn and the highest value of c n reflect the decline of lake conditions and a lack of both primary production and aquatic vegetation in the lake during the amersfoort interstadial a renewed bioproduction and slightly supply of terrestrial plants in organic matter as well as revitalize zooplankton took place

sar monitoring of progressive and seasonal ground deformation using the permanent scatterers technique spaceborne differential radar interferometry has proven a remarkable potential for mapping ground deformation phenomena e g urban subsidence volcano dynamics coseismic and postseismic displacements along faults as well as slope instability however a full operational capability has not been achieved yet due to atmospheric disturbances and phase decorrelation phenomena these drawbacks can often be at least partially overcome by carrying out measurements on a subset of image pixels corresponding to natural or artificial stable reflectors permanent scatterers ps and exploiting long temporal series of interferometric data this approach allows one to push the measurement precision very close to its theoretical limit in the order of spl sim 1 mm for c band european remote sensing ers like sensors in this paper the detection of both time uniform and seasonal deformation phenomena is addressed and a first assessment of the precision achievable by means of the ps technique is discussed results highlighting deformation phenomena occurring in two test sites in california are reported fremont in the southern bay area and san jose in the santa clara valley

slip distribution of the 1997 umbria marche earthquake sequence joint inversion of gps and synthetic aperture radar interferometry data 1 in september october 1997 a seismic sequence struck the southern part of the umbria marche region in central italy the three largest earthquakes took place on 26 september 0033 ut mw 5 7 09 40 gu mw 6 0 near colfiorito and on 14 october mw 5 6 near sellano we present a slip model for the main sequence over two separate fault planes we use previously published global positioning system gps and differential interferometric synthetic aperture radar dinsar data as well as new dinsar data covering the sellano event as input to a simulated annealing solution for the distribution of fault slip the combination of both types of geodetic data as well as the combination of sar interferograms from both ascending and descending satellite tracks allows us to make a more complete geodetically based assessment of the fault slip for these earthquakes for the solutions presented here we required the slip distribution to match the observed surface displacements while maintaining a total scalar seismic moment approximately equal to the seismologically inferred moments for this sequence we find the slip was concentrated in a zone at 4 5 5 km depth with a shallower tail extending to near the surface at the northern end of the colfiorito sequence for both the second colfiorito event and the sellano event individual patches were calculated with maximum slip of 2 m implying stress drops as high as 20 mpa coulomb stress calculations show that the coseismic slip distribution for the colfiorito events facilitated slip of the adjacent sellano earthquake