USGS Science⌕ Search

SEARCH · USGS Science

Results for “earthquake hazards”

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,783 records · Page 99Linked to original sources

The USGS 2023 Conterminous U.S. time‐independent earthquake rupture forecast

We present the 2023 U.S. Geological Survey time‐independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time‐averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19‐member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time‐dependent components, and assess implications with a wider range of hazard and risk metrics.

Bulletin of the Seismological Society of America↗

Map and Database of Probable and Possible Quaternary Faults in Afghanistan

The U.S. Geological Survey (USGS) with support from the U.S. Agency for International Development (USAID) mission in Afghanistan, has prepared a digital map showing the distribution of probable and suspected Quaternary faults in Afghanistan. This map is a key component of a broader effort to assess and map the country's seismic hazards. Our analyses of remote-sensing imagery reveal a complex array of tectonic features that we interpret to be probable and possible active faults within the country and in the surrounding border region. In our compilation, we have mapped previously recognized active faults in greater detail, and have categorized individual features based on their geomorphic expression. We assigned mapped features to eight newly defined domains, each of which contains features that appear to have similar styles of deformation. The styles of deformation associated with each domain provide insight into the kinematics of the modern tectonism, and define a tectonic framework that helps constrain deformational models of the Alpine-Himalayan orogenic belt. The modern fault movements, deformation, and earthquakes in Afghanistan are driven by the collision between the northward-moving Indian subcontinent and Eurasia. The patterns of probable and possible Quaternary faults generally show that much of the modern tectonic activity is related to transfer of plate-boundary deformation across the country. The left-lateral, strike-slip Chaman fault in southeastern Afghanistan probably has the highest slip rate of any fault in the country; to the north, this slip is distributed onto several fault systems. At the southern margin of the Kabul block, the style of faulting changes from mainly strike-slip motion associated with the boundary between the Indian and Eurasian plates, to transpressional and transtensional faulting. North and northeast of the Kabul block, we recognized a complex pattern of potentially active strike-slip, thrust, and normal faults that form a conjugate shear system in a transpressional region of the Trans-Himalayan orogenic belt. The general patterns and orientations of faults and the styles of deformation that we interpret from the imagery are consistent with the styles of faulting determined from focal mechanisms of historical earthquakes. Northwest-trending strike-slip fault zones are cut and displaced by younger, southeast-verging thrust faults; these relations define the interaction between northwest-southeast-oriented contraction and northwest-directed extrusion in the western Himalaya, Pamir, and Hindu Kush regions. Transpression extends into north-central Afghanistan where north-verging contraction along the east-west-trending Alburz-Marmul fault system interacts with northwest-trending strike-slip faults. Pressure ridges related to thrust faulting and extensional basins bounded by normal faults are located at major stepovers in these northwest-trending strike-slip systems. In contrast, young faulting in central and western Afghanistan indicates that the deformation is dominated by extension where strike-slip fault zones transition into regions of normal faults. In addition to these initial observations, our digital map and database provide a foundation that can be expanded, complemented, and modified as future investigations provide more detailed information about the location, characteristics, and history of movement on Quaternary faults in Afghanistan.

Open-File Report↗

Probabilistic seismic-hazard analysis for the western Kingdom of Saudi Arabia

We present a probabilistic seismic-hazard analysis (PSHA) for the west-central part of the Arabian Peninsula. Our study area includes the northern Harrat Rahat volcanic field and the nearby city of Al Madīnah, Kingdom of Saudi Arabia. This young, active volcanic field experienced one historical eruption in 1256 C.E. (654 in the year of the Hijra) that vented 20 to 22 kilometers (km) southeast of the center of Al Madīnah, which has a present population of about 1.4 million. The field also erupted numerous times in the late Pleistocene and possibly in the early Holocene. We used recently developed regional ground-motion prediction equations for Saudi Arabia to calculate the severity of ground shaking as a function of distance from the earthquake source. This information was combined with two key volcanic parameters for this area: (1) the return period of volcanic eruptions and (2) the spatial probability of the next vent opening. The calculated ground-motion levels of peak ground acceleration and peak ground velocity with a 2-percent probability of exceedance in 50 years are expected to be about 0.14 the acceleration due to gravity ( g ) and 10 centimeters per second (cm/s), respectively, at the most probable vent opening location, about 25 km southeast of Al Madīnah’s center, and about 0.07 g and 3 cm/s, respectively, in the city interior. These ground motions are higher than previous estimates that did not consider the nearby earthquakes associated with volcanic activity.

Professional Paper↗

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

Numerical analysis of the mobility of the Palos Verdes debris avalanche, California, and its implication for the generation of tsunamis

Analysis of morphology, failure and post-failure stages of the Palos Verdes debris avalanche reveals that it may have triggered a significant tsunami wave. Our analysis of the failure itself indicates that the slope is stable under aseismic conditions but that a major earthquake (with a magnitude around 7) could have triggered the slide. A post-failure analysis, considering the debris avalanche as a bi-linear flow, shows that peak velocities of up to 45 m/s could have been reached and that the initial movement involved a mass of rock less than 10 km wide, 1 km long and about 50-80 m thick. Initial wave height estimates vary from 10 to 50 m. Tsunami waves propagating to the local shoreline would be significantly smaller. Such a range demonstrates our lack of proper knowledge of the transition from failure to post-failure behavior related to mass movements. Further investigations and analyses of terrestrial and submarine evidence are required for a proper hazard assessment related to tsunami generation in the Los Angeles area. ?? 2003 Elsevier B.V. All rights reserved.

Conference Paper↗

Implications of silicic vent patterns for the presence of large crustal magma chambers

On the basis of the distribution of silicic vents, many volcanic fields can be grouped with (1) igneous systems that may be small and whose vent locations are controlled by regional tectonics, (2) those that include sizable crustal magma bodies which erupt at sites determined by their anomalous local stress fields, or (3) relatively small volume systems that are transitional between categories 1 and 2. Linear vent patterns that are aligned normal to the regional least principal stress (σ 3 ) commonly are associated with absence of evidence for large, shallow magma bodies. The Coso volcanic field and the Inyo‐Mono domes in California and probably the South Sister area in Oregon are examples of such systems. The 1960 dacitic fissure eruption at Cordón Caulle in southern Chile evidently is linked to tectonic stress relaxation associated with the great earthquake that occurred 48 hours before the eruption began. Large, shallow magma chambers are thought to perturb the local stress field so that areal patterns of silicic vents are diffuse, radial, or arcuate. Such systems may erupt great volumes of pyroclastic material catastrophically and produce large calderas. Well‐preserved examples of late precaldera leaks of silicic magma occur at Long Valley, California, and Mount Mazama (Crater Lake), Oregon. Other possible examples are noted. Some of these, those which formed calderas smaller than Crater Lake, apparently were preceded by silicic eruptions from aligned vents. Bearing in mind that there exists a transitional group between tectonically controlled small systems and very large magma chambers, vent distributions can be useful in evaluating potential volcanic hazards for silicic volcanic fields.

California↗

New strategy needed in earthquake, volcano monitoring

Recent advances in space geodesy provide unprecedented opportunities for measuring and understanding processes related to earthquake occurrence and volcanic eruptions in the United States and elsewhere. The Global Positioning System (GPS) uses Earth-orbiting satellites to obtain relative movements of ground points accurate to a few millimeters, either through periodically repeated surveys or by continuous measurements at permanent sites [ Segall and Davis , 1997]. Satellite interferometric synthetic aperture radar (InSAR) uses repeat-pass radar backscatter images of Earth's surface to obtain complete spatial mappings of surface deformation over 100 km × 100 km scenes to centimeter precision [ Massonnet and Feigl , 1998]. During the past 5 years GPS and InSAR increasingly have been applied to local studies of active fault zones and volcanic systems, imaging the sources of deformation buried in Earth's crust and quantifying the hazards they pose to society. The potential now exists to deploy these tools to map and monitor all of the actively deforming western United States, and a new scientific initiative with these goals is under active discussion [ Silver et al , 1998].

Eos Science News↗

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗

Discovery of an active forearc fault in an urban region: Holocene rupture on the XEOLXELEK-Elk Lake fault, Victoria, British Columbia, Canada

Subduction forearcs are subject to seismic hazard from upper plate faults that are often invisible to instrumental monitoring networks. Identifying active faults in forearcs therefore requires integration of geomorphic, geologic, and paleoseismic data. We demonstrate the utility of a combined approach in a densely populated region of Vancouver Island, Canada, by combining remote sensing, historical imagery, field investigations, and shallow geophysical surveys to identify a previously unrecognized active fault, the XEOLXELEK-Elk Lake fault, in the northern Cascadia forearc, ∼10 km north of the city of Victoria. Lidar-derived digital terrain models and historical air photos show a ∼2.5-m-high scarp along the surface of a Quaternary drumlinoid ridge. Paleoseismic trenching and electrical resistivity tomography surveys across the scarp reveal a single reverse-slip earthquake produced a fault-propagation fold above a blind southwest-dipping fault. Five geologically plausible chronological models of radiocarbon dated charcoal constrain the likely earthquake age to between 4.7 and 2.3 ka. Fault-propagation fold modeling indicates ∼3.2 m of reverse slip on a blind, 50° southwest-dipping fault can reproduce the observed deformation. Fault scaling relations suggest a M 6.1–7.6 earthquake with a 13 to 73-km-long surface rupture and 2.3–3.2 m of dip slip may be responsible for the deformation observed in the paleoseismic trench. An earthquake near this magnitude in Greater Victoria could result in major damage, and our results highlight the importance of augmenting instrumental monitoring networks with remote sensing and field studies to identify and characterize active faults in similarily challenging environments.

Victoria↗

Dynamic rupture modeling of the M7.2 2010 El Mayor-Cucapah earthquake: Comparison with a geodetic model

The 2010 M w 7.2 El Mayor-Cucapah earthquake is the largest event recorded in the broader Southern California-Baja California region in the last 18 years. Here we try to analyze primary features of this type of event by using dynamic rupture simulations based on a multifault interface and later compare our results with space geodetic models. Our results show that starting from homogeneous prestress conditions, slip heterogeneity can be achieved as a result of variable dip angle along strike and the modulation imposed by step over segments. We also considered effects from a topographic free surface and find that although this does not produce significant first-order effects for this earthquake, even a low topographic dome such as the Cucapah range can affect the rupture front pattern and fault slip rate. Finally, we inverted available interferometric synthetic aperture radar data, using the same geometry as the dynamic rupture model, and retrieved the space geodetic slip distribution that serves to constrain the dynamic rupture models. The one to one comparison of the final fault slip pattern generated with dynamic rupture models and the space geodetic inversion show good agreement. Our results lead us to the following conclusion: in a possible multifault rupture scenario, and if we have first-order geometry constraints, dynamic rupture models can be very efficient in predicting large-scale slip heterogeneities that are important for the correct assessment of seismic hazard and the magnitude of future events. Our work contributes to understanding the complex nature of multifault systems.

California↗

High-frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington↗

Elders recall an earlier tsunami on Indian Ocean shores

Ten years on, the Indian Ocean tsunami of 26 December 2004 still looms large in efforts to reduce coastal risk. The disaster has spurred worldwide advances in tsunami detection and warning, tsunami-risk assessment, and tsunami awareness [Satake, 2014]. Nearly a lifetime has passed since the northwestern Indian Ocean last produced a devastating tsunami. Documentation of this tsunami, in November 1945, was hindered by international instability in the wake of the Second World War and, in British India, by the approach of independence and partition. The parent earthquake, of magnitude 8.1, was widely recorded, and the tsunami registered on tide gauges, but intelligence reports and newspaper articles say little about inundation limits while permitting a broad range of catalogued death tolls. What has been established about the 1945 tsunami falls short of what's needed today for ground-truthing inundation models, estimating risk to enlarged populations, and anchoring awareness campaigns in local facts. Recent efforts to reduce coastal risk around the Arabian Sea include a project in which eyewitnesses to the 1945 tsunami were found and interviewed (Fig. 1), and related archives were gathered. Results are being made available through UNESCO's Indian Ocean Tsunami Information Center in hopes of increasing scientific understanding and public awareness of the region's tsunami hazards.

Eos, Earth and Space Science News↗

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↗

Meet the people where they are: Assessing user needs for aftershock forecast products in El Salvador, Mexico and the United States

Aftershock forecasts can help communities reduce their seismic risk by conveying how many aftershocks can be expected following a large earthquake, and how the expected number of aftershocks and their corresponding ground shaking evolves over time and space. Prior work finds that graphical forecast products may communicate such information better than only text or numbers. To identify which visual products can serve multiple user groups, we held workshops with members of several professions, including emergency managers, engineers, critical infrastructure operators, public health specialists, science communicators, and more. We conducted these workshops in El Salvador, Mexico and the United States to understand which forecast products may be effective across different countries. In these workshops, users performed small-group activities to elicit the types of aftershock forecast information that would support decisions in their respective roles and how this information would optimally be displayed. Maps of shaking hazards were frequently requested across all professions and countries, even for dissimilar forecast uses. The design of these maps, and other product needs, were differentiated by profession, country and other factors. Other forecast products, including those showing forecasts about the magnitudes and time periods of future aftershocks, also served a variety of users, but for different types of decisions. We found a greater variation in user needs by profession than by country, and that user needs also vary with time, communication channel and other contextual factors. We discuss practical implications for user-centered visual communication of operational aftershock forecasts.

International Journal for Disaster Risk Reduction↗

Characterizing the Kathmandu Valley sediment response through strong motion recordings of the 2015 Gorkha earthquake sequence

We analyze strong motion records and high-rate GPS measurements of the M 7.8 Gorkha mainshock, M 7.3 Dolakha, and two moderate aftershock events recorded at four stations on the Kathmandu basin sediments, and one on rock-outcrop. Recordings on soil from all four events show systematic amplification relative to the rock site at multiple frequencies in the 0.1–2.5 Hz frequency range, and de-amplification of higher frequencies ( >2.5–10 Hz). The soil-to-rock amplification ratios for the M 7.8 and M 7.3 events have lower amplitude and frequency peaks relative to the ratios of the two moderate events, effects that could be suggestive of nonlinear site response. Further, comparisons to ground motion prediction equations show that 1) both soil and rock mainshock recordings were severely depleted of high frequencies, and 2) the depletion at high frequencies is not present in the aftershocks. These observations indicate that the high frequency deamplification is additionally related to characteristics of the source that are not captured by simplified ground motion prediction equations, and allude to seismic hazard analysis models being revised – possibly by treating isolated high frequency radiation sources separately from long period components to capture large magnitude near-source events such as the 2015 Gorkha mainshock.

Kathmandu Valley↗

Rupture passing probabilities at fault bends and steps, with application to rupture length probabilities for earthquake early warning

Earthquake early warning (EEW) systems can quickly identify the beginning of a significant earthquake rupture, but the first seconds of seismic data have not been found to predict the final rupture length. We present two approaches for estimating probabilities of rupture length given the rupture initiation from an EEW system. In the first approach, bends and steps on the fault are interpreted as physical mechanisms for rupture arrest. Arrest probability relations are developed from empirical observations and depend on bend angle and step size. Probability of arrest compounds serially with increasing rupture length as bends or steps are encountered. In the second approach, time‐independent rates among ruptures from the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3), are interpreted to apply to the time‐dependent condition in which rupture grows from a known starting point. Length probabilities from a Gutenberg–Richter magnitude–frequency relation provide a reference of comparison. We illustrate the new approach using the discretized fault model for California developed for UCERF3. For the case of rupture initiating on the southeast end of the San Andreas fault we find the geometric complexity of the Mill Creek section impedes most ruptures, and only ∼5% are predicted to reach to San Bernardino on the eastern edge of the greater Los Angeles region. Conditional probabilities of length can be precompiled in this manner for any initiation point on the fault system and thus are of potential value in seismic hazard and EEW applications.

California↗

Liquefaction and foundation failure of Chevron oil and gasoline tanks at Moss Landing, California

Liquefaction of gray, coarse‐ to medium‐grained, esturine sand during the October 17, 1989 Loma Prieta (Ms=7.1) earthquake led to foundation failures of oil and tanks at the Chevron Marine Station at Moss Landing located about 25 km south‐southeast of the epicenter. Differential settlement of the foundations resulted in 1° to 6° of the tanks towards one another and an associated 2° to 7° inclination of tank piping. In addition, sand was vented around the bases of the tanks that contained the most oil and adjacent cracks. Liquefaction and related ground failure was not restricted to the vicinity of the Chevron tanks but was a widespread phenomenon at Moss Landing. Analysis of geological and geotechnical data at sites such as the Chevron Marine Station will serve to check and to advance current methods of assessing liquefaction potential at a site and of mapping liquefaction hazard on a regional scale.

California↗

Bibliography of literature pertaining to Long Valley Caldera and associated volcanic fields

On May 25-27, 1980, Long Valley caldera was rocked by four M=6 earthquakes that heralded the onset of a wave of seismic activity within the caldera which has continued through the present. Unrest has taken the form of seismic swarms, uplift of the resurgent dome, and areas of vegetation killed by increased CO2 emissions, all interpreted as resulting from magma injection into different levels beneath the caldera, as well as beneath Mammoth Mountain along the southwest rim of the caldera. Continuing economic development in the Mammoth Lakes area has swelled the local population, increasing the risk to people and property if an eruption were to occur. The U.S. Geological Survey (USGS) has been monitoring geophysical activity in the Long Valley area since the mid-1970s and continues to track the unrest in real time with a sophisticated network of geophysical sensors. Hazards information obtained by this monitoring is provided to local, State, and Federal officials and to the public through the Long Valley Observatory. The Long Valley area also was scientifically important before the onset of current unrest. Lying at the eastern foot of the Sierra Nevada, the deposits from this active volcanic system have provided fertile ground for research into Neogene tectonics, Quaternary geology and geomorphology, regional stratigraphy, and volcanology. In the early 1970s, intensive studies of the area began through the USGS Geothermal Investigations Program, owing to the presence of a large young silicic volcanic system. The paroxysmal eruption of Long Valley caldera about 760,000 years ago produced the Bishop Tuff and associated Bishop ash. The Bishop Tuff is a well-preserved ignimbrite deposit that has continued to provide new and developing insights into the dynamics of ignimbrite-forming eruptions. Another extremely important aspect of the Bishop Tuff is that it is the oldest known normally magnetized unit of the Brunhes Chron. Thus, the age of the Bishop Tuff is used to define the beginning of the Brunhes Chron and helps constrain the Brunhes-Matuyama boundary. The Bishop ash, which was dispersed as far east as Nebraska, Kansas, and Texas, provides an important tephrostratigraphic marker throughout the Western United States. The obsidian domes of both the Mono and Inyo Craters, which were produced by rhyolitic eruptions in the past 40,000 years, have been well studied, including extensive scientific drilling through the domes. Exploratory drilling to 3-km depth on the resurgent dome and subsequent instrumentation of the Long Valley Exploratory Well (LVEW) have led to a number of important new insights. Scientific drilling also has been done within the Casa Diablo geothermal field, which, aside from drilling, has been commercially developed and is currently feeding 40 MW of power into the Southern California Edison grid. Studies in all the above-mentioned volcanic fields have contributed to the extensive scientific literature published on the Long Valley region. Although most of this scientific literature has been published since 1970, a significant amount of historical literature extends backward to the late 1800s. The purpose of this bibliography is to compile references pertaining to the Long Valley region from all time periods and all Earth science fields into a single listing, thus providing an easily accessible guide to the published literature for current and future researchers.

Open-File Report↗