surface wave tomography of the western united states from ambient seismic noise rayleigh and love wave phase velocity maps summary we present the results of rayleigh wave and love wave phase velocity tomography in the western united states using ambient seismic noise observed at over 250 broad band stations from the earthscope usarray transportable array and regional networks all available threecomponent time series for the 12 month span between 2005 november 1 and 2006 october 31 have been cross correlated to yield estimated empirical rayleigh and love wave green s functions the love wave signals were observed with higher average signal to noise ratio snr than rayleigh wave signals and hence cannot be fully explained by the scattering of rayleigh waves phase velocity dispersion curves for both rayleigh and love waves between 5 and 40 speriod were measured for each interstation path by applying frequency time analysis the average uncertainty and systematic bias of the measurements are estimated using a method based on analysing thousands of nearly linearly aligned station triplets we find that empirical green s functions can be estimated accurately from the negative time derivative of the symmetric component ambient noise cross correlation without explicit knowledge of the source distribution the average traveltime uncertainty is less than 1 s at periods shorter than 24 s we present rayleigh and love wave phase speed maps at periods of 8 12 16 and 20 s the maps show clear correlations with major geological structures and qualitative agreement with previous results based on rayleigh wave group speeds

internal structure of the san jacinto fault zone in the trifurcation area southeast of anza california from data of dense seismic arrays we image the internal structure of the san jacinto fault zone sjfz in the trifurcation area southeast of anza california with seismic records from dense linear and rectangular arrays the examined data include recordings from more than 20 000 local earthquakes and nine teleseismic events automatic detection algorithms and visual inspection are used to identify p and s body waves along with p and s types fault zone trapped waves fztw the location at depth of the main branch of the sjfz the clark fault is identified from systematic waveform changes across lines of sensors within the dense rectangular array delay times of p arrivals from teleseismic and local events indicate damage asymmetry across the fault with higher damage to the ne producing a local reversal of the velocity contrast in the shallow crust with respect to the large scale structure a portion of the damage zone between the main fault and a second mapped surface trace to the ne generates p and s types fztw inversions of high quality s type fztw indicate that the most likely parameters of the trapping structure are width of 70 m s wave velocity reduction of 60 per cent q value of 60 and depth of 2 km the local reversal of the shallow velocity contrast across the fault with respect to large scale structure is consistent with preferred propagation of earthquake ruptures in the area to the nw

characterization of the san jacinto fault zone near anza california by fault zone trapped waves we installed three 350 m long seismic arrays each array consisting of 12 three component stations across the coyote creek fault ccf clark valley fault cvf and buck ridge fault brf of the san jacinto fault zone sjfz near anza california to record fault zone trapped waves from microearthquakes we observed trapped waves with relatively large amplitudes and long duration at stations close to the fault traces for earthquakes occurring within the fault zone the coda normalized amplitude spectra of trapped waves showed peaks at 4 7 hz which decreased sharply with the distance from the fault trace observations and three dimensional finite difference simulations of trapped waves revealed low velocity and low q waveguides on these active faults with the width of 75 100 m in which shear velocities are reduced by 25 30 from wall rock velocities and q values are 40 90 at depths between the surface and 18 km the locations of earthquakes for which we observed trapped waves delineate the most seismically active fault strands of the sjfz in a region with complicated slip planes near anza the low velocity waveguides inferred from trapped waves extend 15 to 20 km in the length on these active faults and are segmented by the fault discontinuities the waveguide on the brf dips south westward to connect the waveguide on the cvf which dips northeastward this waveguide extends at the seismogenic depth through anza slip gap to another low velocity waveguide on the casa loma fault clf which has been delineated in our previous study of the sjfz using trapped waves li et al 1997 the waveguide on the ccf in coyote mountain is nearly vertical and disconnected from the clf at the south edge of anza gap we interpret the low velocity waveguides on these active strands to partly result from recent prehistoric significant earthquakes on them and evaluate the future earthquake in the anza region

ground motion prediction equations in the san jacinto fault zone significant effects of rupture directivity and fault zone amplification we present a new set of ground motion prediction equations gmpes for horizontal peak ground acceleration peak ground velocity and 5 damped pseudo spectral acceleration psa developed for the san jacinto fault zone sjfz area besides using these equations to quantify seismic shaking in the area the results allow us to examine the physics and local properties controlling the observed ground motions the analyzed dataset includes 30 000 observations from 800 events spanning a magnitude range of 1 5 m 6 0 and recorded by up to 140 stations at epicentral distances ranging from essentially zero to 150 km the local gmpe is developed for the sjfz by applying classical regression techniques with predictive variables that include first distance and magnitude and then site characteristics rupture directivity and fault zone amplification the significance of these effects is determined by measuring the uncertainty reduction of the gmpe due to each factor the results show that in contrast to many regional studies traditional site characteristic has a relatively minor effect on peak amplitudes in our study area however rupture directivity is a significant factor controlling the amplitudes of ground motion even for small events the dense seismic network and newly developed directivity tool enable us to extract efficiently directivity effects with statistical significance using the ground motion dataset during the regression analysis process the obtained rupture directivities are consistent with the main focal mechanism orientations and surface trace orientations known from other studies and predictions for bimaterial ruptures in the trifurcation area of the sjfz fault zone amplification is a second important factor showing strong impact on the peak ground motion values with increasing role for the lower frequency range 10 hz examined in the 5 damped psa values we also observe signatures of large amplitude variances which indicate additional source related control on the distribution of amplitudes besides rupture directivity for aftershocks close in time and location to the m l 5 1 earthquake of march 2013 using the full set of records we present the most complete set of gmpes for the sjfz area including a higher amplitude prediction for regions in the direction of rupture

shallow slip deficit due to large strike slip earthquakes in dynamic rupture simulations with elasto plastic off fault response summary slip inversions of geodetic data from several large magnitude 7 strike slip earthquakes point to coseismic slip deficit at shallow depths 3 4 km that is coseismic slip appears to decrease towards the earth surface while the inferred slip distribution may be consistent with laboratory derived rate and state friction laws suggesting that the uppermost brittle crust may be velocity strengthening there remains a question of how the coseismic slip deficit is accommodated throughout the earthquake cycle the consequence of velocity strengthening fault friction at shallow depths is that the deficit of coseismic slip is relieved by post seismic afterslip and interseismic creep however many seismic events with inferred shallow slip deficit were not associated with either resolvable shallow interseismic creep or robust shallow afterslip hence the origin of shallow slip deficit remains uncertain in this study we investigate whether inelastic failure in the shallow crust due to dynamic earthquake rupture can explain the inferred deficit of shallow slip evidence for such failure is emerging from geologic seismic and geodetic observations we find that the amount of shallow slip deficit is proportional to the amount of inelastic deformation near the earth surface such deformation occurs under a wide range of parameters that characterize rock strength in the upper crust however the largest magnitude of slip deficit in models accounting for off fault yielding is 2 4 times smaller than that inferred from kinematic inversions of geodetic data to explain this discrepancy we further explore to what extent assumptions in the kinematic inversions may bias the inferred slip distributions inelastic deformation in the shallow crust reduces coseismic strain near the fault introducing an additional artificial deficit of up to 10 per cent of the maximum slip in inversions of geodetic data that are based on purely elastic models the largest magnitude of slip deficit in our models combined with the bias in inversions accounts for up to 25 per cent of shallow slip deficit which is comparable but still smaller than 30 60 per cent deficit inferred from kinematic inversions we discuss potential mechanisms that may account for the remaining discrepancy between slip deficit predicted by elasto plastic rupture models and that inferred from inversions of space geodetic data

geological observations of damage asymmetry in the structure of the san jacinto san andreas and punchbowl faults in southern california a possible indicator for preferred rupture propagation direction we present new in situ observations of systematic asymmetry in the pattern of damage expressed by fault zone rocks along sections of the san andreas san jacinto and punchbowl faults in southern california the observed structural asymmetry has consistent manifestations at a fault core scale of millimeters to meters a fault zone scale of meters to tens of meters and related geomorphologic features the observed asymmetric signals are in agreement with other geological and geophysical observations of structural asymmetry in a damage zone scale of tens to hundreds of meters in all of those scales more damage is found on the side of the fault with faster seismic velocities at seismogenic depths the observed correlation between the damage asymmetry and local seismic velocity structure is compatible with theoretical predictions associated with preferred propagation direction of earthquake ruptures along faults that separate different crustal blocks the data are consistent with a preferred northwestward propagation direction for ruptures on all three faults if our results are supported by additional observations asymmetry of structural properties determined in field studies can be utilized to infer preferred propagation direction of large earthquake ruptures along a given fault section the property of a preferred rupture direction can explain anomalous behavior of historic rupture events and may have profound implications for many aspects of earthquake physics on large faults

rupture directivity of small earthquakes at parkfield 1 theoretical modeling of strike slip ruptures along a bimaterial interface suggests that earthquakes initiating on the interface will have a preferred rupture direction we test this model with 450 small earthquakes 2 3 exhibiting southeast rupture characteristics some spatial and temporal variability in rupture directivity is also apparent we observe a higher proportion of northwest directivity ruptures following the 2004 m 6 parkfield earthquake which ruptured toward the northwest our results are generally consistent with the preferred southeast rupture directivity model but suggest that directivity is likely due to several contributing factors

least squares horizons with local slopes and multigrid correlations most seismic horizon extraction methods are based on seismic local reflection slopes that locally follow seismic structural features however these methods often fail to correctly track horizons across discontinuities such faults and noise because the local slopes can only correctly follow laterally continuous reflections in addition seismic amplitude or phase information is not used in these methods to compute horizons that follow consistent phase e g peaks or troughs to solve these problems we propose a novel method to compute horizons that globally fit both local slopes and multi grid correlations of seismic traces in this method we first estimate local reflection slopes by using structure tensors and compute laterally multi grid slopes by using dynamic time warping to correlate seismic traces within multiple laterally coarse grids these coarse grid slopes can correctly correlate reflections that may be significantly dislocated by faults or other discontinuous structures we then compute a

a 21 event 4 000 year history of surface ruptures in the anza seismic gap san jacinto fault and implications for long term earthquake production on a major plate boundary fault paleoseismic work completed at hog lake on the san jacinto fault sjf near anza california indicates that at least 21 surface ruptures have occurred in the anza seismic gap over the past 4 000 years the ages of the ruptures are constrained by 111 radiocarbon dates 97 of which fall in stratigraphic order the average recurrence interval for all ruptures for this period is about 185 105 years although some ruptures such as occurred in the april 1918 earthquake caused only minor displacement we rate the expression of each interpreted event in each of the twelve developed field exposures presented in this work by assigning numeric values for the presence of different criteria that indicate rupture to a paleo ground surface weakly expressed ruptures for example the deformation we interpret to be the result of the historical 1918 earthquake received low scores and are interpreted as smaller earthquakes from this analysis we infer that at least fifteen of the identified ruptures are indicative of large earthquakes similar to the penultimate earthquake inferred to be the m w 7 3 22 november 1800 earthquake the adjusted recurrence interval for large earthquakes lengthens to approximately 254 years comparison with the rupture history at the mystic lake paleoseismic site on the claremont strand indicates that it is plausible that several of the large ruptures identified at hog lake could have jumped the hemet step over at mystic lake and continued on the claremont strand or vice versa but most of the event ages do not match between the two sites indicating that most ruptures do not jump the step finally comparison with san andreas fault ruptures both to the north and south of its juncture with the sjf suggest that some northern sjf ruptures identified at mystic lake may correlate with events identified at wrightwood but that these northern ruptures have no match at hog lake and can not indicate rupture of the entire sjf onto the saf

catastrophic extinction of caribbean rudist bivalves at the cretaceous tertiary boundary strontium isotope ratios 87 sr 86 sr in pristine low mg calcite of shells of rudist bivalves from the titanosarcolites limestones exposed in the central maldon and marchmont inliers of jamaica indicate that species rich rudist coral associations persisted into the latest maastrichtian 66 65 ma this finding contradicts the currently accepted hypothesis of stepwise extinction of rudist bivalves in the middle maastrichtian and argues for a catastrophic impact related demise of caribbean cretaceous reefal ecosystems at the cretaceous tertiary boundary

focal spot imaging based on zero lag cross correlation amplitude fields application to dense array data at the san jacinto fault zone we image the subsurface below a dense seismic array straddling the clark branch of the san jacinto fault zone in southern california the analysis is based on focal spots of surface waves associated with the zero lag amplitudes of noise cross correlations computed between all stations of the dense array local medium properties are inferred from the spatially variable focal spot size and shape based on the first zero crossing of amplitude versus distance distributions the method provides simultaneous estimates of wave speed apparent attenuation and anisotropy without solving a tomographic inverse problem the obtained images of the frequency dependent seismic velocity distributions are consistent with independent estimates from a far field rayleigh wave tomography we observe an anticorrelation between our apparent attenuation coefficient and seismic velocity and a fault parallel alignment of fast propagation directions with greater structural complexity to the southwest of the fault the results imply a complex fault zone structure including a waveguide to the northeast of the fault that is continuous across the observed depth range and a low velocity structure to the southwest associated with a shallow sedimentary basin

pleistocene brawley and ocotillo formations evidence for initial strike slip deformation along the san felipe and san jacinto fault zones southern california abstract we examine the pleistocene tectonic reorganization of the pacific north american plate boundary in the salton trough of southern california with an integrated approach that includes basin analysis magnetostratigraphy and geologic mapping of upper pliocene to pleistocene sedimentary rocks in the san felipe hills these deposits preserve the earliest sedimentary record of movement on the san felipe and san jacinto fault zones that replaced and deactivated the late cenozoic west salton detachment fault sandstone and mudstone of the brawley formation accumulated between 1 1 and 0 6 0 5 ma in a delta on the margin of an arid pleistocene lake which received sediment from alluvial fans of the ocotillo formation to the west southwest our analysis indicates that the ocotillo and brawley formations prograded abruptly to the east northeast across a former mud dominated perennial lake borrego formation at 1 1 ma in response to initiation of the dextral oblique san felipe fault zone the 25 km long

seismic and geodetic evidence for extensive long lived fault damage zones during earthquakes slip is often localized on preexisting faults but it is not well understood how the structure of crustal faults may contribute to slip localization and energetics growing evidence suggests that the crust along active faults undergoes anomalous strain and damage during large earthquakes seismic and geodetic data from the calico fault in the eastern california shear zone reveal a wide zone of reduced seismic velocities and effective elastic moduli using seismic traveltimes trapped waves and interferometric synthetic aperture radar observations we document seismic velocities reduced by 40 50 and shear moduli reduced by 65 compared to wall rock in a 1 5 km wide zone along the calico fault observed velocity reductions likely represent the cumulative mechanical damage from past earthquake ruptures no large earthquake has broken the calico fault in historic time implying that fault damage persists for hundreds or perhaps thousands of years these fi ndings indicate that faults can affect rock properties at substantial distances from primary fault slip surfaces and throughout much of the seismogenic zone a result with implications for the amount of energy expended during rupture to drive cracking and yielding of rock and development of fault systems

stratigraphic coordinates a coordinate system tailored to seismic interpretation in certain seismic data processing and interpretation tasks such as spiking deconvolution tuning analysis impedance inversion and spectral decomposition it is commonly assumed that the vertical direction is normal to reflectors this assumption is false in the case of dipping layers and may therefore lead to inaccurate results to overcome this limitation we propose a coordinate system in which geometry follows the shape of each reflector and the vertical direction corresponds to normal reflectivity we call this coordinate system stratigraphic coordinates we develop a constructive algorithm that transfers seismic images into the stratigraphic coordinate system the algorithm consists of two steps first local slopes of seismic events are estimated by plane wave destruction then structural information is spread along the estimated local slopes and horizons are picked everywhere in the seismic volume by the predictive painting algorithm these picked horizons represent level sets of the first axis of the stratigraphic coordinate system next an upwind finite difference scheme is used to find the two other axes which are perpendicular to the first axis by solving the appropriate gradient equations after seismic data are transformed into stratigraphic coordinates seismic horizons should appear flat and seismic traces should represent the direction normal to the reflectors immediate applications of the stratigraphic coordinate system are in seismic image flattening and spectral decomposition synthetic and real data examples demonstrate the effectiveness of stratigraphic coordinates

properties of seismic fault zone waves and their utility for imaging low velocity structures a two dimensional solution for the scalar wave equation in a model of two vertical layers between two quarter spaces is used to study properties of seismic waves in a laterally heterogeneous low velocity structure the waves referred to as seismic fault zone waves include head waves internal fault zone reflections and trapped waves the analysis aims to clarify the dependency of the waves on media velocities media attenuation coefficients layer widths and source receiver geometry additional calculations with extreme low velocity layers provide examples that may be relevant for volcanic and geothermal domains the interference patterns controlling seismic fault zone waves change with the number of internal reflections in the low velocity structure this number increases with propagation distance along the structure decreases with fault zone width and increases for given length scales with the velocity contrast the relative lateral position of the source within the low velocity layer modifies die length scales associated with internal reflections and influences the resulting interference pattern low values of q affect considerably the dominant period and overall duration of the waves thus there are significant tradeoffs between propagation distance along the structure fault zone width velocity contrast source location within the fault zone and q the lateral and depth receiver coordinates determine the particular point where the interference pattern is sampled and observed motion is a strong function of both coordinates the zone connecting sources generating fault zone waves and observation points with appreciable wave amplitude can be over an order of magnitude larger than the fault zone width calculations for cases with layer p wave velocity of 200 m s 1 modeling a vertical dike or crack with fluid and gas show conspicuous persisting oscillations the results resemble aspects of seismic data in volcanic domains if these waves exist in observed records their explicit recognition and modeling will help to separate source and structural effects and aid in the interpretation of volcano seismology signals although the tradeoffs in parameters governing seismic fault zone waves are significant each variable has its own signature and the parameters may be constrained by additional geophysical data simultaneous modeling of many waveforms with an appropriate solution can resolve the various parameters and provide a high resolution structural image

dynamic rupture on a material interface with spontaneous generation of plastic strain in the bulk abstract we discuss 2d simulations of in plane ruptures on a fault between different solids in a model that includes dynamic generation of plastic strain off the fault the fault is governed by a regularized friction with a gradual response to abrupt changes of normal stress while plastic strain in the bulk is governed by a coulomb yielding criterion the occurrence of off fault damage stabilizes the slip velocity on the fault stable self sustaining ruptures propagate on the material interface as narrow unidirectional pulses in the direction of slip on the more compliant side of the fault plastic strain is generated only on the stiffer side of the fault the tensional quadrant in the direction of rupture propagation in a strip of approximately constant width that is correlated with the slip velocity on the fault simulations for various values of confining normal stress material contrast and cohesion in the bulk suggest that significant generation of off fault damage is limited to conditions that correspond to the top few kilometers of the crust the results are compatible with field observations of shallow asymmetric damage patterns in the structure of several large faults

horizon volumes with interpreted constraints abstracthorizons are geologically significant surfaces that can be extracted from seismic images color coding of horizons based on amplitude or other attributes can help reveal ancient sedimentary environments and structural features extracted horizons are also used for building structure models and stratigraphic interpretations we propose two methods for constructing seismic horizons aligned with reflectors in a 3d seismic image the first method generates horizons one at a time the second method generates an entire volume of horizons at once by first computing a relative geologic time volume from seismic normal vectors rather than gradually building a horizon by extending one or more seed points to a surface along seismic reflectors both of our methods generate whole horizons at once by solving partial differential equations derived from seismic normal vectors the most significant new aspect of both methods is the ability to specify perhaps interactively during interpretation a small number of

structural properties of the san jacinto fault zone at blackburn saddle from seismic data of a dense linear array we image the san jacinto fault zone at blackburn saddle using earthquake waveforms recorded by a 2 km across fault linear array with 108 three component sensors separated by 10 30 m the length and spatiotemporal sampling of the array allow us to derive high resolution information on the internal fault zone structure with spatial extent that can be merged with regional tomography models across fault variations in polarization amplitude and arrival time of teleseismic p waves indicate abrupt changes in subsurface structure near the surface trace of the fault sensor bs55 and 270 m to the northeast sensor bs34 analysis of fault zone head waves from local events reveals the existence of a deep bimaterial interface that extends from the array to at least 50 km southeast and has a section with 10 velocity contrast this analysis also corroborates the teleseismic results and indicates a broad damage zone primarily northeast of the fault bounded by a shallow bimaterial interface near bs34 that merges with the deep interface detection and waveform inversions of love type fault zone trapped waves generated by local events indicate a trapping structure within the broader damage zone with width of 150 m velocity reduction of 55 from the surrounding rock and depth extent of 2 km the performed analyses provide consistent results on the subsurface location of the main seismogenic fault and properties of a major bimaterial interface and damage structure the imaged fault zone properties are consistent with preferred propagation direction of earthquake ruptures in the area to the northwest

late quaternary rate of slip along the san jacinto fault zone near anza southern california the anza section of the san jacinto fault southern california has not experienced a groundbreaking earthquake since at least 1918 and probably since before 1899 leading previous workers to designate this 20 km long fault segment as a seismic gap study of displaced pleistocene and holocene alluvial fan and fluvial deposits across the san jacinto fault near anza california dated using 14c and soil development age control indicates a minimum right lateral slip rate of 9 2 2 mm yr since 9 5 ka and a slip rate of 11 9 5 mm yr since 14 ka 12 9 5 mm yr since 17 ka and 13 10 6 mm yr since 50 ka these estimates agree with a previously determined 700 ka to present minimum slip rate of 10 2 mm yr based on an estimated average slip rate of 12 mm yr from the longer term slip rate estimates about 0 8 m of potential slip has accumulated if the 1918 earthquake broke the anza segment over 1 1 m of potential slip has accumulated if the last earthquake occurred prior to 1899 these data attest to the potential for an earthquake m6 5 for this section of the fault

seismic radiation from an sh line source in a laterally heterogeneous planar fault zone a solution to the class of problems in which an sh line source is operating in a medium consisting of two quarter spaces with an arbitrary number of semi infinite vertical layers in between is given in a frequency wavenumber domain the general result is applied to a model with a single vertical layer and synthetic seismograms are calculated for a source specified as an sh line dislocation at the interface between the vertical layer and the faster quarter space the numerical examples demonstrate that moderate fault normal material heterogeneity can have significant effects on the seismic response of a medium the overall distribution of travel times amplitudes and motion polarities in the presence of fault zone material heterogeneity is different from a corresponding distribution in a homogeneous half space the effects are most prominent near the fault zone where head waves surface waves and trapped modes result in waveform complexities large amplifications over an order of magnitude in the fault zone itself and first motion polarities which are reversed from homogeneous medium predictions the results imply that material heterogeneity should be included in fault zone models which are used for the interpretation of observed seismic data and theoretical calculations of near fault responses such models can be utilized to obtain important fault zone and earthquake source parameters

shallow seismic trapping structure in the san jacinto fault zone near anza california summary we analyse fault zone trapped waves generated by 500 small earthquakes for high resolution imaging of the subsurface structure of the coyote creek clark valley and buck ridge branches of the san jacinto fault zone near anza california based on a small number of selected trapped waves within this data set a previous study concluded on the existence of a lowvelocity waveguide that is continuous to a depth of 15 20 km in contrast our systematic analysis of the larger data set indicates a shallow trapping structure that extends only to a depth of 3 5 km this is based on the following lines of evidence 1 earthquakes clearly outside these fault branches generate fault zone trapped waves that are recorded by stations within the fault zones 2 a traveltime analysis of the difference between the direct s arrivals and trapped wave groups shows no systematic increase moveout with increasing hypocentral distance or event depth estimates based on the observed average moveout values indicate that the propagation distances within the low velocity fault zone layers are 3 5 km 3 quantitative waveform inversions of trapped wave data indicate similar short propagation distances within the low velocity fault zone layers the results are compatible with recent inferences on shallow trapping structures along several other faults and rupture zones the waveform inversions also indicate that the shallow trapping structures are offset to the northeast from the surface trace of each fault branch this may result from a preferred propagation direction of large earthquake ruptures on the san jacinto fault

global campanian upper cretaceous ostracod palaeobiogeography abstract the distribution of 218 campanian marine ostracod genera from 38 geographical areas was used for a quantitative palaeobiogeographical analysis the data were analysed using multivariate statistical methods ordination methods and parsimony analysis of endemicity these methods resulted in distinct geographical area groupings of the studied regions with these groupings seven palaeobiogeographical units based on campanian ostracod occurrences are distinguished the composition of the area groupings and the distribution of the palaeobiogeographical units are discussed their development and spatial distribution are broadly controlled by campanian palaeoclimates partly modified by palaeogeographical conditions

surface wave array tomography in se tibet from ambient seismic noise and two station analysis i phase velocity maps summary empirical green s functions egfs between pairs of seismographs can be estimated from the time derivative of the long time cross correlation of ambient seismic noise these egfs reveal velocity dispersion at relatively short periods which can be used to resolve structures in the crust and uppermost mantle better than with traditional surface wave tomography we combine rayleigh wave dispersion estimates from egfs and from traditional two station ts analysis into a new approach to surface wave array tomography with data from dense receiver arrays we illustrate the methodology with continuous broad band recordings from a temporary seismographic network on the southeastern part of the tibetan plateau in sichuan and yunnan provinces sw china the egfs are robust under temporal changes in regional seismicity and the use of either ambient noise approximated by records without signal from events with magnitude mb 5 or 4 or surface wave coda produces similar results the egfs do not strongly depend on the presence of large earthquakes but they are not reciprocal for stations aligned in the n s direction this directionality reflects the paucity of seismicity to the north of the array using a far field representation of the surface wave green s function and an image transformation technique we infer from the egfs the rayleigh wave phase velocity dispersion in the period band from 10 30 s a classical ts approach is used to determine rayleigh wave phase velocity dispersion between 20 120 s together they constrain phase velocity variations for t 10 120 s which can be used to study the structure from the crust to the upper mantle beneath se tibet short and intermediate period 10 80 s phase velocities are prominently low suggesting that the crust and upper mantle beneath se tibet is characterized by slow shear wave propagation

quantitative analysis of seismic fault zone waves in the rupture zone of the 1992 landers california earthquake evidence for a shallow trapping structure summary we analyse quantitatively a waveform data set of 238 earthquakes recorded by a dense seismic array across and along the rupture zone of the 1992 landers earthquake a grid search method with station delay corrections is used to locate events that do not have catalogue locations the quality of fault zone trapped waves generated by each event is determined from the ratios of seismic energy in time windows corresponding to trapped waves and direct s waves at stations close to and off the fault zone approximately 70 per cent of the events with s p times of less than 2 s including many clearly off the fault produce considerable trapped wave energy this distribution is in marked contrast with previous claims that trapped waves are generated only by sources close to or inside the landers rupture zone the time difference between the s arrival and trapped waves group does not grow systematically with increasing hypocentral distance and depth the dispersion measured from the trapped waves is weak these results imply that the seismic trapping structure at the landers rupture zone is shallow and does not extend continuously along strike by more than a few kilometres synthetic waveform modelling indicates that the fault zone waveguide has depth of approximately 2 4 km a width of approximately 200 m an s wave velocity reduction relative to the host rock of approximately 30 40 per cent and an s wave attenuation coefficient of approximately 20 30 the fault zone waveguide north of the array appears to be shallower and weaker than that south of the array the waveform modelling also indicates that the seismic trapping structure below the array is centred approximately 100 m east of the surface break

body and surface wave reconstruction from seismic noise correlations between arrays at piton de la fournaise volcano body wave reconstruction from ambient seismic noise correlations is an important step toward improving volcano imaging and monitoring here we extract body and surface waves that propagate in piton de la fournaise volcano on la reunion island using ambient noise cross correlation and array processing techniques ambient noise was continuously recorded at three dense arrays each comprising 49 geophones to identify and enhance the green s function from the ambient noise correlation we apply a double beamforming dbf technique between the array pairs the dbf allows us to separate surface and body waves direct and reflected waves and multipathing waves based on their azimuths and slownesses we successfully extract body waves between all the combinations of arrays including the wave that propagates through the active magmatic system of the volcano additionally we identify the effects of uneven noise source distribution and interpret the surface wave reflections

pulverized fault rocks and damage asymmetry along the arima takatsuki tectonic line japan abstract we present field and laboratory data on pulverized rocks at the hakusui kyo outcrop of the arima takatsuki tectonic line attl which is a dextral strike slip fault with 17 km displacement juxtaposing granite to the south against rhyolite to the north the majority of slip at the surface is localized to a clay rich gouge fault core 8 10 cm in width surrounded by a coarsening outwards fault breccias up to 3 m wide fault damage is highly asymmetric with respect to the slipping zone the granite south of the fault has a pulverized damage zone up to 200 m wide while the rhyolite to the north has only about 3 m wide non pulverized fault breccia the degree of pulverization in the granite decreases approximately logarithmically with normal distance from the slip zone the highly fractured pulverized rocks exhibit several distinct textural characteristics in thin section grains appear to be highly comminuted but the original grain shapes and margins are recognizable microfractures tend to be tensile in no preferred orientation grain fragments display little to no rotation and lack evidence of in situ shear consequently at macroscale the rocks appear to preserve original granitic textures despite being highly fractured and friable the observed pulverization and rock damage asymmetry are most consistent with generation mechanism involving ruptures on a bimaterial interface with statistically preferred propagation direction leading to damage primarily on the side with higher seismic velocity at depth this is supported by laboratory measurements of p wave ultrasonic velocities on intact samples which indicate that the granites have consistently higher velocity than the rhyolite with increasing confining pressure

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

application of high resolution dem data to detect rock damage from geomorphic signals along the central san jacinto fault abstract we analyze geomorphic properties extracted from lidar and srtm shuttle radar topography mission data to test whether the damage zone along the central san jacinto fault sjf zone can be resolved with remotely sensed data in a quantitative fashion the sjf is one of the most active faults in southern california with well expressed geomorphology and a fast slip rate as seen in the geology and by gps we use arcmap and the taudem toolbox to compare several morphometric parameters including drainage density dd on both sides of the fault using a 1 km and a 5 km buffer for the lidar and srtm data respectively we also analyze the spatial patterns of dd near the fault using two different definitions of spatial dd the high resolution of the lidar data allows us to focus on a single fault eliminating the effects of parallel nearby faults from the lidar data we find that the highest dd values occur in areas between two fault strands followed generally by rocks on the northeast side of the fault with the lowest dd values occurring on the southwest side of the fault the srtm data shows a band of high dd values centered on the main fault trace with 1 km width our results indicate that there is a strong correlation between drainage density and proximity to the fault with zones of structural complexity along the fault displaying the highest dd we interpret this to largely be an effect of degree of rock damage as these are areas that are expected to be more damaged and field observations support this contention if we are correct then it appears that the northeast side of the sjf is generally more damaged south of the trifurcation area there is evidence that the signal is reversed on the larger scale with more damage on the southwest side of the fault inferred from the srtm data possibly caused by extension between the coyote creek and clark faults the implications of the observed asymmetry could be geological evidence for rupture propagation direction because a preferred propagation direction is predicted to produce asymmetric damage structure that would be recorded in the volume of rock surrounding a fault

spatial variations in slip deficit on the central san andreas fault from insar summary we use ers insar measurements to record spatial variations in creep rate along the creeping segment of the san andreas fault saf california between 1992 and 2001 inversion of geodetic data yields a slip rate distribution along the creeping segment which is used for firstorder moment release and deficit calculations we present a time averaged spatial picture of surface deformation and associated subsurface creep an interferometric stack is constructed from 12 interferograms that show good coherence for the decade of observation the total right lateral offset spanned by the data is 34 mm yr 1 along most of the length of the creeping segment this offset occurs within a narrow 2 km zone close to the fault trace in the northern part a minor part of the offset is taken up by the nearby calaveras paicines fault in general the observed rates of surface creep are consistent with those obtained by several other studies for a longer and or earlier period of time using different geodetic methods this suggests that the average creep rate has been constant over a period of almost four decades a joint gps insar inversion implies that the shallow creep rate is variable along strike reaching up to 31 5 1m m yr 1 in the central section of the creeping segment tapering off along strike to the south and becoming partitioned across two subparallel faults in the north the deep slip rate beneath the seismogenic layer is 33 3m m yr 1 the difference between shallow and deep slip rates suggests that there is a shallow slip deficit on the creeping segment of the saf csaf moment release rate due to aseismic slip is approximately three orders of magnitude greater than seismic moment release the annual creep on the csaf is equivalent to the moment released in a m 6 earthquake the equivalent moment of the slip deficit relative to the deep slip rate is between 4 1 10 17 and 8 4 10 17 nm yr 1 which is equivalent to a magnitude 5 7 5 9 earthquake over a 150 yr period the deficit is equivalent to a m 7 2 7 4 earthquake this slip deficit may be compensated by intermittent moderate sized events or by accelerated afterslip following large events on nearby fault segments

generating a relative geologic time volume by 3d graph cut phase unwrapping method with horizon and unconformity constraints abstractconstruction of a relative geologic time rgt volume is vital to seismic geomorphological and sedimentological interpretation seismic instantaneous phase unwrapping provides an excellent approach for generating an rgt volume although several 2d or 3d seismic phase unwrapping results have been published there is a clear need for discussions on concrete methods for seismic phase unwrapping we have developed the graph cut phase unwrapping method which performs well in the interferometric synthetic aperture radar image processing it has advantages of strong discontinuity preserving ability and high computing efficiency to make it suitable for 3d seismic phase unwrapping the method is improved by extending it from 2d to 3d and by introducing the seismic horizon and unconformity constraints the strong and continuous conformable seismic events which can be easily tracked by certain autopicking methods are introduced as horizon constraints for guiding the phase unwrapping to ensure a constant

internal structure of the san jacinto fault zone at jackass flat from data recorded by a dense linear array the internal structure of the clark fault in the trifurcation area of the san jacinto fault zone is imaged using seismograms recorded by a dense linear array jackass flat jf crossing the surface trace of the fault and an adjacent array tr to the sw delay times between phase arrivals associated with 3500 local earthquakes and nine teleseismic events are used to estimate velocity variations within the arrays the teleseismic p waves travel faster beneath the tr than the jf array in contrast to larger scale tomographic results statistical analysis of local p wave delay times indicates that the entire jf array with an aperture of 400 m is inside a low velocity damage zone this low velocity zone is bounded on the ne side by a shallow bimaterial interface generating fault zone head waves and it contains an inner zone of more intense damage generating fault zone trapped waves the p wave velocity contrast across the local bounding bimaterial interface is 10 15 per cent the trapping structure is associated with a width of 200 m s wave velocity reduction of 35 per cent with respect to the surrounding rock q value of 20 and depth of 3 5 km the imaging results suggest that the main seismogenic fault is near the sw end of the jf array in agreement with a prominent geomorphologic feature the existence of intense local damage on the crustal block with faster larger scale velocity at depth is consistent with common propagation of earthquake ruptures in the area to the nw

fault zone guided wave generation on the locked late interseismic alpine fault new zealand fault zone guided waves fzgws have been observed for the first time within new zealand s transpressional continental plate boundary the alpine fault which is late in its typical seismic cycle ongoing study of these phases provides the opportunity to monitor interseismic conditions in the fault zone distinctive dispersive seismic codas 7 35 hz have been recorded on shallow borehole seismometers installed within 20 m of the principal slip zone near the central alpine fault known for low background seismicity fzgw generating microseismic events are located beyond the catchment scale partitioning of the fault indicating lateral connectivity of the low velocity zone immediately below the near surface segmentation initial modeling of the low velocity zone indicates a waveguide width of 60 200 m with a 10 40 reduction in s wave velocity similar to that inferred for the fault core of other mature plate boundary faults such as the san andreas and north anatolian faults

low velocity zones along the san jacinto fault southern california from body waves recorded in dense linear arrays we derive high resolution information on low velocity fault zone fz structures along the san jacinto fault zone sjfz southern california using waveforms of local earthquakes that are recorded at multiple linear cross fault arrays we observe clear across fault delays of direct p and s waves indicating damage zones at different segments of the sjfz we then compute synthetic traveltimes and waveforms using generalized ray theory and perform forward modeling to constrain the fz parameters at the southern section near the trifurcation area the low velocity zone lvz of the clark branch has a width of 200 m 30 45 reduction in vp and 50 reduction in vs from array data across the anza seismic gap we find a lvz with 200 m width and 50 reduction in both vp and vs nearly as prominent as that on the southern section we only find prominent lvzs beneath three out of the five arrays indicating along strike variations of the fault damage fz reflected phases are considerably less clear than those observed above the rupture zone of the 1992 landers earthquake shortly after the event this may reflect partially healed lvzs with less sharp boundaries at the sjfz given the relatively long lapse time from the last large surface rupturing event alternatively the lack of observed fz reflected phases could be partially due to the relatively small aperture of the arrays nevertheless the clear signatures of damage zones at anza and other locations indicate very slow healing process at least in the top few kilometers of the crust

palaeoecology of corals and rudists in mixed volcaniclastic carbonate small scale rhythms upper cretaceous jamaica abstract the guinea corn formation central jamaica west indies is represented by alternating platform carbonates and volcaniclastic sandstones and siltstones and represents the most landward shift of facies in an upper campanian to maastrichtian transgressive regressive cycle the succession consists of 2 30 m thick rhythms comprising a lower volcaniclastic siltstone sandstone division a middle volcaniclastic siltstone limestone with abundant corals or the rudist antillocaprina stellata and an upper division of rudist rudstones and floatstones the lower division consists of either bioturbated sandy siltstones with lignite and abundant gastropods and infaunal bivalves lower rhythms in the guinea corn formation or interbedded volcaniclastic siltstones and normally graded sandstones the rudist limestones are parallel bedded and consist of high abundance mono to paucispecific assemblages the boundary between the middle and upper parts of rhythms contains a concentration of oncoids with coral nuclei and a mantle of microbial laminates red algae and serpulids four coral assemblages are recognised based on diversity and coral morphology the paracycloseris dasmosmilia assemblage occurs in the lower division of the rhythms and is interpreted as a soft substrate assemblage that was adapted to elevated nutrient levels the ovalastrea actinacis actinacis multicolumnastrea and actinacis calamophyllia gyrodendron assemblages occur in the middle parts of the rhythms and consist of bedded rudstones rudstone mounds and rarer platestones pillarstones and mixstones polyparia arrangements colony morphology and low diversity suggest the ovalastrea actinacis and actinacis multicolumnastrea assemblages were adapted to high sedimentation rates the more diverse actinacis calamophyllia gyrodendron assemblage is interpreted to have grown under more optimum conditions low sedimentation rates and high light intensities rudists are classified on their growth orientation elevator clinger recumbent and on their communal relationships isolated or clustered the cluster elevator biradiolites mooretownensis occurs in the lower division of rhythms and the recumbent a stellata occurs in the middle part of rhythms the upper part of rhythms contains abundant rudists including isolated elevators antillocaprina occidentalis cluster elevators bournonia spp biradiolites jamaicensis thyrastylon spp and chiapasella radiolitiformis and clingers plagioptychus spp with the topmost part contains the large recumbent titanosarcolites

earthquake rupture extents and coseismic slips promoted by damaged fault zones here we investigate the effects of damage fault zones on rupture propagation by conducting a series of 3 d dynamic rupture simulations on a planar vertical strike slip fault we find that damage fault zones can promote rupture extent and increase earthquake potency the waves reflected from the bottom of shallow damage fault zones can increase shear stress on the fault and thus promote rupture propagation in addition the promotional effects increase with the width and depth extent of damage fault zones the overall effects of the waves reflected from the fault parallel side boundaries of damage fault zones are unfavorable for rupture propagation therefore rupture propagation is promoted with the increased width of fault zones due to geometrical spreading effects moreover nonground breaking ruptures may reach the ground surface with the effects of damage fault zones furthermore along strike segmented fault zones as suggested by observations could also promote ruptures and may lead to preferred rupture directions if epicenters are close to fault zones the effects of damage fault zones on rupture propagation hold important implications on assessing earthquake risk

seismic tomography of the southern california plate boundary region from noise based rayleigh and love waves we use cross correlations of ambient seismic noise between pairs of 158 broadband and short period sensors to investigate velocity structure over the top 5 10 km of the crust in the southern california plate boundary region around the san jacinto fault zone sjfz from the 9 component correlation tensors associated with all station pairs we derive dispersion curves of rayleigh and love wave group velocities the dispersion results are inverted first for rayleigh and love waves group velocity maps on a 1 5 1 5 km2 grid that includes portions of the sjfz the san andreas fault saf and the elsinore fault we then invert these maps to 3d shear wave velocities in the top 7 km of the crust the distributions of the rayleigh and love group velocities exhibit 2 azimuthal anisotropy with fast directions parallel to the main faults and rotations in complex areas the reconstructed 3d shear velocity model reveals complex shallow structures correlated with the main geological units and strong velocity contrasts across various fault sections along with low velocity damage zones and basins the sjfz is marked by a clear velocity contrast with higher v s values on the ne block for the section se of the san jacinto basin and a reversed contrast across the section between the san jacinto basin and the saf velocity contrasts are also observed along the southern parts on the saf and the elsinore fault with a faster southwest block in both cases the region around the salton trough is associated with a significant low velocity zone strong velocity reductions following flower shape with depth are observed extensively around both the sjfz and the saf and are especially prominent in areas of geometrical complexity in particular the area between the sjfz and the saf is associated with an extensive low velocity zone correlated with diffuse seismicity at depth and a similar pattern including correlation with deep diffuse seismicity is observed on a smaller scale in the trifurcation area of the sjfz these results augment local earthquake tomography images that have low resolution in the top few km of the crust and provide important constraints for studies concerned with behavior of earthquake ruptures generation of rock damage and seismic shaking hazard in the region

distinct asymmetry in rupture induced inelastic strain across dipping faults an off fault yielding model 1 off fault stresses associated with rupture propagation cause material yielding under a mohr coulomb condition leaving irrecoverable deformation in the fault zone 2d dynamic rupture simulations for a 30 reverse fault and a 60 normal fault in a depth dependent stress environment show that the yielding inelastic zone off the fault broadens as it nears the surface with decreasing confining pressure forming a skewed flower like structure bounded at the top by the free surface the inelastic strain in the hanging wall is significantly larger and broader than the footwall for both reverse and normal faults the occurrence of inelastic strain however significantly reduces ground motion especially on the hanging wall and gives rise to a reduced asymmetry in ground motion on the hanging wall and footwall compared to elastic solutions these results provide theoretical predictions for fault zone structure of dipping faults that can be tested by future observational experiments

quantifying near field and off fault deformation patterns of the 1992 mw 7 3 landers earthquake coseismic surface deformation in large earthquakes is typically measured using field mapping and with a range of geodetic methods e g insar lidar differencing and gps current methods however either fail to capture patterns of near field coseismic surface deformation or lack preevent data consequently the characteristics of off fault deformation and the parameters that control it remain poorly understood we develop a standardized method to fully measure the surface near field coseismic deformation patterns at high resolution using the cosi corr program by correlating pairs of aerial photographs taken before and after the 1992 m w 7 3 landers earthquake cosi corr offers the advantage of measuring displacement across the entire zone of surface deformation and over a wider aperture than that available to field geologists for the landers earthquake our measured displacements are systematically larger than the field measurements indicating the presence of off fault deformation we show that 46 of the total surface displacement occurred as off fault deformation over a mean deformation width of 154 m the magnitude and width of off fault deformation along the rupture is primarily controlled by the macroscopic structural complexity of the fault system with a weak correlation with the type of near surface materials through which the rupture propagated both the magnitude and width of distributed deformation are largest in stepovers bends and at the southern termination of the surface rupture we find that slip along the surface rupture exhibits a consistent degree of variability at all observable length scales and that the slip distribution is self affine fractal with dimension of 1 56

seismic velocity structures in the southern california plate boundary environment from double difference tomography summary we present tomographic images of crustal structures in the southern california plate boundary area with a focus on the san jacinto fault zone sjfz based on double difference inversions of earthquake arrival times absolute arrival times of 247 472 p and 105 448 s wave phase picks for 5493 earthquakes recorded at 139 stations in southern california are used starting with a layered 1 d model and continuing in later iterations with various updated initial models we invert the data for vp and vs in a 270 km long 105 km wide and 35 km deep volume around the sjfz using a spatially variable grid with higher density around the sjfz the examined volume stretches from cajon pass to the northernmost imperial fault zone and includes portions of the southern san andreas fault saf the elsinore fault and the brawley seismic zone in the salton trough because differential traveltimes used in the double difference inversions are most sensitive to near source structures we obtain high resolution around the earthquake sources after 30 iterations we improve the average traveltime misfit by a factor of 16 though ray coverage is limited at shallow depths we obtain detailed images of seismic velocities from 3 to 20 km throughout much of the study area our final velocity results show zones of low velocity and anomalous vp vs ratios associated with various fault strands and sedimentary basins along with clear velocity contrasts across the sjfz and the southern saf the velocity reductions in fault zone regions are generally highest in geometrically complex areas up to 30 50 per cent in the top few kilometres are higher for vs than for vp and follow a flower type pattern with depth in the central section of the sjfz from the san jacinto valley to the trifurcation area the northeast side of the fault has generally higher seismic velocities than the southwest block the obtained contrasts of vp are more persistent and higher up to 20 per cent than the contrasts of vs up to 15 per cent although the differences may stem at least partially from the higher resolution of vp images in the sjfz sections to the northwest and southeast there are patches with reversed velocities contrasts especially in the shallow crust near the san jacinto valley and other basins along the banning fault there is no clear velocity contrast for the southern saf the northeast side has generally lower seismic velocities in the seismogenic zone with patches of contrast reversals in the shallow crust in the brawley seismic zone the northeast side has somewhat lower velocities in the 20 km section near the southern saf and higher velocities farther to the southwest the imaged features have important implications for various aspects of earthquake and crustal dynamics in the region

rupture directivity of microearthquake sequences near parkfield california 1 the direction of propagation is an important factor that affects the pattern of ground motion generated by an earthquake characterizing factors favoring a potential rupture propagation direction is thus an important task here we analyze the earthquake directivity of repeating earthquake sequences located on the san andreas fault near parkfield california all earthquakes of a sequence have very similar waveforms and have overlapping surface ruptures we show that subtle variations of the transfer function between earthquakes of a common sequence can be interpreted as a change of apparent rupture duration relative apparent rupture durations are computed for all pairs of events at all available stations and for each sequence we invert these measurements to obtain an estimation of the apparent rupture duration for each individual event of the sequence relative to a reference event variation of apparent rupture duration with azimuth attests for the rupture directivity we show that the majority of analyzed microearthquakes presents a rupture in the south east direction we also show that on a given repeating sequence most earthquakes tend to show the same rupture direction

dynamic ruptures on a frictional interface with off fault brittle damage feedback mechanisms and effects on slip and near fault motion the spontaneous generation of brittle rock damage near and behind the tip of a propagating rupture can produce dynamic feedback mechanisms that modify significantly the rupture properties seismic radiation and generated fault zone structure in this work we study such feedback mechanisms for single rupture events and their consequences for earthquake physics and various possible observations this is done through numerical simulations of in plane dynamic ruptures on a frictional fault with bulk behavior governed by a brittle damage rheology that incorporates reduction of elastic moduli in off fault yielding regions the model simulations produce several features that modify key properties of the ruptures local wave propagation and fault zone damage these include 1 dynamic generation of near fault regions with lower elastic properties 2 dynamic changes of normal stress on the fault 3 rupture transition from crack like to a detached pulse 4 emergence of a rupture mode consisting of a train of pulses 5 quasi periodic modulation of slip rate on the fault and 6 asymmetric near fault ground motion with higher amplitude and longer duration on the side with reduced elastic moduli the results can have significant implications to multiple topics ranging from rupture directivity and local amplification of seismic motion to near fault tremor like signals

shallow low velocity zone of the san jacinto fault from local earthquake waveform modelling summary we developed a method to determine the depth extent of low velocity zone lvz associated with a fault zone fz using s wave precursors from local earthquakes the precursors are diffracted s waves around the edges of lvz and their relative amplitudes to the direct s waves are sensitive to the lvz depth we applied the method to data recorded by three temporary arrays across three branches of the san jacinto fz the fz dip was constrained by differential traveltimes of p waves between stations at two side of the fz other fz parameters width and velocity contrast were determined by modelling waveforms of direct and fz reflected p and s waves we found that the lvz of the buck ridge fault branch has a width of 150 m with a 30 40 per cent reduction in vp and a 50 60 per cent reduction in vs the fault dips 70 5 to southwest and its lvz extends only to 2 1 km in depth the lvz of the clark valley fault branch has a width of 200 m with 40 per cent reduction in vp and 50 per cent reduction in vs the coyote creek branch is nearly vertical and has a lvz of 150 m in width and of 25 per cent reduction in vp and 50 per cent reduction in vs the lvzs of these three branches are not centred at the surface fault trace but are located to their northeast indicating asymmetric damage during earthquakes

waveform relocated earthquake catalog for southern california 1981 to june 2011 we determine a new relocated catalog hys catalog 2011 for southern california from 1981 through june 2011 about 75 3 of the hypocenters are calculated with absolute and differential travel time picks and 24 7 could be relocated only by using absolute travel time picks with 3d or 1d velocity models the total catalog consists of more than 502 000 earthquakes in the region extending from baja california in the south to coalinga and owens valley in the north the catalog consists of three m 7 1 m 7 2 and m 7 3 mainshocks their foreshocks and aftershocks and background seismicity caused by tectonic and other processes in the southern california crust hypocenters in the new relocated catalog exhibit tighter spatial clustering of seismicity than does the routinely generated catalog and the depth distribution is tighter and reflects the thickness of the seismogenic zone more accurately compared to the standard catalog the relocated hypocenters are more easily related to other data sets such as mapped late quaternary faults

a methodological approach towards high resolution surface wave imaging of the san jacinto fault zone using ambient noise recordings at a spatially dense array we present a new technique for deriving detailed information on seismic velocities of the subsurface material from continuous ambient noise recorded by spatially dense seismic arrays this method uses iterative double beamforming between various subarrays to extract surface wave contributions from the ambient noise data in complex environments with unfavourable noise source distributions the iterative double beamforming extraction makes it possible to retrieve large amounts of rayleigh wave traveltime information in a wide frequency band the method is applied to data recorded by a highly dense nodal array with 1108 vertical geophones centred on the damage zone of the clark branch of the san jacinto fault zone south of anza california the array covers a region of 650 700 m 2 with instrument spacing of 10 30 m and continuous recording at 500 samples s 1 over 30 d in 2014 using this iterative double beamforming on subarrays of 25 sensors and cross correlations between all of the station pairs we separate surface waves from body waves that are abundant in the raw cross correlation data focusing solely on surface waves maps of traveltimes are obtained at different frequencies with unprecedented accuracy at each point of a 15 m spacing grid group velocity inversions at 2 4 hz reveal depth and lateral variations in the structural properties within and around the san jacinto fault zone in the study area this method can be used over wider frequency ranges and can be combined with other imaging techniques such as eikonal tomography to provide unprecedented detailed structural images of the subsurface material

lidar and field observations of slip distribution for the most recent surface ruptures along the central san jacinto fault we measured offsets on tectonically displaced geomorphic features along 80 km of the clark strand of the san jacinto fault sjf to estimate slip per event for the past several surface ruptures we identify 168 offset features from which we make over 490 measurements using b4 light detection and ranging lidar imagery and field observations our results suggest that lidar technology is an exemplary supplement to traditional field methods in slip per event studies displacement estimates indicate that the most recent surface rupturing event mre produced an average of 2 5 2 9 m of right lateral slip with maximum slip of nearly 4 m at anza a m w 7 2 7 5 earthquake average multiple event offsets for the same 80 kms are 5 5 m with maximum values of 3 m at anza for the penultimate event cumulative displacements of 9 10 m through anza suggest the third event was also similar in size paleoseismic work at hog lake dates the most recent surface rupture event at ca 1790 a poorly located large earthquake occurred in southern california on 22 november 1800 we relocate this event to the clark fault based on the mre at hog lake we also recognize the occurrence of a younger rupture along 15 20 km of the fault in blackburn canyon with 1 25 m of average displacement we attribute these offsets to the 21 april 1918 m w 6 9 event these data argue that much or all of the clark fault and possibly also the casa loma fault fail together in large earthquakes but that shorter sections may fail in smaller events online material topographic contour maps and hillshades generated from b4 lidar data corresponding field photographs and data tables comparing lidar based and field based slip measurements for individual geomorphic features along the clark fault

predictive painting of 3d seismic volumes predictive painting is a numerical algorithm that spreads information in 3d volumes according to the local structure of seismic events the algorithm consists of two steps first local spatially variable inline and crossline slopes of seismic events are estimated by the plane wave destruction method next a seed trace is inserted in the volume and the information contained in that trace is spread inside the volume thus automatically painting the data space immediate applications of this technique include automatic horizon picking and flattening in applications to both prestack and poststack seismic data analysis synthetic and field data tests demonstrate the effectiveness of predictive painting

thickness of the seismic slip zone this article reviews geologic and other evidence constraining the thickness of the principal slip zone psz that accommodates the bulk of coseismic shear displacement during an individual rupture event surface deformation from rupturing may occupy swaths tens of meters or more in width but trenches across active faults generally reveal that incremental slip is accommodated by a psz that is tens of centimeters or less in thickness geomorphic evidence coupled with the observations from trenching suggest a psz may stay well localized for distances of several kilometers through many rupture episodes mine exposures and exhumed fault zones demonstrate that pszs separating different lithologies within the fault core although contained within damage zones of variably fractured rock ranging up to hundreds of meters in thickness often comprise just a few centimeters of gouge ultracataclasite that have accommodated large finite displacements 1 km microstructural studies demonstrate incremental slip localized still further down to 1 10 mm as do other fault rock assemblages slickensides and slickenfibers fault veins of pseudotachylyte friction melt intravein septa in hydrothermal fault infills the accumulated evidence indicates that localization of coseismic shearing to less than 10 cm on planar faults is widespread throughout the crustal seismogenic zone with extreme localization to less than 1 cm not uncommon however some distributed coseismic shear may also develop especially at rupture irregularities coseismic reduction of shear resistance from friction melting t 1000 c or from transient thermal pressurization of aqueous fluids t 100 c requires slip during moderate to large earthquakes u 1 m to be restricted to narrow zones respectively a few centimeters or tens of centimeters in thickness given the evidence for slip localization the apparent scarcity of pseudotachylyte suggests either that seismic friction melting is a rare phenomenon or that pseudotachylyte is only rarely preserved in recognizable form within mature hydrated fault zones

stratigraphic surfaces in the depositional and chronostratigraphic wheeler transformed domain one of the controversial and still debated issues of sequence stratigraphy is the assessment of stratigraphic surfaces in a chronostratigraphic framework catuneanu 2006 to assess whether the bounding surfaces of systems tracts or sequences are isochronous or diachronous is of paramount importance for stratigraphic correlation new software developments enable us to study stratigraphic surfaces within such a chronostratigraphic framework and allow us to predict something about their time attributes

relative geologic time age volumes relating every seismic sample to a geologically reasonable horizon the age volume described in this paper can significantly increase the amount of useful geologic information extracted from a seismic data volume this paper has three goals the first is to introduce and describe the relative geologic time volume referred to herein as either an rgt volume or an age volume every seismic sample in an age volume contains an estimate of the relative geologic time for that sample which automatically relates every sample to a geologic horizon the second goal is to discuss some of the applications and advantages of using a rgt volume to store and access seismic interpretation for example it can be used to sculpt a seismic volume along a desired horizon or view depositional changes as a function of geologic time the third goal is to describe two methods for generating such a volume the first employing standard interpretation techniques and the second relying on unwrapping instantaneous phase unwrapping the instantaneous phase is a new holistic approach to interpreting seismic data that can produce a very detailed seismic interpretation the goal of building and using an rgt volume is to extract significantly more geologic information from the seismic data than is available from conventional methods the rgt volume is both a way of storing and retrieving seismic horizons to generate such a volume each seismic sample needs to be assigned an estimate of relative geologic time which in essence connects it to a geologically reasonable horizon in normal sections these relative geologic time values will increase with depth or traveltime high gradients in the rgt volume could be indicative of faults unconformities fluid contacts or interpretation errors while a constant value surface in a rgt volume is a seismic horizon figure 1 contains a portion of the seismic data used to illustrate the

rupture and ground motion models on the northern san jacinto fault incorporating realistic complexity we use the 3d finite element method to conduct dynamic models of rupture and resulting ground motion on the claremont casa loma stepover of the northern san jacinto fault we incorporate complex fault geometry from the u s geological survey usgs quaternary faults database see data and resources 1 a realistic velocity structure the southern california earthquake center community velocity model s a realistic regional stress field with an orientation taken from seismicity relocation literature and several stochastic self similar shear stress distributions as we incorporate more types of complexity the specific effects of any individual factor become less apparent within the overall rupture behavior we also find that the distribution of high and low shear stress that arises from combining regional and stochastic stress fields has the strongest control over where the rupture terminates using a regional stress field alone as well as with the combined regional and stochastic stress realization we find that the stepover presents a significant barrier to rupture regardless of our choice of initial nucleation point and that it is difficult for rupture to propagate the full length of either fault segment greater heterogeneity of stresses tends to produce shorter ruptures within this result we find that the claremont strand is more favorable for long ruptures than the casa loma clark strand low frequency ground motion intensity and distribution are controlled largely by the velocity structure rather than by stress heterogeneity the strongest motions produced in these models are in the san bernardino basin although directivity effects do contribute to the low frequency ground motion distribution particularly in the near field they are secondary to the effects of the velocity structure online material figures of ground motions from models used to calibrate the stress conditions for dynamic rupture propagation 1 sec 18

three dimensional calculations of fault zone guided waves in various irregular structures summary the detailed structure of fault zones fzs plays an important role in problems related to fault mechanics earthquake rupture wave propagation and seismic hazard fzs are thought to consist of an o 10 100 m wide region of decreased seismic velocity but structural details such as their depth extent lateral and vertical variations etc are elusive the small spatial scales involved make such structures difficult to image with ray theoretical methods such as tomography however seismic energy trapped inside fz layers can provide dispersive wave trains that carry information on the fz structure these waves can travel many kilometres inside the fz before reaching the surface and are therefore strongly altered by its properties candidate trapped waves have been observed above several active faults inversion algorithms exist that can model these observations in terms of planar fault zone structures however at present it is not clear how reliable these estimates are as the effects of even small 3 d variations on trapped waves are not well understood the goal of this study is to distinguish 3 d structures that do and do not significantly affect fz waves to achieve this we perform numerical calculations of wave propagation in various fz geometries and analyse the waveforms spectra and envelopes of the synthetic seismograms the main results are that 1 moderate changes of the shape of fz or 2 small scale heterogeneities or 3 depth dependent properties do not strongly affect the observed fz waves in contrast strong effects are to be expected from 4 breaks in the continuity of fz structure e g offsets which may at some point allow imaging of such features at depth

location of largest earthquake slip and fast rupture controlled by along strike change in fault structural maturity due to fault growth earthquake slip distributions are asymmetric along strike but the reasons for the asymmetry are unknown we address this question by establishing empirical relations between earthquake slip profiles and fault properties we analyze the slip distributions of 27 large continental earthquakes in the context of available information on their causative faults in particular on the directions of their long term lengthening we find that the largest slips during each earthquake systematically occurred on that half of the ruptured fault sections most distant from the long term fault propagating tips i e on the most mature half of the broken fault sections meanwhile slip decreased linearly over most of the rupture length in the direction of long term fault propagation i e of decreasing structural maturity along strike we suggest that this earthquake slip asymmetry is governed by along strike changes in fault properties including fault zone compliance and fault strength induced by the evolution of off fault damage fault segmentation and fault planarity with increasing structural maturity we also find higher rupture speeds in more mature rupture sections consistent with predicted effects of low velocity damage zones on rupture dynamics since the direction s of long term fault propagation can be determined from geological evidence it might be possible to anticipate in which direction earthquake slip once nucleated may increase accelerate and possibly lead to a large earthquake our results could thus contribute to earthquake hazard assessment and earthquake early warning

internal structure of the san jacinto fault zone at blackburn saddle from seismic data of a linear array local and teleseismic earthquake waveforms recorded by a 180 m long linear array bb with seven seismometers crossing the clark fault of the san jacinto fault zone northwest of anza are used to image a deep bimaterial interface and core damage structure of the fault delay times of p waves across the array indicate an increase in slowness from the southwest most bb01 to the northeast most bb07 station automatic algorithms combined with visual inspection and additional analyses are used to identify local events generating fault zone head and trapped waves the observed fault zone head waves imply that the clark fault in the area is a sharp bimaterial interface with lower seismic velocity on the southwest side the moveout between the head and direct p arrivals for events within 40 km epicentral distance indicates an average velocity contrast across the fault over that section and the top 20 km of 3 2 per cent a constant moveout for events beyond 40 km to the southeast is due to off fault locations of these events or because the imaged deep bimaterial interface is discontinuous or ends at that distance the lack of head waves from events beyond 20 km to the northwest is associated with structural complexity near the hemet stepover events located in a broad region generate fault zone trapped waves at stations bb04 bb07 waveform inversions indicate that the most likely parameters of the trapping structure are width of 200 m s velocity reduction of 30 40 per cent with respect to the bounding blocks q value of 10 20 and depth of 3 5 km the trapping structure and zone with largest slowness are on the northeast side of the fault the observed sense of velocity contrast and asymmetric damage across the fault suggest preferred rupture direction of earthquakes to the northwest this inference is consistent with results of other geological and seismological studies

nonplanar faults mechanics of slip and off fault damage stress interactions and sliding characteristics of faults with random fractal waviness in a purely elastic medium differ both qualitatively and quantitatively from those of faults with planar surfaces with nonplanar fault models solutions for slip diverge as resolution of the fractal features increases and the scaling of fault slip with fault rupture dimension becomes nonlinear we show that the nonlinear scaling of slip and divergence of solutions arise because stresses from geometric interactions at irregularities along nonplanar faults grow with increasing slip and produce backstresses that progressively impede slip however in real materials with finite strength yielding will halt the growth of the interaction stresses which will profoundly affect slip of nonplanar faults we infer that in the brittle seismogenic portion of the earth s crust off fault yielding occurs on pervasive secondary faults predicted rates of stress relaxation with distance from major faults with random fractal roughness follow a power law relationship that is consistent with reported clustering of background seismicity up to 15 kilometers from faults

fully automated carbonate petrography using deep convolutional neural networks abstract carbonate rocks are important archives of past ocean conditions as well as hosts of economic resources such as hydrocarbons water and minerals geologists typically perform compositional analysis of grain matrix cement and pore types in order to interpret depositional environments diagenetic modification and reservoir quality of carbonate strata such information can be obtained primarily from petrographic analysis a task that is costly labor intensive and requires in depth knowledge of carbonate petrology and micropaleontology recent studies have leveraged machine learning based image analysis including deep convolutional neural networks dcnn to automate description classification and interpretation of thin sections subsurface core images and seismic facies which would accelerate data acquisition and reproducibility for these tasks in carbonate rocks this approach has been applied primarily to recognize carbonate lithofacies and no attempt has been made to individually identify and quantify various types of carbonate grains matrix and cement in this study the applicability and performance of dcnn based object detection and image classification approaches are assessed with respect to carbonate compositional analysis the training data comprised of more than 13 000 individually labelled objects from nearly 4000 carbonate petrographic images the dataset is grouped into six and nine different classes for the image classification and object detection tasks respectively even with a small and relatively imbalanced training set the dcnn was able to achieve an f1 score of 92 for image classification and 84 mean precision for object detection by combining one cycle policy class weight and label mixup smoothing methods this study highlights the inefficiency of image classification as an approach to replicating human description and classification of carbonate petrography by contrast dcnn based object detection appears capable of approaching human speed and accuracy in the area of carbonate petrography because it is able to individually locate and identify different carbonate components with greater cost efficiency speed and reproducibility than conventional human petrographic analysis

can seismic waves be trapped inside an inactive fault zone the case study of nocera umbra central italy in march 1998 a small aperture 350 m linear array of five broadband stations was installed across a pliocene strike slip fault crossing the town of nocera umbra central italy three component seismograms of local earthquakes recorded during the umbria marche seismic sequence have shown that the degradation of mechanical properties in damaged rock causes a significant amplification of ground motion within the fault zone marra et al 2000 during the operation of the array an m w 5 3 subcrustal earthquake occurred at a depth of 48 km and an epicentral distance of 10 km from nocera umbra the near vertical incidence of this event offers a unique opportunity for constraining the 2d geometry of the laterally heterogeneous low velocity structure responsible for the large amplification across the fault based on the observed symmetry of the amplification pattern a 2d wedge shaped model of the fault zone is used for the computation of antiplane sh synthetic waveforms to be compared with observations along the array two parameters of the wedge model can be constrained by independent geophysical and geological data 1 shear wave velocity in pristine rock outside the fault zone v s 2200 m sec and 2 the lateral extent of the damaged rock zone l 160 m the other two model parameters namely the wedge depth h and shear wave velocity of the fault zone v f s are varied computer simulations show that trapped waves which have been observed so far for earthquakes occurring inside the fault zone can develop inside a 2d wedge model when struck by external vertically incident seismic radiation results of numerical modeling indicate that in spite of its extreme simplicity such an idealized model is successful in reproducing many of the essential characters of ground motion within the fault zone of nocera umbra even though the best fit solution is not unique and the model parameters trade off each other a study of parameter sensitivity comparing synthetics with recorded data indicates that the variability of h and v f s can be restricted to a limited range of values the quality of fit is optimal when the shear velocity reduction lies in the range of 40 50 and the depth of the trapping zone is comprised between 1 and 2 km