Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 379 records · Page 21Linked to original sources

Hydrogeologic controls imposed by mechanical stratigraphy in layered rocks of the Chateauguay River Basin, a U.S.-Canada transborder aquifer

[1] The Châteauguay River Basin delineates a transborder watershed with roughly half of its surface area located in northern New York State and half in southern Québec Province, Canada. As part of a multidisciplinary study designed to characterize the hydrogeologic properties of this basin, geophysical logs were obtained in 12 wells strategically located to penetrate the four major sedimentary rock formations that constitute the regional aquifers. The layered rocks were classified according to their elastic properties into three primary units: soft sandstone, hard sandstone, and dolostone. Downhole measurements were analyzed to identify fracture patterns associated with each unit and to evaluate their role in controlling groundwater flow. Fracture networks are composed of orthogonal sets of laterally extensive, subhorizontal bedding plane partings and bed-delimited, subvertical joints with spacings that are consistent with rock mechanics principles and stress models. The vertical distribution of transmissive zones is confined to a few select bedding plane fractures, with soft sandstone having the fewest (one per 70-m depth) and hard sandstone the most (five per 70-m depth). Bed-normal permeability is examined using a probabilistic model that considers the lengths of flow paths winding along joints and bedding plane fractures. Soft sandstone has the smallest bed-normal permeability primarily because of its wide, geomechanically undersaturated joint spacing. Results indicate that the three formations have similar values of bulk transmissivity, within roughly an order of magnitude, but that each rock unit has its own unique system of groundwater flow paths that constitute that transmissivity.

Chateauguay River Basin↗

Elastic moduli of rock glasses under pressure to 8 kilobars and geophysical implications

Shear and longitudinal velocities were measured by the ultrasonic phase comparison method as a function of pressure to 8 kbar on synthetic glasses of basalt, andesite, rhyolite, and quartz composition and on natural obsidian. Velocities of most of the glasses decrease anomalously with pressure, but increasingly more-normal behavior occurs with decrease in SiO 2 content. The pressure derivatives of rigidity and bulk modulus increase linearly, from −3.39 to −0.26 and from −5.91 to +2.09, respectively, with decrease in SiO 2 content from 100 to 49%. The change from negative to positive in the pressure derivatives of both moduli and observed at Poisson's ratio of about 0.25 is consistent with the Smyth model for the anomalous elastic behavior of glass. If the temperature in the upper mantle is about 1500°C, tholeiitic basalt would be molten in accordance with the partial melt explanation for the low-velocity zone; at 1300°C and below, basalt would be in the glassy state, especially if more felsic than tholeiite. At a temperature of 1370°C and at 30 kbar, reasonable values for the upper mantle at 100 km depth, the basalt glass of this study would have a viscosity of about 10 13 P. On the basis of the theory of viscoelasticity the glass would support shear wave propagation at frequencies above 0.1 Hz. Under such conditions of PT, 10 to 30% basalt glass in a matrix either of eclogite or olivine would reduce the seismic velocities by 3 to 9% and could also account for the values observed in the low-velocity zone.

Journal of Geophysical Research Solid Earth↗

Fracture characterization by means of attenuation and generation of tube waves in fractured crystalline rock at Mirror Lake, New Hampshire

Results are presented from experiments carried out in conjunction with the U. S. Geological Survey at the Hubbard Brook Experimental Forest near Mirror Lake, New Hampshire. The study focuses on our ability to obtain orientation and transmissivity estimates of naturally occurring fractures. The collected data set includes a four-offset hydrophone vertical seismic profile, full waveform acoustic logs at 5, 15, and 34 kHz, borehole televiewer, temperature, resistivity, and self-potential logs, and borehole-to-borehole pump test data. Borehole televiewer and other geophysical logs indicate that permeable fractures intersect the Mirror Lake boreholes at numerous depths, but less than half of these fractures appear to have significant permeability beyond the annulus of drilling disturbance on the basis of acoustic waveform log analysis. The vertical seismic profiling (VSP) data indicate a single major permeable fracture near a depth of 44 m, corresponding to one of the most permeable fractures identified in the acoustic waveform log analysis. VSP data also indicate a somewhat less permeable fracture at 220 m and possible fractures at depths of 103 and 135 m; all correspond to major permeable fractures in the acoustic waveform data set. Pump test data confirm the presence of a hydraulic connection between the Mirror Lake boreholes through a shallow dipping zone of permeability at 44 m in depth. Effective fracture apertures calculated from modeled transmissivities correspond to those estimated for the largest fractures indicated on acoustic waveform logs but are over an order of magnitude larger than effective apertures calculated from tube waves in the VSP data set. This discrepancy is attributed to the effect of fracture stiffness. A new model is presented to account for the mechanical strength of asperities in resisting fracture closure during the passage of seismic waves during the generation of VSPs.

New Hampshire↗

Bayesian inversion of data from effusive volcanic eruptions using physics-based models: Application to Mount St. Helens 2004--2008

Physics-based models of volcanic eruptions can directly link magmatic processes with diverse, time-varying geophysical observations, and when used in an inverse procedure make it possible to bring all available information to bear on estimating properties of the volcanic system. We develop a technique for inverting geodetic, extrusive flux, and other types of data using a physics-based model of an effusive silicic volcanic eruption to estimate the geometry, pressure, depth, and volatile content of a magma chamber, and properties of the conduit linking the chamber to the surface. A Bayesian inverse formulation makes it possible to easily incorporate independent information into the inversion, such as petrologic estimates of melt water content, and yields probabilistic estimates for model parameters and other properties of the volcano. Probability distributions are sampled using a Markov-Chain Monte Carlo algorithm. We apply the technique using GPS and extrusion data from the 2004–2008 eruption of Mount St. Helens. In contrast to more traditional inversions such as those involving geodetic data alone in combination with kinematic forward models, this technique is able to provide constraint on properties of the magma, including its volatile content, and on the absolute volume and pressure of the magma chamber. Results suggest a large chamber of >40 km 3 with a centroid depth of 11–18 km and a dissolved water content at the top of the chamber of 2.6–4.9 wt%.

Washington↗

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗

Combined seismic and geodetic analysis before, during and after the 2018 Mt. Etna eruption

In December 2018, Etna volcano experienced one of the largest episodes of unrest since the installation of geophysical monitoring networks in 1970. The unrest culminated in a short eruption with a small volume of lava erupted, a significant seismic crisis and deformation of the entire volcanic edifice of magnitude never recorded before at Mount Etna. Here we describe the evolution of the 2018 eruptive cycle from the analysis of seismic and geodetic data collected in the months preceding, during, and following the intrusion. We model the space-time evolution of high-rate deformation data starting from the active source previously identified from deformation data and the propagation of seismicity in a 3-D velocity model. The intrusion model suggests emplacement of two dikes: a smaller dike located beneath the eruptive fissure and a second, deeper dike between 1 and 5 km below sea level that opened ~2 m. The rise and eruption of magma from the shallower dike did not interrupt the pressurization of a long-lasting deeper reservoir (~6 km) that induced continuous inflation and intense deformation of the eastern flank. Shortly after the intrusion, on 26 December 2018, a M L 4.8 earthquake occurred near Pisano, destroying buildings and roads in two villages. We propose a time-dependent intrusion model that supports the hypothesis of the inflation inducing flank deformation and that this process has been active since September 2018.

Mt. Etna↗

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

Grain-size segregation and levee formation in geophysical mass flows

Data from large-scale debris-flow experiments are combined with modeling of particle-size segregation to explain the formation of lateral levees enriched in coarse grains. The experimental flows consisted of 10 m 3 of water-saturated sand and gravel, which traveled ∼80 m down a steeply inclined flume before forming an elongated leveed deposit 10 m long on a nearly horizontal runout surface. We measured the surface velocity field and observed the sequence of deposition by seeding tracers onto the flow surface and tracking them in video footage. Levees formed by progressive downslope accretion approximately 3.5 m behind the flow front, which advanced steadily at ∼2 m s −1 during most of the runout. Segregation was measured by placing ∼600 coarse tracer pebbles on the bed, which, when entrained into the flow, segregated upwards at ∼6–7.5 cm s −1 . When excavated from the deposit these were distributed in a horseshoe-shaped pattern that became increasingly elevated closer to the deposit termination. Although there was clear evidence for inverse grading during the flow, transect sampling revealed that the resulting leveed deposit was strongly graded laterally, with only weak vertical grading. We construct an empirical, three-dimensional velocity field resembling the experimental observations, and use this with a particle-size segregation model to predict the segregation and transport of material through the flow. We infer that coarse material segregates to the flow surface and is transported to the flow front by shear. Within the flow head, coarse material is overridden, then recirculates in spiral trajectories due to size-segregation, before being advected to the flow edges and deposited to form coarse-particle-enriched levees.

Journal of Geophysical Research↗

The hydrothermal system at Newberry Volcano, Oregon

Results of recent geological and geophysical studies at Newberry Volcano have been incorporated into conceptual and numerical models of a magma-based hydrothermal system. Numerical simulations begin with emplacement of a small magma body, the presumed source of silicic eruptions at Newberry that began about 10,000 B.P., into a thermal regime representing 100,000 years of cooling of a large underlying intrusion. Simulated flow patterns and thermal histories for three sets of hypothetical permeability values are compatible with data from four geothermal drill holes on the volcano. Meteoric recharge cools the caldera-fill deposits, but thermal water moving up a central conduit representing a permeable volcanic vent produces temperatures close to those observed in drill holes within the caldera. Meteoric recharge from the caldera moves down the flanks and creates a near-isothermal zone that extends several hundred meters below the water table, producing temperature profiles similar to those observed in drill holes on the flanks. The temperatures observed in drill holes on the flanks are not influenced by the postulated Holocene magma body. The elevated temperature gradients measured in the lower portions of these holes may be related to the cumulative effect of older intrusions. The models also indicate that meteoric recharge to the deep hydrothermal system probably originates within or near the caldera. Relatively low fluid velocities at depth suggest that at least a significant fraction of the thermal fluid may be very old.

Journal of Geophysical Research Solid Earth↗

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗

A continuous record of intereruption velocity change at Mount St. Helens from coda wave interferometry

In September 2004, Mount St. Helens volcano erupted after nearly 18&thinsp;years of quiescence. However, it is unclear from the limited geophysical observations when or if the magma chamber replenished following the 1980&ndash;1986 eruptions in the years before the 2004&ndash;2008 extrusive eruption. We use coda wave interferometry with repeating earthquakes to measure small changes in the velocity structure of Mount St. Helens volcano that might indicate magmatic intrusion. By combining observations of relative velocity changes from many closely located earthquake sources, we solve for a continuous function of velocity changes with time. We find that seasonal effects dominate the relative velocity changes. Seismicity rates and repeating earthquake occurrence also vary seasonally; therefore, velocity changes and seismicity are likely modulated by snow loading, fluid saturation, and/or changes in groundwater level. We estimate hydrologic effects impart stress changes on the order of tens of kilopascals within the upper 4&thinsp;km, resulting in annual velocity variations of 0.5 to 1%. The largest nonseasonal change is a decrease in velocity at the time of the deep M w &thinsp;=&thinsp;6.8 Nisqually earthquake. We find no systematic velocity changes during the most likely times of intrusions, consistent with a lack of observable surface deformation. We conclude that if replenishing intrusions occurred, they did not alter seismic velocities where this technique is sensitive due to either their small size or the finite compressibility of the magma chamber. We interpret the observed velocity changes and shallow seasonal seismicity as a response to small stress changes in a shallow, pressurized system.

Journal of Geophysical Research B: Solid Earth↗

Constraining the magmatic system at Mount St. Helens (2004–2008) using Bayesian inversion with physics-based models including gas escape and crystallization

Physics-based models of volcanic eruptions track conduit processes as functions of depth and time. When used in inversions, these models permit integration of diverse geological and geophysical data sets to constrain important parameters of magmatic systems. We develop a 1-D steady state conduit model for effusive eruptions including equilibrium crystallization and gas transport through the conduit and compare with the quasi-steady dome growth phase of Mount St. Helens in 2005. Viscosity increase resulting from pressure-dependent crystallization leads to a natural transition from viscous flow to frictional sliding on the conduit margin. Erupted mass flux depends strongly on wall rock and magma permeabilities due to their impact on magma density. Including both lateral and vertical gas transport reveals competing effects that produce nonmonotonic behavior in the mass flux when increasing magma permeability. Using this physics-based model in a Bayesian inversion, we link data sets from Mount St. Helens such as extrusion flux and earthquake depths with petrological data to estimate unknown model parameters, including magma chamber pressure and water content, magma permeability constants, conduit radius, and friction along the conduit walls. Even with this relatively simple model and limited data, we obtain improved constraints on important model parameters. We find that the magma chamber had low (<5wt%) total volatiles and that the magma permeability scale is well constrained at ~10-11.4 m2 to reproduce observed dome rock porosities. Compared with previous results, higher magma overpressure and lower wall friction are required to compensate for increased viscous resistance while keeping extrusion rate at the observed value.

Washington↗

Topographic stress perturbations in southern Davis Mountains, west Texas 2. Hydrogeologic implications

As part of a regional groundwater investigation, geophysical logs were obtained in two municipal water wells located near the west Texas city of Alpine. These boreholes are 252 and 285 m deep and penetrate extrusive rocks of Tertiary age. The deeper well was drilled in the central valley and the other along the northern flank of an east-west trending valley-ridge setting. Analysis and interpretation of the logs reveal that the two wells are subjected to significantly different stress environments because of topographic effects and exhibit significantly different hydrogeologic properties. Water production is associated with two specific types of features common to both wells: (1) the upper and lower contacts of a dense trachyte unit located in the shallow part of the wells and (2) deeper zones of highly fractured rocks within the interior of a basalt formation. The transmissivity of the trachyte boundaries is twice as large in the central valley well as it is in the ridge flank well, whereas the transmissivity of the deeper basalts is an order of magnitude greater in the flank well than it is in the central well. This discrepancy is examined from the perspective of rock failure, fracture opening, and flow enhancement by computing values for a Drucker-Prager stability factor that is based on the magnitudes of the normal and deviatoric stress invariants as a function of depth. Thus the field measurements and subsequent stress analysis offer evidence of a coupled tectonic-hydrologic interaction at this site.

Texas↗

Tectonic history of the north portion of the San Andreas fault system, California, inferred from gravity and magnetic anomalies

Geologic and geophysical data for the San Andreas fault system north of San Francisco suggest that the eastern boundary of the Pacific plate migrated eastward from its presumed original position at the base of the continental slope to its present position along the San Andreas transform fault by means of a series of eastward jumps of the Mendocino triple junction. These eastward jumps total a distance of about 150 km since 29 Ma. Correlation of right-laterally displaced gravity and magnetic anomalies that now have components at San Francisco and on the shelf north of Point Arena indicates that the presently active strand of the San Andreas fault north of the San Francisco peninsula formed recently at about 5 Ma when the triple junction jumped eastward a minimum of 100 km to its present location at the north end of the San Andreas fault. Prior to 5 Ma the triple junction was located at the north end of a proposed northwesterly extension of the Pilarcitos fault. This jump took place at what is now about latitude 38°20′N on the North American plate and may have occurred during part of the time when the relative motion between the North American and Pacific plates rotated 20° clockwise, resulting in transpression along the earlier transcurrent San Andreas fault system. The proposed 150-km eastward movement of the triple junction explains the submarine topography near Cape Mendocino where the continental shelf south of the Mendocino fault extends about 130 km farther west than does the shelf directly north of the fault.

Journal of Geophysical Research Solid Earth↗

Imaging the magmatic system of Mono Basin, California with magnetotellurics in three--dimensions

A three&ndash;dimensional (3D) electrical resistivity model of Mono Basin in eastern California unveils a complex subsurface filled with zones of partial melt, fluid&ndash;filled fracture networks, cold plutons, and regional faults. In 2013, 62 broadband magnetotelluric (MT) stations were collected in an array around southeastern Mono Basin from which a 3D electrical resistivity model was created with a resolvable depth of 35 km. Multiple robust electrical resistivity features were found that correlate with existing geophysical observations. The most robust features are two 300 &plusmn; 50 km3 near-vertical conductive bodies (3&ndash;10 Ω&middot;m) that underlie the southeast and north-eastern margin of Mono Craters below 10 km depth. These features are interpreted as magmatic crystal&ndash;melt mush zones of 15 &plusmn; 5% interstitial melt surrounded by hydrothermal fluids and are likely sources for Holocene eruptions. Two conductive east&ndash;dipping structures appear to connect each magma source region to the surface. A conductive arc&ndash;like structure (< 0.9 Ω&middot;m) links the northernmost mush column at 10 km depth to just below vents near Panum Crater, where the high conductivity suggests the presence of hydrothermal fluids. The connection from the southernmost mush column at 10 km depth to below South Coul&eacute;e is less obvious with higher resistivity (200 Ω&middot;m) suggestive of a cooled connection. A third, less constrained conductive feature (4&ndash;10 Ω&middot;m) 15 km deep extending to 35 km is located west of Mono Craters near the eastern front of the Sierra Nevada escarpment, and is coincident with a zone of sporadic, long&ndash;period earthquakes that are characteristic of a fluid-filled (magmatic or metamorphic) fracture network. A resistive feature (103&ndash;105 Ω&middot;m) located under Aeolian Buttes contains a deep root down to 25 km. The eastern edge of this resistor appears to structurally control the arcuate shape of Mono Craters. These observations have been combined to form a new conceptual model of the magmatic system beneath Mono Craters to a depth of 30 km.

California↗

Geophysical evidence for the evolution of the California Inner Continental Borderland as a metamorphic core complex

We use new seismic and gravity data collected during the 1994 Los Angeles Region Seismic Experiment (LARSE) to discuss the origin of the California Inner Continental Borderland (ICB) as an extended terrain possibly in a metamorphic core complex mode. The data provide detailed crustal structure of the Borderland and its transition to mainland southern California. Using tomographic inversion as well as traditional forward ray tracing to model the wide-angle seismic data, we find little or no sediments, low (≤6.6 km/s) P wave velocity extending down to the crust-mantle boundary, and a thin crust (19 to 23 km thick). Coincident multichannel seismic reflection data show a reflective lower crust under Catalina Ridge. Contrary to other parts of coastal California, we do not find evidence for an underplated fossil oceanic layer at the base of the crust. Coincident gravity data suggest an abrupt increase in crustal thickness under the shelf edge, which represents the transition to the western Transverse Ranges. On the shelf the Palos Verdes Fault merges downward into a landward dipping surface which separates "basement" from low-velocity sediments, but interpretation of this surface as a detachment fault is inconclusive. The seismic velocity structure is interpreted to represent Catalina Schist rocks extending from top to bottom of the crust. This interpretation is compatible with a model for the origin of the ICB as an autochthonous formerly hot highly extended region that was filled with the exhumed metamorphic rocks. The basin and ridge topography and the protracted volcanism probably represent continued extension as a wide rift until ∼13 m.y. ago. Subduction of the young and hot Monterey and Arguello microplates under the Continental Borderland, followed by rotation and translation of the western Transverse Ranges, may have provided the necessary thermomechanical conditions for this extension and crustal inflow.

California↗

The role of farfield tectonic stress in oceanic intraplate deformation, Gulf of Alaska

An integration of geophysical data from the Pacific Plate reveals plate bending anomalies, massive intraplate shearing and deformation, and a lack of oceanic crust magnetic lineaments in different regions across the Gulf of Alaska. We argue that farfield stress from the Yakutat Terrane collision with North America is the major driver for these unusual features. Similar plate motion vectors indicate that the Pacific plate and Yakutat Terrane are largely coupled along their boundary, the Transition Fault, with minimal translation. Our study shows that the Pacific Plate subduction angle shallows toward the Yakutat Terrane and supports the theory that the Pacific Plate and Yakutat Terranemaintain coupling along the subducted region of the Transition Fault. We argue that the outboard transfer of collisional stress to the Pacific Plate could have resulted in significant strain in the NE corner of the Pacific Plate, which created pathways for igneous sill formation just above the Pacific Plate crust in the Surveyor Fan. A shift in Pacific Plate motion during the late Miocene altered the Yakutat collision with North America, changing the stress transfer regime and potentially terminating associated strain in the NE corner of the Pacific Plate. The collision further intensified as the thickest portion of the Yakutat Terrane began to subduct during the Pleistocene, possibly providing the impetus for the creation of the Gulf of Alaska Shear Zone, a>200 km zone of intraplate strike-slip faults that extend from the Transition Fault out into the Pacific Plate. This study highlights the importance of farfield stress from complex tectonic regimes in consideration of large-scale oceanic intraplate deformation.

Gulf of Alaska↗

Crustal structure and composition of the southern Foothills Metamorphic Belt, Sierra Nevada, California, from seismic data

The Foothills Metamorphic Belt is an accreted terrane consisting of Paleozoic and Mesozoic metamorphic rocks that separates the Great Valley from the Sierra Nevada batholith in northern and central California. Until recently, the only available geophysical data for this area were reconnaissance refraction surveys, and gravity and magnetic data. New insights into the structure of the deep crust are provided by the interpretation of a seismic reflection profile (CC‐2), acquired in 1984 by the U.S. Geological Survey at the southern end of the Foothills Metamorphic Belt. Our interpretation is constrained by a new seismic velocity model derived from coincident microearthquake data. Earthquake hypocenters that occur at unusually great depths of 12 to 30 km make the data set particularly useful for obtaining deep crustal velocity information. The velocity model shows velocities of 5.2 to 6.3 km s −1 for the upper 12 km of the crust, and 6.7 to 6.8 km s −1 from 12 km to an estimated Moho at 32 km. The upper crustal velocities correspond to metamorphic rocks and serpentinites of the Foothills Metamorphic Belt as well as to diorites and granodiorites of the Sierra Nevada batholith, while the lower crustal velocities are interpreted to represent intermediate to mafic granulites. The majority of the earthquake hypocenters as well as a 6.7 km s −1 layer in the velocity model corresponds in depth to thick zones of west dipping midcrustal reflections that may represent major shear zones formed during the late Jurassic Nevadan orogeny or synbatholithic ductile shear zones that accommodated crustal extension associated with batholith intrusion. These reflections are truncated updip by an inferred subvertical contact that coincides with the western edge of the Sierra Nevada batholith and the southward trace of the Bear Mountains fault zone. The updip truncation of midcrustal shear zones and high lower crustal velocities indicate that strike‐slip faulting and magmatic underplating can be important processes during the docking and welding of an accreted terrane.

California↗