SEARCH · USGS Science
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.
The U.S. Geological Survey is studying some of the features of the earth's crust in the San Juan Basin region, northwest New Mexico. As a part of this study seismic refraction-reflection measurements were made in and around the basin using explosions and earthquakes as energy sources. Record sections and traveltime tables were derived from the measurements.
Although situated between two States (California and Washington) that have has many violent earthquakes, Oregon is noticeably less active seismically. the greatest damage experienced resulted from a major shock near Olympia, Wash., in 1949. During the short history record available (since 1841), 34 earthquakes of intensity V, Modified Mercalli Scale, or greater have centered within Oregon or near its borders. Only 13 of the earthquakes had an intensity above V, and many of the shocks were local. However, a 1936 earthquake in the eastern Oregon-Washington region caused extensive damage and was felt over an area of 272,000 square kilometers.
We derive a detailed earthquake catalogue and V p , V s and V p / V s models for the region around the 2019 M w 6.4 and M w 7.1 Ridgecrest, California, earthquake sequence using data recorded by rapid-response, densely deployed sensors following the Ridgecrest main shock and the regional network. The new catalogue spans a 4-month period, starting on 1 June 2019, and it includes nearly 95 000 events detected and located with iterative updates to our velocity models. The final V p and V s models correlate well with surface geology in the top 4 km of the crust and spatial seismicity patterns at depth. Joint interpretation of the derived catalogue, velocity models, and surface geology suggests that (i) a compliant low-velocity zone near the Garlock Fault arrested the M w 7.1 rupture at the southeast end; (ii) a stiff high-velocity zone beneath the Coso Mountains acted as a strong barrier that arrested the rupture at the northwest end and (iii) isolated seismicity on the Garlock Fault accommodated transtensional-stepover strain triggered by the main events. The derived catalogue and velocity models can be useful for multiple future studies, including further analysis of seismicity patterns, derivations of accurate source properties (e.g. focal mechanisms) and simulations of earthquake processes and radiated seismic wavefields.
The midcontinent United States between the Appalachian and Rocky Mountains contains 40 known faults or other potentially tectonic features for which published geologic information shows or suggests Quaternary tectonic faulting. We report results of a systematic evaluation of published and other publicly available geologic evidence of Quaternary faulting. These results benefit seismic-hazard assessments by (1) providing some constraints on the recurrence intervals and magnitudes of large, prehistoric earthquakes, and (2) identifying features that warrant additional study. For some features, suggested Quaternary tectonic faulting has been disproved, whereas, for others, the suggested faulting remains questionable. Of the 40 features, nine have clear geologic evidence of Quaternary tectonic faulting associated with prehistoric earthquakes, and another six features have evidence of nontectonic origins. An additional 12 faults, uplifts, or historical seismic zones lack reported paleoseismological evidence of large. Quaternary earthquakes. The remaining 13 features require further paleoseismological study to determine if they have had Quaternary earthquakes that were larger than any known from local historical records; seven of these 13 features are in or near urbanized areas where their study could affect urban hazard estimates. These seven are: (1) the belt of normal faults that rings the Gulf of Mexico from Florida to Texas. (2) the Northeast Ohio seismic zone, (3) the Valmont and (4) Goodpasture faults of Colorado. (5) the Champlain lowlands normal faults of New York State and Vermont, and (6) the Lexington and (7) Kentucky River fault systems of eastern Kentucky. Published by Elsevier Science B.V.
The 2009 eruption of Redoubt produced several very-long-period (VLP) signals associated with explosions. We invert for the source location and mechanism of an explosion at Redoubt volcano using waveform methods applied to broadband recordings. Such characterization of the source carries information on the geometry of the conduit and the physics of the explosion process. Inversions are carried out assuming the volcanic source can be modeled as a point source, with mechanisms described by a) a set of 3 orthogonal forces, b) a moment tensor consisting of force couples, and c) both forces and moment tensor components. We find that the source of the VLP seismic waves during the explosion is well-described by either a combined moment/force source located northeast of the crater and at an elevation of 1.6 km ASL or a moment source at an elevation of 800 m to the southwest of the crater. The moment tensors for the solutions with moment and force and moment-only share similar characteristics. The source time functions for both moment tensors begin with inflation (pressurization) and execute two cycles of deflation-reinflation (depressurization–repressurization). Although the moment/force source provides a better fit to the data, we find that owing to the limited coverage of the broadband stations at Redoubt the moment-only source is the more robust and reliable solution. Based on the moment-only solution, we estimate a volume change of 19,000 m 3 and a pressure change of 7 MPa in a dominant sill and an out-of-phase volume change of 5000 m 3 and pressure change of 1.8 MPa in a subdominant dike at the source location. These results shed new light on the magmatic plumbing system beneath Redoubt and complement previous studies on Vulcanian explosions at other volcanoes.
Southwestern British Columbia and northwestern Washington State are subject to megathrust earthquakes, deep intraslab events, and earthquakes in the continental crust. Of the three types of earthquakes, the most poorly understood are the crustal events. Despite a high level of seismicity, there is no obvious correlation between the historical crustal earthquakes and the mapped surface faults of the region. On 24 June 1997, a ML = 4.6 earthquake occurred 3-4 km beneath the Strait of Georgia, 30 km to the west of Vancouver, British Columbia. This well-recorded earthquake was preceded by 11 days by a felt foreshock (ML = 3.4) and was followed by numerous small aftershocks. This earthquake sequence occurred in one of the few regions of persistent shallow seismic activity in southwestern British Columbia, thus providing an ideal opportunity to attempt to characterize an active near-surface fault. We have computed focal mechanisms and utilized a waveform cross-correlation and joint hypocentral determination routine to obtain accurate relative hypocenters of the mainshock, foreshock, and 53 small aftershocks in an attempt to image the active fault and the extent of rupture associated with this earthquake sequence. Both P-nodal and CMT focal mechanisms show thrust faulting for the mainshock and the foreshock. The relocated hypocenters delineate a north-dipping plane at 2-4 km depth, dipping at 53??, in good agreement with the focal mechanism nodal plane dipping to the north at 47??. The rupture area is estimated to be a 1.3-km-diameter circular area, comparable to that estimated using a Brune rupture model with the estimated seismic moment of 3.17 ?? 1015 N m and the stress drop of 45 bars. The temporal sequence indicates a downdip migration of the seismicity along the fault plane. The results of this study provide the first unambiguous evidence for the orientation and sense of motion for active faulting in the Georgia Strait area of British Columbia.
Mount Rainier has had at least 11 large lahars over the last 6,000 years, including one occurring without evidence of eruptive activity. This prompted the creation of a lahar detection system that uses a combination of seismic, infrasound, and tripwires. We test a laser rangefinder placed on a river channel bank for detecting and confirming mass movements flowing past a station as an alternative to the physical tripwires. After testing the device at an experimental debris-flow flume, the laser rangefinder successfully captured a small debris flow on Mount Rainier in 2023, confirming its effectiveness as a lahar detection and monitoring tool. Over the 2-month deployment at Mount Rainier, we find that spurious recordings in the laser rangefinder data (noise) tend to correlate with high humidity, and that periods of noise do not correlate with increased co-located seismic amplitude. Therefore, the impact of the noise on future alarms can be mitigated by coupling a laser rangefinder alarm with that of independent datasets.
We utilize a two-step process to validate 0-4 Hz ground motion simulations using the 1989 Loma Prieta earthquake. In the first step we run multiple realizations using the Graves and Pitarka hybrid method as implemented on the SCEC Broadband Simulation Platform and compare these with near-fault (R < 40 km) recorded motions. A total of 648 rupture scenarios are examined and from these results we select the best fitting ruptures, which are then used in the second step of the validation. This second step consists of running 3D simulations using a finely sampled seismic velocity mesh. The base 3D structure is obtained from the USGS SF Bay Area velocity model, which is also modified to include small-scale 3D stochastic perturbations. The 3D results are compared with the recorded motions and show very good agreement over the frequency band 0-4 Hz for distances out to 40 km. Most importantly, the newly added features to the simulation process reduce the coherency of the radiated higher frequency (f > 1 Hz) ground motions, and homogenize radiation-pattern effects in this same bandwidth, bringing the simulations into closer agreement with the very near-fault records.
The Calaveras fault lies within a low velocity zone (LVZ) 1-2 km wide near Gilroy, California. Accelerographs G06, located in the LVZ 1.2 km from the Calaveras fault, and G07, 4 km from G06, recorded both the M 6.2 1984 Morgan Hill and the M 6.9 1989 Loma Prieta earthquakes. Comparison of the ground motions shows that a large 0.6-1.0 Hz velocity pulse observed at G06 during the Morgan Hill event may be amplified by focussing caused by the LVZ. Such amplified waves might be a mappable seismic hazard, and the zone of increased hazard can extend as much as 1.2 km from the surface trace of the fault. Finite-difference simulations of ground motions in a simplified LVZ model show a zone of amplified motion similar to the observations.
Blind reverse faults are challenging to detect, and earthquake records can be elusive because deep fault slip does not break the surface along readily recognized scarps. The blind Reelfoot fault in the New Madrid seismic zone in the central United States has been the subject of extensive prior investigation; however, the extent of slip at the southern portion of the fault remains unconstrained. In this study, we use lidar to map terraces and lacustrine landforms in the Obion River valley and investigate apparent broad folding resulting from slip on the buried Reelfoot fault. We compare remote surface mapping results with three auger boreholes in the ∼24 ka Finley terrace and interpret apparent warping as due to tectonic folding and not stratigraphic thickening. We combine our results with historical records of coseismic lake formation that indicate surface deformation dammed the Obion River in the 1812 CE earthquake. Older terraces (deposited at least 35–55 ka) record progressive fold scarps ≥1, ≥2, and ≥8 m high indicating a long record of earthquakes predating the existing paleoseismic record. Broad, distributed folding above the Reelfoot fault into the Obion River valley is consistent with a deep active fault tip along the southern reaches of the fault. Our analyses indicate the entire length of the fault (≥70 km) is capable of rupture and is more consistent with longer rupture scenarios.
The 3 November 2002 Mw 7.9 Denali fault earthquake triggered deformational offsets and microseismicity under Mammoth Mountain (MM) on the rim of Long Valley caldera, California, some 3460 km from the earthquake. Such strain offsets and microseismicity were not recorded at other borehole strain sites along the San Andreas fault system in California. The Long Valley offsets were recorded on borehole strainmeters at three sites around the western part of the caldera that includes Mammoth Mountain - a young volcano on the southwestern rim of the caldera. The largest recorded strain offsets were -0.1 microstrain at PO on the west side of MM, 0.05 microstrain at MX to the southeast of MM, and -0.025 microstrain at BS to the northeast of MM with negative strain extensional. High sample rate strain data show initial triggering of the offsets began at 22:30 UTC during the arrival of the first Rayleigh waves from the Alaskan earthquake with peak-to-peak dynamic strain amplitudes of about 2 microstrain corresponding to a stress amplitude of about 0.06 MPa. The strain offsets grew to their final values in the next 10 min. The associated triggered seismicity occurred beneath the south flank of MM and also began at 22:30 UTC and died away over the next 15 min. This relatively weak seismicity burst included some 60 small events with magnitude all less than M = 1. While poorly constrained, these strain observations are consistent with triggered slip and intrusive opening on a north-striking normal fault centered at a depth of 8 km with a moment of l016 N m, or the equivalent of a M 4.3 earthquake. The cumulative seismic moment for the associated seismicity burst was more than three orders of magnitude smaller. These observations and this model resemble those for the triggered deformation and slip that occurred beneath the north side of MM following the 16 October 1999 M 7.1 Hector Mine, California, earthquake. However, in this case, we see little post-event slip decay reflected in the strain data after the Rayleigh-wave arrivals from the Denali fault earthquake and onset of triggered seismicity did not lag the triggered deformation by 20 min. These observations are also distinctly different from the more widespread and energetic seismicity and deformation triggered by the 1992 M 7.3 Landers earthquake in the Long Valley caldera. Thus, each of the three instances of remotely triggered unrest in Long Valley caldera recorded to date differ. In each case, however, the deformation moment inferred from the strain meter data was more than an order of magnitude larger than the cumulative moment for the associated triggered seismicity.
A reconnaissance shallow seismic reflection/refraction investigation in and around the city of Galena, Alaska suggests that Tertiary and/or Cretaceous bedrock, and possible coalbeds within the Cretaceous, is at least as deep as 550 feet in the immediate vicinity of town. Rock could be deeper than 1000 feet under alternate interpretations. Reflections recorded in these data are believed to be from the sediment/bedrock interface. Analysis of these reflections and associated refractions indicates that this interface, interpreted at most of the six profile locations, has a high seismic velocity, possibly indicating non-sedimentary rock (e.g. volcanic or igneous).
An episode of unrest began at Kīlauea in April 2018 that produced both significant volcanic output and high rates of seismicity, including a M w 6.9 earthquake on 4 May 2018. In this study, we image the rupture process of this earthquake using a genetic algorithm-based back-projection technique. The dominant feature of the earthquake is a slowly propagating western rupture, which shares similar characteristics with the region's largest recorded event in 1975 ( M w 7.7). The location of this western segment suggests that small asperities on this section of the décollement that frequently fail as slow slip events may achieve seismic slip rates when rupture is initiated on adjacent sections of the fault. Given the interaction between volcanic and seismic activity in this region, imaging the rupture properties of these events can improve our understanding of future geologic hazards in this region.
A 180-km-long seismic refraction transect from the eastern High Cascades, across Newberry Volcano, to the eastern High Lava Plains is used to investigate the subvolcanic crustal and upper mantle velocity structure there. Near-surface volcanic flows and sedimentary debris (1.6-4.7 km/s), ranging from 3 to 5 km in thickness, overlie subvolcanic Basin and Range structures. East and west of Newberry Volcano, the subvolcanic basement (5.6 km/s) has been downwarped, producing 5-km-deep basins. The midcrust (8- to 28-km depth) is characterized by velocities ranging from 6.1 to 6.5 km/s and varies laterally in thicknesses. The lower crust is characterized by an unusually high velocity (about 7.4 km/s), and its geometry mirrors the subvolcanic basement geometry. The Moho is located at a depth of 37 km and represents a transition to an upper mantle velocity of 8.1 km/s. The shallow subsurface (1.2 km) beneath Newberry Volcano is characterized by high-velocity 5.6 km/s, versus 4.1 km/s for the surrounding area) intrusions and appears to be located on a basement high. Beneath the seismic refraction array at Newberry Volcano, an absence of low-velocity anomalies suggests that large silicic magma chambers do not exist in the upper crust, but apparent high attenuation of the seismic wave field may be consistent with either partial melts in small volumes, elevated crustal temperatures, and/or poor geophone-recording site coupling.
The observed ground motions from five large aftershocks of the 1999 Chi-Chi, Taiwan, earthquake are compared with predictions from four equations based primarily on data from California. The four equations for active tectonic regions are those developed by Abrahamson and Silva (1997), Boore et al. (1997), Campbell (1997, 2001), and Sadigh et al. (1997). Comparisons are made for horizontal-component peak ground accelerations and 5%-damped pseudoacceleration response spectra at periods between 0.02 sec and 5 sec. The observed motions are in reasonable agreement with the predictions, particularly for distances from 10 to 30 km. This is in marked contrast to the motions from the Chi-Chi mainshock, which are much lower than the predicted motions for periods less than about 1 sec. The results indicate that the low motions in the mainshock are not due to unusual, localized absorption of seismic energy, because waves from the mainshock and the aftershocks generally traverse the same section of the crust and are recorded at the same stations. The aftershock motions at distances of 30-60 km are somewhat lower than the predictions (but not nearly by as small a factor as those for the mainshock), suggesting that the ground motion attenuates more rapidly in this region of Taiwan than it does in the areas we compare with it. We provide equations for the regional attenuation of response spectra, which show increasing decay of motion with distance for decreasing oscillator periods. This observational study also demonstrates that ground motions have large earthquake-location-dependent variability for a specific site. This variability reduces the accuracy with which an earthquake-specific prediction of site response can be predicted. Online Material: PGAs and PSAs from the 1999 Chi-Chi earthquake and five aftershocks.
To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( M w ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.
The area designated for possible oil and gas lease sale as modified from BLM memorandum 3310 #42 (722) and referred to therein as the North Atlantic Outer Continental Shelf (OCS) contains about 58,300 sq km of shelf beneath water depths of less than 200 m and lies chiefly within the Georges Bank basin. The oldest sediments drilled or dredged on the bordering Continental Slope are sandstone, clay, and silt of Upper Cretaceous age. In Upper Cretaceous exposures, on Marthas Vineyard and nearby New England islands, the predominant lithology appears to be clay. About 125 km northeast of the eastern tip of Georges Bank, the Shell B-93 well penetrated clays and silts of Upper and Lower Cretaceous age above dense Jurassic carbonate rocks which overlie a basement of lower Paleozoic slate, schist, quartzite, and granite. Structurally, the Georges Bank basin is a westerly trending trough which opens to the west-southwest. Post-Paleozoic sediments are more than 8 km thick in parts of the basin. Major structural features appear to be directly related to basement structures. Local anticlines, probably caused by differential compaction over basement flexures and horst blocks or by later uplift along basement faults are reflected principally in Lower Cretaceous and older sediments, though some of these features continue upward to within 0.7 of a second (about 650 m) of the seafloor. Tertiary deposits in the Georges Bank basin are probably up to a kilometre thick and are made up of poorly consolidated sand, silt, and clay. The Cretaceous section is inferred to be up to 3.5 km thick and to be mainly clastics -- shale, siltstone, calcareous shale, changing to limestone in the lowest part of the system. Jurassic rocks in the deepest part of the basin appear to be about 3.6 to 4.0 km thick and probably consist mainly of dense carbonates. Potential source rocks in the Georges Bank basin may include organic-rich Cretaceous shale and carbonaceous Jurassic limestone. By analogy with the Scotian Shelf, Cretaceous sandstones are considered to be potential reservoir rocks. Local zones of porous dolomite are believed to be present in carbonate rocks of Jurassic age and should not be overlooked as potential reservoirs. Structural highs related to draping and differential compaction over basement blocks could be important potential petroleum traps. Additional traps may include reef structures near the shelf edge, updip pinchouts, and stratigraphic traps in both clastic and carbonate sediment. A statistical mean for the undiscovered recoverable petroleum resources is calculated to be 0.9 billion barrels of oil and 4.2 trillion cubic feet of gas. At the 5 percent probability level (1 in 20 chance) the undiscovered recoverable petroleum resources are calculated to be 2.4 billion barrels of oil and 12.5 trillion cubic feet of gas. These undiscovered recoverable petroleum resources are those quantities of oil and gas that may be reasonably expected to exist in favorable settings, but which have not yet been identified by drilling. Such estimates, therefore, carry a high degree of uncertainty. Environmental studies of Georges Bank indicated a low-moderate risk from petroleum development. However, the risk estimate is based on very limited data. Drift bottle returns used to infer oil spill trajectories show about a 2% overall recovery rate. Meteorologic data comes mainly from nearby land areas and from ships attempting to avoid storms. Seismicity on Georges Bank is low. This may reflect, in large part, the difficulty of land-based stations in recording earthquakes far from the coast. Direct data on the engineering properties of shallow buried sediment comes mainly from two Texas Tower surveys of limited areas on Georges Shoal and Nantucket Shoals made in the early 1950's. The 17 holes (most less than 30 m deep) reveal some silty layers below loose sand and much lateral variability in sediment type over short distances. The technology for exploration at the required water depths (20 m - 200 m) is available. Mobil drilling units are in great demand around the world and will have to be brought in from other areas along with skilled manpower. Our highest estimates indicate 50 platforms, 800 producing wells, 1,100 km of pipeline, and 5 onshore terminals may be needed. The time frame for production, using our high estimates (5% probability) for the undiscovered recoverable resources, could include 4-5 years for significant development, 6-7 years until production commences, and 18 years until peak production.
Although the probability of a major earthquake in the central U.S. is only one-tenth of the Californian events, any such event is expected to result in damage ten times greater than that anticipated in California. The city of Memphis, which is situated very close to the inferred epicenter of one of the three major 1811-1812 earthquakes, is in a potentially hazardous zone which will be susceptible to the usual seismic hazards. By recognizing the high level of seismicity in the New Madrid area, this study attempts to microzone the potential hazards in the city of Memphis. This study examines the pertinent criteria necessary for any microzonation study. The subjects considered include: (i) The seismicity of the central United States, (ii) Design earthquakes, (iii) Response analysis which allows us to construct the necessary microzonation maps. The seismicity of the region is evaluated according to the state-of-the-art literature available as there is no recorded strong motion data available for the central U.S. The maximum credible earthquakes which are likely to affect Memphis are evaluated as the one in a thousand year occurrence. However, earthquakes of a lower intensity (and lower recurrence rates) are selected as design earthquakes to permit a more realistic microzonation to be performed. For the response analysis, these earthquakes are simulated using synthetically generated accelerograms which exhibit the relevant parameters. These accelerograms displaying the significant characteristics of amplitude, predominant frequency and duration were selected to display, as accurately as possible, the anticipated nature of the horizontal bedrock-motions at Memphis. The horizontal motions (SH-waves) were applied at a depth of 45m below ground surface at numerous sites in the city of Memphis, where the soil stratigraphy had been conceptualized from borehole data. The soils data made available by local sources was the only information used for this study as any laboratory testing was outside the scope of this study. The dynamic soil properties were thus established from the available Standard Penetration Resistances and soil classifications. The results of the response analysis were transformed into microzonation maps depicting (i) zones showing qualitative estimates of ground response, (ii) zones showing the natural frequency of the soils, (iii) zones showing the peak spectral acceleration for 2% damping ratio, (iv) zones of liquefaction potential. These maps are essentially useful for preliminary investigation and we do not expect them to be used on a quantitative basis. However, further investigation is necessary in determining the stratigraphy and soil properties to a more accurate level if an earthquake hazard is anticipated.