stress drop for three m 4 3 4 7 1992 1994 parkfield ca earthquakes in a previous paper fletcher and spudich 1998 we determined the slip distribution for the three m 4 3 to 4 7 parkfield events that occurred between oct 20 1992 and dec 20 1994 the rupture characteristics of these events are of particular interest because they are the largest events to occur at parkfield since the 1966 mainshock to further investigate the rupture process of these events we compute the static stress change or stress drop caused by the slip using a boundary integral method from quin and das 1989 adapted for the forward modeling case of computing stresses from a static slip distribution we find that the initial phase at or very near the hypocenter has the largest stress change in all cases this peak has a value of about 3 8 mpa for the oct 20 1992 event which also has the smallest magnitude of the three at 4 3 and about 45 mpa for the other nov 14 1993 and 60 mpa for the dec 20 1994 event the rest of the stress release is comparatively small and smooth for the nov 14 1993 event which is the deepest event and the most impulsive the pattern of stress release of the other two events oct 20 1992 and dec 20 1994 is more complex with a greater spatial extent the spatial distribution of stress is approximately similar to the spatial distribution of slip but with the larger initial peak in stress drop at the hypocenter

the seismic cycle in the area of the 2011 mw9 0 tohoku oki earthquake we model seismic and aseismic slip on the japan megathrust in the area of the m w9 0 2011 tohoku oki earthquake based on daily time series from 400 gps stations of the geonet network and campaign measurements of six sea floor displacements the coseismic and postseismic slip distributions are inverted simultaneously using principal component analysis based inversion method pcaim exploring a wide range of boundary conditions and regularization constraints we found the coseismic slip distribution to be quite compact with a peak slip between 30 and 50 m near the trench our model shows deep afterslip fringing the downdip edge of the coseismic rupture but also a dominant zone of shallow afterslip afterslip over the first 279 days following the main shock represents about 40 of the coseismic moment we compare the coseismic and postseismic models with an interseismic coupling model derived from inland and sea bottom measurements determined in a self consistent manner assuming that seismic and aseismic slip had to match the long term slip rate along the megathrust the recurrence time of m w9 0 earthquakes is estimated to 100 300 years while historical and paleotsunami records suggest a return period more of the order of 1000 years the discrepancy is smaller if the shallower portion of the megathrust is assumed to produce both aseismic slip as the afterslip model suggests and seismic slip during occasional large tsunamigenic earthquakes

control of rupture by fault geometry during the 1966 parkfield earthquake a reanalysis of the available data for the 1966 parkfleld california earthquake mr 5 suggests that although the ground breakage and aftershocks extended about 40 km along the san andreas fault the initial dynamic rupture was only 20 to 25 km in length the foreshocks and the point of initiation of the main event locate at a small bend in the mapped trace of the fault detailed analysis of the p wave first motions from these events at the gold hill station 20 km southeast indicates that the bend in the fault extends to depth and apparently represents a physical discontinuity on the fault plane other evidence suggests that this discontinuity plays an important part in the recurrence of similar magnitude 5 to 6 earthquakes at parkfield analysis of the strong motion records suggests that the rupture stopped at another discontinuity in the fault plane an en echeion offset near gold hill that lies at the boundary on the san andreas fault between the zone of aseismic slip and the locked zone on which the great 1857 earthquake occurred foreshocks to the 1857 earthquake occurred in this area sieh 1978 and the epicenter of the main shock may have coincided with the offset zone if it did a detailed study of the geological and geophysical character of the region might be rewarding in terms of understanding how and why great earthquakes initiate where they do

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

petroleum systems of the adiyaman region southeastern anatolia turkey abstract the petroliferous cretaceous sequence in the adiyaman region has been identified as a continental platform type system several source rocks exist within the sequence including shales mudstones and carbonates of the derdere karababa and karabogaz formations these source rocks contain mostly an oil prone type ii kerogen in the adiyaman region most oil is produced primarily from pores largely secondary in the karababa dolomitic wackestones and derdere dolomites oil is preserved in structural traps formed during the late cretaceous late miocene deformations and in stratigraphic traps prior to re migrating into late structural traps critical processes controlling the hydrocarbon accumulations occurred from the upper cenomanian lower turonian to aquitanian 93 22 ma the mass of petroleum expelled from the source rocks ranges from 4 0 to 10 3 mg hc g rock calculated fractional conversion and expulsion efficiency values of these source rocks have been found to be about 0 60 and 80 respectively the data presented in this study point to apparently low expelled efficiency values in the area since they also indicate that a significant amount of hydrocarbons have been lost due to exposure as the result of uplift or have escaped been expelled from the source rocks i e generated and are preserved somewhere awaiting discovery

the earthquake prediction experiment at parkfield california since 1985 a focused earthquake prediction experiment has been in progress along the san andreas fault near the town of parkfield in central california parkfield has experienced six moderate earthquakes since 1857 at average intervals of 22 years the most recent a magnitude 6 event in 1966 the probability of another moderate earthquake soon appears high but studies assigning it a 95 chance of occurring before 1993 now appear to have been over simplified the identification of a parkfield fault segment was initially based on geometric features in the surface trace of the san andreas fault but more recent microearthquake studies have demonstrated that those features do not extend to seismogenic depths on the other hand geodetic measurements are consistent with the existence of a locked patch on the fault beneath parkfield that has presently accumulated a slip deficit equal to the slip in the 1966 earthquake a magnitude 4 7 earthquake in october 1992 brought the parkfield experiment to its highest level of alert with a 72 hour public warning that there was a 37 chance of a magnitude 6 event however this warning proved to be a false alarm most data collected at parkfield indicate that strain is accumulating at a constant rate on this part of the san andreas fault but some interesting departures from this behavior have been recorded here we outline the scientific arguments bearing on when the next parkfield earthquake is likely to occur and summarize geophysical observations to date

block modeling with connected fault network geometries and a linear elastic coupling estimator in spherical coordinates geodetic observations of interseismic deformation provide constraints on the partitioning of fault slip across plate boundary zones the spatial distribution of both elastic and inelastic strain accumulation and the nature of the fault system evolution here we describe linear block theory which decomposes surface velocity fields into four components 1 plate rotations 2 elastic deformation from faults with kinematically consistent slip rates 3 elastic deformation from faults with spatially variable coupling and 4 homogeneous intrablock strain elastic deformation rates are computed for each fault segment in a homogeneous elastic half space using multiple optimal planar cartesian coordinate systems to minimize areal distortion and triangular dislocation elements to accurately represent complex fault system geometry block motions fault slip rates elastic coupling and internal block strain rates are determined simultaneously using a linear estimator with constraints from both geodetically determined velocity fields and geologic fault slip rate estimates we also introduce algorithms for efficiently implementing alternative fault network geometries to quantify parameter sensitivity to nonlinear perturbations in model geometry

asperities and barriers on the seismogenic zone in north chile state of the art after the 2007 mw 7 7 tocopilla earthquake inferred by gps and insar data the mw 7 7 2007 november 14 earthquake had an epicentre located close to the city of tocopilla at the southern end of a known seismic gap in north chile through modelling of global positioning system gps and radar interferometry insar data we show that this event ruptured the deeper part of the seismogenic interface 30 50 km and did not reach the surface the earthquake initiated at the hypocentre and was arrested 150 km south beneath the mejillones peninsula an area already identified as an important structural barrier between two segments of the peru chile subduction zone our preferred models for the tocopilla main shock show slip concentrated in two main asperities consistent with previous inversions of seismological data slip appears to have propagated towards relatively shallow depths at its southern extremity under the mejillones peninsula our analysis of post seismic deformation suggests that small but still significant post seismic slip occurred within the first 10 d after the main shock and that it was mostly concentrated at the southern end of the rupture the post seismic deformation occurring in this period represents 12 19 per cent of the coseismic deformation of which 30 55 per cent has been released aseismically post seismic slip appears to concentrate within regions that exhibit low coseismic slip suggesting that the afterslip distribution during the first month of the post seismic interval complements the coseismic slip the 2007 tocopilla earthquake released only 2 5 per cent of the moment deficit accumulated on the interface during the past 130 yr and may be regarded as a possible precursor of a larger subduction earthquake rupturing partially or completely the 500 km long north chile seismic gap

temporal changes in microseismicity and creep near parkfield california the 25 km long section of the san andreas fault near parkfield california ruptured in similar magnitude 6 earthquakes in 1881 1901 1922 1934 and 1966 on the basis of a number of lines of seismological evidence a section of the san andreas fault now termed the parkfield preparation zone has been identified as the locus of the next parkfield earthquake1 3 here we describe coincident changes in surface creep rates and deep seismicity near the parkfield preparation zone following the 2 may 1983 coalinga earthquake and suggest that both respond to the same stimuli these changes were concentrated near the point of initiation of the magnitude 6 characteristic parkfield earthquakes lending credence to the hypothesis that this section of the san andreas fault is characterized by a unique set of physical properties4 6 which make it ideal for earthquake prediction studies and which may also be useful for identifying hypocentral regions in other areas

subduction earthquake deformation associated with 14 november 2007 mw 7 8 tocopilla earthquake in chile results from insar and aftershocks abstract on 14 november 2007 a subduction thrust earthquake magnitude mw 7 8 occurred in the coastal region of northern chile causing substantial damage to the city of tocopilla we investigate the source fault of the earthquake slip distribution and fault interaction by integrating aftershock locations satellite interferometry data and stress model simulations aftershock measurements allow us to locate the area and geometry of the rupture plane in the coastal region between the cities of tocopilla and antofagasta combining two satellite viewing geometries acquired in envisat s wide swath and image modes we observe decimetre scale coseismic deformation the maximum line of sight displacement is found to be about 40 cm located at the mejillones peninsula slip inversions using elastic half space models with geometry constrained by aftershocks suggest rupturing of an area of 160 km by 50 km along the nazca south america convergent margin between latitudes 22 s and 23 5 s the main slip is concentrated on two asperities the largest being located in the southern part of the rupture area at a depth of approximately 30 50 km with a magnitude of about 2 5 m because aftershock distribution may also suggest a region of shallow crustal deformation activity located offshore we investigate whether the 2007 tocopilla earthquake also involved shallow crustal fault slip offshore although we find that the latter assumption is supported by coulomb stress modelling and geologic inferences our geodetic and seismic data provide insufficient constraints to resolve the exact geometry and kinematics of dislocation on this structure

spatiotemporal model of kīlauea s summit magmatic system inferred from insar time series and geometry free time dependent source inversion k lauea volcano hawai i island has a complex magmatic system including summit reservoirs and rift zones kinematic models of the summit reservoir have so far been limited to first order analytical solutions with pre determined geometry to explore the complex geometry and kinematics of the summit reservoir we apply a multitrack wavelet based insar interferometric synthetic aperture radar algorithm and a novel geometry free time dependent modeling scheme to map spatiotemporally distributed surface deformation signals over k lauea s summit we process sar data sets from two overlapping tracks of the envisat satellite including 100 images during the period 2003 2010 following validation against global positioning system gps data we invert the surface deformation time series to constrain the spatiotemporal evolution of the magmatic system without any prior knowledge of the source geometry the optimum model is characterized by a spheroidal and a tube like zone of volume change beneath the summit and the southwest rift zone at 2 3 km depth respectively to reduce the model dimension we apply a principal component analysis pca scheme which allows for the identification of independent reservoirs the first 3 pcs explaining 99 63 8 28 5 and 6 6 respectively of the model include six independent reservoirs with a complex interaction suggested by temporal analysis the data and model presented here in agreement with earlier studies improve the understanding of k lauea s plumbing system through enhancing the knowledge of temporally variable magma supply storage and transport beneath the summit and verify the link between summit magmatic activity seismicity and rift intrusions

automated fault model discretization for inversions for coseismic slip distributions 1 geoscientists increasingly rely on coseismic slip distributions inferred from geodetic observations to drive sophisticated models of the seismic cycle to date little work has been done on optimizing the parameterization of these fault models so that they reflect the resolving power of observed surface displacements the locations of noisy surface displacement observations are often widely scattered far from features we wish to analyze in the subsurface and result in highly variable model resolution with depth the few attempts to produce variably discretized fault planes are arduously constructed by hand and may not correctly reflect the ability of the data to resolve slip features motivated by the increasing size of geodetic data pools and the need for distributed slip models that accurately represent features the data can resolve we present a fully automated algorithm that iteratively adjusts the sizes of dislocations in a fault model we use the concept of smoothing scales derived from the model resolution matrix to resize dislocations so that each dislocation is sized appropriately given the area over which slip in that region of the fault would be smoothed we present a series of synthetic tests that utilize both sparse and dense data sets and we compare our variably discretized inversions to traditional regularly discretized inversions we also use our approach to invert for slip from geodetic observations of the 2004 mw 6 0 parkfield california earthquake and the 1995 mw 8 1 antofagasta chile earthquake

distribution of aseismic slip rate on the hayward fault inferred from seismic and geodetic data 1 we solve for the slip rate distribution on the hayward fault by performing a least squares inversion of geodetic and seismic data sets our analysis focuses on the northern 60 km of the fault interferometric synthetic aperture radar insar data from 13 independent ers interferograms are stacked to obtain range change rates from 1992 to 2000 horizontal surface displacement rates at 141 bench marks are measured using gps from 1994 to 2003 surface creep observations and estimates of deep slip rates determined from characteristic repeating earthquake sequences are also incorporated in the inversion the fault is discretized into 283 triangular dislocation elements that approximate the nonplanar attributes of the fault surface south of the city of hayward a steeply east dipping fault geometry accommodates the divergence of the surface trace and the microseismicity at depth the inferred slip rate distribution is consistent with a fault that creeps aseismically at a rate of 5 mm yr to a depth of 4 6 km the interferometric synthetic aperture radar insar data require an aseismic slip rate that approaches the geologic slip rate on the northernmost fault segment beneath point pinole although the insar data might be complicated by a small dip slip component at this location a low slip rate patch of 1 mm yr is inferred beneath san leandro consistent with the source location of the 1868 earthquake we calculate that the entire fault is accumulating a slip rate deficit equivalent to a mw 6 77 0 05 per century however this estimate of potential coseismic moment represents an upper bound because we do not know how much of the accumulated strain will be released through aseismic processes such as afterslip

geodetic and seismic constraints on some seismogenic zone processes in costa rica new seismic and geodetic data from costa rica provide insight into seismogenic zone processes in central america where the cocos and caribbean plates converge seismic data are from combined land and ocean bottom deployments in the nicoya peninsula in northern costa rica and near the osa peninsula in southern costa rica in nicoya inversion of gps data suggests two locked patches centered at 14 2 and 39 6 km depth interplate microseismicity is concentrated in the more freely slipping intermediate zone suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone the rupture zone for future large earthquakes at least over the short 1 year observation period we also estimate northwest motion of a coastal sliver block at 8 3 mm yr probably related to oblique convergence in the osa region to the south convergence is orthogonal to the trench cocos caribbean relative motion is partitioned here with 8 cm yr on the cocos panama block boundary including a component of permanent shortening across the fila costena fold and thrust belt and 1 cm yr on the panama block caribbean boundary the gps data suggest that the cocos plate panama block boundary is completely locked from 10 50 km depth this large locked zone as well as associated forearc and back arc deformation may be related to subduction of the shallow cocos ridge and or younger lithosphere compared to nicoya with consequent higher coupling and compressive stress in the direction of plate convergence

characteristics on fault coupling along the solomon megathrust based on gps observations from 2011 to 2014 the solomon megathrust along the western solomon arc generated two megathrust earthquakes in the past decade mw 8 1 in 2007 and mw 7 1 in 2010 to investigate the interseismic deformation and inferred coupling on the megathrust we deployed the first continuous gps network in the western solomon islands our 2011 2014 gps data and the back slip inversion model show coupling ratio as high as 73 along the southeastern 2007 rupture segment but only 10 on average along the segment of 2010 event based on the spatial distribution of coseismic slip aftershock clusters derived coupling pattern and paleo geodetic records we discovered the former as a semi permanent asperity and the latter as a potential megathrust barrier we propose that a characteristic earthquake of magnitude not less than mw 8 will recur in an interval of 100 or more years by either single or doublet earthquake

kinematic inversion of physically plausible earthquake source models obtained from dynamic rupture simulations one approach to investigate earthquake source processes is to produce kinematic source models from inversion of seismic records and geodetic data the setup of the inversion requires a variety of assumptions and constraints to restrict the range of possible models here we evaluate to what extent physically plausible earth quake scenarios are reliably restituted in spite of these restrictions we study which characteristics of ruptures such as rupture velocity slip distribution stress drop rise time and slip function can be reliably determined from the inversion of near field seismic and geodetic data using spontaneous dynamic rupture simulations we gen erate five earthquake scenarios each of which has different characteristics of the source process then we conduct a blind test by modeling the synthetic near source data using a standard inversion scheme that optimizes the fit to the observations while searching for solutions with minimum roughness the inversion procedure assumes a rupture front propagating away from the hypocenter with variable rupture velocity and a simple cosine slip time function our results show that overall slip distribution and stress drop are reasonably well determined even for input models with relatively complex histories such as a subshear rupture transitioning to supershear speeds depth averaged rupture velocities are also reasonably well resolved although their estimate progressively deteriorates away from the hypocenter the local rise time and slip function are not well resolved but there is some sensitivity to the rupture pulse width which can be used to differentiate between pulse like and crack like ruptures our test for understanding the inaccuracies in green s functions shows that random 3d perturbations of 5 standard deviation do not lead to significant degradation of the estimation of earthquake source parameters as remedies to the current limitations we propose smoothing slip function parameters and using more complicated inversion schemes only if data necessitates them online material figures showing snapshots of forward and inverse modeling of rupture l curves slip models and waveform fits

effects of layered crust on the coseismic slip inversion and related cff variations hints from the 2012 emilia romagna earthquake abstract the 2012 emilia romagna italy seismic sequence has been extensively studied given the occurrence of two mainshocks both temporally and spatially close to each other the recent literature accounts for several fault models obtained with different inversion methods and different datasets several authors investigated the possibility that the second event was triggered by the first mainshock with elusive results in this work we consider all the available insar and gps datasets and two planar fault geometries which are based on both seismological and geological constraints we account for a layered elastic half space hosting the dislocation and compare the slip distribution resulting from the inversion and the related changes in coulomb failure function cff obtained with both a homogeneous and layered half space finally we focus on the interaction between the two main events discriminating the contributions of coseismic and early postseismic slip of the mainshock on the generation of the second event and discuss the spatio temporal distribution of the seismic sequence when accounting for both insar and gps geodetic data we are able to reproduce a detailed coseismic slip distribution for the two mainshocks that is in accordance with the overall aftershock seismicity distribution furthermore we see that an elastic medium with depth dependent rigidity better accounts for the lack of the shallow seismicity amplifying with respect to the homogeneous case the mechanical interaction of the two mainshocks

detection of a locked zone at depth on the parkfield california segment of the san andreas fault the parkfield california segment of the san andreas fault is transitional in character between the creeping segment of the fault to the northwest and the locked carrizo plain segment to the southeast the rate of shallow fault slip decreases from 25 30 mm yr northwest of the epicenter of the 1966 parkfield earthquake to zero at the southeastern end of the 1966 rupture zone data from a network of trilateration lines spanning the san andreas fault near parkfield and extending to the pacific coast near san luis obispo shed light on the rate of fault slip at depth since the 1966 earthquake in this study average rates of line length change and shallow fault slip were inverted to determine the slip rate at depth on the parkfield fault segment the fault is taken to be a vertical surface with unknown distribution of strike slip displacement in an elastic half space a striking result of the inversions is that all solutions providing acceptable fits to the data exhibit a locked zone essentially coincident with the rupture surface of the 1966 parkfield earthquake the data require that the locked zone extend nearly as far north as the 1966 epicenter however the vertical extent of the locked zone is not well resolved over much of the parkfield segment the fault is slipping faster at the earth s surface than it is at seismogenic depths in order to fit the trilateration measurements it is necessary to include a component of contraction normal to the trend of the san andreas the inversion results suggest a spatially uniform normal strain of 0 06 strain yr the orientation of the contraction is compatible with geologic and seismic evidence of active folding and reverse faulting in the region the magnitude of the contraction is consistent with convergence rates inferred from global plate motion models

source parameters for the 1952 kern county earthquake california a joint inversion of leveling and triangulation observations coseismic leveling and triangulation observations are used to determine the faulting geometry and slip distribution of the july 21 1952 mw 7 3 kern county earthquake on the white wolf fault a singular value decomposition inversion is used to assess the ability of the geodetic network to resolve slip along a multisegment fault and shows that the network is sufficient to resolve slip along the surface rupture to a depth of 10 km below 10 km the network can only resolve dip slip near the fault ends the preferred source model is a two segment right stepping fault with a strike of 51 and a dip of 75 sw the epicentral patch has deep 6 27 km left lateral oblique slip while the northeastern patch has shallow 1 12 5 km reverse slip there is nearly uniform reverse slip epicentral 1 6 m northeast 1 9 m with 3 6 m of left lateral strike slip limited to the epicentral patch the seismic moment is m0 9 2 0 5 1019 n m mw 7 2 the signal to noise ratio of the leveling and triangulation data is reduced by 96 and 49 respectively the slip distribution from the preferred model matches regional geomorphic features and may provide a driving mechanism for regional shortening across the comanche thrust and structural continuity with the scodie seismic lineament to the northeast

slip budget and potential for a m7 earthquake in central california the slip rate budget of the san andreas fault saf in central california which is approximately 33 mm yr is accounted for by a change in the slip release mechanism along the fault in the nw section of the fault between bear valley and monarch peak creep apparently accounts for the slip budget with the seismicity contributing negligibly the section at parkfield marks the transition from a creeping to a locked fault trace since the m8 1857 earthquake five m6 earthquakes have occurred but have not completely accounted for the slip budget southeast of parkfield the saf has been locked since 1857 from cholame to bitterwater valley this section now lags the deep slip by the amount of slip released in 1857 consequently faulting in this section could occur at the time of the next parkfield earthquake if this earthquake releases the slip deficit accumulated in the transition zone and in the locked zone the earthquake will have a moment magnitude m7 2

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

insights into orogenic processes from drab schists and minor intrusions southern sao francisco craton brazil abstract minor altered intrusions and drab retrogressed schists can easily be overlooked in geological studies but this contribution explores these rocks within the archaean and palaeoproterozoic southern sao francisco craton ssfc brazil using geological relationships and accessory mineral in situ analyses in the context of cratonic assembly three magmatic pulses are documented i archaean and ii palaeoproterozoic felsic intrusions both hosted by archaean protoliths and iii palaeoproterozoic felsic intrusions in palaeoproterozoic supracrustal sequences archaean felsic intrusions confirm the palaeoarchaean age of the mafic ultramafic sequence of the rio das velhas greenstone belt and rhyacian intrusions mark the collisional stage of the mineiro belt with the sfc at c 2130 ma greenstone belt schists show a wide distribution of rounded soccer ball archaean detrital and metamorphic zircon grains ranging in age from 3200 to 2750 ma with an interpreted overprinting high grade metamorphic event at c 2700 to 2680 ma most high grade metamorphic rims have th u 0 1 negative ehf t values and ree pattern consistent with eclogite granulite metamorphic facies reinforcing the hypothesis of a dehydrated refractory crust formed during the stabilization of the ssfc even though no such protoliths are preserved this event links crustal thickening and partial melting of archaean lower crust archaean rutile crystals from the greenstone belt schist were reset during the palaeoproterozoic event but still preserve the early archaean high grade metamorphic signature the presence of unstable ilmenite replaced by rutile in the schist associated to felsic intrusions with the same age at c 2130 ma suggest high pressure low temperature prograde metamorphism during the collisional stage of the palaeoproterozoic orogen elongate and prismatic zircon grains from the rhyacian intrusions have low ehf t signature and crystallised from partial melting of sedimentary protoliths accretionary events produced thicker and more differentiated crust by the end of rhyacian time easily overlooked rocks in this study when studied have revealed a rich multi event history of cratonic evolution

frictional constraints on crustal faulting we consider how variations in fault frictional properties affect the phenomenology of earthquake faulting in particular we propose that lateral variations in fault friction produce the marked heterogeneity of slip observed in large earthquakes we model these variations using a rate and state dependent friction law where we differentiate velocity weakening behavior into two fields the strong seismic field is very velocity weakening and the weak seismic field is slightly velocity weakening similarly we differentiate velocity strengthening behavior into two fields the compliant field is slightly velocity strengthening and the viscous field is very velocity strengthening the strong seismic field comprises the seismic slip concentrations or asperities the two intermediate fields weak seismic and compliant have frictional velocity dependences that are close to velocity neutral these fields modulate both the tectonic loading and the dynamic rupture process during the interseismic period the weak seismic and compliant regions slip aseismically while the strong seismic regions remain locked evolving into stress concentrations that fail only in main shocks the weak seismic areas exhibit most of the interseismic activity and aftershocks but can also creep seismically this mixed frictional behavior can be obtained from a sufficiently heterogeneous distribution of the critical slip distance the model also provides a mechanism for rupture arrest dynamic rupture fronts decelerate as they penetrate into unloaded complaint or weak seismic areas producing broad areas of accelerated afterslip aftershocks occur on both the weak seismic and compliant areas around a fault but most of the stress is diffused through aseismic slip rapid afterslip on these peripheral areas can also produce aftershocks within the main shock rupture area by reloading weak fault areas that slipped in the main shock and then healed we test this frictional model by comparing the seismicity and the coseismic slip for the 1966 parkfield 1979 coyote lake and 1984 morgan hill earthquakes the interevent seismicity and aftershocks appear to occur on fault areas outside the regions of significant slip these regions are interpreted as either weak seismic or compliant depending on whether or not they manifest interevent seismicity

coseismic and postseismic slip of the 2004 parkfield earthquake from space geodetic data we invert interferometric synthetic aperture radar insar data jointly with campaign and continuous global positioning system gps data for slip in the coseismic and postseismic periods of the 2004 parkfield earthquake the insar dataset consists of eight interferograms from data collected by the envisat and radarsat satellites spanning the time of the earthquake and variable amounts of the postseismic period the two datasets complement each other with the insar providing dense sampling of motion in the range direction of the satellite and the gps providing more sparse but three dimensional measurements of ground motion the model assumes exponential decay of the postseismic slip with a decay time constant of 0 087 years determined from time series modeling of continuous gps and creepmeter data we find a geodetic moment magnitude of m 6 2 for a 1 day coseismic model and m w 6 1 for the entire postseismic period the coseismic rupture occurred mainly in two slip asperities one near the hypocenter and the other 15 20 km north postseismic slip occurred on the shallow portions of the fault and near the rupture areas of two m 5 0 aftershocks a comparison of the geodetic slip models with seismic moment estimates suggests that the coseismic moment release of the parkfield earthquake is as little as 25 of the total this underlines the importance of aseismic slip in the slip budget for the parkfield segment online material complete data tables and supplemental tables

the determination of displacement fields from geodetic data along a strike slip fault observations of angles or distances between stations of a geodetic network are commonly used to infer information about the movement of the surface of the earth the absence of any observations external to the network leads to an ambiguous displacement field existing techniques for eliminating this ambiguity are all unsatisfactory in some respect the best technique an inner coordinate solution is not appropriate for networks located in a strike slip fault environment the inner coordinate solution zeros the rotation of all stations about their center of mass along a strike slip fault like the san andreas however motion normal to the fault is less likely than motion parallel to the fault the solution presented here an outer coordinate solution finds the rotation of the network that minimizes the components of displacement normal to the fault since motion along a strike slip fault is generally expected to be parallel to the fault the displacements obtained with the outer coordinate solution are more reasonable than those obtained with other techniques examination of a trilateration network near san francisco bay demonstrates the large effect that the choice of adjustment technique can have on the inferred relative motion of the two sides of the fault the inner coordinate solution gave a rate of about 1 mm yr whereas the preferred outer coordinate solution rate was 36 mm yr

forecast model for moderate earthquakes near parkfield california earthquake instability models have possible application to earthquake forecasting because the models simulate both preseismic and coseismic changes of fault slip and ground deformation in the forecast procedure proposed here repeated measurements of preseismic fault slip and ground deformation constrain the values of model parameters the early part of the model simulation corresponds to the available field data and the subsequent part constitutes an estimate of future faulting and ground deformation in particular the time location and size of unstable faulting are estimates of the pending earthquake parameters the forecast accuracy depends on the model realism and parameter resolution the forecast procedure is applied to fault creep and trilateration data measured near parkfield california where at least five magnitude 5 5 to 6 earthquakes have occurred regularly since 1881 the last in 1966 the quasi static model consists of a flat vertical plane embedded in an elastic half space spacially variable fault slip of strike slip sense is driven by an increasing regional shear stress but is impeded by a relatively strong patch of brittle strain softening fault the field data are consistent with these approximate values of patch parameters radius of 3 km patch center 5 km deep and 8 km southeast of the 1966 epicenter and maximum brittle strength of 26 bars fluctuations in the available field data prevent estimating the earthquake time with any more precision than use of the 21 8 year recurrence interval however the model may later give a more precise estimate of the earthquake time if the fault slip rate near the inferred patch increases before the earthquake as predicted by the model

coseismic deformation revealed by inversion of strong motion and gps data the 2003 chengkung earthquake in eastern taiwan summary a moderate earthquake of mw 6 8 occurred on 2003 december 10 it ruptured the chihshang fault in eastern taiwan which is the most active segment of the longitudinal fault as a plate suture fault between the luzon arc of the philippine sea plate and the eurasian plate the largest coseismic displacements were 13 cm horizontal and 26 cm vertical we analyse 40 strong motion and 91 gps data to model the fault geometry and coseismic dislocations the most realistic shape of the chihshang fault surface is listric in type the dipping angle of the seismic zone is steep about 60 70 at depths shallower than 10 km and then gradually decreases to 40 50 at depths of 20 30 km thus the polygonal elements in poly3d are well suited for modelling complex surfaces with curving boundaries using the strong motion data the displacement reaches 1 2 m dip slip on the chihshang fault and decreases to 0 1 m near surface the slip averages 0 34 m releasing a scalar moment of 1 6e26 dyne cm for gps data our model reveals that the maximal dislocation is 1 8 m dip slip the dislocations decrease to 0 1 m near the surface the average slip is 0 48 m giving a scalar moment of 2 2e26 dyne cm regarding post seismic deformation a displacements of 0 5 m were observed near the chihshang fault indicating the strain had not been totally released as a probable result of near surface locking of the fault zone

earthquake potential along the northern hayward fault california the hayward fault slips in large earthquakes and by aseismic creep observed along its surface trace dislocation models of the surface deformation adjacent to the hayward fault measured with the global positioning system and interferometric synthetic aperture radar favor creep at 7 millimeters per year to the bottom of the seismogenic zone along a 20 kilometer long northern fault segment microearthquakes with the same waveform repeatedly occur at 4 to 10 kilometer depths and indicate deep creep at 5 to 7 millimeters per year the difference between current creep rates and the long term slip rate of 10 millimeters per year can be reconciled in a mechanical model of a freely slipping northern hayward fault adjacent to the locked 1868 earthquake rupture which broke the southern 40 to 50 kilometers of the fault the potential for a major independent earthquake of the northern hayward fault might be less than previously thought

preliminary determination of pacfic north america relative motion in the southern gulf of calfornia using the global positioning system global positioning system gps data from experiments conducted in 1985 and 1989 in the southern gulf of california mexico allow a determination of relative motion between the pacific and north american plates the data indicate motion of cabo san lucas on the pacific plate relative to north america at a rate of 47 7 mrn yr and azimuth of 57 6 west of north 1 errors equivalent within uncertainties to the nuvel 1 global plate motion model

recurrence models and parkfield california earthquakes main shocks of the earthquake sequences that occurred on the parkfield section of the san andreas fault in central california in 1922 1934 and 1966 are characterized by southeast rupture expansion over the same 20 to 30 km long section of the fault whereas the seismic moments for the 1922 and 1934 events are identical to within a precision of 10 the seismic moment for 1966 is 20 greater than for the earlier events to within a precision of 20 the parkfield area seismicity in general seems well described by a recurring moderate size characteristic earthquake repeating the same epicenter magnitude seismic moment rupture area and southeast direction of rupture expansion an unexplained 10 year advance of the 1934 event is the only discrepancy in the hypothesis that the parkfield earthquakes in 1857 1881 1901 1922 1934 and 1966 represent a strictly periodic process assuming the strictly periodic model and the absence since 1966 of the perturbations hypothesized for the 1922 to 1934 period the next characteristic parkfield earthquake should occur between 1983 and 1993

coulomb stress accumulation along the san andreas fault system 1 stress accumulation rates along the primary segments of the san andreas fault system are computed using a three dimensional 3 d elastic half space model with realistic fault geometry the model is developed in the fourier domain by solving for the response of an elastic half space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth this approach is then applied to the san andreas fault system using published slip rates along 18 major fault strands of the fault zone gps derived horizontal velocity measurements spanning the entire 1700 200 km region are then used to solve for apparent locking depth along each primary fault segment this simple model fits remarkably well 2 43 mm yr rms misfit although some discrepancies occur in the eastern california shear zone the model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates in addition shear and normal stresses along the major fault strands are used to compute coulomb stress accumulation rate as a result we find earthquake recurrence intervals along the san andreas fault system to be inversely proportional to coulomb stress accumulation rate in agreement with typical coseismic stress drops of 1 10 mpa this 3 d deformation model can ultimately be extended to include both time dependent forcing and viscoelastic response

change of apparent segmentation of the san andreas fault around parkfield from space geodetic observations across multiple periods sequences of earthquakes are commonly represented as a succession of periods of interseismic stress accumulation followed by coseismic and postseismic phases of stress release because the recurrence time of large earthquakes is often greater than the available span of space geodetic data it has been challenging to monitor the evolution of interseismic loading in its entire duration here we analyze large data sets of surface deformation at different key episodes around the cholame parkfield and creeping segments of the san andreas fault that show evidence of significant deceleration of fault slip during the interseismic period we compare the average fault slip rates before and after the 2004 mw6 parkfield earthquake in the 1986 2004 and 2006 2012 periods respectively avoiding 2 years of postseismic deformation after 2004 using a combination of gps data from the plate boundary observatory the southern california earthquake center crustal motion map and the bay area velocity unification networks and interferometric synthetic aperture radar from the advanced land observing satellite alos and envisat satellites we show that the area of coupling at the transition between the parkfield and cholame segments appears larger later in the interseismic period than it does earlier on while strong plate coupling is uniform across the parkfield and cholame segments in the 1986 2004 period creep occurs south of the 2004 epicenter after 2006 making segmentation of the san andreas fault south of parkfield more clearly apparent these observations indicate that analyses of surface deformation late in the earthquake cycle may overestimate the area of plate coupling a fault surface creeping much below plate rate may in some case be a region that does not promote earthquake nucleation but rather just be at a slower stage of its evolution our analysis also shows signs of large variation of slip velocity above and below plate rate in the creeping segment indicating that cycles of weakening and hardening can also be at play in dominantly aseismic areas

impact of megathrust geometry on inversion of coseismic slip from geodetic data application to the 1960 chile earthquake 1 we analyze the role of megathrust geometry on slip estimation using the 1960 chile earthquake mw 9 5 as an example a variable slip distribution for this earthquake has been derived by barrientos and ward 1990 applying an elastic dislocation model with a planar fault geometry their model shows slip patches at 80 110 km depth isolated from the seismogenic zone interpreted as aseismic slip we invert the same geodetic data set using a finite element model fem with precise geometry derived from geophysical data isoparametric fem is implemented to constrain the slip distribution of curve shaped elements slip resolved by our precise geometry model is limited to the shallow region of the plate interface suggesting that the deep patches of moment were most likely an artifact of the planar geometry our study emphasizes the importance of fault geometry on slip estimation of large earthquakes

seismicity and fault geometry of the san andreas fault around parkfield california and their implications abstract fault geometry is a consequence of tectonic evolution and it provides important information on potential seismic hazards we investigated fault geometry and its properties in parkfield california on the basis of local seismicity and seismic velocity residuals refined by an adaptive velocity hypocentral parameter inversion method the station correction terms from the hypocentral parameter inversion present characteristic seismic velocity changes around the fault suggesting low seismic velocities in the region east of the fault and high seismic velocities in the region to the west large seismic velocity anomalies are observed at shallow depths along the whole fault zone at depths of 3 8 km seismic velocity anomalies are small in the central fault zone but are large in the northern and southern fault zones at depths 8 km low seismic velocities are observed in the northern fault zone high seismicity is observed in the southwest fracture zone which has developed beside the creeping segment of the san andreas fault the vertical distribution of seismicity suggests that the fault has spiral geometry dipping ne in the northern region nearly vertical in the central region and sw in the southern region the rapid twisting of the fault plane occurs in a short distance of approximately 50 km the seismic velocity anomalies and fault geometry suggest location dependent piecewise faulting which may cause the periodic m6 events in the parkfield region

slip deficit on the san andreas fault at parkfield california as revealed by inversion of geodetic data a network of geodetic lines spanning the san andreas fault near the rupture zone of the 1966 parkfield california earthquake magnitude m 6 has been repeatedly surveyed since 1959 in the study reported here the average rates of line length change since 1966 were inverted to determine the distribution of interseismic slip rate on the fault these results indicate that the parkfield rupture surface has not slipped significantly since 1966 comparison of the geodetically determined seismic moment of the 1966 earthquake with the interseismic slip deficit rate suggests that the strain released by the latest shock will most likely be restored between 1984 and 1989 although this may not occur until 1995 these results lend independent support to the earlier forecast of an m 6 earthquake near parkfield within 5 years of 1988

plate coupling along the manila subduction zone between taiwan and northern luzon abstract we use gps data trench parallel gravity anomaly tpga and bathymetry to infer plate coupling patterns along the manila subduction zone using a block model and a fault geometry constrained by seismicity we simultaneously solve for the location of euler pole and angular velocity between the sunda and luzon blocks as well as the slip deficit rate on plate interface our estimates show that the euler pole between the sunda and luzon blocks is situated at southern palawan near 8 3 n and 119 4 e with the angular velocity of 4 6 myr 1 the estimated convergence rate along the manila trench continuously decreases southward from 91 mm yr at the northern tip of luzon to 55 mm yr north of mindoro the inversion of gps data reveals partially locked fault patches extending from the west luzon trough to the east of scarborough seamount chain the slip deficit rate in this region is in the range of 20 30 mm yr corresponding to a coupling ratio of 0 4 however the fault slip behavior is not well resolved near the north luzon trough based on a good correlation between locations of large subduction zone earthquakes and areas possessing gravity low we investigate a variety of tpga based plate coupling models assuming different scaling between tpga values and plate coupling ratios the tpga based plate coupling models offer plausible rupture scenarios which are not constrained by current gps data the partially locked fault zone near 15 16 5 n may be associated with the subducted scarborough seamount wherein oceanic floor is highly fractured the great subduction zone earthquake propagates beneath the scarborough seamount seems to be unlikely the densification of gps network in central luzon and seafloor geodetic observations close to trench axis are crucial to distinguish the detailed fault coupling patterns and the role of subducted seamounts

long term fault creep observations in central california the u s geological survey usgs has been monitoring aseismic fault slip in central california for more than 10 years as part of an earthquake prediction experiment since 1968 the usgs creep network has grown from one creep meter at the cienega winery south of hollister to a 44 station network that stretches from hayward east of san francisco bay to palmdale in southern california in general the long term slip pattern is most variable on sections of the faults where several magnitude 4 and larger earthquakes occurred during the recording period e g calaveras fault near hollister and san andreas fault between san juan bautista and bear valley these fault sections are the most difficult to characterize with a single long term slip rate in contrast sections of the faults that are seismically relatively quiet e g san andreas fault between bear valley and parkfield produce the steadiest creep records and are easiest to fit with a single long term slip rate appendix is available with entire article on microfiche order from the american geophysical union 2000 florida avenue n w washington d c 20009 document j82 004 1 00 payment must accompany order

surface deformation associated with the 2015 m w 8 3 illapel earthquake revealed by satellite based geodetic observations and its implications for the seismic cycle abstract in this study we present inter co and post seismic displacements observed in the 2015 illapel earthquake area by global positioning system gps and synthetic aperture radar interferometry insar radarsat 2 alos 2 and sentinel 1a interferograms capture the co and post seismic displacements due to the illapel earthquake based on a layered earth structure we modeled both co and post seismic faulting behaviors on the subduction interface of central chile the best fit model shows that the coseismic rupture broke a 200 km 200 km area with a maximum slip of 10 m at a depth of 20 km two distinct slip centers likely controlled by local ramp flat structure are revealed the total coseismic geodetic moment is 2 76 10 21 n m equivalent to a moment magnitude 8 3 the accumulated afterslip in the first two months after the mainshock is observed on both sides of the coseismic rupture zone with both ascending and descending sentinel 1a interferograms a limited overlap zone between co and post seismic slip models can be observed suggesting partitioning of the frictional properties within the illapel earthquake rupture zone the total afterslip releases 5 0 10 20 n m geodetic moment which is equivalent to an earthquake of m w 7 7 the 2010 m w 8 8 maule earthquake that occurred 400 km away from the illapel earthquake epicenter could have exerted certain effects on the seismic cycle of the illapel earthquake area the seismic records from 2000 to 2015 imply that the rate of annual seismic moment release in the illapel earthquake area dropped from 0 4 to 0 2 10 19 n m yr after the maule earthquake based on the forward modeling with the best fit slip models determined in this study we reproduce the local surface displacements before during and after the illapel earthquake a rough deformation cycle 105 29 yr calculated by using the coseismic displacements and interseismic rate is basically identical with the revisit interval of m8 events in the adjacent areas of the illapel earthquake suggesting that elastic rebound theory is applicable for the long term prediction in this region

the 2016 m7 kumamoto japan earthquake slip field derived from a joint inversion of differential lidar topography optical correlation and insar surface displacements

deep and shallow sources for the lusi mud eruption revealed by surface deformation the lusi mud eruption in east java indonesia began in may 2006 and continues to the present previous analyses of surface deformation data suggested an exponential decay of the pressure in the mud source but did not constrain the location geometry and evolution of the possible source s of the erupting mud and fluids to map the surface deformation we employ multitemporal interferometric synthetic aperture radar and analyze a well populated l band data set acquired by the advanced land observing satellite alos between may 2006 and april 2011 we then apply a time dependent inverse modeling scheme volume changes occur in two regions beneath lusi at 0 3 2 0 km and 3 5 4 75 km depth the cumulative volume change within the shallow source is 2 3 times larger than that of the deep source the observation and model suggest that a shallow source plays a key role by supplying the erupting mud but that additional fluids do ascend from depths 4 km on eruptive timescales

extent and duration of the 2003 cascadia slow earthquake 1 inversion of continuous gps measurements from the pacific northwest show the 2003 cascadia slow earthquake to be among the largest of ten transients recognized here twelve stations bracketing slow slip indicate transient slip propagated bi directionally from initiation in the southern puget basin reaching 300 km along strike over a period of seven weeks this event produced for the first time resolvable vertical subsidence and horizontal displacement reaching six mm in southern washington state inverted for non negative thrust slip a maximum of 3 8 cm of slip is inferred centered at 28 km depth near the sharp arch in the subducting juan de fuca plate nearly all slip lies shallower than 38 km inverted slip shows a total moment release equal to mw 6 6 and a high degree of spatial localization rather than near uniform slip this suggests rupture concentrated along asperities holds for slow earthquakes as well as conventional events

geologic co2 input into groundwater and the atmosphere soda springs id usa abstract a set of co 2 flux geochemical and hydrologic measurement techniques was used to characterize the source of and quantify gaseous and dissolved co 2 discharges from the area of soda springs southeastern idaho an eddy covariance system was deployed for one month near a bubbling spring and measured net co 2 fluxes from 74 to 1147 g m 2 d 1 an inversion of measured eddy covariance co 2 fluxes and corresponding modeled source weight functions mapped the surface co 2 flux distribution within and quantified co 2 emission rate 24 9 t d 1 from a 0 05 km 2 area surrounding the spring soil co 2 fluxes 2 d 1 were measured within a 0 05 km 2 area of diffuse degassing using the accumulation chamber method the estimated co 2 emission rate from this area was 49 t d 1 a carbon mass balance approach was used to estimate dissolved co 2 discharges from contributing sources at nine springs and the soda springs geyser total dissolved inorganic carbon as co 2 discharge for all sampled groundwater features was 57 1 t d 1 of this quantity approximately 3 was derived from biogenic carbon dissolved in infiltrating groundwater 35 was derived from carbonate mineral dissolution within the aquifer s and 62 was derived from deep source s isotopic compositions of helium 1 74 2 37 r a and deeply derived carbon 13 c 3 suggested contribution of volatiles from mantle and carbonate sources assuming that the deeply derived co 2 discharge estimated for sampled groundwater features 35 t d 1 is representative of springs throughout the study area the total rate of deeply derived co 2 input into the groundwater system within this area could be 350 t d 1 similar to co 2 emission rates from a number of quiescent volcanoes

the co seismic slip distribution of the landers earthquake we derived a model for the co seismic slip distribution on the faults which ruptured during the landers earthquake sequence of 28 june 1992 the model is based on the inversion of surface geodetic measurements primarily vector displacements measured using the global positioning system gps the inversion procedure assumes that the slip distribution is to some extent smooth and purely right lateral strike slip for a given fault geometry a family of so lutions of varying smoothness can be generated we choose the optimal model from this family based on cross validation which measures the predictive power of the data and the trade off of misfit and roughness solutions which give roughly equal weight to misfit and smoothness are preferred and have certain features in common 1 there are two main patches of slip on the johnson valley fault and on the homestead valley emerson and camp rock faults 2 virtually all slip is in the upper 10 to 12 km and 3 the model reproduces the general features of the geologically measured sur face displacements without prior constraints on the surface slip in all models regardless of smoothing very little slip is required on the fault that represents the big bear event and the total moment of the landers event is 9 1019 n m the nearly simultaneous rupture of multiple distinct faults suggests that much of the crust in this region must have been close to failure prior to the earthquake

reverse drag revisited why footwall deformation may be the key to inferring listric fault geometry abstract although reverse drag the down warping of hanging wall strata toward a normal fault is widely accepted as an indicator of listric fault geometry previous studies have shown that similar folding may form in response to slip on faults of finite vertical extent with listric or planar geometry in this study we therefore seek more general criteria for inferring subsurface fault geometry from observations of near surface deformation by directly comparing patterns of displacement stress and strain around planar and listric faults as predicted by elastic boundary element models in agreement with previous work we find that models with finite planar planar detached and listric detached faults all develop hanging wall reverse drag folds all of these model geometries also predict a region of tension and elevated maximum coulomb stress in the hanging wall suggesting that the distribution and orientation of near surface joints and secondary faults may also be of limited utility in predicting subsurface fault geometry the most notable difference between the three models however is the magnitude of footwall uplift footwall uplift decreases slightly with introduction of a detachment and more significantly with the addition of a listric fault shape a parametric investigation of faults with constant slip ranging from nearly planar to strongly listric over depths from 1 to 15 km reveals that footwall fold width is sensitive to fault geometry while hanging wall fold width largely reflects fault depth application of a graphical approach based on these results as well as more complete inverse modeling illustrates how patterns of combined hanging wall and footwall deformation may be used to constrain subsurface fault geometry

an interpretation of episodic slip on the calaveras fault near hollister california data from a precision multi wavelength distance measuring mwdm instrument located near hollister california have been analyzed in terms of strike slip faulting in the region covered by the network four episodes of deformation that are readily identifiable in the mwdm data for the year following september 1975 can be modeled as slip on the calaveras fault creepmeter data tend to support this interpretation since the creep events and the observed episodes of deformation appear to be correlated in a temporal sense lower bounds on the moment for each episode of slip can be determined from the mwdm data by the technique of linear programing the results indicate that aseismic slip is the dominant mechanism for strain release since the combined moment of 1 16 1024 dyn cm for the four episodes far exceeds the moment of all of the earthquakes that occurred during the year creepmeter data and geologic evidence that are taken in conjunction with the lower bound of the moment indicate that the depth of aseismic slip for some of the episodes of fault slip could extend to below the greatest depth of earthquake hypocenters for this region of central california

interactions between the landers and hector mine california earthquakes from space geodesy boundary element modeling and time dependent friction the m w 7 1 hector mine earthquake happened 7 years after the 1992 m w 7 3 landers earthquake within a system of right lateral strike slip faults comprising the eastern california shear zone because estimates of recurrence intervals of m 7 earthquakes on these faults range anywhere from 1500 to 50 000 years the close temporal spacing of these two earthquakes suggests that there is some interaction between the two current models of simple coulomb static stress interactions between the earthquake ruptures do not predict an obvious cause and effect relationship the interaction between the computed normal and right lateral shear stress reductions induced by the landers earthquake at the hector mine hypocenter is unclear we use a combination of space geodesy boundary element modeling and computer modeling of time dependent fault friction to investigate the interaction between the two earthquakes we compare stress changes on the hector mine rupture plane induced by the landers earthquake with the detailed slip distribution inferred using a combination of gps and insar data the slip distributions of both earthquakes are also used to infer the magnitude of the shear stress drop on each earthquake rupture and the orientation of the remote background stress consistent with these two recent events for each earthquake the azimuth of the remote maximum compressive stress was approximately 17 6 and the magnitude of the shear stress drop was 8 1 mpa on the landers rupture and 10 2 mpa on the hector mine rupture using a simple spring and slider model as a proxy for the rate and state frictional response of the hector mine faults to landers earthquake imposed stress steps we find that while a decrease in normal stress at the hector mine hypocenter would by itself have caused the hector mine earthquake to nucleate immediately the simultaneous decrease in shear stress may have caused a delay in the peak shear loading consideration of an acceptable range of rate and state friction parameters secular strain rates and crustal stiffnesses leads us to conclude that this delay could have been 0 to 40 years

locking depths estimated from geodesy and seismology along the san andreas fault system implications for seismic moment release 1 the depth of the seismogenic zone is a critical parameter for earthquake hazard models independent observations from seismology and geodesy can provide insight into the depths of faulting but these depths do not always agree here we inspect variations in fault depths of 12 segments of the southern san andreas fault system derived from over 1000 gps velocities and 66 000 relocated earthquake hypocenters geodetically determined locking depths range from 6 to 22 km while seismogenic thicknesses are largely limited to depths of 11 20km theseseismogenicdepthsbestmatchthegeodeticlockingdepthswhenestimated at the 95 cutoff depth in seismicity and most fault segment depths agree to within 2 km however the imperial coyote creek and borrego segments have significant discrepancies in these cases the geodetically inferred locking depths are much shallower than the seismogenic depths we also examine variations in seismic moment accumulation rate per unit fault length as suggested by seismicity and geodesy and find that both approaches yield high rates 1 5 1 8 1 0 13 nm yr km along the mojave and carrizo segments and low rates 0 2 10 13 nm yr km along several san jacinto segments the largest difference in seismic moment between models is calculated for the imperial segment where the moment rate from seismic depths is a factor of 2 5 larger than that from geodetic depths such variability has important implications for the accuracy to which future major earthquake magnitudes can be estimated

deformation near the junction of the creeping and locked segments of the san andreas fault cholame valley california 1970 1980 we have constructed a simple dislocation model of the san andreas fault near cholame california that reproduces the observed deformation of nearby geodetic networks as well as the transition from quasi continuous slip fault creep to a no slip locked condition on the surface trace the model is specified by the slip rates imposed on the fault surface slip on the surface trace is given by the observed fault creep rates and slip at depths greater than 12 6 km is taken to 30 mm yr along the entire length of the fault i e beneath creeping and locked sections alike slip at intermediate depths is assigned by a rather arbitrary extrapolation of the surface creep values this postulated distribution of slip on the fault is tested by using it in a dislocation model to predict the deformation in two geodetic networks a trilateration network 20 km aperture and a 20 km long level line the calculated deformation of the trilateration network agrees with the observed deformation within the observational error agreement between the predicted and observed elevation changes along the level line was satisfactory in the sense that the predicted values were below the level of detection and the observed changes were not significant

application of real time gps to earthquake early warning 1 we explore the use of real time high rate gps displacement data for earthquake early warning using 1 hz displacement waveforms from the april 4 2010 mw 7 2 el mayor cucapah earthquake we compare these data to those provided by the broadband velocity and accelerometer instrumentation of the southern california seismic network the unique information provided by the gps based displacement timeseries is the permanent static displacement using a simple algorithm that can be applied in real time we extract the static offset shortly after the s wave arrival around the time of the observed peak shaking at the same site and before shaking at more distant locations these data can be used as they become available to provide a robust estimate of the earthquake magnitude which ranges from 6 8 to 7 0 in this case we therefore conclude that real time high rate gps can provide a useful and independent assessment of earthquake magnitude for the purpose of earthquake early warning and real time earthquake information systems in general including tsunami warning systems

recent geodetic unrest at santorini caldera greece 1 after approximately 60 years of seismic quiescence within santorini caldera in january 2011 the volcano reawakened with a significant seismic swarm and rapidly expanding radial deformation the deformation is imaged by a dense network of 19 survey and 5 continuous gps stations showing that as of 21 january 2012 the volcano has extended laterally from a point inside the northern segment of the caldera by about 140 mm and is expanding at 180 mm yr a series of spherical source models show the source is not migrating significantly but remains about 4 km depth and has expanded by 14 million m3since inflation began a distributed sill model is also tested which shows a possible n s elongation of the volumetric source while observations of the current deformation sequence are unprecedented at santorini it is not certain that an eruption is imminent as other similar calderas have experienced comparable activity without eruption

assessment of stress coupling among the inter co and post seismic phases related to the 2004 m6 parkfield earthquake due to large uncertainties and non uniqueness in fault slip inversion the investigation of stress coupling based on the direct comparison of independent slip inversions for example between the coseismic slip distribution and the interseismic slip deficit may lead to ambiguous conclusions in this study we therefore adopt the stress constrained joint inversion in the bayesian approach of wang et al and implement the physical hypothesis of stress coupling as a prior we test the hypothesis that interseismic locking is coupled with the coseismic rupture and the early post seismic deformation is a stress relaxation process in response to the coseismic stress perturbation we characterize the role of stress coupling in the seismic cycle by evaluating the efficiency of the model to explain the available data taking the 2004 m6 parkfield earthquake as a study case we find that the stress coupling hypothesis is in agreement with the data the coseismic rupture zone is found to be strongly locked during the interseismic phase and the post seismic slip zone is indicated to be weakly creeping the post seismic deformation plays an important role to rebuild stress in the coseismic rupture zone based on our results for the stress accumulation during both inter and post seismic phase in the coseismic rupture zone together with the coseismic stress drop we estimate a recurrence time of m6 earthquake in parkfield around 23 41 yr suggesting that the duration of 38 yr between the two recent m6 events in parkfield is not a surprise

subsidence derived volumetric strain models for mapping extensional fissures and constraining rock mechanical properties in the san joaquin valley california large scale subsidence due to aquifer overdraft is an ongoing hazard in the san joaquin valley subsidence continues to cause damage to infrastructure and increases the risk of extensional fissures here we use insar derived vertical land motion vlm to model the volumetric strain rate due to groundwater storage change during the 2007 2010 drought in the san joaquin valley central valley california we then use this volumetric strain rate model to calculate surface tensile stress in order to predict regions that are at the highest risk for hazardous tensile surface fissures we find a maximum volumetric strain rate of 232 microstrain yr at a depth of 0 to 200 m in tulare and kings county california the highest risk of tensile fissure development occurs at the periphery of the largest subsiding zones particularly in tulare county and merced county finally we assume that subsidence is likely due to aquifer pressure change which is calculated using groundwater level changes observed at 300 wells during this drought we combine pressure data from selected wells with our volumetric strain maps to estimate the quasi static bulk modulus k a poroelastic parameter applicable when pressure change within the aquifer is inducing volumetric strain this parameter is reflective of a slow deformation process and is one to two orders of magnitude lower than typical values for the bulk modulus found using seismic velocity data the results of this study highlight the importance of large scale high resolution vlm measurements in evaluating aquifer system dynamics hazards associated with overdraft and in estimating important poroelastic parameters

eddy covariance imaging of diffuse volcanic co2 emissions at mammoth mountain ca usa use of eddy covariance ec techniques to map the spatial distribution of diffuse volcanic co2 fluxes and quantify co2 emission rate was tested at the horseshoe lake tree kill area on mammoth mountain california usa ec measurements of co2 flux were made during september october 2010 and ranged from 85 to 1 766 g m 2 day 1 comparative maps of soil co2 flux were simulated and co2 emission rates estimated from three accumulation chamber ac co2 flux surveys least squares inversion of measured eddy covariance co2 fluxes and corresponding modeled source weight functions recovered 58 77 of the co2 emission rates estimated based on simulated ac soil co2 fluxes spatial distributions of modeled surface co2 fluxes based on ec and ac observations showed moderate to good correspondence r2 0 36 to 0 70 results provide a framework for automated monitoring of volcanic co2 emissions over relatively large areas

estimation of interplate coupling in the nankai trough japan using gps data from 1996 to 2006 we used three component surface velocities in southwest japan to estimate plate coupling on the subducting plate interface at the nankai trough we analyzed continuous gps data from the japanese geonet network from 1996 to 2006 using a consistent analysis strategy that generates bias fixed solutions for the entire network we applied systematic time series analysis methods to estimate common mode error which improved position solutions for the entire network to allow for differences in regional deformation sources we modelled the plate coupling on the plate interface beneath shikoku island to kii peninsula and the tokai suruga trough separately the results show strong coupling at a depth of 10 30 km off shikoku and kii peninsula the spatial variation in plate coupling coincides well with the coseismic rupture zones of the past large earthquakes maximum slip deficit rates of 2 3 cm yr 1 at the depth of 5 25 km are found beneath the tokai area consistent with results from other studies the downdip limits of the highly coupled areas and transition zones beneath shikoku and the kii peninsula correspond approximately to estimates of the 450 c isotherms good correlation is observed between the lateral variations of the slip deficit distribution low frequency earthquakes and coseismic slip this correlation suggests that temperature and possibly fluid variations contribute to such correlation in space the interplate slip deficit derived from the gps velocities over the 10 yr of observations is generally compatible with the results over shorter time spans suggesting that plate coupling in sw japan does not change significantly over the period of these gps measurements

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

geodetic measurements near parkfield california 1959 1984 trilateration measurements have been made by various agencies since 1959 near parkfield california in the transition zone between the central creeping section and the southern locked section of the san andreas fault within the transition zone there are several creep meters and alignment arrays four small aperture trilateration networks and a broad trilateration network extending 70 km southwest from the san andreas fault to the pacific coast line length changes observed after the 1966 ml 5 6 parkfield earthquake indicate a geodetic moment of 4 5 1 0 1025 dyn cm the data are consistent with either 91 8 cm slip on a rupture zone which is 30 km long and extends from 3 to 8 km below the surface or 59 5 cm on a 30 km long zone 2 10 km deep this moment estimate includes some postseismic slip and is a factor of 2 4 greater than the moment determined from seismic surface waves the post 1966 data set is well suited to discriminate between the effects of surface and deep slip surface slip rates determined by creep meters alinement arrays and small aperture trilateration networks decrease from 29 mm yr at slack canyon to zero south of highway 46 a gradient of roughly 0 48 0 01 mm yr km deformation of a subset of the trilateration network between cholame and san luis obispo is consistent with shallow slip of 2 9 1 2 mm yr from the surface to a depth of 16 km and 33 4 5 5 mm yr below 16 km the trilateration data alone do not uniquely constrain both the transition depth from shallow to deep slip and the deep slip rate a transition depth near 16 km is preferred because it provides the best fit to the data and because the predicted slip rate is consistent with independent estimates of the total slip rate across the san andreas fault this simple two dimensional model neglects along strike variations in the deformation the data show no evidence for significant strike slip deformation across faults between the san andreas and the pacific coast

an increment of volcano collapse kinematics of the 1975 kalapana hawaii earthquake abstract earthquakes on kilauea volcano can be even more hazardous than volcanic eruptions the most recent large earthquake was the 1975 kalapana quake which generated a locally damaging tsunami over 8 m of horizontal surface displacement and 3 m of subsidence at the coast the seismic magnitude estimates for this earthquake range from ms 7 2 to mw 7 7 the complexity of the seismic data is matched by the complex surface deformation that was observed with a combination of leveling tilt and electronic distance measurements edm the geodetic data shows evidence of collapse of the caldera due to magma withdrawal localized subsidence over the rift zone normal faulting associated with shallow slumping along the hilina pali and displacement on the basal detachment fault we model the summit caldera rift zones and basal detachment fault as dislocations in an elastic half space to better quantify the sources of deformation related to this devastating event using inversion techniques that allow us to let the data constrain the geometry and magnitude of these sources we find 0 04 km3 of magma withdrawal at the summit between 3 5 m of opening along the rift zone and 7 1 m of slip along the basal detachment at 8 3 km depth models that allow for finer spatial resolution of slip on the detachment show that the largest slip occurred west of the seismically recorded earthquake hypocenter along the region of the fault below the coast and that the majority of the fault slip occurred south of the region of microseismicity our best fitting model has a geodetic moment of 4 1 1020 nm mw 7 7 which is consistent with tsunami models and recent analysis of long period seismic data the residual displacements of sites located in the hanging wall of the hilina pali slump system which were not used in the inversion of the other sources suggest that triggered shallow slumping contributed several meters of horizontal displacement to the coseismic displacement at these sites

slip patterns and earthquake populations along different classes of faults in elastic solids numerical simulations of slip instabilities on a vertical strike slip fault in an elastic half space are performed for various models belonging to two different categories the first category consists of inherently discrete cellular fault models such are used to represent fault systems made of segments modeled by numerical cells that can fail independently of one another their quasi independence is assumed to provide an approximate representation of strong fault heterogeneity due to geometric or material property disorder that can arrest ruptures at segment boundaries the second category consists of models having a well defined continuum limit these involve a fault governed by rate and state dependent friction and are used to evaluate what types of property heterogeneity could lead to the quasi independent behavior of neighboring fault segments assumed in the first category the cases examined include models of a cellular fault subjected to various complex spatial distributions of static to kinetic strength drops and models incorporating rate and state dependent friction subjected to various spatial distributions of effective stress normal stress minus pore pressure the results indicate that gradual effective stress variations do not provide a sufficient mechanism for the generation of observed seismic response strong and abrupt fault heterogeneity as envisioned in the inherently discrete category is required for the generation of complex slip patterns and a wide spectrum of event sizes strong fault heterogeneity also facilitates the generation of rough rupture fronts capable of radiating high frequency seismic waves the large earthquakes in both categories of models occur on a quasi periodic basis the degree of periodicity increases with event size and decreases with model complexity however in all discrete segmented cases the models generate nonrepeating sequences of earthquakes and the nature of the large quasi periodic events is highly variable the results indicate that expectations for regular sequences of earthquakes and or simple repetitive precursory slip patterns are unrealistic the frequency size fs statistics of the small failure episodes simulated by the cellular fault models are approximately self similar with b 1 2 and ba 1 where b and ba are b values based on magnitude and rupture area respectively for failure episodes larger than a critical size however the simulated statistics are strongly enhanced with respect to self similar distributions defined by the small events this is due to the fact that the stress concentrated at the edge of a rupture expanding in an elastic solid grows with the rupture size when the fault properties e g geometric irregularities are characterized by a narrow range of size scales the scaling of stress concentrations with the size of the failure zone creates a critical rupture area terminating the self similar earthquake statistics in such systems events reaching the critical size become on the average unstoppable and they continue to grow to a size limited by a characteristic model dimension when however the system is characterized by a broad spectrum of size scales the above phenomena are suppressed and the range of apparent self similar fs statistics is broad and characterized by average b and ba values of about 1 the simulations indicate that power law extrapolations of low magnitude seismicity will often underestimate the rate of occurrence of moderate and large earthquakes the models establish connections between features of fs statistics of earthquakes range of self similar regimes local maxima and structural properties of faults dominant size scales of heterogeneities dimensions of coherent brittle zones the results suggest that observed fs statistics can be used to obtain information on crustal thickness and fault zone structure

inferring fault characteristics using fold geometry constrained by airborne laser swath mapping at raplee ridge utah 1 we used airborne laser swath mapping alsm data to define the geometry of strata in raplee ridge in southeastern utah that were likely deformed by slip on an underlying fault the fold geometry was used to infer the fault geometry and loading responsible for the folding using the boundary element model poly3d monte carlo markov chain mcmc methods created probability density functions pdfs that reveal the most likely fault geometry and loading conditions and their uncertainties despite the restrictive assumptions of poly3d in which body deformations occur from slip on displacement discontinuities in a homogeneous linear elastic half space the strata deflections are matched remarkably well and the resulting inferences of fault geometry and stress field appear reasonable we propose that similar methods may be used with models that include more complicated rheologies and heterogeneous and anisotropic material properties to estimate fault geometries and loading conditions where only deflection of strata record these parameters

earthquake failure sequences along a cellular fault zone in a three dimensional elastic solid containing asperity and nonasperity regions numerical simulations of earthquake failure sequences along a discrete cellular fault zone are performed for a three dimensional 3 d model representing approximately the central san andreas fault the model consists of an upper crust overlying a lower crust and mantle region together defining an elastic half space with a vertical half plane fault the fault contains a region where slip is calculated on a uniform grid of cells governed by a static kinetic friction law and regions where slip is prescribed so as to represent tectonic loading aseismic fault creep and adjacent great earthquakes the computational region models a 70 km long and 17 5 km deep section of the san andreas fault to the nw of the great 1857 rupture zone different distributions of stress drops on failing computational cells are used to model asperity parkfield asperity and nonasperity fault regions the model is inherently discrete and corresponds to a situation in which a characteristic size of geometric disorder within the fault i e cell size here a few hundreds of meters is much larger than the nucleation size of the order of tens of centimeters to tens of meters based on slip weakening or state evolution slip distances the computational grid is loaded by a constant plate motion imposed at the lower crust upper mantle and creeping fault regions and by a staircase slip history imposed at the 1857 and 1906 rupture zones stress transfer along and outside the fault due to the imposed loadings and failure episodes along the computational grid is calculated using 3 d elastic dislocation theory the resulting displacement field in the computational region is compatible with geodetic and seismological observations only when the asperity and nonasperity regions are characterized by significantly different average stress drops the frequency magnitude statistics of the simulated failure episodes are approximately self similar for small events with b 1 2 the b value of statistics based on rupture area ba is about 1 but are strongly enhanced with respect to self similarity for events larger than a critical size this is interpreted as a direct manifestation of our 3 d elastic stress transfer calculations beyond certain rupture area and potency seismic moment divided by rigidity release values the event is usually unstoppable and it continues to grow to a size limited by a characteristic model dimension this effect is not accounted for by cellular automata and block spring models in which the adopted simplified stress transfer laws fail to scale properly with increasing rupture size the simulations suggest that local maxima in observed frequency magnitude statistics correspond to dimensions of coherent brittle zones such as the width of the seismogenic layer or the length of a fault segment bounded by barriers the analysis indicates that a single cell size representing approximately a single scale of geometric disorder cannot induce self similarity in a 3 d elastic model over a broad range of magnitudes a representation of geometric disorder covering a range of scales may thus be required to generate a wide domain of self similar gutenberg richter statistics our simulations show a great diversity in the mode of failure of the parkfield asperity the earthquakes themselves define an irregular sequence of events the modeling like many other discrete fault models suggests that expectations for periodic parkfield earthquakes and or simple precursory patterns repeating from one event to the other are unrealistic

lithological control on the deformation mechanism and the mode of fault slip on the longitudinal valley fault taiwan the longitudinal valley fault lvf in taiwan is creeping at shallow depth along its southern half where it is bounded by the lichi melange by contrast the northern segment of the lvf is locked where it is bounded by forearc sedimentary and volcanoclastic formations structural and petrographic investigations show that the lichi melange most probably formed as a result of internal deformation of the forearc when the continental shelf of south china collided with the luzon arc as a result of the subduction of the south china sea beneath the philippine sea plate the forearc formations constitute the protolith of the lichi melange it seems improbable that the mechanical properties of the minerals of the matrix illite chorite kaolinite in themselves explain the aseismic behavior of the lvf microstructural investigations show that deformation within the fault zone must have resulted from a combination of frictional sliding at grain boundaries cataclasis responsible for grain size comminution and pressure solution creep responsible for the development of the scaly foliation and favored by the mixing of soluble and insoluble minerals the microstructure of the gouge formed in the lichi melange favors effective pressure solution creep which inhibits strain weakening brittle mechanisms and is probably responsible for the dominantly aseismic mode of fault slip since the lichi melange is analogous to any unlithified subduction melanges this study sheds light on the mechanisms which favor aseismic creep on subduction megathrust

variations in fault slip and strain accumulation at parkfield california initial results using two color geodimeter measurements 1984 1988 repeated length measurements of several geodetic baselines near parkfield california have revealed significant variations in the local rates of shallow fault slip and strain this network of baselines surveyed several times each week straddles the san andreas fault in the transition zone between the creeping section to the northwest and the locked section to the southeast the length measurements characterized by a precision approaching 0 1 ppm reveal large fluctuations in the rates of baseline extension principal mode analysis of the length change data indicates that the two largest components of the signal are 1 secular extensions and contractions consistent with surface slip on the main strand of the san andreas fault and 2 a large seasonal oscillation with no obvious spatial coherence on most of these baselines the second component appears to be in phase with seasonal rainfall when data from the baselines with the largest amplitude of the seasonal signal are excluded the remaining data can be modeled in terms of both spatial and temporal variations in surface slip variations in the components of the changes in uniform strain and the possible displacement of the central monument in this radial network in parameterizing this model the spatial variation of slip beneath the near surface is reflected by changes in shear strain although the computed secular shear is highly dependent upon the specified parameterization of surface slip the data are consistent with the hypothesis that slip at intermediate depths lags behind the surface slip rate however the range in models that fit the data does not necessarily imply that there is a deficit in slip at depth relative to the surface comparison of the inferred values of surface slip from the model with the observed fault slip measured by very short baseline creep meters indicates close agreement in secular rates but the short term variations observed with the creep meters are either highly attenuated or nonexistent in the modeled slip since the modeled slip is a spatial average which smooths out possible short wavelength variations in the surface slip for which the creep instruments are most sensitive an interesting conclusion from the two color data is that surface slip on the san andreas fault appears to be spread over a 2 km wide zone on the south flank of middle mountain but is confined to a very narrow zone to the south as the fault passes through the center of the network this conclusion is dependent upon the assumption that a few critical monuments are stable and track tectonic displacements in the long term finally the largest observed strain change is an extensional strain coincident with the kettleman hills earthquake m5 5 in august 1985

time dependent triggered afterslip following the 1989 loma prieta earthquake we employ a modified version of the network inversion filter to investigate time dependent slip following the 1989 loma prieta earthquake previous analysis of global positioning system gps and leveling data suggests afterslip on the loma prieta rupture as well as aseismic slip on a thrust fault northeast of the san andreas fault which we identify with the foothills thrust belt we analyzed 173 daily gps solution files at 62 stations collected from 1989 8 to 1998 3 a total of 1 134 three dimensional relative baseline determinations the observed position changes are assumed to result from secular deformation random benchmark motions and temporally varying fault slip the data reveal temporal variations in slip rate but poorly resolve spatial variations in fault slip the amount of temporal smoothing is estimated by maximum likelihood conditional on this estimate reverse slip on the foothills thrust decays from 45 12 mm yr immediately after the earthquake to zero by 1992 reverse slip on the loma prieta rupture surface decays from 57 11 mm yr to zero by 1994 right lateral slip on the loma prieta rupture surface decays monotonically from 30 10 mm yr to zero by 1994 these results suggest that 1 triggered afterslip can occur off the main rupture zone on adjacent faults 2 shallow afterslip dominated the postseismic deformation for the 8 years following the earthquake and 3 postseismic slip on the foothills thrust may account for a significant portion of its total slip budget

fault slip rates at depth from recurrence intervals of repeating microearthquakes unique attributes in sequences of recurring similar microearthquakes at parkfield california provide a means for inferring slip rate at depth throughout the active fault surface from the time intervals between sequence events application of the method using an 11 year microseismicity record revealed systematic spatial and temporal changes in the slip rate that were synchronous with earthquake activity and other independent measures of fault zone slip if this phenomenon is found to be generally common behavior in active faults it forms the basis for a method to monitor the changing strain field throughout a seismogenic fault zone

anomalous shear strain at parkfield during 1993 94 three borehole tensor strainmeters btsm installed near parkfield and separated by more than 5km have established a good baseline of data spanning the past eight years a significant change in the accumulation rate of shear strain 0 5 and 1 0 microstrain per year at two of the sites commencing approximately in mid 1993 and continuing to 1995 is evident in these data some supporting evidence for a change in the strain rate at that time is seen in dilatometer data geodetic measurements fault zone properties and microearthquake locations there was also an increase in medium level seismic activity m 4 7 on 10 20 92 m 3 9 on 10 26 92 m 3 5 on 3 13 93 m 4 4 on 4 3 93 m 4 8 on 11 14 93 in the middle mountain area while a rainfall induced strain may have caused the observed strain anomaly we propose a more likely explanation is that aseismic slip at depth on a section of the fault south of middle mountain produced the observed strain anomaly this anomaly is unique in the current data set

fault modeling of the 2012 wutai taiwan earthquake and its tectonic implications abstract the mw 5 9 wutai earthquake of 26 february 2012 occurred at a depth of 26 km in southern taiwan where the rupture is not related to any known geologic structures to illustrate the rupture source of the mainshock we employ an elastic half space model and gps coseismic displacements to invert for optimal fault geometry and coseismic slip distribution with observations of both coseismic horizontal and vertical displacements less than 10 mm our preferred fault model strikes 312 and dips 30 to the northeast and exhibits a reverse slip of 28 112 mm and left lateral slip of 9 45 mm estimated geodetic moment of the wutai earthquake is 1 3 1018 n m equivalent to an mw 6 0 earthquake the wutai epicentral area is characterized by a ne sw compression as evidenced by the slaty cleavage orientations and the interpretation of stress tensor inversion of earthquake focal mechanisms using the stress drops of the wutai and the nearby 2010 mw 6 4 jiashian earthquakes we obtain a lower bound of 0 002 for the coefficient of friction on the fault on the other hand studying the crustal thickness contrast in southern taiwan provides an upper bound of the average horizontal compressive force of 1 67 1012 n m transmitted through the taiwan mountain belt and gives an estimate of the maximum friction coefficient for 0 03 the deviation of an order of magnitude difference between the upper and lower bounds for the coefficient of friction suggests that other fault systems may support substantial differential stress in the lithosphere as well

the 2010 mw 7 8 mentawai earthquake very shallow source of a rare tsunami earthquake determined from tsunami field survey and near field gps data 1 the mw 7 8 october 2010mentawai indonesia earthquake was a tsunami earthquake a rare type of earthquake that generates a tsunami much larger than expected based on the seismicmagnitude itproducedalocallydevastatingtsunami withrunupcommonlyinexcess of 6 m we examine this event using a combination of high rate gps data from instruments located on the nearby islands and a tsunami field survey the gps displacement time series are deficient in high frequency energy and show small coseismic displacements 16 m our modeling results show that the combination of the small gps displacements and large tsunami can only be explained by high fault slip at very shallow depths far from the islands and close to the oceanic trench inelastic uplift of trench sediments likely contributed to the size of the tsunami recent results for the 2011 mw 9 0 tohoko oki earthquake have also shown shallow fault slip but the results from our study which involves a smaller earthquake provide much stronger constraints on how shallow the rupture can be with the majority of slip for the mentawai earthquake occurring at depths of 6 km this result challenges the conventional wisdom that the shallow tips of subduction megathrusts are aseismic and therefore raises important questions both about the mechanical properties of the shallow fault zone and the potential seismic and tsunami hazard of this shallow region

comparison of various inversion techniques as applied to the determination of a geophysical deformation model for the 1983 borah peak earthquake abstract a number of techniques are employed to overcome nonuniqueness and instability inherent in linear inverse problems to test the factors that enter into the selection of an inversion technique for fault slip distribution we used a penalty function with smoothness pf s a damped least squares method dls damped least squares method with a positivity constraint dls p and a penalty function with smoothness and a positivity constraint pf s p for inverting the elevation changes for slip associated with the 1983 borah peak earthquake unlike solving an ill posed inverse problem using a gradient technique ward and barrientos 1986 we have restored the well posed character between the elevation changes and normal slip distribution studies showed that the constraints based on sound understanding of the physical nature of the problem are crucial in the derivation of a meaningful solution and dictates primarily the selection of a particular inversion technique all available geological and geophysical information were used to determine a geophysical deformation model for the earthquake it is suggested that the pf s p solution for a fault length of 75 km is the preferred model the long wavelength features in the estimated slip distribution are similar to those obtained by ward and barrientos 1986 whereas the shorter wavelength features differ between two solutions the fault dips 49 to the southwest the slipped zones deepen from the surface at the northwest to about 20 km downdip depth at the southeast the fault extends to the southeast beyond the epicenter of the mainshock it is also shown that only the long wavelength features of the slip distribution are well resolved the resolution is better at shallower levels than at deeper levels the resolution deteriorates when the deformation sources are away from the leveling lines smoothness constraints provide better resolution than damping does at depth the addition of a positivity constraint significantly improves the model resolution

quasi static slip rate shielding by locked and creeping zones as an explanation for small repeating earthquakes at parkfield in a recent study of microearthquakes along the parkfield segment of the san andreas fault nadeau et al 1995 have found that much of the seismicity in the region is characterized by quasi periodic repeating sequences of small earthquakes that are essentially identical in waveform size and location nadeau and johnson 1998 interpreted these as repeated slip on a given asperity driven by a steady slip rate of 2 3 cm yr and concluded that the stress drops needed to be extremely high of the order of 20 kilobars we propose another explanation for these small repeating events namely that an inner asperity is surrounded by a larger creeping zone which in turn is surrounded by a still larger locked zone this geometry produces a local slip velocity much less than the overall creep velocity observed on a still larger scale slip velocity shielding we have constructed a foam rubber model to illustrate the phenomenon the time sequences of small events at the asperity punctuated by large events which rupture the whole block look very similar to the cumulative moment plots of nadeau and johnson the actual dynamic stress drops are of the same order as for the large events thus the results of the model correspond to the observations of nadeau and johnson and suggest that the model may be appropriate to explain their observations without requiring super strong asperities

the parkfield california earthquake prediction experiment five moderate magnitude 6 earthquakes with similar features have occurred on the parkfield section of the san andreas fault in central california since 1857 the next moderate parkfield earthquake is expected to occur before 1993 the parkfield prediction experiment is designed to monitor the details of the final stages of the earthquake preparation process observations and reports of seismicity and aseismic slip associated with the last moderate parkfield earthquake in 1966 constitute much of the basis of the design of the experiment

gps measurement of postseismic deformation following the 1999 chi chi taiwan earthquake 1 in the 15 month period september 1999 to december 2000 after the 1999 chi chi taiwan earthquake mw 7 6 about 80 campaign surveyed gps stations in central taiwan were occupied up to seven times furthermore seven continuous gps stations were put into operation in the epicentral region mostly within three weeks after the main shock gps data from these near field stations and some pre existing far field stations were utilized to study the postseismic deformation following the chi chi earthquake the postseismic displacements of the gps stations were approximated by a combination of an exponential transient decay function with a relaxation time of 86 days and a postseismic linear rate change stations on the hanging wall displayed west northwesterly horizontal displacements of up to 252 mm and uplift as large as 229 mm postseismic displacements are much larger to the south near the epicenter of chi chi main shock than in the north where the maximum coseismic displacement occurred stations on the footwall show only small displacements the postseismic gps data were inverted to infer the deeper fault geometry and afterslip distribution based on a four segment fault model with its shallow part constrained by coseismic fault geometry deduced from gps and seismological studies the optimal model requires the lower fault segment to be a horizontal plane at a depth of 10 4 km consistent with a model based on gps data taken during the first three months after the main shock and other geological or geophysical studies the maximum afterslip 459 mm occurs in the hypocentral region and in the northern part of the decollement slip of about 430 mm also occurs in the southern part of the decollement the afterslip moment inferred from 15 months of gps data is 4 7 1019 n m approximately 2 3 times the seismic moment released by aftershocks this implies that a major part of the postseismic deformation is aseismic

patterns and mechanisms of coseismic and postseismic slips of the 2011 mw 7 1 van turkey earthquake revealed by multi platform synthetic aperture radar interferometry on 23rd october 2011 a mw 7 1 reverse slip earthquake occurred in the bardakci saray thrust fault zone in the van region eastern turkey earlier geodetic studies have found different slip distributions in terms of both magnitude and pattern in this paper we present several cosmo skymed csk envisat asar and radarsat 2 interferograms spanning different time intervals showing that significant postseismic signals can be observed in the first three days after the mainshock using observations that combine coseismic and postseismic signals is shown to significantly underestimate coseismic slip we hence employed the csk pair with the minimum postseismic signals to generate one conventional interferogram and one along track interferogram for further coseismic modelling our best fit coseismic slip model suggests that 1 this event is associated with a buried nnw dipping fault with a preferable dip angle of 49 and a maximum slip of 6 5 m at a depth of 12 km and 2 two unequal asperities can be observed consistent with previous seismic solutions significant oblique aseismic slip with predominant left lateral slip components above the coseismic rupture zone within the first 3 days after the mainshock is also revealed by a postseismic csk interferogram indicating that the greatest principal stress axis might have rotated due to a significant stress drop during the coseismic rupture

a bayesian method for incorporating self similarity into earthquake slip inversions distributions of coseismic slip help illuminate many properties of earthquakes including fault geometry stress changes frictional properties and potential future hazard slip inversions take observations and calculate slip at depth but there are a number of commonly adopted assumptions such as minimizing the second spatial derivative of slip the laplacian that have little physical basis and potentially bias the result in light of growing evidence that fault slip shows fractal properties we suggest that this information should be incorporated into slip inversions as a form of regularization instead of laplacian smoothing we have developed a bayesian approach to efficiently solve for slip incorporating von karman regularization in synthetic tests our approach retrieves fractal slip better than laplacian regularization as expected but even performs comparably or better when the input slip is not fractal we apply this to the 2014 mw 6 0 napa valley earthquake on a two segment fault using interferometric synthetic aperture radar insar and global positioning system gps data we find the von karman and laplacian inversions give similar slip magnitude but in different locations and the von karman solution has much tighter confidence bounds on slip than the laplacian solution differences in earthquake slip due to the regularization technique could have important implications for the interpretation and modeling of stress changes on the causative and neighboring faults we therefore recommend that choice of regularization method should be routinely made explicit and justified and that von karman regularization is a better default than laplacian