Geology Reports⌕ Search

Geology topics

W.L. Ellsworth

Publications and source records attributed to W.L. Ellsworth.

At least 19 recordsLinked to original sources
USGS2010

Locating non-volcanic tremor along the San Andreas Fault using a multiple array source imaging technique

Non-volcanic tremor (NVT) has been observed at several subduction zones and at the San Andreas Fault (SAF). Tremor locations are commonly derived by cross-correlating envelope-transformed seismic traces in combination with source-scanning techniques. Recently, they have also been located by using relative relocations with master events, that is low-frequency earthquakes that are part of the tremor; locations are derived by conventional traveltime-based methods. Here we present a method to locate the sources of NVT using an imaging approach for multiple array data. The performance of the method is checked with synthetic tests and the relocation of earthquakes. We also applied the method to tremor occurring near Cholame, California. A set of small-aperture arrays (i.e. an array consisting of arrays) installed around Cholame provided the data set for this study. We observed several tremor episodes and located tremor sources in the vicinity of SAF. During individual tremor episodes, we observed a systematic change of source location, indicating rapid migration of the tremor source along SAF. ?? 2010 The Authors Geophysical Journal International ?? 2010 RAS.

Geophysical Journal International↗
USGS2010

Precise estimation of repeating earthquake moment: Example from parkfield, california

We offer a new method for estimating the relative size of repeating earthquakes using the singular value decomposition (SVD). This method takes advantage of the highly coherent waveforms of repeating earthquakes and arrives at far more precise and accurate descriptions of earthquake size than standard catalog techniques allow. We demonstrate that uncertainty in relative moment estimates is reduced from ??75% for standard coda-duration techniques employed by the network to an uncertainty of ??6.6% when the SVD method is used. This implies that a single-station estimate of moment using the SVD method has far less uncertainty than the whole-network estimates of moment based on coda duration. The SVD method offers a significant improvement in our ability to describe the size of repeating earthquakes and thus an opportunity to better understand how they accommodate slip as a function of time.

Bulletin of the Seismological Society of America↗
USGS2009

Precise location of San Andreas Fault tremors near Cholame, California using seismometer clusters: Slip on the deep extension of the fault?

We examine a 24-hour period of active San Andreas Fault (SAF) tremor and show that this tremor is largely composed of repeated similar events. Utilizing this similarity, we locate the subset of the tremor with waveforms similar to an identified low frequency earthquake (LFE) "master template," located using P and S wave arrivals to be ???26 km deep. To compensate for low signal-to-noise, we estimate event-pair differential times at "clusters" of nearby stations rather than at single stations. We find that the locations form a near-linear structure in map view, striking parallel to the SAF and near the surface trace. Therefore, we suggest that at least a portion of the tremor occurs on the deep extension of the fault, likely reflecting shear slip, similar to subduction zone tremor. If so, the SAF may extend to the base of the crust, ???10 km below the deepest regular earthquakes on the fault. ?? 2009 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2007

Source parameters of a M4.8 and its accompanying repeating earthquakes off Kamaishi, NE Japan: Implications for the hierarchical structure of asperities and earthquake cycle

We determine the source parameters of a M4.9 ?? 0.1 'characteristic earthquake' sequence and its accompanying microearthquakes at ???50 km depth on the subduction plate boundary offshore of Kamaishi, NE Japan. The microearthquakes tend to occur more frequently in the latter half of the recurrence intervals of the M4.9 ?? 0.1 events. Our results show that the microearthquakes are repeating events and they are located not only around but also within the slip area for the 2001 M4.8 event. From the hierarchical structure of slip areas and smaller stress drops for the microearthquakes compared to the M4.8 event, we infer the small repeating earthquakes rupture relatively weak patches in and around the slip area for the M4.8 event and their activity reflects a stress concentration process and/or change in frictional property (healing) at the area. We also infer the patches for the M4.9 ?? 0.1 and other repeating earthquakes undergo aseismic slip during their interseismic period. Copyright 2007 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2007

Monitoring microearthquakes with the San Andreas fault observatory at depth

In 2005, the San Andreas Fault Observatory at Depth (SAFOD) was drilled through the San Andreas Fault zone at a depth of about 3.1 km. The borehole has subsequently been instrumented with high-frequency geophones in order to better constrain locations and source processes of nearby microearthquakes that will be targeted in the upcoming phase of SAFOD. The microseismic monitoring software MIMO, developed by NORSAR, has been installed at SAFOD to provide near-real time locations and magnitude estimates using the high sampling rate (4000 Hz) waveform data. To improve the detection and location accuracy, we incorporate data from the nearby, shallow borehole (???250 m) seismometers of the High Resolution Seismic Network (HRSN). The event association algorithm of the MIMO software incorporates HRSN detections provided by the USGS real time earthworm software. The concept of the new event association is based on the generalized beam forming, primarily used in array seismology. The method requires the pre-computation of theoretical travel times in a 3D grid of potential microearthquake locations to the seismometers of the current station network. By minimizing the differences between theoretical and observed detection times an event is associated and the location accuracy is significantly improved.

Conference Paper↗
USGS2004

Streaks, multiplets, and holes: High-resolution spatio-temporal behavior of Parkfield seismicity

Double-difference locations of ???8000 earthquakes from 1969-2002 on the Parkfield section of the San Andreas Fault reveal detailed fault structures and seismicity that is, although complex, highly organized in both space and time. Distinctive features of the seismicity include: 1) multiple recurrence of earthquakes of the same size at precisely the same location on the fault (multiplets), implying frictional or geometric controls on their location and size; 2) sub-horizontal alignments of hypocenters along the fault plane (streaks), suggestive of rheological transitions within the fault zone and/or stress concentrations between locked and creeping areas; 3) regions devoid of microearthquakes with typical dimensions of 1-5 km (holes), one of which contains the M6 1966 Parkfield earthquake hypocenter. These features represent long lived structures that persist through many cycles of individual event. Copyright 2004 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2003

Detailed fault structure of the 2000 Western Tottori, Japan, earthquake sequence

We investigate the faulting process of the aftershock region of the 2000 western Tottori earthquake (Mw 6.6) by combining aftershock hypocenters and moment tensor solutions. Aftershock locations were precisely determined by the double difference method using P- and S-phase arrival data of the Japan Meteorological Agency unified catalog. By combining the relocated hypocenters and moment tensor solutions of aftershocks by broadband waveform inversion of FREESIA (F-net), we successfully resolved very detailed fault structures activated by the mainshock. The estimated fault model resolves 15 individual fault segments that are consistent with both aftershock distribution and focal mechanism solutions. Rupture in the mainshock was principally confined to the three fault elements in the southern half of the zone, which is also where the earliest aftershocks concentrate. With time, the northern part of the zone becomes activated, which is also reflected in the postseismic deformation field. From the stress tensor analysis of aftershock focal mechanisms, we found a rather uniform stress field in the aftershock region, although fault strikes were scattered. The maximum stress direction is N107??E, which is consistent with the tectonic stress field in this region. In the northern part of the fault, where no slip occurred during the mainshock but postseismic slip was observed, the maximum stress direction of N130??E was possible as an alternative solution of stress tensor inversion.

Tottori↗
USGS2003

Imaging the complexity of an active normal fault system: The 1997 Colfiorito (central Italy) case study

Six moderate magnitude earthquakes (5 < Mw < 6) ruptured normal fault segments of the southern sector of the North Apennine belt (central Italy) in the 1997 Colfiorito earthquake sequence. We study the progressive activation of adjacent and nearby parallel faults of this complex normal fault system using ???1650 earthquake locations obtained by applying a double-difference location method, using travel time picks and waveform cross-correlation measurements. The lateral extent of the fault segments range from 5 to 10 km and make up a broad, ???45 km long, NW trending fault system. The geometry of each segment is quite simple and consists of planar faults gently dipping toward SW with an average dip of 40??-45??. The fault planes are not listric but maintain a constant dip through the entire seismogenic volume, down to 8 km depth. We observe the activation of faults on the hanging wall and the absence of seismicity in the footwall of the structure. The observed fault segmentation appears to be due to the lateral heterogeneity of the upper crust: preexisting thrusts inherited from Neogene's compressional tectonic intersect the active normal faults and control their maximum length. The stress tensor obtained by inverting the six main shock focal mechanisms of the sequence is in agreement with the tectonic stress active in the inner chain of the Apennine, revealing a clear NE trending extension direction. Aftershock focal mechanisms show a consistent extensional kinematics, 70% of which are mechanically consistent with the main shock stress field.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Apparent break in earthquake scaling due to path and site effects on deep borehole recordings

We reexamine the scaling of stress drop and apparent stress, rigidity times the ratio between seismically radiated energy to seismic moment, with earthquake size for a set of microearthquakes recorded in a deep borehole in Long Valley, California. In the first set of calculations, we assume a constant Q and solve for the corner frequency and seismic moment. In the second set of calculations, we model the spectral ratio of nearby events to determine the same quantities. We find that the spectral ratio technique, which can account for path and site effects or nonconstant Q, yields higher stress drops, particularly for the smaller events in the data set. The measurements determined from spectral ratios indicate no departure from constant stress drop scaling down to the smallest events in our data set (Mw 0.8). Our results indicate that propagation effects can contaminate measurements of source parameters even in the relatively clean recording environment of a deep borehole, just as they do at the Earth's surface. The scaling of source properties of microearthquakes made from deep borehole recordings may need to be reevaluated.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

A Brownian model for recurrent earthquakes

We construct a probability model for rupture times on a recurrent earthquake source. Adding Brownian perturbations to steady tectonic loading produces a stochastic load-state process. Rupture is assumed to occur when this process reaches a critical-failure threshold. An earthquake relaxes the load state to a characteristic ground level and begins a new failure cycle. The load-state process is a Brownian relaxation oscillator. Intervals between events have a Brownian passage-time distribution that may serve as a temporal model for time-dependent, long-term seismic forecasting. This distribution has the following noteworthy properties: (1) the probability of immediate rerupture is zero; (2) the hazard rate increases steadily from zero at t = 0 to a finite maximum near the mean recurrence time and then decreases asymptotically to a quasi-stationary level, in which the conditional probability of an event becomes time independent; and (3) the quasi-stationary failure rate is greater than, equal to, or less than the mean failure rate because the coefficient of variation is less than, equal to, or greater than 1/???2 ??? 0.707. In addition, the model provides expressions for the hazard rate and probability of rupture on faults for which only a bound can be placed on the time of the last rupture. The Brownian relaxation oscillator provides a connection between observable event times and a formal state variable that reflects the macromechanics of stress and strain accumulation. Analysis of this process reveals that the quasi-stationary distance to failure has a gamma distribution, and residual life has a related exponential distribution. It also enables calculation of "interaction" effects due to external perturbations to the state, such as stress-transfer effects from earthquakes outside the target source. The influence of interaction effects on recurrence times is transient and strongly dependent on when in the loading cycle step pertubations occur. Transient effects may be much stronger than would be predicted by the "clock change" method and characteristically decay inversely with elapsed time after the perturbation.

Bulletin of the Seismological Society of America↗
USGS2002

Comments on potential geologic and seismic hazards affecting Mare Island, Solano County, California

This report was prepared in response to a written request from the City of Vallejo, California, to the U.S. Geological Survey (USGS). By letter of October 4, 2002, the City requested that the USGS "provide advice to the City’s LNG Health and Safety Committee on its review of a potential liquid natural gas project" on the southern portion of Mare Island. The City specifically requested that the USGS advise the committee on potential hazards including fault rupture, earthquake ground motion, soil failure during earthquakes, tsunami and seiche, and landslides. The City requested that the USGS: (1) comment on these hazards, (2) describe its degree of confidence in its opinions, and (3) describe the scope of additional studies that will be needed if the City enters into an agreement with project sponsors. Advice was also requested on the selection of the safe shutdown and operating basis earthquakes as specified in the NFPA 59A standard (NFPA, 2001). This review of published reports and other publicly available information indicates that all of the hazards on which the USGS was asked to comment should be considered for the proposed project on the southern portion of Mare Island. Available information differs greatly for each of these potential hazards, and adequate understanding for design will require detailed site-specific investigations.

California↗
USGS2000

A Double-difference Earthquake location algorithm: Method and application to the Northern Hayward Fault, California

We have developed an efficient method to determine high-resolution hypocenter locations over large distances. The location method incorporates ordinary absolute travel-time measurements and/or cross-correlation P-and S-wave differential travel-time measurements. Residuals between observed and theoretical travel-time differences (or double-differences) are minimized for pairs of earthquakes at each station while linking together all observed event-station pairs. A least-squares solution is found by iteratively adjusting the vector difference between hypocentral pairs. The double-difference algorithm minimizes errors due to unmodeled velocity structure without the use of station corrections. Because catalog and cross-correlation data are combined into one system of equations, interevent distances within multiplets are determined to the accuracy of the cross-correlation data, while the relative locations between multiplets and uncorrelated events are simultaneously determined to the accuracy of the absolute travel-time data. Statistical resampling methods are used to estimate data accuracy and location errors. Uncertainties in double-difference locations are improved by more than an order of magnitude compared to catalog locations. The algorithm is tested, and its performance is demonstrated on two clusters of earthquakes located on the northern Hayward fault, California. There it colapses the diffuse catalog locations into sharp images of seismicity and reveals horizontal lineations of hypocenter that define the narrow regions on the fault where stress is released by brittle failure.

Bulletin of the Seismological Society of America↗
USGS1999

Slip-parallel seismic lineations on the Northern Hayward Fault, California

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

Geophysical Research Letters↗
USGS1998

Observation of the seismic nucleation phase in the Ridgecrest, California, earthquake sequence

Near-source observations of five M 3.8–5.2 earthquakes near Ridgecrest, California are consistent with the presence of a seismic nucleation phase. These earthquakes start abruptly, but then slow or stop before rapidly growing again toward their maximum rate of moment release. Deconvolution of instrument and path effects by empirical Green's functions demonstrates that the initial complexity at the start of the earthquake is a source effect. The rapid growth of the P -wave arrival at the start of the seismic nucleation phase supports the conclusion of Mori and Kanamori [1996] that these earthquakes begin without a magnitude-scaled slow initial phase of the type observed by Iio [1992, 1995].

Geophysical Research Letters↗
USGS1996

Detailed observations of California foreshock sequences: Implications for the earthquake initiation process

We find that foreshocks provide clear evidence for an extended nucleation process before some earthquakes. In this study, we examine in detail the evolution of six California foreshock sequences, the 1986 Mount Lewis ( M L = 5.5), the 1986 Chalfant ( M L = 6.4), the 1986 Stone Canyon ( M L = 4.7), the 1990 Upland ( M L = 5.2), the 1992 Joshua Tree ( M W = 6.1), and the 1992 Landers ( M W = 7.3) sequence. Typically, uncertainties in hypocentral parameters are too large to establish the geometry of foreshock sequences and hence to understand their evolution. However, the similarity of location and focal mechanisms for the events in these sequences leads to similar foreshock waveforms that we cross correlate to obtain extremely accurate relative locations. We use these results to identify small-scale fault zone structures that could influence nucleation and to determine the stress evolution leading up to the mainshock. In general, these foreshock sequences are not compatible with a cascading failure nucleation model in which the foreshocks all occur on a single fault plane and trigger the mainshock by static stress transfer. Instead, the foreshocks seem to concentrate near structural discontinuities in the fault and may themselves be a product of an aseismic nucleation process. Fault zone heterogeneity may also be important in controlling the number of foreshocks, i.e., the stronger the heterogeneity, the greater the number of foreshocks. The size of the nucleation region, as measured by the extent of the foreshock sequence, appears to scale with mainshock moment in the same manner as determined independently by measurements of the seismic nucleation phase. We also find evidence for slip localization as predicted by some models of earthquake nucleation.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Properties of the seismic nucleation phase

Near-source observations show that earthquakes begin abruptly at the P-wave arrival, but that this beginning is weak, with a low moment rate relative to the rest of the main shock. We term this initial phase of low moment rate the seismic nucleation phase. We have observed the seismic nucleation phase for a set of 48 earthquakes ranging in magnitude from 1.1–8.1 The size and duration of the seismic nucleation phase scale with the total seismic moment of the earthquake, suggesting that the process responsible for the seismic nucleation phase carries information about the eventual size of the earthquake. The seismic nucleation phase is characteristically followed by quadratic growth in the moment rate, consistent with self-similar rupture at constant stress drop. In this paper we quantify the properties of the seismic nucleation phase and offer several possible explanations for it.

California↗
USGS1995

Seismic evidence for an earthquake nucleation phase

Near-source observations show that earthquakes initiate with a distinctive seismic nucleation phase that is characterized by a low rate of moment release relative to the rest of the event. This phase was observed for the 30 earthquakes having moment magnitudes 2.6 to 8.1, and the size and duration of this phase scale with the eventual size of the earthquake. During the nucleation phase, moment release was irregular and appears to have been confined to a limited region of the fault. It was characteristically followed by quadratic growth in the moment rate as rupture began to propagate away from the nucleation zone. These observations suggest that the nucleation process exerts a strong influence on the size of the eventual earthquake.

Science↗