slip parallel seismic lineations on the northern hayward fault california a high resolution relative earthquake location procedure is used to image the fine scale seismicity structure of the northern hayward fault california the seismicity defines a narrow near vertical fault zone containing horizontal alignments of hypocenters extending along the fault zone the lineations persist over the 15 year observation interval implying the localization of conditions on the fault where brittle failure conditions are met the horizontal orientation of the lineations parallels the slip direction of the fault suggesting that they are the result of the smearing of frictionally weak material along the fault plane over thousands of years

seismotectonics of the loma prieta california region determined from three dimensional v p v p v s and seismicity three dimensional vp and vp vs velocity models for the loma prieta region were developed from the inversion of local travel time data 21 925 p arrivals and 1 116 s arrivals from earthquakes refraction shots and blasts recorded on 1700 stations from the northern california seismic network and numerous portable seismograph deployments the velocity and density models and microearthquake hypocenters reveal a complex structure that includes a san andreas fault extending to the base of the seismogenic layer a body with high vp extends the length of the rupture and fills the 5 km wide volume between the loma prieta mainshock rupture and the san andreas and sargent faults we suggest that this body controls both the pattern of background seismicity on the san andreas and sargent faults and the extent of rupture during the mainshock thus explaining how the background seismicity outlined the along strike and depth extent of the mainshock rupture on a different fault plane 5 km away new aftershock focal mechanisms based on three dimensional ray tracing through the velocity model support a heterogeneous postseismic stress field and can not resolve a uniform fault normal compression the subvertical or steeply dipping san andreas fault and the fault surfaces that ruptured in the 1989 loma prieta earthquake are both parts of the san andreas fault zone and this section of the fault zone does not have a single type of characteristic event

changing seas in the late miocene northern aegean a paratethyan approach to mediterranean basin evolution abstract the northern aegean region evolved during the miocene as a restricted land locked basin with small ephemeral connections to both the eastern paratethys former black sea and mediterranean its biostratigraphic data show mixed paratethys mediterranean components but the paratethys fauna has generally been neglected for chronologic reconstructions here we review this biostratigraphic data from a paratethyan perspective and present revised paleogeographic reconstructions of the northern aegean throughout the late miocene in the tortonian all sub basins show mainly fluvio deltaic terrestrial environments with a series of scattered lakes that are predominantly fed by local rivers and short lived paratethys connections the first persisting marine conditions still alternating with brackish paratethyan environments indicate a middle messinian late maeotian age 6 9 6 1 ma when the region formed a semi isolated egemar sea with multiple marine influxes the termination of marine conditions is very well documented by a marked palaeoenvironmental change to the brackish water environments that correlate to the maeotian pontian boundary 6 1 ma in eastern paratethys during the messinian salinity crisis 5 97 5 33 ma the northern aegean was a brackish water system lake egemar that formed a passageway for paratethyan overspill waters towards the mediterranean we conclude that the thick evaporites of the northern aegean domain do not reflect the classic mediterranean msc sequences but are more likely related to older badenian or maeotian salinity incursions

deep fault plane geometry inferred from multiplet relative relocation beneath the south flank of kilauea dense microearthquake swarms occur in the upper south flank of kilauea providing multiplets composed of hundreds of events the similarity of their waveforms and the quality of the data have been sufficient to provide accurate relative relocations of their hypocenters a simple and efficient method has been developed which allowed the relative relocation of more than 250 events with an average precision of about 50 m horizontally and 75 m vertically relocation of these events greatly improves the definition of the seismic image of the fault that generates them indeed relative relocations define a plane dipping about 6 northward although corresponding absolute locations are widely dispersed in the swarm a composite focal mechanism built from events providing a correct spatial sampling of the multiplet also gives a well constrained northward dip of about 5 to the near horizontal plane this technique thus collapses the clouds of hypocenters of single event locations to a plane coinciding with the slip plane revealed by previous focal mechanism studies we cannot conclude that all south flank earthquakes collapse to a single plane there may locally be several planes perhaps with different dips and depths throughout the south flank volume the 6 northward dipping plane we found is too steep to represent the overall flexure of the oceanic crust under the load of the island of hawaii this plane is probably an important feature that characterizes the basal slip layer below the upper south flank of kilauea volcano differences in seismicity rate and surface deformations between the upper and lower south flank could be related to the geometry of this deep fault plane the present work illustrates how high precision relative relocations of similar events in dense swarms combined with the analysis of geodetic measurements can help to describe deep fault plane geometry systematic selection and extensive relative relocation of similar earthquakes could be attempted in other well instrumented highly seismic areas to provide reliable basic information especially useful for understanding of earthquake generation processes

monitoring velocity variations in the crust using earthquake doublets an application to the calaveras fault california we present a technique that greatly improves the precision in measuring temporal variations of crustal velocities using an earthquake doublet or pair of microearthquakes that have nearly identical waveforms and the same hypocenter and magnitude but occur on different dates we compute differences in arrival times between seismograms recorded at the same station in the freqency domain by cross correlation of short windows of signal a moving window analysis of the entire seismograms including the coda gives t the difference in arrival times versus running time along the seismogram the time resolution of the method is an order of magnitude better than the digitization interval the t technique is illustrated with a pair of microearthquakes m 1 7 and 2 0 that occurred before and after the coyote lake california earthquake m 5 9 of august 6 1979 and on the same segment of the calaveras fault that ruptured during the earthquake the coda wave arrivals for some stations are progressively delayed for the second earthquake in the doublet so that its seismogram appears as a stretched version of the earlier event we interpret this systematic variation in t along the coda as a change in the average s velocity in the upper crust in the time interval between the two doublets s wave velocities appear to have decreased by 0 2 in an oblong region 5 10 km in radius at the south end of the aftershock zone

measurements of the strain field associated with episodic creep events on the san andreas fault at san juan bautista california over the past 6 years repeated strain steps recorded on a tensor strainmeter installed close to the san andreas fault trace at san juan bautista in california correlated closely with episodic creep events registered on a nearby creepmeter the strain events were remarkably similar in character were of about 1 hour duration and were followed within hours to days by creep events of 2 6 days duration the episodic slip of a source region from 200 m to 500 m in depth and at most a few kilometers in length is consistent with the observed strain and creep data we postulate that the additional stress transferred to this region by the loma prieta earthquake is giving rise to an increased rate of accumulation of shear strain around a locked patch beneath the creep source region thus accounting for the increased frequency of the creep strain events following the loma prieta earthquake

a slow earthquake sequence on the san andreas fault earthquakes typically release stored strain energy on timescales of the order of seconds limited by the velocity of sound in rock over the past 20 years observations1 13 and laboratory experiments14 have indicated that rupture can also occur more slowly with durations up to hours such events may be important in earthquake nucleation15 and in accounting for the excess of plate convergence over seismic slip in subduction zones the detection of events with larger timescales requires near field deformation measurements in december 1992 two borehole strainmeters close to the san andreas fault in california recorded a slow strain event of about a week in duration and we show here that the strain changes were produced by a slow earthquake sequence equivalent magnitude 4 8 with complexity similar to that of regular earthquakes the largest earthquakes associated with these slow events were small local magnitude 3 7 and contributed negligible strain release the importance of slow earthquakes in the seismogenic process remains an open question but these observations extend the observed timescale for slow events by two orders of magnitude

earthquake source scaling relationships from 1 to 5 ml using seismograms recorded at 2 5 km depth the scaling relationships of earthquake sources less than about magnitude 3 have been the subject of considerable controversy over the last two decades studies of such events have shown a tendency for the constant stress drop self similarity of larger earthquakes to breakdown at small magnitudes and an apparent minimum source dimension of about 100 m has been observed other studies showed that this apparent breakdown in scaling could be an artifact of severe near surface attenuation limiting the spatial resolution of surface data in this study source parameters are determined for over 100 nearby tectonic earthquakes from recordings at a depth of 2 5 km in granite in the cajon pass scientific drill hole southern california comparison of the seismograms recorded at this depth with those at the wellhead clearly demonstrates the effect of the severe attenuation in the upper kilometers of the earth s crust source parameters are calculated by spectral modeling of three component p and s waves assuming four source models based on the brune 2 n 2 model in model l n 2 is fixed and q of p and s waves is determined to be 912 581 1433 and 1078 879 1323 respectively the numbers in parentheses are 1 standard deviation in model 2 qp qs 1000 is assumed and n is allowed to vary the 2 model is a good average for the data set but there is some real scatter supported by the data in model 3 qp qs 1000 is also assumed and 2 is constrained and in model 4 attenuation is ignored and n is allowed to vary source dimensions of less than 10m are observed for all four models 10 times smaller than the proposed minimum no breakdown in constant stress drop scaling is seen in the downhole data approximately ml 1 to 5 5 m0 109 1016 nm the ratio between radiated seismic energy estimated by integrating the velocity squared spectra with adequate signal bandwidth and seismic moment appears to decrease gradually with decreasing moment in the magnitude range 1 to 7 this is not incompatible with a constant stress drop but could result from errors in calculating energy the ratio of the s wave energy to that radiated by the p waves is about 14 after correction for attenuation this low value is consistent with the corner frequency shift of about 1 3 this corner frequency shift is observed for all four source models and therefore is interpreted as being source controlled

source array analysis of coda waves near the 1989 loma prieta california mainshock implications for the mechanism of coseismic velocity changes coda waves are often considered to be generated by backscattering of primary waves from randomly distributed heterogeneities in the crust however modeling and experimental work indicates that much of the coda may be due to scattering near the receiver knowing where the coda is generated is important for several reasons coda waves have been used to characterize site amplification scattering and attenuation throughout the crust and velocity changes associated with large earthquakes however a lack of knowledge of where coda waves are generated can lead to ambiguous interpretations in this study we analyze 26 s of coda waves recorded at up to 78 stations using slowness stacking on two source arrays and find that at nearly all stations regardless of distance from the source nearly all the arrivals are within one hemisphere most are strongly clustered near the direct arrival and the clusters are in general symmetrically disposed about the array axis fault plane this finding suggests that the coda consists primarily of waves scattered near the stations rather than waves scattered throughout the crust this result holds even at stations very close to the sources where the amount of coda examined is about 6 times the direct s travel time and it suggests that much of the coseismic velocity decrease associated with the 1989 loma prieta earthquake ellsworth et al 1992 occurred in the shallow crust near the stations one possible mechanism for inducing such a change is a decrease in the mean stress from the mainshock slip however on the basis of static stress modeling we find that this mechanism would produce a velocity increase in the region where velocity was observed to decrease we conclude that the velocity decrease was more likely caused by crack opening or crack coalescence due to strong shaking from the mainshock

seismological studies at parkfield i simultaneous inversion for velocity structure and hypocenters using cubic b splines parameterization a model for the 3d p and s wave velocity structure at parkfield has been developed by the systematic inversion of arrival times for microearthquakes and surface vibroseis sources recorded by nearby calnet stations and by the local borehole installed 10 station high resolution seismic network there a modification of the linear b splines model parameterization method of thurber 1983 to a cubic b splines formulation through tests inverting synthetic data sets for the parkfield configuration is found to produce smooth velocity distributions that still retain some small scale features application of the method to 953 p wave and 288 s wave arrivals from 111 microearthquakes and 8 vibroseis locations throughout the study area yields a velocity structure having variable resolution a low velocity zone appears to coincide with the creeping segment of the san andreas fault zone northwest of the 1966 parkfield earthquake epicenter an increase of about 15 per cent in seismic velocities is seen from the northeast to southwest across this segment of the fault zone the higher velocities to the southwest apparently correspond to the granitic basement the contrast across the fault is subdued to the southeast along the 1966 rupture zone a high vp vs ratio of 1 9 is seen in a 2 km wide volume near the location of the 1966 main shock hypocenter presumably the nucleation zone for the expected m 6 parkfield earthquake

seismic trapped modes in the oroville and san andreas fault zones three component borehole seismic profiling of the recently active oroville california normal fault and microearthquake event recording with a near fault three component borehole seismometer on the san andreas fault at parkfield california have shown numerous instances of pronounced dispersive wave trains following the shear wave arrivals these wave trains are interpreted as fault zone trapped seismic modes parkfield earthquakes exciting trapped modes have been located as deep as 10 kilometers as shallow as 4 kilometers and extend 12 kilometers along the fault on either side of the recording station selected oroville and parkfield wave forms are modeled as the fundamental and first higher trapped sh modes of a narrow low velocity layer at the fault modeling results suggest that the oroville fault zone is 18 meters wide at depth and has a shear wave velocity of 1 kilometer per second whereas at parkfield the fault gouge is 100 to 150 meters wide and has a shear wave velocity of 1 1 to 1 8 kilometers per second these low velocity layers are probably the rupture planes on which earthquakes occur

clustering and periodic recurrence of microearthquakes on the san andreas fault at parkfield california the san andreas fault at parkfield california apparently late in an interval between repeating magnitude 6 earthquakes is yielding to tectonic loading partly by seismic slip concentrated in a relatively sparse distribution of small clusters 20 meter radius of microearthquakes within these clusters which account for 63 of the earthquakes in a 1987 92 study interval virtually identical small earthquakes occurred with a regularity that can be described by the statistical model used previously in forecasting large characteristic earthquakes sympathetic occurrence of microearthquakes in nearby clusters was observed within a range of about 200 meters at communication speeds of 10 to 100 centimeters per second the rate of earthquake occurrence particularly at depth increased significantly during the study period but the fraction of earthquakes that were cluster members decreased

source scaling of earthquakes in the shumagin region alaska time domain inversions of regional waveforms the scaling of pulse duration to seismic moment is estimated for earthquakes along an interplate thrust zone from digital waveforms recorded by short period and broad band instruments of the east aleutian shumagin seismic network we measure pulse duration using an empirical green s function technique based on damped time domain deconvolution from several thousand events 22 earthquakes with magnitudes 3 0 7 0 and depths 23 56 km are found to give reliable estimates of pulse duration durations are also determined directly from one parameter nonlinear inversions for a variety of simple functional forms of source time functions symmetric source pulses boxcar or triangle shapes fit waveforms better than an asymmetric model t exp 2t d for most 62 per cent of the waveform pairs while the asymmetric model fits best for only 8 per cent of the data pulse duration increases with the size of events from 0 1 to 10s over the seismic moment m0 range of 1014 to 3 1019 nm when normalized by the cube root of seismic moment pulse durations show 8 x variation comparable static stress drop estimates range from 0 2 to 135 mpa contrary to predictions of some laboratory and theoretical studies earthquakes at the deepest part of the thrust zone do not show significantly higher stress drops than do shallower events rupture properties however show a strong dependence on earthquake size the three largest events m0 5 1018nm have the three longest normalized durations on average 3 8 times longer than those for smaller events the durations require smaller events to have 10 100 x larger static stress drops or 4 x faster rupture velocities or some combination of the two possibly the largest events rupture both strong and weak patches while smaller events just rupture strong patches on the fault surface the characteristic dimension that separates large from small events 3 15 km is comparable to characteristic wavelengths of pacific basin bathymetry and may reflect the influence of the subducted sea floor upon fault zone heterogeneity

the magnitude frequency distribution of earthquakes recorded with deep seismometers at cajon pass southern california abstract the cumulative b value the slope of the gutenberg richter relationship between earthquake occurrence rate and magnitude is commonly found to be constant 1 network catalogues however reveal a decrease at small magnitudes b value is not just an artifact of catalogue incompleteness but that small earthquakes are really not as numerous as a constant b value extrapolated from larger events would predict in the cajon pass area southern california the b value of seismicity recorded by the local network scsn appears to decrease below about m l 1 6 in order to investigate whether this decrease is real or simply represents the network detection threshold we use seismicity recorded by the deep 1 5 and 2 5 km seismometers deployed in the cajon pass scientific drillhole between april 1992 and october 1994 the maximum amplitudes recorded downhole are compared to scsn magnitudes for events recorded by the network to determine the relationship between amplitude and m l as a function of hypocentral distance from the borehole magnitudes are then calculated for 1300 earthquakes which occurred within 40 km of the borehole magnitude frequency curves are calculated for those events within 18 km of the borehole and a constant b value is observed to m l 0 5

resolution of fault slip along the 470 km long rupture of the great 1906 san francisco earthquake and its implications data from all available triangulation networks affected by the 1906 earthquake have been combined to assess the trade off between slip resolution and its uncertainty and to construct a conservative image of coseismic slip along the rupture because of varying network aperture and station density slip resolution is very uneven although slip is determined within uncertainties of 1 0 m along 60 of the fault constraints are poor on the remaining mostly offshore portions of the rupture slip decreases from maxima of 8 6 and 7 5 m at shelter cove and tomales bay to 4 5 m near mount tamalpais and 2 7 m at loma prieta the geodetically derived slip distribution is in poor agreement with estimates based on analysis of s wave seismograms probably because these waves register only 20 30 of the total seismic moment obtained from longer period surface waves consideration of a range of fault geometries for 1906 slip near loma prieta indicates right lateral motions lie between 2 3 and 3 1 m these values are considerably greater than the 1 5 m of measured surface slip on which several assessments of high earthquake hazard for this fault segment were based this factor along with the absence of 1989 slippage where 1906 surface slip was used to make the forecasts casts doubt on some claims of success in predicting the 1989 m 6 9 loma prieta earthquake

modeling dynamic rupture in a 3d earthquake fault model we propose a fourth order staggered grid finite difference method to study dynamic faulting in three dimensions the method uses an implementation of the boundary conditions on the fault that allows the use of general friction models including slip weakening and rate dependence because the staggered grid method defines stresses and particle velocities at different grid points we preserve symmetry by implementing a two grid row thick fault zone slip is computed between points located at the borders of the fault zone while the two components of shear traction on the fault are forced to be symmetric inside the fault zone we study the properties of the numerical method comparing our simulations with well known properties of seismic ruptures in 3d among the properties that are well modeled by our method are full elastic wave interactions frictional instability rupture initiation from a finite initial patch spontaneous rupture growth at subsonic and supersonic speeds as well as healing by either stopping phases or rate dependent friction we use this method for simulating spontaneous rupture propagation along an arbitrarily loaded planar fault starting from a localized asperity on circular and rectangular faults the shape of the rupture front is close to elliptical and is systematically elongated in the in plane direction of traction drop this elongation is due to the presence of a strong shear stress peak that moves ahead of the rupture in the in plane direction at high initial stresses the rupture front becomes unstable and jumps to super shear speeds in the direction of in plane shear another interesting effect is the development of relatively narrow rupture fronts due to the presence of rate weakening friction the solutions for the thick fault boundary conditions scale with the slip weakening distance do and are stable and reproducible for do greater than about 4 in terms of 2t t ax finally a comparison of scalar and vector boundary conditions for the friction shows that slip is dominant along the direction of the prestress with the largest deviations in slip rate direction occurring near the rupture front and the edges of the fault

seismological studies at parkfield vi moment release rates and estimates of source parameters for small repeating earthquakes abstract waveform data from a borehole network of broadband seismographic stations have been used to study microearthquakes along the parkfield segment of the san andreas fault saf analysis of almost 10 years of such data demonstrates that much of the seismicity in this region consists of repeating sequences quasiperiodic sequences of earthquakes that are essentially identical in terms of waveform size and location scalar seismic moments have been estimated for 53 of these repeating sequences and combined with equivalent estimates from 8 similar but larger event sequences from the stone canyon section of the fault and the main parkfield sequence these estimates show that seismic moment is being released as a function of time in a very regular manner measurements of the moment release rate combined with an assumed tectonic loading rate lead to estimates of the seismic parameters source area slip and recurrence interval such parameters exhibit a systematic dependence upon source size over a range of 10 10 in seismic moment that can be described by three simple scaling relationships several implications of these scaling relationships are explored including the repeat time of earthquakes average stress drop strength of the fault and heat generated by earthquakes what emerges from this analysis of moment release rates is a quantitative description of an earthquake process that is controlled by small strong asperities that occupy less than 1 of the fault area this means that the fault is highly heterogeneous with respect to stress strength and heat generation such heterogeneity helps to explain many of the apparent contradictions that are encountered in the study of earthquakes such as why faults appear weak why significant heat flow is not observed how significant high frequencies can be generated by large earthquakes and how various geologic features such as pseudotachylytes might form

mechanism diversity of the loma prieta aftershocks and the mechanics of mainshock aftershock interaction the diverse aftershock sequence of the 1989 loma prieta earthquake is inconsistent with conventional models of mainshock aftershock interaction because the aftershocks do not accommodate mainshock induced stress changes instead the sense of slip of the aftershocks is consistent with failure in response to a nearly uniaxial stress field in which the maximum principal stress acts almost normal to the mainshock fault plane this orientation implies that i stress drop in the mainshock was nearly complete ii mainshock induced decreases of fault strength helped were important in controlling the occurrence of after shocks and iii mainshock rupture was limited to those sections of the fault with preexisting shear stress available to drive fault slip

effects of a low velocity zone on a dynamic rupture dynamic crack earthquake simulations generally assume that the crustal material surrounding faults is laterally homogeneous tomographic and near fault seismic studies indicate that the crust near faults is instead comprised of rocks of varying material velocities we have tested the effects of adding material velocity variation to simulations of spontaneously propagating earthquakes we used two dimensional plane strain conditions coupled with a slip weakening fracture criterion and examined earthquakes on faults that bisect finite width low velocity zones embedded in country rock and earthquakes on faults that bound two different velocity materials when a fault bisects a low velocity zone the normal stress remains unchanged but both the rupture velocity and slip velocity pulse shape are perturbed the presence of the low velocity zone induces high frequency oscillations in the slip function near the rupture front when the fault is on the edge of the low velocity zone the oscillations are more pronounced and repeated sticking and slipping can occur near the rupture front for the slip weakening velocity independent friction model however the temporary sticking does not lead to permanent arrest of slip and slip duration is still controlled by the overall rupture dimension when an earthquake ruptures a fault juxtaposing a lower velocity material against a higher velocity material the normal stress across the fault near the crack tip is perturbed the sign of the normal stress perturbation depends on the direction of rupture leading in some cases to a directional dependence of rupture velocity when slip is accompanied by stress reduction a positive feedback develops between the normal and shear stress changes as previously noted by andrews and ben zion 1997 resulting in an apparently unavoidable grid size dependence in computation of stress change near the rupture front numerical experiments indicate however that the rupture velocity is insensitive to this zone size dependence which is highly localized immediately behind the crack tip the factors controlling the rupture velocity in the simulations including directional dependence are further elucidated by a new analytical solution for rupture of an asperity on a frictionless interface

discriminating characteristics of tectonic and human induced seismicity we analyze statistical features of background and clustered subpopulations of earthquakes in different regions in an effort to distinguish between human induced and natural seismicity analysis of end member areas known to be dominated by human induced earthquakes the geyser geothermal field in northern california and tautona gold mine in south africa and regular tectonic activity the san jacinto fault zone in southern california and the coso region excluding the coso geothermal field in eastern central california reveals several distinguishing characteristics induced seismicity is shown to have 1 higher rate of background events both absolute and relative to the total rate 2 faster temporal offspring decay 3 higher rate of repeating events 4 larger proportion of small clusters and 5 larger spatial separation between parent and offspring compared to regular tectonic activity these differences also successfully discriminate seismicity within the coso and salton sea geothermal fields in california before and after the expansion of geothermal production during the 1980s online material figures of examined regions earthquake epicenters and joint distributions of the time and space components of the nearest neighbor distance

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

implications for mechanical properties of brittle faults from observations of the punchbowl fault zone california field observations of the punchbowl fault zone an inactive trace of the san andreas are integrated with results from experimental deformation of naturally deformed punchbowl fault rocks for a qualitative description of the mechanical properties of the fault and additional information for conceptual models of crustal faulting the punchbowl fault zone consists of a single continuous gouge layer bounded by zones of extensively damaged host rock fault displacements were not only localized to the gouge layer but also to discrete shear surfaces within the gouge deformation in the exposure studied probably occurred at depths of 2 to 4 km and was dominated by cataclastic mechanisms textural data also suggest that significant amounts of pore fluids were present during faulting and that fluid assisted mechanisms such as dissolution diffusion and precipitation were operative the experimental data on specimens collected from the fault zone suggest that there is a gradual decrease in strength and elastic modulus and an increase in relative ductility and permeability toward the main gouge zone the gouge layer has fairly uniform mechanical properites and it has significantly lower strength elastic modulus and permeability than both the damaged and the undeformed host rock for the punchbowl fault and possibly other brittle faults the variations in loading of the gouge zone with time are primarily governed by the morphology of the fault and the mechanical properties of the damaged host rock in addition the damaged zone acts as the permeable unit of the fault zone and surrounding rock it appears that the gouge primarily governs whether displacements are localized and it therefore may have a significant influence on the mode of slip

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

the vorticity balance in the tropical upper troposphere of a general circulation model abstract the time mean vorticity balance in the summertime tropical upper troposphere of an atmospheric general circulation model constructed at the geophysical fluid dynamics laboratory is examined with particular emphasis on the detailed balance in the tibetan anticyclone the model produces a reasonable simulation of the large scale features of the northern summer 200 mb flow in the tropics without the inclusion of subgrid scale processes that strongly damp the upper tropospheric vorticity the vorticity balance is essentially nonlinear and nearly inviscid there is considerable cancellation between the stretching and horizontal advection of vorticity by the time mean flow in the vicinity of the tibetan anticyclone with much of the remainder balanced by vertical advection and twisting mixing by the resolved transients is not negligible in some regions but considerably smaller than the horizontal advection overall and less well correlated with the stretching subgrid scale mixing consisting only o

postseismic response of repeating aftershocks the recurrence intervals of repeating earthquakes on the san andreas fault in the loma prieta aftershock zone follow the characteristic 1 t decay of omori s law a model in which these earthquakes occur on isolated patches of the fault that fail in stick slip with creep around them can explain this observation in this model the recurrence interval is inversely proportional to the loading rate due to creep logarithmic velocity strengthening friction predicts 1 t decay in creep rate following the mainshock the time dependence of recurrence is inconsistent with a simple viscous constitutive relationship which predicts an exponential decay of loading rate thus our observations imply postseismic slip at seismogenic depth under a power law rheology the time dependence of postseismic deformation measured geodetically may be diagnostic of whether postseismic deformation is caused by creep or possible viscoelastic deformation at greater depths

variations in rupture process with recurrence interval in a repeated small earthquake in theory and in laboratory experiments friction on sliding surfaces such as rock glass and metal increases with time since the previous episode of slip1 this time dependence is a central pillar of the friction laws widely used to model earthquake phenomena2 3 on natural faults other properties such as rupture velocity4 5 porosity and fluid pressure6 11 may also vary with the recurrence interval eighteen repetitions of the same small earthquake separated by intervals ranging from a few days to several years allow us to test these laboratory predictions in situ the events with the longest time since the previous earthquake tend to have about 15 larger seismic moment than those with the shortest intervals although this trend is weak in addition the rupture durations of the events with the longest recurrence intervals are more than a factor of two shorter than for the events with the shortest intervals both decreased duration and increased friction are consistent with progressive fault healing during the time of stationary contact