Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research B: 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 721 records · Page 40Linked to original sources
USGS2012

Laboratory studies of volcanic jets

The study of the fluid dynamics of violent volcanic eruptions by laboratory experiment is described, and the important fluid-dynamic processes that can be examined in laboratory models are discussed in detail. In preliminary experiments, pure gases are erupted from small reservoirs. The gases used are Freon 12 and Freon 22, two gases of high molecular weight and high density that are good analogs of heavy and particulate-laden volcanic gases; nitrogen, a moderate molecular weight, moderate density gas for which the thermodynamic properties are well known; and helium, a low molecular weight, lowdensity gas that is used as a basis for comparison with the behavior of the heavier gases and as an analog of steam, the gas that dominates many volcanic eruptions. Transient jets erupt from the reservoir into the laboratory upon rupture of a thin diaphragm at the exit of a convergent nozzle. The gas accelerates from rest in the reservoir to high velocity in the jet. Reservoir pressures and geometries are such that the fluid velocity in the jets is initially supersonic and later decays to subsonic. The measured reservoir pressure decreases as the fluid expands through repetitively reflecting rarefaction waves, but for the conditions of these experiments, a simple steady-discharge model is sufficient to explain the pressure decay and to predict the duration of the flow. Density variations in the flow field have been visualized with schlieren and shadowgraph photography. The observed structure of the jet is correlated with the measured pressure history. The starting vortex generated when the diaphragm ruptures becomes the head of the jet. Though the exit velocity is sonic, the flow head in the helium jet decelerates to about one-third of sonic velocity in the first few nozzle diameters, the nitrogen head decelerates to about three-fourths of sonic velocity, while Freon maintains nearly sonic velocity. The impulsive acceleration of reservoir fluid into the surrounding atmosphere produces a compression wave. The strength of this wave depends primarily on the sound speed of the fluid in the reservoir but also, secondarily with opposite effect, on the density: helium produces a relatively strong atmospheric shock while the Freons do not produce any optically observable wave front. Well-formed N waves are detected with a microphone far from the reservoir. Barrel shocks, Mach disks, and other familiar features of steady underexpanded supersonic jets form inside the jet almost immediately after passage of the flow head. These features are maintained until the pressure in the reservoir decays to sonic conditions. At low pressures the jets are relatively structureless. Gas-particle jets from volcanic eruptions may behave as pseudogases if particle concentrations and mass and momentum exchange between the components are sufficiently small. The sound speed of volcanic pseudogases can be as large as 1000 m s −1 or as small as a few tens of meters per second depending on the mass loading and initial temperature. Fluids of high sound speed produce stronger atmospheric shock waves than do those of low sound speed. Therefore eruption of a hot gas lightly laden with particulates should produce a stronger shock than eruption of a cooler or heavily laden fluid. An empirical expression suggests that the initial velocity of the head of supersonic volcanic jets is controlled by the sound speed and the ratio of the density of the erupting fluid to that of the atmosphere. The duration of gas or pseudogas eruptions is controlled by the sound speed of the fluid and the ratio of reservoir volume to vent area.

Journal of Geophysical Research Solid Earth↗
USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗
USGS2012

The growth of geological structures by repeated earthquakes: 2, Field examples of continental dip-slip faults

A strong test of our understanding of the earthquake cycle is the ability to reproduce extant fault-bounded geological structures, such as basins and ranges, which are built by repeated cycles of deformation. Along strike-slip faults, the coseismic and interseismic deformation can be nearly equal in magnitude and opposite in sign, resulting in little permanent deformation except for the fault offset. For dip-slip faults, portions of the crust are lifted and dropped, and so buoyancy forces are exerted. The seismic and interseismic deformations do not balance, and structures grow and become subject to erosion and deposition. We consider three examples for which the structure and fault geometry are well known: the White Wolf reverse fault in California, site of the 1952 Kern County M=7.3 earthquake, the Lost River normal fault in Idaho, site of the 1983 Borah Peak M=7.0 earthquake, and the Cricket Mountain normal fault in Utah, site of Quaternary slip events. Basin stratigraphy and seismic reflection records are used to profile the structure, and coseismic deformation measured by leveling surveys is used to estimate the fault geometry. To reproduce these structures, we add the deformation associated with the earthquake cycle (the coseismic slip and postseismic relaxation) to the flexure caused by the observed sediment load, treating the crust as a thin elastic plate overlying a fluid substrate. The cumulative deformation is principally dependent on the elastic plate thickness, modestly sensitive to the sediment-substrate density difference, and insensitive to the fluid viscosity for the 4- to 8-Ma structures. We deduce a longterm flexural rigidity of 2–15 × 10 19 Nm; this is equivalent to an elastic plate thickness of 2–4 km for a Young's modulus of 2.5 × 10 10 Nm −2 . This value is found where independent estimates of the elastic thickness from the coherence between surface topography and gravity yield values of about 4 km, but where coseismic fault slip extends to a depth of 10–15 km. Thus much of the seismogenic crust must weaken substantially during the life of active faults, causing the fault-bounded basins to narrow over time.

Journal of Geophysical Research Solid Earth↗
USGS2012

Geochemistry of some gases in hydrothermal fluids from the southern Juan de Fuca Ridge

Five samples of hydrothermal fluids from two vent areas on the southern Juan de Fuca Ridge were analyzed for dissolved gases. Concentrations in the end-member hydrothermal fluid of H 2 (270–527 μmol/kg), CH 4 (82–118 μmol/kg), and CO 2 (3920–4460 μmol/kg) are well above values in ambient seawater and are similar to concentrations reported for other ridge crest hydrothermal systems. The carbon isotopic ratios of the CH 4 (δ 13 C = −17.8 to −20.8) and CO 2 (δ 13 C = −3.6 to −4.7) suggest that at least some of the CH 4 and CO 2 in the fluids is basalt-derived. The range of δ 13 C values for the basalt-derived CO 2 is −6.8 to −9.7, calculated by assuming conservation of recharge ΣCO 2 during hydrothermal circulation. Apparent temperatures of equilibration between the CH 4 and the basalt-derived CO 2 range from 640°C to 750°C. Small amounts of ethane (C 2 H 6 /CH 4 ≅ 0.9 × 10 −3 −2.2 × 10 −3 ), propane, and butane detected in the samples may also have formed in the basalt. One sample of almost pure (95.5%) hydrothermal fluid contained a significant fraction, up to 63% and 74%, respectively, of the recharge Ar and N 2 . This suggests that the fluid has not undergone extensive vapor-liquid phase separation.

Journal of Geophysical Research Solid Earth↗
USGS2012

Paleomagnetism of the Oligocene Kalamazoo Tuff: implications for middle Tertiary extension in east central Nevada

The Oligocene Kalamazoo Tuff (≃35 Ma) was sampled for paleomagnetic analysis across a 100-km-wide zone of highly extended crust in east central Nevada to estimate between-site vertical axis rotations and thus the relative importance of strike-slip faulting to the mechanism of extension. Subordinate sampling was also undertaken in a younger hornblende dacite intrusion (≃35 Ma) and the overlying tuff of North Creek (≃35 Ma). In general, ash flow tuffs cool and are magnetized within a short period of time relative to geomagnetic secular variation and their eutaxitic structures typically indicate the paleohorizontal. These attributes of ash flow tuffs, in addition to their generally widespread occurrence, make them excellent candidates for combined paleomagnetic and structural studies of regional deformation. However, in this study a number of problems associated with the paleomagnetism of the Kalamazoo Tuff are documented; they include topography-related variations of the eutaxitic structures, secular variation during slow cooling, postblocking rheomorphic flow, and local remagnetization of the tuff by younger igneous activity. The tilt-corrected data, with these sources of error reduced or eliminated, exhibit a 28° ± 12° clockwise rotation of the Schell Creek Range relative to the Kern Mountains region. This rotation implies differential extension accommodated by strike-slip faulting or N-S shortening. The paleomagnetic results also suggest that large changes in strike of layered units near faults with presumed strike-slip movement need not be the result of oroclinal bending, but could result from superimposed sets of orthogonal normal faults.

Journal of Geophysical Research Solid Earth↗
USGS2012

Columbia Glacier, Alaska: Changes in velocity 1977-1986

The Columbia Glacier, a grounded, iceberg-calving tidewater glacier near Valdez, Alaska, began to retreat about 1977. Drastic retreat occurred in 1984, and by early 1986, retreat amounted to 2 km. The glacier has thinned more than 100 m since 1974 at a point 4 km behind the 1974 terminus position. Between 1977 and 1985 the lower glacier ice velocity increased from 3–8 m/d to 10–15 m/d. Ice velocity in the region 0.5 km above the terminus was highest near the time the glacier was most receded (late fall), and lowest near the time of maximum length (early summer), for years 1977–1982. Velocity in the region 52–57 from the head of the glacier was highest in mid-spring, and lowest in early fall from 1977 to 1985. Through the years 1983–1985, the dates of maximum and minimum velocities within 0.5 km of the receding terminus tended toward the dates of the 52–57 km maximum and minimums. This occurred because as the terminus receded, it was no longer strongly influenced by the reverse slope of the terminal moraine shoal. Velocities near the terminus fluctuated by 2–3 m/d during summer and fall, when liquid water input was variable, and were relatively constant during winter. Hourly variations in ice velocities are controlled by liquid water input to the glacier hydraulic system and tide stage. Velocity increases near periods of high surface water input and decreases during periods of high tide as a result of hydrostatic back pressure.

Journal of Geophysical Research Solid Earth↗
USGS2012

Usefulness of weak bands in midinfrared remote sensing of particulate planetary surfaces

Midinfrared (2.5–25 μm) reflectance spectra of minerals are often used to predict emittance qualitatively. These spectra display weak overtone and combination tone bands, which may be as diagnostic of composition as the strong fundamental molecular vibration bands usually considered for remote sensing applications, but which have been widely ignored. However, unlike the strong bands, the contrast of weak bands relative to the continuum in the midinfrared usually does not decrease with decreasing particle size, but typically increases. To illustrate this behavior, transmittance and reflectance spectra of calcite and quartz are presented for the wavelength range from 4000 to 400 cm −1 (2.5–25 μm). It is the purpose of this paper to point out that these weak bands are potentially useful for compositional remote sensing of particulate planetary surfaces; this will require completion of supporting laboratory studies to document the occurrence of weak bands.

Journal of Geophysical Research Solid Earth↗
USGS2020

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

Wyoming↗
USGS2012

Interpretation of gravity data in a complex volcano-tectonic setting, southwestern Nevada

This regional gravity study, based on an irregular 2-km data grid, was conducted during the past few years at Yucca Mountain, southern Nye County, Nevada, as part of a program to locate a suitable repository for high-level nuclear waste. About 100 surface rock samples, three borehole gamma-gamma logs, and one borehole gravity study provide excellent density control. A nearly linear increase in density of 0.26 g/cm 3 per kilometer of depth is indicated in the thick tuff sequences that underlie the mountain. Isostatic and 2.0-g/cm 3 Bouguer corrections were applied to the observed gravity values to remove regional gradients and topographic effects, respectively. The Bare Mountain gravity high, with an isostatic anomaly maximum of 48 mGal, is connected with a greater gravity high over the Funeral Mountains, to the southwest; together, these highs result from a continuous block of dense, metamorphosed Precambrian and Paleozoic rocks that stretches across much of the Walker Lane from the east edge of Death Valley to Bare Mountain. The Calico Hills gravity high appears more likely to originate from a northeast trending buried ridge of Paleozoic rocks that extends southwestward beneath Busted Butte, 5 km southeast of the proposed repository, where two- and three-dimensional modeling indicates that the pre-Cenozoic rocks lie less than 1000 m beneath the surface. Tuff, at least 4000 m thick, fills a large steep-sided depression in the pretuff rocks beneath Yucca Mountain and Crater Flat. The gravity low and the thick tuff section lie within a large collapse area that includes the Crater Flat-Timber Mountain-Silent Canyon caldera complexes. Gravity lows in Crater Flat itself are interpreted to coincide with the source areas of the Prow Pass Member, the Bullfrog Member, and the Tram Member of the Crater Flat Tuff; these source areas add nearly 350 km 2 to the previously recognized extent of the local caldera complexes. Southward extension of the broad gravity low associated with Crater Flat into the Amargosa Desert is evidence for sector graben-type collapse segments related to the formation of the Timber Mountain caldera and superimposed on the other volcanic and extensional structures within Crater Flat.

Journal of Geophysical Research Solid Earth↗
USGS2012

Age of the Peach Springs Tuff, southeastern California and western Arizona

Sanidine separates from pumice of the early Miocene Peach Springs Tuff are concordantly dated at 18.5±0.2 Ma by two isotopic techniques. The Peach Springs Tuff is the only known unit that can be correlated between isolated outcrops of Miocene strata from the central Mojave Desert of southeastern California to the western Colorado Plateau in Arizona, across five structural provinces, a distance of 350 km. Thus the age of the Peach Springs Tuff is important to structural and paleogeographic reconstructions of a large region. Biotite and sanidine separates from bulk samples of the Peach Springs Tuff from zones of welding and vapor-phase alteration have not produced consistent ages by the K-Ar method. Published ages of mineral separates from 17 localities ranged from 16.2 to 20.5 Ma. Discordant 40 Ar/ 39 Ar incremental release spectra were obtained for one biotite and two of the sanidine separates. Ages that correspond to the last gas increments are as old as 27 Ma. The 40 Ar/ 39 Ar incremental release determinations on sanidine separated from blocks of Peach Springs Tuff pumice yield ages of 18.3±0.3 and 18.6±0.4 Ma. Laser fusion measurements yield a mean age of 18.51±0.10. The results suggest that sanidine and biotite K-Ar ages older than about 18.5 Ma are due to inherited Ar from pre-Tertiary contaminants, which likely were incorporated into the tuff during deposition. Sanidine K-Ar ages younger than 18 Ma probably indicate incomplete extraction of radiogenic 40 Ar, whereas laser fusion dates of biotite and hornblende younger than 18 Ma likely are due to postdepositional alteration. Laser fusion ages as high as 19.01 Ma on biotite grains from pumice suggest that minerals from pre-Tertiary country rocks also were incorporated in the magma chamber.

Journal of Geophysical Research Solid Earth↗
USGS2012

An assessment of the accuracy of the geodetic measurements that define the southern California uplift

Examination of the charge that the geodetic measurements which define the southern California uplift are seriously flawed by height-dependent systematic errors indicates that this charge is unfounded. Our review of those factors designed to minimize the effects of rod error shows that systematic errors attributable to rod miscalibration are trivial and would tend to cancel during the particularly pertinent period 1955–1965. Comparisons between the results of measurements in which the correlation between topography and signal is generally poor reveal large aseismic tilts in a number of places within and around the margins of the uplift. Especially significant in this context are the results of preuplift and postuplift levelings over routes characterized by diverse length, topography, and atmospheric conditions that produce closely matching, temporally equivalent heights for a representative bench mark within the uplift. Statistical analyses of the short-wavelength components of uplift signal and elevation along a frequently repeated survey line indicate variable correlations between the two. However, the spatial and temporal patterns of these correlations are inconsistent with their attribution to height-dependent systematic errors. Statistically significant short-wavelength correlations along this survey line are reasonably explained in part as due to real movement associated with differentially subsiding bench marks, and there is a strong likelihood that this movement tends to dominate many of the correlations.

Journal of Geophysical Research Solid Earth↗
USGS2020

Plastic faulting in ice

Plastic faulting is a brittle‐like failure phenomenon exhibited by water ice and several other rock types under confinement. It is suspected to be the mechanism of deep earthquakes and extreme cases of shear localization in shallow rocks. Unlike ordinary Coulombic failure, plastic faulting is characterized by a pressure‐independent failure strength and fault plane oriented 45° to maximum principal stress. To research the question of how the instability initiates, we conducted over 50 constant‐displacement‐rate experiments on polycrystalline ice (phases Ih and II) near the brittle‐to‐ductile (B‐D) transition, at confining pressures P = 0–300 MPa, applied strain rates = 5 × 10 −5 – 7 × 10 −3 s −1 , temperatures T = 105–233 K, and mean grain sizes d = 0.25–1.18 mm. We find that (1) the width of the B‐D transition in variable space is vanishingly narrow, to the point of appearing as a crossover, (2) a plastic fault plane, once formed, is not a zone of subsequent weakness, (3) distributed ice I→II phase transformation in small amounts (<1 vol%) shows no causal relationship to subsequent failure, and (4) plastic faulting also occurs in ice II. We hypothesize that the elusive nucleating “trigger” parallels that of metals and ceramics undergoing severe plastic deformation, wherein transient local structural rearrangement occurs, in turn causing material strength to drop to a level sufficiently low, in a volume sufficiently large, that adiabatic instability is nucleated. Our results do not require and often are inconsistent with phase transformation. Plastic faulting may therefore be available to all solids undergoing severe deformation, and its appearance in so few is simply the result of insufficiently extreme conditions.

Journal of Geophysical Research- Solid Earth↗
USGS2023

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma↗
USGS2012

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗
USGS2021

Cataloging tectonic tremor energy radiation in the Cascadia subduction zone

For the past ∼12 years the Pacific Northwest Seismic Network has been automatically detecting and locating tectonic tremor across the Cascadia subduction zone, resulting in a catalog of more than 500,000 tremor epicenters to date, which has served as a valuable resource for tremor and slip research. This manuscript presents an updated methodology for routine tremor detection in Cascadia and a new catalog of over 180,000 tremor epicenters including amplitudes detected along the subduction zone margin from 2017 to 2021. The events are detected via cross-correlation of continuous vertical envelope data of 128 stations from northern California to northern Vancouver Island. The modified approach results in less scatter and a 55% increase in detected epicenters than previously observed, as well as a newly identified tremor source offset updip from the main tremor and slip region at the southern edge of the subduction zone. Radiated seismic energy in the 1.5–5 Hz band is used to assign epicenters an energy magnitude ( M eL ), which is calibrated to the M L of local earthquakes. Southern Cascadia is most active, but the highest tremor energy rates occur in northern Cascadia. Tremor in central Cascadia is systematically weaker and less frequent. Individual epicenter magnitudes range from ∼0.5–2 and spatiotemporally cluster into 1,060 swarms with cumulative M eL ranging from ∼0.8 to 3.7. The swarms reflect underlying slow slip events and occur with an earthquake-like energy distribution with a b value ∼1. Tremor epicenters, however, follow a tapered Gutenberg-Richter distribution with high b values, suggesting individual tremor bursts and their constituent low-frequency earthquakes are fault-dimension limited.

California, Oregon, Washington↗
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↗
USGS2022

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗
USGS2025

New paleomagnetic constraints on the eruption timing, stratigraphy, and post-emplacement deformation of the Picture Gorge Basalt within the Columbia River Basalt Group

The Picture Gorge Basalt (PGB) is part of the Miocene Columbia River Basalt Group (CRBG). The PGB, which outcrops in eastern Oregon, is considered coincident in time with the voluminous Grande Ronde Basalt. New radiometric ages have expanded the age‐range of the PGB, suggesting it began erupting prior to the Steens Basalt to its south and continued until after cessation of the Grande Ronde Basalt eruptions, an interval of 1.5 Ma. However, the existing paleomagnetism of the PGB implies this eruption timeline is an overestimate. To reconcile the radiometric and paleomagnetic timescales for the PGB, we conducted a paleomagnetic study on sections of the PGB to construct a detailed, high‐quality magnetostratigraphy. Our data indicate the stratigraphically lowest lava flows in the PGB are of reversed polarity, revealing a new paleomagnetic transition with the PGB and a reversed (R)–normal (N)–reversed (R) sequence. This suggests one of two timeline possibilities for PGB volcanism: (a) eruptions began and during through CRBG polarity chrons R0–N0–R1, penecontemporaneous with Steens Basalt, or (b) eruptions began and persisted during CRBG polarity chrons R1–N1–R2. Our work supports a longer interval of PGB volcanism than was suggested by previous paleomagnetic data but is at odds with the suggestion that PGB eruptions lasted through the entire main CRBG. We favor a scenario wherein PGB eruptions begin with R0 and continue into the R1 paleomagnetic interval. The paleomagnetic results also record a ∼18° vertical‐axis rotation of east‐central Oregon after ∼16 Ma with respect to stable North America.

Oregon↗