USGS Science⌕ Search

SEARCH · USGS Science

Results for “Tectonics”

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,387 records · Page 77Linked to original sources
USGS2002

The geohydrologic setting of Yucca Mountain, Nevada

This paper provides a geologic and hydrologic framework of the Yucca Mountain region for the geochemical papers in this volume. The regional geologic units, which range in age from late Precambrian through Holocene, are briefly described. Yucca Mountain is composed of dominantly pyroclastic units that range in age from 11.4 to 15.2 Ma. The principal focus of study has been on the Paintbrush Group, which includes two major zoned and welded ash-flow tuffs separated by an important hydrogeologic unit referred to as the Paintbrush non-welded (PTn). The regional structural setting is currently one of extension, and the major local tectonic domains are presented together with a tectonic model that is consistent with the known structures at Yucca Mountain. Streamflow in this arid to semi-arid region occurs principally in intermittent or ephemeral channels. Near Yucca Mountain, the channels of Fortymile Wash and Amargosa River collect infrequent runoff from tributary basins, ultimately draining to Death Valley. Beneath the surface, large-scale interbasin flow of groundwater from one valley to another occurs commonly in the region. Regional groundwater flow beneath Yucca Mountain originates in the high mesas to the north and returns to the surface either in southern Amargosa Desert or in Death Valley, where it is consumed by evapotranspiration. The water table is very deep beneath the upland areas such as Yucca Mountain, where it is 500-750 m below the land surface, providing a large thickness of unsaturated rocks that are potentially suitable to host a nuclear-waste repository. The nature of unsaturated flow processes, which are important for assessing radionuclide migration, are inferred mainly from hydrochemical or isotopic evidence, from pneumatic tests of the fracture systems, and from the results of in situ experiments. Water seeping down through the unsaturated zone flows rapidly through fractures and more slowly through the pores of the rock matrix. Although capillary forces are expected to divert much of the flow around repository openings, some may drip onto waste packages, ultimately causing release of radionuclides, followed by transport down to the water table. ?? 2002 Elsevier Science Ltd. All rights reserved.

Applied Geochemistry↗
USGS2016

Mississippian coral latitudinal diversity gradients (western interior United States): Testing the limits of high resolution diversity data

Analysis of high resolution diversity data for Mississippian corals in the western interior United States yielded mild latitudinal diversity gradients despite the small geographic area covered by samples and a large influence on diversity patterns by geographic sampling intensity (sample bias). Three competing plate tectonic reconstructions were tested using the diversity patterns. Although none could be forcefully rejected, one reconstruction proved less consistent with diversity patterns than the other two and additional coral diversity data from farther north in Canada would better discriminate the two equivalent reconstructions. Despite the relatively high sampling intensity represented by the analyzed database, diversity patterns were greatly affected by sample abundance and distribution. Hence, some effort at recognizing and accounting for sample bias should be undertaken in any study of latitudinal diversity gradients. Small-scale geographic lumping of sample localities had only small effects on geographic diversity patterns. However, large-scale (e.g., regional) geographic lumping of diversity data may not yield latitudinally sensitive diversity patterns. Temporal changes in coral diversity in this region reflect changes in eustacy, local tectonism, and terrigenous sediment flux, far more than they do shifting latitude. Highest regional diversity occurred during the interval when the studied region occupied the highest latitude. Therefore, diversity data from different regions may not be comparable, in terms of latitudinal inference. Small-scale stratigraphic lumping of the data caused a nearly complete loss of the latitudinal diversity patterns apparent prior to lumping. Hence, the narrowest possible stratigraphic resolution should be maintained in analyzing latitudinal diversity gradients.

Journal of Paleontology↗
USGS1997

The rupture process of the Manjil, Iran earthquake of 20 June 1990 and implications for intraplate strike-slip earthquakes

In terms of seismically radiated energy or moment release, the earthquake of 20 January 1990 in the Manjil Basin-Alborz Mountain region of Iran is the second largest strike-slip earthquake to have occurred in an intracontinental setting in the past decade. It caused enormous loss of life and the virtual destruction of several cities. Despite a very large meizoseismal area, the identification of the causative faults has been hampered by the lack of reliable earthquake locations and conflicting field reports of surface displacement. Using broadband data from global networks of digitally recording seismographs, we analyse broadband seismic waveforms to derive characteristics of the rupture process. Complexities in waveforms generated by the earthquake indicate that the main shock consisted of a tiny precursory subevent followed in the next 20 seconds by a series of four major subevents with depths ranging from 10 to 15 km. The focal mechanisms of the major subevents, which are predominantly strike-slip, have a common nodal plane striking about 285°–295°. Based on the coincidence of this strike with the dominant tectonic fabric of the region we presume that the EW striking planes are the fault planes. The first major subevent nucleated slightly south of the initial precursor. The second subevent occurred northwest of the initial precursor. The last two subevents moved progressively southeastward of the first subevent in a direction collinear with the predominant strike of the fault planes. The offsets in the relative locations and the temporal delays of the rupture subevents indicate heterogeneous distribution of fracture strength and the involvement of multiple faults. The spatial distribution of teleseismic aftershocks, which at first appears uncorrelated with meizoseismal contours, can be decomposed into stages. The initial activity, being within and on the periphery of the rupture zone, correlates in shape and length with meizoseismal lines. In the second stage of activity the aftershock zone expands and appears to cluster about the geomorphic and geologic features several tens of kilometres from the rupture zone. The activity is interpreted as a regional response to quasistatic stress migration along zones of tectonic weakness. The radiated energy of the main shock and the estimate of seismic moment yields an apparent stress of 20 bars. High apparent stress may be typical of strike slip earthquakes occurring in intracontinental environments undergoing continental collision.

Studia Geophysica et Geodaetica↗
USGS2013

Paleoseismicity and neotectonics of the Aleutian subduction zone — An overview

The Aleutian subduction zone is one of the most seismically active plate boundaries and the source of several of the world’s largest historic earthquakes. The structural architecture of the subduction zone varies considerably along its length. At the eastern end is a tectonically complex collision zone where the allochthonous Yakutat terrane is moving northwest into mainland Alaska. West of the collision zone a shallow-dipping subducted plate beneath a wide forearc, nearly orthogonal convergence, and a continental-type subduction regime characterizes the eastern part of the subduction zone. In the central part of the subduction zone, convergence becomes increasingly right oblique and the forearc is divided into a series of large clockwise-rotated fault-bounded blocks. Highly oblique convergence and island arc tectonics characterize the western part of the subduction zone. At the extreme western end of the arc, the relative plate motion is nearly pure strike-slip. A series of great subduction earthquakes ruptured most of the 4000-km length of the subduction zone during a period of several decades in the mid 1900s. The majority of these earthquakes broke multiple segments as defined by the large-scale structure of the overriding plate margin and patterns of historic seismicity. Several of these earthquakes generated Pacific-wide tsunamis and significant damage in the southwestern and south-central regions of Alaska. Characterization of previous subduction earthquakes is important in assessing future seismic and tsunami hazards. However, at present such information is available only for the eastern part of the subduction zone. The 1964 Alaska earthquake (M 9.2) ruptured about ~950 km of the plate boundary that encompassed the Kodiak and Prince William Sound (PWS) segments. Within this region, nine paleosubduction earthquakes in the past ~5000 years are recognized on the basis of geologic evidence of sudden land level change and, at some sites, coeval tsunami deposits. Carbon 14-based chronologies indicate recurrence intervals between median calibrated ages for these paleoearthquakes range from 333 to 875 years. The most recent occurred about 489 years ago and broke only the Kodiak segment. During the previous three cycles, both the Kodiak and PWS segments were involved in either multiple-segment ruptures or closely timed pairs of single segment ruptures. evidence for the earlier paleosubduction earthquakes has been found only at sites in the PWS segment. Thus, future work on the paleoseismicity of other segments would by particular valuable in defining the seismic behavior of the subduction zone.

Alaska↗
USGS2013

Does a boundary of the Wrangell Block extend through southern Cook Inlet and Shelikof Strait, Alaska?

In southcentral Alaska, the boundaries of two different tectonic blocks extend southwestward from the Denali Fault toward Cook Inlet and Shelikof Strait. We use offshore multichannel seismic reflection data and oil-well stratigraphy to evaluate whether local geologic structures are compatible with boundaries of either tectonic block and with the relative motion expected across the block boundaries. Our main conclusion is that a block boundary does not extend southwestward the entire length of Shelikof Strait, as was proposed for one of the blocks. Furthermore, below southern Cook Inlet, no high-strain extensional structures that might be related to either proposed boundary are evident. Small normal faults below southern Cook Inlet could have been caused by block rotation, but they represent only minor strain. One way to explain the lack of larger structures is that the rotation began recently so that indicative boundary structures have not yet formed. Alternatively, deformation associated with the block boundaries could be distributed through onshore areas.

Alaska↗
USGS1999

On relating apparent stress to the stress causing earthquake fault slip

Apparent stress τ a is defined as , where is the average shear stress loading the fault plane to cause slip and η is the seismic efficiency, defined as E a / W , where E a is the energy radiated seismically and W is the total energy released by the earthquake. The results of a recent study in which apparent stresses of mining-induced earthquakes were compared to those measured for laboratory stick-slip friction events led to the hypothesis that . This hypothesis is tested here against a substantially augmented data set of earthquakes for which can be estimated, mostly from in situ stress measurements, for comparison with τ a . The expanded data set, which includes earthquakes artificially triggered at a depth of 9 km in the German Kontinentales Tiefbohrprogramm der Bundesrepublik Deutschland (KTB) borehole and natural tectonic earthquakes, covers a broad range of hypocentral depths, rock types, pore pressures, and tectonic settings. Nonetheless, over ∼14 orders of magnitude in seismic moment, apparent stresses exhibit distinct upper bounds defined by a maximum seismic efficiency of ∼0.06, consistent with the hypothesis proposed before. This behavior of τ a and η can be expressed in terms of two parameters measured for stick-slip friction events in the laboratory: the ratio of the static to the dynamic coefficient of friction and the fault slip overshoot. Typical values for these two parameters yield seismic efficiencies of ∼0.06. In contrast to efficiencies for laboratory events for which η is always near 0.06, those for earthquakes tend to be less than this bounding value because E a for earthquakes is usually underestimated due to factors such as band-limited recording. Thus upper bounds on appear to be controlled by just a few fundamental aspects of frictional stick-slip behavior that are common to shallow earthquakes everywhere. Estimates of from measurements of τ a for suites of earthquakes, using , are found to be comparable in magnitude to estimates of shear stress on the basis of extrapolating in situ stress data to seismogenic depths.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Alternating seismic uplift and subsidence in the late Holocene at Madang, Papua New Guinea: Evidence from raised reefs

Well-preserved mid-late Holocene coral reefs are exposed in low coastal cliffs in the vicinity of the Madang lagoon on the north coast of Papua New Guinea. Results from U/Th and 14 C dating of corals, surveying, and field mapping indicate several major changes in relative sea level over this period. Specifically, there is evidence for a relative sea level fall of ≥4.5 m about 3000 calendar years B.P., followed by relative sea level rises of ∼1.5 m about 2400 calendar years B.P. and ≥0.5 m about 1200 calendar years B.P. and a subsequent relative sea level fall of ≥3 m some time in the past 1000 years. Since regional eustatic sea levels are believed to have been dropping gradually over this time frame, these observed changes in relative sea level are interpreted as reflecting alternating tectonic uplift and subsidence. Furthermore, the detailed structure and age relationships of the coral deposits indicate that both uplift and subsidence occurred rapidly, most probably as coseismic events with vertical displacements of 0.5 to 4.5 m. These events may be related to rupture on NW-SE trending reverse faults which have been mapped in the nearby Adelbert Range and possibly on NE trending cross faults which have been inferred from seismicity. This interpretation implies a much greater degree of tectonic instability and potential seismic hazard in the region than previously recognized, although the inferred coseismic vertical displacements are shown to be consistent with present-day local seismicity. In a broader context, the study illustrates how detailed analysis of vertical changes in coral reef structure and assemblages may be used as a sensitive indicator of changing relative sea level, capable of resolving century timescale events and reversals.

Madang↗
USGS2000

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Ground deformation associated with the March 1996 earthquake swarm at Akutan volcano, Alaska, revealed by satellite radar interferometry

In March 1996 an intense swarm of volcano-tectonic earthquakes (???3000 felt by local residents, Mmax = 5.1, cumulative moment of 2.7 ??1018 N m) beneath Akutan Island in the Aleutian volcanic arc, Alaska, produced extensive ground cracks but no eruption of Akutan volcano. Synthetic aperture radar interferograms that span the time of the swarm reveal complex island-wide deformation: the western part of the island including Akutan volcano moved upward, while the eastern part moved downward. The axis of the deformation approximately aligns with new ground cracks on the western part of the island and with Holocene normal faults that were reactivated during the swarm on the eastern part of the island. The axis is also roughly parallel to the direction of greatest compressional stress in the region. No ground movements greater than 2.83 cm were observed outside the volcano's summit caldera for periods of 4 years before or 2 years after the swarm. We modeled the deformation primarily as the emplacement of a shallow, east-west trending, north dipping dike plus inflation of a deep, Mogi-type magma body beneath the volcano. The pattern of subsidence on the eastern part of the island is poorly constrained. It might have been produced by extensional tectonic strain that both reactivated preexisting faults on the eastern part of the island and facilitated magma movement beneath the western part. Alternatively, magma intrusion beneath the volcano might have been the cause of extension and subsidence in the eastern part of the island. We attribute localized subsidence in an area of active fumaroles within the Akutan caldera, by as much as 10 cm during 1992-1993 and 1996-1998, to fluid withdrawal or depressurization of the shallow hydrothermal system. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses

We have applied a previously developed analytical stress model to interpret subsurface stress conditions inferred from acoustic televiewer logs obtained in two municipal water wells located in a valley in the southern Davis Mountains near Alpine, Texas. The appearance of stress-induced breakouts with orientations that shift by 90° at two different depths in one of the wells is explained by results from exact solutions for the effects of valleys on gravity and tectonically induced subsurface stresses. The theoretical results demonstrate that above a reference depth termed the hinge point, a location that is dependent on Poisson's ratio, valley shape, and magnitude of the maximum horizontal tectonic stress normal to the long axis of the valley, horizontal stresses parallel to the valley axis are greater than those normal to it. At depths below this hinge point the situation reverses and horizontal stresses normal to the valley axis are greater than those parallel to it. Application of the theoretical model at Alpine is accommodated by the fact that nearby earthquake focal mechanisms establish an extensional stress regime with the regional maximum horizontal principal stress aligned perpendicular to the valley axis. We conclude that the localized stress field associated with a valley setting can be highly variable and that breakouts need to be examined in this context when estimating the orientations and magnitudes of regional principal stresses.

Texas↗
USGS2006

Seismic characteristics of central Brazil crust and upper mantle: A deep seismic refraction study

[1] A two‐dimensional model of the Brazilian central crust and upper mantle was obtained from the traveltime interpretation of deep seismic refraction data from the Porangatu and Cavalcante lines, each approximately 300 km long. When the lines were deployed, they overlapped by 50 km, forming an E‐W transect approximately 530 km long across the Tocantins Province and western São Francisco Craton. The Tocantins Province formed during the Neoproterozoic when the São Francisco, the Paranapanema, and the Amazon cratons collided, following the subduction of the former Goiás ocean basin. Average crustal V P and V P /V S ratios, Moho topography, and lateral discontinuities within crustal layers suggest that the crust beneath central Brazil can be associated with major geological domains recognized at the surface. The Moho is an irregular interface, between 36 and 44 km deep, that shows evidences of first‐order tectonic structures. The 8.05 and 8.23 km s −1 P wave velocities identify the upper mantle beneath the Porangatu and Cavalcante lines, respectively. The observed seismic features allow for the identification of (1) the crust has largely felsic composition in the studied region, (2) the absence of the mafic‐ultramafic root beneath the Goiás magmatic arc, and (3) block tectonics in the foreland fold‐and‐thrust belt of the northern Brasília Belt during the Neoproterozoic. Seismic data also suggested that the Bouguer gravimetric discontinuities are mainly compensated by differences in mass distribution within the lithospheric mantle. Finally, the Goiás‐Tocantins seismic belt can be interpreted as a natural seismic alignment related to the Neoproterozoic mantle domain.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Modeling the evolution of channel shape: Balancing computational efficiency with hydraulic fidelity

The cross-sectional shape of a natural river channel controls the capacity of the system to carry water off a landscape, to convey sediment derived from hillslopes, and to erode its bed and banks. Numerical models that describe the response of a landscape to changes in climate or tectonics therefore require formulations that can accommodate evolution of channel cross-sectional geometry. However, fully two-dimensional (2-D) flow models are too computationally expensive to implement in large-scale landscape evolution models, while available simple empirical relationships between width and discharge do not adequately capture the dynamics of channel adjustment. We have developed a simplified 2-D numerical model of channel evolution in a cohesive, detachment-limited substrate subject to steady, unidirectional flow. Erosion is assumed to be proportional to boundary shear stress, which is calculated using an approximation of the flow field in which log-velocity profiles are assumed to apply along vectors that are perpendicular to the local channel bed. Model predictions of the velocity structure, peak boundary shear stress, and equilibrium channel shape compare well with predictions of a more sophisticated but more computationally demanding ray-isovel model. For example, the mean velocities computed by the two models are consistent to within ???3%, and the predicted peak shear stress is consistent to within ???7%. Furthermore, the shear stress distributions predicted by our model compare favorably with available laboratory measurements for prescribed channel shapes. A modification to our simplified code in which the flow includes a high-velocity core allows the model to be extended to estimate shear stress distributions in channels with large width-to-depth ratios. Our model is efficient enough to incorporate into large-scale landscape evolution codes and can be used to examine how channels adjust both cross-sectional shape and slope in response to tectonic and climatic forcing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2009

Anatomy of the dead sea transform from lithospheric to microscopic scale

Fault zones are the locations where motion of tectonic plates, often associated with earthquakes, is accommodated. Despite a rapid increase in the understanding of faults in the last decades, our knowledge of their geometry, petrophysical properties, and controlling processes remains incomplete. The central questions addressed here in our study of the Dead Sea Transform (DST) in the Middle East are as follows: (1) What are the structure and kinematics of a large fault zone? (2) What controls its structure and kinematics? (3) How does the DST compare to other plate boundary fault zones? The DST has accommodated a total of 105 km of leftlateral transform motion between the African and Arabian plates since early Miocene (???20 Ma). The DST segment between the Dead Sea and the Red Sea, called the Arava/ Araba Fault (AF), is studied here using a multidisciplinary and multiscale approach from the ??m to the plate tectonic scale. We observe that under the DST a narrow, subvertical zone cuts through crust and lithosphere. First, from west to east the crustal thickness increases smoothly from 26 to 39 km, and a subhorizontal lower crustal reflector is detected east of the AF. Second, several faults exist in the upper crust in a 40 km wide zone centered on the AF, but none have kilometer-size zones of decreased seismic velocities or zones of high electrical conductivities in the upper crust expected for large damage zones. Third, the AF is the main branch of the DST system, even though it has accommodated only a part (up to 60 km) of the overall 105 km of sinistral plate motion. Fourth, the AF acts as a barrier to fluids to a depth of 4 km, and the lithology changes abruptly across it. Fifth, in the top few hundred meters of the AF a locally transpressional regime is observed in a 100-300 m wide zone of deformed and displaced material, bordered by subparallel faults forming a positive flower structure. Other segments of the AF have a transtensional character with small pull-aparts along them. The damage zones of the individual faults are only 5-20 m wide at this depth range. Sixth, two areas on the AF show mesoscale to microscale faulting and veining in limestone sequences with faulting depths between 2 and 5 km. Seventh, fluids in the AF are carried downward into the fault zone. Only a minor fraction of fluids is derived from ascending hydrothermal fluids. However, we found that on the kilometer scale the AF does not act as an important fluid conduit. Most of these findings are corroborated using thermomechanical modeling where shear deformation in the upper crust is localized in one or two major faults; at larger depth, shear deformation occurs in a 20-40 km wide zone with a mechanically weak decoupling zone extending subvertically through the entire lithosphere. Copyright 2009 by the American Geophysical Union.

Reviews of Geophysics↗
USGS2010

Glimpses of East Antarctica: Aeromagnetic and satellite magnetic view from the central Transantarctic Mountains of East Antarctica

Aeromagnetic and satellite magnetic data provide glimpses of the crustal architecture within the Ross Sea sector of the enigmatic, ice-covered East Antarctic shield critical for understanding both global tectonic and climate history. In the central Transantarctic Mountains (CTAM), exposures of Precambrian basement, coupled with new high-resolution magnetic data, other recent aeromagnetic transects, and satellite magnetic and seismic tomography data, show that the shield in this region comprises an Archean craton modified both by Proterozoic magmatism and early Paleozoic orogenic basement reactivation. CTAM basement structures linked to the Ross Orogeny are imaged 50–100 km farther west than previously mapped, bounded by inboard upper crustal Proterozoic granites of the Nimrod igneous province. Magnetic contrasts between craton and rift margin sediments define the Neoproterozoic rift margin, likely reactivated during Ross orogenesis and Jurassic extension. Interpretation of satellite magnetic and aeromagnetic patterns suggests that the Neoproterozoic rift margin of East Antarctica is offset by transfer zones to form a stepwise series of salients tracing from the CTAM northward through the western margin of the Wilkes Subglacial Basin to the coast at Terre Adélie. Thinned Precambrian crust inferred to lie east of the rift margin cannot be imaged magnetically because of modification by Neoproterozoic and younger tectonic events.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Connecting the Yakima fold and thrust belt to active faults in the Puget Lowland, Washington

High-resolution aeromagnetic surveys of the Cascade Range and Yakima fold and thrust belt (YFTB), Washington, provide insights on tectonic connections between forearc and back-arc regions of the Cascadia convergent margin. Magnetic surveys were measured at a nominal altitude of 250 m above terrain and along flight lines spaced 400 m apart. Upper crustal rocks in this region have diverse magnetic properties, ranging from highly magnetic rocks of the Miocene Columbia River Basalt Group to weakly magnetic sedimentary rocks of various ages. These distinctive magnetic properties permit mapping of important faults and folds from exposures to covered areas. Magnetic lineaments correspond with mapped Quaternary faults and with scarps identified in lidar (light detection and ranging) topographic data and aerial photography. A two-dimensional model of the northwest striking Umtanum Ridge fault zone, based on magnetic and gravity data and constrained by geologic mapping and three deep wells, suggests that thrust faults extend through the Tertiary section and into underlying pre-Tertiary basement. Excavation of two trenches across a prominent scarp at the base of Umtanum Ridge uncovered evidence for bending moment faulting possibly caused by a blind thrust. Using aeromagnetic, gravity, and paleoseismic evidence, we postulate possible tectonic connections between the YFTB in eastern Washington and active faults of the Puget Lowland. We suggest that faults and folds of Umtanum Ridge extend northwestward through the Cascade Range and merge with the Southern Whidbey Island and Seattle faults near Snoqualmie Pass 35 km east of Seattle. Recent earthquakes (M W ≤ 5.3) suggest that this confluence of faults may be seismically active today.

Washington↗
USGS2019

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth↗
USGS2020

Structural controls on slope failure within the western Santa Barbara Channel based on 2D and 3D seismic imaging

The Santa Barbara Channel, offshore California, contains several submarine landslides and ample evidence for incipient failure. This region hosts active thrust and reverse faults that accommodate several mm/yr of convergence, yet the relationships between tectonic deformation and slope failure remain unclear. We present 3‐D and 2‐D multichannel seismic reflection (MCS) data sets, multibeam bathymetry, and chronostratigraphic constraints to investigate the controls on slope failure. Splay faulting along the North Channel Deformation Trend (NCDT) coincides with a distinct zone of compressional uplift and onlapping of steeply dipping Quaternary strata. The NCDT is spatially correlated with seafloor fissures, and 3‐D seismic analyses reveal an intricate system of en echelon reverse faults that offset sediments younger than ~25 ka. Localized uplift zones are located between faults, one of which underlies the Gaviota landslide headscarp. We observe a direct relationship between slope failure and along‐strike variations in the tectonostratigraphic framework. Based on geophysical properties at Ocean Drilling Program (ODP) Site 893, we predict a trend in compaction and porosity reduction in the basin that drives pore fluids up‐dip, toward the zone of onlap above the NCDT, thus reducing slope stability. This interplay between tectonic, sedimentary, and fluid‐flow processes along the NCDT has created a confluence of preconditioning factors, with Gaviota and Goleta landslides being distinguished from the surrounding slopes by their position above the NCDT. The distribution of seafloor fissures suggests sections of the slope remain unstable and are prone to future landsliding. These results provide insights into the processes and 3‐D feedbacks that lead to slope instability along other convergent margins.

California↗
USGS2021

Indicators of volcanic eruptions revealed by global M4+ earthquakes

Determining whether seismicity near volcanoes is due primarily to tectonic or magmatic processes is a challenging but critical endeavor for volcanic eruption forecasting and detection, especially at poorly monitored volcanoes. Global statistics on the occurrence and timing of earthquakes near volcanoes both within and outside of eruptive periods reveal patterns in eruptive seismicity that may improve our ability to discern magmatically driven seismicity from purely tectonic seismicity. In this paper, we catalog magnitude four and greater (M4+) earthquakes near volcanoes globally and compute statistics on their occurrence with respect to various eruptive and volcanic attributes, evaluating their utility as diagnostic indicators of eruptions. Using a 2‐week time window and a 30 km radius around the volcanoes, we find that 11% of eruptions are preceded by at least one M4+ earthquake, but only 1% of such earthquakes is followed by eruption. However, earthquakes located 5–15 km from the volcano, those with normal faulting mechanisms and/or large nondouble‐couple components, and those occurring as groups are more commonly associated with eruptions, providing significant forecasting utility in some cases. Similarly, certain volcanoes are more likely to exhibit such precursors, such as those with long repose periods. We illustrate the use of these data in eruption forecasting scenarios, including rapid identification of analogous earthquake sequences at other volcanoes. When integrated within the context of multiparametric, multidisciplinary probabilistic assessments of volcanic activity, global earthquake statistics can improve eruption forecasts, and our work provides a model for use on other rapidly expanding global volcanological databases.

Journal of Geophysical Research↗