Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geological Society of America Bulletin”

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 829 records · Page 46Linked to original sources

Late Holocene slip rate of the Mojave section of the San Andreas Fault near Palmdale, California

The geologic slip rate on the Mojave section of the San Andreas fault is poorly constrained, despite its importance for understanding earthquake hazard, apparent discrepancies between geologic and geodetic slip rates along this fault section, and long‐term fault interactions in southern California. Here, we use surficial geologic mapping, excavations, and radiocarbon and luminescence dating to quantify the displacements and ages of late Holocene landforms offset by the fault at three sites. At the Ranch Center site, the slip rate is determined using the base of a fan marking incision and deflection of an ephemeral channel. At the adjacent Key Slide site, the margin of a landslide deposited on indigenous fire hearths provides a minimum rate. At the X‐12 site, the slip rate is determined from a channel that incised into a broad fan surface, and is deflected and beheaded by the fault. We use maximum–minimum bounds on both the displacement and age of each offset feature to calculate slip rate for each site independently. Overlap of the three independent rate ranges yields a rate of 33–39 mm/yr over the last 3 ka, under the assumption that the sites share a common history, given their proximity. Considered in sequence, site‐level epistemic uncertainties in the data permit but do not require a rate increase since ∼1200 cal B.P. Modest rate changes can be explained by aleatory variability in earthquake timing and magnitude; larger changes could suggest a shared regional variation with the Garlock and other faults. The new late Holocene slip rates are consistent with geodetic model estimates that include a viscoelastic crust and earthquake cycle effects. The geologic slip rates also provide average slip over dozens of earthquake cycles—a key constraint for long‐term earthquake rupture forecasts.

California↗

Magnetic field variations in Alaska: Recording space weather events on seismic stations in Alaska

Seismometers are highly sensitive instruments to not only ground motion but also many other nonseismic noise sources (e.g., temperature, pressure, and magnetic field variations). We show that the Alaska component of the Transportable Array is particularly susceptible to recording magnetic storms and other space weather events because the sensors used in this network are unshielded and magnetic flux variations are stronger at higher latitudes. We also show that vertical‐component seismic records across Alaska are directly recording magnetic field variations between 40 and 800 s period as opposed to actual ground motion during geomagnetic events with sensitivities ranging from 0.004 to 0.48    ( m / s 2 ) / T "> 0.48 ( m / s 2 ) / T 0.48 (m/s2)/T ⁠ . These sensitivities were found on a day where the root mean square variation in the magnetic field was 225 nT. Using a method developed by Forbriger (2007 , his section 3.1), we show that improving vertical seismic resolution of an unshielded sensor by as much as 10 dB in the 100–400 s period band using magnetic data from a collocated three‐component magnetometer is possible. However, due to large spatial variations in Earth’s magnetic field, this methodology becomes increasingly ineffective as the distance between the seismometer and magnetometer increases (no more than 200 km separation). A potential solution to this issue may be to incorporate relatively low‐cost magnetometers as an additional environmental data stream at high‐latitude seismic stations. We demonstrate that the Bartington Mag‐690 sensors currently deployed at Global Seismographic Network sites are not only acceptable for performing corrections to seismic data, but are also capable of recording many magnetic field signals with similar signal‐to‐noise ratios, in the 20–1000 s period band, as the observatory grade magnetometers operated by the U.S. Geological Survey Geomagnetism Program. This approach would densify magnetic field observations and could also contribute to space weather monitoring by supplementing highly calibrated magnetometers with additional sensors.

Alaska↗

Finite-fault source inversion using teleseismic P waves: Simple parameterization and rapid analysis

We examine the ability of teleseismic P waves to provide a timely image of the rupture history for large earthquakes using a simple, 2D finite‐fault source parameterization. We analyze the broadband displacement waveforms recorded for the 2010 M w ∼7 Darfield (New Zealand) and El Mayor‐Cucapah (Baja California) earthquakes using a single planar fault with a fixed rake. Both of these earthquakes were observed to have complicated fault geometries following detailed source studies conducted by other investigators using various data types. Our kinematic, finite‐fault analysis of the events yields rupture models that similarly identify the principal areas of large coseismic slip along the fault. The results also indicate that the amount of stabilization required to spatially smooth the slip across the fault and minimize the seismic moment is related to the amplitudes of the observed P waveforms and can be estimated from the absolute values of the elements of the coefficient matrix. This empirical relationship persists for earthquakes of different magnitudes and is consistent with the stabilization constraint obtained from the L‐curve in Tikhonov regularization. We use the relation to estimate the smoothing parameters for the 2011 M w 7.1 East Turkey, 2012 M w 8.6 Northern Sumatra, and 2011 M w 9.0 Tohoku, Japan, earthquakes and invert the teleseismic P waves in a single step to recover timely, preliminary slip models that identify the principal source features observed in finite‐fault solutions obtained by the U.S. Geological Survey National Earthquake Information Center (USGS/NEIC) from the analysis of body‐ and surface‐wave data. These results indicate that smoothing constraints can be estimated a priori to derive a preliminary, first‐order image of the coseismic slip using teleseismic records.

Darfield, El Mayor-Cucapah↗

Developing a map of geologically defined site-condition categories for California

Consideration of site conditions is a vital step in analyzing and predicting earthquake ground motion. The importance of amplification by soil conditions has long been recognized, but though many seismic-instrument sites have been characterized by their geologic conditions, there has been no consistent, simple classification applied to all sites. As classification of sites by shear-wave velocity has become more common, the need to go back and provide a simple uniform classification for all stations has become apparent. Within the Pacific Earthquake Engineering Research Center's Next Generation Attenuation equation project, developers of attenuation equations recognized the need to consider site conditions and asked that the California Geological Survey provide site conditions information for all stations that have recorded earthquake ground motion in California. To provide these estimates, we sorted the available shear-wave velocity data by geologic unit, generalized the geologic units, and prepared a map so that we could use the extent of the map units to transfer the velocity characteristics from the sites where they were measured to sites on the same or similar materials. This new map is different from the California Geological Survey "preliminary site-conditions map of California" in that 19 geologically defined categories are used, rather than National Earthquake Hazards Reduction Program categories. Although this map does not yet cover all of California, when completed it may provide a basis for more precise consideration of site conditions in ground-motion calculations.

Bulletin of the Seismological Society of America↗

Length-scale-dependent relationships between VS30 and topographic slopes in southern California

In engineering seismology, the time‐averaged shear‐wave velocity ( ⁠ V S ⁠ ) of the upper 30 m of the crust ( ⁠V S30⁠ ) is the primary parameter used in ground‐motion models to predict seismic site effects. ⁠V S30⁠ is typically derived from in situ recordings of V S ⁠ , although proxy‐based approaches (using geologic and/or geomorphometric classifications) are provisionally adopted when measurement‐based ⁠V S30⁠ are sparse or not readily available. Despite the acceptance of proxy approaches, there are limited studies that examine the empirical relationships between ⁠V S30⁠ and topographic attributes measured from various length scales and different resolutions of the digital elevation model. In this study, we examine the relationships between compiled ⁠V S30⁠ measurements from 218 sites in southern California and topographic metrics of slope and relief measured over various length scales. We find that the correlations between topographic metrics and ⁠V S30⁠ are weak but statistically significant. The correlations are improved when topographic slopes and relief are measured over length scales longer than typical hillslopes and ⁠V S30⁠ sites are separated by different geologic groups. This is likely because ⁠V S30⁠ ⁠ , especially on the rock sites, is better reflected in topographic metrics that capture large‐scale topographic relief, as well as landscape positions such as hilltops and valley bottoms.

California↗

The seismic response of the Los Angeles basin, California

Using strong-motion data recorded in the Los Angeles region from the 1992 ( M w 7.3) Landers earthquake, we have tested the accuracy of existing three-dimensional (3D) velocity models on the simulation of long-period (≧2 sec) ground motions in the Los Angeles basin and surrounding San Fernando and San Gabriel Valleys. First, the overall pattern and degree of long-period excitation of the basins were identified in the observations. Within the Los Angeles basin, the recorded amplitudes are about three to four times larger than at sites outside the basins; amplitudes within the San Fernando and San Gabriel Valleys are nearly a factor of 3 greater than surrounding bedrock sites. Then, using a 3D finite-difference numerical modeling approach, we analyzed how variations in 3D earth structure affect simulated waveforms, amplitudes, and the fit to the observed patterns of amplification. Significant differences exist in the 3D velocity models of southern California that we tested (Magistrale et al. , 1996; Graves, 1996a; Hauksson and Haase, 1997). Major differences in the models include the velocity of the assumed background models; the depth of the Los Angeles basin; and the depth, location, and geometry of smaller basins. The largest disparities in the response of the models are seen for the San Fernando Valley and the deepest portion of the Los Angeles basin. These arise in large part from variations in the structure of the basins, particularly the effective depth extent, which is mainly due to alternative assumptions about the nature of the basin sediment fill. The general ground-motion characteristics are matched by the 3D model simulations, validating the use of 3D modeling with geologically based velocity-structure models. However, significant shortcomings exist in the overall patterns of amplification and the duration of the long-period response. The successes and limitations of the models for reproducing the recorded ground motions as discussed provide the basis and direction for necessary improvements to earth structure models, whether geologically or tomographically derived. The differences in the response of the earth models tested also translate to variable success in the ability to successfully model the data and add uncertainty to estimates of the basin response given input “scenario” earthquake source models.

California↗

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↗

A study of possible ground-motion amplification at the Coyote Lake Dam, California

The abutment site at the Coyote Lake Dam recorded an unusually large peak acceleration of 1.29g during the 1984 Morgan Hill earthquake. Following this earthquake another strong-motion station was installed about 700 m downstream from the abutment station. We study all events (seven) recorded on these stations, using ratios of peak accelerations, spectral ratios, and particle motion polarization (using holograms) to investigate the relative ground motion at the two sites. We find that in all but one case the motion at the abutment site is larger than the downstream site over a broad frequency band. The polarizations are similar for the two sites for a given event, but can vary from one event to another. This suggests that the dam itself is not strongly influencing the records. Although we can be sure that the relative motion is usually larger at the abutment site, we cannot conclude that there is anomalous site amplification at the abutment site. The downstream site could have lower-than-usual near-surface amplifications. On the other hand, the geology near the abutment site is extremely complex and includes fault slivers, with rapid lateral changes in materials and presumably seismic velocities. For this reason alone, the abutment site should not be considered a normal free-field site.

Bulletin of the Seismological Society of America↗

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↗

Maximum magnitude ( M max ) in the central and eastern United States for the 2014 U.S. Geological Survey Hazard Model

Probabilistic seismic‐hazard assessment (PSHA) requires an estimate of M max , the moment magnitude M of the largest earthquake that could occur within a specified area. Sparse seismicity hinders M max estimation in the central and eastern United States (CEUS) and tectonically similar regions worldwide (stable continental regions [SCRs]). A new global catalog of moderate‐to‐large SCR earthquakes is analyzed with minimal assumptions about enigmatic geologic controls on SCR M max . An earlier observation that SCR earthquakes of M 7.0 and larger occur in young (250–23 Ma) passive continental margins and associated rifts but not in cratons is not strongly supported by the new catalog. SCR earthquakes of M 7.5 and larger are slightly more numerous and reach slightly higher M in young passive margins and rifts than in cratons. However, overall histograms of M from young margins and rifts and from cratons are statistically indistinguishable. This conclusion is robust under uncertainties in M , the locations of SCR boundaries, and which of two available global SCR catalogs is used. The conclusion stems largely from recent findings that (1) large southeast Asian earthquakes once thought to be SCR were in actively deforming crust and (2) long escarpments in cratonic Australia were formed by prehistoric faulting. The 2014 seismic‐hazard model of the U.S. Geological Survey represents CEUS M max as four‐point probability distributions. The distributions have weighted averages of M 7.0 in cratons and M 7.4 in passive margins and rifts. These weighted averages are consistent with M max estimates of other SCR PSHAs of the CEUS, southeastern Canada, Australia, and India.

Bulletin of the Seismological Society of America↗

Coda duration magnitudes in central California: An empirical approach

A new empirical coda magnitude M D is presented for the central California seismic network (CALNET), in which station corrections are introduced to reduce the influence of hypocenter and station distributions. The new magnitude scale is obtained by relating the lapse-time τ^ to local magnitudes M L using data from 55 earthquakes with 1.1 ≦ M L ≤ 5.6 from June 1977 through 1981. Magnitude estimates are significantly improved by introducing corrections for the seismograph instrument attenuation setting α ij and a site correction δ j . The magnitude of the i th event 〈 M D 〉, is the median value of individual M D ij = -1.03 + 2.10log 10 τˆ ij + 0.00268 τˆ ij + α ij + δ j , where τ^ ij is the lapse-time measured at the j th station. Local site effects influence the duration of the coda and cause stations to either overestimate or underestimate 〈 M D 〉, typically up to 0.5 magnitude units. The distribution of the site corrections is spatially correlated, changing coherently across major tectonic and geologic features such as the San Andreas fault. It appears that the site corrections are influenced by the physical attenuation properties near the station, seemingly contradicting the concept that the coda is related to back-scattering and attenuation in a volume with dimensions on the order of the hypocenter to receiver distance.

California↗

MMI attenuation and historical earthquakes in the basin and range province of western North America

Earthquakes in central Nevada (1932-1959) were used to develop a modified Mercalli intensity (MMI) attenuation model for estimating moment magnitude M for earthquakes in the Basin and Range province of interior western North America. M is 7.4-7.5 for the 26 March 1872 Owens Valley, California, earthquake, in agreement with Beanland and Clark's (1994) M 7.6 that was estimated from geologic field observations. M is 7.5 for the 3 May 1887 Sonora, Mexico, earthquake, in agreement with Natali and Sbar's (1982) M 7.4 and Suter's (2006) M 7.5, both estimated from geologic field observations. MMI at sites in California for earthquakes in the Nevada Basin and Range apparently are not much affected by the Sierra Nevada except at sites near the Sierra Nevada where MMI is reduced. This reduction in MMI is consistent with a shadow zone produced by the root of the Sierra Nevada. In contrast, MMI assignments for earthquakes located in the eastern Sierra Nevada near the west margin of the Basin and Range are greater than predicted at sites in California. These higher MMI values may result from critical reflections due to layering near the base of the Sierra Nevada.

Bulletin of the Seismological Society of America↗

How to determine an earthquake rate budget? A comparison of geodetic, geologic, and seismologic moment rates

How much tectonic geodetic strain is accommodated as earthquakes? Can estimates of tectonic geodetic strain be used in probabilistic seismic hazard analyses? Understanding the balance of moment accumulation (i.e., geodetic moment) and accommodation (i.e., geologic and seismologic moment) is a necessary step toward leveraging corrected geodetic strain rates to quantify earthquake budgets. We investigate these questions using updated geodetic, geologic, and seismologic datasets compiled for the western United States (WUS) portion of the 2023 U.S. National Seismic Hazard Model (NSHM23‐WUS) to provide new views of moment rates across the WUS. We compare moment rates using a floating footprint approach across the extent of the NSHM23‐WUS fault sections database. We make these comparisons using “moment rate deviation (MRD),” which is defined as a surplus or deficit of geodetic moment rate relative to the sum of geologic and seismologic moment rates. We find that, in general, higher strain rate areas of the WUS (i.e., plate boundary) have smaller values of MRD, meaning that geodetic moment rates approximate the earthquake record preserved by geologic and seismologic moment rates well. In contrast, the lower strain rate, eastern portions of our analysis domain (i.e., intraplate) have very high values of MRD, with geodetic moment rates exceeding measured geologic and seismologic moment rates. We find that geodetic moment rates appear to be representative of the earthquake budget in relatively well‐characterized, interplate regions, and geodetic estimates of off‐fault deformation could be considered in seismic hazard analyses where MRD is minimized.

western United States↗

Normal-faulting slip maxima and stress-drop variability: A geological perspective

We present an empirical estimate of maximum slip in continental normal-faulting earthquakes and present evidence that stress drop in intraplate extensional environments is dependent on fault maturity. A survey of reported slip in historical earthquakes globally and in latest Quaternary paleoearthquakes in the Western Cordillera of the United States indicates maximum vertical displacements as large as 6–6.5 m. A difference in the ratio of maximum-to-mean displacements between data sets of prehistoric and historical earthquakes, together with constraints on bias in estimates of mean paleodisplacement, suggest that applying a correction factor of 1.4±0.3 to the largest observed displacement along a paleorupture may provide a reasonable estimate of the maximum displacement. Adjusting the largest paleodisplacements in our regional data set (∼6 m) by a factor of 1.4 yields a possible upper-bound vertical displacement for the Western Cordillera of about 8.4 m, although a smaller correction factor may be more appropriate for the longest ruptures. Because maximum slip is highly localized along strike, if such large displacements occur, they are extremely rare. Static stress drop in surface-rupturing earthquakes in the Western Cordillera, as represented by maximum reported displacement as a fraction of modeled rupture length, appears to be larger on normal faults with low cumulative geologic displacement (<2 km) and larger in regions such as the Rocky Mountains, where immature, low-throw faults are concentrated. This conclusion is consistent with a growing recognition that structural development influences stress drop and indicates that this influence is significant enough to be evident among faults within a single intraplate environment.

Bulletin of the Seismological Society of America↗

Localized surface disruptions observed by InSAR during strong earthquakes in Java and Hawai'i

Interferometric Synthetic Aperture Radar data spanning strong earthquakes on the islands of Java and Hawai‘i in 2006 reveal patches of subsidence and incoherence indicative of localized ground failure. Interferograms spanning the 26 May 2006 Java earthquake suggest an area of about 7.5 km 2 of subsidence (~2 cm) and incoherence south of the city of Yogyakarta that correlates with significant damage to housing, high modeled peak ground accelerations, and poorly consolidated geologic deposits. The subsidence and incoherence is inferred to be a result of intense shaking and/or damage. At least five subsidence patches on the west side of the Island of Hawai‘i, ranging 0.3–2.2 km 2 in area and 3–8 cm in magnitude, occurred as a result of a pair of strong earthquakes on 15 October 2006. Although no felt reports or seismic data are available from the areas in Hawai‘i, the Java example suggests that the subsidence patches indicate areas of amplified earthquake shaking. Surprisingly, all subsidence areas in Hawai‘i were limited to recent, and supposedly stable, lava flows and may reflect geological conditions not detectable at the surface. In addition, two ‘a‘ā lava flows in Hawai‘i were partially incoherent in interferograms spanning the earthquakes, indicating surface disruption as a result of the earthquake shaking. Coearthquake incoherence of rubbly deposits, like ‘a‘ā flows, should be explored as a potential indicator of earthquake intensity and past strong seismic activity.

Bulletin of the Seismological Society of America↗

Slip rate on the San Diego trough fault zone, inner California Borderland, and the 1986 Oceanside earthquake swarm revisited

The San Diego trough fault zone (SDTFZ) is part of a 90-km-wide zone of faults within the inner California Borderland that accommodates motion between the Pacific and North American plates. Along with most faults offshore southern California, the slip rate and paleoseismic history of the SDTFZ are unknown. We present new seismic reflection data that show that the fault zone steps across a 5-km-wide stepover to continue for an additional 60 km north of its previously mapped extent. The 1986 Oceanside earthquake swarm is located within the 20-km-long restraining stepover. Farther north, at the latitude of Santa Catalina Island, the SDTFZ bends 20&deg; to the west and may be linked via a complex zone of folds with the San Pedro basin fault zone (SPBFZ). In a cooperative program between the U.S. Geological Survey (USGS) and the Monterey Bay Aquarium Research Institute (MBARI), we measure and date the coseismic offset of a submarine channel that intersects the fault zone near the SDTFZ&ndash;SPBFZ junction. We estimate a horizontal slip rate of about 1:5  0:3 mm=yr over the past 12,270 yr.

California↗

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗

Changes in crustal seismic deformation rates associated with the 1964 Great Alaska earthquake

We calculated seismic moment rates from crustal earthquake information for the upper Cook Inlet region, including Anchorage, Alaska, for the 30 yr prior to and 36 yr following the 1964 Great Alaska earthquake. Our results suggest over a factor of 1000 decrease in seismic moment rate (in units of dyne centimeters per year) following the 1964 mainshock. We used geologic information on structures within the Cook Inlet basin to estimate a regional geologic moment rate, assuming the structures extend to 30 km depth and have near-vertical dips. The geologic moment rates could underestimate the true rates by up to 70% since it is difficult determine the amount of horizontal offset that has occurred along many structures within the basin. Nevertheless, the geologic moment rate is only 3-7 times lower than the pre-1964 seismic moment rate, suggesting the 1964 mainshock has significantly slowed regional crustal deformation. If we compare the geologic moment rate to the post-1964 seismic moment rate, the moment rate deficit over the past 36 yr is equivalent to a moment magnitude 6.6-7.0 earthquake. These observed differences in moment rates highlight the difficulty in using seismicity in the decades following a large megathrust earthquake to adequately characterize long-term crustal deformation.

Alaska↗