scale dependence of seismic energy to moment ratio for strike slip earthquakes in japan we analyzed four pairs of a large m w 6 and a small m w 3 5 to 4 shallow strike slip earthquakes to investigate the scale dependence of the seismic energy to moment ratio an important macroscopic parameter which reflects the basic physical process of seismic slip these earthquakes occurred in the south western part of japan and high quality close in records epicentral distance 50 km are available for both the small and large earthquakes the paired events have almost the same focal mechanism and hypocenter location we used the spectral ratio of the paired events in order to remove the effects of attenuation along the wave propagation path and the station site response we then estimated the seismic energy from the source spectra estimated from the spectral ratio the energy to moment ratio increases with the earthquake size this scale dependence is very similar to that found earlier for earthquakes in southern california

complex source processes and the interaction of moderate earthquakes during the earthquake swarm in the hida mountains japan 1998 abstract this paper presents a set of detailed source models of 18 moderate m4 5 earthquakes that were part of the swarm activity beneath the hida mountains in central japan in 1998 to reveal the detailed rupture process this study uses a waveform inversion procedure with the empirical green function egf most of our models show a complicated rupture process that includes several separate rupture areas and that involves ruptures that in some events seem to repeat at almost the same location in our models we measured as characteristic values the maximums of slip slip rate stress drop and average rupture propagation velocity in comparison with the models for larger events determined by similar methods the maximum slip scales with the cubic squares of the seismic moment for large stress drop events only the slip rate has an upper limit of about 1 2 m s and the stress drop is distributed widely from 1 to 50 mpa the average stress drops estimated using the corner frequency are almost uniform at about 1 mpa the first event in each swarm cluster tends to have a large maximum stress drop which reflects the change in the stress field during swarm activity the average rupture velocity determined by the least squares method is relatively slow and shows a negative relationship with the maximum stress drop this finding cannot be explained by simple crack models and instead suggests a very complicated mechanism like the interactive fault system the interaction between the slip areas of neighbor events resulted in a non overlapping slip distribution this finding indicates that previous rupture areas restrict the size and location of later events even in such complex swarm activities

crustal structure of the diablo and gabilan ranges central california a reinterpretation of existing data the compressional wave velocity structure of the crust of the coast ranges of central california was modeled from seismic refraction data reported previously for both the diablo range east of the san andreas fault and the gabilan range west of the fault two alternative velocity models were derived for each range by iterative two dimensional ray tracing each pair of models shows a sedimentary layer and an underlying crust composed of three or four layers for the diablo range the velocities in the sediments range from 2 9 to 4 8 km sec and the depth to basement varies from 0 2 km within the central diablo range to a maximum of 3 4 km north of the livermore valley in the gabilan range the velocities in the sediments overlying the granitic basement are higher north of the gabilan range 3 6 to 4 6 km sec than to the south 2 2 to 3 8 km sec and the depth to basement both north and south of the range is less than 2 km below the sediments and fractured near surface material a resolvable difference in the crustal velocity structure on opposite sides of the san andreas fault indicates compositional differences at depth the upper crust has an average velocity of 5 7 km sec at depths between 3 and 16 3 km in the diablo range and an average velocity of 6 1 km sec at depths between 2 and 9 1 km in the gabilan range the lower crust has an average velocity of 6 9 km sec at depths between 16 3 and 30 km in the diablo range and an average velocity of 6 5 km sec at depths between 10 and 24 km in the gabilan range the models also show that the depth to moho differs by several kilometers between the two ranges in the diablo range models the crust mantle boundary becomes shallower from south to north rising from a depth of 30 to 26 km in the gabilan range models the moho is at a depth of 24 to 26 km depending on the velocity assumed at the base of the crust laboratory measurements of rock velocities and the mapped surface geology allow us to interpret the velocity models in terms of crustal composition we conclude that the diablo range probably consists largely of metagraywacke to a depth of 16 3 km and gabbroic material below this depth the crust of the gabilan range probably consists of granitic material to a depth of 9 1 km and gneissic material below this depth franciscan rocks are not regionally present in the crust of the gabilan range if the gneissic lower crust consists of the same rocks as are found at the surface in the gabilan range then the presumed large scale lateral motion of the salinian block has taken place at or below the crust mantle boundary rather than along a mid crustal slip plane

source rupture plane determination from directivity doppler effect for small earthquakes recorded by local networks abstract fault plane solutions from first motion polarity of p waves give two nodal planes among which it is not possible to distinguish the fault plane from the auxiliary one we found that for earthquakes recorded by a dense local network it is possible to identify the rupture direction from the doppler effect produced by directivity this principle was successfully applied using three different simple techniques direct inspection and time and frequency domains analyses to one m l 3 6 and six smaller 1 2 m l 2 0 well recorded earthquakes from northern baja california mexico the results indicate activity in faults conjugate to those currently mapped in this region

reflection imaging with earthquake sources and dense arrays the mw 5 8 earthquake that occurred in louisa county virginia on 23 august 2011 provided an opportunity to record with several high density seismic arrays in addition to traditional sparse temporary seismic networks traditional aftershock networks consist of a few dozen stations spread over tens of kilometers as a result the recorded seismic waveforms suffer from spatial aliasing that is so severe that many types of waveform processing are not applicable here we report the results of recording with a large number of oil industry type instruments deployed at a spacing that is an order of magnitude closer than in traditional deployments the objective was to image subsurface structure with array methods using the aftershocks as sources the dense array recorded continuously for 12 days and consisted of 172 vertical component seismometers that were placed at 200 400 m and a 60 km long three component regional profile with stations every 2 km we demonstrate how processing techniques from vertical seismic profiling can produce high resolution 3 d reflection images of structure beneath the array these images display reflectivity that correlates with that observed on a nearby deep reflection survey collected by the u s geological survey of particular interest is a strong reflector imaged across multiple profiles our analysis demonstrates how a surface array of seismometers can provide 3 d images of structure using microearthquake sources when wavefields are sampled sufficiently densely

high resolution imaging of the bear valley section of the san andreas fault at seismogenic depths with fault zone head waves and relocated seismicity summary detailed imaging of fault zone fz material properties at seismogenic depths is a difficult seismological problem owing to the short length scales of the structural features seismic energy trapped within a low velocity damage zone has been utilized to image the fault core at shallow depths but these phases appear to lack sensitivity to structure in the depth range where earthquakes nucleate major faults that juxtapose rocks of significantly different elastic properties generate a related phase termed a fault zone head wave fzhw that spends the majority of its path refracting along the fault we utilize data from a dense temporary array of seismometers in the bear valley region of the san andreas fault to demonstrate that head waves have sensitivity to fz structure throughout the seismogenic zone measured differential arrival times between the head waves and direct p arrivals and waveform modelling of these phases provide high resolution information on the velocity contrast across the fault the obtained values document along strike fault normal and downdip variations in the strength of the velocity contrast ranging from 20 to 50 per cent depending on the regions being averaged by the ray paths the complexity of the fz waveforms increases dramatically in a region of the fault that has two active strands producing two separate bands of seismicity synthetic waveform calculations indicate that geological observations of the thickness and rock type of the layer between the two strands are valid also for the subsurface structure of the fault the results show that joint analysis of fzhws and direct p arrivals can resolve important smallscale elements of the fz structure at seismogenic depths detailed characterization of material contrasts across faults and their relation to earthquake ruptures is necessary for evaluating theoretical predictions of the effects that these structures have on rupture propagation

empirical green s function study of the january 17 1994 northridge california earthquake the january 17 1994 northridge main shock was analyzed using an empirical green s function approach to estimate the spatial and temporal distribution of fault slip far field source time functions obtained from seismograms recorded by the berkeley digital seismic network bdsn and terrascope stations reveal the distributed nature of the source and the presence of a second large subevent the location of the fault beneath the san fernando valley as well as the updip directivity of the rupture contributed to the locally heavy damage due to this earthquake

aftershocks of microearthquakes as probes of the mechanics of rupture 1 using a waveform cross correlation technique rubin and gillard 2000 obtained precise relative locations for 4300 0 5 m 3 5 earthquakes occurring along 50 km of the san andreas fault this study adds to that another 5000 earthquakes distributed along 10 km of the san andreas fault and 20 km of the calaveras fault errors in relative location are typically tens of meters for earthquakes separated by hundreds of meters and after correcting for time dependent station delays meters for earthquakes separated by tens of meters along both faults the minimum separation between consecutive earthquakes scales with magnitude in a manner consistent with a magnitude independent stress drop by treating each earthquake on the san andreas as if it were a main shock scaling the distances to all subsequent earthquakes by main shock size and stacking the results a composite aftershock sequence is produced that has many of the characteristics predicted by rate and state friction models projected onto the fault surface these aftershocks outline a quasi elliptical roughly 4 mpa stress drop main shock elongate in the slip parallel direction by 40 at the ends of the major axes of this ellipse over twice as many aftershocks occur to the nw than to the se an asymmetry attributed to the contrast in material properties across the fault unlike the san andreas the calaveras fault exhibits little p wave velocity contrast and no discernible aftershock asymmetry however the earliest part of the aftershock sequence on the calaveras might be truncated by the 30 s blind time of the network following a triggering event

resolving fault plane ambiguity for small earthquakes summary we have developed an automated procedure to resolve fault plane ambiguity for small to medium sized earthquakes 2 5 ml 5 using synthetic green s tensors computed in a 3 d earth structure model and applied this procedure to 35 earthquakes in the los angeles area for 69 per cent of the events we resolved fault plane ambiguity of our cmt solutions at 70 per cent or higher probability for some earthquakes the fault planes selected by our automated procedure were confirmed by the distributions of relocated aftershock hypocentres in regions where there are no precisely relocated aftershocks or for earthquakes with few aftershocks we expect our method to provide the most convenient means for resolving fault plane ambiguity our procedure does not rely on detecting directivity effects therefore it is applicable to any types of earthquakes

the rupture extent of low frequency earthquakes near parkfield ca the low frequency earthquakes lfes that constitute tectonic tremor are often inferred to be slow to have durations of 0 2 to 0 5 s a factor of 10 to 100 longer than those of typical m w 1 2 earthquakes here we examine lfes near parkfield ca in order to assess several proposed explanations for lfes long durations we determine lfe rupture areas and location distributions using a new approach similar to directivity analysis where we examine how signals coming from various locations within lfes finite rupture extents create differences in the apparent source time functions recorded at various stations we use synthetic ruptures to determine how much the lfe signals recorded at each station would be modified by spatial variations of the source station travel time within the rupture area given various possible rupture diameters and then compare those synthetics with the data our synthetics show that the methodology can identify inter station variations created by heterogeneous slip distributions or complex rupture edges and thus lets us estimate lfe rupture extents for unilateral or bilateral ruptures to obtain robust estimates of the sources similarity across stations we stack signals from thousands of lfes using an empirical green s function approach to isolate the lfes apparent source time functions from the path effects our analysis of lfes in parkfield implies that lfes apparent source time functions are similar across stations at frequencies up to 8 to 16 hz depending on the family the inter station coherence observed at these relatively high frequencies or short wave lengths down to 0 2 to 0 5 km suggest that lfes in each of the 7 families examined occur on asperities they are clustered in patches with sub 1 km diameters the individual lfes rupture diameters are estimated to be smaller than 1 1 km for all families and smaller than 0 5 km and 1 km for the two shallowest families which were previously found to have 0 2 s durations coupling the diameters with the durations suggests that it is possible to model these m w 1 2 lfes with earthquake like rupture speeds around 70 of the shear wave speed however that rupture speed matches the data only at the edge of our uncertainty estimates for the family with highest coherence the data for that family are better matched if lfes have rupture velocities smaller than 40 of the shear wave speed or if lfes have different rupture dynamics they could have long rise times contain composite sub ruptures or have slip distributions that persist from event to event

estimating finite source properties of small earthquake ruptures many of the most fundamental questions in earthquake science are currently limited by a lack of knowledge about small earthquake ruptures small earthquakes are difficult to study owing to the poor constraints placed on many of the interesting physical parameters by bandlimited far field seismic data traditionally dynamic models such as an expanding circular crack have been utilized to bridge the gap between the easily measurable quantities for small earthquakes and more interesting physical parameters such as stress drop and rupture velocity here i present a method for estimating the basic finite source properties of a rupture that is independent of any a priori model and utilizes the description of a finite source the second moments that far field waves are inherently sensitive to application to two magnitude 5 events in southern california demonstrates the ability of an empirical green9s function approach to estimating the second moments to resolve the fault plane ambiguity rupture length and overall directivity additional results are presented for two example m 2 7 events from the creeping section of the san andreas fault to examine the likely lower bound on event size that can be studied with surface seismometers the creeping section earthquakes have very similar rupture areas but would be improperly interpreted as significantly different using the traditional methodology one of these events presents a relatively clear interpretation of the velocity of rupture front propagation which is about 0 8 of the rayleigh speed suggesting little difference in rupture velocity between it and typical large earthquakes

source parameters and rupture velocity of small m 2 1 reservoir induced earthquakes summary we calculate stress drop and rupture speed for ml 2 1 shallow reservoir induced earthquakes and find them to be similar to those of large natural earthquakes previous studies have suggested that hydrofractures mining and reservoir induced earthquakes have lower average stress drop than natural tectonic earthquakes this difference might result from the different tectonic setting or the shallower hypocentral depths of induced earthquakes alternatively difficulties in correcting for attenuation and site effects in earlier studies may lead to underestimation of stress drop in addition most studies assume the rupture velocity of small reservoir induced earthquakes to be the same as for the large earthquakes we analyse a set of 101 ml 2 1 earthquakes induced by changing water level in the acu reservoir ne brazil the earthquakes are shallow depth 5 km and the region has negligible natural seismicity we use three different approaches to calculate the source parameters of the six largest 1 9 ml 2 1 earthquakes we model the individual spectra to find corner frequency frequency independent q and long period amplitude we use collocated small earthquakes as empirical green s functions to calculate the spectral ratios and determine the relative source time functions estimates of the source duration and corner frequency imply stress drops in the range of 26 179 mpa these are similar to or higher than tectonic earthquakes suggesting that the shallow hypocentral depth and the presence of water do not affect stress drop we observe clear directivity for one of the earthquakes and use the azimuthal variation in pulse width to estimate a rupture velocity of 0 6

fine structure of the rupture zone of the april 26 and 27 1997 northridge aftershocks we investigated the rupture geometry of two m w 4 5 earthquakes april 26 1997 and april 27 1997 that occurred on the western edge of the aftershock zone of the 1994 m w 6 7 northridge earthquake both events have thrust mechanisms with a steep plane dipping 75 se and a shallow plane dipping 45 ne an empirical green s function deconvolution followed by a waveform inversion was used to determine the slip distribution of the two events the inversion results show that the steep plane fits the data slightly better than the shallow plane the background seismicity aftershocks of the 1994 northridge event in the epicentral region shows the existence of a north dipping fault plane similar to that of the 1971 san fernando earthquake however the spatial trend of the two april 1997 events and their aftershocks reveals tightly clustered seismicity on a steep plane dipping south relative relocation of the april 26 and 27 events shows that the april 27 event ruptured 1 4 km n70 e of the april 26 event and at a slightly shallower depth i e almost along strike and on a steep plane dipping south these observations suggest that the steep plane is the fault plane of both events the two events ruptured on a plane which is almost perpendicular to the trend of the regional background seismicity thus the seismogenic structure beneath the transverse ranges exhibits complexity on scales of a few kilometers the rupture area for both the april 26 and 27 events is 1 km 2 with a stress drop of at least 20 to 30 bars

predominance of unilateral rupture for a global catalog of large earthquakes the manner in which an earthquake rupture propagates across a fault reflects both the initial properties of the fault and the dynamical stresses produced by the rupture we quantify the propagation of an earthquake rupture using the second moments of the earthquake space time distribution in particular the second moments provide a simple way to differentiate between approximately bilateral and predominantly unilateral ruptures we determined the second moments for a catalog of m w 7 earthquakes that have occurred since 1994 the results show that approximately 80 of large shallow ruptures are predominantly unilateral our result agrees well with strong motion inversions from previous studies on moderate and large earthquakes the predominance of unilateral propagation in large earthquakes is explained by a simple characteristic earthquake model in which large events rupture an entire structurally defined fault segment and nucleation points are uniformly distributed along the fault the predominance of unilateral rupture may also be enhanced by the dynamic stress field produced by contrasts in the elastic properties between the two sides of plate boundary faults manuscript received 16 november 2001

an inversion method to analyze rupture processes of small earthquakes using stopping phases 1 we propose a method to obtain source parameters most importantly the rupture velocity of small earthquakes on the basis of stopping phases it is assumed that a rupture begins at one focus of an ellipse and spreads radially at a constant rupture velocity on the fault surface eventually terminating on the elliptical boundary this model allows us to investigate a wide variety of rupture configurations ranging from a circular to a unilaterally propagating rupture in this model we find that two high frequency stopping phases hilbert transformations of each other are radiated and that the difference in arrival times between the two stopping phases is dependent on the average value of the rupture velocity the radius of the major semiaxes of the ellipse and its ellipticity these parameters can be estimated by a nonlinear least squares inversion method we perform simulation tests which allow us to detect a hilbert transform pair whose arrival times coincide with ones predicted by an assumed source model from noisy data these simulation tests show that our inversion method is useful since our method does not require an a priori assumption of the average rupture velocity the orientation of the fault plane and the circular crack model it provides us with a tool for estimating accurate source dimensions as well as rupture processes of small earthquakes

disputable non double couple mechanisms of several strong earthquakes second degree moment approach the earthquake source mechanism description of the geometry of the source and its strength in terms of forces equivalent to rupturing of the rock mass is routinely modeled as a moment tensor mt describing the earthquake focus as a point however sometimes this approximation is not satisfied and seismic radiation keeps the directivity due to the finite extent of the source inversion into the mt then may yield a biased mechanism synthetic study in the adamova and silen 2010 demonstrated the appearance of spurious non double couple components in the mechanism even for a pure double couple dc source their method was designed to reduce the spurious source components by using the second degree moments to evaluate their contribution into the records and subtract it from the data here we applied the procedure to five moderate to large regional events with large nonshear components they are mostly located on large tectonic faults where predominantly pure shear slip is expected we studied one event on the prominent north anatolian fault three events in the pacific area and one event in bolivia in most cases the non dc components essentially were reduced and the geometry of the mechanism remained largely unchanged this confirms the hypothesis that part of the non dc components in regional mt solutions may be spurious due to the neglect of the source finiteness in the routine procedure of the mt retrieval in addition the geometrical and kinematical characteristics of the foci provided by the second degree moments the source ellipsoid and rupture velocity vector are mostly consistent with the fault geometry aftershock distribution and estimates of rupture speed from available previous studies

finite moment tensor of the 3 september 2002 yorba linda earthquake we have developed procedures for inverting broadband waveforms for the finite moment tensors fmts of regional earthquakes the fmt is defined in terms of second order polynomial moments of the source space time function it removes the fault plane ambiguity of the centroid moment tensor cmt and yields several additional parameters of seismological interest the characteristic length lc width wc and duration tc of the faulting as well as the directivity vector v d of the fault slip following mcguire et al 2001 we represent the observed waveform relative to the synthetic in terms of two frequency dependent differential times a phase delay p and an amplitude reduction time q which we measure using the generalized seismological data functional gsdf method gee and jordan 1992 we numerically calculate the fmt partial derivatives which allows us to use synthetics computed by using any forward modeling tools we have tested our methodology on southern california seismic network scsn recordings of the 03 september 2002 yorba linda earthquake m w 4 3 using 1d synthetic green s functions we determined the cmt and resolved fault plane ambiguity to resolve the details of source finiteness we employed a joint inversion technique that recovers the cmt parameters of the aftershocks and the cmt and fmt parameters of the mainshock the joint system of equations relating the data to the source parameters of the mainshock aftershock cluster is denuisanced for path anomalies in both observables this projection operation effectively corrects the mainshock data for path related anomalies in a way similar to but more flexible than empirical green s function egf techniques our results indicate that the yorba linda rupture occurred as left lateral slip on a fault patch conjugate to the nearby right lateral whittier fault we obtained source dimensions of lc 0 7 0 1 km wc 0 4 0 1 km and tc 0 2 0 05 sec implying a stress drop of about 3 2 mpa and we found a directivity of 0 8 0 2 oriented up and to the northeast the inferred fault plane is consistent with the mainshock aftershock distribution relocated by hauksson et al 2002

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

crustal stress field in southern california and its implications for fault mechanics we present a new high spatial resolution image of stress orientation in southern california based on the inversion of earthquake focal mechanisms we use this image to study the mechanics of faulting in the plate boundary region the stress field contains significant spatial heterogeneity which in some cases appears to be a result of the complexity of faulting and in other cases appears to be a cause temporal changes in the stress field are also observed primarily related to major earthquakes the observed 15 10 rotation of the stress axes due to the 1992 m7 3 landers mainshock implies that the deviatoric stress magnitude in the crust is low of the order of 10 mpa this suggests that active faults in southern california are weak the maximum principal stress axis near the san andreas fault is often at 50 to the fault strike indicating that the shear stress on the fault is comparable to the deviatoric stress the san andreas in southern california may therefore be a weak fault in a low strength crust

aftershock asymmetry rupture directivity among central san andreas fault microearthquakes using a waveform cross correlation technique we have obtained precise relative locations for nearly 75 of the northern california seismic network catalog 4300 earthquakes occurring between 1984 and 1997 along 50 km of the san andreas fault errors in relative location are meters to tens of meters for events separated by tens to hundreds of meters we find that consecutive earthquakes in the relocated catalog occur no closer than a distance approximately equal to the radius of the first rupture as estimated from the moment magnitude relationship of abercrombie 1996 assuming a 10 mpa stress drop when the relative position vectors between consecutive events are normalized by this distance and projected onto the fault surface they define a hole whose shape suggests that typical microearthquakes are elongate in the mode ii slip parallel direction by several tens of percent moreover of the 100 immediate aftershocks occurring closest to the mode ii edges of the prior rupture more than twice as many occur to the northwest than to the southeast we interpret this asymmetry as resulting from the large contrast in material properties across the fault models of dynamic rupture between dissimilar media predict that ruptures in this region may run preferentially to the southeast in the direction of motion of the lower velocity material if so then the barriers that stop rupture fronts moving to the southeast should initially be farther from failure on average than the barriers that stop rupture fronts moving to the northwest once the rupture stops the induced stress change is more symmetric but the fault remains farther from failure on average to the southeast this interpretation receives some support from pulse width measurements on a localized set of 72 magnitude 0 6 to 3 6 earthquakes

source process and tectonic implication of the january 20 2007 odaesan earthquake south korea abstract the source process for the 20th of january 2007 mw 4 5 odaesan earthquake in south korea is investigated in the low and high frequency bands using velocity and acceleration waveform data recorded by the korea meteorological administration seismographic network at distances less than 70 km from the epicenter synthetic green functions are adopted for the low frequency band of 0 1 0 3 hz by using the wave number integration technique and the one dimensional velocity model beneath the epicentral area an iterative technique was performed by a grid search across the strike dip rake and focal depth of rupture nucleation parameters to find the best fit double couple mechanism to resolve the nodal plane ambiguity the spatiotemporal slip distribution on the fault surface was recovered using a non negative least square algorithm for each set of the grid searched parameters the focal depth of 10 km was determined through the grid search for depths in the range of 6 14 km the best fit double couple mechanism obtained from the finite source model indicates a vertical strike slip faulting mechanism the nw faulting plane gives comparatively smaller root mean squares rms error than its auxiliary plane slip pattern event provides simple source process due to the effect of low frequency that acted as a point source model three empirical green functions are adopted to investigate the source process in the high frequency band a set of slip models was recovered on both nodal planes of the focal mechanism with various rupture velocities in the range of 2 0 4 0 km s although there is a small difference between the rms errors produced by the two orthogonal nodal planes the sw dipping plane gives a smaller rms error than its auxiliary plane the slip distribution is relatively assessable by the oblique pattern recovered around the hypocenter in the high frequency analysis indicating a complex rupture scenario for such moderate sized earthquake similar to those reported for large earthquakes

can observations of earthquake scaling constrain slip weakening summary we use observations of earthquake source parameters over a wide magnitude range m w 0 7 to place constraints on constitutive fault weakening the data suggest a scale dependence of apparent stress and stress drop both may increase slightly with earthquake size we show that this scale dependence need not imply any difference in fault zone properties for different sized earthquakes we select 30 earthquakes well recorded at 2 5 km depth at cajon pass california we use individual and empirical green s function spectral analysis to improve the resolution of source parameters including static stress drop and total slip s we also measure radiated energy e s w ecompare the cajon pass results with those from larger california earthquakes including aftershocks of the 1994 northridge earthquake and confirm the results of abercrombie 1995 e s m 0 where rigidity and both e s m 0 and increase as m 0 and s increases uncertainties remain large due to model assumptions and variations between possible models and earthquake scale independence is possible within the resolution assuming that the average trends are real we define a quantity g 2 es m 0 s 2w hich is the total energy dissipation in friction and fracture minus 1 s w here 1 is the final static stress if 1 d the dynamic shear strength during the last increments of seismic slip then g g the fracture energy in a slip weakening interpretation of dissipation

teleseismic inversion for the second degree moments of earthquake space time distributions summary determining detailed descriptions of the space time distribution of an earthquake s rupture using only teleseismic data is often a difficult and non unique process however a description of an earthquake that goes beyond a point source is often needed for many seismotectonic problems we have developed a method that inverts measurements of the shifts in traveltime and amplitude of various seismic phases at global stations for the second degree polynomial moments of an event s space time distribution our values for the second moments of the 1995 mw 8 0 jalisco mexico thrust event correspond to a characteristic rupture length of 121t10 km along strike a characteristic rupture width of 76t11 km downdip a characteristic duration of 29t2 s and a direction of rupture propagation along the strike of the subduction zone to the nw for the 1995 mw 7 2 gulf of aqaba earthquake our estimates of the second moments correspond to a characteristic rupture length of 53t2 km a characteristic duration of 12t1 s and rupture propagation to the ne along the strike of the dead sea transform our estimates of the first and second moments agree well with values from local seismic and geodetic networks for both the jalisco and aqaba events we are able to resolve the fault plane ambiguity associated with the events moment tensors using our estimates of the second spatial moment

two dimensional seismic image of the san andreas fault in the northern gabilan range central california evidence for fluids in the fault zone a joint inversion for two dimensional p wave velocity vp p to s velocity ratio vp vs and earthquake locations along the san andreas fault saf in central california reveals a complex relationship among seismicity fault zone structure and the surface fault trace a zone of low vp and high vp vs lies beneath the saf surface trace safst extending to a depth of about 6 km most of the seismic activity along the saf occurs at depths of 3 to 7 km in a southwest dipping zone that roughly intersects the safst and lies near the southwest edge of the low vp and high vp vs zones tests indicate that models in which this seismic zone is significantly closer to vertical can be confidently rejected a second high vp vs zone extends to the northeast apparently dipping beneath the diablo range another zone of seismicity underlies the northeast portion of this vp vs high the high vp vs zones cut across areas of very different vp values indicating that the high vp vs values are due to the presence of fluids not just lithology the close association between the zones of high vp vs and seismicity suggests a direct involvement of fluids in the faulting process

determination of rupture duration and stress drop for earthquakes in southern california a simple technique is developed for determining the rupture duration and stress drop of earthquakes between magnitudes 3 5 and 4 0 using the time between the p wave onset and the first zero crossing 1 2 on seismograms from local seismic networks this method is applied to 10 main shocks in southern california to investigate regional variations in stress drop the initial pulse widths of 65 foreshocks or aftershocks of these events were measured values of 1 2 for small earthquakes below about magnitude 2 2 are generally observed to remain constant with decreasing magnitude in four sequences studied the relative pulse width of a particular main shock m 3 5 at a given station is found to be correlated with the relative pulse width of its aftershocks recorded at that station these observations are interpreted to signify that the waveforms of these small events m 2 2 are essentially the impulse response of the path between the source and receiver values of 1 2 determined from small foreshocks and aftershocks are therefore subtracted in effect deconvolved from those of each main shock to obtain an estimate of the rupture duration of the main shock which is corrected for path effects significant variations in rupture duration and stress drop are observed for the main shocks studied aftershock locations and azimuthal variations in 1 2 both indicate that the rupture zone of one earthquake expanded unilaterally a factor of 10 variation in stress drop is calculated for two adjacent events of similar seismic moments occurring 1 hr apart on the san jacinto fault system the first event in this pair had the highest stress drop of the events studied 860 bars and was followed within 8 months by a magnitude 5 5 earthquake 2 km away

total moment spectra of fourteen large earthquakes total moment spectra mt m where m is the moment rate tensor are computed for 14 large earthquakes recorded by the international deployment of accelerometers ida network using the scalar moment retrieval method proposed by silver and jordan 1982 for each event we obtain estimates of mt averaged over the 10 disjunct 1 mhz intervals in the low frequency band 1 11 mhz typical ida record sets from events with mt 0 2 a 1 a 1027 dyne cm yield standard errors on the 1 mhz averages that are generally less than 20 our multiple band estimates of mt are usually consistent with comparable single band values found by other investigators from the total moment spectra we derive the zero frequency static moment mt 0 and the characteristic source duration c 2 var f t where f t is the time function of the moment rate tensor assumed to be the same for all components the parameter c formally depends on the second degree temporal spatial and mixed spatial temporal moments of an extended source but calculations with realistic source geometries indicate that the interpretation of tc strictly in terms of the second temporal moment leads to very little error a plot of c against shows considerable scatter some events lie significantly above the empirical scaling curve of kanamori and given 1981 slow earthquakes and some below fast earthquakes examples of the former include all three deep focus events analyzed here as well as the 1970 colombia and 1963 peru bolivia deep focus earthquakes whose total moment spectra are calculated from the moment tensor solutions of gilbert and dziewonski 1975 an example of the latter is the great sumbawa earthquake of august 9 1977 24 3 a which is inferred to have a relatively shallow spatial centroid 20 km and a nearly flat moment spectrum suggesting that the faulting was confined to the upper more brittle portion of the oceanic lithosphere three other large earthquakes tonga june 22 1977 23 2 a kuril islands december 6 1978 3 6 0 2 a and tumaco december 12 1979 25 2 a exhibit moment spectra which significantly decrease toward higher frequencies indicating larger values of c the aftershocks of both the tonga and kuril islands events are distributed over a considerable range of depths consistent with rupture into the lower more ductile portions of the oceanic lithosphere we speculate that the characteristic source duration may be causally related to the depth of coseismic rupturing within the lithosphere