Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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 847 records · Page 47Linked to original sources

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Using tracer variance decay to quantify variability of salinity mixing in the Hudson River Estuary

The salinity structure in an estuary is controlled by time‐dependent mixing processes. However, the locations and temporal variability of where significant mixing occurs is not well‐understood. Here we utilize a tracer variance approach to demonstrate the spatial and temporal structure of salinity mixing in the Hudson River Estuary. We run a 4‐month hydrodynamic simulation of the tides, currents, and salinity that captures the spring‐neap tidal variability as well as wind‐driven and freshwater flow events. On a spring‐neap time scale, salinity variance dissipation (mixing) occurs predominantly during the transition from neap to spring tides. On a tidal time scale, 60% of the salinity variance dissipation occurs during ebb tides and 40% during flood tides. Spatially, mixing during ebbs occurs primarily where lateral bottom salinity fronts intersect the bed at the transition from the main channel to adjacent shoals. During ebbs, these lateral fronts form seaward of constrictions located at multiple locations along the estuary. During floods, mixing is generated by a shear layer elevated in the water column at the top of the mixed bottom boundary layer, where variations in the along channel density gradients locally enhance the baroclinic pressure gradient leading to stronger vertical shear and more mixing. For both ebb and flood, the mixing occurs at the location of overlap of strong vertical stratification and eddy diffusivity, not at the maximum of either of those quantities. This understanding lends a new insight to the spatial and time dependence of the estuarine salinity structure.

Hudson River

Seasonal controls on sediment delivery and hydrodynamics in a vegetated tidally influenced interdistributary island

River deltas are maintained by a continuous supply of terrestrial sediments that provide critical land building material to help sustain and protect vulnerable ecological communities and serve as natural storm protection barriers. Local hydrodynamics are important in determining the degree to which fluvial sediments are removed from the water column and retained on the delta complex. During 2014, we measured hydrodynamics and sediment transport characteristics at one of the world's most rapidly prograding deltas, the Wax Lake delta in Louisiana, USA. We observed waves to be the dominant source of bottom stress for 70% of our observations. Sediment concentration tended to increase with shear stress, but only after stresses exceeded 0.01–0.02 Pa. Significant wave height and bottom stress were substantially reduced after June, when the emergence of American lotus ( Nelumbo lutea ) formed a dense canopy over the intertidal regions of the island splay. Hydrodynamics during these summer vegetated conditions were much more favorable to floc formation, and by extension particle settling, as shown by trends in the Kolmogorov microscale parameter over the course of the measurement campaign. Together, these findings suggest that the timing between peak river discharge and the emergence of vegetation may have a strong influence on rates of progradation in seasonally vegetated delta splays, whereby sediments delivered by flood events that extend late into summer may be governed by hydrodynamics that favor particle deposition, whereas those delivered prior to the summer may be more prone to remain in suspension and bypass the delta complex.

Louisiana

Persistent multidecadal variability since the 15th century in the southern Barents Sea derived from annually resolved shell-based records

In the North Atlantic Ocean, multidecadal variability in sea surface temperatures (SSTs) over the past several centuries has largely been inferred through terrestrial proxies and decadally resolved marine proxies. Annually resolved proxy records from marine archives provide valuable insight into this variability, but are especially rare from high latitude environments, particularly for centennial timescales. We constructed continuous, absolutely dated records of shell growth (1449–2014 CE; 564 years) and oxygen isotope ratios (δ 18 O shell ; 1539–2014 CE; 476 years) from shells of the bivalve Arctica islandica from coastal northern Norway, a location sensitive to large-scale North Atlantic Ocean dynamics. An annual (January–December) SST reconstruction derived from δ 18 O shell for the past five centuries suggests an increase of at least 2°C from the mid-18th century to 2014. The SST reconstruction correlates significantly with instrumental records and with other proxy reconstructions in the southern Barents Sea region. Spectral analysis of the shell growth and isotope records supports evidence for Atlantic multidecadal variability (65–80 year periodicity) extending into polar and subpolar latitudes for the past five centuries. These results provide additional evidence that multidecadal variability in SSTs are a persistent feature of the North Atlantic marine system.

Greenland Sea, Iceland Sea, Norwegian Sea, Barents

Compositional layering in Io driven by magmatic segregation and volcanism

The compositional evolution of volcanic bodies like Io is not well understood. Magmatic segregation and volcanic eruptions transport tidal heat from Io's interior to its surface. Several observed eruptions appear to be extremely high temperature (≥ 1600 K), suggesting either very high degrees of melting, refractory source regions, or intensive viscous heating on ascent. To address this ambiguity, we develop a model that couples crust and mantle dynamics to a simple compositional system. We analyse the model to investigate chemical structure and evolution. We demonstrate that magmatic segregation and volcanic eruptions lead to stratification of the mantle, the extent of which depends on how easily high temperature melts from the more refractory lower mantle can migrate upwards. We propose that Io's highest temperature eruptions originate from this lower mantle region, and that such eruptions act to limit the degree of compositional stratification.

Journal of Geophysical Research