Search USGS⌕ Search

SEARCH · Search USGS

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,369 records · Page 76Linked to original sources
USGS2021

Response of an asymmetrical five-story building in Fairbanks, Alaska during the November 30, 2018 M7.1 Anchorage, Alaska earthquake

A recently constructed, five-story, asymmetrical steel building on the campus of the University of Alaska, Fairbanks was equipped with a strong-motion array that recorded the M7.1 Anchorage earthquake of November 30, 2018 at an epicentral distance of 408 km. The largest recorded peak accelerations at the basement and top of the building are 0.021g and 0.071g, respectively. The steel building is designed with several bays that utilize K-shaped buckling restrained braces. The building response records allow identification of fundamental periods (frequencies) as 0.73s (1.4 Hz), 0.63s (1.60 Hz), and 0.56s (1.78 Hz) in the NS, EW, and torsional directions, respectively. System identification computations resulted in estimated critical damping percentages as 7.7% and 4.7 % in the NS and EW directions, respectively. At this low-level of shaking, the building is not expected to (and did not) experience observable damage, which is confirmed with very small average drift ratios. This is the first time a seismic response from this structural array has been analyzed.

Alaska↗
USGS1995

A mechanism for decoupling within the oceanic lithosphere revealed in the Troodos ophiolite

Contrasting kinematic histories recorded in the sheeted dykes and underlying plutonic rocks of the Troodos ophiolite provide a new perspective on the mechanical evolution of oceanic spreading centres. The kinematic framework of the decoupling zone that partitions deformation between the sheeted dykes and plutonics contrasts with low-angle detachment models for slow-spreading ridges based on continental-rift analogues. A model for the generation of multiple, horizontal decoupling horizons, linked by planar normal faults, demonstrates new possibilities for the kinematic and rheological significance of seismic reflectors in oceanic lithosphere.

Nature↗
USGS2001

Geophysics: Silent slip on the Cascadia subduction interface

G eodetic satellites of the Global Positioning System (GPS) now permit continuous recording of surface motions around earthquake faults and volcanoes with millimeter precision. Data from distinct points on Earth's surface can be combined to infer the locations of the sources of deformation at depth in the crust. These data are radically altering our understanding of earthquake processes, justifying ambitious new sensor arrays to image active deformation sources in Earth's crust. In recent years, large continuous GPS arrays have been deployed in Japan and southern California, and smaller networks have been installed in other seismically active regions. On page 1525 of this issue ( 1 ), Dragert et al. provide a glimpse of the kinds of insights we can expect from these arrays.

Cascadia subduction zone↗
USGS2007

Shallow P- and S-wave velocities and site resonances in the St. Louis region, Missouri-Illinois

As part of the seismic hazard-mapping efforts in the St. Louis metropolitan area we determined the compressional and shear-wave velocities (Vp and Vs) to about a 40-m depth at 17 locations in this area. The Vs measurements were made using high-resolution seismic refraction and reflection methods. We find a clear difference in the Vs profiles between sites located on the river floodplains and those located in the upland urban areas of St. Louis. Vs30 (average Vs to 30-m depth) values in floodplain areas range from 200 to 290 m/s (NEHRP category D) and contrast with sites on the upland areas of St. Louis, which have Vs30 values ranging from 410 to 785 m/s (NEHRP categories C and B). The lower Vs30 values and earthquake recordings in the floodplains suggest a greater potential for stronger and more prolonged ground shaking in an earthquake. Spectral analysis of a M3.6 earthquake recorded on the St. Louis-area ANSS seismograph network indicates stronger shaking and potentially damaging S-wave resonant frequencies at NEHRP category D sites compared to ground motions at a rock site located on the Saint Louis University campus. ?? 2007, Earthquake Engineering Research Institute.

Earthquake Spectra↗
USGS2004

Catalog of earthquake hypocenters at Alaskan volcanoes: January 1 through December 31, 2003

The Alaska Volcano Observatory (AVO), a cooperative program of the U.S. Geological Survey, the Geophysical Institute of the University of Alaska Fairbanks, and the Alaska Division of Geological and Geophysical Surveys, has maintained seismic monitoring networks at historically active volcanoes in Alaska since 1988. The primary objectives of this program are the near real time seismic monitoring of active, potentially hazardous, Alaskan volcanoes and the investigation of seismic processes associated with active volcanism. This catalog presents the calculated earthquake hypocenter and phase arrival data, and changes in the seismic monitoring program for the period January 1 through December 31, 2003. The AVO seismograph network was used to monitor the seismic activity at twenty-seven volcanoes within Alaska in 2003. These include Mount Wrangell, Mount Spurr, Redoubt Volcano, Iliamna Volcano, Augustine Volcano, Katmai volcanic cluster (Snowy Mountain, Mount Griggs, Mount Katmai, Novarupta, Trident Volcano, Mount Mageik, Mount Martin), Aniakchak Crater, Mount Veniaminof, Pavlof Volcano, Mount Dutton, Isanotski Peaks, Shishaldin Volcano, Fisher Caldera, Westdahl Peak, Akutan Peak, Makushin Volcano, Okmok Caldera, Great Sitkin Volcano, Kanaga Volcano, Tanaga Volcano, and Mount Gareloi. Monitoring highlights in 2003 include: continuing elevated seismicity at Mount Veniaminof in January-April (volcanic unrest began in August 2002), volcanogenic seismic swarms at Shishaldin Volcano throughout the year, and low-level tremor at Okmok Caldera throughout the year. Instrumentation and data acquisition highlights in 2003 were the installation of subnetworks on Tanaga and Gareloi Islands, the installation of broadband installations on Akutan Volcano and Okmok Caldera, and the establishment of telemetry for the Okmok Caldera subnetwork. AVO located 3911 earthquakes in 2003. This catalog includes: (1) a description of instruments deployed in the field and their locations; (2) a description of earthquake detection, recording, analysis, and data archival systems; (3) a description of velocity models used for earthquake locations; (4) a summary of earthquakes located in 2003; and (5) an accompanying UNIX tar-file with a summary of earthquake origin times, hypocenters, magnitudes, phase arrival times, and location quality statistics; daily station usage statistics; and all HYPOELLIPSE files used to determine the earthquake locations in 2003.

Alaska↗
USGS2026

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗
USGS2017

Geologic field-trip guide to the volcanic and hydrothermal landscape of the Yellowstone Plateau

Yellowstone National Park, a nearly 9,000 km 2 (~3,468 mi 2 ) area, was preserved in 1872 as the world’s first national park for its unique, extraordinary, and magnificent natural features. Rimmed by a crescent of older mountainous terrain, Yellowstone National Park has at its core the Quaternary Yellowstone Plateau, an undulating landscape shaped by forces of late Cenozoic explosive and effusive volcanism, on-going tectonism, glaciation, and hydrothermal activity. The Yellowstone Caldera is the centerpiece of the Yellowstone Plateau. The Yellowstone Plateau lies at the most northeastern front of the 17-Ma Yellowstone hot spot track, one of the few places on Earth where time-transgressive processes on continental crust can be observed in the volcanic and tectonic (faulting and uplift) record at the rate and direction predicted by plate motion. Over six days, this field trip presents an intensive overview into volcanism, tectonism, and hydrothermal activity on the Yellowstone Plateau (fig. 1). Field stops are linked directly to conceptual models related to monitoring of the various volcanic, geochemical, hydrothermal, and tectonic aspects of the greater Yellowstone system. Recent interest in young and possible future volcanism at Yellowstone as well as new discoveries and synthesis of previous studies, (for example, tomographic, deformation, gas, aeromagnetic, bathymetric, and seismic surveys), provide a framework in which to discuss volcanic, hydrothermal, and seismic activity in this dynamic region.

Yellowstone Plateau↗
USGS2014

Holocene and latest Pleistocene paleoseismology of the Salt Lake City segment of the Wasatch Fault Zone, Utah, at the Penrose Drive Trench Site

The Salt Lake City segment (SLCS) of the Wasatch fault zone (WFZ) and the West Valley fault zone (WVFZ) compromise Holocene-active normal faults that bound a large intrabasin graben in northern Salt Lake Valley and have evidence of recurrent, large-magnitude (M ~6-7) surface-faulting earthquakes. However, at the time of this investigation, questions remained regarding the timing, displacement, and recurrence of latest Pleistocene and Holocene earthquakes on the northern SLCS and WVFZ , and whether the WVFZ is seismically independent of, or moves coseismically with, the SLCS. To improve paleoseismic data for the SLCS, we conducted a fault-trench investigation at the Penrose Drive site on the northern SLCS. Two trenches, excavated across an 11-m-high scarp near the northern end of the East Bench fault, exposed colluvial-wedge evidence for fize of six (preferred) surface-faulting earthquakes postdating to Provo-phase shoreline of Lake Bonneville (~14-18 ka). Radiocarbon and luminescence ages support earthquake times at 4.0 ± 0.5 ka (2σ) (PD1), 5.9 ± 0.7 ka (PD2), 7.5 ± 0.8 ka (PD3a), 9.7 ± 1.1 ka (PD3b), 10.9 ± 0.2 ka (PD4), and 12.1 ± 1.6 ka (PD5). At least one additional earthquake occurred at 16.5 ± 1.9 ka (PD6) based on an erosional unconformity that separates deformed Lake Bonneville sily and flat-lying Provo-phase shoreline gravel. Earthquakes PD5-PD1 yield latest Pleistocene (post-Provo) and Holocene mean recurrence intervals of ~1.6 kyr and ~1.7-1.9 kyr, respectively. Using 1.0-1.4 m of per-event vertical displacement for PD5-PD3b corroborate previously identified SLCS earthquakes at 4-10 ka. PD4 and PD5 occurred within an ~8-kyr *17-9 ka) time interval on the SLCS previously interpreted as a period of seismic quiescence, and PD6 possibly corresponds with a previously identified earthquake at ~17 ka (although both events have large timing uncertainties). The Penrose data, when combined with previous paleoseismic results, improve the latest Pleistocene-Holocene earthquake chronology of the SLCS, and demonstrate that the SLCS has been a consistently active source of large-magnitude earthquakes since the latest Pleistocene. At least nine surface-faulting earthquakes (S1-S9) have occurred since the highstand of Lake Bonneville (~18 ka). Where the SLCS earthquake record is most complete (since ~14 ka), per-site estimates of mean recurrence are similar for the latest Pleistocene (post-Provo) (~1.6 kyr), Holocene (~1.6-1.9 kyr), and late Holocene (~1.2-1.4 kyr). These SLCS paleoearthquake data indicate an essentially stable rate of earthquake recurrence since the latest Pleistocene and are important for understanding the earthquake potential of the SLCS, clarifying the seismogenic relation between the SLCS and WVFZ, and forecasting the probabilities of future large-magnitude earthquake in the Wasatch Front region.

Utah↗
USGS2004

Radiated energy and the rupture process of the Denali fault earthquake sequence of 2002 from broadband teleseismic body waves

Displacement, velocity, and velocity-squared records of P and SH body waves recorded at teleseismic distances are analyzed to determine the rupture characteristics of the Denali fault, Alaska, earthquake of 3 November 2002 ( M W 7.9, M e 8.1 ). Three episodes of rupture can be identified from broadband ( ∼0.1–5.0 Hz ) waveforms. The Denali fault earthquake started as a M W 7.3 thrust event. Subsequent right-lateral strike-slip rupture events with centroid depths of 9 km occurred about 22 and 49 sec later. The teleseismic P waves are dominated by energy at intermediate frequencies (0.1-1 Hz) radiated by the thrust event, while the SH waves are dominated by energy at lower frequencies (0.05-0.2 Hz) radiated by the strike-slip events. The strike-slip events exhibit strong directivity in the teleseismic SH waves. Correcting the recorded P -wave acceleration spectra for the effect of the free surface yields an estimate of 2.8 × 10 15 N m for the energy radiated by the thrust event. Correcting the recorded SH -wave acceleration spectra similarly yields an estimate of 3.3 × 10 16 N m for the energy radiated by the two strike-slip events. The average rupture velocity for the strike-slip rupture process is 1.1 β –1.2 β . The strike-slip events were located 90 and 188 km east of the epicenter. The rupture length over which significant or resolvable energy is radiated is, thus, far shorter than the 340-km fault length over which surface displacements were observed. However, the seismic moment released by these three events, 4 × 10 20 N m, was approximately half the seismic moment determined from very low-frequency analyses of the earthquake. The difference in seismic moment can be reasonably attributed to slip on fault segments that did not radiate significant or coherent seismic energy. These results suggest that very large and great strike-slip earthquakes can generate stress pulses that rapidly produce substantial slip with negligible stress drop and little discernible radiated energy on fault segments distant from the initial point of nucleation. The existence of this energy-deficient rupture mode has important implications for the evaluation of the seismic hazard of very large strike-slip earthquakes.

Alaska↗
USGS2012

Strain accumulation in western Washington

The Juan de Fuca plate is subducted beneath the North American plate off the coast of Washington at a rate of about 40 mm/yr N68°E. The average principal strain rates (extension reckoned positive) measured in northwestern Washington are as follows: Olympic peninsula 25 km south of Port Angeles from 1982 through 1990, and and near Seattle from 1972 through 1985, and . Both strain measurements are consistent with uniaxial contraction in the direction of plate convergence. Uplift rates inferred from tide gage recordings are about 4 mm/yr on the Pacific coast and near 0 mm/yr farther inland near Seattle. These deformation rates are consistent with a model of the Cascadia subduction zone in which the plate interface beneath the continental slope and outer continental shelf is locked but free to slip farther landward. The limited downdip extent of the locked segment of the plate interface is consistent with a shallow depth (∼20 km) of the isotherm (∼450°C) that defines the brittle‐ductile transition. Small thrust events diagnostic of seismic subduction should then occur only offshore and at shallow depths. The principal strain rates measured from 1972 through 1983 in the back arc region near Richland, Washington, are and .

Washington↗
USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗
USGS2005

Preliminary report on the 28 September 2004, M 6.0 Parkfield, California earthquake

The M w 6.0 Parkfield earthquake struck central California at 17:15:14 UTC on 28 September 2004. The epicenter was located 11 km southeast of the rural town of Parkfield, adjacent to Gold Hill and on the San Andreas Fault ( Figure 1 ). The California Integrated Seismic Network (CISN) reported that the hypocenter was located at 35.819°N, 120.364°W at a depth of 8.8 km. From the distribution of aftershocks and from models of seismograms, strain changes, and geodetic displacements from the earthquake, it appears that the rupture propagated to the northwest along the San Andreas Fault from its hypocenter beneath Gold Hill to Middle Mountain. Because of the earthquake's moderate size and the low population density, only minimal damage was reported, but strong ground motions of approximately 1 g were recorded at a few isolated points (Shakal et al. , 2005, this issue).

California↗
USGS2008

A geophysical investigation of shallow deformation along an anomalous section of the Wasatch fault zone, Utah, USA

We report the results of a geophysical study of the Wasatch fault zone near the Provo and Salt Lake City segment boundary. This area is anomalous because the fault zone strikes more east-west than north-south. Vibroseis was used to record a common mid-point (CMP) profile that provides information to depths of ???500 m. A tomographic velocity model, derived from first breaks, constrained source and receiver static corrections; this was required due to complex terrain and significant lateral velocity contrasts. The profile reveals an ???250-m-wide graben in the hanging wall of the main fault that is associated with both synthetic and antithetic faults. Faults defined by apparent reflector offsets propagate upward toward topographic gradients. Faults mapped from a nearby trench and the seismic profile also appear to correlate with topographic alignments on LiDAR gradient maps. The faults as measured in the trench show a wide range of apparent dips, 20??-90??, and appear to steepen with depth on the seismic section. Although the fault zone is likely composed of numerous small faults, the broad asymmetric structure in the hanging wall is fairly simple and dominated by two inward-facing ruptures. Our results indicate the feasibility of mapping fault zones in rugged terrain and complex near-surface geology using low-frequency vibroseis. Further, the integration of geologic mapping and seismic reflection can extend surface observations in areas where structural deformation is obscured by poorly stratified or otherwise unmappable deposits. Therefore, the vibroseis technique, when integrated with geological information, provides constraints for assessing geologic hazards in areas of potential development.

Conference Paper↗
USGS2014

Stratigraphic and microfossil evidence for a 4500-year history of Cascadia subduction zone earthquakes and tsunamis at Yaquina River estuary, Oregon, USA

The Sallys Bend swamp and marsh area on the central Oregon coast onshore of the Cascadia subduction zone contains a sequence of buried coastal wetland soils that extends back ∼4500 yr B.P. The upper 10 of the 12 soils are represented in multiple cores. Each soil is abruptly overlain by a sandy deposit and then, in most cases, by greater than 10 cm of mud. For eight of the 10 buried soils, times of soil burial are constrained through radiocarbon ages on fine, delicate detritus from the top of the buried soil; for two of the buried soils, diatom and foraminifera data constrain paleoenvironment at the time of soil burial. We infer that each buried soil represents a Cascadia subduction zone earthquake because the soils are laterally extensive and abruptly overlain by sandy deposits and mud. Preservation of coseismically buried soils occurred from 4500 yr ago until ∼500–600 yr ago, after which preservation was compromised by cessation of gradual relative sea-level rise, which in turn precluded drowning of marsh soils during instances of coseismic subsidence. Based on grain-size and microfossil data, sandy deposits overlying buried soils accumulated immediately after a subduction zone earthquake, during tsunami incursion into Sallys Bend. The possibility that the sandy deposits were sourced directly from landslides triggered upstream in the Yaquina River basin by seismic shaking was discounted based on sedimentologic, microfossil, and depositional site characteristics of the sandy deposits, which were inconsistent with a fluvial origin. Biostratigraphic analyses of sediment above two buried soils—in the case of two earthquakes, one occurring shortly after 1541–1708 cal. yr B.P. and the other occurring shortly after 3227–3444 cal. yr B.P.—provide estimates that coseismic subsidence was a minimum of 0.4 m. The average recurrence interval of subduction zone earthquakes is 420–580 yr, based on an ∼3750–4050-yr-long record and seven to nine interearthquake intervals. The comparison of the Yaquina Bay earthquake record to similar records at other Cascadia coastal sites helps to define potential patterns of rupture for different earthquakes, although inherent uncertainty in dating precludes definitive statements about rupture length during earthquakes. We infer that in the first half of the last millennia, the northern Oregon part of the subduction zone had a different rupture history than the southern Oregon part of the subduction zone, and we also infer that at ca. 1.6 ka, two earthquakes closely spaced in time together ruptured a length of the megathrust that extends at least from southwestern Washington to southern Oregon.

Oregon↗
USGS2003

The Mw 7.7 Bhuj earthquake: Global lessons for earthquake hazard in intra-plate regions

The Mw 7.7 Bhuj earthquake occurred in the Kachchh District of the State of Gujarat, India on 26 January 2001, and was one of the most damaging intraplate earthquakes ever recorded. This earthquake is in many ways similar to the three great New Madrid earthquakes that occurred in the central United States in 1811-1812, An Indo-US team is studying the similarities and differences of these sequences in order to learn lessons for earthquake hazard in intraplate regions. Herein we present some preliminary conclusions from that study. Both the Kutch and New Madrid regions have rift type geotectonic setting. In both regions the strain rates are of the order of 10-9/yr and attenuation of seismic waves as inferred from observations of intensity and liquefaction are low. These strain rates predict recurrence intervals for Bhuj or New Madrid sized earthquakes of several thousand years or more. In contrast, intervals estimated from paleoseismic studies and from other independent data are significantly shorter, probably hundreds of years. All these observations together may suggest that earthquakes relax high ambient stresses that are locally concentrated by rheologic heterogeneities, rather than loading by plate-tectonic forces. The latter model generally underlies basic assumptions made in earthquake hazard assessment, that the long-term average rate of energy released by earthquakes is determined by the tectonic loading rate, which thus implies an inherent average periodicity of earthquake occurrence. Interpreting the observations in terms of the former model therefore may require re-examining the basic assumptions of hazard assessment.

Gujarat↗
USGS2018

A collection of historic seismic instrumentation photographs at the Albuquerque Seismological Laboratory

The Albuquerque Seismological Laboratory (ASL) has preserved a collection of photographs of seismographic equipment, stations, and drawings used by the United States Coast and Geodetic Survey (USC&GS) in the early-to-mid-twentieth century. The photographs were transferred to ASL from the US Department of Commerce building in Washington DC after ASL became established as a USC&GS facility for seismological research and instrumentation testing in 1961. The photographs were variously used in publications to document USC&GS activities or used internally to document instrument installation, operation, and/or components. The historic photograph collection is made up of hundreds of photographs that were originally printed in the 1930s and 1940s. Our hope is that the archival of these records, as well as brief descriptions about their significance, will help preserve some of the knowledge about seismological advances in the United States during this time period. We provide digital scans of the complete ASL historic photo collection and discuss the historical context of select images in the Electronic Supplement of this document.

Seismological Research Letters↗
USGS2005

An efficient algorithm for double-difference tomography and location in heterogeneous media, with an application to the Kilauea volcano

Improving our understanding of crustal processes requires a better knowledge of the geometry and the position of geological bodies. In this study we have designed a method based upon double-difference relocation and tomography to image, as accurately as possible, a heterogeneous medium containing seismogenic objects. Our approach consisted not only of incorporating double difference in tomography but also partly in revisiting tomographic schemes for choosing accurate and stable numerical strategies, adapted to the use of cross-spectral time delays. We used a finite difference solution to the eikonal equation for travel time computation and a Tarantola-Valette approach for both the classical and double-difference three-dimensional tomographic inversion to find accurate earthquake locations and seismic velocity estimates. We estimated efficiently the square root of the inverse model's covariance matrix in the case of a Gaussian correlation function. It allows the use of correlation length and a priori model variance criteria to determine the optimal solution. Double-difference relocation of similar earthquakes is performed in the optimal velocity model, making absolute and relative locations less biased by the velocity model. Double-difference tomography is achieved by using high-accuracy time delay measurements. These algorithms have been applied to earthquake data recorded in the vicinity of Kilauea and Mauna Loa volcanoes for imaging the volcanic structures. Stable and detailed velocity models are obtained: the regional tomography unambiguously highlights the structure of the island of Hawaii and the double-difference tomography shows a detailed image of the southern Kilauea caldera-upper east rift zone magmatic complex. Copyright 2005 by the American Geophysical Union.

Hawaii↗
USGS2016

Compartmentalization of the Coso East Flank geothermal field imaged by 3-D full-tensor MT inversion

Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2–5 Ohm m) and appreciable lateral extent (>1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Z xx and Z yy data components leads to a robust east-dip (60†) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.

Coso Range↗