Geology Reports⌕ Search

Geology topics

S. E. Hough

Publications and source records attributed to S. E. Hough.

At least 37 records · Page 2Linked to original sources
USGS2002

Source properties of earthquakes near the Salton Sea triggered by the 16 October 1999 M 7.1 Hector Mine, California, earthquake

We analyze the source properties of a sequence of triggered earthquakes that occurred near the Salton Sea in southern California in the immediate aftermath of the M 7.1 Hector Mine earthquake of 16 October 1999. The sequence produced a number of early events that were not initially located by the regional network, including two moderate earthquakes: the first within 30 sec of the P-wave arrival and a second approximately 10 minutes after the mainshock. We use available amplitude and waveform data from these events to estimate magnitudes to be approximately 4.7 and 4.4, respectively, and to obtain crude estimates of their locations. The sequence of small events following the initial M 4.7 earthquake is clustered and suggestive of a local aftershock sequence. Using both broadband TriNet data and analog data from the Southern California Seismic Network (SCSN), we also investigate the spectral characteristics of the M 4.4 event and other triggered earthquakes using empirical Green's function (EGF) analysis. We find that the source spectra of the events are consistent with expectations for tectonic (brittle shear failure) earthquakes, and infer stress drop values of 0.1 to 6 MPa for six M 2.1 to M 4.4 events. The estimated stress drop values are within the range observed for tectonic earthquakes elsewhere. They are relatively low compared to typically observed stress drop values, which is consistent with expectations for faulting in an extensional, high heat flow regime. The results therefore suggest that, at least in this case, triggered earthquakes are associated with a brittle shear failure mechanism. This further suggests that triggered earthquakes may tend to occur in geothermal-volcanic regions because shear failure occurs at, and can be triggered by, relatively low stresses in extensional regimes.

Bulletin of the Seismological Society of America↗
USGS2002

Magnitude estimates of two large aftershocks of the 16 December 1811 New Madrid earthquake

The three principal New Madrid mainshocks of 1811-1812 were followed by extensive aftershock sequences that included numerous felt events. Although no instrumental data are available for either the mainshocks or the aftershocks, available historical accounts do provide information that can be used to estimate magnitudes and locations for the large events. In this article we investigate two of the largest aftershocks: one near dawn following the first mainshock on 16 December 1811, and one near midday on 17 December 1811. We reinterpret original felt reports to obtain a set of 48 and 20 modified Mercalli intensity values of the two aftershocks, respectively. For the dawn aftershock, we infer a Mw of approximately 7.0 based on a comparison of its intensities with those of the smallest New Madrid mainshock. Based on a detailed account that appears to describe near-field ground motions, we further propose a new fault rupture scenario for the dawn aftershock. We suggest that the aftershock had a thrust mechanism and occurred on a southeastern limb of the Reelfoot fault. For the 17 December 1811 aftershock, we infer a M W of approximately 6.1 ± 0.2 . This value is determined using the method of Bakun et al. (2002), which is based on a new calibration of intensity versus distance for earthquakes in central and eastern North America. The location of this event is not well constrained, but the available accounts suggest an epicenter beyond the southern end of the New Madrid Seismic Zone.

Arkansas, Illinois, Missouri, Tennessee↗
USGS2002

The 26 January 2001 M 7.6 Bhuj, India, earthquake: Observed and predicted ground motions

Although local and regional instrumental recordings of the devastating 26, January 2001, Bhuj earthquake are sparse, the distribution of macroseismic effects can provide important constraints on the mainshock ground motions. We compiled available news accounts describing damage and other effects and interpreted them to obtain modified Mercalli intensities (MMIs) at >200 locations throughout the Indian subcontinent. These values are then used to map the intensity distribution throughout the subcontinent using a simple mathematical interpolation method. Although preliminary, the maps reveal several interesting features. Within the Kachchh region, the most heavily damaged villages are concentrated toward the western edge of the inferred fault, consistent with western directivity. Significant sediment-induced amplification is also suggested at a number of locations around the Gulf of Kachchh to the south of the epicenter. Away from the Kachchh region, intensities were clearly amplified significantly in areas that are along rivers, within deltas, or on coastal alluvium, such as mudflats and salt pans. In addition, we use fault-rupture parameters inferred from teleseismic data to predict shaking intensity at distances of 0-1000 km. We then convert the predicted hard-rock ground-motion parameters to MMI by using a relationship (derived from Internet-based intensity surveys) that assigns MMI based on the average effects in a region. The predicted MMIs are typically lower by 1-3 units than those estimated from news accounts, although they do predict near-field ground motions of approximately 80%g and potentially damaging ground motions on hard-rock sites to distances of approximately 300 km. For the most part, this discrepancy is consistent with the expected effect of sediment response, but it could also reflect other factors, such as unusually high building vulnerability in the Bhuj region and a tendency for media accounts to focus on the most dramatic damage, rather than the average effects. The discrepancy may also be partly attributable to the inadequacy of the empirical relationship between MMI and peak ground acceleration (PGA), when applied to India. The MMI-PGA relationship was developed using data from California earthquakes, which might have a systematically different stress drop and therefore, a different frequency content than intraplate events. When a relationship between response spectra and MMI is used, we obtain larger predicted MMI values, in better agreement with the observations.

Bulletin of the Seismological Society of America↗
USGS2001

Empirical Green's function analysis of recent moderate events in California

I use seismic data from portable digital stations and the broadband Terrascope network in southern California to investigate radiated earthquake source spectra and discuss the results in light of previous studies on both static stress drop and apparent stress. Applying the empirical Green's function (EGF) method to two sets of M 4-6.1 events, I obtain deconvolved source-spectra estimates and corner frequencies. The results are consistent with an ??2 source model and constant Brune stress drop. However, consideration of the raw spectral shapes of the largest events provides evidence for a high-frequency decay more shallow than ??2. The intermediate (???f-1) slope cannot be explained plausibly with attenuation or site effects and is qualitatively consistent with a model incorporating directivity effects and a fractional stress-drop rupture process, as suggested by Haddon (1996). However, the results obtained in this study are not consistent with the model of Haddon (1996) in that the intermediate slope is not revealed with EGF analysis. This could reflect either bandwidth limitations inherent in EGF analysis or perhaps a rupture process that is not self-similar. I show that a model with an intermediate spectral decay can also reconcile the apparent discrepancy between the scaling of static stress drop and that of apparent stress drop for moderate-to-large events.

Bulletin of the Seismological Society of America↗
USGS2001

Triggered earthquakes and the 1811-1812 New Madrid, central United States, earthquake sequence

The 1811-1812 New Madrid, central United States, earthquake sequence included at least three events with magnitudes estimated at well above M 7.0. I discuss evidence that the sequence also produced at least three substantial triggered events well outside the New Madrid Seismic Zone, most likely in the vicinity of Cincinnati, Ohio. The largest of these events is estimated to have a magnitude in the low to mid M 5 range. Events of this size are large enough to cause damage, especially in regions with low levels of preparedness. Remotely triggered earthquakes have been observed in tectonically active regions in recent years, but not previously in stable continental regions. The results of this study suggest, however, that potentially damaging triggered earthquakes may be common following large mainshocks in stable continental regions. Thus, in areas of low seismic activity such as central/ eastern North America, the hazard associated with localized source zones might be more far reaching than previously recognized. The results also provide additional evidence that intraplate crust is critically stressed, such that small stress changes are especially effective at triggering earthquakes.

Bulletin of the Seismological Society of America↗
USGS2000

On the modified Mercalli intensities and magnitudes of the 1811-1812 New Madrid earthquakes

We reexamine original felt reports from the 1811–1812 New Madrid earthquakes and determine revised isoseismal maps for the three principal mainshocks. In many cases we interpret lower values than those assigned by earlier studies. In some cases the revisions result from an interpretation of original felt reports with an appreciation for site response issues. Additionally, earlier studies had assigned modified Mercalli intensity (MMI) values of V-VII to a substantial number of reports that we conclude do not describe damage commensurate with intensities this high. We investigate several approaches to contouring the MMI values using both analytical and subjective methods. For the first mainshock on 02∶15 LT December 16, 1811, our preferred contouring yields M ω7.2–7.3 using the area-moment regressions of Johnston [1996]. For the 08∶00 LT on January 23, 1812, and 03∶45 LT on February 7, 1812, mainshocks, we obtain M ω7.0 and M ω7.4–7.5, respectively. Our magnitude for the February mainshock is consistent with the established geometry of the Reelfoot fault, which all evidence suggests to have been the causative structure for this event. We note that the inference of lower magnitudes for the New Madrid events implies that site response plays a significant role in controlling seismic hazard at alluvial sites in the central and eastern United States. We also note that our results suggest that thrusting may have been the dominant mechanism of faulting associated with the 1811–1812 sequence.

Arkansas, Illinois, Kentucky, Missouri, Tennessee↗
USGS2000

The 1998 earthquake sequence south of Long Valley Caldera, California: Hints of magmatic involvement

A significant episode of seismic and geodetic unrest took place at Long Valley Caldera, California, beginning in the summer of 1997. Activity through late May of 1998 was concentrated in and around the south moat and the south margin of the resurgent dome. The Sierran Nevada block (SNB) region to the south/southeast remained relatively quiet until a M 5.1 event occurred there on 9 June 1998 (UT). A second M 5.1 event followed on 15 July (UT); both events were followed by appreciable aftershock sequences. An additional, distinct burst of activity began on 1 August 1998. The number of events in the August sequence (over the first week or two) was similar to the aftershock sequence of the 15 July 1998 M 5.1 event, but the later sequence was not associated with any events larger than M 4.3. All of the summer 1998 SNB activity was considered tectonic rather than magmatic; in general the SNB is considered an unlikely location for future eruptions. However, the August sequence-an 'aftershock sequence without a mainshock'- is suggestive of a strain event larger than the cumulative seismotectonic strain release. Moreover, a careful examination of waveforms from the August sequence reveals a small handful of events whose spectral signature is strikingly harmonic. We investigate the waveforms of these events using spectral, autocorrelation, and empirical Green's function techniques and conclude that they were most likely associated with a fluid-controlled source. Our observations suggest that there may have been some degree of magma or magma-derived fluid involvement in the 1998 SNB sequence.

California↗
USGS1999

Attenuation and source properties at the Coso Geothermal area, California

We use a multiple-empirical Green's function method to determine source properties of small ( M −0.4 to 1.3) earthquakes and P - and S -wave attenuation at the Coso Geothermal Field, California. Source properties of a previously identified set of clustered events from the Coso geothermal region are first analyzed using an empirical Green's function ( EGF ) method. Stress-drop values of at least 0.5-1 MPa are inferred for all of the events; in many cases, the corner frequency is outside the usable bandwidth, and the stress drop can only be constrained as being higher than 3 MPa. P - and S -wave stress-drop estimates are identical to the resolution limits of the data. These results are indistinguishable from numerous EGF studies of M 2-5 earthquakes, suggesting a similarity in rupture processes that extends to events that are both tiny and induced, providing further support for Byerlee's Law. Whole-path Q estimates for P and S waves are determined using the multiple-empirical Green's function ( MEGF ) method of Hough (1997), whereby spectra from clusters of colocated events at a given station are inverted for a single attenuation parameter, κ, with source parameters constrained from EGF analysis. The κ estimates, which we infer to be resolved to within 0.01 sec or better, exhibit almost as much scatter as a function of hypocentral distance as do values from previous single-spectrum studies for which much higher uncertainties in individual κ estimates are expected. The variability in κ estimates determined here therefore suggests real lateral variability in Q structure. Although the ray-path coverage is too sparse to yield a complete three-dimensional attenuation tomographic image, we invert the inferred κ value for three-dimensional structure using a damped least-squares method, and the results do reveal significant lateral variability in Q structure. The inferred attenuation variability corresponds to the heat-flow variations within the geothermal region. A central low- Q region corresponds well with the central high-heat flow region; additional detailed structure is also suggested.

California↗
USGS1997

Empirical Green's function analysis: Taking the next step

An extension of the empirical Green's function (EGF) method is presented that involves determination of source parameters using standard EGF deconvolution, followed by inversion for a common attenuation parameter for a set of colocated events. Recordings of three or more colocated events can thus be used to constrain a single path attenuation estimate. I apply this method to recordings from the 1995–1996 Ridgecrest, California, earthquake sequence; I analyze four clusters consisting of 13 total events with magnitudes between 2.6 and 4.9. I first obtain corner frequencies, which are used to infer Brune stress drop estimates. I obtain stress drop values of 0.3–53 MPa (with all but one between 0.3 and 11 MPa), with no resolved increase of stress drop with moment. With the corner frequencies constrained, the inferred attenuation parameters are very consistent; they imply an average shear wave quality factor of approximately 20–25 for alluvial sediments within the Indian Wells Valley. Although the resultant spectral fitting (using corner frequency and κ) is good, the residuals are consistent among the clusters analyzed. Their spectral shape is similar to the the theoretical one-dimensional response of a layered low-velocity structure in the valley (an absolute site response cannot be determined by this method, because of an ambiguity between absolute response and source spectral amplitudes). I show that even this subtle site response can significantly bias estimates of corner frequency and κ, if it is ignored in an inversion for only source and path effects. The multiple-EGF method presented in this paper is analogous to a joint inversion for source, path, and site effects; the use of colocated sets of earthquakes appears to offer significant advantages in improving resolution of all three estimates, especially if data are from a single site or sites with similar site response.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

On the coherence of ground motion in the San Fernando valley

We present an analysis of the coherence of seismic ground motion recorded on alluvial sediments in the San Fernando Valley, California. Using aftershocks of the 17 January 1994 M w 6.7 earthquake recorded at a quasi-dense array of portable stations, we analyze the coherence of three well-recorded magnitude 3.7 to 4.0 events over the frequency range 0.5 to 15 Hz and a distance range of 0.5 to 5.3 km. All stations are located at sites with broadly similar near-site geology, characterized by medium to fine-grain Quaternary alluvial sediments. On average, relatively high values of coherence are observed for distances up to 3 to 4 km and frequencies up to 2 to 3 Hz; coherence drops sharply at frequencies near and above 3 Hz. Although average coherence functions are described reasonably well by a log-linear relationship with frequency, the curves at all distances exhibit a flattening at low frequencies that is not consistent with previous observations of coherence at hardrock sites. The distance decay of coherence is also markedly less strong, with high coherence values observed over station separations corresponding to multiple wavelengths. This may reflect fundamental differences in shallow-wave propagation in the two environments, with high-frequency scattering relatively more dominant in regions of hard-rock near-surface geology. Within a sedimentary basin or valley, the site response itself generally reflects a resonance phenomenon that may tend to give rise to more uniform ground motions. However, previous studies have demonstrated the existence of pathological focusing and amplification effects within complex sedimentary basin environments such as the greater Los Angeles region; our results undoubtedly do not quantify the full range of ground-motion variability at all sites, but rather represent the level of that variability that can be expected, and quantified, for typical source/receiver paths.

California↗
USGS1996

Observational constraints on earthquake source scaling: Understanding the limits in resolution

I examine the resolution of the type of stress drop estimates that have been used to place observational constraints on the scaling of earthquake source processes. I first show that apparent stress and Brune stress drop are equivalent to within a constant given any source spectral decay between ω 1.5 and ω 3 (i.e., any plausible value) and so consistent scaling is expected for the two estimates. I then discuss the resolution and scaling of Brune stress drop estimates, in the context of empirical Green's function results from recent earthquake sequences, including the 1992 Joshua Tree, California, mainshock and its aftershocks. I show that no definitive scaling of stress drop with moment is revealed over the moment range 10 19 –10 25 ; within this sequence, however, there is a tendency for moderate-sized ( M 4–5) events to be characterized by high stress drops. However, well-resolved results for recent M > 6 events are inconsistent with any extrapolated stress increase with moment for the aftershocks. Focusing on corner frequency estimates for smaller ( M < 3.5) events, I show that resolution is extremely limited even after empirical Green's function deconvolutions. A fundamental limitation to resolution is the paucity of good signal-to-noise at frequencies above 60 Hz, a limitation that will affect nearly all surficial recordings of ground motion in California and many other regions. Thus, while the best available observational results support a constant stress drop for moderate- to large-sized events, very little robust observational evidence exists to constrain the quantities that bear most critically on our understanding of source processes: stress drop values and stress drop scaling for small events.

California↗
USGS1995

Earthquakes in the Los Angeles metropolitan region: A possible fractal distribution of rupture size

Although there is debate on the maximum size of earthquake that is possible on any of several known fault systems in the greater Los Angeles metropolitan region, it is reasonable to assume that the distribution of earthquakes will follow a fractal distribution of rupture areas. For this assumption and an overall slip-rate for the region of approximately 1 centimeter per year, roughly one magnitude 7.4 to 7.5 event is expected to occur every 245 to 325 years. A model in which the earthquake distribution is fractal predicts that, additionally, there should be approximately six events in the range of magnitude 6.6 in this same span of time, a higher rate than has occurred in the historic record.

California↗
USGS1994

Southern surface rupture associated with the M 7.3 1992 Landers, California, earthquake

Although most evidence suggests that the 28 June 1992 M 7.3 Landers earthquake ruptured unilaterally north, significant surface rupture was mapped on the Eureka Peak and Burnt Mountain faults, to the south of the Landers epicenter. An eyewitness account reports that surface rupture occurred on the northern Eureka Peak fault within approximately 35 sec of the mainshock initiation. Array analysis of the Landers mainshock provides evidence in support of this report; a significant southern subevent in the early mainshock coda. I also analyze dense array recordings of a M 5.6 aftershock that occurred 3 min after the mainshock at 34°7.65′N, 116°23.82′W and show that there is strong evidence that this event was also associated with significant rupture on the Eureka Peak fault. This analysis thus suggests that the Eureka Peak fault rupture was not caused by direct bilateral mainshock rupture but instead was initially triggered less than a minute after the mainshock and reruptured by the M 5.6 aftershock. Results for the evolution of the Landers sequence suggest that mainshock subevents may in some cases be accurately described as aftershocks (i.e., disjoint triggered events) that occur within the duration of mainshock strong ground motion.

California↗
USGS1994

Fault-zone waves observed at the southern Joshua Tree earthquake rupture zone

Waveform and spectral characteristics of several aftershocks of the M 6.1 22 April 1992 Joshua Tree earthquake recorded at stations just north of the Indio Hills in the Coachella Valley can be interpreted in terms of waves propagating within narrow, low-velocity, high-attenuation, vertical zones. Evidence for our interpretation consists of: (1) emergent P arrivals prior to and opposite in polarity to the impulsive direct phase; these arrivals can be modeled as headwaves indicative of a transfault velocity contrast; (2) spectral peaks in the S wave train that can be interpreted as internally reflected, low-velocity fault-zone wave energy; and (3) spatial selectivity of event-station pairs at which these data are observed, suggesting a long, narrow geologic structure. The observed waveforms are modeled using the analytical solution of Ben-Zion and Aki (1990) for a plane-parallel layered fault-zone structure. Synthetic waveform fits to the observed data indicate the presence of NS-trending vertical fault-zone layers characterized by a thickness of 50 to 100 m, a velocity decrease of 10 to 15% relative to the surrounding rock, and a P -wave quality factor in the range 25 to 50.

Bulletin of the Seismological Society of America↗