USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Observer”

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 1,225 records · Page 68Linked to original sources

Self-similar rupture implied by scaling properties of volcanic earthquakes occurring during the 2004-2008 eruption of Mount St. Helens, Washington

We analyze a group of 6073 low-frequency earthquakes recorded during a week-long temporary deployment of broadband seismometers at distances of less than 3 km from the crater at Mount St. Helens in September of 2006. We estimate the seismic moment ( M 0 ) and spectral corner frequency ( f 0 ) using a spectral ratio approach for events with a high signal-to-noise (SNR) ratio that have a cross-correlation coefficient of 0.8 or greater with at least five other events. A cluster analysis of cross-correlation values indicates that the group of 421 events meeting the SNR and cross-correlation criteria forms eight event families that exhibit largely self-similar scaling. We estimate the M 0 and f 0 values of the 421 events and calculate their static stress drop and scaled energy ( E R / M 0 ) values. The estimated values suggest self-similar scaling within families, as well as between five of eight families (i.e., and constant). We speculate that differences in scaled energy values for the two families with variable scaling may result from a lack of resolution in the velocity model. The observation of self-similar scaling is the first of its kind for such a large group of low-frequency volcanic tectonic events occurring during a single active dome extrusion eruption.

Journal of Geophysical Research B: Solid Earth↗

Thresholds for post-wildfire debris flows: Insights from the Pinal Fire, Arizona, USA

Wildfire significantly alters the hydrologic properties of a burned area, leading to increases in overland flow, erosion, and the potential for runoff-generated debris flows. The initiation of debris flows in recently burned areas is well-characterized by rainfall intensity-duration (ID) thresholds. However, there is currently a paucity of data quantifying the rainfall intensities required to trigger post-wildfire debris flows, which limits our understanding of how and why rainfall ID thresholds vary in different climatic and geologic settings. In this study, we monitored debris-flow activity following the Pinal Fire in central Arizona, which differs from both a climatic and hydrogeomorphic perspective from other regions in the western U.S. where ID thresholds for post-wildfire debris flows are well-established, namely the Transverse Ranges of southern CA. Since the peak rainfall intensity within a rainstorm may exceed the rainfall intensity required to trigger a debris flow, the development of robust rainfall ID thresholds requires knowledge of the timing of debris flows within rainstorms. Existing post-wildfire debris-flow studies in Arizona only constrain the peak rainfall intensity within debris-flow-producing storms, which may far exceed the intensity that actually triggered the observed debris flow. In this study, we used pressure transducers within 5 burned drainage basins to constrain the timing of debris flows within rainstorms. Rainfall ID thresholds derived here from triggering rainfall intensities are, on average, 22 mm/h lower than ID thresholds derived under the assumption that the triggering intensity is equal to the maximum rainfall intensity recorded during a rainstorm. We then use a hydrologic model to demonstrate that the magnitude of the 15-minute rainfall ID threshold at the Pinal Fire site is associated with the rainfall intensity required to exceed a recently proposed dimensionless discharge threshold for debris-flow initiation. Model results further suggest that previously observed differences in regional ID thresholds between Arizona and the San Gabriel Mountains of southern CA may be attributed, in large part, to differences in the hydraulic properties of burned soils.

Arizona↗

Runoff-generated debris flows: observations and modeling of surge initiation, magnitude, and frequency

Runoff during intense rainstorms plays a major role in generating debris flows in many alpine areas and burned steeplands. Yet compared to debris flow initiation from shallow landslides, the mechanics by which runoff generates a debris flow are less understood. To better understand debris flow initiation by surface water runoff, we monitored flow stage and rainfall associated with debris flows in the headwaters of two small catchments: a bedrock-dominated alpine basin in central Colorado (0.06 km 2 ) and a recently burned area in southern California (0.01 km 2 ). We also obtained video footage of debris flow initiation and flow dynamics from three cameras at the Colorado site. Stage observations at both sites display distinct patterns in debris flow surge characteristics relative to rainfall intensity (I). We observe small, quasiperiodic surges at low I; large, quasiperiodic surges at intermediate I; and a single large surge followed by small-amplitude fluctuations about a more steady high flow at high I. Video observations of surge formation lead us to the hypothesis that these flow patterns are controlled by upstream variations in channel slope, in which low-gradient sections act as “sediment capacitors,” temporarily storing incoming bed load transported by water flow and periodically releasing the accumulated sediment as a debris flow surge. To explore this hypothesis, we develop a simple one-dimensional morphodynamic model of a sediment capacitor that consists of a system of coupled equations for water flow, bed load transport, slope stability, and mass flow. This model reproduces the essential patterns in surge magnitude and frequency with rainfall intensity observed at the two field sites and provides a new framework for predicting the runoff threshold for debris flow initiation in a burned or alpine setting.

California;Colorado↗

Identifying calcium sources at an acid deposition-impacted spruce forest: A strontium isotope, alkaline earth element multi-tracer approach

Depletion of calcium from forest soils has important implications for forest productivity and health. Ca is available to fine feeder roots from a number of soil organic and mineral sources, but identifying the primary source or changes of sources in response to environmental change is problematic. We used strontium isotope and alkaline earth element concentration ratios of trees and soils to discern the record of Ca sources for red spruce at a base-poor, acid deposition-impacted watershed. We measured 87 Sr/ 86 Sr and chemical compositions of cross-sectional stemwood cores of red spruce, other spruce tissues and sequential extracts of co-located soil samples. 87 Sr/ 86 Sr and Sr/Ba ratios together provide a tracer of alkaline earth element sources that distinguishes the plant-available fraction of the shallow organic soils from those of deeper organic and mineral soils. Ca/Sr ratios proved less diagnostic, due to within-tree processes that fractionate these elements from each other. Over the growth period from 1870 to 1960, 87 Sr/ 86 Sr and Sr/Ba ratios of stemwood samples became progressively more variable and on average trended toward values that considered together are characteristic of the uppermost forest floor. In detail the stemwood chemistry revealed an episode of simultaneous enhanced uptake of all alkaline earth elements during the growth period from 1930 to 1960, coincident with reported local and regional increases in atmospheric inputs of inorganic acidity. We attribute the temporal trends in stemwood chemistry to progressive shallowing of the effective depth of alkaline earth element uptake by fine roots over this growth period, due to preferential concentration of fine roots in the upper forest floor coupled with reduced nutrient uptake by roots in the lower organic and upper mineral soils in response to acid-induced aluminum toxicity. Although both increased atmospheric deposition and selective weathering of Ca-rich minerals such as apatite provide possible alternative explanations of aspects of the observed trends, the chemical buffering capacity of the forest floor-biomass pool limits their effectiveness as causal mechanisms.

Biogeochemistry↗

Modeling the height of young forests regenerating from recent disturbances in Mississippi using Landsat and ICESat data

Many forestry and earth science applications require spatially detailed forest height data sets. Among the various remote sensing technologies, lidar offers the most potential for obtaining reliable height measurement. However, existing and planned spaceborne lidar systems do not have the capability to produce spatially contiguous, fine resolution forest height maps over large areas. This paper describes a Landsat–lidar fusion approach for modeling the height of young forests by integrating historical Landsat observations with lidar data acquired by the Geoscience Laser Altimeter System (GLAS) instrument onboard the Ice, Cloud, and land Elevation (ICESat) satellite. In this approach, “young” forests refer to forests reestablished following recent disturbances mapped using Landsat time-series stacks (LTSS) and a vegetation change tracker (VCT) algorithm. The GLAS lidar data is used to retrieve forest height at sample locations represented by the footprints of the lidar data. These samples are used to establish relationships between lidar-based forest height measurements and LTSS–VCT disturbance products. The height of “young” forest is then mapped based on the derived relationships and the LTSS–VCT disturbance products. This approach was developed and tested over the state of Mississippi. Of the various models evaluated, a regression tree model predicting forest height from age since disturbance and three cumulative indices produced by the LTSS–VCT method yielded the lowest cross validation error. The R 2 and root mean square difference (RMSD) between predicted and GLAS-based height measurements were 0.91 and 1.97 m, respectively. Predictions of this model had much higher errors than indicated by cross validation analysis when evaluated using field plot data collected through the Forest Inventory and Analysis Program of USDA Forest Service. Much of these errors were due to a lack of separation between stand clearing and non-stand clearing disturbances in current LTSS–VCT products and difficulty in deriving reliable forest height measurements using GLAS samples when terrain relief was present within their footprints. In addition, a systematic underestimation of about 5 m by the developed model was also observed, half of which could be explained by forest growth that occurred between field measurement year and model target year. The remaining difference suggests that tree height measurements derived using waveform lidar data could be significantly underestimated, especially for young pine forests. Options for improving the height modeling approach developed in this study were discussed.

Remote Sensing of Environment↗

Reconstructing the velocity and deformation of a rapid landslide using multiview video

Noncontact measurements of spatially varied ground surface deformation during landslide motion can provide important constraints on landslide mechanics. Here, we present and test a new method for extracting measurements of rapid landslide surface displacement and velocity (accelerations of approximately 1 m/s 2 ) using sequences of stereo images obtained from a pair of inexpensive, stationary 4K video cameras with nominal frame rates of 29.97 Hz. The method combines elements of Structure from Motion with those of optical flow to extract data on 3‐D evolution of the ground surface during slope failure. We apply the method to an experiment at the U.S. Geological Survey debris‐flow flume in which a high‐speed, liquefying landslide was triggered by gradually adding water to a 6‐m 3 prism of loosely packed sediment on a 31° slope. Strip‐scanning lidar measurements made during the experiment corroborate our video‐based measurements, but the latter encompassed the entire landslide surface and were much lower in cost. Our video‐based measurements enabled computation of depth‐integrated landslide dilation/contraction rates. The range of computed rates was within the ranges inferred from independent measurements of evolving pore water pressures and reasonable estimates of the hydraulic permeability of the sediment. Dilation and contraction rates play a crucial role in landslide mechanics. The dilation and contraction we observe contradict the incompressible flow assumption used in many studies that have employed noncontact methods to infer landslide properties.

Oregon↗

A high-resolution seismic reflection/refraction study of the Chugach-Peninsular terrane boundary, southern Alaska

We present results from a high-resolution seismic refraction analysis of the shallow (approximately 2 km) crustal structure along the 107-km-long Trans-Alaska Crustal Transect Chugach reflection line in southern Alaska and a comparison with laboratory measurements of field samples. The refraction analysis includes the two-dimensional interpretation of several thousand first- and secondary-arrival travel times digitized from 1024-channel split-spread common shot gathers. The velocity model derived from this analysis better defines the location and geometry of terrane boundaries than does the normal incidence reflection section and agrees well with surface mapping of lithologies. Furthermore, the model predicts travel times within 100 ms of the reflection times recorded from the base of the Quaternary on the Chugach reflection section. Thicknesses of Quaternary deposits, with velocities between 1.1 and 2.0 km/s, correlate inversely with the quantity of observed lower crustal reflections on the Chugach section, suggesting that the presence or absence of these sediments in sufficient thickness exerted primary control on the quality of the deeper portion of the section. There is a significant velocity contrast between crystalline rocks across the Border Ranges fault (5.0 versus 5.6 km/s), the major contact between the Chugach and Peninsular terranes, in agreement with laboratory measurements of field specimens. In the Peninsular terrane the modeling indicates that an unnamed fault delimiting the southern flank of the Copper River Basin dips steeply northward at 50° and has about 1300 m of vertical offset. Laboratory measurements document a maximum velocity anisotropy of 20% for phyllitic schists of the Valdez Group in the Chugach terrane. In agreement with the observed E-W strike and near-vertical dip of the Valdez Group, we determined a significant (14%) velocity anisotropy for ray paths oriented N-S versus NE-SW.

Journal of Geophysical Research Solid Earth↗

Interpreting and reporting 40Ar/39Ar geochronologic data

The 40 Ar/ 39 Ar dating method is among the most versatile of geochronometers, having the potential to date a broad variety of K-bearing materials spanning from the time of Earth’s formation into the historical realm. Measurements using modern noble-gas mass spectrometers are now producing 40 Ar/ 39 Ar dates with analytical uncertainties of ∼0.1%, thereby providing precise time constraints for a wide range of geologic and extraterrestrial processes. Analyses of increasingly smaller subsamples have revealed age dispersion in many materials, including some minerals used as neutron fluence monitors. Accordingly, interpretive strategies are evolving to address observed dispersion in dates from a single sample. Moreover, inferring a geologically meaningful “age” from a measured “date” or set of dates is dependent on the geological problem being addressed and the salient assumptions associated with each set of data. We highlight requirements for collateral information that will better constrain the interpretation of 40 Ar/ 39 Ar data sets, including those associated with single-crystal fusion analyses, incremental heating experiments, and in situ analyses of microsampled domains. To ensure the utility and viability of published results, we emphasize previous recommendations for reporting 40 Ar/ 39 Ar data and the related essential metadata, with the amendment that data conform to evolving standards of being findable, accessible, interoperable, and reusable (FAIR) by both humans and computers. Our examples provide guidance for the presentation and interpretation of 40 Ar/ 39 Ar dates to maximize their interdisciplinary usage, reproducibility, and longevity.

GSA Bulletin↗

Sedimentary response to orogenic exhumation in the northern Rocky Mountain Basin and Range province, Flint Creek basin, west-central Montana

Middle Eocene through Upper Miocene sedimentary and volcanic rocks of the Flint Creek basin in western Montana accumulated during a period of significant paleoclimatic change and extension across the northern Rocky Mountain Basin and Range province. Gravity modelling, borehole data, and geologic mapping from the Flint Creek basin indicate that subsidence was focused along an extensionally reactivated Sevier thrust fault, which accommodated up to 800 m of basin fill while relaying stress between the dextral transtensional Lewis and Clark lineament to the north and the Anaconda core complex to the south. Northwesterly paleocurrent indicators, foliated metamorphic lithics, 64 Ma (40Ar/39Ar) muscovite grains, and 76 Ma (U-Pb) zircons in a ca. 27 Ma arkosic sandstone are consistent with Oligocene exhumation and erosion of the Anaconda core complex. The core complex and volcanic and magmatic rocks in its hangingwall created an important drainage divide during the Paleogene shedding detritus to the NNW and ESE. Following a major period of Early Miocene tectonism and erosion, regional drainage networks were reorganized such that paleoflow in the Flint Creek basin flowed east into an internally drained saline lake system. Renewed tectonism during Middle to Late Miocene time reestablished a west-directed drainage that is recorded by fluvial strata within a Late Miocene paleovalley. These tectonic reorganizations and associated drainage divide explain observed discrepancies in provenance studies across the province. Regional correlation of unconformities and lithofacies mapping in the Flint Creek basin suggest that localized tectonism and relative base level fluctuations controlled lithostratigraphic architecture.

Canadian Journal of Earth Sciences↗

Rayleigh wave amplitude uncertainty across the Global Seismographic Network and potential implications for global tomography

The Global Seismographic Network (GSN) is a multiuse, globally distributed seismic network used by seismologists, to both characterize earthquakes and study the Earth’s interior. Most stations in the network have two collocated broadband seismometers, which enable network operators to identify potential metadata and sensor issues. In this study, we investigate the accuracy with which surface waves can be measured across the GSN, by comparing waveforms of vertical‐component Rayleigh waves from M w "> M w Mw 6 and larger events between collocated sensor pairs. We calculate both the amplitude deviation and correlation coefficient between waveforms at sensor pairs. In total, we make measurements on over 670,000 event–station pairs from events that occurred from 1 January 2010 to 1 January 2020. We find that the average sensor‐pair amplitude deviation, and, therefore, GSN calibration level, is, approximately, 4% in the 25–250 s period band. Although, we find little difference in sensor‐pair amplitude deviations as a function of period across the entire network, the amount of useable data decreases rapidly as a function of increasing period. For instance, we determined that just over 12% of records at 250 s period provided useable recordings (e.g., sensor‐pair amplitude deviations of less than 20% and sensor‐pair correlation greater than 0.95). We then use these amplitude‐estimate deviations to identify how data coverage and quality could be limiting our ability to invert for whole Earth 3D attenuation models. We find an increase in the variance of our attenuation models with increasing period. For example, our degree 12 attenuation inversion at 250 s period shows 32% more variance than our degree 12 attenuation model at 25 s. This indicates that discrepancies of deep‐mantle tomography between studies could be the result of these large uncertainties. Because these high uncertainties arise from limited, high‐quality observations of long‐period ( ⁠ > 100    s ">

Bulletin of the Seismological Society of America↗

Teleseismic waves reveal anisotropic poroelastic response of wastewater disposal reservoir

Connecting earthquake nucleation in basement rock to fluid injection in basal, sedimentary reservoirs, depends heavily on choices related to the poroelastic properties of the fluid-rock system, thermo-chemical effects notwithstanding. Direct constraints on these parameters outside of laboratory settings are rare, and it is commonly assumed that the rock layers are isotropic. With the Arbuckle wastewater disposal reservoir in Osage County, Oklahoma, high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic; this includes evidence of static shifts in pressure that presumably relate to changes in local permeability. The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures. Using dynamic strains from a nearby borehole strainmeter, we show that the ratio of shear to volumetric strain coupling is ∼ 0.41 "> ∼ 0.41 ∼0.41 which implies a mean Skempton's coefficient of A = 0.24 "> A = 0.24 A=0.24 over the plausible range of the undrained Poisson's ratio. Since these observations are made at relatively low confining pressure and differential stress, we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.

Earth and Planetary Physics↗

New insights into Kilauea's volcano dynamics brought by large-scale relative relocation of microearthquakes

We investigated the microseismicity recorded in an active volcano to infer information concerning the volcano structure and long-term dynamics, by using relative relocations and focal mechanisms of microearthquakes. There were 32,000 earthquakes of the Mauna Loa and Kilauea volcanoes recorded by more than eight stations of the Hawaiian Volcano Observatory seismic network between 1988 and 1999. We studied 17,000 of these events and relocated more than 70%, with an accuracy ranging from 10 to 500 m. About 75% of these relocated events are located in the vicinity of subhorizontal decollement planes, at a depth of 8-11 km. However, the striking features revealed by these relocation results are steep southeast dipping fault planes working as reverse faults, clearly located below the decollement plane and which intersect it. If this decollement plane coincides with the pre-Mauna Loa seafloor, as hypothesized by numerous authors, such reverse faults rupture the pre-Mauna Loa oceanic crust. The weight of the volcano and pressure in the magma storage system are possible causes of these ruptures, fully compatible with the local stress tensor computed by Gillard et al. [1996]. Reverse faults are suspected of producing scarps revealed by kilometer-long horizontal slip-perpendicular lineations along the decollement surface and therefore large-scale roughness, asperities, and normal stress variations. These are capable of generating stick-slip, large-magnitude earthquakes, the spatial microseismic pattern observed in the south flank of Kilauea volcano, and Hilina-type instabilities. Rupture intersecting the decollement surface, causing its large-scale roughness, may be an important parameter controlling the growth of Hawaiian volcanoes.

Journal of Geophysical Research B: Solid Earth↗

Rise of a variable-viscosity fluid in a steadily spreading wedge-shaped conduit with accreting walls

Relatively rigid plates making up the outer 50 to 100 km of the Earth are steadily separating from one another along narrow globe-circling zones of submarine volcanism, the oceanic spreading centers. Continuity requires that the viscous underlying material rise beneath spreading centers and accrete onto the steadily diverging plates. It is likely that during the rise the viscosity changes systematically and that the viscous tractions exerted on the plates contribute to the unique pattern of submarine mountains and earthquake faults observed at spreading centers. The process is modeled by viscous creep in a wedge-shaped conduit (with apex at the sea floor) in which the viscosity varies as r m where r is distance from the apex and m is a parameter. For these conditions, the governing differential equations take a simple form. The solution for the velocity is independent of r and of the sign of m . As viscous stresses vary as r m-1 , the pattern of stress on the conduit wall is sensitive to viscosity variation. For negative m , the viscous pressure along the base of the conduit is quite uniform; for positive m , it falls toward zero in the axial region as the conduit base widens. For small opening angles, viscous forces push the plates apart, and for large ones, they oppose plate separation. Though highly idealized, the solution provides a tool for investigating tectonic processes at spreading centers.

Journal of Research of the U.S. Geological Survey↗

Failure in laboratory fault models in triaxial tests

A model of a fault in the Earth is a sand-filled saw cut in a granite cylinder subjected to a triaxial test. The saw cut is inclined at an angle α to the cylinder axis, and the sand filling is intended to represent gouge. The triaxial test subjects the granite cylinder to a constant confining pressure and increasing axial stress to maintain a constant rate of shortening of the cylinder. The required axial stress increases at a decreasing rate to a maximum, beyond which a roughly constant axial stress is sufficient to maintain the constant rate of shortening. Such triaxial tests were run for saw cuts inclined at angles α of 20°, 25°, 30°, 35°, 40°, 45°, and 50° to the cylinder axis, and the apparent coefficient of friction μ a (ratio of the shear stress to the normal stress, both stresses resolved onto the saw cut) at failure was determined. Subject to the assumption that the observed failure involves slip on Coulomb shears (orientation unspecified), the orientation of the principal compression axis within the gouge can be calculated as a function of μ a for a given value of the coefficient of internal friction μ i . The rotation of the principal stress axes within the gouge in a triaxial test can then be followed as the shear strain across the gouge layer increases. For μ i ∼ 0.8, an appropriate value for highly sheared sand, the observed values μ a imply that the principal axis of compression within the gouge rotates so as to approach being parallel to the cylinder axis for all saw cut angles (20° < α < 50°). In the limiting state (principal compression axis parallel to cylinder axis) the stress state in the gouge layer would be the same as that in the granite cylinder, and the failure criterion would be independent of the saw cut angle.

Journal of Geophysical Research B: Solid Earth↗

Subsurface energy storage and transport for solar-powered geysers on Triton

The location of active geyser-like eruptions and related features close to the current subsolar latitude on Triton suggests a solar energy source for these phenomena. Solidstate greenhouse calculations have shown that sunlight can generate substantially elevated subsurface temperatures. A variety of models for the storage of solar energy in a sub-greenhouse layer and for the supply of gas and energy to a geyser are examined. "Leaky greenhouse" models with only vertical gas transport are inconsistent with the observed upper limit on geyser radius of ∼ 1.5 kilometers. However, lateral transport of energy by gas flow in a porous N 2 layer with a block size on the order of a meter can supply the required amount of gas to a source region ∼1 kilometer in radius. The decline of gas output to steady state may occur over a period comparable with the inferred active geyser lifetime of five Earth years. The required subsurface permeability may be maintained by thermal fracturing of the residual N2 polar cap. A lower limit on geyser source radius of ∼50 to 100 meters predicted by a theory of negatively buoyant jets is not readily attained.

Science↗

Controls on debris‐flow initiation on burned and unburned hillslopes during an exceptional rainstorm in southern New Mexico, USA

AbstractUsing observations from 688 debris flows, we analyse the hydrologic and landscape characteristics that influenced debris‐flow initiation mechanisms and locations in a watershed that had been partially burned by the 2012 Whitewater‐Baldy Complex Fire in the Gila Mountains, southern New Mexico. Debris flows can initiate due to different processes. Slopes can fail as discrete landslides and then become fluidized and move downstream as debris flows (landslide initiated) or progressive bulking of sediment from a distributed area can become channelized and concentrated as it moves downslope (runoff generated). In this study, we have an unusual opportunity to investigate both types of debris‐flow initiation mechanisms in our observations of debris flows, triggered by an exceptional rainstorm in the autumn of 2013. Additionally, we compare our observations with those of a dataset of 1138 debris flows in the Colorado Front Range, triggered during the same weather system. We found that runoff‐generated debris flows dominated in burn areas, and runoff required to start these flows could be well characterized by the Shields stress. Landslide‐initiated debris flows were dominant in unburned areas. Debris‐flow densities were tied to total rainfall and precipitation intensities. Like the observations in the Colorado Front Range, debris‐flow initiation locations were found primarily in areas of relatively sparse vegetation on south‐facing slopes between 25 and 40°, and with upslope contributing areas less than 1000 m 2 . In terms of preferential locations for debris‐flow initiations, 2013 vegetation coverage, approximated by Green–Red Vegetation Index metrics, proved to be more influential than the 2012 burn‐severity designation. The uniformity of observations between our study area and those in the Colorado Front Range indicate that the underlying hydrologic and landscape patterns of the debris‐flow initiation locations documented in these studies could be applicable to the wider southwest and Rocky Mountain regions.

New Mexico↗

Experimental observations of pressure oscillations and flow regimes in an analogue volcanic system

Gas-liquid flows, designed to be analogous to those in volcanic conduits, are generated in the laboratory using organic gas-gum rosin mixtures expanding in a vertically mounted tube. The expanding fluid shows a range of both flow and pressure oscillation behaviors. Weakly supersaturated source liquids produce a low Reynolds number flow with foam expanding from the top surface of a liquid that exhibits zero fluid velocity at the tube wall; i.e., the conventional "no-slip" boundary condition. Pressure oscillations, often with strong long-period characteristics and consistent with longitudinal and radial resonant oscillation modes, are detected in these fluids. Strongly supersaturated source liquids generate more energetic flows that display a number of flow regimes. These regimes include a static liquid source, viscous flow, detached flow (comprising gas-pockets-at-wall and foam-in-gas annular flow, therefore demonstrating strong radial heterogeneity), and a fully turbulent transonic fragmented or mist flow. Each of these flow regimes displays characteristic pressure oscillations that can be related to resonance of flow features or wall impact phenomena. The pressure oscillations are produced by the degassing processes without the need of elastic coupling to the confining medium or flow restrictors and valvelike features. The oscillatory behavior of the experimental flows is compared to seismoacoustic data from a range of volcanoes where resonant oscillation of the fluid within the conduit is also often invoked as controlling the observed oscillation frequencies. On the basis of the experimental data we postulate on the nature of seismic signals that may be measured during large-scale explosive activity.

Journal of Geophysical Research B: Solid Earth↗

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗