Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Atmospheres”

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 973 records · Page 54Linked to original sources
USGS2020

Hydrologically induced deformation in Long Valley Caldera and adjacent Sierra Nevada

Vertical and horizontal components of GNSS displacements in the Long Valley Caldera and adjacent Sierra Nevada range show a clear correlation with hydrological trends at both multiyear and seasonal time scales. We observe a clear vertical and horizontal seasonal deformation pattern primarily attributable to the solid earth response to hydrological surface loading at large-to-regional (Sierra Nevada range) scales. Several GNSS sites, located at the eastern edge of the Sierra Nevada along the southwestern rim of Long Valley Caldera, also show significant horizontal deformation that cannot be explained by elastic deformation from surface loading. Due to the location of these sites and the strong correlation between their horizontal displacements and spring discharge, we hypothesize that this deformation reflects poroelastic processes related to snowmelt runoff water infiltrating into the Sierra Nevada slopes around Long Valley Caldera. Interestingly, this is also an area where water infiltrates to feed the local hydrothermal system, and where snowmelt-induced earthquake swarms have been recently detected.

California, Nevada↗
USGS2020

The role of pre-magmatic rifting in shaping a volcanic continental margin: An example from the Eastern North American Margin

Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2‐D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along‐strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post‐rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn‐rift structures, and the SDR there extend ~ 200 km seaward of the hinge line. Early post‐rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn‐rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.

Journal of Geophysical Research-- Solid Earth↗
USGS2020

Shallow basin structure and attenuation are key to predicting long shaking duration in Los Angeles Basin

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three-dimensional simulations using the CVM-S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM-H do not match the timing of the initial arrivals as well as CVM-S4.26.M01, and the strong late arrivals in the CVM-H simulation travel with an apparent velocity slower than observed. A higher-quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off-azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long-period ground motions in basins.

California↗
USGS2020

Origin and properties of hydrothermal tremor at Lone Star Geyser, Yellowstone National Park, USA

Geysers are rare geologic features that intermittently discharge liquid water and steam driven by heating and decompression boiling. The cause of variability in eruptive styles and the associated seismic signals are not well understood. Data collected from five broadband seismometers at Lone Star Geyser, Yellowstone National Park are used to determine the properties, location, and temporal patterns of hydrothermal tremor. The tremor is harmonic at some stages of the eruption cycle and is caused by near‐periodic repetition of discrete seismic events. Using the polarization of ground motion, we identify the location of tremor sources throughout several eruption cycles. During preplay episodes (smaller eruptions preceding the more vigorous major eruption), tremor occurs at depths of 7–10 m and is laterally offset from the geyser's cone by ~5 m. At the onset of the main eruption, tremor sources migrate laterally and become shallower. As the eruption progresses, tremor sources migrate along the same path but in the opposite direction, ending where preplay tremor originates. The upward and then downward migration of tremor sources during eruptions are consistent with warming of the conduit followed by evacuation of water during the main eruption. We identify systematic relations among the two types of preplays, discharge, and the main eruption. A point‐source moment tensor fit to low‐frequency waveforms of an individual tremor event using half‐space velocity models indicates average V S ≳ 0.8 km/s, source depths ~4–20 m, and moment tensors with primarily positive isotropic and compensated linear vector dipole moments.

Wyoming↗
USGS2020

Along-margin variations in breakup volcanism at the Eastern North American Margin

We model the magnetic signature of rift-related volcanism to understand the distribution and volumeofmagmatic activity that occurred during the breakup of Pangaea and early Atlantic opening at the Eastern North American Margin (ENAM).Along-strike variations in the amplitude and character of the prominent East Coast Magnetic Anomaly (ECMA) suggest that the emplacement of the volcanic layers producing this anomaly similarly varied along the margin. We use three-dimensional magnetic forward modeling constrained by seismic interpretationsto identify along-margin variations in volcanic thickness and width that can explain the observed amplitude and character of the ECMA. Our model results suggest that the ECMA is produced by a combination of both first-order (~600-1000 km)and second-order (~50-31100 km) magmatic segmentation. The first-order magmatic segmentation could have resulted from preexisting variations in crustal thickness and rheology developed during the tectonic amalgamation of Pangaea. The second-order magmatic segmentation developed during continental breakup and likely influenced the segmentation and transform fault spacing of the initial, and modern, Mid-Atlantic Ridge. These variations in magmatism showhow extension and thermal weakening was distributed at the ENAM during continental breakup and how this breakup magmatism was related to both previous and subsequent Wilson Cycle stages.

East Coast of United States, Atlantic Ocean↗
USGS2020

The induced Mw 5.0 March 2020 west Texas seismic sequence

On March 26, 2020, a M 5.0 earthquake occurred in the Delaware Basin, Texas, near the border between Reeves and Culberson Counties. This was the third largest earthquake recorded in Texas and the largest earthquake in the Central and Eastern United States since the three M 5.0–5.8 induced events in Oklahoma during 2016. Using multistation waveform template matching, we detect 3,940 earthquakes in the sequence with the first event in the area occurring in May 2018. The M 5.0 earthquake sequence occurred on a ENE (∼082°) normal fault dipping ∼37° toward the south. The earthquake caused 6 mm of oblique surface deformation, and geodetic slip inversion suggests slip was isolated above 6 km depth. We find that the sequence was most likely induced by nearby wastewater disposal operations, and seismicity rates in the region surrounding the M 5.0 will likely continue to increase in the future if disposal operations continue unaltered.

Texas↗
USGS2021

Airborne geophysical imaging of weak zones on Iliamna Volcano, Alaska: Implications for slope stability

Water‐saturated, hydrothermally altered rocks reduce the strength of volcanic edifices and increase the potential for sector collapses and far‐traveled mass flows of unconsolidated debris. Iliamna Volcano is an andesitic stratovolcano located on the western side of the Cook Inlet, ∼225 km southwest of Anchorage and is a source of repeated avalanches. The widespread snow and ice cover on Iliamna Volcano make surface alteration difficult to identify. However, intense hydrothermal alteration significantly reduces both the electrical resistivity and magnetization of volcanic rock and can therefore be identified with airborne geophysical measurements. We use airborne electromagnetic and magnetic data to map snow and ice thickness and identify underlying alteration zones at Iliamna Volcano, Alaska. Resistivities were calculated to an average depth of >300 m, and a 3‐D susceptibility model extends from the surface to the base of the volcano, about 3,000 m below the summit. Geophysical models image low resistivity (<30 ohm‐m) and low susceptibilities near the summit of Iliamna and below its older vent complex, with the low susceptibilities indicating alteration up to ∼800 m in thickness. Thin conductors (∼50–100 m thick) on the edifice slopes coincide with recorded locations of repeated debris avalanches over the past ∼60 years and are attributed to saturated zones at high elevation. Three‐dimensional slope stability models based upon the geophysically constrained alteration distribution suggest the edifice of Iliamna is unstable and could lead to collapse scars ∼400 m deep near the current and former vent complexes.

Alaska↗
USGS2021

Characteristics of frequent dynamic triggering of microearthquakes in Southern California

Dynamic triggering of earthquakes has been reported at various fault systems. The triggered earthquakes are thought to be caused either directly by dynamic stress changes due to the passing seismic waves, or indirectly by other nonlinear processes that are initiated by the passing waves. Distinguishing these physical mechanisms is difficult because of the general lack of high‐resolution earthquake catalogs and robust means to quantitatively evaluate triggering responses, particularly, delayed responses. Here we use the high‐resolution Quake Template Matching catalog in Southern California to systematically evaluate teleseismic dynamic triggering patterns in the San Jacinto Fault Zone and the Salton Sea Geothermal Field from 2008 to 2017. We develop a new statistical approach to identify triggered cases, finding that approximately 1 out of every 5 global M w ≥ 6 earthquakes dynamically trigger microearthquakes in Southern California. The triggering responses include both instantaneous and delayed triggering, showing a highly heterogeneous pattern and indicating possible evolving triggering thresholds. We do not observe a clear peak ground velocity triggering threshold that can differentiate triggering earthquakes from non‐triggering events, but there are subtle differences in the frequency content that may possibly differentiate the earthquakes. In contrast to the depth distribution of background seismicity, the identified triggered earthquakes tend to concentrate at the edges of the seismogenic zones. Although instantaneously triggered earthquakes are likely a result of dynamic Coulomb stress changes, the cases of delayed dynamic triggering are best explained by nonlinear triggering processes, including cyclic material fatigue, accelerated transient creep, and stochastic frictional heterogeneities.

California↗
USGS2021

Comparative rhyolite systems: Inferences from vent patterns and eruptive episodicities: Eastern California and Laguna del Maule

Distilling my experience in having field mapped in detail the volcanic fields at Laguna del Maule and Long Valley and having worked out their time-volume-composition magmatic histories, I compare and contrast the postglacial rhyolites of the former with six multi-vent eruptive sequences of rhyolite in California. Compilations and discussions are made of volcanic-field areas and longevities, their compositions, vent distributions, individual batch and total volumes, eruptive episodicities, and tectonic influences. Growth of long-lived pluton-scale reservoirs of granitic crystal mush, from which the rhyolite melts separated, are interpreted in terms of conceptual models I published previously—(1) fundamentally basaltic transcrustal magmatism, 1981; (2) the deep-crustal MASH zone model, 1988; and (3) the rhyolite-melt crystal-mush model, 2001. Inferences and speculations are advanced concerning processes and timescales of rhyolite-melt separation from granitic mush and of prompt or long-delayed subsequent eruption.

California↗
USGS2021

The birth of a Hawaiian fissure eruption

Most basaltic explosive eruptions intensify abruptly, allowing little time to document processes at the start of eruption. One opportunity came with the initiation of activity from fissure 8 (F8) during the 2018 eruption on the lower East Rift Zone of Kīlauea, Hawaii. F8 erupted in four episodes. We recorded 28 min of high‐definition video during a 51‐min period, capturing the onset of the second episode on 5 May. From the videos, we were able to analyze the following in‐flight parameters: frequency and duration of explosions; ejecta heights; pyroclast exit velocities; in‐flight total mass and estimated mass eruption rates; and the in‐flight total grain size distributions. The videos record a transition from initial pulsating outgassing, via spaced, but increasingly rapid, discrete explosions, to quasisustained, unsteady fountaining. This transition accompanied waxing intensity (mass flux) of the F8 eruption. We infer that all activity was driven by a combination of the ascent of a coupled mixture of small bubbles and melt, and the buoyant rise of decoupled gas slugs and/or pockets. The balance between these two types of concurrent flow determined the exact form of the eruptive activity at any point in time, and changes to their relative contributions drove the transition we observed at early F8. Qualitative observations of other Hawaiian fountains at Kīlauea suggest that this physical model may apply more generally. This study demonstrates the value of in‐flight parameters derived from high‐resolution videos, which offer a rapid and highly time‐sensitive alternative to measurements based on sampling of deposits posteruption.

Island of Hawai'i↗
USGS2021

The seismo-acoustics of submarine volcanic eruptions

Many of the world’s volcanoes are hidden beneath the ocean’s surface where eruptions are difficult to observe. However, seismo‐acoustic signals produced by these eruptions provide a useful means of identifying active submarine volcanism. A literature survey revealed reports of 119 seismo‐acoustically recorded submarine eruptions since 1939. Submarine eruptions have been recorded in all major tectonic settings, with a range of geochemistries, and at a variety of water depths, but the reports are dominated by eruptions in the Pacific and at only a few locations. Many of the reports offer little detail, with over half of the observations made from distances >500 km, and only about half were confirmed as eruptions by non‐seismo‐acoustic evidence. The reported seismo‐acoustic signals cover a wide variety of processes, including earthquakes, explosions, various types of tremor, signals related to lava extrusion, and landslides. Recorded signals can sometimes be difficult to classify or confidently associate with an eruption, although there has been progress in this regard. Real‐time monitoring of submarine eruptions has been on‐going for several decades on regional and global scales with growing interest and effort in local networks. Real‐time networks are complemented by short‐term instrument deployments that often give more detailed insights into the dynamics and processes of submarine eruptions. Thorough seismo‐acoustic monitoring and study has increased our understanding of submarine eruptions, especially of deep‐sea volcanoes and spreading centers. Despite this, there are still many outstanding questions that need to be addressed for submarine volcanoes to be as well understood and monitored as their terrestrial counterparts.

Journal of Geophysical Research↗
USGS2021

Limited mantle hydration by bending faults at the Middle America Trench

Seismic anisotropy measurements show that upper mantle hydration at the Middle America Trench (MAT) is limited to serpentinization and/or water in fault zones, rather than distributed uniformly. Subduction of hydrated oceanic lithosphere recycles water back into the deep mantle, drives arc volcanism, and affects seismicity at subduction zones. Constraining the extent of upper mantle hydration is an important part of understanding many fundamental processes on Earth. Substantially reduced seismic velocities in tomography suggest that outer rise plate‐bending faults provide a pathway for seawater to rehydrate the slab mantle just prior to subduction. Estimates of outer‐rise hydration based on tomograms vary significantly, with some large enough to imply that, globally, subduction has consumed more than two oceans worth of water during the Phanerozoic. We found that, while the mean upper mantle wavespeed is reduced at the MAT outer rise, the amplitude and orientation of inherited anisotropy are preserved at depths >1 km below the Moho. At shallower depths, relict anisotropy is replaced by slowing in the fault‐normal direction. These observations are incompatible with pervasive hydration but consistent with models of wave propagation through serpentinized fault zones that thin to <100‐m in width at depths >1 km below Moho. Confining hydration to fault zones reduces water storage estimates for the MAT upper mantle from ∼3.5 wt% to <0.9 wt% H 2 0. Since the intermediate thermal structure in the ∼24 Myr‐old MAT slab favors serpentinization, limited hydration suggests that fault mechanics are the limiting factor, not temperatures. Subducting mantle may be similarly dry globally.

Journal of Geophysical Research↗
USGS2021

B-positive: A robust estimator of aftershock magnitude distribution in transiently incomplete catalogs

The earthquake magnitude-frequency distribution is characterized by the b -value, which describes the relative frequency of large versus small earthquakes. It has been suggested that changes in b -value after an earthquake can be used to discriminate whether that earthquake is part of a foreshock sequence or a more typical mainshock-aftershock sequence, with a decrease in b -value heralding a larger earthquake to come. However, the measurement of b -value during an active aftershock sequence is strongly biased by short-term incompleteness of the earthquake catalog and by data-windowing, and these biases have the same direction as the proposed signal. Here I develop a new estimator of the b- value that is insensitive to transient changes in catalog completeness and that does not require data windowing. The new estimator “ b -positive” is based on the positive-only subset of the differences in magnitude between successive earthquakes, which are described by a double-exponential (Laplace) distribution with the same b -value as the magnitude distribution itself. The b -positive estimator greatly improves the robustness of continuous b -value measurements during active earthquake sequences, as well as in historical catalogs with unknown or variable completeness. The new estimator confirms some of the observations of Gulia and Wiemer (2019), although at a reduced level, showing a decrease and recovery of the b -value during several recent foreshock sequences that cannot be attributed simply to measurement bias. However, the unbiased b -value changes may be too subtle to use in a real-time earthquake alarm system.

Journal of Geophysical Research: Solid Earth↗
USGS2021

Heat flux from a vapor-dominated hydrothermal field beneath Yellowstone Lake

We report results from 149 heat flux measurements made over n ∼2-year interval at sites in and around a vapor-dominated geothermal field located at water depths of ∼100–120 m in Yellowstone Lake, Wyoming. Measurements of both in situ temperature and thermal conductivity as a function of depth were made with a 1 m probe via a remotely operated vehicle, and are combined to compute the vertical conductive heat flux. Inside the ∼55.5 × 10 3 m 2 bathymetric depression demarcating the vapor-dominated field, the median conductive flux is 13 W m −2 , with a conductive output of 0.72 MW. Outside the thermal field, the median conductive flux is 3.5 W m −2 . We observed 49 active vents inside the thermal field, with an estimated mass discharge rate of 56 kg s −1 , a median exit-fluid temperature of 132°C, and a total heat output of 29 MW. We find evidence for relatively weak secondary convection with a total output of 0.09 MW in thermal area lake floor sediments. Our data indicate that vapor beneath the thermal field is trapped by a low-permeability cap at a temperature of ∼189°C and a depth of ∼15 m below the lake floor. The thermal output of the Deep Hole is among the highest of any vapor-dominated field in Yellowstone, due in part to the high boiling temperatures associated with the elevated lake floor pressures.

Wyoming↗
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↗
USGS2021

Geometry of the décollement below eastern Bangladesh and implications for seismic hazard

Eastern Bangladesh sits on the seismically active Chittagong-Myanmar fold and thrust belt (CMFB), a north-trending accretionary wedge on the eastern side of the India-Eurasia collision. Earthquakes on the basal décollement and associated thrusts within the CMFB present a hazard to this densely populated region. In this study, we interpret 28 seismic reflection profiles from both published and unpublished sources to constrain the depth of the basal décollement. To convert profiles from the time domain to the depth domain, we integrate sonic log and seismic stacking velocity data to generate time-velocity relationships for different parts of the CMFB. Our analysis reveals that the décollement is ∼9 km deep in northeast and southeast Bangladesh, but shallows to ∼5 km in east-central Bangladesh. The décollement has an area of 7.25 × 10 4 km 2 (∼150 × 450 km), making it capable of an Mw 8.5 earthquake. However, the warped geometry of this fault might act as a rupture barrier were a large earthquake to occur on the décollement. Our combined velocity and fault model lay the groundwork for future studies to address seismic segmentation, ground shaking, and rupture modeling in the CMFB. Finally, we use our compiled data set to analyze the evolution of fold kinematics in the CMFB. We observe that folding style and failure mode varies, from mainly ductile deformation in the foreland to mainly brittle in the hinterland. The dual-failure modes within the CMFB support the hypothesis that a region with ductile deformation may still be capable of seismic behavior.

Journal of Geophysical Research: Solid Earth↗
USGS2021

Sources of volcanic tremor associated with the summit caldera collapse during the 2018 east rift eruption of Kīlauea Volcano, Hawai'i

Volcanic tremor occurring at the beginning of the 2018 Kīlauea eruption is characterized using both seismic and tilt data recorded at the Kīlauea summit. An automatic seismic network-based approach detects several types of tremor including (a) 0.5–1 Hz long-period tremor preceding the eruption, located at the south-southwest edge of Halema'uma'u Crater and attributed to the quasi-steady radiation from a shallow hydrothermal system and (b) two sequences of gliding tremor at the beginning of the eruption, both with locations on the edges of the crater and within it. The first sequence is attributed to two swarms of low-amplitude regularly repeating earthquakes induced by the jerky motions of a cylindrical rock piston with radius of 325 m, height of 250 m, and mass of 2.07 × 10 11 kg, progressively intruding 12.3 m into the shallow hydrothermal system with volume of 10 8 m 3 and depth extent of 300 m. The second sequence is attributed to a gradual evolution in the properties of a bubbly magma within an east-striking dike below Halema'uma'u Crater, impacted by repeated roof collapses. A fluid-filled crack model points to a decrease in gas volume fraction from 4.22% to 1.6 × 10 −2 % in the magma filling the dike, and a model of gas retro-diffusion within the melt suggests a two orders of magnitude decrease in bubble number density from 7 × 10 8 m −3 down to 4 × 10 6 m −3 . Both models feature a quasi to totally degassed magma by May 26.

Hawaii↗
USGS2021

Earthquake source mechanisms and stress field variations associated with wastewater-induced seismicity in southern Kansas, USA

The strong increase of seismicity rates in the contiguous USA over the last 10 years is linked to the injection of huge amounts of wastewater from oil and gas production in unconventional hydrocarbon reservoirs. We calculated 549 moment tensors of induced earthquakes ( M W ≤ 4.9) in southern Kansas to study their source mechanisms and their relation to injection activity. Seventeen percent of the events analyzed contained significant volumetric (ISO%) components, and these events mostly occurred near the two largest local earthquakes during the 4 months of largest active wastewater disposal. Mapping the local stress field, we determined that most of the region lies within a transtensional stress regime, with a maximum horizontal stress σ Hmax trending N75°E. In the epicentral area of the M W 4.9 Milan earthquake, the σ Hmax trend is rotated to about S80°E. Locally, two areas display a change in the stress field orientation with depth, from transtensional above 5.5 km depth to strike slip deeper in the basement. Relating the resolved fault geometries to the obtained local stress field orientation, we find that most of the activated fault planes were optimally oriented to the current stress field and thus small stress perturbations caused by the water injection could lead to failure.

Kansas↗