fault interactions in the sea of marmara pull apart north anatolian fault earthquake clustering and propagating earthquake sequences summary knowledge on large earthquakes m 7 0 geology and fault kinematics is used to analyse conditions that favour isolated seismicity clustered earthquakes or propagating sequences along the north anatolian fault naf and the sea of marmara pull apart the overall naf marmara fault system is one of the most appropriate on earth to document fault interactions because reliable information covers almost completely two seismic cycles the past 500 yr coulomb stress analysis is used to characterize loading evolution in well identified fault segments including secular loading from below and lateral loading imposed by the occurrence of previous earthquakes earthquakes along the naf tend to occur where previous events have increased the stress but significant isolated events in the sea of marmara region 1894 1912 have occurred suggesting the secular loading has been the determining factor present day loading appears to be particularly high along the 70 km long segment located in the central marmara sea southwest of istanbul for the 18th century m 7 0 earthquake clusters we construct scenarios consistent with the tectonic and historical data we find that scenarios consistent with slip deficit and secular loading distributions from below clearly involve a sequence that propagates westward through the sea of marmara despite the structural complexity however the inference of a propagating sequence implies that each event has occurred in a segment previously stressed by lateral coulomb stress interactions the most likely scenarios for the propagating sequence are also consistent with coulomb stress interactions between faults with significant normal slip across the cinarcik basin propagating earthquake sequences do not occur every seismic cycle along the naf the loading has to be in a particular state of stress close to failure and uniform all along the fault segments to experience propagating earthquake sequences non uniform stress relief during the 18th century sequence explains the occurrence of isolated events in marmara in 1894 and 1912 as a consequence the wellknown 20th century sequence along the naf has not propagated as a sequence across the sea of marmara the most linear part of the naf across northern turkey behaves as a single fault segment accumulating stress during hundreds of years and rupturing entirely during very short periods the marmara pull apart fault system behaves as a major geometric complexity stopping or delaying the progression of earthquake clustering and propagating sequences fault zones interact with each other at a very large scale

block models of crustal motion in southern california constrained by gps measurements 1 we estimate slip rates on major active faults in southern california using a block model constrained by global positioning system measurements of interseismic deformation the block model includes the effects of block rotation and elastic strain accumulation consistent with a simple model of the earthquake cycle our estimates of the right lateral strike slip rate on the san andreas fault vary by at least a factor of 5 from a high of 35 9 0 5 mm yr in the carrizo plain to a low of 5 1 1 5 mm yr through the san bernadino segment shortening across the puente hills thrust and left lateral slip on the raymond hill fault are consistent with both thickening and escape tectonics in the los angeles basin discrepancies between geodetic and geologic slip rate estimates along the san andreas and san jacinto faults as well as in the eastern california shear zone may be explained by a temporal change in fault system behavior we find no substantial evidence for long term postseismic relaxation and infer that the viscosity of the lower crust upper mantle may be relatively high 1019 pa s

the 2013 mw 7 7 balochistan earthquake seismic potential of an accretionary wedge great earthquakes rarely occur within active accretionary prisms despite the intense long term deformation associated with the formation of these geologic structures this paucity of earthquakes is often attributed to partitioning of deformation across multiple structures as well as aseismic deformation within and at the base of the prism davis et al 1983 we use teleseismic data and satellite optical and radar imaging of the 2013 m w 7 7 earthquake that occurred on the southeastern edge of the makran plate boundary zone to study this unexpected earthquake we first compute a multiple point source solution from w phase waveforms to estimate fault geometry and rupture duration and timing we then derive the distribution of subsurface fault slip from geodetic coseismic offsets we sample for the slip posterior probability density function using a bayesian approach including a full description of the data covariance and accounting for errors in the elastic structure of the crust the rupture nucleated on a subvertical segment branching out of the chaman fault system and grew into a major earthquake along a 50 north dipping thrust fault with significant along strike curvature fault slip propagated at an average speed of 3 0 km s for about 180 km and is concentrated in the top 10 km with no displacement on the underlying decollement this earthquake does not exhibit significant slip deficit near the surface nor is there significant segmentation of the rupture we propose that complex interaction between the subduction accommodating the arabia eurasia convergence to the south and the ornach nal fault plate boundary between india and eurasia resulted in the significant strain gradient observed prior to this earthquake convergence in this region is accommodated both along the subduction megathrust and as internal deformation of the accretionary wedge

geological evidence for the last two faulting events on the north anatolian fault zone in the mudurnu valley western turkey abstract recurrence intervals of large earthquakes produced by slip on the north anatolian fault zone have been inferred mainly from historical records and included much ambiguity it is therefore important to date individual earthquakes using geological and geomorphological methods we excavated an exploratory trench across the surface rupture zone associated with the mudurnu valley western turkey earthquake of july 22 1967 and revealed evidence for the last two faulting events including the 1967 earthquake on the basis of the radiocarbon dates of sediments in the trench we estimate most conservatively that the penultimate event occurred after 1480 a d but believe it likely that the event occurred shortly after 1650 20 a d the penultimate event is likely to be correlated with the anatolian earthquake of august 17 1668 which caused severe damage along the north anatolian fault zone for about 600 km from bolu to erzincan

horizontal coseismic deformation of the 1999 chi chi earthquake measured from spot satellite images implications for the seismic cycle along the western foothills of central taiwan 1 the 1999 chi chi earthquake mw 7 6 broke a major thrust fault along the western foothills of the central range of taiwan we have measured the horizontal coseismic displacement field by correlating optical satellite images acquired before and after the earthquake these data reveal the fault trace and a clockwise rotation of surface displacements toward the north with much larger displacements and strain in the hanging wall at the surface coseismic slip increases from 5 6 m near the epicenter to 10 11 m to the north in the epicentral area we observe a left lateral strike slip zone trending n125 e and farther north a fault zone trending n s with a right lateral component the data were modeled using elastic dislocations the fault geometry consists of a shallow

seismic variability of subduction thrust faults insights from laboratory models 1 laboratory models are realized to investigate the role of interface roughness driving rate and pressure on friction dynamics the setup consists of a gelatin block driven at constant velocity over sand paper the interface roughness is quantified in terms of amplitude and wavelength of protrusions jointly expressed by a reference roughness parameter obtained by their product frictional behavior shows a systematic dependence on system parameters both stick slip and stable sliding occur depending on driving rate and interface roughness stress drop and frequency of slip episodes vary directly and inversely respectively with the reference roughness parameter reflecting the fundamental role for the amplitude of protrusions an increase in pressure tends to favor stick slip static friction is a steeply decreasing function of the reference roughness parameter the velocity strengthening weakening parameter in the state and rate dependent dynamic friction law becomes negative for specific values of the reference roughness parameter which are intermediate with respect to the explored range despite the simplifications of the adopted setup which does not address the problem of off fault fracturing a comparison of the experimental results with the depth distribution of seismic energy release along subduction thrust faults leads to the hypothesis that their behavior is primarily controlled by the depth and time dependent distribution of protrusions a rough subduction fault at shallow depths unable to produce significant seismicity because of low lithostatic pressure evolves into a moderately rough velocity weakening fault at intermediate depths the magnitude of events in this range is calibrated by the interplay between surface roughness and subduction rate at larger depths the roughness further decreases and stable sliding becomes gradually more predominant thus although interplate seismicity is ultimately controlled by tectonic parameters velocity of the plates trench and the thermal regime the direct control is exercised by the resulting frictional properties of the plate interface

the active main marmara fault abstract this paper presents selections from and a synthesis of a high resolution bathymetric sparker and deep towed seismic reflection data set recently acquired by the french ifremer r v le suroit in an e w deep trough that forms the northern half of the sea of marmara in nw turkey it includes the first high resolution complete bathymetric map of this area a single throughgoing dextral strike slip fault system which is the western continuation of the northern branch of the north anatolian fault cuts this trough lengthwise and joins the 1999 8 17 kocaeli earthquake fault with the 1912 8 09 arkoy murefte earthquake fault both of which display strike slip offset in its eastern one fourth the structure follows closely the northern margin of the deep trough whereas in the west it hugs its southern margin the eastern one fourth of the structure has a minor component of its displacement distributed across the deep trough owing to a possible original bend in the course of the dextral structure the present course of the north anatolian fault in the sea of marmara originated some 2 105 a ago by cutting across the older basin fabric generated by a dominant nne ssw extension before it began taking up major motion in the pliocene

breakdown in power law scaling in an analogue model of earthquake rupture and stick slip the gutenberg richter relationship for earthquakes in both global and fault specific seismicity has been interpreted as an indication of self similarity for events of all magnitudes it is observed however that the linear extrapolation of this relationship to high magnitude often fails to predict the frequency of occurrence of large events earthquake hazard cannot generally be predicted solely by examination of microseismicity it has been postulated that this break in power law scaling is due to a change in the rupture mechanism as the size of the rupture becomes larger than the thickness of the seismogenic layer here we examine the scaling of earthquake like failures in an analogue model consisting of an elastic solid confined and driven by a steel stressing apparatus sensors embedded in the solid close to the interface with a rough substrate measure local strain and allow the examination of rupture scaling and moment release the frequency magnitude distribution for events is a power law at low magnitude but exhibits enhancement for high energy events which is consistent with a characteristic earthquake model we show unambiguously that this break in scaling corresponds to ruptures whose diameter is equal to the smaller dimension of the model fault

postseismic stress transfer explains time clustering of large earthquakes in mongolia abstract three m 8 earthquakes have occurred in mongolia during a 52 year period 1905 1957 since these earthquakes were well separated in space 400 km the coseismic stress change is far too low 0 001 bar to explain mutual earthquake triggering by contrast postseismic relaxation gradually causes a significant stress change 0 1 0 9 bar over large distances using a spring slider model to simulate earthquake interaction we find that viscoelastic stress transfer may be responsible for the earthquake time clustering observed in active tectonic areas therefore revealing postseismic strain by satellite geodesy and modeling earthquake clusters could improve our understanding of earthquake occurrence especially in zones where large earthquakes have already struck in past decades

a dislocation model for the earthquake cycle at convergent plate boundaries summary a kinematic model for the earthquake cycle at convergent plate boundaries has been constructed on the basis of dislocation theory we model the lithosphere asthenosphere system by a stratified semi infinite medium under gravity consisting of an elastic surface layer an intervening layer with maxwell viscoelasticity and an elastic substratum the steady motion of plate convergence is naturally represented by uniform slip at a constant rate on the upper boundary of the descending oceanic plate the occurrence of interplate earthquakes is regarded as a perturbation of the steady state plate motion and represented by a periodic sequence of step slips on a finite seismic zone of the plate interface based on dislocation theory we can show that the steady slip on the interface deeper than the lithosphereasthenosphere boundary does not contribute to surface deformation in the steady state accordingly the surface deformation associated with the earthquake cycle is given by the superposition of viscoelastic responses to the steady slip on the interface shallower than the lithosphere asthenosphere boundary the steady backslip on the seismic zone and the periodic sequence of seismic slips we have computed cyclic patterns of vertical displacements at the free surface for two representative cases taking account of gravity effects the patterns of vertical displacements differ notably in the latter stage of the cycle depending on whether or not the seismic zone extends through the entire thickness of the lithosphere after the completion of one earthquake cycle our model yields a certain amount of permanent deformation resulting from the steady plate convergence the pattern of the permanent deformation that is characterized by steep uplift on the continental side sharp subsidence at the plate boundary and gentle uplift on the oceanic side is irrespective of the cyclic process of stress accumulation and release repeated on the seismic zone the earthquake cycle model developed here provides a possible explanation for the formation of earthquake origin marine terraces

postseismic motion after the 2001 mw 7 8 kokoxili earthquake in tibet observed by insar time series 1 on november 14th 2001 a mw 7 8 earthquake occurred in the kokoxili region of northern tibet the earthquake ruptured more than 400 km along the western part of the kunlun fault with a maximum of 8 m left lateral slip in this paper we use a multitemporal interferometric sar insar time series technique to map the postseismic motion following the large kokoxili event sar data from envisat descending orbits along five adjacent tracks covering almost the entire ruptured fault length are used to calculate the displacement time series for a period between 2 and 6 years after the earthquake a peak to trough signal of 8 cm in the radar line of sight is observed during the period between 2003 and 2008 two different mechanisms are employed to explain the observed surface displacements namely afterslip and viscoelastic relaxation the observations inverted for afterslip on and below the coseismic rupture plane shows that the maximum slip in the afterslip model is 0 6 m the position of the maximum postseismic slip is located in the middle of two relatively high coseismic slip patches which suggests that afterslip is a plausible mechanism models of viscoelastic stress relaxation in a maxwell half space give a best fitting viscosity for the mid to lower crust of 2 5 1019 pa s and the principal postseismic relaxation process is due to viscous flow in the lower crust to upper mantle however the insar observations are incapable of distinguishing between localized afterslip and distributed viscoelastic relaxation deformation and the lowest misfits are produced by mixed models of viscoelastic relaxation in the mantle below 70 km and afterslip in the crust modeling of viscoelastic relaxation in a maxwell half space and also a mixed mechanism model enables us to place an effective viscosity of 2 1019 pa s on the lower crust to mantle of northern tibet

self healing pulse like shear ruptures in the laboratory models predict that dynamic shear ruptures during earthquake faulting occur as either sliding cracks where a large section of the interface slides behind a fast moving rupture front or self healing slip pulses where the fault relocks shortly behind the rupture front we report experimental visualizations of crack like pulse like and mixed rupture modes propagating along frictionally held incoherent interfaces separating identical solids and we describe the conditions under which those modes develop a combination of simultaneously performed measurements via dynamic photoelasticity and laser interferometry reveals the rupture mode type the exact point of rupture initiation the sliding velocity history and the rupture propagation speed

crustal deformation in great california earthquake cycles periodic crustal deformation associated with repeated strike slip earthquakes is computed for the following model a depth l h extending downward from the earth s surface at a transform boundary between uniform elastic lithospheric plates of thickness h is locked between earthquakes it slips an amount consistent with remote plate velocity vpl after each lapse of earthquake cycle time tcy lower portions of the fault zone at the boundary slip continuously so as to maintain constant resistive shear stress the plates are coupled at their base to a maxwellian viscoelastic asthenosphere through which steady deep seated mantle motions compatible with plate velocity are transmitted to the surface plates the coupling is described approximately through a generalized elsasser model we argue that the model gives a more realistic physical description of tectonic loading including the time dependence of deep slip and crustal stress buildup throughout the earthquake cycle than do simpler kinematic models in which loading is represented as imposed uniform dislocation slip on the fault below the locked zone parameters of the model are chosen to fit seismic and geologic constraints and the apparent time dependence of surface strain rates along presently locked traces of the 1857 and 1906 san andreas ruptures we fix vpl 35 mm yr tcy 160 years and l 9 11 km based on earthquake nucleation depths the geodetic data are then found to be described reasonably within the context of a model that is locally uniform along strike and symmetric about a single san andreas fault strand by lithosphere thickness h 20 30 km and elsasser relaxation time tr 10 16 years we therefore suggest that the asthenosphere appropriate to describe crustal deformation on the earthquake cycle time scale lies in the lower crust and perhaps crust mantle transition zone and has an effective viscosity between about 2 1018 and 1019 pa s depending on the thickness assigned to the asthenospheric layer predictions based on the chosen set of parameters are also consistent with data on variations of contemporary surface strain and displacement rates as a function of distance from the 1857 and 1906 rupture traces although the fit is degraded by asymmetry relative to the fault and by slip on adjacent fault strands

a physical model of the effect of a shallow weak layer on strong ground motion for strike slip ruptures we report results of foam rubber modeling of the effect of a shallow weak layer on ground motion from strike slip ruptures computer modeling of strong ground motion from strike slip earthquakes has involved somewhat arbitrary assumptions about the nature of slip along the shallow part of the fault e g fixing the slip to be zero along the upper 2 kilometers of the fault plane in order to match certain strong motion accelerograms most modeling studies of earthquake strong ground motion have used what is termed kinematic dislocation modeling in kinematic modeling the time function for slip on the fault is prescribed and the response of the layered medium is calculated unfortunately there is no guarantee that the model and the prescribed slip are physically reasonable unless the true nature of the medium and its motions are known ahead of time there is good reason to believe that in many cases faults are weak along the upper few kilometers of the fault zone and may not be able to maintain high levels of shear strain required for high dynamic energy release during earthquakes physical models of faulting as distinct from numerical or mathematical models are guaranteed to obey static and dynamic mechanical laws foam rubber modeling studies have been reported in a number of publications the object of this paper is to present results of physical modeling using a shallow weak layer in order to verify the physical basis for assuming a long rise time and a reduced high frequency pulse for the slip on the shallow part of faults it appears a 2 kilometer deep weak zone along strike slip faults could indeed reduce the high frequency energy radiated from shallow slip and that this effect can best be represented by superimposing a small amplitude short rise time pulse at the onset of a much longer rise time slip a weak zone was modeled by inserting weak plastic layers of a few inches in thickness into the foam rubber model for the 15 cm weak zone the average pulse is reduced by a factor of 0 46 the factor for the 20 cm case reduction is 0 11 for the 30 cm case it is 0 045 from these results we can see that the thicker the weak layer the more difficult it is for a short rise time acceleration pulse to push its way through the weak layer to the surface this is thus an approximate justification for reducing the high frequency radiation from shallower parts of strike slip faults if it is known that the shallow part of the fault is weak or has not stored up shear stress

thin plate modeling of interseismic deformation and asymmetry across the altyn tagh fault zone 1 creeping dislocations in an elastic half space are commonly used to model interseismic deformation however this semi infinite conventional model can lead to biased inferences of the slip rate and the fault locking depth we therefore favor the use of a thin plate model sheared at its base in agreement with the current knowledge of the seismogenic thickness of the lithosphere we apply these two approaches to interpret insar data across the altyn tagh fault system near longitude 94 e in a region of clear asymmetry of interseismic velocity distribution our results suggest a locking depth of 7 9 km and a present day geodetic slip rate of 8 10 mm yr consistent with the rate based on gps measurements at this longitude we interpret the asymmetric pattern as the joint effect of a rigidity decrease from the tarim basin north of the fault to the qaidam basin south of it and of a 5 7 km southward offset of strain concentration from the quaternary northern active fault this suggests that the rigidity contrast as well as most of the creep at depth occurs on the southern geological fault

a dislocation model of strain accumulation and release at a subduction zone strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface this process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip because steady state subduction does not contribute to the deformation at the free surface deformation observed there is completely equivalent to that produced by the slip cycle alone the response to that slip is simply the response of a particular earth model to embedded dislocations for a purely elastic earth model the deformation cycle consists of a coseismic offset followed by a linear in time recovery to the initial value during the interval between earthquakes for an elastic viscoelastic earth model elastic lithosphere over a viscoelastic asthenosphere the postearthquake recovery is not linear in time records of local uplift as a function of time indicate that the long term postseimic recovery is approximately linear suggesting that elastic earth models are adequate to describe the deformation cycle however the deformation predicted for amore simple elastic half space earth does not reproduce the deformation observed along the subduction zones in japna at all well if stick slip is restricted to the main thrust zone as recognized earlier by shimazaki seno and kato the uplift profiles could be explained if stick slips were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km but independent evidence suggests that stick slip at such depths if unlikely less

earthquake nucleation on faults with rate and state dependent strength abstract dieterich j h 1992 earthquake nucleation on faults with rate and state dependent strength in t mikumo k aki m ohnaka l j ruff and p k p spudich editors earthquake source physics and earthquake precursors tectonophysics 211 115 134 faults with rate and state dependent constitutive properties reproduce a range of observed fault slip phenomena including spontaneous nucleation of slip instabilities at stresses above some critical stress level and recovery of strength following slip instability calculations with a plane strain fault model with spatially varying properties demonstrate that accelerating slip precedes instability and becomes localized to a fault patch the dimensions of the fault patch follow scaling relations for the minimum critical length for unstable fault slip the critical length is a function of normal stress loading conditions and constitutive parameters which include dc the characteristic slip distance if slip starts on a patch that exceeds the critical size the length of the rapidly accelerating zone tends to shrink to the characteristic size as the time of instability approaches solutions have been obtained for a uniform fixed patch model that are in good agreement with results from the plane strain model over a wide range of conditions above the steady state stress the logarithm of the time to instability linearly decreases as the initial stress increases because nucleation patch length and premonitory displacement are proportional to dc the moment of premonitory slip scales by d3c the scaling of dc is currently an open question unless dc for earthquake faults is significantly greater than that observed on laboratory faults premonitory strain arising from the nucleation process for earthquakes may by too small to detect using current observation methods excluding the possibility that dc in the nucleation zone controls the magnitude of the subsequent earthquake then the source dimensions of the smallest earthquakes in a region provide an upper limit for the size of the nucleation patch

stress evolution before and after the 2008 wenchuan china earthquake abstract we examine the distribution of coulomb stress before and following the m 7 9 wenchuan earthquake and model the time dependent components of the stress field into the future our modelling shows that the combined coseismic plus postseismic viscoelastic stress change due to the preceding large events m 6 7 at the hypocentral area of the wenchuan earthquake is negative which may have led to a delay on its occurrence between 16 and 95 yr depending on the rake of the rupture on the other hand coseismic stress modelling of the wenchuan earthquake shows strong interaction particularly with the faults close to its sw and ne ends such interactions have previously been shown to trigger further large earthquakes in addition redistribution of stresses due to continuing postseismic deformation will potentially increase the stress levels not only on these faults but also on the distant faults such as the xianshuihe kunlun and tazang faults postseismic viscoelastic stress load on these distant faults over the next 22 yr since the earthquake could be up to 2 5 times larger than that of the coseismic stress increase alone this would potentially advance these faults in their seismic cycle and thereby encourage future earthquakes resulting in clustering of large earthquakes in the region this behaviour however strongly depends on the stage that these faults are in their seismic cycle which is not known currently given the high density of population in the area with several cities chengdu guangyuan and hanzhong of more than 1 million inhabitants our results indicate a need for detailed paleoseismic studies to help elucidate the state of the stressed faults

evidence for postseismic deformation of the lower crust following the 2004 mw6 0 parkfield earthquake previous studies have shown that postseismic relaxation following the 2004 mw6 0 parkfield ca earthquake is dominated by afterslip however we show that some fraction of the afterslip inferred from kinematic inversion to have occurred immediately below the seismically ruptured area may in fact be a substitute for viscous postseismic deformation of the lower crust using continuous gps and synthetic aperture radar interferometry we estimate the relative contribution of shallow afterslip at depth less than 20km and deeper seated deformation required to account for observed postseismic surface displacements exploiting the possible separation in space and time of the time series of displacements predicted from viscoelastic relaxation we devise a linear inversion scheme that allows inverting jointly for the contribution of afterslip and viscoelastic flow as a function of time we find that a wide range of models involving variable amounts of viscoelastic deformation can fit the observations equally well provided that they allow some fraction of deep seated deformation at depth larger than 20 km these models require that the moment released by postseismic relaxation over 5 years following the earthquake reached nearly as much as 200 of the coseismic moment all the models show a remarkable complementarity of coseismic and shallow afterslip distributions some significant deformation at lower crustal depth 20 26 km is required to fit the geodetic data the condition that postseismic deformation cannot exceed complete relaxation places a constraint on the amount of deep seated deformation the analysis requires an effective viscosity of at least 10 18 pa s of the lower crust assuming a semi infinite homogeneous viscous domain this deep seated deformation is consistent with the depth range of tremors which also show a transient postseismic response and could explain as much as 50 of the total postseismic geodetic moment the remaining fraction being due to afterslip at depth shallower than 20 km lower crustal postseismic deformation could reflect a combination of localized ductile deformation and aseismic frictional sliding

a viscoelastic coupling model for the cyclic deformation due to periodically repeated earthquakes at subduction zones vertical displacements due to periodic reverse faulting events in an elastic plate overlying a viscoelastic maxwell half space are obtained and compared with the observed deformation cycle coseismic strain release postseismic transients interseismic strain accumulation from japan the viscoelastic effects including the influence of buoyant restoring forces are obtained using the method developed by rundle and plate convergence and strain accumulation are incorporated following the procedure suggested by savage the resulting deformation cycle is compared with that of an analogous elastic half space dislocation model in which postearthquake effects are due to transient aseismic slip below the coseismic fault cyclic deformation is similar but not identical for the two models and observations from southwest japan suggest the superiority of the viscoelastic coupling model in particular inclusion of the effects of steady state flow in the asthenosphere overcomes a defect of the elastic half space model and results in agreement with the observed interseismic movement pattern several aspects of the postseismic deformation its landward migration and its transition to the interseismic phase of the cycle are explained as well but the short duration of near trench transients relative to those observed farther inland is not matched the success of a buried slip model in explaining early postseismic near trench movements and asthenospheric flow in accounting for cumulative postearthquake transient motions suggests the existence of a transition zone between lithosphere and asthenosphere whose behavior is brittle elastic in the short term and ductile for longer term deformation and such a modification may reconcile remaining discordant observations however reasonable variations in coupling model parameters cannot account for observed differences in the deformation cycle in other parts of japan and these regional differences remain unexplained

optimal fault resolution in geodetic inversion of coseismic data summary with the continued growth in availability of differential interferometry from synthetic aperture radar and gps data space based geodesy has been widely applied to image the coseismic displacement field and to retrieve the static dislocation over the fault plane for almost all the significant earthquakes of the past two decades this is performed by linear data inversion over a set of subfaults generally characterized by a constant and predefined or manually adjusted dimensions in this paper we propose a new algorithm to automatically retrieve an optimized fault subdivision in the linear inversion of coseismic geodetic data the code iteratively keeps the parameter resolution close to a predefined high value we first discuss the rationale supporting ouralgorithmand afteradetaileddescriptionofitsimplementation weanalysetheadvantages of its introduction in the data inversion the algorithm was tested against an exhaustive range of synthetic and real data sets and fault mechanisms among them we present the results for themw 6 2 2009 l aquila central italy earthquake and compare the new and previously published slip distributions showing the disappearance of misleading slip pattern and the increased resolution for shallower zones

photo acoustic study of subshear and supershear ruptures in the laboratory abstract we report supershear ruptures observed in spontaneously nucleating laboratory microearthquakes and describe the signature of the associated mach wavefront radiation transducers detect the wavefield both close and at a distance from the fault the rupture velocities are inferred from either photoelastic high speed imaging or using the acoustic recordings both methods yield compatible estimates of sufficient accuracy to discriminate between sub shear and supershear ruptures the acoustic records allow to characterize the mach wavefront radiated from the supershear rupture front in particular its amplitude and decay velocity functions recovered by integrating high frequency accelerometer signals in the case of supershear ruptures consist in a double pulse function a first pulse traveling at supershear velocity followed by a second pulse traveling at the rayleigh wave velocity conversely the sub shear event is characterized by a single pulse finally we perform numerical simulations of our experiment using a prescribed supershear rupture velocity the synthetic waveforms obtained from these simulations including the mach wave amplitude and phase yield a satisfactory fit to the experimental results

constitutive relations between dynamic physical parameters near a tip of the propagating slip zone during stick slip shear failure abstract constitutive relations between physical parameters in the cohesive zone during stick slip shear failure are experimentally investigated stick slip was generated along a 40 cm long precut fault in tsukuba granite samples using a servocontrolled biaxial loading apparatus dynamic behavior during local breakdown processes near a tip of the slipping zone is revealed the slip velocity and acceleration are given as a function of the slip displacement and the cohesive or breakdown shear stress as a function of the slip velocity a cycle of the breakdown and restrengthening process of stick slip is composed of five phases characterized in terms of the cohesive strength and the slip velocity the cohesive strength can degrade regardless of the slip velocity during slip instabilities the maximum slip acceleration u max and the maximum slip velocity u max are obtained experimentally as u max 2 u c u max 2 and u max b g v where u c is the critical displacement b the breakdown stress drop g the rigidity and v the rupture velocity these relations are consistent with ida s theoretical estimation based on the cohesive zone model the above formula gives good estimates for the maximum slip acceleration of actual earthquakes the cutoff frequency max of the power spectral density of the slip acceleration increases with increasing normal stress in particular max is found to be directly proportional to the normal stress n within the normal stress range less than 17 mpa as max khz 4 0 n mpa n max increase with an increase in u max or u max all these results lead to the conclusion that u max u max and max increase with increasing normal stress this is consistent with a previous observation that b increases with increasing normal stress the above empirical linear relation between max and n can be explained by a linear dependence of b on n the size scale dependence of physical parameters is discussed and such size scale dependent parameters as u max and apparent fracture energy are scaled to the breakdown zone size regarded as a characteristic length to extend the results obtained in the laboratory to an earthquake failure in the earth

new empirical relationships among magnitude rupture length rupture width rupture area and surface displacement source parameters for historical earthquakes worldwide are compiled to develop a series of empirical relationships among moment magnitude m surface rupture length subsurface rupture length downdip rupture width rupture area and maximum and average displacement per event the resulting data base is a significant update of previous compilations and includes the additional source parameters of seismic moment moment magnitude subsurface rupture length downdip rupture width and average surface displacement each source parameter is classified as reliable or unreliable based on our evaluation of the accuracy of individual values only the reliable source parameters are used in the final analyses in comparing source parameters we note the following trends 1 generally the length of rupture at the surface is equal to 75 of the subsurface rupture length however the ratio of surface rupture length to subsurface rupture length increases with magnitude 2 the average surface displacement per event is about one half the maximum surface displacement per event and 3 the average subsurface displacement on the fault plane is less than the maximum surface displacement but more than the average surface displacement thus for most earthquakes in this data base slip on the fault plane at seismogenic depths is manifested by similar displacements at the surface log linear regressions between earthquake magnitude and surface rupture length subsurface rupture length and rupture area are especially well correlated showing standard deviations of 0 25 to 0 35 magnitude units most relationships are not statistically different at a 95 significance level as a function of the style of faulting thus we consider the regressions for all slip types to be appropriate for most applications regressions between magnitude and displacement magnitude and rupture width and between displacement and rupture length are less well correlated and have larger standard deviation than regressions between magnitude and length or area the large number of data points in most of these regressions and their statistical stability suggest that they are unlikely to change significantly in response to additional data separating the data according to extensional and compressional tectonic environments neither provides statistically different results nor improves the statistical significance of the regressions regressions for cases in which earthquake magnitude is either the independent or the dependent parameter can be used to estimate maximum earthquake magnitudes both for surface faults and for subsurface seismic sources such as blind faults and to estimate the expected surface displacement along a fault for a given size earthquake

probing the mechanical properties of seismically active crust with space geodesy study of the coseismic deformation due to the 1992 mw7 3 landers southern california earthquake 1 the coseismic deformation due to the 1992 m w 7 3 landers earthquake southern california is investigated using synthetic aperture radar sar and global positioning system gps measurements the ers 1 satellite data from the ascending and descending orbits are used to generate contiguous maps of three orthogonal components east north up of the coseismic surface displacement field the coseismic displacement field exhibits symmetries with respect to the rupture plane that are suggestive of a linear relationship between stress and strain in the crust interferometric synthetic aperture radar insar data show small scale deformation on nearby faults of the eastern california shear zone some of these faults in particular the calico rodman and pinto mountain faults were also subsequently strained by the 1999 m w 7 1 hector mine earthquake i test the hypothesis that the anomalous fault strain represents essentially an elastic response of kilometer scale compliant fault zones to stressing by nearby earthquakes fialko et al 2002 the coseismic stress perturbations due to the landers earthquake are computed using a slip model derived from inversions of the insar and gps data calculations are performed for both homogeneous and transversely isotropic half space models the compliant zone model that best explains the deformation on the calico and pinto mountain faults due to the hector mine earthquake successfully predicts the coseismic displacements on these faults induced by the landers earthquake deformation on the calico and pinto mountain faults implies about a factor of 2 reduction in the effective shear modulus within the similar to 2 km wide fault zones the depth extent of the low rigidity zones is poorly constrained but is likely in excess of a few kilometers the same type of structure is able to explain high gradients in the radar line of sight displacements observed on other faults adjacent to the landers rupture in particular the lenwood fault north of the soggy lake has likely experienced a few centimeters of left lateral motion across 1 km wide compliant fault zone having the rigidity reduction of more than a factor of 2 the inferred compliant fault zones are interpreted to be a result of extensive damage due to past earthquakes

controls of shear zone rheology and tectonic loading on postseismic creep 1 postseismic deformation is well documented in geodetic data collected in the aftermath of large earthquakes in the postseismic time interval gps is most sensitive to creep in the lower crust or upper mantle activated by earthquake generated stress perturbations in these regions deformation may be localized on an aseismic frictional surface or on a ductile shear zone these two hypotheses imply specific rheologies and therefore time dependence of postseismic creep hence postseismic creep constitutes a potential probe into the rheology of aseismic regions of the lithosphere i present a simple shear zone model of postseismic creep in which the rheology of the creeping element can be varied in the absence of tectonic loading during the postseismic time interval the displacement history of the shear zone obeying a power law rheology with stress exponent n follows an analytical relaxation curve parameterized by 1 n for a frictional surface postseismic creep follows the same relaxation law in the limit 1 n 0 a rough estimate of the apparent stress exponent can be obtained from continuous gps records application to data collected after the 1994 sanriku earthquake yields 1 n 0 1 which is consistent with dislocation creep mechanisms however the records of two other subduction zone events the 2001 peru event and the 1997 kronotski earthquake and a continental strike slip earthquake the 1999 izmit earthquake require negative 1 n rather than characterizing the shear zone rheology these negative exponents indicate that reloading of the shear zone by tectonic forces is important numerical simulations of postseismic deformation with nonnegligible reloading produce curves that are well fit by the generalized relaxation law with 1 n 0 although the actual stress exponent of the rheology is positive while this prevents rheology from being tightly constrained by the studied gps records it indicates that reloading is important in the postseismic time interval in other words the stress perturbation induced by an earthquake is comparable to the stress supported by ductile shear zones in the interseismic period

interaction of the san andreas fault creeping segment with adjacent great rupture zones and earthquake recurrence at parkfield three dimensional finite element calculations are employed to study interactions in space and time between the creeping segment of the san andreas fault in central california and the adjacent currently locked zones of the 1857 and 1906 great earthquakes vertically the model consists of an elastic upper crust over a maxwell viscoelastic region representing the entire lower crust or a narrower horizontal detachment layer and a stiffer and more viscous upper mantle the crust has a single vertical fault extending to the top of the mantle at 25 km depth in zones along strike corresponding to the 1857 and 1906 events the top 12 5 km of the fault is locked against slip except in great earthquakes below the locked zones and everywhere along the creeping region between them the fault is freely slipping the model parameters are compatible with seismological and geological observations and with a ratio of maxwell relaxation time to the relaxing layer thickness in the range 1 to 2 yr km as established by li and rice 1987 and fares and rice 1988 based on fits to geodetic data along the san andreas fault an imposed constant far field shear motion and periodic 1857 and 1906 type earthquakes generate slip rates along the creeping fault segment that evolve in time throughout the entire earthquake cycle shortly after an adjacent great earthquake slip rates in the creeping zone are higher than the far field velocity while later in the cycle they are lower hence time dependency should be accounted for when measurements of fault slip are used to estimate the plate motion if parkfield earthquakes are a response to a time dependent loading of the type simulated here their recurrence interval would tend to lengthen with time since the 1857 event thus the hypothesis of characteristic periodic earthquakes at parkfield may not provide the best estimate of the occurrence time of the next event using for example the statistics of past events and assuming that parkfield earthquakes are a response to a slip deficit near middle mountain and that the elastic crustal layer is 17 5 km thick we find that the next event is predicted for about 1992 9 years if the lower crust is a 7 5 km thick layer having a material relaxation time of 15 years and 1995 11 years if the 7 5 km thick lower crust is characterized by a relaxation time of 7 5 years these values may be compared to the 1988 7 years estimate based on periodicity in time the modeling results also indicate that the interaction between the 1857 and 1906 rupture zones is small

earthquake scaling fault segmentation and structural maturity abstract slip and length measurements on earthquakes suggest large stress drop variability we analyze an extended set of slip length measurements for large earthquakes m 6 to seek for the possible origin s of this apparent variability we propose that such variability arises from earthquakes breaking a variable number of major fault segments that number depends on the strength of the inter segment zones which itself depends on the structural maturity of the faults we propose new d max l parameterizations based on that idea of multiple segment ruptures in such parameterizations each broken segment roughly scales as a crack while the total multi segment rupture does not stress drop on individual segments is roughly constant only varying between 3 5 to 9 mpa the slight variation that is still observed depends on fault structural maturity more mature faults have lower stress drops than immature ones the new d max l functions that we propose reduce uncertainties with respect to available relationships they thus provide a more solid basis to estimate seismic hazard by integrating fault properties revealed by geological studies

identifying the unique ground motion signatures of supershear earthquakes theory and experiments the near field ground motion signatures associated with sub rayleigh and supershear ruptures are investigated using the laboratory earthquake experiment originally developed by rosakis and co workers xia et al 2004 2005a lu et al 2007 rosakis et al 2007 heterodyne laser interferometers enable continuous high bandwidth measurements of fault normal fn fault parallel fp and vertical v particle velocity ground motion records at discrete locations on the surface of a homalite 100 test specimen as a sub rayleigh or a supershear rupture sweeps along the frictional fault photoelastic interference fringes acquired using high speed digital photography provide a synchronized spatially resolved whole field view of the advancing rupture tip and surrounding maximum shear stress field the first phase of experimental investigations examine and verify the ground motion signatures of supershear ruptures experimental results demonstrate that a shear mach front produced by a stable supershear rupture is characterized by a dominant fp velocity component the situation is shown to reverse in the sub rayleigh rupture speed regime whereby the fn particle velocity component dominates the ground motion record additional distinguishing particle velocity signatures consistent with theoretical and numerical predictions and repeatedly observed in experimental records are 1 a pronounced peak in the fp velocity record induced by the leading dilatational field which sweeps the measurement station in advance of the shear mach front and 2 a pronounced velocity swing in the fn record associated with the arrival of a trailing rayleigh sub rayleigh secondary rupture which follows the arrival of the shear mach front analysis of the particle velocity records also confirms 2d steady state theoretical predictions pertaining to the separation attenuation and radiation partitioning of the shear and dilatational portions of the rupture velocity field components the second phase of our experimental investigations re examine the 2002 mw7 9 denali fault earthquake and the remarkable set of near source ground motion records obtained at ps10 located approximately 85 km east of the epicenter and just 3 km north of the fault along the alaska pipeline motivated by the analysis and interpretation of these records by ellsworth et al 2004 dunham and archuleta 2004 2005 we attempt to mimic the denali strike slip rupture scenario and replicate the ps10 ground motion signatures using a laboratory earthquake experiment the experiments feature a left to right west to east propagating right lateral rupture within a homalite 100 test specimen with particle velocity data collected at a near field station situated just above north of the fault both sub rayleigh and supershear laboratory earthquake experiments are conducted using the denali ps10 configuration in order to compare and contrast the resulting particle velocity signatures supershear laboratory records capture all of the prominent features displayed within the ps10 ground motion records noted velocity signatures are correlated to the location of the rupture fronts and their noted arrival times in the synchronized photoelastic image sequence scaling relationships are also presented which transform the laboratory records through six orders of magnitude in time to match the scale of the ps10 ground motion records the strong correlation between the scaled experimental records and the actual ps10 ground motion records support the hypothesis that the denali strike slip fault exhibited a supershear burst finally we present a 2d steady state stress velocity formulation that relates the fp and fn particle velocity records measured close to the fault to the evolution of the stress tensor at the same location a locally steady state condition is assumed within a restricted time interval in order to invoke these relationships and estimate the dynamic stresses xx t and t at the near fault station dynamic stress measurements enable a new class of friction investigations using the laboratory earthquake configuration experimental findings are presented which capture the temporal and spatial distributions of xx and evolution of the dynamic friction coefficient and velocity weakening behavior of a supershear slip pulse

the seismic cycle at subduction thrusts insights from seismo thermo mechanical models 1 the underestimation of the size of recent megathrust earthquakes illustrates our limited understanding of their spatiotemporal occurrence and governing physics to unravel their relation to associated subduction dynamics and long term deformation we developed a 2 d continuum viscoelastoplastic model that uses an eulerian lagrangian finite difference framework with similar on and off fault physics we extend the validation of this numerical tool to a realistic subduction zone setting that resembles southern chile the resulting quasi periodic pattern of quasi characteristic m8 m9 megathrust events compares quantitatively with observed recurrence and earthquake source parameters albeit at very slow coseismic speeds without any data fitting surface displacements agree with gps data recorded before and during the 2010 m8 8 maule earthquake including the presence of a second order flexural bulge these surface displacements show cycle to cycle variations of slip deficits which overall accommodate 5 of permanent internal shortening we find that thermally and stress driven creep governs a spontaneous conditionally stable downdip transition zone between temperatures of 350 c and 450 c ruptures initiate above it and below the forearc moho propagate within it interspersed by small intermittent events and arrest below it as ductile shearing relaxes stresses ruptures typically propagate upward along lithological boundaries and widen as pressures drop the main thrust is constrained to be weak due to fluid induced weakening required to sustain regular subduction and to generate events with natural characteristics fluid pressures of 75 99 of solid pressures the agreement with a range of seismological geodetic and geological observations demonstrates the validity and strength of this physically consistent seismo thermo mechanical approach

shocks in a box an analogue model of subduction earthquake cycles with application to seismotectonic forearc evolution 1 we introduce and test an experimental approach to simulate elastoplastic megathrust earthquake cycles using an analogue model and apply it to study the seismotectonic evolution of subduction zones the quasi two dimensional analogue model features rate and state dependent elastic frictional plastic and viscoelastic material properties and is scaled for gravity inertia elasticity friction and viscosity the experiments are monitored with a high resolution strain analysis tool based on digital image correlation particle imaging velocimetry piv providing deformation time series comparable to seismologic geodetic and geologic observations in order to separate elastic and nonelastic effects inherent the experimental deformation patterns we integrate elastic dislocation modeling edm into a hybrid approach we use the analogue earthquake slip and interseismic locking distribution as edm dislocation input and forward model the coseismic and interseismic elastic response the residual which remains when the edm prediction is subtracted from the experimental deformation pattern highlights the accumulation of permanent deformation in the model the setup generates analogue earthquake sequences with realistic source mechanisms and elastic forearc response and recurrence patterns and reproduces principal earthquake scaling relations by applying the model to an accretionary type plate margin we demonstrate how strain localization at the rupture peripheries may lead to a seismotectonically segmented forearc including a tectonically stable shelf and coastal high 20 plate convergence accommodated by internal shortening overlying the area of large megathrust earthquake slip fifty to 75 of plate convergence is accommodated by internal shortening in the slope region where earthquake slip tapers out toward the trench the inner forearc region remains undeformed and represents a basin

the 2013 mw 7 7 balochistan earthquake energetic strike slip reactivation of a thrust fault we analyse the mw 7 7 balochistan earthquake of 09 24 2013 based on ground surface deformation measured from sub pixel correlation of landsat 8 images combined with back projection and finite source modeling of teleseismic waveforms the earthquake nucleated south of the chaman strike slip fault and propagated southwestward along the hoshab fault at the front of the kech band the rupture was mostly unilateral propagated at 3 km s on average and produced a 200 km surface fault trace with purely strike slip displacement peaking to 10 m and averaging around 6 m the finite source model shows that slip was maximum near the surface although the hoshab fault is dipping by 45 to the north in accordance with its origin as a thrust fault within the makran accretionary prism slip was nearly purely strike slip during that earthquake large seismic slip on such a non optimally oriented fault was enhanced possibly due to the influence of the free surface on dynamic stresses or to particular properties of the fault zone allowing for strong dynamic weakening strike slip faulting on thrust fault within the eastern makran is interpreted as due to eastward extrusion of the accretionary prism as it bulges out over the indian plate portions of the makran megathrust some thrust faults in the kirthar range and strike slip faults within the chaman fault system have been brought closer to failure by this earthquake aftershocks cluster within the chaman fault system north of the epicenter opposite to the direction of rupture propagation by contrast few aftershocks were detected in the area of maximum moment release in this example aftershocks cannot be used to infer earthquake characteristics

modelling hanging wall accommodation above rigid thrust ramps abstract experimental models are used to study the role of material rheology in hanging wall accommodation above rigid flat ramp flat thrust footwalls the deformation in the hanging wall was accomplished by forwards sliding along a rigid basal staircase trajectory with a variable ramp angle ranging from 15 to 60 we model different ramp angles to examine hanging wall accommodation styles above thrust ramps of overthrust faults ranging from 15 to 30 as well as above pre existing normal faults ranging from 45 to 60 for the hanging walls we used stratified frictional sand and viscous silicone putty materials in this paper we study three types of models type 1 models represent purely frictional hanging walls where accommodation above thrust ramps was by layer parallel thickening and by generating a series of back thrusts type 2 and 3 models represent stratified frictional viscous hanging walls in these models accommodation was by a complex association of reverse and normal faults mainly controlled by the rheological anisotropy as well as by the ramp inclination angle in type 2 models the silicone covered only the lower flat while in type 3 models it also covered the rigid ramp for 30 in type 2 models and 45 in type 3 models the viscous layer inhibited the development of back thrusts in the frictional hanging wall instead the silicone thickened to develop a ductile ramp for values higher than 30 in type 2 models and 45 in type 3 models back thrusts develop in response to the bulk compression the experiments simulate many structures observed above natural thrust ramps with 30 and pre existing normal faults with 45 the models emphasise the importance of a basal ductile layer which allows the hanging wall to step up over the rigid ramp by building up its own ductile ramp the models also emphasise that foreland directed normal faulting can develop at a thrust front in the case that the vertical stress due to gravity exceeds the horizontal stress due to end loading within a thrust wedge

deformation due to the 17 august 1999 izmit turkey earthquake measured from spot images the geometry of the ruptured areas and the coseismic slip distribution data are key to highlighting the behavior of seismic faults this information is generally retrieved from field investigations and geodetic measurements or synthetic aperture radar sar interferometry here we show that spot images can also be used to accurately map the fault zone and to determine the slip distribution by subpixel correlation of images acquired before and after an earthquake the measured slip includes the contribution of possible distributed shear that might not be clearly expressed in surface ruptures and smoothes out possible along strike variability due to near surface fault complexities we apply the technique to the m s 7 4 1999 izmit earthquake our results reveal a 100 m wide and very linear fault zone that can be traced for 70 km from golcuk to akyazi along which supershear rupture has been inferred the obtained slip distribution compares well with the field measurements and is consistent with ground deformation measured at some distance from the fault zone using sar images very little deformation was accommodated off the main fault plane maximum slip is observed near sapanca lake at a small fault jog that has probably influenced rupture propagation

present day kinematics at the india asia collision zone the collision of the indian subcontinent with asia drives the growth and evolution of the greater tibetan plateau region fault slip rates resulting from the relative motion between crustal blocks can provide a kinematic description of the distribution of presentday deformation i construct a three dimensional regional scale elastic block model of the india asia collision zone that is consistent with geodetic observations of interseismic deforma tion mapped fault system geometry historical seismicity and the mechanics of the earthquake cycle this mechanical model of the elastic upper crust yields a set of kinematically consistent fault slip rates and block motions that may serve to constrain dynamic models of continental crustal dynamics

experiments on shortening of a 4 layer model of the continental lithosphere davy ph and cobbold p r 1991 experiments on shortening of a 4layer model of the continental lithosphere in p r cobbold editor experimental and numerical modelling of continental deformation tecronophysics 188 l 25 we present an experimental method for modelhng lithospheric deformation the method is based on a simplified profile of lithosphere rheology where only four layers are considered brittle crust ductile crust brittle mantle and ductile mantle the strengths of the four layers calculated from laboratory determined flow laws depend strongly on geothenns and slightly on strain rates and deformation regimes the rheological profiles of the lithosphere have been recreated in experimental models built with analog materials sand silicone putty and golden syrup which flow under their own weight during lithospheric shortening the style of deformation depends on the relative strengths of the four layers especially those in the mantle of special importance is the existence or non existence of a brittle mantle the initial stages of compression often produce flexures or buckles without local isostatic compensation buckling exerts a strong control on the localisation of thrust faults in later stages the thrusts tend to form at the inflexion points of buckles the style of crustal thickening depends on the relative thicknesses and resistances of the lithospheric layers in 2 layer lithospheric models without resistant mantle shortening was concentrated upon a small number of thrusts widely spaced in 3 layer models shortening was concentrated into a large number of reverse faults more closely spaced and more steeply dipping than before in clayer models a single thrust formed in the brittle mantle motion on this thrust dragged the crustal layers down forming an asymmetric crustal root with steep imbricate slices

mantle flow beneath a continental strike slip fault postseismic deformation after the 1999 hector mine earthquake two recent large earthquakes in the mojave desert california the magnitude 7 3 1992 landers and magnitude 7 1 1999 hector mine earthquakes have each been followed by elevated crustal strain rates over periods of months and years geodetic data collected after the hector mine earthquake exhibit a temporally decaying horizontal velocity field and a quadrant uplift pattern opposite to that expected for localized shear beneath the earthquake rupture we interpret the origin of this accelerated crustal deformation to be vigorous flow in the upper mantle in response to the stress changes generated by the earthquake our results suggest that transient flow in the upper mantle is a fundamental component of the earthquake cycle and that the lower crust is a coherent stress guide coupling the upper crust with the upper mantle

on the mechanics of earthquake afterslip we propose a model for earthquake afterslip based on rate and state variable friction laws in the model afterslip is attributed to the interaction of a velocity weakening region at depth within which earthquakes nucleate with an upper region of velocity strengthening frictional behavior the existence of this upper region is supported by independent seismologic observations and the results of laboratory friction experiments in our model afterslip is the result of relaxation of a stress perturbation within the velocity strengthening region which arises when an earthquake propagates into that region from below we derive the stress perturbation and its decay from the friction constitutive law using a simple 1 degree of freedom approximation for the elastic interaction between the fault and its surroundings this approximation is based on thickness averaged displacements and slip velocities within the velocity strengthening region which is assumed to slip as a rigid block coseismic and postseismic slip are coupled through the thickness averaged stiffness k of the velocity strengthening region we assume k to be inversely proportional to the thickness of this region which means that thicker velocity strengthening regions have a greater tendency to arrest coseismic slip we model the afterslip time histories of the 1966 parkfield and 1987 superstition hills earthquakes and relate the model parameters to physical parameters which may govern the rheologic behavior of the faults in accord with field observations our model predicts 1 that afterslip on some faults scales with the thickness of the unconsolidated sedimentary cover and 2 that proportionally more afterslip occurs for earthquakes in which coseismic surface slip is small compared with coseismic slip at depth velocity strengthening frictional behavior is to be expected for faults within poorly consolidated sediments and for those that contain significant gouge zones about 500 m within their shallow regions 3 5 km combining our results with those of recent laboratory friction studies indicates that relatively young faults with little accumulated fault gouge should exhibit little afterslip

centrifuge modelling of the evolution of low angle detachment faults from high angle normal faults abstract centrifuge models composed of a ductile layer overlain by a semi brittle layer are used to study how deformation localised by a high angle normal fault promotes detachment faulting during lateral extension driven by centrifugal force localised extension along a pre existing fault initiated localised isostatic upwelling of the denser lower layer where the lower ductile layer was significantly less dense than the semi brittle upper layer localised extension along the prescribed cut initiated upwelling of the ductile lower layer based on model results we argue that the transition from high angle normal faulting to low angle detachment faulting is an inevitable consequence of localising extension provided that there is viscous coupling between the extending upper layer and the upwelling lower layer in models with a lower layer of equal density or a denser lower layer this rotation takes place at the later stages of localised thinning in the upper semi brittle layer whereas in models with a less dense lower layer the rotation takes place earlier due to the buoyant rise of the ductile lower layer in areas of distributed crustal stretching e g rift basins where extension of the upper layer is accommodated by numerous steep faults distributed over a wide area and upwelling of lower ductile materials is distributed across the area normal faults remain more planar despite a large amount of extension models show that distributed extension along several closely spaced normal faults encourages rotation of blocks rather than their distortion to form listric faults we further conclude that the only configuration whereby localisation of extension would not result in detachment style faulting is when the upper and lower layers were completely decoupled

submarine fault scarps in the sea of marmara pull apart north anatolian fault implications for seismic hazard in istanbul earthquake scarps associated with recent historical events have been found on the floor of the sea of marmara along the north anatolian fault naf the marmarascarps cruise using an unmanned submersible rov provides direct observations to study the fine scale morphology and geology of those scarps their distribution and geometry the observations are consistent with the diversity of fault mechanisms and the fault segmentation within the north marmara extensional step over between the strike slip ganos and izmit faults smaller strike slip segments and pull apart basins alternate within the main step over commonly combining strike slip and extension rapid sedimentation rates of 1 3 mm yr appear to compete with normal faulting components of up to 6 mm yr at the pull apart margins in spite of the fast sedimentation rates the submarine scarps are preserved and accumulate relief sets of youthful earthquake scarps extend offshore from the ganos and izmit faults on land into the sea of marmara our observations suggest that they correspond to the submarine ruptures of the 1999 izmit mw 7 4 and the 1912 ganos ms 7 4 earthquakes while the 1999 rupture ends at the immediate eastern entrance of the extensional cinarcik basin the 1912 rupture appears to have crossed the ganos restraining bend into the sea of marmara floor for 60 km with a right lateral slip of 5 m ending in the central basin step over from the gulf of saros to marmara the total 1912 rupture length is probably about 140 km not 50 km as previously thought the direct observations of submarine scarps in marmara are critical to defining barriers that have arrested past earthquakes as well as defining a possible segmentation of the contemporary state of loading incorporating the submarine scarp evidence modifies substantially our understanding of the current state of loading along the naf next to istanbul coulomb stress modeling shows a zone of maximum loading with at least 4 5 m of slip deficit encompassing the strike slip segment 70 km long between the cinarcik and central basins that segment alone would be capable of generating a large magnitude earthquake mw 7 2 other segments in marmara appear less loaded

the distribution of 234u 238u activity ratios in river sediments abstract uranium u isotopes can be used to estimate the comminution age of sediments i e the time elapsed from sediment production on continents via weathering and physical erosion to deposition in the sedimentary record the calculation of this comminution age is based on measured 234u 238u activity ratios in river sediments and inferred time dependent recoil effect which leads to the preferential release of 234u from mineral lattices during erosion processes in this study we report on a large scale 234u 238u investigation of modern river sediments worldwide with the aim to determine the extent to which parameters such as grain size lithology weathering climate and geomorphology may influence the distribution of u isotopes in fine grained sediments our extensive dataset n 64 includes u isotopic measurements for many of the world s largest rivers but also rivers draining particular climatic and geological settings our results indicate that sediments collected from river basins draining mostly igneous metamorphic or volcanic rocks often display 234u 238u ratios 1 with clay size fractions

stress transfer and strain rate variations during the seismic cycle 1 the balance of forces implies stress transfers during the seismic cycle between the elastobrittle upper crust and the viscoelastic lower crust this could induce observable time variations of crustal straining in the interseismic period we simulate these variations using a one dimensional system of springs sliders and dashpot loaded by a constant force the seismogenic zone and the zone of afterslip below are modeled from rate and state friction the ductile deeper fault zone is modeled from a viscous slider with newtonian viscosity n the force per unit length f must exceed a critical value fc to overcome friction resistance of the fault system this simple system produces periodic earthquakes the recurrence period tcycle and the duration of the postseismic relaxation phase which is driven dominantly by afterslip then both scale linearly with n between two earthquakes interseismic strain buildup across the whole system is nonstationary with the convergence rates vi just after each earthquake being systematically higher than the value vf at the end of the interseismic period we show

case for very low coupling stress on the cascadia ssubduction fault a fundamental problem in plate tectonics is the shear strength of major plate boundary faults this translates to the question whether the generally observed small earthquake stress drops of 3 10 mpa on major faults release most of the accumulated stress or only a small fraction of it there is strong evidence that the san andreas fault a major transform plate boundary is weak 20 mpa shear resistance it is not yet clear whether subduction thrust faults are also weak we present two types of evidence from the northern cascadia subduction zone that indicate very low coupling shear stress on that plate interface and hence very low strength of the subduction thrust fault comparable to that estimated for the san andreas fault first the well defined surface heat flow and heat generation allow negligible frictional heating on the plate interface the average shear stress on the fault must thus be very low over a time scale of a few million years second focal mechanism solutions for small crustal earthquakes in the southern vancouver island area indicate that the horizontal stress in the direction of plate convergence has a similar magnitude to the vertical stress this inferred stress state requires the present tectonic stress coupled across the subduction thrust fault to be very low one explanation for the weakness of the fault is the presence of near lithostatic pore fluid pressure in the region of the fault zone for which there is independent evidence the conclusion of a weak subduction thrust fault does not conflict with geodetic observations of contemporary surface deformation which indicate that the fault is currently locked accumulating strain energy toward a future great earthquake the surface deformation responds to the small 20 mpa temporal changes of the stress field associated with the subduction earthquake cycle this transient stress is superimposed on the larger background regional stress field in which the maximum compression is parallel to the margin the weakness of the cascadia subduction thrust fault and the unusual stress state of the forearc region have important implications for earthquake hazards for example a subduction earthquake may induce large strike slip earthquakes in the forearc that affect a large area

effective friction law for small scale fault heterogeneity in 3d dynamic rupture 1 we address the problem of modeling dynamic rupture on multiscale heterogeneous faults in 3d under the assumption of slip weakening friction we numerically construct effective friction laws that integrate the effects of small scale heterogeneity during the rupture this homogenization process is based on the description of the initial phase of the rupture by the dominant unstable spectral mode its dynamics is influenced by the geometry of the fault the static friction heterogeneities and the friction law we first define a periodic small scale heterogeneous model introducing heterogeneity in the distribution of the static friction coefficient we then describe a method for constructing this effective friction law applying this new law homogeneously on the fault permits to reproduce the dynamic evolution of the heterogeneous fault furthermore we show that the effective friction law can be used to replace small scale heterogeneities in two scale heterogeneous models while preserving their effects we study three kinds of two scale models with growing complexity first periodic at both scales then periodic only at small scale and finally irregular at both scales this homogenization method can be adapted to the case where the heterogeneity is introduced in the initial stress rather than in the static friction value finally we show in a simple example that the effective friction law permits to reproduce the transition between subshear and supershear rupture propagation originally produced by heterogeneities on the fault

constraining earthquake source inversions with gps data 1 resolution based removal of artifacts 1 we present a resolution analysis of an inversion of gps data from the 2004 mw 6 0 parkfield earthquake this earthquake was recorded at thirteen 1 hz gps receivers which provides for a truly coseismic data set that can be used to infer the static slip field we find that the resolution of our inverted slip model is poor at depth and near the edges of the modeled fault plane that are far from gps receivers the spatial heterogeneity of the model resolution in the static field inversion leads to artifacts in poorly resolved areas of the fault plane these artifacts look qualitatively similar to asperities commonly seen in the final slip models of earthquake source inversions but in this inversion they are caused by a surplus of free parameters the location of the artifacts depends on the station geometry and the assumed velocity structure we demonstrate that a nonuniform gridding of model parameters on the fault can remove these artifacts from the inversion we generate a nonuniform grid with a grid spacing that matches the local resolution length on the fault and show that it outperforms uniform grids which either generate spurious structure in poorly resolved regions or lose recoverable information in well resolved areas of the fault in a synthetic test the nonuniform grid correctly averages slip in poorly resolved areas of the fault while recovering small scale structure near the surface finally we present an inversion of the parkfield gps data set on the nonuniform grid and analyze the errors in the final model

faulting characteristics of supershear earthquakes article i nfo recent observations of large continental strike slip earthquakes show that different fault segments may rupture at different speeds and that rupture may propagate faster than the shear wave velocity of surrounding rocks we report that all the documented occurrences of supershear rupture are associated with faults which have simple geometry the surface trace of these faults is described in the field or seen on satellite images as remarkably linear continuous and narrow segmentation features along these segments are small or absent and the deformation is highly localized as theoretically predicted supershear is clearly associated with mode ii rupture

coseismic deformation derived from analyses of c and l band sar data and fault slip inversion of the yushu ms7 1 earthquake china in 2010 abstract we obtained the coseismic deformation field of the ms7 1 earthquake in yushu qinghai china on 14 april 2010 using l band and c band sar images the results show that the deformation fields derived from l and c bands separately are consistent in general they are all characterized by concentric elliptic interference fringes surrounding a nw trending fault with two local intensively deformed regions meanwhile they have small differences in area coverage of the deformation field magnitude of los displacements and other details of deformation nearby the fault the deformation from c band and l band covers 89 59 km square and 77 43 km square with maximum displacements of about 45 cm and 65 cm in the radar line of sight respectively using these two kinds of insar measurements separately and jointly we constructed a three segment fault model and inverted the coseismic fault slip of the yushu event the inversion results show that the slip distributions constrained by the two kinds of insar data are in agreement on the whole they all indicate two concentration regions of slip distribution one is located near the jiegu town on the southeastern segment of the yushu fault with a large area and the maximum slip of 2 4 m and the other is at the epicenter on the northwestern segment of the fault with a relatively small area there are also discrepancies between the two kinds of slip distribution models in general the scale magnitude and depth of slip distribution constrained by c band asar data are larger than that by l band palsar data the slip distribution from a joint inversion of the two kinds of sar data is considered to be the most robust model compared to the results of that constrained by l band sar data or c band sar data alone

coseismic and postseismic fault slip for the 17 august 1999 m 7 5 izmit turkey earthquake we use global positioning system gps observations and elastic half space models to estimate the distribution of coseismic and postseismic slip along the izmit earthquake rupture our results indicate that large coseismic slip reaching 5 7 meters is confined to the upper 10 kilometers of the crust correlates with structurally distinct fault segments and is relatively low near the hypocenter continued surface deformation during the first 75 days after the earthquake indicates an aseismic fault slip of as much as 0 43 meters on and below the coseismic rupture these observations are consistent with a transition from unstable episodic large earthquakes to stable fault creep sliding at the base of the seismogenic zone

fault stability inferred from granite sliding experiments at hydrothermal conditions seismicity on crustal faults is concentrated in the depth interval 1 3 to 12 15km tse and rice 1986 suggested that the lower bound on seismicity is due to a switch with increasing temperature from velocity weakening destabilizing to velocity strengthening stabilizing friction they inferred this transition from friction data for dry granite however pore fluids exist at elevated temperatures throughout the crust and may strongly influence strength and sliding behavior we present new data from sliding experiments on granite at elevated t 23 to 600 plus elevated ph2o 100 mpa our results show velocity strengthening at room temperature but velocity weakening from 100 to 350 c except at 250 from 350 to 600 there are systematic trends from velocity weakening to strong velocity strengthening and from high to low friction neither trend was seen in tests on dry granite the velocity dependence data imply the potential for unstable slip in the interval 100 to 350 using a geotherm to map temperature to depth this interval closely matches the observed earthquake distribution

strain accumulation across the carrizo segment of the san andreas fault california impact of laterally varying crustal properties 1 major strike slip faults juxtapose geologically dissimilar terrain which may vary in mechanical properties leading to an asymmetric pattern of strain accumulation we present new gps data on the carrizo segment of the san andreas fault separating the salinian block southwest of the fault from franciscan terrane northeast of the fault to better quantify asymmetric strain accumulation we also present a series of finite element models to investigate the possible role of variable elastic layer thickness and material properties of the upper crust the geodetic data are well fit with a simple model comprising a weak upper crustal zone 10 25 km wide northeast of the fault this model is also consistent with geologic data on the distribution of major rock types and corresponding laboratory data on their material properties as well as paleoseismic seismic and magnetotelluric data using this model we estimate a long term average over several seismic cycles slip rate for the san andreas fault of 36 1 5 2 mm yr in agreement with the known holocene rate within uncertainties and a viscosity for the combined lower crust upper mantle of 2 5 1019 pa s

experimental observation of stick slip instability fronts summary stick slip dynamic instability is a key mechanism governing frictional processes from microscale physics to earthquake faults and landslides yet challenging questions are still open about its nucleation and propagation dynamics here we present novel observations on laboratory experimental faults where spontaneously nucleating fractures are produced describing 1 an initial quasistatic stable rupture front propagating at about 5 per cent the shear wave velocity vs accelerating to subshear and then to intersonic velocity 2 the arisal of a higher degree of complexity when the friction to prestress ratio is increased on the sliding surface the complex behaviour includes stop and go sequences irregular propagation and rerupturing episodes within short time intervals implying rapid restrengthening of the surface and the formation of self healing pulses reproducing experimentally for the first time a behaviour observed on seismic faults

gelatins as rock analogs a systematic study of their rheological and physical properties abstract in a laboratory model the behavior of the selected analog material must be properly scaled to nature here we systematically investigate both rheological and physical properties of a wide range of gelatins as functions of temperature composition concentration ageing and applied strain rate results show that gelatins behavior changes gradually from purely elastic to visco elasto brittle to purely viscous nonlinear rheology going from the gel to the sol state the rheological variability of gelatins appears promising for the potential use of these analog materials to simulate crustal and lithospheric rheological behavior in particular we found that pig skin 2 5 wt at 10 c has the required rheological properties for a suitable experimental set up to model crustal deformation

space and time evolution of rupture and faulting during the 1999 izmit turkey earthquake we use the records of the ground motion obtained at near fault accelerometers to study the space and time evolution of rupture and faulting during the izmit earthquake we find that the rupture propagated at the sub rayleigh speed of about 3 km sec on the western and eastern segments of the fault but that the central segment izmit sapanca lake sakarya nearly 50 km long broke at the supershear speed of about 4 8 km sec this value within the range of uncertainties is the one theoretically predicted sqrt 2 v mathrm s in fracture dynamics for stable shear crack growth at intersonic speed we infer an average fault slip of about 2 9 m over a total rupture length of about 150 km with the largest values of up to 6 m occurring in the golcuk area to the west and in the sakarya region to the east the strong motion data also indicate that the slip diminished gradually to the west beyond the hersek peninsula over about 30 km whereas it stopped abruptly at depth at the termination of the eastern karadere segment the slip duration is between 2 and 4 sec except in the hypocentral area which slipped in about 1 sec the seismic moment inferred is about 2 5 10 20 n m

the 2005 mw 7 6 kashmir earthquake sub pixel correlation of aster images and seismic waveforms analysis we analyze the m w 7 6 kashmir earthquake of october 8 2005 using sub pixel correlation of aster images to measure ground deformation and modeling seismic waveforms the surface rupture is continuous over a distance of 75 km and cuts across the hazara syntaxis reactivating the tanda and the muzaffarabad faults north of muzaffarabad the surface rupture coincides approximately with the mbt on the southwestern flank of the syntaxis although the two faults have opposite dip angles the rupture terminates abruptly at the hairpin turn of the mbt showing a strong structural control the fault offset is 4 m on average and peaks to 7 m northwest of muzaffarabad the rupture lasted about 25 s and propagated updip and bi laterally 2 km s with a rise time of 2 5 s the shallowness and compactness of the rupture both in time and space provide an explanation for the intensity of destructions this kind of analysis could be achieved as soon as a post earthquake image is available and would provide key information for early assessment of damages the study sheds some light on seismic hazard in the himalaya and raises concern regarding the possibility of a repetition of the 1555 event which presumably ruptured the himalayan front south of the kashmir basin and may have reached a magnitude m w 8

numerical modeling of landquakes the thurwieser landslide that occurred in italy in 2004 and the seismic waves it generated have been simulated and compared to the seismic signal recorded a few tens of kilometers from the landslide source i e landquake the main features of the low frequency seismic signal are reproduced by the simulation topography effects on the flowing mass have a major impact on the generated seismic signal whereas they weakly affect low frequency wave propagation simulation of the seismic signal makes it possible to discriminate between possible alternative scenarios for flow dynamics and to provide first estimates of the rheological parameters during the flow as landquakes are continuously recorded by seismic networks our results provide a new way to collect data on the dynamics and rheology of natural flows

predictive modeling of the seismic cycle of the greater san francisco bay region the seismic cycle for the san francisco bay region is synthesized by a model combining the pre and post 1906 seismic histories the long term acceleration of seismic release seismic moment benioff strain release or event count in the seismic cycle and the shorter term accelerations preceding the larger earthquakes within that cycle are modeled using an empirical predictive technique called time to failure analysis in which rate of seismic release is proportional to an inverse power of the remaining time to failure the exponent of time to failure in the accelerating sequences appears to be scale invariant and the length of the full cycle is estimated at 269 50 years the 1989 loma prieta earthquake which is the culmination of the first subcycle in the present long term seismic cycle should have been predictable with an uncertainty of 2 years in time and 0 5 in magnitude although the specific location at loma prieta was not predictable by this technique if our model is correct and if the loma prieta earthquake is the culmination of a subcycle the san francisco bay region should be entering a relatively long 20 50 years period of seismic quiescence above magnitude 6 a great earthquake such as the 1906 san francisco event would appear to be more than a century in the future

laboratory earthquakes the sub rayleigh to supershear rupture transition we report on the experimental observation of spontaneously nucleated supershear rupture and on the visualization of sub rayleigh to supershear rupture transitions in frictionally held interfaces the laboratory experiments mimic natural earthquakes the results suggest that under certain conditions supershear rupture propagation can be facilitated during large earthquake events

evidence of nonlinear elasticity of the crust from the mw7 6 manyi tibet earthquake satellite synthetic aperture radar sar interferometry shows that the magnitude 7 6 manyi earthquake of 8 november 1997 produced a 170 kilometer long surface break with up to 7 meters of left lateral slip reactivating a n76 e quaternary fault in western tibet the radar interferometric map reveals asymmetric along strike displacement profiles between the two sides of the surface rupture a pattern that cannot be explained with linear elastic theory this observation suggests that the elastic moduli of the crust in tension and in compression are different because of the presence of cracks in the crust at shallow depth a model indicates that a ratio of 2 between compressive and tensile elastic moduli can account for the observed asymmetry a ratio that is consistent with laboratory and borehole measurements

earthquake supercycles inferred from sea level changes recorded in the corals of west sumatra records of relative sea level change extracted from corals of the mentawai islands sumatra imply that this 700 kilometer long section of the sunda megathrust has generated broadly similar sequences of great earthquakes about every two centuries for at least the past 700 years the moment magnitude 8 4 earthquake of september 2007 represents the first in a series of large partial failures of the mentawai section that will probably be completed within the next several decades

sm nd dating and cooling history of scourian granulites sutherland the scourian gneisses of the lewisian complex of north west scotland have been the subject of many geochronological investigations sm nd whole rock measurements suggest that the protoliths of the lewisian complex differentiated 2 92 0 05 gyr ago from an approximately chrondritic mantle1 rb sr and pb pb whole rock ages between 2 8 and 2 6 gyr reflect subsequent depletion in rb and u refs 2 3 which has generally been associated with granulite facies metamorphism zircons which are thought to have crystallized during this metamorphism have also given u pb ages of 2 66 gyr ref 4 this granulite facies metamorphism involving pressures in excess of 10 kbar and peak temperatures estimated between 820 50 c ref 5 and 1 250 c refs 6 7 affects a terrain including supracrustal rocks7 and thus implies a major tectonic event we report here the first geochronological data measured on mineral assemblages that have also provided estimates of the conditions of metamorphism such measurements are possible because the major minerals of granulite facies assemblages show an unusually favourable spread in sm nd ratios garnet especially is strongly enriched in sm ref 8 data for three garnetiferous basic gneiss assemblages suggest that closure of the sm nd system occurred significantly later than the peak of the metamorphism

from sub rayleigh to supershear ruptures during stick slip experiments on crustal rocks supershear earthquake ruptures propagate faster than the shear wave velocity although there is evidence that this occurs in nature it has not been experimentally demonstrated with the use of crustal rocks we performed stick slip experiments with westerly granite under controlled upper crustal stress conditions supershear ruptures systematically occur when the normal stress exceeds 43 megapascals mpa with resulting stress drops on the order of 3 to 25 mpa comparable to the stress drops inferred by seismology for crustal earthquakes in our experiments the sub rayleigh to supershear transition length is a few centimeters at most suggesting that the rupture of asperities along a fault may propagate locally at supershear velocities in turn these sudden accelerations and decelerations could play an important role in the generation of high frequency radiation and the overall rupture energy budget

evidence for self similar triangular slip distributions on earthquakes implications for earthquake and fault mechanics we characterize average slip distributions on earthquakes beyond their individual heterogeneity for that we analyze a large number of seismic slip distributions both measured at the surface after earthquakes 44 profiles and derived from slip inversion models 76 models investigating the overall shape of these slip profiles we find that they are roughly triangular both along strike and dip and most of them 70 80 are asymmetric long linear slopes and high slip gradients therefore are the key ingredients to describe earthquake slip profiles the scaling relations between maximum displacement and length or width suggest furthermore that the triangular slip profiles are self similar such slip patterns make earthquakes dominated by one major zone of maximum slip hence one major asperity analyzing the position of hypocenters with respect to these asperities we find that earthquakes nucleate at a distance from them that averages 20 30 of their total length compiling observations on 56 earthquakes we show that this distance i e the asperity size is structurally defined we then compare the earthquake slip profiles to cumulative slip profiles measured on long term faults of various ages and sizes and find that all profiles have a similar shape triangular and asymmetric hence combining data for a large number of earthquakes leads to point out average generic characteristics of the coseismic slip that are similar to those that emerge from the accumulation of events with time on a single fault this suggests that these characteristics result from robust physical properties

the seismic cycle at subduction thrusts 1 insights from laboratory models 1 subduction megathrust earthquakes occur at the interface between the subducting and overriding plates these hazardous phenomena are only partially understood because of the absence of direct observations the restriction of the instrumental seismic record to the past century and the limited resolution completeness of historical to geological archives to overcome these restrictions modeling has become a key tool to study megathrust earthquakes we present a novel model to investigate the seismic cycle at subduction thrusts using complementary analog paper 1 and numerical paper 2 approaches here we introduce a simple scaled gelatin on sandpaper setup including realistic tectonic loading spontaneous rupture nucleation and viscoelastic response of the lithosphere particle image velocimetry allows to derive model deformation and earthquake source parameters analog earthquakes are characterized by quasi periodic recurrence consistent with elastic theory the interseismic stage shows rearward motion subsidence in the outer wedge and uplift of the coastal area as a response of locked plate interface at shallow depth the coseismic stage exhibits order of magnitude higher velocities and reversal of the interseismic deformation pattern in the seaward direction subsidence of the coastal area and uplift in the outer wedge like natural earthquakes analog earthquakes generally nucleate in the deeper portion of the rupture area and preferentially propagate upward in a crack like fashion scaled rupture width slip proportionality and seismic moment duration scaling verifies dynamic similarities with earthquakes experimental repeatability is statistically verified comparing analog results with natural observations we conclude that this technique is suitable for investigating the parameter space influencing the subduction interplate seismic cycle

recent plate motions and crustal deformation reports by u s workers on geodetic measurements of recent plate motions or crustal deformation published in 1987 1990 are reviewed the review begins with global plate motions proceeds through plate boundaries in california alaska and the pacific northwest and finishes with volcanic phenomena monument stability and longevity and gps relative position measurements 184 refs

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

crustal earthquake instability in relation to the depth variation of frictional slip properties recent stability studies using constitutive relations of the type found by dieterich ruina and others to describe frictional slip of rocks in the laboratory have provided a new explanation of the depth cutoff of shallow crustal earthquakes the class of friction laws discussed has the property that the sliding stress depends on normal stress temperature slip rate and slip history for sliding at a fixed slip rate v and fixed environment e g normal stress temperature etc the shear strength evolves toward a steady state value ss v stability analysis show that for d ss v dv 0 velocity strengthening implies stable steady state sliding at least for a broad class of constitutive laws more fully described in the paper experiments by dieterich and by tullis and weeks on westerly granite with mature sliding surfaces indicate that d ss v dv is negative at room temperature whereas higher temperature experiments by stesky show that d ss v dv becomes positive above approximately 300 c therefore in the case of the earth where temperature increases with depth the above observations seem to suggest that the depth cutoff of crustal earthquake activity can be understood in terms of the variation of the frictional response with depth from a regime with d ss v dv 0 this is not inconsistent with but rather refines the suggestion by sibson and others that the depth cutoff is due to a transition from brittle friction to ductile flow further our results show definitively that the shallow depth confinement of seismicity is compatible with a model in which deformation is localized to a fault zone extending well below the seismogenic depth the depth confinement does not require a model showing transition to a zone of broadly distributed creep flow beneath the seismogenic zone following mavko a two dimensional quasi static strike slip fault model is analyzed but using different numerical procedures and a depth variation of frictional properties based on the laboratory data mentioned above and the lachenbruch sass depth variation of temperature for the san andreas fault the resulting predictions of such features as the confinement of crustal earthquakes to shallow depths the development of locked patches the recurrence time for the seismic cycle the seismic stress drop and displacement etc are generally in agreement with the observed characteristics of large scale strike slip earthquakes along the san andreas fault calculations based on the laboratory data predict for example that stick slip nucleates around 5 7 km that large seismic motion occurs around the nucleation depth and above but diminishes gradually to zero at 13 15 km depths that rapid postseismic creep occurs over another 3 4 km depth and that at greater depths steady slip consistent with the average plate velocity is only modestly perturbed by the earthquakes occurring above the most uncertain parameter is the slip weakening distance l for evolution of fault surface state we illustrate how features of the predicted earthquake cycles vary with l over the range for which calculations are feasible

on the use of dislocations to model interseismic strain and stress build up at intracontinental thrust faults creeping dislocations in an elastic half space are commonly used to model interseismic deformation at subduction zones and might also apply to major intracontinental thrust faults such as the main himalayan thrust here we compare such models with a more realistic 2 d finite element model that accounts for the mechanical layering of the continental lithosphere and surface processes and that was found to fit all available constraints on interseismic and long term surface displacements these can also be fitted satisfactorily from dislocation models the conventional back slip model commonly used for subduction zones may however lead to a biased inference about the geometry of the locked portion of the thrust fault we therefore favour the use of a creeping buried dislocation that simulates the ductile shear zone in the lower crust a limitation of dislocation models is that the mechanical response of the lithosphere to the growth of the topography by bending of the elastic cores and ductile flow in the lower crust cannot be easily introduced fortunately these effects can be neglected because we may assume to first order a stationary topography moreover we show that not only can dislocation models be used to adjust surface displacements but with some caution they can also provide a physically sound rationale to interpret interseismic microseismicity in terms of stress variations

stress coupling between earthquakes in northwest turkey and the north aegean sea we have investigated the coulomb stress interactions of 29 earthquakes m s greater than or equal to 6 0 that have occurred in the region of northwest turkey and north aegean sea since 1912 of these events 23 may be related to earlier events and 16 are clearly related to earlier events all events after 1967 are related to previous events events in the early part of our time interval that show no correlation could be related to historical events as yet unidentified in some cases faults that have received a stress reduction from earlier events are prepared for an event by an earthquake occurring a few years before that creates a local coulomb stress rise thus regions of coulomb stress shadow can become regions where a damaging earthquake may occur the relation between smaller events and the coulomb stress distribution is less clear but may be related to poor data quality and practical limitations of our modeling technique nonetheless there are 4 times as many events per unit area in regions of enhanced stress than where stress is reduced we discuss the contemporary distribution of coulomb stress and argue that it is possible to identify the likely locations of future damaging earthquakes including identifying the most likely candidate faults

creeping along the ismetpasa section of the north anatolian fault western turkey rate and extent from insar creeping along the north anatolian fault naf at ismetpasa turkey was discovered some thirty years ago about a decade after the first observations of the phenomenon along the san andreas fault in california however little is known about its lateral extent and rate in order to study its three dimensional nature and rupture characteristics we use synthetic aperture radar interferometry insar and elastic dislocation models compared also with field observations interferograms with temporal baselines ranging between 1 25 and 5 years show that the creeping section starts at the western termination of the 1943 m 7 6 earthquake rupture it continues about 70 km to the west overlapping with the eastern part of the 1944 m 7 3 earthquake rupture offsets along strike indicate a maximum creep rate of 11f3 mm year near the mid point of the creeping section decreasing gradually towards the edges near ismetpasa insar data yield 8f3 mm year of creep rate consistent with recent instrumental triangulation and creepmeter measurements modeling of the insar and gps data suggests that the faultcreep occurs most probably at a shallow depth 0 7 km our analysis combined with previous studies suggests that creeping might have commenced following the 1944 earthquake and thus may be a long lasting but transient slip episode d 2005 elsevier b v all rights reserved

nonlinear strain buildup and the earthquake cycle on the san andreas fault two contrasting models of the earthquake deformation cycle on strike slip faults predict significant temporal declines in shear strain rate near the fault accompanied by a progressive broadening of the zone of deformation adjacent to it in the thin lithosphere model transient deformation results from flow in the asthenosphere due to stress relaxation following faulting through most or all of the lithosphere for an earth model with a thick elastic lithosphere plate thickness depth of seismic slip transient motions are due to postearthquake aseisrnic slip below the coseismic fault plane data from the san andreas fault indicate a long term temporal decrease in strain rate that persists for at least 30 years and may extend through the entire earthquake cycle observations support a cycle long rate decrease and a temporal spreading of the deformation profile only if movement cycles on the northern and southern locked sections of the fault are basically similar if so the usually lower strain rates and broader deformation zone currently observed on the southern san andreas represent a later evolutionary stage of the northern locked section where a great earthquake is a more recent occurrence although the data allow some extreme models to be discarded no sufficiently strong constraints exist to decide between the thin and thick lithosphere models regardless of the appropriate model the geodetic observations themselves indicate that strain buildup is sufficiently nonlinear to cause significant departures from recurrence estimates based on linear strain accumulation and the time predictable model

joint inversion of insar gps teleseismic and strong motion data for the spatial and temporal distribution of earthquake slip application to the 1999 izmit mainshock the space time distribution of slip of the 17 august 1999 izmit earthquake is investigated by inverting synthetic aperture radar sar interferometry and global positioning system gps data together with teleseismic broadband and near field strong motion records surface offsets are used as an added constraint special emphasis is given to analysis of the resolution of the different data sets we use a four segment finite fault model and a nonlinear inversion scheme allowing slip to vary in amplitude direction and duration as well as variable rupture velocity from the inversion of synthetic data we find that the best spatial resolution can be expected in the upper half of the fault model above 12 km where coverage of the interferometric sar data is good western half of the rupture and near the gps and strong motion stations teleseismic data are found to have a lower resolution that is more evenly distributed over the fault model the joint inversion of all the data sets has an increased resolving power compared with the separate inversions and gives a more robust description of the space and time distribution of slip our study shows the importance of resolution tests in evaluating the reliability of earthquake kinematic models and it confirms that an excellent fit of a single kind of data does not necessarily imply a good retrieval of the kinematic properties of an earthquake the izmit rupture which is almost pure right lateral strike slip faulting is dominated by the bilateral breaking of a central asperity located between 29 7 e about 10 km west of the city of golcuk and 30 4 e eastern margin of sapanka lake with slip reaching 6 8 m in the depth range 6 12 km the western termination of the rupture is found near the city of yalova but large slip ends around 29 7 e about 10 km east of hersek delta a second area of large slip is required by all the data sets further east toward the city of duzce between 30 7 e and 31 1 e karadere and duzce faults this eastern slip zone which is separated from the main central asperity by an area of greatly reduced slip is less well constrained by the data however a strong motion station near the city of duzce helps to locate a high slip patch near 31 1 e in the depth range 6 12 km the total seismic moment resulting from the joint inversion is 2 4 1027 dyne cm most of the energy release occurred in a short time less than 15 sec corresponding to the bilateral breaking of the central asperity rupture propagation is relatively uniform and fast toward the west with a rupture velocity close to 3 5 km sec propagation of large slip toward the east is initially slower but it accelerates during a short time interval about 10 sec after rupture nucleation eastward progression then slows down to less than 2 km sec after 15 sec and rupture almost vanishes in amplitude ca 20 sec after initiation rupture propagation then proceeds on the easternmost karadere and duzce fault segments east of 30 7 e from 22 to ca 50 sec supershear rupture propagation is not required for modeling the waveforms considered in this study the hypocenter of the duzce earthquake which occurred 3 months later 12 november 1999 m w 7 2 is located in the immediate vicinity of the easternmost slip patch of the izmit earthquake manuscript received 28 july 2000

deformation cycles of subduction earthquakes in a viscoelastic earth subduction zones produce the largest earthquakes over the past two decades space geodesy has revolutionized our view of crustal deformation between consecutive earthquakes the short time span of modern measurements necessitates comparative studies of subduction zones that are at different stages of the deformation cycle piecing together geodetic snapshots from different subduction zones leads to a unifying picture in which the deformation is controlled by both the short term years and long term decades and centuries viscous behaviour of the mantle traditional views based on elastic models such as coseismic deformation being a mirror image of interseismic deformation are being thoroughly revised snapshots of subduction zones using space geodesy reveal that the viscous behaviour of the mantle controls crustal deformation requiring the revision of traditional elastic models for earthquake risk assessment space geodesy the use of satellites to monitor earth has revolutionized our view of crustal deformation between consecutive earthquakes however the brevity of such measurements means that studies must be done by comparing multiple subduction zones at different stages of the earthquake cycle by piecing together geodetic snapshots from the sumatra chile and cascadia subduction zones this review presents a unifying picture in which deformation is controlled by the viscous behaviour of the mantle in both the short term years and the long term decades to centuries

fault behavior and characteristic earthquakes examples from the wasatch and san andreas fault zones paleoseismological data for the wasatch and san andreas fault zones have led to the formulation of the characteristic earthquake model which postulates that individual faults and fault segments tend to generate essentially same size or characteristic earthquakes having a relatively narrow range of magnitudes near the maximum analysis of scarp derived colluvium in trench exposures across the wasatch fault provides estimates of the timing and displacement associated with individual surface faulting earthquakes at all of the sites studied the displacement per event has been consistently large measured values range from 1 6 to 2 6 m and the average is about 2 m on the basis of variability in the timing of individual events as well as changes in scarp morphology and fault geometry six major segments are recognized along the wasatch fault on the basis of the most likely number of surface faulting events 18 that have occurred on segments of the wasatch fault zone during the past 8000 years an average recurrence interval of 400 666 years with a preferred average of 444 years is calculated for the entire zone geologic data on the distribution of slip associated with prehistoric earthquakes and slip rates along the south central segment of the san andreas fault suggest that the m 8 1857 earthquake is a characteristic earthquake for this segment comparisons of earthquake recurrence relationships on both the wasatch and san andreas faults based on historical seismicity data and geologic data show that a linear constant b value extrapolation of the cumulative recurrence curve from the smaller magnitudes leads to gross underestimates of the frequency of occurrence of the large or characteristic earthquakes only by assuming a low b value in the moderate magnitude range can the seismicity data on small earthquakes be reconciled with geologic data on large earthquakes the characteristic earthquake appears to be a fundamental aspect of the behavior of the wasatch and san andreas faults and may apply to many other faults as well

a new multilayered visco elasto plastic experimental model to study strike slip fault seismic cycle nowadays technological advances in satellite imagery measurements as well as the development of dense geodetic and seismologic networks allow for a detailed analysis of surface deformation associated with active fault seismic cycle however the study of earthquake dynamics faces several limiting factors related to the difficulty to access the deep source of earthquake and to integrate the characteristic time scales of deformation processes that extend from seconds to thousands of years to overcome part of these limitations and better constrain the role and couplings between kinematic and mechanical parameters we have developed a new experimental approach allowing for the simulation of strike slip fault earthquakes and analyze in detail hundreds of successive seismic cycle model rheology is made of multilayered visco elasto plastic analog materials to account for the mechanical behavior of the upper and lower crust and to allow simulating brittle ductile coupling postseismic deformation phase and far field stress transfers the kinematic evolution of the model surface is monitored using an optical system based on subpixel spectral correlation of high resolution digital images first results show that the model succeed in reproducing the deformation mechanisms and surface kinematics associated to the main phases of the seismic cycle indicating that model scaling is satisfactory these results are comforted by using numerical algorithms to study the strain and stress distribution at the surface and at depth along the fault plane our analog modeling approach appears then as an efficient complementary approach to investigate earthquake dynamics

power law viscous materials for analogue experiments new data on the rheology of highly filled silicone polymers abstract the selection of appropriate analogue materials is a central consideration in the design of realistic physical models we investigate the rheology of highly filled silicone polymers in order to find materials with a power law strain rate softening rheology suitable for modelling rock deformation by dislocation creep and report the rheological properties of the materials as functions of the filler content the mixtures exhibit strain rate softening behaviour but with increasing amounts of filler become strain dependent for the strain independent viscous materials flow laws are presented while for strain dependent materials the relative importance of strain and strain rate softening hardening is reported if the stress or strain rate is above a threshold value some highly filled silicone polymers may be considered linear visco elastic strain independent and power law strain rate softening the power law exponent can be raised from 1 to 3 by using mixtures of high viscosity silicone and plasticine however the need for high shear strain rates to obtain the power law rheology imposes some restrictions on the usage of such materials for geodynamic modelling two simple shear experiments are presented that use newtonian and power law strain rate softening materials the results demonstrate how materials with power law rheology result in better strain localization in analogue experiments

elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate and state dependent friction we present an efficient and rigorous numerical procedure for calculating the elastodynamic response of a fault subjected to slow tectonic loading processes of long duration within which there are episodes of rapid earthquake failure this is done for a general class of rate and state dependent friction laws with positive direct velocity effect the algorithm allows us to treat accurately within a single computational procedure loading intervals of thousands of years and to calculate for each earthquake episode initially aseismic accelerating slip prior to dynamic rupture the rupture propagation itself rapid post seismic deformation which follows and also ongoing creep slippage throughout the loading period in velocity strengthening fault regions the methodology is presented using the two dimensional 2 d antiplane spectral formulation and can be readily extended to the 2 d in plane and 3 d spectral formulations and with certain modifications to the space time boundary integral formulations as well as to their discretized development using finite difference or finite element methods the methodology can be used to address a number of important issues such as fault operation under low overall stress interaction of dynamic rupture propagation with pore pressure development patterns of rupture propagation in events nucleated naturally as a part of a sequence the earthquake nucleation process earthquake sequences on faults with heterogeneous frictional properties and or normal stress and others the procedure is illustrated for a 2 d crustal strike slip fault model with depth variable properties for lower values of the state evolution distance of the friction law small events appear the nucleation phases of the small and large events are very similar suggesting that the size of an event is determined by the conditions on the fault segments the event is propagating into rather than by the nucleation process itself we demonstrate the importance of incorporating slow tectonic loading with elastodynamics by evaluating two simplified approaches one with the slow tectonic loading but no wave effects and the other with all dynamic effects included but much higher loading rate

earthquake supercycles in central italy inferred from 36cl exposure dating we use 36cl surface exposure dating to determine the slip release pattern over the last 15 kyr for the velino magnola fault a major active normal fault in central italy we sampled the fault at five well separated sites along its length and modeled the 36cl concentrations measured in the 376 samples we find that the fault broke in at least 9 large earthquakes that occurred in two 5 6 ka long supercycles each cycle included a 4 5 ka long phase of relative quiescence followed by a cluster of at least 3 large earthquakes or earthquake sequences that released most of the accumulated strain in 1 ka all 9 identified events broke the entire fault and produced maximum surface slips of 2 3 m though the velino magnola fault seems presently in a stage of relative quiescence it may re enter a phase of paroxysmal seismic activity in a few hundred of years

modeling afterslip and aftershocks following the 1992 landers earthquake one way to probe the rheology of the lithosphere and fault zones is to analyze the temporal evolution of deformation following a large earthquake in such a case the lithosphere responds to a known stress change that can be assessed from earthquake slip models constrained from seismology and geodesy here we model the postseismic response of a fault zone that is assumed to obey a rate strengthening rheology where the frictional stress varies as a ln e e being the deformation rate and a 0 a rheological parameter the model is simple enough that these parameters can be estimated by inversion of postseismic geodetic data we apply this approach to the analysis of geodetic displacements following the m w 7 3 1992 landers earthquake the model adjusts well the measured displacements and implies a 0 47 0 53 mpa in addition we show that aftershocks and afterslip follow the same temporal evolution and that the spatiotemporal distribution of aftershocks is consistent with the idea that they are driven by reloading of the seismogenic zone resulting from frictional afterslip