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 541 records · Page 30Linked to original sources

A great tsunami earthquake component of the 1957 Aleutian Islands earthquake

The great 1957 Aleutian Islands earthquake ruptured ∼1200 km of the plate boundary along the Aleutian subduction zone and produced a destructive tsunami across Hawaiʻi. Early seismic and tsunami analyses indicated that large megathrust fault slip was concentrated in the western Aleutian Islands, but tsunami waves generated by slip in the west cannot explain the large observed runup in Hawaiʻi far to the southeast. Recently mapped 1957 geologic deposits on eastern Aleutian Islands suggest occurrence of very large nearby slip. Jointly modeling tsunami runup along the eastern Aleutian and Hawaiian Islands together with tide gauge recordings across the Pacific resolves 12-26 m shallow slip along 600 km of the eastern Aleutian Islands in addition to modest, deeper western slip inferred from seismic records. The eastern near-trench slip results in an M W 8.3-8.6 tsunami earthquake component of the M W 8.6-8.8 rupture, comparable in size to the adjacent 1946 Aleutian tsunami earthquake to the east. The reexamination of the 1957 rupture confirms the tsunami hazards posed by the eastern Aleutian subduction zone to Hawaiʻi and lays the groundwork for investigation of large prehistoric earthquakes through modeling tsunami runup inferred from stratigraphic observations to constrain their rupture processes.

Alaska↗

Scattered surface wave energy in the seismic coda

One of the many important contributions that Aki has made to seismology pertains to the origin of coda waves (Aki, 1969; Aki and Chouet, 1975). In this paper, I revisit Aki's original idea of the role of scattered surface waves in the seismic coda. Based on the radiative transfer theory, I developed a new set of scattered wave energy equations by including scattered surface waves and body wave to surface wave scattering conversions. The work is an extended study of Zeng et al. (1991), Zeng (1993) and Sato (1994a) on multiple isotropic-scattering, and may shed new insight into the seismic coda wave interpretation. The scattering equations are solved numerically by first discretizing the model at regular grids and then solving the linear integral equations iteratively. The results show that scattered wave energy can be well approximated by body-wave to body wave scattering at earlier arrival times and short distances. At long distances from the source, scattered surface waves dominate scattered body waves at surface stations. Since surface waves are 2-D propagating waves, their scattered energies should in theory follow a common decay curve. The observed common decay trends on seismic coda of local earthquake recordings particular at long lapse times suggest that perhaps later seismic codas are dominated by scattered surface waves. When efficient body wave to surface wave conversion mechanisms are present in the shallow crustal layers, such as soft sediment layers, the scattered surface waves dominate the seismic coda at even early arrival times for shallow sources and at later arrival times for deeper events.

Pure and Applied Geophysics↗

Seismic evidence for a tilted mantle plume and north-south mantle flow beneath Iceland

Shear waves converted from compressional waves at mantle discontinuities near 410- and 660-km depth recorded by two broadband seismic experiments in Iceland reveal that the center of an area of anomalously thin mantle transition zone lies at least 100 km south of the upper-mantle low-velocity anomaly imaged tomographically beneath the hotspot. This offset is evidence for a tilted plume conduit in the upper mantle, the result of either northward flow of the Icelandic asthenosphere or southward flow of the upper part of the lower mantle in a no-net-rotation reference frame. ?? 2002 Elsevier Science B.V. All rights reserved.

Earth and Planetary Science Letters↗

Quantitative analysis of seismic fault zone waves in the rupture zone of the 1992 Landers, California, earthquake: Evidence for a shallow trapping structure

We analyse quantitatively a waveform data set of 238 earthquakes recorded by a dense seismic array across and along the rupture zone of the 1992 Landers earthquake. A grid-search method with station delay corrections is used to locate events that do not have catalogue locations. The quality of fault zone trapped waves generated by each event is determined from the ratios of seismic energy in time windows corresponding to trapped waves and direct S waves at stations close to and off the fault zone. Approximately 70 per cent of the events with S-P times of less than 2 s, including many clearly off the fault, produce considerable trapped wave energy. This distribution is in marked contrast with previous claims that trapped waves are generated only by sources close to or inside the Landers rupture zone. The time difference between the S arrival and trapped waves group does not grow systematically with increasing hypocentral distance and depth. The dispersion measured from the trapped waves is weak. These results imply that the seismic trapping structure at the Landers rupture zone is shallow and does not extend continuously along-strike by more than a few kilometres. Synthetic waveform modelling indicates that the fault zone waveguide has depth of approximately 2-4 km, a width of approximately 200 m, an S-wave velocity reduction relative to the host rock of approximately 30-40 per cent and an S-wave attenuation coefficient of approximately 20-30. The fault zone waveguide north of the array appears to be shallower and weaker than that south of the array. The waveform modelling also indicates that the seismic trapping structure below the array is centred approximately 100 m east of the surface break.

Geophysical Journal International↗

A seismic refraction survey of the Imperial Valley Region, California

The U.S. Geological Survey conducted an extensive seismic refraction survey in the Imperial Valley region of California in 1979. The Imperial Valley is located in the Salton Trough, an active rift between the Pacific and North American plates. Forty shots fired at seven shot points were recorded by 100 portable seismic instruments at typical spacing of 0.5–1 km. More than 1300 recording locations were occupied, and more than 3000 usable seismograms were obtained. We analyzed five profiles using a standard ray‐tracing program, constructed a contour map of reduced travel times from our most widely recorded shot point, and modeled an existing gravity profile across the Salton Trough. Results are itemized: (1) All models have in common a sedimentary layer ( V p = 1.8–5.0 km/s), a “transition zone” ( V p = 5.0–5.65 km/s), a basement ( V p = 5.65 km/s in the Imperial Valley, 5.9 km/s on the bordering mesas), and subbasement ( V p = 7.2 km/s). (2) The sedimentary layer ranges in thickness along the axis of the Salton Trough from 3.7 km (Salton Sea) to 4.8 km (U.S.‐Mexican border). On the bordering mesas it is quite variable in thickness. (3) The “transition” zone is about 1 km thick in most places. In the Imperial Valley there are no marked velocity discontinuities in this zone between the sedimentary layer and basement. On the bordering mesas, however, there is a discontinuity at the top of this zone. (4) There are apparently two types of basement. On the bordering mesas, basement is crystalline igneous and metamorphic rocks. In the Imperial Valley, basement is mostly lower‐greenshist‐facies sedimentary rocks, based primarily on the smooth transition in character from sediment to basement arrivals, the low value of basement velocity, and the fact that deep (4 km) wells in the valley penetrate only the upper part of the known Cenozoic stratigraphic column for the Salton Trough. (5) The subbasement, or intermediate crustal layer, ranges in depth along the axis of the Salton Trough from 16 km (Salton Sea) to 10 km (U.S.‐Mexican border). Gravity modeling requires that this layer deepen and/or pinch out beneath the bordering mesas and mountain ranges. Based on its high velocity and the presence of intrusive basaltic rocks in the sedimentary section in the Imperial Valley, the subbasement is thought to be a mafic intrusive complex similar to oceanic middle crust. (6) Several structures are seen that affect basement, transition zone, and deeper parts of the sedimentary layer. They include a scarp along the Imperial fault, as much as 1 km down to the northeast, and a scarp passing roughly along the topographic boundary between the Imperial Valley and the bordering mesa to the west, as much as 3½ km down to the east. We interpret the latter scarp to be the suture, or rift boundary, between the older crystalline basement on the mesa and the younger metasedimentary basement in the Imperial Valley. (7) On a contour map of reduced travel time from our most widely recorded shot point, subtle patches of early arrivals among otherwise late arrivals in the central Imperial Valley correlate well with known geothermal resource areas having reservoir temperatures of more than 150°C. Apparently the Salton Trough is a location where new crust is being generated. As the rift opens, mafic intrusive rocks fill it from below as sedimentary rocks fill it from above. Rifting and intrusion produce high heat flow that metamorphoses the sedimentary rocks to shallow depth (metasedimentary basement in the Imperial Valley) and thus consolidates the new crust.

California↗

Minimal clustering of injection-induced earthquakes observed with a large-n seismic array

The clustering behavior of injection‐induced earthquakes is examined using one month of data recorded by the LArge‐ n Seismic Survey in Oklahoma (LASSO) array. The 1829‐node seismic array was deployed in a 25 km×32 km area of active saltwater disposal in northern Oklahoma between 14 April and 10 May 2016. Injection rates in the study area are nearly constant around the time of the deployment. We develop a local magnitude ( ⁠ M L ⁠ ) equation for the region and estimate magnitudes for 1104 earthquakes recorded by the deployment. The determined earthquake magnitudes range from M L 0.01 to 3.0. The majority of earthquakes occurred between 1.5 and 5.5 km depth, and the shallowest earthquake depths overlap with the base of injection wells at depths between 1.5 and 2.5 km. We compute focal mechanisms of the largest events ( ⁠ M L > 2.0 ⁠ ), and find a mix of normal‐ and strike‐slip‐faulting types. Earthquakes occur regularly in time during the deployment, but are not evenly distributed in space across the study area, that is, they are spatially clustered. Analysis of the nearest‐neighbor distances in the space–time–magnitude domain shows the seismicity is dominated by single‐event clusters (i.e., independent events). This high proportion of single‐event clusters compared with multievent clusters has been previously noted for induced events at geothermal sites. When clustering occurs, the number of events in a cluster is typically small. We observe only four clusters with 10 or more events. For these larger clusters, we find equivalent numbers of foreshocks and aftershocks; however, the foreshock sequences are significantly longer in duration lasting days to tens of days, while aftershock sequences are observed only on the order of one day. The minimal clustering observed for events in the LASSO array suggests that the majority of events are being directly driven by stress changes due to local saltwater disposal.

Oklahoma↗

High-resolution seismic-reflection investigation of the northern Gulf of Mexico gas-hydrate-stability zone

We recorded high-resolution seismic-reflection data in the northern Gulf of Mexico to study gas and gas-hydrate distribution and their relation to seafloor slides. Gas hydrate is widely reported near the seafloor, but is described at only one deep drill site. Our data show high-reflectivity zones (HRZs) near faults, diapirs, and gas vents and interbedded within sedimentary sections at shallow depth (<1 km). The HRZs lie below the gas-hydrate-stability zone (GHSZ) as well as within the zone (less common), and they coincide with zones of shallow water-flows. Bottom simulating reflections are rare in the Gulf, and not documented in our data. We infer HRZs result largely from free gas in sandy beds, with gas hydrate within the GHSZ. Our estimates for the base BHSZ correlate reasonably with the top of HRZs in some thick well-layered basin sections, but poorly where shallow sediments are thin and strongly deformed. The equivocal correlation results from large natural variability of parameters that are used to calculate the base of the GHSZ. The HRZs may, however, be potential indicators of nearby gas hydrate. The HRZs also lie at the base of at least two large seafloor slides (e.g. up to 250 km2) that may be actively moving along decollement faults that sole within the GHSZ or close to the estimated base of the GHSZ. We suspect that water/gas flow along these and other faults such as 'chimney' features provide gas to permit crystallization of gas hydrate in the GHSZ. Such flows weaken sediment that slide down salt-oversteepened slopes when triggered by earthquakes. Published by Elsevier Science Ltd.

Marine and Petroleum Geology↗

Broadband records of earthquakes in deep gold mines and a comparison with results from SAFOD, California

For one week during September 2007, we deployed a temporary network of field recorders and accelerometers at four sites within two deep, seismically active mines. The ground-motion data, recorded at 200 samples/sec, are well suited to determining source and ground-motion parameters for the mining-induced earthquakes within and adjacent to our network. Four earthquakes with magnitudes close to 2 were recorded with high signal/noise at all four sites. Analysis of seismic moments and peak velocities, in conjunction with the results of laboratory stick-slip friction experiments, were used to estimate source processes that are key to understanding source physics and to assessing underground seismic hazard. The maximum displacements on the rupture surfaces can be estimated from the parameter , where is the peak ground velocity at a given recording site, and R is the hypocentral distance. For each earthquake, the maximum slip and seismic moment can be combined with results from laboratory friction experiments to estimate the maximum slip rate within the rupture zone. Analysis of the four M 2 earthquakes recorded during our deployment and one of special interest recorded by the in-mine seismic network in 2004 revealed maximum slips ranging from 4 to 27 mm and maximum slip rates from 1.1 to 6.3 m/sec. Applying the same analyses to an M 2.1 earthquake within a cluster of repeating earthquakes near the San Andreas Fault Observatory at Depth site, California, yielded similar results for maximum slip and slip rate, 14 mm and 4.0 m/sec.

Bulletin of the Seismological Society of America↗

Tidal and seasonal variations in calving flux observed with passive seismology

The seismic signatures of calving events, i.e., calving icequakes, offer an opportunity to examine calving variability with greater precision than is available with other methods. Here using observations from Yahtse Glacier, Alaska, we describe methods to detect, locate, and characterize calving icequakes. We combine these icequake records with a coincident, manually generated record of observed calving events to develop and validate a statistical model through which we can infer iceberg sizes from the properties of calving icequakes. We find that the icequake duration is the single most significant predictor of an iceberg's size. We then apply this model to 18 months of seismic recordings and find elevated iceberg calving flux during the summer and fall and a pronounced lull in calving during midwinter. Calving flux is sensitive to semidiurnal tidal stage. Large calving events are tens of percent more likely during falling and low tides than during rising and high tides, consistent with a view that deeper water has a stabilizing influence on glacier termini. Multiple factors affect the occurrence of mechanical fractures that ultimately lead to iceberg calving. At Yahtse Glacier, seismology allows us to demonstrate that variations in the rate of submarine melt are a dominant control on iceberg calving rates at seasonal timescales. On hourly to daily timescales, tidal modulation of the normal stress against the glacier terminus reveals the nonlinear glacier response to changes in the near-terminus stress field.

Alaska↗

Hydrologic effects of the earthquake of March 27, 1964, outside Alaska, with sections on Hydroseismograms from the Nunn-Bush Shoe Co. well, Wisconsin, and Alaska earthquake effects on ground water in Iowa: Chapter C in The Alaska earthquakes, March 27, 1964: effects on hydrologic regimen

The Alaska earthquake of March 27, 1964, had widespread hydrologic effects throughout practically all of the United States. More than 1,450 water-level recorders, scattered throughout all the 50 States except Connecticut, Delaware, and Rhode Island, registered the earthquake. Half of the water-level records were obtained from ground-water observation wells and half at surface-water gaging stations. The earthquake is also known to have registered on water-level recorders on wells in Canada, England, Denmark, Belgium, Egypt, Israel, Libya, Philippine Islands, South-West Africa, South Africa, and Northern Territory of Australia. The Alaska earthquake is the first for which widespread surface-water effects are known. The effects were recorded at stations on flowing streams, rivers, reservoirs, lakes, and ponds. The 755 surface-water stations recording effects are spread through 38 States, but are most numerous in the south-central and southeastern States, especially in Florida and Louisiana. Most of the fluctuations recorded can be referred to more precisely as seismic seiches; however, a few stations recorded the quake as a minor change in stage. The largest recorded seiche outside Alaska was 1.83 feet on a reservoir in Michigan. The next largest was 1.45 feet on Lake Ouachita in Arkansas. The largest fluctuation in a well was 23 feet registered by a pressure recorder near Belle Fourche, S. Dak. Fluctuations of more than 10 feet were reported from wells in Alabama, Florida, Georgia, Illinois, Missouri, and Pennsylvania. A 3.40-foot fluctuation was recorded in a well in Puerto Rico. The Alaska earthquake was registered by about seven times as many water-level recorders as recorded the Hebgen Lake, Mont., earthquake of August 19, 1959.

Professional Paper↗

Evolution of the 2015 Cotopaxi eruption revealed by combined geochemical & seismic observations

Through integration of multiple data streams to monitor volcanic unrest scientists are able to make more robust eruption forecast and to obtain a more holistic interpretation of volcanic systems. We examined gas emission and gas geochemistry, seismic and petrologic data recorded during the 2015 unrest of Cotopaxi (Ecuador) in order to decipher the origin and temporal evolution of this eruption. Identification of families of similar seismic events and the use of seismic amplitude ratios reveals temporal changes in volcanic processes. SO2 (300 to 24000 t/d), BrO/SO2 (5-10 x10-5), SO2/HCl (5.8 ± 4.8 and 6.6 ± 3.0) and CO2/SO2 (0.6 to 2.1) measured throughout the eruption indicate a shallow magmatic source. Bulk ash and glass chemistry indicate a homogenous andesitic (SiO2 wt%=56.94 ± 0.25) magma having undergone extensive S-exsolution and degassing during ascent. These data lead us to interpret this eruption as a magma intrusion and ascend to shallow levels. The intrusion progressively interacted with the hydrothermal system, boiled off water, and produced hydromagmatic explosions. A small volume of this intrusion continued to fragment and produced episodic ash emissions until it was sufficiently degassed and rheologically stiff. Based on the 470 kt of measured SO2 we estimate that ~ 65.3 x106 m3 of magma were required to supply the emitted gases. This volume exceeds the volume of erupted juvenile material by a factor of 50. This result emphasizes the importance of careful monitoring of Cotopaxi to identify the intrusion of a new batch of magma, which could rejuvenate the non-erupted material.

Geochemistry, Geophysics, Geosystems↗

Teleseismic and near-field analysis of the Nahanni earthquakes in the Northwest Territories, Canada

The analysis of the Nahanni earthquakes of October 5, 1985 (MS 6.6), and December 23, 1985 (MS 6.9), will have important implications for the assessment of seismic hazards in intraplate environments. To maximize the information available to seismic engineers, broadband data recorded teleseismically are analyzed jointly with strong-motion data recorded in the near field. The time-domain analysis of teleseismic data yields the source mechanisms, depths, and complexities of rupture of each earthquake. Both earthquakes occurred as shallow thrusts with centroid depths (6 to 7 km) and shallowly dipping fault planes that correspond well with the aftershock distributions obtained from a local survey run by the Canadian Geological Survey. The shallow nodal plane for the October 5 earthquake dips 30° to the WSW, while the shallow nodal planes of the subevents for the December 23 earthquake dip an average of 23° to the WSW. The October 5 earthquake has an impulsive initial rupture, followed by a weak subevent of longer duration but smaller moment release. The December 23 earthquake exhibits more complexity, being comprised of three subevents of similar size. The subevent delays derived from the teleseismic analysis are used to help interpret arrivals in records of ground velocity recorded in the near field of the December 23 earthquake. The rupture geometries inferred from the joint near- and far-field analysis suggest that the rupture processes were unusually complicated and that the 2g peak that occurs late in one of the near-field records could be a localized phenomenon. Spectral analyses of the teleseismic P waves yield the following source parameters for the October 5 and December 23 earthquakes, respectively: the seismic moments are 1.2 and 1.8 × 1026 dyne-cm, the radiated energies are 1.8 and 2.8 × 1021 dyne-cm, and the dynamic stress drops are 65 and 50 bar. The acceleration source spectra of both earthquakes exhibit an intermediate slope (| üα(ω) | ∝ω) from 0.03 to 0.3 Hz, suggesting that the earthquakes represent the failure of asperities. Extrapolating the teleseismic P-wave spectra to estimate the near-field S-wave spectra yields good fits to the acceleration spectra from two strong motion records, but underestimates the spectra from a third strong motion record with the strongest, but possibly localized, accelerations.

Bulletin of the Seismological Society of America↗

Probabilistic seismic demand analysis using advanced ground motion intensity measures

One of the objectives in performance-based earthquake engineering is to quantify the seismic reliability of a structure at a site. For that purpose, probabilistic seismic demand analysis (PSDA) is used as a tool to estimate the mean annual frequency of exceeding a specified value of a structural demand parameter (e.g. interstorey drift). This paper compares and contrasts the use, in PSDA, of certain advanced scalar versus vector and conventional scalar ground motion intensity measures (IMs). One of the benefits of using a well-chosen IM is that more accurate evaluations of seismic performance are achieved without the need to perform detailed ground motion record selection for the nonlinear dynamic structural analyses involved in PSDA (e.g. record selection with respect to seismic parameters such as earthquake magnitude, source-to-site distance, and ground motion epsilon). For structural demands that are dominated by a first mode of vibration, using inelastic spectral displacement (Sdi) can be advantageous relative to the conventionally used elastic spectral acceleration (Sa) and the vector IM consisting of Sa and epsilon (??). This paper demonstrates that this is true for ordinary and for near-source pulse-like earthquake records. The latter ground motions cannot be adequately characterized by either Sa alone or the vector of Sa and ??. For structural demands with significant higher-mode contributions (under either of the two types of ground motions), even Sdi (alone) is not sufficient, so an advanced scalar IM that additionally incorporates higher modes is used.

Earthquake Engineering and Structural Dynamics↗

Analysis of high-resolution single channel seismic data for use in sediment resource evaluation, eastern Texas and western Louisiana Continental Shelf, Gulf of Mexico

Shallow subsurface geologic data recorded as high-resolution seismic profiles are used to interpret the geology of coastal and marine systems. These data were originally recorded on paper rolls that are stored in geophysical archives. Data collection has since converted to entirely digital formats, yet the analog data are still useful for geologic interpretation. This report describes the process of recovering analog copies of seismic profiles from physical archives, electronically scanning, and converting them to industry-standard digital format. The recovered data are also reviewed and assessed for potential sediment resources. The data recovered in this study were collected from the Gulf of Mexico continental shelf offshore of East Texas and West Louisiana. The project is a collaborative study between the U.S. Geological Survey and the Bureau of Ocean Energy Management.

Louisiana, Texas↗

Shallow seismic structure of Kunlun fault zone in northern Tibetan Plateau, China: Implications for the 2001 M s8.1 Kunlun earthquake

The shallow seismic velocity structure of the Kunlun fault zone (KLFZ) was jointly deduced from seismic refraction profiling and the records of trapped waves that were excited by five explosions. The data were collected after the 2001 Kunlun M s8.1 earthquake in the northern Tibetan Plateau. Seismic phases for the in-line record sections (26 records up to a distance of 15 km) along the fault zone were analysed, and 1-D P - and S -wave velocity models of shallow crust within the fault zone were determined by using the seismic refraction method. Sixteen seismic stations were deployed along the off-line profile perpendicular to the fault zone. Fault-zone trapped waves appear clearly on the record sections, which were simulated with a 3-D finite difference algorithm. Quantitative analysis of the correlation coefficients of the synthetic and observed trapped waveforms indicates that the Kunlun fault-zone width is 300 m, and S -wave quality factor Q within the fault zone is 15. Significantly, S -wave velocities within the fault zone are reduced by 30–45 per cent from surrounding rocks to a depth of at least 1–2 km, while P -wave velocities are reduced by 7–20 per cent. A fault-zone with such P - and S -low velocities is an indication of high fluid pressure because Vs is affected more than Vp . The low-velocity and low- Q zone in the KLFZ model is the effect of multiple ruptures along the fault trace of the 2001 M s8.1 Kunlun earthquake.

Tibetan Plateau↗

Volcano dome dynamics at Mount St. Helens: Deformation and intermittent subsidence monitored by seismicity and camera imagery pixel offsets

The surface deformation field measured at volcanic domes provides insights into the effects of magmatic processes, gravity- and gas-driven processes, and the development and distribution of internal dome structures. Here we study short-term dome deformation associated with earthquakes at Mount St. Helens, recorded by a permanent optical camera and seismic monitoring network. We use Digital Image Correlation (DIC) to compute the displacement field between successive images and compare the results to the occurrence and characteristics of seismic events during a 6 week period of dome growth in 2006. The results reveal that dome growth at Mount St. Helens was repeatedly interrupted by short-term meter-scale downward displacements at the dome surface, which were associated in time with low-frequency, large-magnitude seismic events followed by a tremor-like signal. The tremor was only recorded by the seismic stations closest to the dome. We find a correlation between the magnitudes of the camera-derived displacements and the spectral amplitudes of the associated tremor. We use the DIC results from two cameras and a high-resolution topographic model to derive full 3-D displacement maps, which reveals internal dome structures and the effect of the seismic activity on daily surface velocities. We postulate that the tremor is recording the gravity-driven response of the upper dome due to mechanical collapse or depressurization and fault-controlled slumping. Our results highlight the different scales and structural expressions during growth and disintegration of lava domes and the relationships between seismic and deformation signals.

Washington↗

Seismic hazard and risk assessment in the intraplate environment: The New Madrid seismic zone of the central United States

Although the causes of large intraplate earthquakes are still not fully understood, they pose certain hazard and risk to societies. Estimating hazard and risk in these regions is difficult because of lack of earthquake records. The New Madrid seismic zone is one such region where large and rare intraplate earthquakes (M = 7.0 or greater) pose significant hazard and risk. Many different definitions of hazard and risk have been used, and the resulting estimates differ dramatically. In this paper, seismic hazard is defined as the natural phenomenon generated by earthquakes, such as ground motion, and is quantified by two parameters: a level of hazard and its occurrence frequency or mean recurrence interval; seismic risk is defined as the probability of occurrence of a specific level of seismic hazard over a certain time and is quantified by three parameters: probability, a level of hazard, and exposure time. Probabilistic seismic hazard analysis (PSHA), a commonly used method for estimating seismic hazard and risk, derives a relationship between a ground motion parameter and its return period (hazard curve). The return period is not an independent temporal parameter but a mathematical extrapolation of the recurrence interval of earthquakes and the uncertainty of ground motion. Therefore, it is difficult to understand and use PSHA. A new method is proposed and applied here for estimating seismic hazard in the New Madrid seismic zone. This method provides hazard estimates that are consistent with the state of our knowledge and can be easily applied to other intraplate regions. ?? 2007 The Geological Society of America.

Special Paper of the Geological Society of America↗

Upper crustal structure in Puget Lowland, Washington: Results from the 1998 Seismic Hazards Investigation in Puget Sound

A new three-dimensional (3-D) model shows seismic velocities beneath the Puget Lowland to a depth of 11 km. The model is based on a tomographic inversion of nearly one million first-arrival travel times recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS), allowing higher-resolution mapping of subsurface structures than previously possible. The model allows us to refine the subsurface geometry of previously proposed faults (e.g., Seattle, Hood Canal, southern Whidbey Island, and Devils Mountain fault zones) as well as to identify structures (Tacoma, Lofall, and Sequim fault zones) that warrant additional study. The largest and most important of these newly identified structures lies along the northern boundary of the Tacoma basin; we informally refer to this structure here as the Tacoma fault zone. Although tomography cannot provide information on the recency of motion on any structure, Holocene earthquake activity on the Tacoma fault zone is suggested by seismicity along it and paleoseismic evidence for abrupt uplift of tidal marsh deposits to its north. The tomography reveals four large, west to northwest trending low-velocity basins (Tacoma, Seattle, Everett, and Port Townsend) separated by regions of higher velocity ridges that are coincident with fault-bounded uplifts of Eocene Crescent Formation basalt and pre-Tertiary basement. The shapes of the basins and uplifts are similar to those observed in gravity data; gravity anomalies calculated from the 3-D tomography model are in close agreement with the observed anomalies. In velocity cross sections the Tacoma and Seattle basins are asymmetric: the basin floor dips gently toward a steep boundary with the adjacent high-velocity uplift, locally with a velocity "overhang" that suggests a basin vergent thrust fault boundary. Crustal fault zones grow from minor folds into much larger structures along strike. Inferred structural relief across the Tacoma fault zone increases by several kilometers westward along the fault zone to Lynch Cove, where we interpret it as a zone of south vergent faulting overthrusting Tacoma basin. In contrast, structural relief along the Seattle fault zone decreases west of Seattle, which we interpret as evidence that the N-S directed compression is being accommodated by slip transfer between the Seattle and Tacoma fault zones. Together, the Tacoma and Seattle fault zones raise the Seattle uplift, one of a series of east-west trending, pop-up structures underlying Puget Lowland from the Black Hills to the San Juan Islands.

Washington↗