USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Seismic Record”

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 1,711 records · Page 95Linked to original sources
USGS2012

Source parameters of the 1980 Mammoth Lakes, California, earthquake sequence

From the more than 1500 Mammoth Lakes earthquakes recorded on three-component digital seismographs (Spudich et al., 1981), 150 were used in an analysis of the locations, mechanism, and source parameters. A composite fault plane solution of nine earthquakes 3.9 ≤ M ≤ 5.1 defines a right-lateral strike slip mechanism on a steeply dipping nearly east-west plane striking S75°E or left-lateral strike slip on a nearly north-south plane striking N10°E. Vertical cross sections of well-located aftershocks indicate possibly three east-west planes that coincide with the locations of the four largest earthquakes with M L ≥ 6.0. Using the spectral analysis of S waves (Brune, 1970), source parameters for 67 earthquakes were determined. Forty-eight had magnitudes greater than or equal to 3.0. Seismic moments ranged from 9.20×10 18 dyn cm to 2.33×10 24 dyn cm. Earthquakes with seismic moment greater than about 1.0×10 21 dyn cm had nearly constant stress drops (≃ 50 bars); earthquakes with seismic moment less than about 1.0×10 21 dyn cm had stress drops that apparently decrease as seismic moment decreases.

California↗
USGS2022

Using a grid-search approach to validate the Graves-Pitarka broadband simulation method

This work assesses the ability of the Graves–Pitarka simulation approach to reproduce observed ground motions for 12 California and Baja California earthquakes. A total of 240 realizations are computed for each earthquake and compared with recorded strong motions from near-fault sites. In addition to spatial variability in slip, each realization samples from discrete combinations of average rupture speed, rise time, and down-dip fault width. Ground motions for each realization of a given earthquake are compared to the observations using a pseudo-spectral acceleration goodness-of-fit (GoF) metric. Parameters for the lowest misfit cases are then tabulated to develop relations for estimating median values and ranges for future applications. The results are generally consistent with published scaling relations for estimating rupture dimensions as a function of magnitude. Additionally, I find the average rise time scales as 2.3 × α T × 10 − 9 × M o 1 / 3 "> 2.3 × α T × 10 - 9 × M o 1 / 3 ( M o is seismic moment in dyne-cm) and average rupture speed as ( 0.765 ± 0.075 ) × α T − 1 × V s "> ( 0.765 ± 0.075 ) × α T - 1 × V s ( V s is local shear wave velocity), where 2.3 × α T × 10 − 9 × M o 1 / 3 "> α T is a mechanism adjustment ranging from a value of 1.0 for pure strike-slip to 0.9 for pure reverse-slip cases. There are cases where slightly different combinations of parameters produce equally good fits to the observations, demonstrating the non-uniqueness of using a single GoF metric in this approach. Nonetheless, the results reinforce the importance of adequately sampling ranges of rupture parameters when performing validations, as well as when simulating ground motions for future events.

Earth, Planets and Space↗
USGS2024

Data-driven adjustments for combined use of NGA-East hard-rock ground motion and site amplification models

Model development in the Next Generation Attenuation-East (NGA-East) project included two components developed concurrently and independently: (1) earthquake ground-motion models (GMMs) that predict the median and aleatory variability of various intensity measures conditioned on magnitude and distance, derived for a reference hard-rock site condition with an average shear-wave velocity in the upper 30 m ( V S30 ) = 3000 m/s; and (2) a site amplification model that modifies intensity measures for softer site conditions. We investigate whether these models, when used in tandem, are compatible with ground-motion recordings in central and eastern North America (CENA) using an expanded version of the NGA-East database that includes new events from November 2011 (end date of NGA-East data curation) to April 2022. Following this expansion, the data set has 187 events, 2096 sites, and 16,272 three-component recordings, although the magnitude range remains limited (∼4 to 5.8). We compute residuals using 17 NGA-East GMMs and three data selection criteria that reflect within-CENA regional variations in ground-motion attributes. Mixed-effects regression of the residuals reveals a persistent pattern in which ground motions are overpredicted at short periods (0.01–0.6 s, including peak ground acceleration (PGA)) and underpredicted at longer periods. These misfits are regionally variable, with the Texas–Oklahoma–Kansas region having larger absolute misfits than other parts of CENA. Two factors potentially influencing these misfits are (1) differences in the site amplification models used to adjust the data to the reference condition during NGA-East GMM development relative to CENA amplification models applied since the 2018 National Seismic Hazard Model (NSHM), and (2) potential bias in simulation-based factors used to adjust ground motions from the hard-rock reference condition to a V S30 = 760 m/s condition. We provide adjustment factors and their epistemic uncertainties and discuss implications for applications.

Kansas, Oklahoma, Texas↗
USGS2013

Site Response and Basin Waves in the Sacramento–San Joaquin Delta, California

The Sacramento–San Joaquin Delta is an inland delta at the western extent of the Central Valley. Levees were built around swampy islands starting after the Civil War to reclaim these lands for farming. Various studies show that these levees could fail in concert from shaking from a major local or regional earthquake resulting in salty water from the San Francisco Bay contaminating the water in the Delta. We installed seismographs around the Delta and on levees to assess the contribution of site response to the seismic hazard of the levees. Cone penetrometer testing shows that the upper 10 s of meters of soil in the Delta have shear‐wave velocities of about 200 m/s, which would give a strong site response. Seismographs were sited following two strategies: pairs of stations to compare the response of the levees to nearby sites, and a more regional deployment in the Delta. Site response was determined in two different ways: a traditional spectral ratio (TSR) approach of S waves using station BDM of the Berkeley Digital Seismic Net as a reference site, and using SH / SV ratios of noise (or Nakamura’s method). Both estimates usually agree in spectral character for stations whose response is dominated by a resonant peak, but the most obvious peaks in the SH / SV ratios usually are about two‐thirds as large as the main peaks in the TSRs. Levee sites typically have large narrow resonances in the site response function compared to sites in the farmland of the Delta. These resonances, at a frequency of about 1–3 Hz, have amplitudes of about 15 with TSR and 10–12 with Nakamura’s method. Sites on farmland in the Delta also have amplifications, but these are typically broader and not as resonant in appearance. Late (slow) Rayleigh waves were recorded at stations in the Delta, have a dominant period of about one second, and are highly monochromatic. Results from a three‐station array at the Holland Marina suggest that they have a phase velocity of about 600 m/s and arrive at about the same azimuth as the straight‐line back azimuth to the source. A dispersion curve determined for the basin or valley waves yields a shallow velocity profile that increases from about 350 m/s in the upper 0.2 km to about 1.1 km/s at a depth of about 2 km.

California↗
USGS1981

Northwest margin of California continental borderland: Marine geology and tectonic evolution

The northwest margin of the California continental borderland consists mainly of two northwest-trending pre-Neogene lithologic belts blanketed by Miocene and younger strata. These belts, which are lithologically and structurally correlated with the Franciscan Complex and Great Valley sequence of northern California, are interpreted to represent facies corresponding to the subduction complex and forearc-basin deposits of a late Mesozoic and Paleogene continental-margin arc-trench system. The outer belt, which forms the acoustic basement underlying the continental slope and Patton Ridge, is characterized by: (1) moderately high compressional velocities (5.1 to 6.2 km/sec), (2) discordant and discontinuous seismic reflectors, and (3) nonfossiliferous and diverse rock types t at range in metamorphic grade from zeolite-bearing arkosic sandstone to blueschist. Landward, the inner Great Valley sequence belt underlies Santa Rosa-Cortes Ridge and Santa Cruz and San Nicolas Basins; in contrast to the outer belt, this belt is characterized by: (1) intermediate compressional velocities (4.5 to 4.6 km/sec), (2) concordant and relatively continuous reflectors, and (3) a thick turbidite sequence of relatively undeformed and unmetamorphosed fossiliferous Cretaceous and Paleogene strata. The inner belt overlies a 6.6-km/sec layer that is probably oceanic crust, a relation consistent with the hypothesis that trapped old oceanic crust underlies the forearc basins along parts of some convergent margins. Well-developed wrench-fault structures in overlying Miocene strata record a transition from Paleogene subduction to Neogene transform faulting. Related events include widespread andesitic to dacitic volcanism and concurrent uplift and erosion of Miocene and older rocks. The proximity of these Miocene silicic volcanic rocks to the paleotrench (< 20 km) is best explained as resulting from passage of the southward-migrating Rivera ridge-fault-trench triple junction along the margin. On the basis of the age of the adjacent sea floor, passage of this triple junction occurred about 17 m.y. ago, a date that corresponds closely to the inception of volcanism and uplift in the northwest margin of the borderland. A widespread unconformity that separates highly deformed lower and lower middle iocene strata from relatively undeformed younger strata is assumed to mark the subsequent onset of wrench faulting. The timing of volcanism and uplift, and of the slightly younger wrenching event, closely matches the plate-tectonic model of Atwater and Molnar if the northwest margin of the borderland is restored to an earlier tectonic setting off Baja California, approximately 260 km southeast of its present position.

California↗
USGS1980

Multichannel seismic-reflection profiles collected along the U.S. continental margin in 1978

During 1978, the U.S.. Geological Survey (USGS) contracted with Geophysical Services, Inc. (GSI) for GSI to collect 4,813 km of 48-channel seismic-reflection profiles along the continental margin between North Carolina and Maine. Fifteen lines were acquired perpendicular to the margin Clines 18 through 32), and six lines were acquired parallel to the margin Clines 33 through 38). The profiles were shot by the Motor Vessel (MV) CARINO and MV CECIL GREEN between May and November 1978. A 3,600-long hydrophone streamer was used; it was composed of 24 groups of 100-m length each near the ship, followed by 24 groups of 50-mlength each. Tuned airgun arrays totaling 2,000 and 1,400 cubic inches of air were used by the MV- CARINO and MV CECIL GREEN, respectively. The profiles were stacked by GSI and displayed in conventional time-variant scaled format to 12 seconds of recording time with vertical scales of .2.5 inches per second and horizontal scales of 1.2 km per inch. Additional protessing was applied to Line 25 across the Baltimore Canyon Trough and Line 32 across the Carolina Trough. Data collected along Line 25 over the Outer Shelf, Slope and upper Rise were restacked with 1 1/2-km-spaced velocity analyses, signature deconvolution, deep-water multiple deconvolution, velocity filtering Con shelf only), and time migration. The restacked part and remaining parts of Line 25 were converted to depth sections having a vertical scale of 1.2 km per inch (vertical exaggeration: 2/1). Data collected along two parts of Line 32 were also restacked with signature and deep-water multiple deconvolution. Line 32 was also converted to a depth section except for the first 80 km near shore. Profiles 18-38 may be inspected at the USGS in Woods Hole, MA 02543. Copies of the profiles may be, purchased only from the National Geophysical and Solar-Terrestrial Data Center, NOAA, Boulder, Co 80303.

Open-File Report↗
USGS2017

Dynamic strains for earthquake source characterization

Strainmeters measure elastodynamic deformation associated with earthquakes over a broad frequency band, with detection characteristics that complement traditional instrumentation, but they are commonly used to study slow transient deformation along active faults and at subduction zones, for example. Here, we analyze dynamic strains at Plate Boundary Observatory (PBO) borehole strainmeters (BSM) associated with 146 local and regional earthquakes from 2004–2014, with magnitudes from M 4.5 to 7.2. We find that peak values in seismic strain can be predicted from a general regression against distance and magnitude, with improvements in accuracy gained by accounting for biases associated with site–station effects and source–path effects, the latter exhibiting the strongest influence on the regression coefficients. To account for the influence of these biases in a general way, we include crustal‐type classifications from the CRUST1.0 global velocity model, which demonstrates that high‐frequency strain data from the PBO BSM network carry information on crustal structure and fault mechanics: earthquakes nucleating offshore on the Blanco fracture zone, for example, generate consistently lower dynamic strains than earthquakes around the Sierra Nevada microplate and in the Salton trough. Finally, we test our dynamic strain prediction equations on the 2011 M 9 Tohoku‐Oki earthquake, specifically continuous strain records derived from triangulation of 137 high‐rate Global Navigation Satellite System Earth Observation Network stations in Japan. Moment magnitudes inferred from these data and the strain model are in agreement when Global Positioning System subnetworks are unaffected by spatial aliasing.

Seismological Research Letters↗
USGS2008

The ShakeOut Scenario

This is the initial publication of the results of a cooperative project to examine the implications of a major earthquake in southern California. The study comprised eight counties: Imperial, Kern, Los Angeles, Orange, Riverside, San Bernardino, San Diego, and Ventura. Its results will be used as the basis of an emergency response and preparedness exercise, the Great Southern California ShakeOut, and for this purpose we defined our earthquake as occurring at 10:00 a.m. on November 13, 2008. As members of the southern California community use the ShakeOut Scenario to plan and execute the exercise, we anticipate discussion and feedback. This community input will be used to refine our assessment and will lead to a formal publication in early 2009. Our goal in the ShakeOut Scenario is to identify the physical, social and economic consequences of a major earthquake in southern California and in so doing, enable the users of our results to identify what they can change now?before the earthquake?to avoid catastrophic impact after the inevitable earthquake occurs. To do so, we had to determine the physical damages (casualties and losses) caused by the earthquake and the impact of those damages on the region?s social and economic systems. To do this, we needed to know about the earthquake ground shaking and fault rupture. So we first constructed an earthquake, taking all available earthquake research information, from trenching and exposed evidence of prehistoric earthquakes, to analysis of instrumental recordings of large earthquakes and the latest theory in earthquake source physics. We modeled a magnitude (M) 7.8 earthquake on the southern San Andreas Fault, a plausible event on the fault most likely to produce a major earthquake. This information was then fed forward into the rest of the ShakeOut Scenario. The damage impacts of the scenario earthquake were estimated using both HAZUS-MH and expert opinion through 13 special studies and 6 expert panels, and fall into four categories: building damages, non-structural damages, damage to lifelines and infrastructure, and fire losses. The magnitude 7.8 ShakeOut earthquake is modeled to cause about 1800 deaths and $213 billion of economic losses. These numbers are as low as they are because of aggressive retrofitting programs that have increased the seismic resistance of buildings, highways and lifelines, and economic resiliency. These numbers are as large as they are because much more retrofitting could still be done. The earthquake modeled here may never happen. Big earthquakes on the San Andreas Fault are inevitable, and by geologic standards extremely common, but probably will not be exactly like this one. The next very damaging earthquake could easily be on another fault. However, lessons learned from this particular event apply to many other events and could provide benefits in many possible future events.

Open-File Report↗
USGS1980

Geophysical observations of Kilauea Volcano, Hawaii, 2. Constraints on the magma supply during November 1975-September 1977

Following a 22-month hiatus in eruptive activity, Kilauea volcano extruded roughly 35 ?? 106 m3 of tholeiitic basalt from vents along its middle east rift zone during 13 September-1 October, 1977. The lengthy prelude to this eruption began with a magnitude 7.2 earthquake on 29 November, 1975, and included rapid summit deflation episodes in June, July, and August 1976 and February 1977. Synthesis of seismic, geodetic, gravimetric, and electrical self-potential observations suggests the following model for this atypical Kilauea eruptive cycle. Rapid summit deflation initiated by the November 1975 earthquake reflected substantial migration of magma from beneath the summit region of Kilauea into the east and southwest rift zones. Simultaneous leveling and microgravity observations suggest that 40-90 ?? 106 m3 of void space was created within the summit magma chamber as a result of the earthquake. If this volume was filled by magma from depth before the east rift zone intrusive event of June 1976, the average rate of supply was 6-13 ?? 106 m3/month, a rate that is consistent with the value of 9 ?? 106 m3/month suggested from observations of long-duration Kilauea eruptions. Essentially zero net vertical change was recorded at the summit during the 15-month period beginning with the June 1976 intrusion and ending with the September 1977 eruption. This fact suggests that most magma supplied from depth during this interval was eventually delivered to the east rift zone, at least in part during four rapid summit deflation episodes. Microearthquake epicenters migrated downrift to the middle east rift zone for the first time during the later stages of the February 1977 intrusion, an occurrence presumably reflecting movement of magma into the eventual eruptive zone. This observation was confirmed by tilt surveys in May 1977 that revealed a major inflation center roughly 30 km east of the summit in an area of anomalous steaming and forest kill first noted in March 1976. ?? 1980.

Journal of Volcanology and Geothermal Research↗
USGS2018

Role of fault gouge during Interaction between hydraulic fracture and a preexisting fracture

Enhanced reservoir connectivity generally requires maximizing the intersection between hydraulic fracture (HF) and preexisting underground natural fractures (NF), while having the hydraulic fracture continue to propagate across the natural fractures. Observations of downhole core samples suggest that these natural fractures are in fact veins filled with minerals such as calcite (Mighani et al., 2016). We study this interaction during the approach of a hydraulic fracture to a smooth saw-cut fracture under triaxial stress conditions. The specimen is Solnhofen limestone, a fine-grained (<5 µm grain), low permeability (<10 nD) carbonate. The differential stress (1-20 MPa) and inclination of the fault which determines the approach angle, θ (30, 60°) are the experimental variables. We conduct the experiments on both bare surface and gouge-filled fault surfaces. The gouge is a 1 mm thick crushed powder of Solnhofen limestone with <106 µm grain size. During the hydraulic fracture, acoustic emissions (AE), inferred slip velocity, axial stress and pore pressure are recorded at a 5 MHz sampling rate. The hydraulic fracture was able to cross the bare surface fault with small induced fault slip. The fault gouge increased the coefficient of friction significantly from 0.12 (bare, polished surface) to > 0.44 (gouged layer). However, the gouge-filled fault arrested the hydraulic fracture and generated a slip event with different characteristics: 1- The stress drop was larger while the generated AE signals had lower magnitude. 2- Slip velocity recorded by the vibrometer was of the same order of magnitude for the bare and gouge-filled faults, but the slip duration increased from 29 µsec for bare surface to ~2.5 msec (~90 times longer rise time) for the gouge-filled fault. The experiments suggest that the gouge-filled fault can accommodate much larger displacement while promoting slow slip on the fault which is harder to detect as AE signals. The observed long duration slip events are similar to the field observations of the long period and long duration (LPLD) events during the stimulation of clay-rich shale formations (Zoback et al., 2012). While the intrinsic low strength, high ductility, and unfavorably oriented natural fractures in shale formations are expected to reduce the occurrence of induced seismicity, our experiments suggest an additional mechanism for the observed LPLD events, i.e. the role of fault gouge. They also suggest that the microseismic detection techniques may under-predict the stimulated volume as the activation of natural gouge-filled fractures may proceed aseismically.

Conference Paper↗
USGS1997

Review of electric and magnetic fields accompanying seismic and volcanic activity

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.

Surveys in Geophysics↗
USGS2025

Introduction to the special section on improving measurements of earthquake source parameters

Earthquake source parameters such as magnitude, seismic moment, source dimension, stress drop, and radiated energy are fundamental to understanding earthquake physics, and are also key ingredients in earthquake ground‐motion modeling, rupture simulation, and statistical seismology. However, the uncertainties in these parameters estimated from the radiated seismic wavefield are large due to variability in approaches, including site and attenuation characterizations, and so estimates for an individual earthquake made by different studies can vary greatly. Estimating spectral source parameters remains a popular topic, due to a combination of their intrinsic importance and their apparent ease of measurement, but also a controversial one, due to many sources of variability and large uncertainty. The available methods coupled with necessary parameter choices and assumptions in the analysis make it challenging for researchers to apply methods or understand the reliability in results or reported source parameters. This Special Section on Improving Measurements of Earthquake Source Parameters showcasing comparisons between methods and studies seeks to alleviate some of these difficulties to help the community identify the important components and trade‐offs of decomposing recorded seismograms into their source, path, and site components.

Bulletin of the Seismological Society of America↗
USGS1994

Slip triggered on southern California faults by the 1992 Joshua Tree, Landers, and big bear earthquakes

Five out of six functioning creepmeters on southern California faults recorded slip triggered at the time of some or all of the three largest events of the 1992 Landers earthquake sequence. Digital creep data indicate that dextral slip was triggered within 1 min of each mainshock and that maximum slip velocities occurred 2 to 3 min later. The duration of triggered slip events ranged from a few hours to several weeks. We note that triggered slip occurs commonly on faults that exhibit fault creep. To account for the observation that slip can be triggered repeatedly on a fault, we propose that the amplitude of triggered slip may be proportional to the depth of slip in the creep event and to the available near-surface tectonic strain that would otherwise eventually be released as fault creep. We advance the notion that seismic surface waves, perhaps amplified by sediments, generate transient local conditions that favor the release of tectonic strain to varying depths. Synthetic strain seismograms are presented that suggest increased pore pressure during periods of fault-normal contraction may be responsible for triggered slip, since maximum dextral shear strain transients correspond to times of maximum fault-normal contraction.

California↗
USGS1981

Landslides from the February 4, 1976, Guatemala earthquake

The M (Richter magnitude) = 7.5 Guatemala earthquake of February 4, 1976, generated more than 10,000 landslides throughout an area of approximately 16,000 km2. These landslides caused hundreds of fatalities as well as extensive property damage. Landslides disrupted both highways and the railroad system and thus severely hindered early rescue efforts. In Guatemala City, extensive property damage and loss of life were due to ground failure beneath dwellings built too close to the edges of steeply incised canyons. We have recorded the distribution of landslides from this earthquake by mapping individual slides at a scale of 1:50,000 for most of the landslide-affected area, using high-altitude aerial photography. The highest density of landslides was in the highlands west of Guatemala City. The predominant types of earthquake-triggered landslides were rock falls and debris slides of less than 15,000 m3 volume; in addition to these smaller landslides, 11 large landslides had volumes of more than 100,000 m3. Several of these large landslides posed special hazards to people and property from lakes impounded by the landslide debris and from the ensuing floods that occurred upon breaching and rapid erosion of the debris. The regional landslide distribution was observed to depend on five major factors: (1) seismic intensity; (2) lithology: 90 percent of all landslides were within Pleistocene pumice deposits; (3) slope steepness; (4) topographic amplification of seismic ground motion; and (5) regional fractures. The presence of preearthquake landslides had no apparent effect on the landslide distribution, and landslide concentration in the Guatemala City area does not correlate with local seismic-intensity data. The landslide concentration, examined at this scale, appears to be governed mainly by lithologic differences within the pumice deposits, preexisting fractures, and amplification of ground motion by topography-all factors related to site conditions.

Professional Paper↗
USGS1992

The Loma Prieta, California, Earthquake of October 17, 1989: Strong ground motion and ground failure

Professional Paper 1551 describes the effects at the land surface caused by the Loma Prieta earthquake. These effects: include the pattern and characteristics of strong ground shaking, liquefaction of both floodplain deposits along the Pajaro and Salinas Rivers in the Monterey Bay region and sandy artificial fills along the margins of San Francisco Bay, landslides in the epicentral region, and increased stream flow. Some significant findings and their impacts were: * Strong shaking that was amplified by a factor of about two by soft soils caused damage at up to 100 kilometers (60 miles) from the epicenter. * Instrumental recordings of the ground shaking have been used to improve how building codes consider site amplification effects from soft soils. * Liquefaction at 134 locations caused $99.2 million of the total earthquake loss of $5.9 billion. Liquefaction of floodplain deposits and sandy artificial fills was similar in nature to that which occurred in the 1906 San Francisco earthquake and indicated that many areas remain susceptible to liquefaction damage in the San Francisco and Monterey Bay regions. * Landslides caused $30 million in earthquake losses, damaging at least 200 residences. Many landslides showed evidence of movement in previous earthquakes. * Recognition of the similarities between liquefaction and landslides in 1906 and 1989 and research in intervening years that established methodologies to map liquefaction and landslide hazards prompted the California legislature to pass in 1990 the Seismic Hazards Mapping Act that required the California Geological Survey to delineate regulatory zones of areas potentially susceptible to these hazards. * The earthquake caused the flow of many streams in the epicentral region to increase. Effects were noted up to 88 km from the epicenter. * Post-earthquake studies of the Marina District of San Francisco provide perhaps the most comprehensive case history of earthquake effects at a specific site developed for any earthquake. Soft soils beneath the Marina amplified ground shaking to damaging levels and caused liquefaction of sandy artificial fills. Liquefaction required 123 repairs of pipelines in the Municipal Water Supply System, more than three times the number of repairs elsewhere in the system. Approximately 13.6 km of gas-distribution lines were replaced, and more than 20% of the wastewater collection lines were repaired or replaced.

Professional Paper↗
USGS2026

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii↗
USGS2001

Leg 188 synthesis: Transitions in the glacial history of the Prydz Bay region, East Antarctica, from ODP drilling

Drilling during Leg 119 (1988) and Leg 188 (2000; Sites 1165–1167) of the Ocean Drilling Program (ODP) provides direct evidence for long- and short-term changes in Cenozoic paleoenvironments in the Prydz Bay region. Cores from across the continental margin reveal that in preglacial times the present shelf was an alluvial plain system with austral conifer woodland in the Late Cretaceous that changed to cooler Nothofagus rainforest scrub by the middle to late Eocene (Site 1166). Earliest recovered evidence of nearby mountain glaciation is seen in late Eocene–age grain textures in fluvial sands. In the late Eocene to early Oligocene, Prydz Bay permanently shifted from being a fluvio-deltaic complex to an exclusively marine continental shelf environment. This transition is marked by a marine flooding surface later covered by overcompacted glacial sediments that denote the first advance of the ice sheet onto the shelf. Cores do not exist for the early Oligocene to early Miocene, and seismic data are used to infer the transition from a shallow to normal depth prograding continental shelf with submarine canyons on the slope and channel/levees on the rise. Cores from the continental rise at Site 1165 show long-term (millions of years) early Miocene and younger decreases in sedimentation rates as well as short-term (Milankovitch periods) cyclicity between principally biogenic and terrigenous sediment supply—resulting from the cyclic presence of onshore glaciers and changes in ocean circulation. Middle Miocene transitions include rapid decreases in sedimentation rates, increased ice-rafted debris, shifts in clays and other minerals, and regional erosion of the slope and rise. These transitions may reflect enhanced glacial erosion and reduced glacial meltwater from progressively colder ice. At this time, seismic data show that depocenters began to shift from the outer continental rise to the base of the continental slope coincident with the initial stages of the glacial erosion and overdeepening of the continental shelf. During the late Miocene to early Pliocene there was a transition to greater subglacial activity on the shelf and more pronounced cyclic facies variations on the continental rise. At this time, severe glacial morphologies initiated on the shelf with the erosion of Prydz Channel and other troughs by fast-moving ice and the deposition of overcompacted glacial diamictons by slow-moving ice on adjacent banks. The Prydz trough-mouth fan also began to form with alternating deposition of debris flows (ice at shelf edge) and muddy units (reduced ice) (Site 1167). The fan also records a transition during the late Pleistocene for times younger than 780 k.y. when short-term glacial variations continued but ice reached the shelf edge only a few times. Both short-term and long-term transitions characterize the Cenozoic evolution of the Prydz Bay region from the Cretaceous nonglacial to late Neogene full-glacial paleoenvironments. These transitions are known only from ODP cores, and further insights will require additional drilling.

Prydz Bay region↗
USGS1995

Finite-fault analysis of the 1979 March 14 Petatlan, Mexico, earthquake using teleseismic P waveforms

Vertical, teleseismic P waves recorded for the 1979 March 14 Petatlan, Mexico, earthquake were used to derive the distribution of coseismic slip using a linear finite-fault inversion scheme that solves for the amount of slip in each of a series of consecutive time windows. Data recorded by six stations of the Global Digital Seismograph Network were inverted in addition to digitized analogue long-period recordings available from nine Worldwide Standardized Seismograph Network stations. The digital data include four broad-band and short-period velocity waveforms reconstructed from the short- and long-period components. The time-window approach allows for a variable rise time on the fault and accounts for the source multiplicity evident in the recorded P waveforms. Synthetic tests conducted using the inversion method on the limited data set, however, reveal that the data are insufficient to identify the exact dislocation duration on the fault. The method is thus implemented by prescribing the fault rise time using five consecutive 1 s time windows. The coseismic slip inferred from the P waves shows a small 70 cm peak near the earthquake hypocentre and a large zone of dislocation (1.2 m maximum) further south-east. The slip pattern covers depths from 3 to 25 km and is located south-east of other recent large interplate ruptures on the Michoacan segment of the Mexican subduction zone. This result indicates that the 1979 Petatlan earthquake broke an independent, adjacent portion of the Cocos-North America plate boundary. The seismic moment of 1.5 × 10 27 dyn cm inferred from the P waves is approximately one-half the long-period moment estimated by other investigators from the observed surface waves. Although the discrepancy is within the uncertainty of the seismic-moment estimates, it may suggest the presence of a component of slow interplate motion that did not radiate significant P -wave energy.

Geophysical Journal International↗