Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research-- Solid Earth”

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 343 records · Page 19Linked to original sources
USGS2017

Low stress drops observed for aftershocks of the 2011 Mw 5.7 Prague, Oklahoma, earthquake

In November 2011, three M w ≥ 4.8 earthquakes and thousands of aftershocks occurred along the structurally complex Wilzetta fault system near Prague, Oklahoma. Previous studies suggest that wastewater injection induced a M w 4.8 foreshock, which subsequently triggered a M w 5.7 mainshock. We examine source properties of aftershocks with a standard Brune-type spectral model and jointly solve for seismic moment ( M 0 ), corner frequency ( f 0 ), and kappa ( κ ) with an iterative Gauss-Newton global downhill optimization method. We examine 934 earthquakes with initial moment magnitudes ( M w ) between 0.33 and 4.99 based on the pseudospectral acceleration and recover reasonable M 0 , f 0 , and κ for 87 earthquakes with M w 1.83–3.51 determined by spectral fit. We use M 0 and f 0 to estimate the Brune-type stress drop, assuming a circular fault and shear-wave velocity at the hypocentral depth of the event. Our observations suggest that stress drops range between 0.005 and 4.8 MPa with a median of 0.2 MPa (0.03–26.4 MPa with a median of 1.1 MPa for Madariaga-type), which is significantly lower than typical eastern United States intraplate events (>10 MPa). We find that stress drops correlate weakly with hypocentral depth and magnitude. Additionally, we find the stress drops increase with time after the mainshock, although temporal variation in stress drop is difficult to separate from spatial heterogeneity and changing event locations. The overall low median stress drop suggests that the fault segments may have been primed to fail as a result of high pore fluid pressures, likely related to nearby wastewater injection.

Oklahoma↗
USGS2016

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Quantitative estimations of the processes that control the mass balance of sulfur between melt, crystals, and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere. The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H 2 O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma, and (3) the transport of sulfur-rich vapor (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed-system processes where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “excess sulfur” problem, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by amount of sulfur originally dissolved in erupted melt, as estimated from melt inclusion sulfur concentrations (the “petrologic estimate”), reflects the challenges in closing the sulfur mass balance between crystals, melt, and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed- and open-system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the excess sulfur observed in crystal-rich arc magmas. This result does not negate the important role of vapor migration in sulfur mass balance but rather points out that second boiling (in situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for crystal-rich systems. In contrast, in crystal-poor systems, magma recharge that releases sulfur-rich bubbles is necessary and most likely the primary contributor to sulfur mass balance. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan supereruption in Argentina (2.08 Ma) and show the potential importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California↗
USGS2017

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↗
USGS2017

Sediment gravity flows triggered by remotely generated earthquake waves

Recent great earthquakes and tsunamis around the world have heightened awareness of the inevitability of similar events occurring within the Cascadia Subduction Zone of the Pacific Northwest. We analyzed seafloor temperature, pressure, and seismic signals, and video stills of sediment-enveloped instruments recorded during the 2011–2015 Cascadia Initiative experiment, and seafloor morphology. Our results led us to suggest that thick accretionary prism sediments amplified and extended seismic wave durations from the 11 April 2012 M w 8.6 Indian Ocean earthquake, located more than 13,500 km away. These waves triggered a sequence of small slope failures on the Cascadia margin that led to sediment gravity flows culminating in turbidity currents. Previous studies have related the triggering of sediment-laden gravity flows and turbidite deposition to local earthquakes, but this is the first study in which the originating seismic event is extremely distant (> 10,000 km). The possibility of remotely triggered slope failures that generate sediment-laden gravity flows should be considered in inferences of recurrence intervals of past great Cascadia earthquakes from turbidite sequences. Future similar studies may provide new understanding of submarine slope failures and turbidity currents and the hazards they pose to seafloor infrastructure and tsunami generation in regions both with and without local earthquakes.

Cascadia subduction zone↗
USGS2017

Integrating puffing and explosions in a general scheme for Strombolian-style activity

Strombolian eruptions are among the most common subaerial styles of explosive volcanism worldwide. Distinctive features of each volcano lead to a correspondingly wide range of variations of magnitude and erupted products, but most papers focus on a single type of event at a single volcano. Here, in order to emphasize the common features underlying this diversity of styles, we scrutinize a database from 35 different erupting vents, including 21 thermal infrared videos from Stromboli (Italy), Etna (Italy), Yasur (Vanuatu), and Batu Tara (Indonesia), from puffing, through rapid explosions to normal explosions, with variable ejection parameters and relative abundance of gas, ash, and bombs. Using field observations and high-speed thermal infrared videos processed by a new algorithm, we identify the distinguishing characteristics of each type of activity and how they may relate and interact. In particular, we record that ash-poor normal explosions may be preceded and followed by the onset or the increase of the puffing activity, while ash-rich explosions are emergent, i.e., with inflation of the free surface followed directly by emission of increasingly large gas pockets. Overall, we see that all Strombolian activities form a continuum arising from a common mechanism and are modulated by the combination of two well-established controls: (1) the length of the bursting gas pocket with respect to the vent diameter and (2) the presence and thickness of a high-viscosity layer in the uppermost part of the volcanic conduit.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Temporal variation of tectonic tremor activity in southern Taiwan around the 2010 ML6.4 Jiashian Earthquake

Deep tectonic tremor, which is extremely sensitive to small stress variations, could be used to monitor fault-zone processes during large earthquake cycles and aseismic processes before large earthquakes. In this study, we develop an algorithm for the automatic detection and location of tectonic tremor beneath the southern Central Range of Taiwan and examine the spatio-temporal relationship between tremor and the 4 March 2010 ML6.4 Jiashian earthquake, located about 20 km from active tremor sources. We find that tremor in this region has a relatively short duration, short recurrence time, and no consistent correlation with surface GPS data. We find a short-term increase in the tremor rate 19 days before the Jiashian mainshock, and around the time when the tremor rate began to rise, one GPS station recorded a flip in its direction of motion. We hypothesize that tremor is driven by a slow-slip event that preceded the occurrence of the shallower Jiashian mainshock, even though the inferred slip is too small to be observed by all GPS stations. Our study shows that tectonic tremor may reflect stress variation during the pre-nucleation process of a nearby earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Remote measurement of high preeruptive water vapor emissions at Sabancaya volcano by passive differential optical absorption spectroscopy

Water (H 2 O) is by far the most abundant volcanic volatile species and plays a predominant role in driving volcanic eruptions. However, numerous difficulties associated with making accurate measurements of water vapor in volcanic plumes have limited their use as a diagnostic tool. Here we present the first detection of water vapor in a volcanic plume using passive visible-light differential optical absorption spectroscopy (DOAS). Ultraviolet and visible-light DOAS measurements were made on 21 May 2016 at Sabancaya Volcano, Peru. We find that Sabancaya's plume contained an exceptionally high relative water vapor abundance 6 months prior to its November 2016 eruption. Our measurements yielded average sulfur dioxide (SO 2 ) emission rates of 800–900 t/d, H 2 O emission rates of around 250,000 t/d, and an H 2 O/SO 2 molecular ratio of 1000 which is about an order of magnitude larger than typically found in high-temperature volcanic gases. We attribute the high water vapor emissions to a boiling-off of Sabancaya's hydrothermal system caused by intrusion of magma to shallow depths. This hypothesis is supported by a significant increase in the thermal output of the volcanic edifice detected in infrared satellite imagery leading up to and after our measurements. Though the measurement conditions encountered at Sabancaya were very favorable for our experiment, we show that visible-light DOAS systems could be used to measure water vapor emissions at numerous other high-elevation volcanoes. Such measurements would provide observatories with additional information particularly useful for forecasting eruptions at volcanoes harboring significant hydrothermal systems.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Rock friction under variable normal stress

This study is to determine the detailed response of shear strength and other fault properties to changes in normal stress at room temperature using dry initially bare rock surfaces of granite at normal stresses between 5 and 7 MPa. Rapid normal stress changes result in gradual, approximately exponential changes in shear resistance with fault slip. The characteristic length of the exponential change is similar for both increases and decreases in normal stress. In contrast, changes in fault normal displacement and the amplitude of small high-frequency elastic waves transmitted across the surface follow a two stage response consisting of a large immediate and a smaller gradual response with slip. The characteristic slip distance of the small gradual response is significantly smaller than that of shear resistance. The stability of sliding in response to large step decreases in normal stress is well predicted using the shear resistance slip length observed in step increases. Analysis of the shear resistance and slip-time histories suggest nearly immediate changes in strength occur in response to rapid changes in normal stress; these are manifested as an immediate change in slip speed. These changes in slip speed can be qualitatively accounted for using a rate-independent strength model. Collectively, the observations and model show that acceleration or deceleration in response to normal stress change depends on the size of the change, the frictional characteristics of the fault surface, and the elastic properties of the loading system.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Connecting crustal seismicity and earthquake-driven stress evolution in Southern California

Tectonic stress in the crust evolves during a seismic cycle, with slow stress accumulation over interseismic periods, episodic stress steps at the time of earthquakes, and transient stress readjustment during a postseismic period that may last months to years. Static stress transfer to surrounding faults has been well documented to alter regional seismicity rates over both short and long time scales. While static stress transfer is instantaneous and long lived, postseismic stress transfer driven by viscoelastic relaxation of the ductile lower crust and mantle leads to additional, slowly varying stress perturbations. Both processes may be tested by comparing a decade-long record of regional seismicity to predicted time-dependent seismicity rates based on a stress evolution model that includes viscoelastic stress transfer. Here we explore crustal stress evolution arising from the seismic cycle in Southern California from 1981 to 2014 using five M ≥6.5 source quakes: the M 7.3 1992 Landers, M6.5 1992 Big Bear, M 6.7 1994 Big Bear, M 7.1 1999 Hector Mine, and M 7.2 2010 El Mayor-Cucapah earthquakes. We relate the stress readjustment in the surrounding crust generated by each quake to regional seismicity using rate-and-state friction theory. Using a log likelihood approach, we quantify the potential to trigger seismicity of both static and viscoelastic stress transfer, finding that both processes have systematically shaped the spatial pattern of Southern California seismicity since 1992.

California↗
USGS2017

Delayed seismicity rate changes controlled by static stress transfer

On 15 June 2010, a M w 5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity.

California↗
USGS2017

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

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

California↗
USGS2017

Refining the formation and early evolution of the Eastern North American Margin: New insights from multiscale magnetic anomaly analyses

To investigate the oceanic lithosphere formation and early seafloor spreading history of the North Atlantic Ocean, we examine multiscale magnetic anomaly data from the Jurassic/Early Cretaceous age Eastern North American Margin (ENAM) between 31 and 40°N. We integrate newly acquired sea surface magnetic anomaly and seismic reflection data with publicly available aeromagnetic and composite magnetic anomaly grids, satellite-derived gravity anomaly, and satellite-derived and shipboard bathymetry data. We evaluate these data sets to (1) refine magnetic anomaly correlations throughout the ENAM and assign updated ages and chron numbers to M0–M25 and eight pre-M25 anomalies; (2) identify five correlatable magnetic anomalies between the East Coast Magnetic Anomaly (ECMA) and Blake Spur Magnetic Anomaly (BSMA), which may document the earliest Atlantic seafloor spreading or synrift magmatism; (3) suggest preexisting margin structure and rifting segmentation may have influenced the seafloor spreading regimes in the Atlantic Jurassic Quiet Zone (JQZ); (4) suggest that, if the BSMA source is oceanic crust, the BSMA may be M series magnetic anomaly M42 (~168.5 Ma); (5) examine the along and across margin variation in seafloor spreading rates and spreading center orientations from the BSMA to M25, suggesting asymmetric crustal accretion accommodated the straightening of the ridge from the bend in the ECMA to the more linear M25; and (6) observe anomalously high-amplitude magnetic anomalies near the Hudson Fan, which may be related to a short-lived propagating rift segment that could have helped accommodate the crustal alignment during the early Atlantic opening.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

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↗
USGS2017

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

A decade of induced slip on the causative fault of the 2015 Mw 4.0 Venus earthquake, northeast Johnson County, Texas

On 7 May 2015, a M w 4.0 earthquake occurred near Venus, northeast Johnson County, Texas, in an area of the Bend Arch-Fort Worth Basin that reports long-term, high-volume wastewater disposal and that has hosted felt earthquakes since 2009. In the weeks following the M w 4.0 earthquake, we deployed a local seismic network and purchased nearby active-source seismic reflection data to capture additional events, characterize the causative fault, and explore potential links between ongoing industry activity and seismicity. Hypocenter relocations of the resulting local earthquake catalog span ~4–6 km depth and indicate a fault striking ~230 ° , dipping to the west, consistent with a nodal plane of the M w 4.0 regional moment tensor. Fault plane solutions indicate normal faulting, with B axes striking parallel to maximum horizontal compressive stress. Seismic reflection data image the reactivated basement fault penetrating the Ordovician disposal layer and Mississippian production layer, but not displacing post-Lower Pennsylvanian units. Template matching at regional seismic stations indicates that low-magnitude earthquakes with similar waveforms began in April 2008, with increasing magnitude over time. Pressure data from five saltwater disposal wells within 5 km of the active fault indicate a disposal formation that is 0.9–4.8 MPa above hydrostatic. We suggest that the injection of 28,000,000 m 3 of wastewater between 2006 and 2015 at these wells led to an increase in subsurface pore fluid pressure that contributed to inducing this long-lived earthquake sequence. The 2015 M w 4.0 event represents the largest event in the continuing evolution of slip on the causative fault.

Texas↗
USGS2017

Magnetotelluric imaging of lower crustal melt and lithospheric hydration in the Rocky Mountain Front transition zone, Colorado, USA

We present an electrical resistivity model of the crust and upper mantle from two‐dimensional (2‐D) anisotropic inversion of magnetotelluric data collected along a 450 km transect of the Rio Grande rift, southern Rocky Mountains, and High Plains in Colorado, USA. Our model provides a window into the modern‐day lithosphere beneath the Rocky Mountain Front to depths in excess of 150 km. Two key features of the 2‐D resistivity model are (1) a broad zone (~200 km wide) of enhanced electrical conductivity (<20 Ωm) in the midcrust to lower crust that is centered beneath the highest elevations of the southern Rocky Mountains and (2) hydrated lithospheric mantle beneath the Great Plains with water content in excess of 100 ppm. We interpret the high conductivity region of the lower crust as a zone of partially molten basalt and associated deep‐crustal fluids that is the result of recent (less than 10 Ma) tectonic activity in the region. The recent supply of volatiles and/or heat to the base of the crust in the late Cenozoic implies that modern‐day tectonic activity in the western United States extends to at least the western margin of the Great Plains. The transition from conductive to resistive upper mantle is caused by a gradient in lithospheric modification, likely including hydration of nominally anhydrous minerals, with maximum hydration occurring beneath the Rocky Mountain Front. This lithospheric “hydration front” has implications for the tectonic evolution of the continental interior and the mechanisms by which water infiltrates the lithosphere.

Colorado↗