USGS Science⌕ Search

SEARCH · USGS Science

Results for “Article”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47Linked to original sources
USGS2006

Surface waves in the western Taiwan coastal plain from an aftershock of the 1999 Chi-Chi, Taiwan, earthquake

Significant surface waves were recorded in the western coastal plain (WCP) of Taiwan during the 1999 Chi-Chi, Taiwan, earthquake and its series of aftershocks. We study in detail the surface waves produced by one aftershock (20 September 1999, 18hr 03m 41.16sec, M 6.2) in this paper. We take the Chelungpu-Chukou fault to be the eastern edge of the WCP because it marks a distinct lateral contrast in seismic wave velocities in the upper few kilometers of the surface. For many records from stations within the WCP, body waves and surface waves separate well in both the time domain and the period domain. Long-period (e.g., >2 sec) ground motions in the plain are dominated by surface waves. Significant prograde Rayleigh wave particle motions were observed in the WCP. The observed peak ground velocities are about 3-5 times larger than standard predictions in the central and western part of the plain. Observed response spectra at 3 sec, 4 sec, and 5 sec at the center of the plain can be 15 times larger than standard predictions and 10 times larger than the predictions of Joyner (2000) based on surface wave data from the Los Angeles basin. The strong surface waves were probably generated at the boundary of the WCP and then propagated toward the west, largely along radial directions relative to the epicenter. The geometry of the boundary may have had a slight effect on propagation directions of surface waves. Group velocities of fundamental mode Rayleigh and Love waves are estimated using the multiple filter analysis (MFA) technique and are refined with phase matched filtering (PMF). Group velocities of fundamental mode surface waves range from about 0.7 km/sec to 1.5 km/sec for the phases at periods from 3 sec to 10 sec. One important observation from this study is that the strongest surface waves were recorded in the center of the plain. The specific location of the strongest motions depends largely on the period of surface waves rather than on specific site conditions or plain structures. Accordingly, we conjecture that surface waves could be generated in a wide area close to boundaries of low-velocity sedimentary wave guides. In the case studied in this article the area can be as wide as 30 km (from the Chelungpu fault to the center of the plain). Surface waves converted by P and S waves at different locations would overlap each other and add constructively along their propagation paths. As a result, the surface waves would get stronger and stronger. Beyond a certain distance to the boundary, no more surface waves would be generated. Consequently, no more local surface waves would be superimposed into the invasive surface waves, and the surface waves would tend to decay in amplitude with distance.

Bulletin of the Seismological Society of America↗
USGS2006

Uncertainty of earthquake losses due to model uncertainty of input ground motions in the Los Angeles area

In a recent study we used the Monte Carlo simulation method to evaluate the ground-motion uncertainty of the 2002 update of the California probabilistic seismic hazard model. The resulting ground-motion distribution is used in this article to evaluate the contribution of the hazard model to the uncertainty in earthquake loss ratio, the ratio of the expected loss to the total value of a structure. We use the Hazards U.S. (HAZUS) methodology for loss estimation because it is a widely used and publicly available risk model and intended for regional studies by public agencies and for use by governmental decision makers. We found that the loss ratio uncertainty depends not only on the ground-motion uncertainty but also on the mean ground-motion level. The ground-motion uncertainty, as measured by the coefficient of variation (COV), is amplified when converting to the loss ratio uncertainty because loss increases concavely with ground motion. By comparing the ground-motion uncertainty with the corresponding loss ratio uncertainty for the structural damage of light wood-frame buildings in Los Angeles area, we show that the COV of loss ratio is almost twice the COV of ground motion with a return period of 475 years around the San Andreas fault and other major faults in the area. The loss ratio for the 2475-year ground-motion maps is about a factor of three higher than for the 475-year maps. However, the uncertainties in ground motion and loss ratio for the longer return periods are lower than for the shorter return periods because the uncertainty parameters in the hazard logic tree are independent of the return period, but the mean ground motion increases with return period.

Bulletin of the Seismological Society of America↗
USGS2006

Orientation-independent measures of ground motion

The geometric mean of the response spectra for two orthogonal horizontal components of motion, commonly used as the response variable in predictions of strong ground motion, depends on the orientation of the sensors as installed in the field. This means that the measure of ground-motion intensity could differ for the same actual ground motion. This dependence on sensor orientation is most pronounced for strongly correlated motion (the extreme example being linearly polarized motion), such as often occurs at periods of 1 sec or longer. We propose two new measures of the geometric mean, GMRotDpp, and GMRotIpp, that are independent of the sensor orientations. Both are based on a set of geometric means computed from the as-recorded orthogonal horizontal motions rotated through all possible non-redundant rotation angles. GMRotDpp is determined as the ppth percentile of the set of geometric means for a given oscillator period. For example, GMRotDOO, GMRotD50, and GMRotD100 correspond to the minimum, median, and maximum values, respectively. The rotations that lead to GMRotDpp depend on period, whereas a single-period-independent rotation is used for GMRotIpp, the angle being chosen to minimize the spread of the rotation-dependent geometric mean (normalized by GMRotDpp) over the usable range of oscillator periods. GMRotI50 is the ground-motion intensity measure being used in the development of new ground-motion prediction equations by the Pacific Earthquake Engineering Center Next Generation Attenuation project. Comparisons with as-recorded geometric means for a large dataset show that the new measures are systematically larger than the geometric-mean response spectra using the as-recorded values of ground acceleration, but only by a small amount (less than 3%). The theoretical advantage of the new measures is that they remove sensor orientation as a contributor to aleatory uncertainty. Whether the reduction is of practical significance awaits detailed studies of large datasets. A preliminary analysis contained in a companion article by Beyer and Bommer finds that the reduction is small-to-nonexistent for equations based on a wide range of magnitudes and distances. The results of Beyer and Bommer do suggest, however, that there is an increasing reduction as period increases. Whether the reduction increases with other subdivisions of the dataset for which strongly correlated motions might be expected (e.g., pulselike motions close to faults) awaits further analysis.

Bulletin of the Seismological Society of America↗
USGS2008

High-resolution seismic images and seismic velocities of the San Andreas fault zone at Burro Flats, Southern California

To better understand the structure of the San Andreas fault (SAF) at Burro Flats in southern California, we acquired a three-dimensional combined set of seismic reflection and refraction profiles centered on the main active trace at Burro Flats. In this article, we discuss the variation in shallow-depth velocities along each seismic profile, with special emphasis on the 1500 m/sec P -wave velocity contour, which can be an indicator of shallow-depth water-saturated unconsolidated sediments. Along the four seismic profiles, minimum depths of the groundwater table, as inferred from 1500 m/sec velocity contour, range from 10 to about 20 m. The largest variations in depth to the top of the groundwater table occur in areas near mapped faults, suggesting that the groundwater flow in Burro Flats is strongly affected by the locations of fault traces. We also used the seismic data to develop seismic reflection images that show multiple strands of the SAF in the upper 60 m. Reflectors above the 10 m depth probably correspond to Holocene alluvial deposits; reflectors below the 15 m depth probably arise from velocity or density variations within the Precambrian gneiss complex, likely due to weathering. Apparent vertical offsets of reflectors are observed along profiles (lines 1 and 2) that are normal to the SAF, indicating minor apparent vertical offsets on the SAF at shallow depths. Along line 2, the apparently vertically offset reflectors correlate with zones of relatively low P -wave velocity. Along the central part of lines 1 and 2, the faults form a flower structure, which is typical of strike-slip faults such as the SAF.

California↗
USGS2009

Suggested notation conventions for rotational seismology

We note substantial inconsistency among authors discussing rotational motions observed with inertial seismic sensors (and much more so in the broader topic of rotational phenomena). Working from physics and other precedents, we propose standard terminology and a preferred reference frame for inertial sensors (Fig. 1) that may be consistently used in discussions of both finite and infinitesimal observed rotational and translational motions in seismology and earthquake engineering. The scope of this article is limited to observations because there are significant differences in the analysis of finite and infinitesimal rotations, though such discussions should remain compatible with those presented here where possible. We recommend the general use of the notation conventions presented in this tutorial, and we recommend that any deviations or alternatives be explicitly defined.

Bulletin of the Seismological Society of America↗
USGS2009

Laboratory and field testing of commercial rotational seismometers

There are a small number of commercially available sensors to measure rotational motion in the frequency and amplitude ranges appropriate for earthquake motions on the ground and in structures. However, the performance of these rotational seismometers has not been rigorously and independently tested and characterized for earthquake monitoring purposes as is done for translational strong- and weak-motion seismometers. Quantities such as sensitivity, frequency response, resolution, and linearity are needed for the understanding of recorded rotational data. To address this need, we, with assistance from colleagues in the United States and Taiwan, have been developing performance test methodologies and equipment for rotational seismometers. In this article the performance testing methodologies are applied to samples of a commonly used commercial rotational seismometer, the eentec model R-1. Several examples were obtained for various test sequences in 2006, 2007, and 2008. Performance testing of these sensors consisted of measuring: (1) sensitivity and frequency response; (2) clip level; (3) self noise and resolution; and (4) cross-axis sensitivity, both rotational and translational. These sensor-specific results will assist in understanding the performance envelope of the R-1 rotational seismometer, and the test methodologies can be applied to other rotational seismometers.

Bulletin of the Seismological Society of America↗
USGS2009

Major earthquakes recorded by Speleothems in Midwestern U.S. caves

Historic earthquakes generated by the New Madrid seismic zone represent some of the largest recorded in the United States, yet prehistoric events are recognized only through deformation in late-Wisconsin to Holocene-age, near surface sediments (liquefaction, monoclinal folding, and changes in river meanders). In this article, we show that speleothems in caves of southwestern Illinois and southeastern Missouri may constitute a previously unrecognized recorder of large earthquakes in the U.S. midcontinent region. The timing of the initiation and regrowth of stalagmites in southwestern Illinois and southeastern Missouri caves is consistent with the historic and prehistoric record of several known seismic events in the U.S. midcontinent region. We conclude that dating the initiation of original stalagmite growth and later postearthquake rejuvenation constitutes a new paleoseismic method that has the potential for being applied to any region around the world in the vicinity of major seismic zones where caves exist. Use of this technique could expand the geographical distribution of paleoseimic data, document prehistoric earthquakes, and help improve interpretations of paleoearthquakes.

Bulletin of the Seismological Society of America↗
USGS2010

Tsunami simulations of the 1867 Virgin Island earthquake: Constraints on epicenter location and fault parameters

The 18 November 1867 Virgin Island earthquake and the tsunami that closely followed caused considerable loss of life and damage in several places in the northeast Caribbean region. The earthquake was likely a manifestation of the complex tectonic deformation of the Anegada Passage, which cuts across the Antilles island arc between the Virgin Islands and the Lesser Antilles. In this article, we attempt to characterize the 1867 earthquake with respect to fault orientation, rake, dip, fault dimensions, and first tsunami wave propagating phase, using tsunami simulations that employ high-resolution multibeam bathymetry. In addition, we present new geophysical and geological observations from the region of the suggested earthquake source. Results of our tsunami simulations based on relative amplitude comparison limit the earthquake source to be along the northern wall of the Virgin Islands basin, as suggested by Reid and Taber (1920) , or on the carbonate platform north of the basin, and not in the Virgin Islands basin, as commonly assumed. The numerical simulations suggest the 1867 fault was striking 120°–135° and had a mixed normal and left-lateral motion. First propagating wave phase analysis suggests a fault striking 300°–315° is also possible. The best-fitting rupture length was found to be relatively small (50 km), probably indicating the earthquake had a moment magnitude of ∼7.2. Detailed multibeam echo sounder surveys of the Anegada Passage bathymetry between St. Croix and St. Thomas reveal a scarp, which cuts the northern wall of the Virgin Islands basin. High-resolution seismic profiles further indicate it to be a reasonable fault candidate. However, the fault orientation and the orientation of other subparallel faults in the area are more compatible with right-lateral motion. For the other possible source region, no clear disruption in the bathymetry or seismic profiles was found on the carbonate platform north of the basin.

Bulletin of the Seismological Society of America↗
USGS2010

Reply to “Comment on ‘Is There a Basis for Preferring Characteristic Earthquakes over a Gutenberg–Richter Distribution in Probabilistic Earthquake Forecasting?’ by Tom Parsons and Eric L. Geist” by Jens-Uwe Klügel

The focus of Parsons and Geist (2009) was to test whether the key observational data used in earthquake forecasting necessitate a characteristic earthquake rupture model. The point of our article was not to suggest that a specific form of the Gutenberg–Richter earthquake distribution is a perfect representation of reality. The uncertainties surrounding event slip estimates, paleoseismic event rates, and observed a and b values in catalog magnitude–frequency distributions are broad. So broad, in fact, that giving full weight to just one model of earthquake rupture behavior in formal forecasting is unjustified. Further, the characteristic earthquake model requires definition of rupture segments, which introduces a series of unquantifiable uncertainties that are seldom addressed in forecasts (e.g., Field et al., 2009).

Bulletin of the Seismological Society of America↗
USGS2015

Revisions to some parameters used in stochastic-method simulations of ground motion

The stochastic method of ground‐motion simulation specifies the amplitude spectrum as a function of magnitude ( M ) and distance ( R ). The manner in which the amplitude spectrum varies with M and R depends on physical‐based parameters that are often constrained by recorded motions for a particular region (e.g., stress parameter, geometrical spreading, quality factor, and crustal amplifications), which we refer to as the seismological model. The remaining ingredient for the stochastic method is the ground‐motion duration. Although the duration obviously affects the character of the ground motion in the time domain, it also significantly affects the response of a single‐degree‐of‐freedom oscillator. Recently published updates to the stochastic method include a new generalized double‐corner‐frequency source model, a new finite‐fault correction, a new parameterization of duration, and a new duration model for active crustal regions. In this article, we augment these updates with a new crustal amplification model and a new duration model for stable continental regions. Random‐vibration theory (RVT) provides a computationally efficient method to compute the peak oscillator response directly from the ground‐motion amplitude spectrum and duration. Because the correction factor used to account for the nonstationarity of the ground motion depends on the ground‐motion amplitude spectrum and duration, we also present new RVT correction factors for both active and stable regions.

Bulletin of the Seismological Society of America↗
USGS2015

A century of induced earthquakes in Oklahoma?

Seismicity rates have increased sharply since 2009 in the central and eastern United States, with especially high rates of activity in the state of Oklahoma. Growing evidence indicates that many of these events are induced, primarily by injection of wastewater in deep disposal wells. The upsurge in activity has raised two questions: What is the background rate of tectonic earthquakes in Oklahoma? How much has the rate varied throughout historical and early instrumental times? In this article, we show that (1) seismicity rates since 2009 surpass previously observed rates throughout the twentieth century; (2) several lines of evidence suggest that most of the significant earthquakes in Oklahoma during the twentieth century were likely induced by oil production activities, as they exhibit statistically significant temporal and spatial correspondence with disposal wells, and intensity measurements for the 1952 El Reno earthquake and possibly the 1956 Tulsa County earthquake follow the pattern observed in other induced earthquakes; and (3) there is evidence for a low level of tectonic seismicity in southeastern Oklahoma associated with the Ouachita structural belt. The 22 October 1882 Choctaw Nation earthquake, for which we estimate Mw 4.8, occurred in this zone.

Oklahoma↗
USGS2016

Uncertainty in V s30 -based site response

Methods that account for site response range in complexity from simple linear categorical adjustment factors to sophisticated nonlinear constitutive models. Seismic‐hazard analysis usually relies on ground‐motion prediction equations (GMPEs); within this framework site response is modeled statistically with simplified site parameters that include the time‐averaged shear‐wave velocity to 30 m ( V S 30 ) and basin depth parameters. Because V S 30 is not known in most locations, it must be interpolated or inferred through secondary information such as geology or topography. In this article, we analyze a subset of stations for which V S 30 has been measured to address effects of V S 30 proxies on the uncertainty in the ground motions as modeled by GMPEs. The stations we analyze also include multiple recordings, which allow us to compute the repeatable site effects (or empirical amplification factors [EAFs]) from the ground motions. Although all methods exhibit similar bias, the proxy methods only reduce the ground‐motion standard deviations at long periods when compared to GMPEs without a site term, whereas measured V S 30 values reduce the standard deviations at all periods. The standard deviation of the ground motions are much lower when the EAFs are used, indicating that future refinements of the site term in GMPEs have the potential to substantially reduce the overall uncertainty in the prediction of ground motions by GMPEs.

Bulletin of the Seismological Society of America↗
USGS2017

Combining multiple earthquake models in real time for earthquake early warning

The ultimate goal of earthquake early warning (EEW) is to provide local shaking information to users before the strong shaking from an earthquake reaches their location. This is accomplished by operating one or more real‐time analyses that attempt to predict shaking intensity, often by estimating the earthquake’s location and magnitude and then predicting the ground motion from that point source. Other EEW algorithms use finite rupture models or may directly estimate ground motion without first solving for an earthquake source. EEW performance could be improved if the information from these diverse and independent prediction models could be combined into one unified, ground‐motion prediction. In this article, we set the forecast shaking at each location as the common ground to combine all these predictions and introduce a Bayesian approach to creating better ground‐motion predictions. We also describe how this methodology could be used to build a new generation of EEW systems that provide optimal decisions customized for each user based on the user’s individual false‐alarm tolerance and the time necessary for that user to react.

Bulletin of the Seismological Society of America↗
USGS2017

Determination of earthquake magnitude for early warning from the time-dependence of P-wave amplitudes

We propose a method that utilizes the time dependence of P ‐wave displacement amplitudes to estimate the final magnitude ( ⁠M⁠ ) for earthquake early warning (EEW) before the arrival of the peak amplitude. A relation between M and P ‐wave displacement amplitude is employed for the method. Its value is set as a function of time from the P arrival, and is determined using a K‐NET dataset recorded in Japan from a scaling relation between M and the time dependence of P ‐wave displacement. A test to check the performance of the proposed equation demonstrates in a statistical sense that this technique enables us to estimate M more rapidly than conventional methods without loss of accuracy. We conclude that the approach proposed in this article effectively gains a longer lead time as well as reduces the blind zone for EEW.

Bulletin of the Seismological Society of America↗
USGS2018

Spatial and spectral interpolation of ground-motion intensity measure observations

Following a significant earthquake, ground‐motion observations are available for a limited set of locations and intensity measures (IMs). Typically, however, it is desirable to know the ground motions for additional IMs and at locations where observations are unavailable. Various interpolation methods are available, but because IMs or their logarithms are normally distributed, spatially correlated, and correlated with each other at a given location, it is possible to apply the conditional multivariate normal (MVN) distribution to the problem of estimating unobserved IMs. In this article, we review the MVN and its application to general estimation problems, and then apply the MVN to the specific problem of ground‐motion IM interpolation. In particular, we present (1) a formulation of the MVN for the simultaneous interpolation of IMs across space and IM type (most commonly, spectral response at different oscillator periods) and (2) the inclusion of uncertain observation data in the MVN formulation. These techniques, in combination with modern empirical ground‐motion models and correlation functions, provide a flexible framework for estimating a variety of IMs at arbitrary locations.

Bulletin of the Seismological Society of America↗
USGS2017

Aftershocks, earthquake effects, and the location of the large 14 December 1872 earthquake near Entiat, central Washington

Reported aftershock durations, earthquake effects, and other observations from the large 14 December 1872 earthquake in central Washington are consistent with an epicenter near Entiat, Washington. Aftershocks were reported for more than 3 months only near Entiat. Modal intensity data described in this article are consistent with an Entiat area epicenter, where the largest modified Mercalli intensities, VIII, were assigned between Lake Chelan and Wenatchee. Although ground failures and water effects were widespread, there is a concentration of these features along the Columbia River and its tributaries in the Entiat area. Assuming linear ray paths, misfits from 23 reports of the directions of horizontal shaking have a local minima at Entiat, assuming the reports are describing surface waves, but the region having comparable misfit is large. Broadband seismograms recorded for comparable ray paths provide insight into the reasons why possible S – P times estimated from felt reports at two locations are several seconds too small to be consistent with an Entiat area epicenter.

Washington↗
USGS2018

The case for mean rupture distance in ground‐motion estimation

This article advocates for the use of mean rupture distances that we contend are more physically representative of the distance to an earthquake and are simpler than minimum distances. Many current ground‐motion models (GMMs) rely on numerous modifications of minimum rupture distances to accurately model near‐source ground motions. These modifications, that include additional distance definitions and saturation terms, result in complicated functional forms and are often not easily understood on a seismological basis, such as the magnitude‐dependent near‐fault saturation term. The use of mean distance represents the location of a station in relation to the entire rupture plane and results in a simpler, more physically meaningful GMM that models near‐source ground motion as accurately as other GMMs that have more inputs and more complex functional forms. We demonstrate the use of mean distance by developing a GMM for shallow‐crustal earthquakes with the Next Generation Attenuation‐West2 (NGA‐West2) project database. Specifically, we use the generalized mean distance, also known as the power mean, in which the power varies with frequency. We show that this new GMM fits the NGA‐West2 database as well as the NGA‐West2 GMMs and exhibits similar near‐source amplitude scaling. An additional benefit of mean distance is that it can provide a mechanism to account for spatially variable slip. We prospectively validate this GMM against the 2016 M "> M 7.8 Kaikōura, New Zealand, earthquake, which was not used in model development. To better understand the magnitude dependence of geometrical spreading, we employ a simple conceptual model based on fundamental principles to show that the GMM is consistent with common seismological understanding.

Bulletin of the Seismological Society of America↗
USGS2020

Earthquake magnitude and Lg Q variations between the Grenville and northern Appalachian geologic provinces of eastern Canada

This article assesses the ability of regionally specific, frequency‐dependent crustal attenuation ( ⁠ 1 / Q "> 1 / Q ⁠ ) to reduce mean magnitude discrepancies between seismic stations in the northern Appalachian and Grenville provinces (NAP and GP) of Canada. Lg Q ( f ) "> Q ( f ) is an important parameter in ground‐motion models used in probabilistic seismic hazard analysis. Discrepancies in regional magnitude estimates have long been noted to exist between stations in the two provinces for common event origins. Such discrepancies could arise from systematic site condition variations between the geologic provinces or from varying crustal attenuative properties. To evaluate the effect of frequency‐dependent anelastic attenuation, Lg Q ( f ) "> Q ( f ) on estimated magnitudes, we analyze Lg amplitudes from > 6000 "> > 6000 waveforms recorded by Grenville and northern Appalachian receivers from 420 natural earthquakes of M N "> M N magnitude 3–5.6. Waveform analysis is strictly limited to analyst‐reviewed, vertical‐component waveforms in which Lg is clearly identified, ensuring that the datasets exhibit dominant, high‐frequency energy in the Lg velocity window. Lg Q ( f ) "> Q ( f ) is found to be higher in the GP than in the northern Appalachians. In the Grenville, Q ( f ) = 761 ( ± 145 ) f 0.25 ( ± 0.014 ) "> Q ( f ) = 761 ( ± 145 ) f 0.25 ( ± 0.014 ) ⁠ , and in the northern Appalachians, attenuation is higher: Q ( f ) = 506 ( ± 172 ) f 0.33 ( ± 0.310 ) "> Q ( f ) = 506 ( ± 172 ) f 0.33 ( ± 0.310 ) . Earthquake magnitude determined using the peak amplitude of the Lg phase ( ⁠ m b L g "> m b L g ) for eastern Canada is corrected to incorporate the frequency‐dependent, regionally specific Lg Q ( f ) "> Q ( f ) determined in this study. Using the new Lg Q ( f ) "> Q ( f ) values diminishes and nearly resolves magnitude discrepancies between the provinces. Correcting regional magnitude discrepancies between provinces is critical for reliable regional seismic hazard estimates because magnitude error in a particular region could lead to increased uncertainty in seismic hazard models.

Grenville Geologic Province, Northern Appalachian ↗