Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Bulletin Seismological Society America”

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 307 records · Page 17Linked to original sources
USGS2010

Broadband ground-motion simulation using a hybrid approach

This paper describes refinements to the hybrid broadband ground-motion simulation methodology of Graves and Pitarka (2004), which combines a deterministic approach at low frequencies (f< 1 Hz) with a semistochastic approach at high frequencies (f> 1 Hz). In our approach, fault rupture is represented kinematically and incorporates spatial heterogeneity in slip, rupture speed, and rise time. The prescribed slip distribution is constrained to follow an inverse wavenumber-squared fall-off and the average rupture speed is set at 80% of the local shear-wave velocity, which is then adjusted such that the rupture propagates faster in regions of high slip and slower in regions of low slip. We use a Kostrov-like slip-rate function having a rise time proportional to the square root of slip, with the average rise time across the entire fault constrained empirically. Recent observations from large surface rupturing earthquakes indicate a reduction of rupture propagation speed and lengthening of rise time in the near surface, which we model by applying a 70% reduction of the rupture speed and increasing the rise time by a factor of 2 in a zone extending from the surface to a depth of 5 km. We demonstrate the fidelity of the technique by modeling the strong-motion recordings from the Imperial Valley, Loma Prieta, Landers, and Northridge earthquakes.

Bulletin of the Seismological Society of America↗
USGS2010

A revised ground-motion and intensity interpolation scheme for shakemap

We describe a weighted-average approach for incorporating various types of data (observed peak ground motions and intensities and estimates from groundmotion prediction equations) into the ShakeMap ground motion and intensity mapping framework. This approach represents a fundamental revision of our existing ShakeMap methodology. In addition, the increased availability of near-real-time macroseismic intensity data, the development of newrelationships between intensity and peak ground motions, and new relationships to directly predict intensity from earthquake source information have facilitated the inclusion of intensity measurements directly into ShakeMap computations. Our approach allows for the combination of (1) direct observations (ground-motion measurements or reported intensities), (2) observations converted from intensity to ground motion (or vice versa), and (3) estimated ground motions and intensities from prediction equations or numerical models. Critically, each of the aforementioned data types must include an estimate of its uncertainties, including those caused by scaling the influence of observations to surrounding grid points and those associated with estimates given an unknown fault geometry. The ShakeMap ground-motion and intensity estimates are an uncertainty-weighted combination of these various data and estimates. A natural by-product of this interpolation process is an estimate of total uncertainty at each point on the map, which can be vital for comprehensive inventory loss calculations. We perform a number of tests to validate this new methodology and find that it produces a substantial improvement in the accuracy of ground-motion predictions over empirical prediction equations alone.

Bulletin of the Seismological Society of America↗
USGS2011

A Comparison of seismic instrument noise coherence analysis techniques

The self-noise of a seismic instrument is a fundamental characteristic used to evaluate the quality of the instrument. It is important to be able to measure this self-noise robustly, to understand how differences among test configurations affect the tests, and to understand how different processing techniques and isolation methods (from nonseismic sources) can contribute to differences in results. We compare two popular coherence methods used for calculating incoherent noise, which is widely used as an estimate of instrument self-noise (incoherent noise and self-noise are not strictly identical but in observatory practice are approximately equivalent; Holcomb, 1989; Sleeman et al., 2006). Beyond directly comparing these two coherence methods on similar models of seismometers, we compare how small changes in test conditions can contribute to incoherent-noise estimates. These conditions include timing errors, signal-to-noise ratio changes (ratios between background noise and instrument incoherent noise), relative sensor locations, misalignment errors, processing techniques, and different configurations of sensor types.

Bulletin of the Seismological Society of America↗
USGS2011

Spatiotemporal earthquake clusters along the North Anatolian fault zone offshore Istanbul

We investigate earthquakes with similar waveforms in order to characterize spatiotemporal microseismicity clusters within the North Anatolian fault zone (NAFZ) in northwest Turkey along the transition between the 1999 İzmit rupture zone and the Marmara Sea seismic gap. Earthquakes within distinct activity clusters are relocated with cross-correlation derived relative travel times using the double-difference method. The spatiotemporal distribution of microearthquakes within individual clusters is resolved with relative location accuracy comparable to or better than the source size. High-precision relative hypocenters define the geometry of individual fault patches, permitting a better understanding of fault kinematics and their role in local-scale seismotectonics along the region of interest. Temporal seismic sequences observed in the eastern Sea of Marmara region suggest progressive failure of mostly nonoverlapping areas on adjacent fault patches and systematic migration of microearthquakes within clusters during the progressive failure of neighboring fault patches. The temporal distributions of magnitudes as well as the number of events follow swarmlike behavior rather than a mainshock/aftershock pattern.

Istanbul↗
USGS2012

Significant motions between GPS sites in the New Madrid region: implications for seismic hazard

Position time series from Global Positioning System (GPS) stations in the New Madrid region were differenced to determine the relative motions between stations. Uncertainties in rates were estimated using a three‐component noise model consisting of white, flicker, and random walk noise, following the methodology of Langbein, 2004. Significant motions of 0.37±0.07 (one standard error) mm/yr were found between sites PTGV and STLE, for which the baseline crosses the inferred deep portion of the Reelfoot fault. Baselines between STLE and three other sites also show significant motion. Site MCTY (adjacent to STLE) also exhibits significant motion with respect to PTGV. These motions are consistent with a model of interseismic slip of about 4 mm/yr on the Reelfoot fault at depths between 12 and 20 km. If constant over time, this rate of slip produces sufficient slip for an M 7.3 earthquake on the shallow portion of the Reelfoot fault, using the geologically derived recurrence time of 500 years. This model assumes that the shallow portion of the fault has been previously loaded by the intraplate stress. A GPS site near Little Rock, Arkansas, shows significant southward motion of 0.3–0.4 mm/yr (±0.08 mm/yr) relative to three sites to the north, indicating strain consistent with focal mechanisms of earthquake swarms in northern Arkansas.

Bulletin of the Seismological Society of America↗
USGS2011

Testing long-period ground-motion simulations of scenario earthquakes using the Mw 7.2 El Mayor-Cucapah mainshock: Evaluation of finite-fault rupture characterization and 3D seismic velocity models

Using a suite of five hypothetical finite-fault rupture models, we test the ability of long-period ( T >2.0 s) ground-motion simulations of scenario earthquakes to produce waveforms throughout southern California consistent with those recorded during the 4 April 2010 M w 7.2 El Mayor-Cucapah earthquake. The hypothetical ruptures are generated using the methodology proposed by Graves and Pitarka (2010) and require, as inputs, only a general description of the fault location and geometry, event magnitude, and hypocenter, as would be done for a scenario event. For each rupture model, two Southern California Earthquake Center three-dimensional community seismic velocity models (CVM-4m and CVM-H62) are used, resulting in a total of 10 ground-motion simulations, which we compare with recorded ground motions. While the details of the motions vary across the simulations, the median levels match the observed peak ground velocities reasonably well, with the standard deviation of the residuals generally within 50% of the median. Simulations with the CVM-4m model yield somewhat lower variance than those with the CVM-H62 model. Both models tend to overpredict motions in the San Diego region and underpredict motions in the Mojave desert. Within the greater Los Angeles basin, the CVM-4m model generally matches the level of observed motions, whereas the CVM-H62 model tends to overpredict the motions, particularly in the southern portion of the basin. The variance in the peak velocity residuals is lowest for a rupture that has significant shallow slip (<5 km depth), whereas the variance in the residuals is greatest for ruptures with large asperities below 10 km depth. Overall, these results are encouraging and provide confidence in the predictive capabilities of the simulation methodology, while also suggesting some regions in which the seismic velocity models may need improvement.

California↗
USGS2012

A terrain-based site characterization map of California with implications for the contiguous United States

We present an approach based on geomorphometry to predict material properties and characterize site conditions using the V S 30 parameter (time‐averaged shear‐wave velocity to a depth of 30 m). Our framework consists of an automated terrain classification scheme based on taxonomic criteria (slope gradient, local convexity, and surface texture) that systematically identifies 16 terrain types from 1‐km spatial resolution (30 arcsec) Shuttle Radar Topography Mission digital elevation models (SRTMDEMs). Using 853 V S 30 values from California, we apply a simulation‐based statistical method to determine the mean V S 30 for each terrain type in California. We then compare the V S 30 values with models based on individual proxies, such as mapped surface geology and topographic slope, and show that our systematic terrain‐based approach consistently performs better than semiempirical estimates based on individual proxies. To further evaluate our model, we apply our California‐based estimates to terrains of the contiguous United States. Comparisons of our estimates with 325 V S 30 measurements outside of California, as well as estimates based on the topographic slope model, indicate our method to be statistically robust and more accurate. Our approach thus provides an objective and robust method for extending estimates of V S 30 for regions where in situ measurements are sparse or not readily available.

California↗
USGS2011

Shear-wave velocity characterization of the USGS Hawaiian strong-motion network on the Island of Hawaii and development of an NEHRP site-class map

To assess the level and nature of ground shaking in Hawaii for the purposes of earthquake hazard mitigation and seismic design, empirical ground-motion prediction models are desired. To develop such empirical relationships, knowledge of the subsurface site conditions beneath strong-motion stations is critical. Thus, as a first step to develop ground-motion prediction models for Hawaii, spectral-analysis-of-surface-waves (SASW) profiling was performed at the 22 free-field U.S. Geological Survey (USGS) strong-motion sites on the Big Island to obtain shear-wave velocity ( V S ) data. Nineteen of these stations recorded the 2006 Kiholo Bay moment magnitude (M) 6.7 earthquake, and 17 stations recorded the triggered M 6.0 Mahukona earthquake. V S profiling was performed to reach depths of more than 100 ft. Most of the USGS stations are situated on sites underlain by basalt, based on surficial geologic maps. However, the sites have varying degrees of weathering and soil development. The remaining strong-motion stations are located on alluvium or volcanic ash. V S 30 (average V S in the top 30 m) values for the stations on basalt ranged from 906 to 1908 ft/s [National Earthquake Hazards Reduction Program (NEHRP) site classes C and D], because most sites were covered with soil of variable thickness. Based on these data, an NEHRP site-class map was developed for the Big Island. These new V S data will be a significant input into an update of the USGS statewide hazard maps and to the operation of ShakeMap on the island of Hawaii.

Hawai'i↗
USGS2012

Including foreshocks and aftershocks in time-independent probabilistic seismic hazard analyses

Time‐independent probabilistic seismic‐hazard analysis treats each source as being temporally and spatially independent; hence foreshocks and aftershocks, which are both spatially and temporally dependent on the mainshock, are removed from earthquake catalogs. Yet, intuitively, these earthquakes should be considered part of the seismic hazard, capable of producing damaging ground motions. In this study, I consider the mainshock and its dependents as a time‐independent cluster, each cluster being temporally and spatially independent from any other. The cluster has a recurrence time of the mainshock; and, by considering the earthquakes in the cluster as a union of events, dependent events have an opportunity to contribute to seismic ground motions and hazard. Based on the methods of the U.S. Geological Survey for a high‐hazard site, the inclusion of dependent events causes ground motions that are exceeded at probability levels of engineering interest to increase by about 10% but could be as high as 20% if variations in aftershock productivity can be accounted for reliably.

Bulletin of the Seismological Society of America↗
USGS2011

Specifying initial stress for dynamic heterogeneous earthquake source models

Dynamic rupture calculations using heterogeneous stress drop that is random and self-similar with a power-law spatial spectrum have great promise of producing realistic ground-motion predictions. We present procedures to specify initial stress for random events with a target rupture length and target magnitude. The stress function is modified in the depth dimension to account for the brittle-ductile transition at the base of the seismogenic zone. Self-similar fluctuations in stress drop are tied in this work to the long-wavelength stress variation that determines rupture length. Heterogeneous stress is related to friction levels in order to relate the model to physical concepts. In a variant of the model, there are high-stress asperities with low background stress. This procedure has a number of advantages: (1) rupture stops naturally, not at artificial barriers; (2) the amplitude of short-wavelength fluctuations of stress drop is not arbitrary: the spectrum is fixed to the long-wavelength fluctuation that determines rupture length; and (3) large stress drop can be confined to asperities occupying a small fraction of the total rupture area, producing slip distributions with enhanced peaks.

Bulletin of the Seismological Society of America↗
USGS2011

Model and parametric uncertainty in source-based kinematic models of earthquake ground motion

Four independent ground-motion simulation codes are used to model the strong ground motion for three earthquakes: 1994 M w 6.7 Northridge, 1989 M w 6.9 Loma Prieta, and 1999 M w 7.5 Izmit. These 12 sets of synthetics are used to make estimates of the variability in ground-motion predictions. In addition, ground-motion predictions over a grid of sites are used to estimate parametric uncertainty for changes in rupture velocity. We find that the combined model uncertainty and random variability of the simulations is in the same range as the variability of regional empirical ground-motion data sets. The majority of the standard deviations lie between 0.5 and 0.7 natural-log units for response spectra and 0.5 and 0.8 for Fourier spectra. The estimate of model epistemic uncertainty, based on the different model predictions, lies between 0.2 and 0.4, which is about one-half of the estimates for the standard deviation of the combined model uncertainty and random variability. Parametric uncertainty, based on variation of just the average rupture velocity, is shown to be consistent in amplitude with previous estimates, showing percentage changes in ground motion from 50% to 300% when rupture velocity changes from 2.5 to 2.9 km/s. In addition, there is some evidence that mean biases can be reduced by averaging ground-motion estimates from different methods.

Bulletin of the Seismological Society of America↗
USGS2012

Long-term creep rates on the Hayward Fault: Evidence for controls on the size and frequency of large earthquakes

The Hayward fault (HF) in California exhibits large (Mw 6.5&ndash;7.1) earthquakes with short recurrence times (161&plusmn;65 yr), probably kept short by a 26%&ndash;78% aseismic release rate (including postseismic). Its interseismic release rate varies locally over time, as we infer from many decades of surface creep data. Earliest estimates of creep rate, primarily from infrequent surveys of offset cultural features, revealed distinct spatial variation in rates along the fault, but no detectable temporal variation. Since the 1989 Mw 6.9 Loma Prieta earthquake (LPE), monitoring on 32 alinement arrays and 5 creepmeters has greatly improved the spatial and temporal resolution of creep rate. We now identify significant temporal variations, mostly associated with local and regional earthquakes. The largest rate change was a 6‐yr cessation of creep along a 5‐km length near the south end of the HF, attributed to a regional stress drop from the LPE, ending in 1996 with a 2‐cm creep event. North of there near Union City starting in 1991, rates apparently increased by 25% above pre‐LPE levels on a 16‐km‐long reach of the fault. Near Oakland in 2007 an Mw 4.2 earthquake initiated a 1&ndash;2 cm creep event extending 10&ndash;15 km along the fault. Using new better‐constrained long‐term creep rates, we updated earlier estimates of depth to locking along the HF. The locking depths outline a single, &sim;50‐km‐long locked or retarded patch with the potential for an Mw&sim;6.8 event equaling the 1868 HF earthquake. We propose that this inferred patch regulates the size and frequency of large earthquakes on HF.

California↗
USGS2012

A new perspective on the geometry of the San Andreas Fault in southern California and its relationship to lithospheric structure

The widely held perception that the San Andreas fault (SAF) is vertical or steeply dipping in most places in southern California may not be correct. From studies of potential‐field data, active‐source imaging, and seismicity, the dip of the SAF is significantly nonvertical in many locations. The direction of dip appears to change in a systematic way through the Transverse Ranges: moderately southwest (55°–75°) in the western bend of the SAF in the Transverse Ranges (Big Bend); vertical to steep in the Mojave Desert; and moderately northeast (37°–65°) in a region extending from San Bernardino to the Salton Sea, spanning the eastern bend of the SAF in the Transverse Ranges. The shape of the modeled SAF is crudely that of a propeller. If confirmed by further studies, the geometry of the modeled SAF would have important implications for tectonics and strong ground motions from SAF earthquakes. The SAF can be traced or projected through the crust to the north side of a well documented high‐velocity body (HVB) in the upper mantle beneath the Transverse Ranges. The north side of this HVB may be an extension of the plate boundary into the mantle, and the HVB would appear to be part of the Pacific plate.

California↗
USGS2012

Characterization of intrabasin faulting and deformation for earthquake hazards in southern Utah Valley, Utah, from high-resolution seismic imaging

We conducted active and passive seismic imaging investigations along a 5.6-km-long, east&ndash;west transect ending at the mapped trace of the Wasatch fault in southern Utah Valley. Using two-dimensional (2D) P-wave seismic reflection data, we imaged basin deformation and faulting to a depth of 1.4 km and developed a detailed interval velocity model for prestack depth migration and 2D ground-motion simulations. Passive-source microtremor data acquired at two sites along the seismic reflection transect resolve S-wave velocities of approximately 200 m/s at the surface to about 900 m/s at 160 m depth and confirm a substantial thickening of low-velocity material westward into the valley. From the P-wave reflection profile, we interpret shallow (100&ndash;600 m) bedrock deformation extending from the surface trace of the Wasatch fault to roughly 1.5 km west into the valley. The bedrock deformation is caused by multiple interpreted fault splays displacing fault blocks downward to the west of the range front. Further west in the valley, the P-wave data reveal subhorizontal horizons from approximately 90 to 900 m depth that vary in thickness and whose dip increases with depth eastward toward the Wasatch fault. Another inferred fault about 4 km west of the mapped Wasatch fault displaces horizons within the valley to as shallow as 100 m depth. The overall deformational pattern imaged in our data is consistent with the Wasatch fault migrating eastward through time and with the abandonment of earlier synextensional faults, as part of the evolution of an inferred 20-km-wide half-graben structure within Utah Valley. Finite-difference 2D modeling suggests the imaged subsurface basin geometry can cause fourfold variation in peak ground velocity over distances of 300 m.

Utah↗
USGS2013

Assessing the location and magnitude of the 20 October 1870 Charlevoix, Quebec, earthquake

The Charlevoix, Quebec, earthquake of 20 October 1870 caused damage to several towns in Quebec and was felt throughout much of southeastern Canada and along the U.S. Atlantic seaboard from Maine to Maryland. Site‐specific damage and felt reports from Canadian and U.S. cities and towns were used in analyses of the location and magnitude of the earthquake. The macroseismic center of the earthquake was very close to Baie‐St‐Paul, where the greatest damage was reported, and the intensity magnitude M I was found to be 5.8, with a 95% probability range of 5.5–6.0. After corrections for epicentral‐distance differences are applied, the modified Mercalli intensity (MMI) data for the 1870 earthquake and for the moment magnitude M 6.2 Charlevoix earthquake of 1925 at common sites show that on average, the MMI readings are about 0.8 intensity units smaller for the 1870 earthquake than for the 1925 earthquake, suggesting that the 1870 earthquake was M I 5.7. A similar comparison of the MMI data for the 1870 earthquake with the corresponding data for the M 5.9 1988 Saguenay event suggests that the 1870 earthquake was M I 6.0. These analyses all suggest that the magnitude of the 1870 Charlevoix earthquake is between M I 5.5 and M I 6.0, with a best estimate of M I 5.8.

Bulletin of the Seismological Society of America↗
USGS2012

Contemporary seismicity in and around the Yakima-Fold-and-Thrust Belt in eastern Washington

We examined characteristics of routinely cataloged seismicity from 1970 to the present in and around the Yakima fold‐and‐thrust belt (YFTB) in eastern Washington to determine if the characteristics of contemporary seismicity provide clues about regional‐scale active tectonics or about more localized, near‐surface processes. We employed new structural and hydrologic models of the Columbia River basalts (CRB) and found that one‐third to one‐half of the cataloged earthquakes occur within the CRB and that these CRB earthquakes exhibit significantly more clustered, and swarmlike, behavior than those outside. These results and inferences from published studies led us to hypothesize that clustered seismicity is likely associated with hydrologic changes in the CRB, which hosts the regional aquifer system. While some general features of the regional groundwater system support this hypothesis, seismicity patterns and mapped long‐term changes in groundwater levels and present‐day irrigation neither support nor refute it. Regional tectonic processes and crustal‐scale structures likely influence the distribution of earthquakes both outside and within the CRB as well. We based this inference on qualitatively assessed alignments between the dominant northwest trends in the geologic structure and the seismicity generally and between specific faults and characteristics of the 2009 Wooded Island swarm and aseismic slip, which is the only cluster studied in detail and the most vigorous since regional monitoring began.

Washington↗