granite ascent in convergent orogenic belts testing a model the common spatial relationship in convergent orogenic belts between a crustal scale shear zone system high grade metamorphic rocks and granites suggests a feedback relation between crustal anatexis and contractional deformation that helps granite extraction and focuses granite ascent such a feedback relation has been proposed for ascent of early devonian granites in west central maine this interpretation requires that deformation metamorphism and plutonism were synchronous we have determined precise u pb zircon and monazite ages that we interpret to record time of crystallization of syntectonic granite in metric to decametric sheets and kilometric plutons and of schlieric granite within migmatites ages are in the range ca 408 404 ma within 1 m y at 95 confidence limits these ages are similar to extant u pb monazite ages of ca 405 399 2 ma for syntectonic regional metamorphism in the same area the coincidence between the age of peak metamorphism and crystallization ages of granite shows tectonics metamorphism and magmatism were contemporaneous in support of the feedback model

determination of the location of faulting beneath kobe during the 1995 hyogo ken nanbu japan earthquake from near source particle motion we constrained the location of faulting beneath kobe during the 1995 hyogo ken nanbu kobe japan earthquake by investigating the particle motion of strong motion records from sites close to the causative fault the intersection of the earth s surface and the extension of buried fault plane is consistent with known active fault lines but not with the severely damaged zone the intersection line also agrees with geodetic data

coseismic fault slip associated with the 1992 m w 6 1 joshua tree california earthquake implications for the joshua tree landers earthquake sequence coseismic surface deformation associated with the mw 6 1 april 23 1992 joshua tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from global positioning system and trilateration measurements the rms signal to noise ratio for these inferred displacements is 1 8 with near fault displacement estimates exceeding 40 mm in order to determine the long wavelength distribution of slip over the plane of rupture a tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right lateral slip and zero surface slip constraints the resulting slip distribution yields a geodetic moment estimate of 1 7 1018 n m with corresponding maximum slip around 0 8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations from empirical green s function analyses a rupture duration of 5 s is obtained which implies a rupture radius of 6 8 km most of the inferred slip lies to the north of the hypocenter consistent with northward rupture propagation stress drop estimates are in the range of 2 4 mpa in addition predicted coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0 1 mpa in contrast predicted stress changes are near zero at the hypocenter of the mw 7 3 june 28 1992 landers earthquake which nucleated about 20 km beyond the northernmost edge of the joshua tree rupture based on aftershock migrations and the predicted static stress field we speculate that redistribution of joshua tree induced stress perturbations played a role in the spatio temporal development of the earthquake sequence culminating in the landers event

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

rupture process of the 1987 superstition hills earthquake from the inversion of strong motion data a pair of significant earthquakes occurred on conjugate faults in the western imperial valley involving the through going superstition hills fault and the elmore ranch cross fault the first event was located on the elmore ranch fault m s 6 2 and the larger event on the superstition hills fault m s 6 6 the latter event is seen as a doublet teleseismically with the amplitudes in the ratio of 1 2 and delayed by about 8 sec this 8 sec delay is also seen in about a dozen strong motion records these strong motion records are used in a constrained least squares inversion scheme to determine the distribution of slip on a 2 d fault upon closer examination the first of the doublets was found to be itself complex requiring two episodes of slip thus the rupture model was allowed to have three separate subevents treated as separate ruptures with independent locations and start times the best fits were obtained when all three events initiated at the northwestern end of the fault near the intersection of the cross fault their respective delays are 2 1 and 8 6 sec relative to the first subevent and their moments are 0 4 0 9 and 3 5 10 25 dyne cm which is about half of that seen teleseismically this slip distribution suggests multi rupturing of a single asperity with stress drops of 60 200 and 15 bars respectively the first two subevents were confined to a small area around the epicenter while the third propagated 18 km southwestward compatible with the teleseismic and afterslip observations

inversion for slip distribution using teleseismic p waveforms north palm springs borah peak and michoacan earthquakes we have inverted the teleseismic p waveforms recorded by stations of the global digital seismograph network for the 8 july 1986 north palm springs california the 28 october 1983 borah peak idaho and the 19 september 1985 michoacan mexico earthquakes to recover the distribution of slip on each of the faults using a point by point inversion method with smoothing and positivity constraints in the inversion procedure a fault plane with fixed strike and dip is placed in the region of the earthquake hypocenter and divided into a large number of subfaults rupture is assumed to propagate at a constant velocity away from the hypocenter and synthetic ground motions for pure strike slip and dip slip dislocations are calculated at the teleseismic stations for each subfault the observed seismograms are then inverted to obtain the distribution of strike slip and dip slip displacement for the earthquake results of the inversion indicate that the global digital seismograph network data are useful for deriving fault dislocation models for moderate to large events however a wide range of frequencies which includes periods shorter than those within the passband of the long period global digital seismograph network instruments is necessary to infer the distribution of slip on the earthquake fault although the long period waveforms define the size dimensions and seismic moment of the earthquake data at shorter periods provide additional constraints on the variation of slip on the fault dislocation models obtained for all three earthquakes are consistent with a heterogeneous rupture process where failure is controlled largely by the size and location of high strength asperity regions

fractals in geology and geophysics the definition of a fractal distribution is that the number of objectsn with a characteristic size greater thanr scales with the relationn r d the frequency size distributions for islands earthquakes fragments ore deposits and oil fields often satisfy this relation fractals were originally introduced by mandelbrot to relate the length of a coastline to the length of the measuring stick this application illustrates a fundamental aspect of fractal distributions scale invariance the requirement of an object to define a scale in photographs of many geological features is one indication of the wide applicability of scale invariance to geological problems scale invariance can lead to fractal clustering geophysical spectra can also be related to fractals these are self affine fractals rather than self similar fractals examples include the earth s topography and geoid

comparison of seismic waveform inversion results for the rupture history of a finite fault application to the 1986 north palm springs california earthquake the july 8 1986 north palm springs earthquake is used as a basis for comparison of several different approaches to the solution for the rupture history of a finite fault the inversion of different waveform data is considered both teleseismic p waveforms and local strong ground motion records linear parametrizations for slip amplitude are compared with nonlinear parametrizations for both slip amplitude and rupture time inversions using both synthetic and empirical green s functions are considered in general accurate green s functions are more readily calculable for the teleseismic problem where simple ray theory and flat layered velocity structures are usually sufficient however uncertainties in the variation in t with frequency most limit the resolution of teleseismic inversions a set of empirical green s functions that are well recorded at teleseismic distances could avoid the uncertainties in attenuation in the inversion of strong motion data the accurate calculation of propagation path effects other than attenuation effects is the limiting factor in the resolution of source parameters the assumption of a laterally homogeneous velocity structure is usually not a good one and the use of empirical green s functions is desirable considering the parametrization of the problem any degree of fault rupture complexity can be described in terms of a linear parametrization for slip amplitudes however a nonlinear parametrization for rupture times and slip amplitudes can have a distinct advantage over a simple linear one by limiting the number of unknown parameters regardless of the choice of data or the type of parametrization the model or solution will be affected by the choice of minimization norm and the type of stabilization used

a spatial random field model to characterize complexity in earthquake slip 1 finite fault source inversions reveal the spatial complexity of earthquake slip over the fault plane we develop a stochastic characterization of earthquake slip complexity based on published finite source rupture models in which we model the distribution of slip as a spatial random field the model most consistent with the data follows a von karman autocorrelation function acf for which the correlation lengths a increase with source dimension for earthquakes with large fault aspect ratios we observe substantial differences of the correlation length in the along strike a x and downdip a z directions increasing correlation length with increasing magnitude can be understood using concepts of dynamic rupture propagation the power spectrum of the slip distribution can also be well described with a power law decay i e a fractal distribution in which the fractal dimension d remains scale invariant with a median value d 2 29 0 23 while the comer wave number k c which is inversely proportional to source size decreases with earthquake magnitude accounting for larger slip patches for large magnitude events our stochastic slip model can be used to generate realizations of scenario earthquakes for near source ground motion simulations

global positioning system measurements of deformations associated with the 1987 superstition hills earthquake evidence for conjugate faulting large station displacements observed from imperial valley global positioning system gps campaigns are attributed to the november 24 1987 superstition hills earthquake sequence thirty sites from a 42 station gps network established in 1986 have been reoccupied during 1988 and or 1990 displacements at three sites within 3 km of the surface rupture approach 0 5 m eight additional stations within 20 km of the seismic zone are displaced at least 10 cm this is the first occurrence of a large earthquake ms 6 6 within a preexisting gps network best fitting uniform slip models of rectangular dislocations in an elastic half space indicate 130 8 cm right lateral displacement along the northwest trending superstition hills fault and 30 10 cm left lateral displacement along the conjugate northeast trending elmore ranch fault the geodetic moments are 9 4 1025 dyn cm and 2 3 1025 dyn cm for the superstition hills and elmore ranch faults respectively consistent with teleseismic source parameters the data also suggest that postseismic slip along the superstition hills fault is concentrated at shallow depths distributed slip solutions using singular value decomposition indicate near uniform displacement along the elmore ranch fault and concentrated slip to the northwest and southeast along the superstition hills fault a significant component of nonseismic secular displacement is observed across the imperial valley which is attributed to interseismic plate boundary deformation

high resolution images of the lower crust deep seismic reflections from 15 to 180 hz abstract we report the results of an unconventional broad band deep seismic reflection survey undertaken on land in southwest england previous marine seismic lines in this area have shown strong layered lower crustal reflectivity the broad band survey shows near normal incidence p wave reflections from the lower continental crust and from the moho with frequencies that range from 15 to 180 hz high frequency lower crustal s wave reflections were also observed these high frequencies were obtained by placing both sources and receivers in boreholes up to 2 km deep drilled within a large homogeneous granite pluton recording down hole avoids loss of high frequency energy during transmission through heterogeneous and highly attenuative near surface layers the high frequencies were generated using explosive sources and were recorded on sub surface accelerometers and hydrophones and on a conventional array of surface geophones seismic attenuation within the granite is low the quality factor q lies above 2000 these new high frequency observations show that the lower crust does not preferentially reflect particular frequencies the data are most readily explained as resulting from the combination of a frequency independent q in the upper crust with a white reflectivity response in the lower crust the data are not consistent with crustal models which are highly tuned to reflect preferentially the frequencies employed in conventional deep seismic reflection surveys the data require relatively sharp boundaries within the lower crust significant variation in physical properties must occur over vertical distances of not more than about 10 m sharp boundaries are most easily explained if the reflections result from lithological contrasts

subevent structure of large earthquakes a ground motion perspective our ability to predict ground motions from future earthquakes hinges on the accurate modeling of the radiation from the earthquake source a successful approach to modeling earthquake radiation has been to represent faults as a series of discrete independently rupturing subfaults although their physical interpretation has remained largely unclear our simulation of ground motions from twenty six well recorded moderate to large earthquakes substantiates the hypothesis that large ruptures are made up of a sequence of smaller subevents the size of a subevent follows a simple linear relationship with the size of the earthquake fault in an apparently deterministic manner the strength of the high frequency radiation is controlled by the maximum slip velocity which varies stochastically over a small range both the characteristic subevent size and the slip velocity appear to be region independent indicating a remarkable uniformity in earthquake source properties this observation helps to reduce the uncertainty in ground motion prediction by constraining its essential and previously ambiguous parameters

linearized inversion for fault rupture behavior application to the 1984 morgan hill california earthquake we present a technique to infer the rupture history of an earthquake from near source records of ground motion unlike most previous studies each point on the fault is assumed to slip only once when the rupture front passes with a spatially variable slip intensity in this parameterization the data are linearly related to slip intensity but nonlinearly related to rupture time we perform a linearized inversion for slip intensity and rupture time by iteratively perturbing an assumed starting model the inverse problem for the model perturbation is solved using a tomographic back projection technique smoothing and inequality constraints are applied to ensure that the resulting solution is stable and feasible asymptotic ray theory is used to calculate theoretical seismograms and partial derivatives with respect to model parameters in test cases with noisy synthetic data sets we found that it is possible to distinguish between rupture models having variable slip amplitude and models having variable rupture velocity provided the station coverage is adequate we apply the technique to recordings of the april 24 1984 morgan hill earthquake the results indicate that slip amplitude on the fault plane was extremely variable and that the rupture front did not propagate uniformly away from the hypocenter in particular rupture was delayed on a 12 km2 section of the fault approximately 14 km to the southeast of the hypocenter the rupture front surrounded the region which subsequently failed with a component of rupture propagation back toward the hypocenter similar behavior has been observed in dynamic rupture models with stress or strength inhomogeneities this segment of the fault ruptured with a large slip amplitude releasing 12 of the total seismic moment from 4 of the total area of the aftershock zone the surface trace of this section of the fault is characterized by a complex left step that could act to increase the normal stress acting across the fault however the distribution of aftershocks suggests that the fault at depth is simpler and that it may bend to the right in either case our rupture model suggests that this segment of the fault represents an asperity which initially resisted rupture but eventually ruptured massively we estimate the shear fracture energy for this earthquake to be 2 106 j m2

spatial variability of seismicity parameters in aftershock zones the spatial variability of the b value of the frequency magnitude relationship and the decay rate of aftershocks as described by the p value of the modified omori law is investigated by using dense spatial grids we map out the distribution of b and p values within the landers northridge morgan hill and kobe aftershock sequences considerable spatial variability is found with b values of independent subvolumes ranging from 0 6 to 1 4 and p values ranging from 0 6 to 1 8 these systematic and statistically highly significant differences argue that it is an oversimplification to assign one single p and b value to an aftershock sequence that extends up to 100 km the spatial distribution of these two parameters is compared with the slip distribution during the mainshock suggesting that the areas of largest slip release correlate with high b value regions we hypothesize that the frictional heat created during the event may influence the p value distribution within an aftershock zone while applied shear stress crack density and pore pressure govern the frequency magnitude distribution by investigating the frequency magnitude distribution separately for preseismic and postseismic periods for the morgan hill mainshock we find that only the volume in the vicinity of the highest slip release shows a significant increase in the b value which decays to premainshock values within a year surrounding areas of the aftershock zone show an approximately constant b value with time because the aftershock hazard after a mainshock depends strongly on both the b and p value we propose that aftershock hazard assessment can be improved by taking into account the spatial distribution of the parameters

simultaneous inversion for earthquake location and velocity structure beneath central costa rica we have imaged the complex crustal and upper mantle structure beneath central costa rica using p wave arrival times from locally recorded earthquakes thurber s 1983 iterative inversion method is used to simultaneously estimate velocities along a three dimensional grid and hypocentral parameters of local earthquakes our data consist of over 12 000 arrival times from more than 1300 earthquakes recorded by stations of a permanent seismographic network in costa rica our resulting velocity model correlates well with mapped geologic units at very shallow depth and with tectonic features at greater depth we find low velocities 4 0 to 4 8 km sec in the shallow crust above 10 km near the active volcanoes and associated with a nw se trending late cretaceous to late tertiary sedimentary basin southeast of herradura peninsula high velocities 5 4 to 5 7 km sec in the shallow crust correlate with outcrops of late jurassic to early tertiary ultramafic ophiolitic units and with basic tertiary volcanic units at depths between 20 and 30 km high velocities 6 8 to 7 2 km sec are associated with the subducting cocos plate under costa rica and two prominent low velocity bodies 6 3 to 6 5 km sec are present about 30 km trenchward of the volcanic arc and along the projection of the aseismic cocos ridge as it subducts beneath costa rica the thickened oceanic crust of the cocos ridge is most likely responsible for its low velocities the deep low velocity anomaly located trenchward of the axis of the volcanoes may indicate the presence of a low density intrusive resulting from an earlier phase of magmatism possibly the late miocene episode that produced the talamanca intrusive complex

how strong is an asperity a recent study of repeating earthquakes on the san andreas fault in central california by nadeau and johnson 1998 found that the smallest events occurred on patches having a linear dimension of the order of 0 5 m displacements of about 2 cm and stress drops of the order of 2000 mpa roughly 10 times larger than rock strengths measured in the laboratory the stress drop for larger events was observed to decrease as a power law of the seismic moment reaching the commonly observed value of 10 mpa at about magnitude 6 these large strengths are shown here to be consistent with laboratory data if the preexisting microcracks are all healed a hierarchical fractal asperity model is presented which is based on recent laboratory observations of contact distributions in sliding friction experiments this cantor dust model is shown to be consistent with the observed power law decrease in stress drop and increase in displacement with increasing event size the spatial distribution of hypocenters in the parkfield area is shown to be consistent with this simple fractal model and with a hierarchical clustering of asperities having a fractal dimension of d 1 and discrete rescaling factor of about 20

a composite source model for computing realistic synthetic strong ground motions a composite source model is presented for convolution with synthetic green s functions in order to synthesize strong ground motions due to a complex rupture process of a large earthquake subevents with a power law distribution of sizes are located randomly on the fault each subevent radiates a displacement pulse with the shape of a brune s pulse in the far field at a time determined by a constant rupture velocity propagating from the hypocenter thus all the input parameters have a physical basis we simulate strong ground motions for event station pairs that correspond to records obtained in mexico by the guerrero accelerograph network the synthetic accelerations velocities and displacements have realistic amplitudes durations and fourier spectra

a hybrid method for calculating near source broadband seismograms application to strong motion prediction abstract we present a hybrid method for computing broadband strong motion seismograms in the near field of large earthquakes we combine complete seismograms at low frequency with ray theory seismograms at high frequency to form a composite broadband seismogram that spans the entire frequency range of interest in our approach the amplitude spectra of the two sets of synthetic seismograms are reconciled at intermediate frequencies where their domain of validity overlaps we demonstrate the method with scenario earthquakes based on the spatial random field model for complex earthquake slip j geophys res 107 b11 2002 2308 the hybrid near source broadband seismograms are useful both for detailed source modeling and for incorporating source effects into probabilistic seismic hazard analysis

source parameters of the 23 april 1992 m 6 1 joshua tree california earthquake and its aftershocks empirical green s function analysis of geos and terrascope data source parameters of the m 6 1 23 april 1992 joshua tree mainshock and 86 m 1 8 to 4 9 aftershocks are determined using an empirical green s function methodology for the aftershocks deconvolved p and s wave spectra are calculated for 126 pairs of closely spaced events recorded on portable geos stations s wave spectra from the two horizontal components are averaged the deconvolved spectra are fit by a ratio of omega square source models yielding an optimal least squares corner frequency for both the large and the small event in each pair we find no resolved difference between the inferred p and s wave corner frequencies using the standard brune model for stress drop we also find no resolved nonconstant scaling of stress drop with moment although we also conclude that detailed scaling systematics would be difficult to resolve in particular a weak increase of stress drop with moment over a limited moment magnitude cannot be ruled out for magnitudes smaller than m 3 to 3 5 the inferred stress drop values will be limited by the maximum observable corner frequency value of 60 hz for the mainshock source time functions are obtained from mainshock recordings at three terrascope stations pfo pas and gsc using an m 4 3 foreshock as an empirical green s function the results indicate a fairly simple single pulse source time function with clear south to north directivity and an inferred rupture radius of 5 to 6 km the deconvolved source time functions are inverted to obtain a finite rupture model that gives a robust estimate of rupture dimension early aftershocks are found to lie along the perimeters of regions with high mainshock slip the inferred mainshock stress drop value 56 bars is within the range determined for the aftershocks our derived mainshock source spectra do not show resolvable deviation from the omega square model

seismic imaging of the 1999 izmit turkey rupture inferred from the near fault recordings we use near fault accelerograms to infer the space time history of rupture on the fault during the izmit earthquake the records show that the ground displacement and velocity near the fault were surprisingly simple rupture propagated toward the west at a velocity of about 3 km s and toward the east at a remarkably high average velocity of 4 7 km s over a distance of about 45 km before decelerating to about 3 1 km s on the eastern segment slip on the fault is particularly large down to a depth of 20 km on the central portion of the fault where it reaches about 7 m slip is large also below 10 km on the eastern fault segment and this may have contributed to the loading of shear stress on the duzce fault on the western fault segment large slip seems confined to shallow depths

fault geometry at the rupture termination of the 1995 hyogo ken nanbu earthquake the source geometry and slip distribution at rupture termination of the 1995 hyogo ken nanbu earthquake were investigated using waveform inversion on the assumption of fault branching in the northeastern part of the rupture model possible branching of the okamoto fault is suggested both by the static displacement distribution and damage extension east of kobe nishinomiya area to exclude data contaminated by the basin edge diffracted wave in the waveform inversion process we examined the spatio temporal variation of its influence and from a comparison with a flat model we determined windows appropriate for the data three subevents were identified as was reported in previous works the largest was in the shallow part on the awaji side whereas the smaller two occurred at great depths 7km on the kobe side we found a smaller subevent in the deep part of the branch total variance reduction was larger and the abic value was smaller when we assumed the branch than when we did not which shows the superiority of the branching fault model resolution checks showed that the slips on proposed branched portions are physical and not caused by random data noise or systematic errors in green s func tions that arise from misestimation of velocity structures calculation of the relative static displacement between two leveling observation stations near the branching fault also showed the greater utility of the branching fault model the effect of slip on the branched fault has been shown in near source ground motion simulation using the 3d finite difference method iwata et al 1999 the characteristic distribution of ground motion in the near source region indicated by the damage distribution is well reproduced by the modeling of both the source process and wave propagation in the realistic 3d velocity structure slip on the branched fault affected the ground motion in eastern kobe nada and higashi nada wards ashiya and nishinomiya cities but its contribution is not dominant even in those regions about 30 to 50 maximum velocity in the frequency range of 0 1 to 1 0 hz we conclude that branching rupture of the okamoto fault during the hyogo ken nanbu earthquake is the preferable interpretation on the source process of this earth quake this supports the idea of a geometrical barrier suggested by geological ob servations

the 1971 san fernando earthquake a double event evidence is presented which suggest that the 1971 san fernando earthquake may have been a double event that occurred on two separate subparallel thrust faults it is postulated that the initial event took place at depth on the sierra madre fault zone which runs along the base of the san gabriel mountains rupture is postulated to have occurred from a depth of about 15 km to a depth of about 3 km a second event is thought to have initiated about 4 sec later on another steeply dipping thrust fault which is located about 4 km south of the sierra madre fault zone the surface trace of this fault coincides with the san fernando fault zone which was the principal fault associated with surface rupture it is postulated that rupture propagated from a depth of 8 km to the free surface the moments of the first and second events are approximately 0 7 10 26 dyne cm and 1 0 10 26 dyne cm respectively this model is found to explain the combined data sets of strong ground motions teleseismic p and s waveforms and static offsets better than previous models which consist of either a single fault plane or a plane having a dip angle which shallows with decreasing depth nevertheless many features of the observed motions remain unexplained and considerable uncertainty still exists regarding the faulting history of the san fernando earthquake

high frequency spectral falloff of earthquakes fractal dimension of complex rupture b value and the scaling of strength on faults the high frequency falloff omega sup minus gamma of earthquake displacement spectra and the b value of aftershock sequences are attributed to the character of spatially varying strength along fault zones the author assumes that the high frequency energy of a main shock is produced by a self similar distribution of subevents where the number of subevents with radii greated than r is proportional to r sup minus d d being the fractal dimension in this model an earthquake is composed of a hierarchical set of smaller earthquakes the static stress drop is parameterized to be proportional to r sup eta and strength is assumed to be proportional to static stress drop he finds that a distribution of subevents with d 2 and stress drop independent of seismic moment eta 0 produces a main shock with an omega sup minus 2 falloff if the subevent areas fill the rupture area of the main shock by equating subevents to islands of high stress of a random self similar stress field on a fault he relates d to the scaling of strength on a fault such that d 2 minus eta thus d 2 corresponds to constant stress drop scaling eta 0 and scale invariant fault strength amore self similar model of aftershock rupture zones on a fault is used to determine the relationship between the b value the size distribution of aftershock rupture zones and the scaling of strength on a fault the b value for aftershock sequences on a fault is found to equal 3 minus 1 5 eta 3 eta therefore this model indicates that the typically observed spectral falloffs of omega sup minus 2 and b values of 1 can be entirely caused by scale invariant strength eta 0 along fault zones less

comparison of the coseismic rupture with the aftershock distribution in the hyuga nada earthquakes of 1996 on october 19 1996 a large underthrusting earthquake ms 6 7 the hyuga nada japan earthquake occurred along the southern end of nankai trough about two months later a second large earthquake ms 6 7 occurred in the adjacent region we study the source process of the two large earthquakes in the hyuga nada region and compare the coseismic rupture area with aftershock distribution the main source parameters obtained for the first mainshock are strike dip rake 210 12 81 the seismic moment mo 2 3 1019 nm mw 6 8 the rupture area s 20 15 km and the source duration t 17 s for the second main shock strike dip rake 210 12 87 the seismic moment mo 1 5 1019 nm mw 6 7 the rupture area s 18 18 km and the source duration t 15 s the coseismic rupture areas do not overlap the aftershock area while the aftershock areas of the two main shocks mutually overlap this implies that the common aftershock region takes a role of barriers to dynamic rupture it is also seen that the aftershock area expanded during the first one day

application of an iterative least squares waveform inversion of strong motion and teleseismic records to the 1978 tabas iran earthquake an iterative least squares technique is used to simultaneously invert the strong motion records and teleseismic p waveforms for the 1978 tabas iran earthquake to deduce the rupture history the method solves for both slip amplitudes and rupture times on a finite fault several inversions of the data are done to determine the best fitting hypocenter average rupture velocity and source time function within the range of our parameter search the effects of using different data sets and different parametrizations of the problem linear versus nonlinear are also considered teleseismic body waves are shown to be essential to obtain a reliable solution due to the sparse coverage of the fault plane by strong motion stations a simple triangular source time function is compared with a modified kostrov source time function for a circular rupture and a constant propagating stress drop although the waveforms used in the inversions can not resolve between these two time functions the kostrov source time function gives a moment estimate and maximum fault displacement more consistent with long period surface waves and field observations the results of the iterative inversion for slip amplitudes and rupture times are compared with the results of a linear inversion for slip amplitudes in which the rupture velocity is fixed flexibility in the sequence of faulting for the linear inversion is obtained by allowing multiple rupture intervals a consensus of all the inversion runs indicates a complex multiple source for the tabas earthquake with four main source regions over a fault length of 90 km and an average rupture velocity of 2 5 km sec

sylvite and halite on particles recovered from 25143 itokawa a preliminary report we observed cross sectional ultra thin sections near the surface of 12 particles recovered from the s type asteroid itokawa by the hayabusa spacecraft in 2010 using spherical aberration corrected stem and conventional tem although their mineralogy is almost identical to the equilibrated ll chondrites and therefore basically anhydrous micrometer to submicrometer sized sylvite was identified on the surface of itokawa particle ra qd02 0034 separately micrometer sized halite was also identified on the surface of itokawa particle ra qd02 0129 detailed inspection of the sample processing procedures at the jaxa s planetary materials sample curation facility and textural observation of the sylvite and halite indicate that they were clearly present on two itokawa particles before they were removed from clean chamber 2 at jaxa however there is no direct evidence for their extraterrestrial origin at present if the sylvite and halite are extraterrestrial their presence suggests that they may be more abundant on the surface of s type asteroids than previously thought

source rupture process of the kocaeli turkey earthquake of august 17 1999 obtained by joint inversion of near field data and teleseismic data the rupture process of the 1999 turkey earthquake is examined using both near field strong motion data and teleseismic body wave data the derived source parameters are as follows strike dip slip 268 86 180 nearly pure strike slip the seismic moment mo 1 7 1020 nm mw 7 4 the source duration 20 sec the fault length 70 km the fault width 15 km the rupture process is characterized by an asymmetric bilateral rupture propagation and smooth slip it consists of two major fault segments a rupture propagating to the west and a second rupture propagating to the east the maximum dislocation and the maximum dislocation velocity are 6 3 m and 2 7 m s respectively both found at the former segment the average dislocation is about 4 m the extent of the coseismic rupture suggests that a considerable part of the anticipated seismic gap remains unruptured

source complexity of the 1987 whittier narrows california earthquake from the inversion of strong motion records strong motion records for the whittier narrows earthquake are inverted to obtain the history of slip both constant rupture velocity models and variable rupture velocity models are considered the results show a complex rupture process within a relatively small source volume with at least four separate concentrations of slip two sources are associated with the hypocenter the larger having a slip of 55 90 cm depending on the rupture model these sources have a radius of approximately 2 3 km and are ringed by a region of reduced slip the aftershocks fall within this low slip annulus other sources with slips from 40 to 70 cm each ring the central source region and the aftershock pattern all the sources are predominantly thrust although some minor right lateral strike slip motion is seen the overall dimensions of the whittier earthquake from the strong notion inversions is 10 km long along the strike and 6 km wide down the dip the preferred dip is 30 degree and the preferred average rupture velocity is 2 5 km s moment estimates range from 7 4 to 10 0 times 10 sup 24 dyn cm depending on the rupture model

fractal geometry in the san andreas fault system it has been noted that the spatial distribution of earthquakes and the mode of strain release in the san andreas fault system is related to the complexity of fault geometry because of their rough appearance over many length scales faults can be regarded as fractal surfaces direct estimates of fractal dimension d of portions of the san andreas fault system between the northern gabilan range and the salton sea including the postulated extent of the great 1857 fort tejon earthquake are obtained from measured fault lengths analogous to the lengths of coastlines as discussed by mandelbrot regions characterized by complicated fault geometry are associated with larger values of d based on fault traces mapped at a scale of 1 750 000 d is 1 3 for this reach of the fault defined as a 30 km wide band about a main fault trace for that part near parkfield which could be associated with the nucleation of the 1857 earthquake d is 1 1 at this same scale d is 1 4 for the san andreas and related faults near san bernardino where the 1857 rupture stopped compared to 1 2 for the san andreas san juan fault segments near the point of arrest of the 1966 parkfield earthquake at finer map scales 1 24 000 and 1 62 500 critical lengths of 500 m and 1 km are identified which might relate to the extent of off san andreas fault offsets the critical lengths also suggest that fault geometry is not self similar if this fractal geometry persists through the seismic cycle it may be possible to use a quantitative measure of complexity to explain the occurrence of great and characteristic earthquakes along a given reach of fault

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

strong motion and broadband teleseismic analysis of the 1991 sierra madre california earthquake short period and broadband teleseismic waveform data and three component strong motion records were analyzed to obtain the source parameters of the 1991 sierra madre earthquake close in strong motion velocity records analyzed from 5 s to 5 hz show two distinct pulses about 0 35 s apart requiring some rupture complexity the near field shear wave displacement pulse from this event has a relatively short duration about 1 s for the magnitude of the event requiring a particularly high average stress drop to further constrain the rupture process the data were used in a finite fault source inversion to determine the temporal and spatial distribution of slip we chose a fault plane orientation striking s 62 w and dipping 50 toward the northwest as required by the distribution of aftershocks the first motion mechanism and the teleseismic bodywave point source inversion in addition to the aftershock locations depth constraints are provided by teleseismic short period and broadband recordings which require a centroid depth of 10 11 km our inverse modeling results indicate that both the teleseismic and strong motion data sets can be fit with a compact rupture area about 12 km2 southwest and up dip from the hypocenter the average slip is approximately 50 cm and the maximum slip is 120 cm the seismic moment obtained from either of the separate data sets or both sets combined is about 2 8 0 3 1024 dyn cm and the potency is 0 01 km3 using the area of significant slip estimated from the finite fault inversion the resulting stress drop is on the order of 150 200 bars

stochastic modeling of seafloor morphology inversion of sea beam data for second order statistics at scale lengths less than 100 km or so statistical descriptions of seafloor morphology can be usefully employed to characterize ridge crest processes off ridge tectonics and vulcanism sedimentation and postdepositional transport we seek to develop methods for the estimation of seafloor statistics that take into account the finite precision resolution and sampling obtained by actual echo sounding systems in this initial paper we restrict our attention to the problem of recovering second order statistics from data sets collected by multibeam devices such as sea beam the seafloor is modeled as a stationary zero mean gaussian random field completely specified by its two point covariance function we introduce an anisotropic two point covariance function that has five free parameters describing the amplitude orientation characteristic wave numbers and hausdorff fractal dimension of seafloor topography we formulate the general forward problem relating this model to the statistics of an ideal multibeam echo sounder in particular the along track autocovariance functions of individual beams and the cross covariance functions between beams of arbitrary separation using these second moments as data functionals we then pose the inverse problem of estimating the seafloor parameters from realistic noisy data sets with finite sampling and beamwidth and we solve this inverse problem by an iterative linearized least squares method the inversion method is applied to sea beam transit data from both the pacific and atlantic oceans sea beam system noise stands out as a sharp spike on the along track autocovariance function and can be modeled as a white noise process whose amplitude generally increases with beam angle the five parameters in our second order model can be estimated from the inversion of data sets comprising 100 200 km of track length in general the cross track wave number is the most poorly determined although uncertainties in the assumed sea beam response may bias the values of the fractal dimension using the assumed beamwidth the measured noise values and the seafloor parameters recovered from the inversion we generate sea beam synthetics whose statistical character can be directly compared with raw sea beam data for most of the track segments we have processed thus far the synthetics are similar to the data in the case of one atlantic profile however the comparison clearly indicates the necessity of incorporating higher order statistics the space domain procedures described in this paper can be extended for this purpose

the slip history of the 1994 northridge california earthquake determined from strong motion teleseismic gps and leveling data we present a rupture model of the northridge earthquake determined from the joint inversion of near source strong ground motion recordings p and sh teleseismic body waves global positioning system gps displacement vectors and permanent uplift measured along leveling lines the fault is defined to strike 122 and dip 40 to the south southwest the average rake vector is determined to be 101 and average slip is 1 3 m the peak slip reaches about 3 m our estimate of the seismic moment is 1 3 0 2 10 26 dyne cm potency of 0 4 km3 the rupture area is small relative to the overall aftershock dimensions and is approximately 15 km along strike nearly 20 km in the dip direction and there is no indication of slip shallower than about 5 to 6 km the up dip strong motion velocity waveforms are dominated by large s wave pulses attributed to source directivity and are comprised of at least 2 to 3 distinct arrivals a few seconds apart stations at southern azimuths indicate two main s wave arrivals separated longer in time about 4 to 5 sec these observations are best modeled with a complex distribution of subevents the initial s wave arrival comes from an asperity that begins at the hypocenter and extends up dip and to the north where a second larger subevent is centered about 12 km away the secondary s arrivals at southern azimuths are best fit with additional energy radiation from another high slip region at a depth of 19 km 8 km west of the hypocenter the resolving power of the individual data sets is examined by predicting the geodetic gps and leveling displacements with the dislocation model determined from the waveform data and vice versa and also by analyzing how well the teleseismic solution predicts the recorded strong motions the general features of the geodetic displacements are not well predicted from the model determined independently from the strong motion data likewise the slip model determined from geodetic data does not adequately reproduce the strong motion characteristics whereas a particularly smooth slip pattern is sufficient to satisfy the geodetic data the strong motion and teleseismic data require a more heterogeneous slip distribution in order to reproduce the velocity amplitudes and frequency content although the teleseismic model can adequately reproduce the overall amplitude and frequency content of the strong motion velocity recordings it does a poor job of predicting the geodetic data consequently a robust representation of the slip history and heterogeneity requires a combined analysis of these data sets

mapping alluvial fans in death valley california using multichannel thermal infrared images we have mapped alluvial fans in death valley california using nasa s 8 12 m six channel airborne thermal infrared multispectral scanner tims we are able to recognize both composition and relative age differences age unit boundaries are generally consistent with those obtained by conventional mapping composition was verified by field investigation and comparison with existing geologic maps bedrock and its young derived fan gravels have similar emissivities the original composition of the fans is modified by differential erosion and weathering permitting relative age mapping with tims

slip distribution of the 1992 landers earthquake and its implications for earthquake source mechanics we use near source 10 to 164 km displacement seismograms to image the slip distribution and rupture history of the 28 june 1992 mw 7 3 landers california earthquake aftershock seismograms from similar distances are modeled to find the velocity model and frequency range 0 05 to 0 25 hz over which theoretical green s functions are most accurate and the measure of fit is used as an upper bound on theoretical error in the mainshock inversion we represent the rupture surface with three planar segments divided into 3 by 3 km elements extending from the surface to 18 km depth and solve for the slip distribution and the rupture time model that minimizes misfit to both recorded seismograms and mapped surface displacement in a least squares sense we investigate two faulting models one with a uniformly short 3 sec rise time everywhere on the fault surface and the other with a variable rise time 2 to 6 sec we perform sensitivity tests using synthetic data from a hypothetical earthquake to evaluate model resolution and solution stability in the sensitivity tests both fault models recover similar slip and rupture features but neither is capable of imaging details of the rise time for the landers earthquake we find an average rupture velocity of 2 5 km sec and use this average for a starting model in a linearized inversion for rupture time in our solutions of slip amplitude distribution the southernmost johnson valley fault segment has 20 of the total seismic moment 6 to 8 1019 n m with small displacements near the hypocenter the homestead valley segment contributes half of the moment with the largest slip amplitudes 25 to 35 km northwest of the hypocenter at 4 to 12 km depth and the camp rock emerson segment contributes the remaining moment with the largest slip amplitude 35 to 50 km northwest of the hypocenter in the shallow crust 9 km there is some evidence that the rupture front is delayed as it encounters high slip regions suggesting that prior to the mainshock these areas were further from failure owing either to greater strength or lower prestress

frequency domain inversion of strong motions application to the 1992 landers earthquake we present a frequency domain inversion in which the observed earthquake strong ground motions are used to constrain the space time dependence of slip on a fault green s functions are numerically evaluated and the parameters describing the rupture are the local slip rupture time and rise time these parameters are simultaneously evaluated without additional constraints this procedure allows for large variations in the local rupture velocity the june 28 1992 landers earthquake mw 7 3 offers an exceptional opportunity to apply this technique to a major strike slip event we model the rupture evolution including local differences in slip durations and variations in rupture velocity our final results are in good agreement with other inversion studies geodetic and surface observations the main discrepancies occurred at depth and at the end of the johnson valley fault we show that a relatively low resolution could be an explanation for these differences rupture velocity and slip are extremely heterogeneous both along strike and with depth a moment of 0 90 1020 n m was found the slip distribution shows that this event consists of a series of regions of high slip subevents separated by regions of relative low slip approximately 50 of the moment was released on the homestead valley fault in this region of large slip the rupture velocity inferred by our inversion is well constrained and is equal to 3 0 km s at depth and 2 5 km s near the surface our inversion favors the hypothesis that the duration of the slip at each point of the fault is of the order of the duration of rupture of each subevent

a faulting model for the 1979 imperial valley earthquake by comparing synthetic particle velocities with the near source strong motion data we have constructed by trial and error a faulting model for the 1979 imperial valley earthquake the calculation of the synthetic seismograms takes into account the vertical inhomogeneity of the elastic parameters in the imperial valley and the spatial variation of the slip rate parameters on the fault plane the independent slip rate parameters are 1 the strike slip rate amplitude 2 the dip slip rate amplitude 3 the duration that slip rate is nonzero the rise time of the slip function and 4 the rupture time which determines when the slip rate is initiated our faulting model has the following principal features 1 faulting occurred on the imperial fault and on the brawley fault rupture on the brawley fault being triggered by rupture on the imperial fault 2 the imperial fault is a plane 35 km long and 13 km wide with a strike of 323 measured clockwise from north and a dip of 80 ne the brawley fault is a 10 km long and 8 km wide plane with a strike of 360 and a dip of 90 3 faulting on the imperial fault is primarily right lateral strike slip with a small component of normal dip slip in the sediments at its northern end the larger strike slip rates are generally confined between depths of 5 and 13 km with maximum values of about 1 0 m s the duration varies on the fault with a maximum of 1 9 s which is considerably shorter than the total time for the rupture to take place 4 the rupture velocity on the imperial fault is highly variable locally it exceeds the shear wave velocity and in one instance the compressional wave velocity the average rupture velocity though is less than the shear wave velocity 5 although the slip on the brawley fault contributes only about 4 of the total moment it greatly affects the ground motion at nearby stations 6 the total seismic moment is 6 7 1018 n m where the imperial fault contributes 6 4 1018 n m and the brawley fault contributes 2 7 1017 n m in the process of trying almost 300 faulting models we found that given the elastic parameters of the medium the synthetic seismograms were most sensitive to the specification of the rupture time although the slip rate amplitudes are linearly related to the data the rupture time and the duration are not the parameterization of the nonlinear variables has a strong effect on the generation of synthetic seismograms from a finite fault

inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 imperial valley california earthquake a least squares point by point inversion of strong ground motion and teleseismic body waves is used to infer the fault rupture history of the 1979 imperial valley california earthquake the imperial fault is represented by a plane embedded in a half space where the elastic properties vary with depth the inversion yields both the spatial and temporal variations in dislocation on the fault plane for both right lateral strike slip and normal dip slip components of motion inversions are run for different fault dips and for both constant and variable rupture velocity models effects of different data sets are also investigated inversions are compared which use the strong ground motions alone the teleseismic body waves alone and simultaneously the strong ground motion and teleseismic records the inversions are stabilized by adding both smoothing and positivity constraints the moment is estimated to be 5 0 10 25 dyne cm and the fault dip 90 5 dislocation in the hypocentral region south of the united states mexican border is relatively small and almost dies out near the border dislocation then increases sharply north of the border to a maximum of about 2 m under interstate 8 dipslip motion is minor compared to strike slip motion and is concentrated in the sediments the best fitting constant rupture velocity is 80 per cent of the local shear wave velocity however there is a suggestion that the rupture front accelerated from the hypocenter northward the 1979 imperial valley earthquake can be characterized as a magnitude 5 earthquake at the hypocenter which then grew into or triggered a magnitude 6 earthquake north of the border

rupture history of the 1984 morgan hill california earthquake from the inversion of strong motion records near source strong motion velocity records and teleseismic short period p waveforms are modeled to obtain the spatial and temporal distribution of slip for the 1984 morgan hill earthquake both forward modeling and constrained least squares inversion techniques are used to interpret the strong motion velocity waveforms in the frequency range of approximately 0 2 to 2 0 hz these data support a nearly unilateral rupture to the southeast with a rupture propagation velocity of nine tenths of the local s wave velocity the majority of the slip occurs over a fault length of 25 km and to a first approximation can be interpreted as two main source regions each with an extent of about 5 km with their centers separated by about 12 km however each of the sources has detailed structure of its own and a simple two point source model is not an accurate representation of the morgan hill earthquake the second source occurs about 4 5 sec after the first and is approximately 3 times larger the maximum dislocation on the fault plane is about 1 m the total moment of the earthquake is estimated to be 2 1 10 25 dyne cm the morgan hill earthquake offers convincing evidence for very inhomogeneous slip and stress distributions on shallow strike slip faults

a stochastic view of lower crustal fabric based on evidence from the ivrea zone despite its complicated history the ivrea zone is considered to be a representative surface exposure of extended continental crust we have digitized two standard 1 25 000 geological maps from this area and evaluated their structural statistics because of the subvertical orientation of the ivrea zone these maps can be considered as small scale cross sections through the lower continental crust the autocorrelation functions of the digitized lithologies measured from these maps show a clear self similar or fractal rather than a gaussian or deterministic trend we found that an anisotropic von karman correlation function with an aspect ratio around 4 and a hurst number of 0 3 corresponding to a fractal dimension of 2 7 matches the observed data our results represent an explicit confirmation of previous indirect evidence for the fractal nature of lithospheric heterogeneities and provide the means to construct realistic crustal scale seismic models of ivrea type lower continental crust

fractal analysis applied to characteristic segments of the san andreas fault fractal theory is applied in a quantitative analysis of the san andreas fault saf geometry the method which directly measures the increase in total fault length with a decrease in ruler size gives the fractal dimension d and scaling properties for the chosen length band 0 5 1000 km a physical interpretation of d is that it measures the irregularity of the fault trace in the selected band the fault is subdivided into six segments of distinctive seismic behavior a main fault trace which shows either maximum coseismic slip or creeping was selected for analysis with three alternative branches examined for the saf system south of san bernardino branches of the fault trace are not considered in this work fractal dimensions calculated for the different segments range from 1 0008 to 1 0191 and are different from 1 0 to the 95 confidence interval these small changes reflect the overall smoothness of the main fault trace significant variations in d among segments indicate heterogeneities in the fault smoothness along strike d also changes significantly from the short length band to the long length band where the demarcation point ranges from 1 to 2 km the short length band has larger d values a slight correlation is obtained between the fractal dimension of the main trace and the extent of subparallel faulting this indicates some correspondence between the main fault trace irregularity and the complexity of subsidiary fault traces in plan view

rupture process of the 1999 kocaeli turkey earthquake estimated from strong motion waveforms the multiple time window linear waveform inversion is applied to study the source process of the 1999 kocaeli turkey earthquake using strong motion data two peculiar observations characterize this event the first one is a conflicting observation of the long surface ruptures and short source duration resulting from the observed strong motion waveforms the second is an anomalously short s p time observed at a station skr considering its hypocentral distance a supershear rupture propagation ellsworth and celebi 1999 and the p wave triggering of an asperity anderson et al 2000 were both proposed to explain the record at skr our waveform inversion study is aimed to examine these possibilities in order to clarify the relation between the spatial distribution of the surface ruptures and coseismic slips we divide the assumed fault plane into three parts and look for the first time window front propagation velocity propagation velocity of triggering front of the first time window at each subfault in each part by performing numerous inversions changing the first time window front propagation velocity separately in each part and selecting the best first time window propagation velocities that give largest variance reduction at the same time appropriate smoothness of spatiotemporal slip distribution is searched for by performing inversions with different strength of smoothing constraints against the observational equations and judging the relative appropriateness of inversion results in terms of the akaike s bayesian information criterion the best source model is characterized by an asperity that is about 35 km east of the epicenter and near a segment boundary at sapanca lake triggered by the first time window front with 5 8 km sec close to p wave velocity the time progression of the rupture suggests a scenario where a p wave arriving from the hypocenter or somewhere around the hypocenter triggered the second largest asperity near the junction of the segments at sapanca lake which is a preferable interpretation the first time window propagation velocities west of the hypocenter and east of skr are expected to be 3 0 km sec manuscript received 20 august 2000