shallow creep on the haiyuan fault gansu china revealed by sar interferometry interferometric synthetic aperture radar data are used to map the interseismic velocity field along the haiyuan fault system hfs at the north eastern boundary of the tibetan plateau two m 8 earthquakes ruptured the hfs in 1920 and 1927 but its 260 km long central section known as the tianzhu seismic gap remains unbroken since 1000 years the envisat sar data spanning the 2003 2009 period cover about 200 300 km2 along three descending and two ascending tracks interferograms are processed using an adapted version of roi pac the signal due to stratified atmospheric phase delay is empirically corrected together with orbital residuals mean line of sight velocity maps are computed using a constrained time series analysis after selection of interferograms with low atmospheric noise these maps show a dominant left lateral motion across the hfs and reveal a narrow 35 km long zone of high velocity gradient across the fault in between the tianzhu gap and the 1920 rupture we model the observed velocity field using a discretized fault creeping at shallow depth and a least squares inversion the inferred shallow slip rate distribution reveals aseismic slip in between two fully locked segments the average creep rate is 5 mm yr 1 comparable in magnitude with the estimated loading rate at depth suggesting no strain accumulation on this segment the modeled creep rate locally exceeds the long term rate reaching 8 mm yr 1 suggesting transient creep episodes the present study emphasizes the need for continuous monitoring of the surface velocity in the vicinity of major seismic gaps in terms of seismic hazard assessment

the displacement field of the landers earthquake mapped by radar interferometry geodetic data obtained by ground or space based techniques can be used to infer the distribution of slip on a fault that has ruptured in an earthquake although most geodetic techniques require a surveyed network to be in place before the earthquake1 3 satellite images when collected at regular intervals can capture co seismic displacements without advance knowledge of the earthquake s location synthetic aperture radar sar interferometry first introduced4 in 1974 for topographic mapping5 8 can also be used to detect changes in the ground surface by removing the signal from the topography9 10 here we use sar interferometry to capture the movements produced by the 1992 earthquake in landers california11 we construct an interferogram by combining topographic information with sar images obtained by the ers 1 satellite before and after the earthquake the observed changes in range from the ground surface to the satellite agree well with the slip measured in the field with the displacements measured by surveying and with the results of an elastic dislocation model as a geodetic tool the sar interferogram provides a denser spatial sampling 100 m per pixel than surveying methods1 3 and a better precision 3 cm than previous space imaging techniques12 13

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

the motion of crustal blocks driven by flow of the lower lithosphere and implications for slip rates of continental strike slip faults geodetic measurements in actively deforming areas of the continents reveal the pattern of deformation in the lithosphere if the dominant forces acting on crustal blocks are tractions at their bases then the long term motion of each block will be given by the average velocity of the underlying lithosphere slip rates between blocks estimated in this way from recent geodetic measurements across fault zones in the south island of new zealand and southern california are in good agreement with slip rates estimated geologically

the off fault deformation on the north anatolian fault zone and assessment of slip rate from carbonate veins abstract discrepancies in earthquake characteristics between geodetically determined and geologically observed slip rates are prominent for the eastern part of the north anatolian fault naf zone we investigated the spatial distribution of deformation across a naf fault segment near erzincan ne turkey through field observations of physiographic features along the principle fault and a carbonate filled fissure system we propose that distributed deformation across the principal fault in a stretching zone caused sigmoidal simple shear rotation of the carbonate filled fissure system displacement data of that sigmoidal fissure using isotopic age data of the carbonate fill indicates an off fault deformation rate of 15 90 3 19 mm per year for the holocene which corresponds to at least 50 of the geodetically estimated annual slip across the fault it is critical to take the distributed deformation into account for determining earthquake characteristics of active fault systems

modeling spatial fracture intensity as a control on flow in fractured rock spatial fracture intensity p 32 fracture area by volume is an important characteristic of a jointed rock mass although it can hardly ever be measured p 32 can be modeled based on available geological information such as spatial data of the fracture network flow in a mass composed of low permeability hard rock is controlled by joints and fractures in this article models were developed from a geological data set of fractured andesite in lanyu island taiwan where a site is investigated for possible disposal of low level and intermediate level radionuclide waste three different types of conceptual models of spatial fracture intensity distribution were generated an enhanced baecher s model ebm a levy lee fractal model llfm and a nearest neighborhood model nnm modeling was conducted on a 10 10 10 m synthetic fractured block simulated flow was forced by a 1 hydraulic gradient between two vertical x z faces of the cube from north to south with other boundaries set to no flow conditions resulting flow vectors are very sensitive to spatial fracture intensity p 32 flow velocity increases with higher fracture intensity p 32 r squared values of regression analysis for the variables velocity v v max and fracture intensity p 32 are 0 293 0 353 and 0 408 in linear fit and 0 028 0 08 and 0 084 in power fit higher r 2 values are positively linked with structural features but the relation between velocity and fracture intensity is non linear possible flow channels are identified by stream traces in the levy leefractal model

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

slip distribution along the north anatolian fault associated with the large earthquakes of the period 1939 to 1967 abstract between 1939 and 1967 six large fault ruptures formed a westward migrating sequence of events along a 900 km long nearly continuous portion of the north anatolian fault for these events the 1939 erzincan 1942 niksar erbaa 1943 tosya 1944 bolu gerede 1957 abant and 1967 mudurnu valley earthquakes i have compiled a record of dextral slip which contains nearly 100 points these data indicate that the amount of slip is irregularly distributed along the 1939 to 1967 rupture zone the maximum slip 7 5 m occurred during the 1939 earthquake in the eastern 150 km of the 900 km long rupture zone dextral offsets diminish very abruptly eastward but very gradually westward the rate of westward decrease in the 1939 to 1967 offsets is only slightly greater than the rate of westward decrease of post miocene displacement along the north anatolian fault this suggests that westward decrease in slip can be expected to be a general characteristic of earthquake ruptures along the north anatolian fault in the future nevertheless within the 1939 to 1967 slip distribution there are three regions that had less slip than the neighboring regions these i interpret as possible sites of large future earthquakes one of these seismic gaps has already experienced another earthquake in 1951 and subsequent creep

equivalent strike slip earthquake cycles in half space and lithosphere asthenosphere earth models by virtue of the images used in the dislocation solution the deformation at the free surface produced throughout the earthquake cycle by slippage on a long strike slip fault in an earth model consisting of an elastic plate lithosphere overlying a viscoelastic half space asthenosphere can be duplicated by prescribed slip on a vertical fault embedded in an elastic half space for the case in which each earthquake ruptures the entire lithosphere thickness h the half space equivalent slip rate is as follows depth interval 0 h slip identical to that in lithosphere asthenosphere model i e abrupt coseismic slip and no subsequent slip depth interval 2n 1 h to 2n 1 h n 1 2 slip rate uniform in space and dependent upon time as fn t exp t a where fn is a n 1 degree polynomial in t a is twice the asthenosphere relaxation time and t is measured from the instant after the preceding earthquake the slip rate averaged over the seismic cycle in each depth interval equals the secular rate of relative plate motion for reasonable values of a slip rates below 5h do not vary much from that mean value and can be treated as constant thus the surface deformation due to the earthquake cycle in the lithosphere asthenosphere model can be calculated very simply from the half space model with time dependent slip in the two depth intervals h 3h and 3h 5h and uniform slip at a rate equal to the secular relative plate velocity below depth 5h inversion of 1973 1988 geodetic measurements of deformation across the segment of the san andreas fault in the transverse ranges north of los angeles for the half space equivalent slip distribution suggests no significant slip on the fault above 30 km and a uniform slip rate of 36 mm yr below 30 km one equivalent lithosphere asthenosphere model would have a 30 km thick lithosphere and an asthenosphere relaxation time greater than 33 years but other models are possible

synthetic aperture radar interferometry to measure earth s surface topography and its deformation synthetic aperture radar interferometry insar from earth orbiting spacecraft provides a new tool to map global topography and deformation of the earth s surface radar images taken from slightly different viewing directions allow the construction of digital elevation models of meter scale accuracy these data sets aid in the analysis and interpretation of tectonic and volcanic landscapes if the earth s surface deformed between two radar image acquisitions a map of the surface dis placement with tens of meters resolution and subcentimeter accuracy can be con structed this review gives a basic overview of insar for earth scientists and presents a selection of geologic applications that demonstrate the unique capabilities of insar for mapping the topography and deformation of the earth

calibrating interferometric synthetic aperture radar insar images with regional gps network atmosphere models summary interferometric synthetic aperture radar insar technology provides a valuable tool for obtaining earth surface deformation and topography at high spatial resolution for crustal deformation studies similar to global positioning system gps insar measurements are affected by the earth s ionospheric and tropospheric layers as the electromagnetic signals significantly refract while propagating through the different layers while gps signals propagating through the neutral atmosphere are affected primarily by the distribution pressure and temperature of atmospheric gases including water vapour the propagation through the ionosphere is mainly affectedbythenumberoffreeelectronsalongthesignalpath here wepresenttheuseofdense regional gps networks for extracting tropospheric zenith delays and ionospheric total electron content tec maps in order to reduce the noise levels in insar images the results show significant reduction in the root mean square rms values when simultaneously combining the two corrections both at short time periods where no surface deformation is expected and at longer periods where imaging of localized subsidence and fault creep is enhanced

continental transforms a view from the alpine fault abstract continental transform faults are dominantly highly localized strike slip shear zones hundreds of kilometers long that accumulate tens to hundreds of kilometers of displacement from work on the alpine fault we pose the questions what is the deep structure of a continental transform and how does the displacement become localized we review research on the alpine fault and propose a model in which the fault partitions at depth into a steep zone extending into the mantle with largely fault parallel motion and a flat ductile decollement in the lower crust the fault localizes around two thirds of the plate motion within a 100 km wide zone of distributed deformation a review of other active continental fault systems suggests that variation between them may reflect their tectonic origins the nature of the crust in which they develop the presence of a significant oblique component of motion and the displacement rate all however have evidence for the development of a single principal fault zone that carries 50 of the total displacement and extends as a localized zone of shear into the upper mantle we review mechanisms of strain weakening and suggest that localization of a principal fault may be initiated in the seismogenic crust and through a series of positive feedbacks eventually extend through the lower crust into the upper mantle

geodetic observations of postseismic creep in the decade after the 1999 izmit earthquake turkey implications for a shallow slip deficit the relationship between aseismic slip and tectonic loading is important for understanding both the pattern of strain accumulation along a fault and its ability to generate large earthquakes we investigate the spatial distribution and temporal evolution of aseismic creep on the western north anatolian fault naf using time series analysis of envisat interferometric synthetic aperture radar insar data covering the full extent of the 1999 izmit and duzce earthquake ruptures and spanning 2002 2010 discontinuities in the line of sight velocity across the fault imply that fault creep reaches the earth s surface at an average fault parallel rate of 5 mm yr along an 80 km section of the naf by combining insar and published gps velocities we are able to extract the vertical and east west components of motion and show that the adapazari basin is subsiding at a rate of 6 mm yr vertical motions have biased previous estimates of creep in this region the displacement time series close to the fault is consistent with an afterslip model based on rate and state friction which predicts a rapid deceleration in fault creep rate after the izmit earthquake to a near steady state 5 mm yr after 5 years projecting our model 200 years into the future we find that the cumulative displacement of 1 1 3 m is insufficient to account for the shallow coseismic slip deficit observed in previous studies distributed off fault deformation in the shallow crust or transient episodes of faster slip are likely required to release some of the long term strain during the earthquake cycle

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

reconciling viscoelastic models of postseismic and interseismic deformation effects of viscous shear zones and finite length ruptures we have developed a suite of earthquake cycle models for strike slip faults to investigate how finite ruptures and lithosphere scale viscous shear zones affect interseismic deformation in particular we assess whether localized and stationary interseismic deformation and large scale rapidly decaying postseismic transients may be explained with models incorporating either or both of these features models incorporating viscous shear zones give more stationary interseismic deformation than layered half space maxwell viscoelastic models producing similar early postseismic deformation in the near field this tendency is accentuated when the effective viscosity per unit width of the shear zone increases either with depth or with interseismic time models with finite 200 km long ruptures produce time dependent deformation similar to that from models with infinitely long ruptures though with smaller magnitude and somewhat more localized surface velocity fluctuations models incorporating a stiff lithosphere a low viscosity mantle asthenosphere and a crust to lithosphere scale viscous shear zone can replicate postseismic and interseismic deformation typical of large strike slip earthquakes such models require depth dependent power law or transient rheologies for the viscous shear zone material in the crust as long as the effective viscosity per unit shear zone width is approximately 1015 pa s m at crustal depths during the postseismic interval below the moho the shear zone effective viscosity per unit width must be higher throughout the seismic cycle at least over a short depth interval for example if shear zone viscosity per unit width is 5 1016 pa s m in the mantle lithosphere a model with a 50 km thick lithosphere and asthenosphere effective viscosities of 1 to 5 1018 pa s can reproduce reference velocity profiles for both postseismic and later interseismic deformation

measurement of interseismic strain accumulation across the north anatolian fault by satellite radar interferometry in recent years interseismic crustal velocities and strains have been determined for a number of tectonically active areas through repeated measurements using the global positioning system the terrain in such areas is often remote and difficult and the density of gps measurements relatively sparse in principle satellite radar interferometry can be used to make millimetric precision measurements of surface displacement over large surface areas in practice the small crustal deformation signal is dominated over short time intervals by errors due to atmospheric topographic and orbital effects here we show that these effects can be overcome by stacking multiple interferograms after screening for atmospheric anomalies effectively creating a new interferogram that covers a longer time interval in this way we have isolated a 70 km wide region of crustal deformation across the eastern end of the north anatolian fault turkey the distribution of deformation is consistent with slip of 17 32 mm yr below 5 33 km on the extension of the surface fault at depth if the gps determined slip rate of 24 1 mm yr is accepted the locking depth is constrained to 18 6 km

evaluation of transport of radioactive contaminant in fractured rock fractures are one of the major travel paths for contaminant transport and dispersion within a rock mass in this paper a methodology that combines the orthogonal fractured network and the mass conservation theorem is developed to allow safety assessment in a fractured rock mass information equivalent to fracture geometric characteristics is used to establish the orthogonal fracture network the transport model is further implemented to investigate radioactive contaminant migrations in the fractured media in addition this paper employs the proposed methodology in accordance with the investigated data on a radioactive waste disposal site in lan yu taiwan results indicated that the attenuation of concentration was due to the radioactive decay and the sorption of fractures comparison of the transport simulations between 137cs and 60co shows that the environmental impact of 137cs is more obvious than that of 60co 137cs requires more time and more distance to reduce the concentration of the contaminant plume to below the maximum permissible activity

present day strain distribution across a segment of the central bend of the north anatolian fault zone from a persistent scatterers insar analysis of the ers and envisat archives

monitoring interseismic activity on the ilan plain ne taiwan using small baseline ps insar gps and leveling measurements partitioning from arc continent collision and backarc extension

application of fracture network model with crack permeability tensor on flow and transport in fractured rock article i nfo very efficient and practical models that use a crack tensor in terms of fracture geometry parameters for field rock mass are proposed a framework of the crack tensor approach combined with the identifying fracture controlled characteristics was successively used in a hypothetical hydro geological field for determining the dynamics of flow and transport of nuclides in fractured rock the modeling results and site specific data of lan yu island in taiwan reveal the characteristics associated with the geometry of fracture networks under the prescribed boundaries the simulated fracture networks were transformed into a grid based simulation procedure under consideration of the delaying effect in matrix owing to particle migration synthetic flow results indicate that the flowing vector is very sensitive to fracture intensity p32 high fracture intensity leads to high velocity heterogeneity generated by fractures field will result in the shortest and the longest transmission distance of 12 fold difference while taking into consideration of the ingrowth decay transport the plumes of the daughter product were relatively tortuous and anomalous to predict

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

correction to surface deformation due to a strike slip fault in an elastic gravitational layer overlying a viscoelastic gravitational half space by ting to yu john b rundle and jose fernandez in the paper surface deformation due to a strike slip fault in an elastic gravitational layer overlying aviscoelastic gravitational half space by ting to yu john b rundle and jos6 fern sndez journal of geophysical research 101 b2 3199 3214 1996 we have found errors while using the methods for calculating deformations produced by rectangular strike slip fault in the described viscoelastic gravitational layered earth

late caledonian dyke swarms of northern britain spatial and temporal intimacy between lamprophyric and granitic magmatism around the ross of mull pluton inner hebrides numerous sheet intrusions in iona and western mull scotland range from basic calc alkaline lamprophyres 45 sio2 mg number 65 75 thornton tuttle differentiation index di 45 through intermediate microdioritic compositions malchites to microgranitic porphyries 70 sio2 di 85 well defined major element trends remarkably consistent normalized multielement patterns and strong linear correlations between homologous elements e g sc v fe ti mg ni cr suggest that these sheets form one cogenetic suite although rb sr systems have been disturbed by archaean country rock contamination isotopic signatures and field evidence suggest that sheet emplacement bracketed that of the 414 ma ross of mull granitoid pluton a few late sheets cut the pluton exceptionally well developed multiple and composite sheets provide clear evidence forsimultaneity of lamprophyric and granitic magmas and an early group of sheets appears to be metamorphosed by the pluton the composite sheets comprise lamprophyric borders around thicker microgranitic quartz feldspar porphyries each rock type carrying xenocrysts or xenoliths of the other mutual contacts mark wide discontinuities in major element composition and sites of local hybridisation yet close trace element similarities strongly imply the lamprophyres and porphyries are cogenetic least squares modelling further shows that fractionation of actual phenocrysts in one lamprophyric margin could reproduce the porphyry with which it is in intimate contact the most primitive mg 70 lamprophyres combine very high k rb th and ba with morb like contents of cr ni ti y and sc implying a substantial mantle derived component and in some cases near primary status but not precluding some crustal participation biotite and hornblende lamprophyres are locally heteromorphic discontinuities in some chemical trends field evidence and initial87sr 86sr ratios 0 706 imply more complex interactions of fractionation and hybridisation with lamprophyric parent magmas plus crustal contamination in the genesis of the intermediate and microgranitic rocks the association between lamprophyric and granitic magmatism in iona western mull is not therefore fortuitous but genetic the lamprophyres have an undoubted role in the genesis not merely of the microgranitic minor intrusions but probably of the granite pluton itself

why are there no clockwise rotations along the north anatolian fault zone mechanical analysis of lithospheric deformation associated with an intracontinental transform fault suggests that it should be associated with a pattern of vertical axis rotations governed by the length of the fault its displacement and the mechanical properties of the lithosphere paleomagnetic results from two transects across the north anatolian fault zone indicate that there have been no measurable clockwise rotations on a regional scale associated with neogene dextral strike slip motion along the fault the data come from two groups of volcanic rocks that have been dated by k ar whole rock geochronology and have yielded ages ranging from 45 3 3 1 m y to 41 8 1 3 m y eocene and 21 7 1 5 m y to 7 5 0 2 m y miocene the absence of clockwise rotations suggests one of the following explanations 1 in contrast to the assumptions underlying current mechanical analyses deformation in the upper crust may not directly reflect the flow in the lithosphere and is localized along a single fault trace 2 deformation is distributed but rotations have been prevented by the presence of elongate fault blocks subparallel to the trace of the main fault 3 displacement along the north anatolian fault may be too small or the rheologic properties of the lithosphere such that rotations about vertical axes cannot be resolved instead of the predicted clockwise rotations measured declinations show that there has been approximately 30 of anticlockwise rotation of both the eocene and miocene volcanic rocks this suggests that no significant rotation occurred between eocene and miocene times and that the anticlockwise rotation occurred either immediately before or concurrent with the plio pleistocene slip on the north anatolian fault

the relationship between seismic and aseismic slip on the philippine fault on leyte island bayesian modeling of fault slip and geothermal subsidence delineating seismic and aseismic slip on faults allows the exploration of the complex relationship between these different modes of slip further quantifying them helps in the assessment of seismogenic potential we present a distributed slip model of the rate of aseismic slip along the leyte island section of the philippine fault and make comparisons with the extent of seismic slip from the latest significant earthquakes in july 2017 we derived both coseismic and aseismic slip distributions from kinematic inversions of synthetic aperture radar interferometric insar observations using a probabilistic bayesian framework velocity maps from stacking and time series analysis of alos interferograms spanning 4 years 2007 2011 show a step change in ground deformation right on the fault inverting for slip at depth reveals along strike variations of aseismic slip rate over 100 km aseismic slip on the surface reaches rates of more than 3 cm yr equivalent to the long term slip rate 3 3 0 2 cm yr over the same period a 20 km segment in tongonan appears to be locked this segment ruptured in mw 6 5 and mw 5 8 earthquakes on 6 and 10 july 2017 respectively as constrained by sentinel 1 and alos 2 insar data seismic slip appears to be restricted within the tongonan segment with up to 152 21 cm of left lateral displacement the slip budget and complementarity between the extents of interseismic and coseismic slip suggest that a seismogenic asperity exists in tongonan the presence of active hydrothermal systems and rate strengthening materials provide physical conditions that can promote aseismic slip

crustal deformation of the marlborough fault zone in the south island of new zealand geodetic constraints over the interval 1982 1994 crustal deformation across the marlborough fault zone in the south island of new zealand has been investigated by resurveying with the global positioning system gps a triangulation and trilateration network across part of the zone the principal strain rates which mostly have errors between 5 and 10 vary systematically across the region the principal axis of horizontal contraction gradually swings from se nw southeast of the hope fault to almost east west in the west the component of velocity parallel to the strike of the faults can account for 95 of the relative motion between the australian and pacific plates this component exhibits little variation west of the alpine wairau fault and shows an almost linear variation from the wairau fault to the east coast the horizontal strain rate northwest of the wairau fault is a small east west uniaxial contraction 0 1 ppm yr 1 if it is assumed that the crustal blocks bounded by the major faults are primarily driven by basal shear tractions with stresses on the faults being relatively insignificant then the slip rates expected on the faults may be derived from the observed velocity variations across the fault zone these predicted slip rates are in good agreement with geologically observed slip rates the larger slip rates on the hope fault are related to the observation that the strain field extends southeast of the fault over a distance large compared with the spacing between the faults the agreement between predicted and geological slip rates lends support to the idea that the observed pattern of deformation reflects that in the lower lithosphere

global positioning system constraints on plate kinematics and dynamics in the eastern mediterranean and caucasus we present and interpret global positioning system gps measurements of crustal motions for the period 1988 1997 at 189 sites extending east west from the caucasus mountains to the adriatic sea and north south from the southern edge of the eurasian plate to the northern edge of the african plate sites on the northern arabian platform move 18 2 mm yr at n25 5 w relative to eurasia less than the nuvel 1a circuit closure rate 25 1 mm yr at n21 7 w preliminary motion estimates 1994 1997 for stations located in egypt on the northeastern part of africa show northward motion at 5 6 2 mm yr also slower than nuvel ia estimates 10 1 mm yr at n2 4 e eastern turkey is characterized by distributed deformation while central turkey is characterized by coherent plate motion internal deformation of 2 mm yr involving westward displacement and counterclockwise rotation of the anatolian plate the anatolian plate is de coupled from eurasia along the right lateral strike slip north anatolian fault naf we derive a best fitting euler vector for anatolia eurasia motion of 30 7 0 8 n 32 6 0 4 e 1 2 0 1 myr the euler vector gives an upper bound for naf slip rate of 24 1 mm yr we determine a preliminary gps arabia anatolia euler vector of 32 9 1 2 n 40 3 1 1 e 0 8 0 2 myr and an upper bound on left lateral slip on the east anatolian fault eaf of 9 1 mm yr the central and southern aegean is characterized by coherent motion internal deformation of 2 mm yr toward the sw at 30 1 mm yr relative to eurasia stations in the se aegean deviate significantly from the overall motion of the southern aegean showing increasing velocities toward the trench and reaching 10 1 mm yr relative to the southern aegean as a whole

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

the erzincan earthquake of 13 march 1992 in eastern turkey the 13 march 1992 erzincan earthquake m 6 8 occurred in the eastern half of the erzincan basin the largest aftershock took place near pulumur on 15 march 1992 no clear surface breaks were observed although teleseismic studies suggested that it was a strike slip earthquake striking parallel to the north anatolian fault with a focus of approximately 10 2 km depth 30 km rupture length 95 cm of slip and a 1 16 1026 dyn cm seismic moment the aftershock distribution concentrated at an area of the intersection between the north anatolian fault and the ovacik fault these results indicate that the previously suggested seismic gap along the north anatolian fault east of erzincan still remains unruptured

geodetic determination of the kinematics of central greece with respect to europe implications for eastern mediterranean tectonics we use a new satellite laser ranging global positioning system slr gps solution at seven sites in anatolia and aegea to obtain a better definition of the extrusion motion of the anatolian aegean block with respect to europe we show that this motion can be described in a first approximation by a counterclockwise rotation which transfers most of the motion of arabia to anatolia we combine 78 displacement vectors obtained at common points of two triangulation nets measured in central greece in 1895 and 1975 with the slr gps measurements to compute the velocity field over greece with respect to europe these measurements indicate that central greece is a zone of extension between the anatolian aegean counterclockwise rotation to the south and the northern greece clockwise rotation to the north this extension is principally localized within the gulf of corinth to the east but is distributed to the west we then extrapolate this velocity field to the whole aegea and western anatolia using recently published gps results as well as the slr results the narrow dextral north anatolian fault which limits the velocity field to the north progressively gives way to a much wider boundary zone where extension becomes dominant we show that the collision between the mediterranean ridge and africa began 3 6 ma and we describe the modifications that this collision has produced on the kinematic pattern both in aegea and on the mediterranean ridge