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,045 records · Page 58Linked to original sources
USGS1994

Seismic images of the Brooks Range fold and thrust belt, Arctic Alaska, from an integrated seismic reflection/refraction experiment

We describe results of an integrated seismic reflection/refraction experiment across the Brooks Range and flanking geologic provinces in Arctic Alaska. The seismic acquisition was unusual in that reflection and refraction data were collected simultaneously with a 700 channel seismograph system deployed numerous times along a 315 km profile. Shot records show continuous Moho reflections from 0–180 km offset, as well as numerous upper- and mid-crustal wide-angle events. Single and low-fold near-vertical incidence common midpoint (CMP) reflection images show complex upper- and middle-crustal structure across the range from the unmetamorphosed Endicott Mountains allochthon (EMA) in the north, to the metamorphic belts in the south. Lower-crustal and Moho reflections are visible across the entire reflection profile. Travel-time inversion of PmP arrivals shows that the Moho, at 33 km depth beneath the North Slope foothills, deepens abruptly beneath the EMA to a maximum of 46 km, and then shallows southward to 35 km at the southern edge of the range. Two zones of upper- and middle-crustal reflections underlie the northern Brooks Range above ~ 12–15 km depth. The upper zone, interpreted as the base of the EMA, lies at a maximum depth of 6 km and extends over 50 km from the range front to the north central Brooks Range where the base of the EMA outcrops above the metasedimentary rocks exposed in the Doonerak window. We interpret the base of the lower zone, at ~ 12 km depth, to be from carbonate rocks above the master detachment upon which the Brooks Range formed. The seismic data suggest that the master detachment is connected to the faults in the EMA by several ramps. In the highly metamorphosed terranes south of the Doonerak window, the CMP section shows numerous south-dipping events which we interpret as a crustal scale duplex involving the Doonerak window rocks. The basal detachment reflections can be traced approximately 100 km, and dip southward from about 10–12 km near the range front, to 14–18 km beneath the Doonerak window, to 26–28 km beneath the metamorphic belts in the central Brooks Range. The section documents middle- and lower-crustal involvement in the formation of the Brooks Range.

Alaska↗
USGS2020

California Historical Intensity Mapping Project (CHIMP): A consistently reinterpreted dataset of seismic intensities for the past 162 years and implications for seismic hazard maps

Historical seismic intensity data are useful for myriad reasons, including assessment of the performance of Probabilistic Seismic Hazard Assessment (PSHA) models and corresponding hazard maps by comparing their predictions to a dataset of historically observed intensities in the region. To assess PSHA models for California, a long and consistently interpreted intensity record is necessary. For this purpose, the California Historical Intensity Mapping Project (CHIMP) has compiled a dataset that combines and reinterprets intensity information that has been stored in disparate and sometimes hard-to-access locations. The CHIMP dataset also includes new observations of intensity from archival research and oral history collection. Version 1 of the dataset includes 46,502 intensity observations for 62 earthquakes with estimated magnitudes ranging from 4.7 to 7.9. The 162 years of shaking data show observed shaking lower than expected from seismic hazard models. This discrepancy is reduced, but persists, if historical intensity data for the largest earthquakes are smoothed to reduce the effects of spatial under-sampling. Possible reasons for this discrepancy include other limitations of the CHIMP dataset, the hazard models, and the possibility that California seismicity throughout the historical period has been lower than the long-term average. Some of these issues may also explain similar discrepancies observed for Italy and Japan.

California↗
USGS2013

Very long period conduit oscillations induced by rockfalls at Kilauea Volcano, Hawaii

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with an 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 August 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–1000 s. The VLP signals associated with the rockfalls originate in a source region ∼1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200–4500 m 3 ), length of magma column (120–210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3–6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400–15,000 kg s −1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on 29 August 2012.

Hawaii↗
USGS2007

A catastrophic meltwater flood event and the formation of the Hudson Shelf Valley

The Hudson Shelf Valley (HSV) is the largest physiographic feature on the U.S. mid-Atlantic continental shelf. The 150-km long valley is the submerged extension of the ancestral Hudson River Valley that connects to the Hudson Canyon. Unlike other incised valleys on the mid-Atlantic shelf, it has not been infilled with sediment during the Holocene. Analyses of multibeam bathymetry, acoustic backscatter intensity, and high-resolution seismic reflection profiles reveal morphologic and stratigraphic evidence for a catastrophic meltwater flood event that formed the modern HSV. The valley and its distal deposits record a discrete flood event that carved 15-m high banks, formed a 120-km 2 field of 3- to 6-m high bedforms, and deposited a subaqueous delta on the outer shelf. The HSV is inferred to have been carved initially by precipitation and meltwater runoff during the advance of the Laurentide Ice Sheet, and later by the drainage of early proglacial lakes through stable spillways. A flood resulting from the failure of the terminal moraine dam at the Narrows between Staten Island and Long Island, New York, allowed glacial lakes in the Hudson and Ontario basins to drain across the continental shelf. Water level changes in the Hudson River basin associated with the catastrophic drainage of glacial lakes Iroquois, Vermont, and Albany around 11,450 14 C year BP (∼ 13,350 cal BP) may have precipitated dam failure at the Narrows. This 3200 km 3 discharge of freshwater entered the North Atlantic proximal to the Gulf Stream and may have affected thermohaline circulation at the onset of the Intra-Allerød Cold Period. Based on bedform characteristics and fluvial morphology in the HSV, the maximum freshwater flux during the flood event is estimated to be ∼ 0.46 Sv for a duration of ∼ 80 days.

New Jersey, New York↗
USGS2023

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska↗
USGS2000

The 1999 Southern California Seismic Network bulletin

The Pasadena office of the U.S. Geological Survey (USGS), together with the Caltech Seismological Laboratory, operates a network of more than 350 remote seismometers in southern California called the S outhern C alifornia S eismic N etwork (SCSN). SCSN is part of TriNet, a cooperative project between the USGS, Caltech, and the California Division of Mines and Geology (CDMG). The TriNet project is halfway completed and is upgrading the existing network to digital, adding new stations, and developing real-time and earthquake-alert capabilities. Signals from the SCSN sites are telemetered to a central processing location at the Caltech Seismological Lab in Pasadena. Computers that detect and record thousands of earthquakes each year continuously monitor these signals. Phase arrival times for these events are picked by analysts and are archived along with digital seismograms. Data acquisition, processing, and archiving are achieved using the Caltech/USGS Seismic Processing (CUSP) system ( Dollar, 1989 ). These data have been compiled into the SCSN Catalog of Earthquakes, a list beginning in 1932 that currently contains more than 356,300 events. Waveform, phase, and catalog data are archived by the Southern California Earthquake Center Data Center (SCEC_DC). This data set is critical to the evaluation of earthquake hazards in California and to the advancement of geoscience as a whole.

California↗
USGS2015

Natural or Induced: Identifying Natural and Induced Swarms from Pre-production and Co-production Microseismic Catalogs at the Coso Geothermal Field

Increased levels of seismicity coinciding with injection of reservoir fluids have prompted interest in methods to distinguish induced from natural seismicity. Discrimination between induced and natural seismicity is especially difficult in areas that have high levels of natural seismicity, such as the geothermal fields at the Salton Sea and Coso, both in California. Both areas show swarm-like sequences that could be related to natural, deep fluid migration as part of the natural hydrothermal system. Therefore, swarms often have spatio-temporal patterns that resemble fluid-induced seismicity, and might possibly share other characteristics. The Coso Geothermal Field and its surroundings is one of the most seismically active areas in California with a large proportion of its activity occurring as seismic swarms. Here we analyze clustered seismicity in and surrounding the currently produced reservoir comparatively for pre-production and co-production periods. We perform a cluster analysis, based on the inter-event distance in a space-time-energy domain to identify notable earthquake sequences. For each event j, the closest previous event i is identified and their relationship categorized. If this nearest neighbor’s distance is below a threshold based on the local minimum of the bimodal distribution of nearest neighbor distances, then the event j is included in the cluster as a child to this parent event i. If it is above the threshold, event j begins a new cluster. This process identifies subsets of events whose nearest neighbor distances and relative timing qualify as a cluster as well as a characterizing the parent-child relationships among events in the cluster. We apply this method to three different catalogs: (1) a two-year microseismic survey of the Coso geothermal area that was acquired before exploration drilling in the area began; (2) the HYS_catalog_2013 that contains 52,000 double-difference relocated events and covers the years 1981 to 2013; and (3) a catalog of 57,000 events with absolute locations from the local network recorded between 2002 and 2007. Using this method we identify 10 clusters of more than 20 events each in the pre-production survey and more than 200 distinct seismicity clusters that each contain at least 20 and up to more than 1000 earthquakes in the more extensive catalogs. The cluster identification method used yields a hierarchy of links between multiple generations of parent and offspring events. We analyze different topological parameters of this hierarchy to better characterize and thus differentiate natural swarms from induced clustered seismicity and also to identify aftershock sequences of notable mainshocks. We find that the branching characteristic given by the average number of child events per parent event is significantly different for clusters below than for clusters around the produced field.

Conference Paper↗
USGS1975

Earthquake shaking and damage to buildings

Ground shaking close to the causative fault of an earthquake is more intense than it was previously believed to be. This raises the possibility that large numbers of buildings and other structures are not sufficiently resistant for the intense levels of shaking that can occur close to the fault. Many structures were built before earthquake codes were adopted; others were built according to codes formulated when less was known about the intensity of near-fault shaking. Although many building types are more resistant than conventional design analyses imply, the margin of safety is difficult to quantify. Many modern structures, such as freeways, have not been subjected to and tested by near-fault shaking in major earthquakes (magnitude 7 or greater). Damage patterns in recent moderate-sized earthquakes occurring in or adjacent to urbanized areas, however, indicate that many structures, including some modern ones designed to meet earthquake code requirements, cannot withstand the severe shaking that can occur close to a fault. It is necessary to review the ground motion assumed and the methods utilized in the design of important existing structures and, if necessary, to strengthen or modify the use of structures that are found to be weak. New structures situated close to active faults should be designed on the basis of ground motion estimates greater than those used in the past. The ultimate balance between risk of earthquake losses and cost for both remedial strengthening and improved earthquake-resistant construction must be decided by the public. Scientists and engineers must inform the public about earthquake shaking and its effect on structures. The exposure to damage from seismic shaking is steadily increasing because of continuing urbanization and the increasing complexity of lifeline systems, such as power, water, transportation, and communication systems. In the near future we should expect additional painful examples of the damage potential of moderate-sized earthquakes in urban areas. Over a longer time span, however, we can significantly reduce the risk to life and property from seismic shaking through better land utilization, improved building codes and construction practices, and at least the gradual replacement of poor buildings by more resistant buildings. Progress toward reducing risk from seismic shaking through better building design is slowed by deficiencies in our knowledge about the nature of damaging ground motion and the failure mechanisms in structures. For example, lacking observational data, seismologists must rely on simplified theoretical and numerical models of the earthquake process to estimate near-fault ground motion, especially for earthquakes as large as magnitude 7 and 8. Because such models have not been adequately tested against data, their reliability is unknown. Engineers lack detailed information about failure processes in structures during an earthquake. Although many structures have been instrumented to measure their response to an earthquake, few records have been obtained from buildings that actually sustained significant structural damage and few structures are properly instrumented to measure all the modes of deformation that are likely to contribute to failure. Moreover, the fact that many structures have withstood ground motion more intense than that assumed in their design indicates that conventional methods of design do not take into account important contributions to earthquake resistance by nonstructural elements and by the ability of structural elements to deform inelastically without necessarily causing failure of the structure. It is fortunate when such reserve resistance exists, but better understanding of the sources of reserve strength is needed to determine how large a margin of safety they confer and how they might be affected by changes in construction practices and materials with time. In the next few years we look forward to significant advances in knowledge and to more effective application of what is already known, largely because of substantial funding of research related to seismic engineering by the National Science Foundation . The increasing number of strong-motion seismographs operating in seismically active regions will likely provide for the first time a number of records of damaging levels of ground motion. Significant effort is being directed toward obtaining near-fault records, although many probable sites of future large earthquakes remain inadequately instrumented, especially outside the conterminous United States. New and more complete information on building response and damage mechanisms will be obtained by improved instrumentation of structures and through laboratory investigations of failure in structures and structural elements. Further developments in computer technology and in computer modeling techniques will permit more realistic simulations of the seismic response of soils and structures that take into account their inelastic behavior and their strain-dependent properties. Earthquake design codes will continually be revised to better utilize existing knowledge concerning the nature of strong ground motion and the dynamic behavior of buildings during earthquakes and to incorporate new knowledge and also experiences gained from future earthquakes. We believe that application of new knowledge, improvements in earthquake-resistant design and construction, and remedial strengthening or replacement of weak existing structures can significantly reduce our current level of exposure to earthquake hazards.

Science↗
USGS1998

Observations of hybrid seismic events at Soufriere Hills Volcano, Montserrat: July 1995 to September 1996

Swarms of small repetitive events with similar waveforms and magnitudes are often observed during the emplacement of lava domes. Over 300,000 such events were recorded in association with the emplacement of the lava dome at Soufriere Hills Volcano, Montserrat, from August 1995 through August 1996. These events originated < 2–3 km deep. They exhibited energy ranging over ˜ 1.5–4.5 Hz and were broader band than typical long-period events. We term the events “hybrid” between long-period and volcano-tectonic. The events were more impulsive and broader band prior to, compared with during and after, periods of inferred increased magma flux rate. Individual swarms contained up to 10,000 events often exhibiting very similar magnitudes and waveforms throughout the swarm. Swarms lasted hours to weeks, during which inter-event intervals generally increased, then decreased, often several times. Long-duration swarms began about every two months starting in late September 1995. We speculate that the events were produced as the magma column degassed into adjacent cracks.

Geophysical Research Letters↗
USGS2012

Anomalous chemical changes in well waters and possible relation to earthquakes

Water level, temperature, salinity, electric conductivity, and pH have been measured periodically for several years at three water wells located along a 17-km segment of the San Andreas fault between San Juan Bautista and Cienega Winery in central California. Water samples were collected at the same time for subsequent chemical analyses in the laboratory. Some sudden large changes in salinity and conductivity were recorded in early March 1980 at the two wells near San Juan Bautista. These changes coincided approximately with the beginning of an episode of increased local seismicity, including a magnitude 4.8 earthquake on April 13. Analyses of water samples revealed corresponding changes in ion concentrations, especially of Na + , Ca ++ , Mg ++ , SO 4 −− , HCO 3 − , F − , and Cl − . The observed changes may be the result of mixing of waters from different aquifers through cracks developed in the water barriers by a possible crustal strain episode that may have occurred in the study area.

Geophysical Research Letters↗
USGS2019

Characterizing large earthquakes before rupture is complete

Whether large and very large earthquakes are distinguishable from each other early on in the rupture process has been a subject often debated over the past several decades. Studies have shown that the frequency content of radiated seismic energy in the first few seconds of an earthquake scales with the final magnitude of the event, implying determinism. Other studies have shown that the recordings of ground displacement from small-to-moderate sized earthquakes are indistinguishable, and thus earthquakes share a universal early rupture process. Regardless of how earthquakes start, however, at some point in the rupture process events of different sizes must be distinguishable from one another. If that difference occurs early - i.e., before the rupture duration of the smaller event - this implies that earthquakes demonstrate some level of determinism. Here we show through analysis of a large database of source time functions and near-source displacement records that after an initiation phase, ruptures of M7-9 earthquakes organize into a slip pulse, the kinematic properties of which scale with magnitude. As such, early in the rupture process - after about 10s - large and very large earthquakes demonstrate different properties and can thus be distinguished.

Science Advances↗
USGS2021

Subducting oceanic basement roughness impacts on upper plate tectonic structure and a backstop splay fault zone activated in the southern Kodiak aftershock region of the Mw 9.2, 1964 megathrust rupture, Alaska

In 1964, the Alaska margin ruptured in a giant Mw 9.2 megathrust earthquake, the 2nd largest during worldwide instrumental recording. The coseismic slip and aftershock region offshore Kodiak Island was surveyed in 1977 – 1981 to understand the region’s tectonics. We re-processed multichannel seismic (MCS) field data using current standard Kirchhoff depth migration and/or MCS traveltime tomography. Further surveys in 1994 added P-wave velocity structure from wide-angle seismic lines and multibeam bathymetry. Published regional gravity, backscatter, and earthquake compilations also became available at this time. Beneath the trench, rough oceanic crust is covered by ~3 to 5 km thick sediment. Sediment on the subducting plate modulates the plate interface relief. The accreted prism’s imbricate thrust faults have a complex P-wave velocity structure. Landward, an accelerated increase in P-wave velocities is marked by a backstop splay fault zone (BSFZ) that marks a transition from the prism to the higher rigidity rock beneath the middle and upper slope. Structures associated with this feature may indicate fluid flow. Further upslope, another fault extends >100 km along-strike across the middle slope. Erosion from subducting seamounts leaves embayments in the frontal prism. Plate interface roughness varies along the subduction zone. Beneath the lower and middle slope, 2.5 D plate interface images show modest relief whereas the oceanic basement image is rougher. The 1964 earthquake slip maximum coincides with the leading/landward flank of a subducting seamount and the BSFZ. The BSFZ is a potentially active structure and should be considered in tsunami hazard assessments.

Alaska↗
USGS2010

Documentation for assessment of modal pushover-based scaling procedure for nonlinear response history analysis of "ordinary standard" bridges

The earthquake engineering profession is increasingly utilizing nonlinear response history analyses (RHA) to evaluate seismic performance of existing structures and proposed designs of new structures. One of the main ingredients of nonlinear RHA is a set of ground-motion records representing the expected hazard environment for the structure. When recorded motions do not exist (as is the case for the central United States), or when high-intensity records are needed (as is the case for San Francisco and Los Angeles), ground motions from other tectonically similar regions need to be selected and scaled. The modal-pushover-based scaling (MPS) procedure recently was developed to determine scale factors for a small number of records, such that the scaled records provide accurate and efficient estimates of 'true' median structural responses. The adjective 'accurate' refers to the discrepancy between the benchmark responses and those computed from the MPS procedure. The adjective 'efficient' refers to the record-to-record variability of responses. Herein, the accuracy and efficiency of the MPS procedure are evaluated by applying it to four types of existing 'ordinary standard' bridges typical of reinforced-concrete bridge construction in California. These bridges are the single-bent overpass, multi span bridge, curved-bridge, and skew-bridge. As compared to benchmark analyses of unscaled records using a larger catalog of ground motions, it is demonstrated that the MPS procedure provided an accurate estimate of the engineering demand parameters (EDPs) accompanied by significantly reduced record-to-record variability of the responses. Thus, the MPS procedure is a useful tool for scaling ground motions as input to nonlinear RHAs of 'ordinary standard' bridges.

Open-File Report↗
USGS2022

The dynamic floor of Yellowstone Lake, Wyoming, USA: The last 14 k.y. of hydrothermal explosions, venting, doming, and faulting

Hydrothermal explosions are significant potential hazards in Yellowstone National Park, Wyoming, USA. The northern Yellowstone Lake area hosts the three largest hydrothermal explosion craters known on Earth empowered by the highest heat flow values in Yellowstone and active seismicity and deformation. Geological and geochemical studies of eighteen sublacustrine cores provide the first detailed synthesis of the age, sedimentary facies, and origin of multiple hydrothermal explosion deposits. New tephrochronology and radiocarbon results provide a four-dimensional view of recent geologic activity since recession at ca. 15–14.5 ka of the >1-km-thick Pinedale ice sheet. The sedimentary record in Yellowstone Lake contains multiple hydrothermal explosion deposits ranging in age from ca. 13 ka to ~1860 CE. Hydrothermal explosions require a sudden drop in pressure resulting in rapid expansion of high-temperature fluids causing fragmentation, ejection, and crater formation; explosions may be initiated by seismicity, faulting, deformation, or rapid lake-level changes. Fallout and transport of ejecta produces distinct facies of subaqueous hydrothermal explosion deposits. Yellowstone hydrothermal systems are characterized by alkaline-Cl and/or vapor-dominated fluids that, respectively, produce alteration dominated by silica-smectite-chlorite or by kaolinite. Alkaline-Cl liquids flash to steam during hydrothermal explosions, producing much more energetic events than simple vapor expansion in vapor-dominated systems. Two enormous explosion events in Yellowstone Lake were triggered quite differently: Elliott’s Crater explosion resulted from a major seismic event (8 ka) that ruptured an impervious hydrothermal dome, whereas the Mary Bay explosion (13 ka) was triggered by a sudden drop in lake level stimulated by a seismic event, tsunami, and outlet channel erosion.

Wyoming↗
USGS1994

In search of earthquake-related hydrologic and chemical changes along Hayward Fault

Flow and chemical measurements have been made about once a month, and more frequently when required, since 1976 at two springs in Alum Rock Park in eastern San Jose, California, and since 1980 at two shallow wells in eastern Oakland in search of earthquake-related changes. All sites are on or near the Hayward Fault and are about 55 km apart. Temperature, electric conductivity, and water level or flow rate were measured in situ with portable instruments. Water samples were collected for later chemical and isotopic analyses in the laboratory. The measured flow rate at one of the springs showed a long-term decrease of about 40% since 1987, when a multi-year drought began in California. It also showed several increases that lasted a few days to a few months with amplitudes of 2.4 to 8.6 times the standard deviations above the background rate. Five of these increases were recorded shortly after nearby earthquakes of magnitude 5.0 or larger, and may have resulted from unclogging of the flow path and increase of permeability caused by strong seismic shaking. Two other flow increases were possibly induced by exceptionally heavy rainfalls. The water in both wells showed seasonal temperature and chemical variations, largely in response to rainfall. In 1980 the water also showed some clear chemical changes unrelated to rainfall that lasted a few months; these changes were followed by a magnitude 4 earthquake 37 km away. The chemical composition at one of the wells and at the springs also showed some longer-term variations that were not correlated with rainfall but possibly correlated with the five earthquakes mentioned above. These correlations suggest a common tectonic origin for the earthquakes and the anomalies. The last variation at the affected well occurred abruptly in 1989, shortly before a magnitude 5.0 earthquake 54 km away. ?? 1993.

California↗
USGS2000

SMSIM--Fortran programs for simulating ground motions from earthquakes: Version 2.0.--a revision of OFR 96-80-A

A simple and powerful method for simulating ground motions is based on the assumption that the amplitude of ground motion at a site can be specified in a deterministic way, with a random phase spectrum modified such that the motion is distributed over a duration related to the earthquake magnitude and to distance from the source. This method of simulating ground motions often goes by the name "the stochastic method." It is particularly useful for simulating the higher-frequency ground motions of most interest to engineers, and it is widely used to predict ground motions for regions of the world in which recordings of motion from damaging earthquakes are not available. This simple method has been successful in matching a variety of ground-motion measures for earthquakes with seismic moments spanning more than 12 orders of magnitude. One of the essential characteristics of the method is that it distills what is known about the various factors affecting ground motions (source, path, and site) into simple functional forms that can be used to predict ground motions. SMSIM is a set of programs for simulating ground motions based on the stochastic method. This Open-File Report is a revision of an earlier report (Boore, 1996) describing a set of programs for simulating ground motions from earthquakes. The programs are based on modifications I have made to the stochastic method first introduced by Hanks and McGuire (1981). The report contains source codes, written in Fortran, and executables that can be used on a PC. Programs are included both for time-domain and for random vibration simulations. In addition, programs are included to produce Fourier amplitude spectra for the models used in the simulations and to convert shear velocity vs. depth into frequency-dependent amplification. The revision to the previous report is needed because the input and output files have changed significantly, and a number of new programs have been included in the set.

Open-File Report↗
USGS2019

Geologic map of the Ferncliff and Louisa quadrangles, Louisa, Fluvanna, and Goochland Counties, Virginia

The area encompassed by the geologic map of the Ferncliff and Louisa, Va., 7.5-minute quadrangles includes the hypothetical surface projection of the Quail fault, which is the subsurface fault that was responsible for the 2011 magnitude 5.8 (M5.8) Mineral, Va., earthquake. The mapping shows that the Quail fault appears to have reactivated the Harris Creek fault, a Paleozoic fault that has been mapped and named in the study area and marks the boundary between the Ellisville pluton neck and Chopawamsic Formation. The Harris Creek fault was also reactivated in the early Mesozoic. Another result of the mapping is a well-defined, southwest to northeast, narrow zone of metagraywacke and ultramafic rocks (both part of the informal Shores complex) that marks the closure of a small ocean basin and the accretion of the 468- to 460-Ma (mega-annum) Ordovician Chopawamsic volcanic arc (part of the Carolina terrane) onto Laurentia. The accretion zone is truncated by the Ordovician-Silurian (444 Ma) Ellisville pluton; the 444-Ma age of the pluton therefore represents the minimum age of the accretion zone and indicates likely closure of the ocean basin during the Taconic orogeny. Across the map area, the metamorphic grade ranges from lower-greenschist facies in the northwest to amphibolite facies in the southeast. 40 Ar/ 39 Ar age-dating across this metamorphic gradient indicates that Taconic metamorphism was overprinted by late-Paleozoic Alleghanian metamorphism that was accompanied by refolding and faulting of Taconic structures. Quaternary terraces mapped along the South Anna River record a long history of incision and downcutting. A continuing question is how much of this downcutting was a result of neotectonic uplift in the central Virginia seismic zone. This report consists of a single geologic map sheet and an online geographic information systems database that includes bedrock geologic unit contacts and polygons, surficial geologic polygons, faults, and structural geologic information.

Virginia↗
USGS1969

A seismic-refraction survey of crustal structure in central Arizona

The U.S. Geological Survey conducted a seismic-refraction study of the earth's crust and upper mantle near the Tonto Forest Seismological Observatory (TFO), located 10miles south of the Mogollon Rim near Payson in central Arizona. Two recording lines 400 km long intersect in the approximate form of a cross at TFO; one line trends southeast and the other northeast. The sedimentary layer at most places southwest of the rim is less than 1 km thick, but north of the rim it is 2 to 3 km thick. The velocity in this uppermost layer ranges from 2.6 to 4.7 km/sec, with the higher limit measured near or north of the rim. Arrivals refracted in the upper crust (P g ) can be attributed to two layers for all the shot points south of the rim. The velocity in the upper layer is about 5.9 km/sec with thickness ranging from 2 to 8 km; beneath the upper layer the velocity is about 6.1 km/sec. The upper layer seems to be absent northeast of the rim, where two shot points generated P g arrivals that show only a velocity of 6.2 km/sec. A Poisson ratio of 0.22 for the upper crustal layers was measured from shear and compressional arrivals. The lower crust could not be identified from the first and later refraction arrivals; however, minimum depths to the intermediate layer were determined. An average crustal velocity of 6.2 km/sec was measured from wide-angle reflection alignments. A thin intermediate layer would explain the seismic measurements. A delay-time method was used to map the configuration of the M-discontinuity. The depth below sea level is about 36 km along the northwest-trending line. The northeast-trending line shows a shallow depth of 21 km near Gila Bend, increasing depth to about 34 km under TFO, and a flat M-discontinuity at 40 km depth under the Mogollon Mesa northeast to Sunrise Springs. There is evidence of an abrupt depth change of about 4 km on the M-discontinuity in the vicinity of TFO. The velocity in the upper mantle is 7.85 km/sec. The relation of topographic elevation to crustal thickness suggests an approach to isostatic equilibrium, which is deduced from a near-zero regional free-air gravity anomaly. However, lateral density change in the upper mantle is required to make the crustal-refraction model fit the observed gravity-anomaly values, provided that velocity and density are linearly related.

Arizona↗