SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23Linked to original sources

USGS2012

The glass transition in basalt

The glass transition has been experimentally detected in basalt as (1) an increase in the aggregate linear thermal expansion coefficient α L , (2) an abrupt change in the temperature dependence of Young's modulus dE / dT , and (3) a change in stress relaxation behavior that effectively separates the T > T G and T < T G creep regimes. Transition temperatures determined by the respective experimental methods are (1) 730° ± 15°C, by expansion dilatometry, (2) 725° ± 10°C, by acoustic spectroscopy, and (3) 725° ± 25°C, by stress relaxation. For olivine tholeiite from Kilauea Iki lava lake, Hawaii, thermal expansivity rises through the glass transition interval, attaining α υ = 39.0 × 10 −6 °C −1 at 850°C. For olivine tholeiites with moderate glass contents, the temperature dependence of the Young's modulus changes on either side of the inferred T G of 725°C, and dE / dT = −0.071 kbar °C −1 for T < T G , while dE / dT = − 0.21 kbar °C −1 for T > T G . Collectively, the mechanical results suggest that for Hawaiian olivine tholeiite at 1-atm pressure, the principal material responses are (1) elastic ( T ≤ 600°C), (2) reduced creep (600 < T < 725°C), (3) glass transition ( T ≃ 725°C), (4) enhanced creep (725 < T < 980°C), and (5) partial melt ( T > 980°C). Fracture surface morphologies developed during solidification suggest that the presence of the supercooled melt grain boundary phase may participate in the regulation of the thermal cracking process. Well-preserved fracture surfaces formed by incremental crack growth are found to be covered with striations that correspond to the inferred sequential stopping positions of the advancing fracture front. These striae may be rationalized in terms of experimental analogues that have been produced in viscoelastic polymers by cyclic tension-tension loading. The fracture surface morphology, mechanical data, and the controlled crack growth analogues suggest that thermal fracture in solidifying basalt is an incremental and cyclic process, involving three steps: (1) the accumulation of elastic strain energy in cooling rock at temperatures below that required for stress relaxation due to viscous flow in the intercrystalline liquid phase, (2) fracture at a Δ T determined primarily by the aggregate thermal expansion coefficient α υ and Young's modulus E , (3) the penetration by the advancing crack tip, of the thermal horizon capable of relaxing stress due to the creep of intercrystalline supercooled melt, producing the rough surface texture associated with the termination of a striation. Further crack growth must now await the migration of the solidus. The cycle then repeats. Striations measured in deep Hawaiian lava lakes have been compared with the crack advance increments expected in the vicinity of the glass transition, based on two tests: (1) thermal gradients measured in Kilauea Iki, combined with the mechanical properties of olivine tholeiite evaluated near T G , and (2) the crack advance required to match the recorded seismic stopping phases for prexisting cracks of the dimensions expected for Kilauea Iki. The observed versus predicted comparisons are (1) 31 versus 36 cm; and (2) 31 versus 30 cm. We envision this incremental crack growth process as contributing to the control on the downward movement of the thermal cracking front—and its associated hydrothermal circulation zone—in the upper portions of solidifying subaerial and submarine ponded basalt.

Journal of Geophysical Research Solid Earth↗

USGS2005

Distribution of aseismic slip rate on the Hayward fault inferred from seismic and geodetic data

We solve for the slip rate distribution on the Hayward fault by performing a least squares inversion,of geodetic and seismic data sets. Our analysis focuses on the northern 60 km of the fault. Interferometric synthetic aperture radar (InSAR) data from 13 independent ERS interferograms are stacked to obtain range change rates from 1992 to 2000. Horizontal surface displacement rates at 141 bench marks are measured using GPS from 1994 to 2003. Surface creep observations and estimates of deep slip rates determined from characteristic repeating earthquake sequences are also incorporated in the inversion. The fault is discretized into 283 triangular dislocation elements that approximate the nonplanar attributes of the fault surface. South of the city of Hayward, a steeply, east dipping fault geometry accommodates the divergence of the surface trace and the microseismicity at depth. The inferred slip rate distribution is consistent with a fault that creeps aseismically at a rate of ???5 mm/yr to a depth of 4-6 km. The interferometric synthetic aperture radar (InSAR) data require an aseismic slip rate that approaches the geologic slip rate on the northernmost fault segment beneath Point Pinole, although the InSAR data might be complicated by a small dip-slip component at this location. A low slip rate patch of <1 mm/yr is inferred beneath San Leandro consistent with the source location of the 1868 earthquake. We calculate that the entire fault is accumulating a slip rate deficit equivalent to a Mw = 6.77 ?? 0.05 per century. However, this estimate of potential coseismic moment represents an upper bound because we do not know how much of the accumulated strain will be released through aseismic processes such as afterslip. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

A model for seasonal changes in GPS positions and seismic wave speeds due to thermoelastic and hydrologic variations

It is known that GPS time series contain a seasonal variation that is not due to tectonic motions, and it has recently been shown that crustal seismic velocities may also vary seasonally. In order to explain these changes, a number of hypotheses have been given, among which thermoelastic and hydrology-induced stresses and strains are leading candidates. Unfortunately, though, since a general framework does not exist for understanding such seasonal variations, it is currently not possible to quickly evaluate the plausibility of these hypotheses. To fill this gap in the literature, I generalize a two-dimensional thermoelastic strain model to provide an analytic solution for the displacements and wave speed changes due to either thermoelastic stresses or hydrologic loading, which consists of poroelastic stresses and purely elastic stresses. The thermoelastic model assumes a periodic surface temperature, and the hydrologic models similarly assume a periodic near-surface water load. Since all three models are two-dimensional and periodic, they are expected to only approximate any realistic scenario; but the models nonetheless provide a quantitative framework for estimating the effects of thermoelastic and hydrologic variations. Quantitative comparison between the models and observations is further complicated by the large uncertainty in some of the relevant parameters. Despite this uncertainty, though, I find that maximum realistic thermoelastic effects are unlikely to explain a large fraction of the observed annual variation in a typical GPS displacement time series or of the observed annual variations in seismic wave speeds in southern California. Hydrologic loading, on the other hand, may be able to explain a larger fraction of both the annual variations in displacements and seismic wave speeds. Neither model is likely to explain all of the seismic wave speed variations inferred from observations. However, more definitive conclusions cannot be made until the model parameters are better constrained. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

USGS2012

Conductive heat flux in VC-1 and the thermal regime of Valles caldera, Jemez Mountains, New Mexico

Over 5% of heat in the western United States is lost through Quaternary silicic volcanic centers, including the Valles caldera in north central New Mexico. These centers are the sites of major hydrothermal activity and upper crustal metamorphism, metasomatism, and mineralization, producing associated geothermal resources. We present new heat flow data from Valles caldera core hole 1 (VC-1), drilled in the southwestern margin of the Valles caldera. Thermal conductivities were measured on 55 segments of core from VC-1, waxed and wrapped to preserve fluids. These values were combined with temperature gradient data to calculate heat flow. Above 335 m, which is probably unsaturated, heat flow is 247±16 mW m −2 . The only deep temperature information available is from an uncalibrated commercial log made 19 months after drilling. Gradients, derived from uncalibrated temperature logs, and conductivities are inversely correlated between 335 and 737 m, indicating a conductive thermal regime, and component heat fluxes over three depth intervals (335–539 m, 549–628 m, and 628–737 m) are in excellent agreement with each other with an average of 504±15 mW m −2 . Temperature logs to 518 m depth with well-calibrated temperature sensors result in a revised heat flow of 463±15 mW m. We use shallow thermal gradient data from 75 other sites in and around the caldera to interpret the thermal regime at the VC-1 site. A critical review of published thermal conductivity data from the Valles caldera yields an average thermal conductivity of ≥1 W m −1 K −1 for the near-surface tuffaceous material, and we assume that shallow gradient values (°C km −1 ) are approximately numerically equal to heat flow (mW m −2 ). Heat loss from the caldera is asymmetrically distributed, with higher values (400 mW m −2 or higher) concentrated in the west-southwestern quadrant of the caldera. This quadrant also contains the main drainage from the caldera and the youngest volcanism associated with the caldera. We interpret the shallow thermal gradient data and the thermal regime at VC-1 to indicate a long-lived hydrothermal (and magmatic) system in the southwestern Valles caldera that has been maintained through the generation of shallow magma bodies during the long postcollapse history of the caldera. High heat flow at the VC-1 site is interpreted to result from hot water circulating below the base of the core hole, and we attribute the lower heat flow in the unsaturated zone to hydrologic recharge.

Journal of Geophysical Research Solid Earth↗

USGS2000

Crustal deformation associated with glacial fluctuations in the eastern Chugach Mountains, Alaska

The changes of the solid Earth in south central Alaska in response to two major glacial fluctuations on different temporal and spatial scales have been estimated and we evaluated their influence on the stress state and ongoing tectonic deformation of the region. During the recent (1993–1995) Bering Glacier surge, a large transfer of ice from the Bagley Ice Field to the Bering Glacier terminus region occurred. We estimated the elastic displacement of the solid Earth due to ice mass redistribution from Global Positioning System (GPS) measurements at sites near the surging glacier. We can account for these displacements by transfer of an ice volume of about 14 km 3 from the surge reservoir area to the terminus region. We examined the background seismicity (M L ≥ 2.5) before, during, and after the surge. We found that the occurrence of small earthquakes (M L ≤ 4.0) in the surge reservoir region increased during the surge time interval possibly in response to a decrease in ice mass. This suggests that a small decrease in the vertical stress, σ 3 , could be enough to modulate the occurrence of small, shallow earthquakes in this dominantly thrust fault setting. During this century the southern Alaska coastal glaciers have been undergoing an overall decrease in volume. Based on our compilation of changes in the extent and thickness of the coastal glaciers between the Malaspina and Bering, we calculated surface displacements due to the Earth's viscoelastic response to annual thinning and to the cumulative retreat over the last 100 years. The uplift of the region due to an average annual thinning rate of 1–6 m/yr in the ablation region is 1–12 mm/yr. For our reference model with a viscosity of 5×10 19 Pa s for depths between ≈ 40 and 200 km the total viscoelastic response due to the retreat over the last century may be as much as a couple of meters within the coastal ablation zone near Icy Bay. The maximum decrease in σv between 0 and 10 km was ≈ 1.0 MPa, which is significant in relation to the stress drops in recent earthquakes (≈ 2 to 10 MPa) but small in relation to the estimated tectonic stress magnitude. Therefore the occurrence of an earthquake such as the St. Elias (1979, M S = 7.2) may have been advanced in time; however, most of the ongoing stress accumulation would be primarily due to tectonic forces.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.

Wyoming↗

USGS2015

Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California

Northwest directed slip from the southern San Andreas Fault is transferred to the Mission Creek, Banning, and Garnet Hill fault strands in the northwestern Coachella Valley. How slip is partitioned between these three faults is critical to southern California seismic hazard estimates but is poorly understood. In this paper, we report the first slip rate measured for the Banning fault strand. We constrain the depositional age of an alluvial fan offset 25 ± 5 m from its source by the Banning strand to between 5.1 ± 0.4 ka (95% confidence interval (CI)) and 6.4 + 3.7/−2.1 ka (95% CI) using U-series dating of pedogenic carbonate clast coatings and 10 Be cosmogenic nuclide exposure dating of surface clasts. We calculate a Holocene geologic slip rate for the Banning strand of 3.9 + 2.3/−1.6 mm/yr (median, 95% CI) to 4.9 + 1.0/−0.9 mm/yr (median, 95% CI). This rate represents only 25–35% of the total slip accommodated by this section of the southern San Andreas Fault, suggesting a model in which slip is less concentrated on the Banning strand than previously thought. In rejecting the possibility that the Banning strand is the dominant structure, our results highlight an even greater need for slip rate and paleoseismic measurements along faults in the northwestern Coachella Valley in order to test the validity of current earthquake hazard models. In addition, our comparison of ages measured with U-series and 10 Be exposure dating demonstrates the importance of using multiple geochronometers when estimating the depositional age of alluvial landforms.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Nucleation and growth of strike slip faults in granite

Fractures within granodiorite of the central Sierra Nevada, California, were studied to elucidate the mechanics of faulting in crystalline rocks, with emphasis on the nucleation of new fault surfaces and their subsequent propagation and growth. Within the study area the fractures form a single, subparallel array which strikes N50°–70°E and dips steeply to the south. Some of these fractures are identified as joints because displacements across the fracture surfaces exhibit dilation but no slip. The joints are filled with undeformed minerals, including epidote and chlorite. Other fractures are identified as small faults because they display left-lateral strike slip separations of up to 2 m. Slickensides, developed on fault surfaces, plunge 0°–20° to the east. The faults occur parallel to, and in the same outcrop with, the joints. The faults are filled with epidote, chlorite, and quartz, which exhibit textural evidence of shear deformation. These observations indicate that the strike slip faults nucleated on earlier formed, mineral-filled joints. Secondary, dilational fractures propagated from near the ends of some small faults contemporaneously with the left-lateral slip on the faults. These fractures trend 25°±10° from the fault planes, parallel to the direction of inferred local maximum compressive stress. The faults did not propagate into intact rock in their own planes as shear fractures. Rather, adjacent faults were linked together by secondary, dilational fractures. Extensive secondary fracturing between faults produced larger fault zones that accommodate 10–100 m of left-lateral slip. As deformation progressed, faulting evolved from relatively short, closely spaced faults to longer, more widely spaced fault zones.

Journal of Geophysical Research Solid Earth↗

USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗

USGS2009

Moon meteoritic seismic hum: Steady state prediction

We use three different statistical models describing the frequency of meteoroid impacts on Earth to estimate the seismic background noise due to impacts on the lunar surface. Because of diffraction, seismic events on the Moon are typically characterized by long codas, lasting 1 h or more. We find that the small but frequent impacts generate seismic signals whose codas overlap in time, resulting in a permanent seismic noise that we term the "lunar hum" by analogy with the Earth's continuous seismic background seismic hum. We find that the Apollo era impact detection rates and amplitudes are well explained by a model that parameterizes (1) the net seismic impulse due to the impactor and resulting ejecta and (2) the effects of diffraction and attenuation. The formulation permits the calculation of a composite waveform at any point on the Moon due to simulated impacts at any epicentral distance. The root-mean-square amplitude of this waveform yields a background noise level that is about 100 times lower than the resolution of the Apollo long-period seismometers. At 2 s periods, this noise level is more than 1000 times lower than the low noise model prediction for Earth's microseismic noise. Sufficiently sensitive seismometers will allow the future detection of several impacts per day at body wave frequencies. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

USGS2001

Vertical tectonics in northern Escanaba Trough as recorded by thick late Quaternary turbidites

Escanaba Trough, the southernmost segment of the Gorda Ridge, is filled by as much as 500 m of late Quaternary turbidite and hemipelagic sediment. Coring at Deep Sea Drilling Project Site 35 and Ocean Drilling Program (ODP) Sites 1037 and 1038 together with 4.5-kHz deep-tow and 3.5-kHz surface-ship seismic reflection profiles enable a distinct pattern of reflections to be mapped throughout Escanaba Trough in the upper part of this sediment fill. The uppermost 80 m of turbidite sediment, which includes at least 11 turbidity current events, were deposited in <7 kyr. Nine of these turbidites are found throughout Escanaba Trough at water depths >3200 m. The turbidity currents were trapped upon entering Escanaba Trough, resulting in all of the sediment in suspension in the flows being deposited. The thickness of the turbidite layers reflects both the flow thickness and the vertical grain concentration within the flow that deposited the layer. Variations in the turbidite thickness with respect to water depth can be used to estimate the degree of relative vertical movement within the floor of Escanaba Trough. In the area of hydrothermal activity near ODP Site 1038, uplift of as much as 140 m has occurred over the past 8 kyr.

California↗

USGS2012

Modeling crater topography and albedo from monoscopic Viking orbiter images: 1. Methodology

A new photoclinometric technique for extraction of topographic data from single planetary images is presented that overcomes many previous limitations of photoclinometry. The procedure fully compensates for oblique viewing geometry prevalent in spacecraft images. Albedo variations have been one of the most serious obstacles in the application of photoclinometry to planetary surfaces. This problem is overcome in the topographic solution by simultaneously utilizing brightness data from a pair of profiles; both segments are assumed to have the same topographic and albedo variations along their lengths. Profile directions are chosen where the orientation of downslope or upslope is obvious, thus resolving a major ambiguity in photoclinometry. This requirement is particularly easy to satisfy for craters and not very difficult for many irregular features. An additional procedure is presented that eliminates even the requirement of topographic symmetry along the pair of profiles. If two profiles have the same relief but their shapes are very different, another method can be used in an iterative process to derive topographic profiles; however, this procedure does, require that the albedo not vary along the profiles. Test results indicate that both procedures have an accuracy and precision of approximately 2° for slopes of typical bowl-shaped craters, which translates to approximately 5% for depths.

Journal of Geophysical Research Solid Earth↗

USGS2000

Three-dimensional seismic velocity structure of the San Francisco Bay area

Seismic travel times from the northern California earthquake catalogue and from the 1991 Bay Area Seismic Imaging Experiment (BASIX) refraction survey were used to obtain a three-dimensional model of the seismic velocity structure of the San Francisco Bay area. Nonlinear tomography was used to simultaneously invert for both velocity and hypocenters. The new hypocenter inversion algorithm uses finite difference travel times and is an extension of an existing velocity tomography algorithm. Numerous inversions were performed with different parameters to test the reliability of the resulting velocity model. Most hypocenters were relocated <2 km from their catalogue locations. Large lateral velocity variations at shallow (<4 km) depth correlate with known surface geology, including low-velocity Cenozoic sedimentary basins, high-velocity Cenozoic volcanic rocks, and outcrop patterns of the major Mesozoic geologic terranes. Salinian arc rocks have higher velocities than the Franciscan melange, which in turn are faster than Great Valley Sequence forearc rocks. The thickess of low-velocity sediment is defined, including >12 km under the Sacramento River Delta, 6 km beneath Livermore Valley, 5 km beneath the Santa Clara Valley, and 4 km beneath eastern San Pablo Bay. The Great Valley Sequence east of San Francisco Bay is 4-6 km thick. A relatively high velocity body exists in the upper 10 km beneath the Sonoma volcanic field, but no evidence for a large intrusion or magma chamber exists in the crust under The Geysers or the Clear Lake volcanic center. Lateral velocity contrasts indicate that the major strike-slip faults extend subvertically beneath their surface locations through most of the crust. Strong lateral velocity contrasts of 0.3-0.6 km/s are observed across the San Andreas Fault in the middle crust and across the Hayward, Rogers Creek, Calaveras, and Greenville Faults at shallow depth. Weaker velocity contrasts (0.1-0.3 km/s) exist across the San Andreas, Hayward, and Rogers Creek Faults at all other depths. Low spatial resolution evidence in the lower crust suggests that the top of high-velocity mafic rocks gets deeper from west to east and may be offset under the major faults. The data suggest that the major strike-slip faults extend subvertically through the middle and perhaps the lower crust and juxtapose differing lithology due to accumulated strike-slip motion. The extent and physical properties of the major geologic units as constrained by the model should be used to improve studies of seismicity, strong ground motion, and regional stress.

Journal of Geophysical Research B: Solid Earth↗

USGS2007

Generation of fine hydromagmatic ash by growth and disintegration of glassy rinds

The deposits of mafic hydromagmatic eruptions are more fine grained and variable in vesicularity than dry magmatic deposits. Blocky, equant shapes of many hydromagmatic clasts also contrast with droplet, thread, and bubble wall morphology of dry magmatic fragments. Small (<∼180 μ m), blocky hydromagmatic pyroclasts have traditionally been interpreted to result from discrete vapor explosions, although such explosions tend to occur only under certain conditions. This paper considers a process of hydromagmatic ash formation that involves repeated growth and disintegration of glassy rinds on pyroclast surfaces as they deform within turbulent flows. This process, termed “turbulent shedding”, may occur during the expansion phase of vapor explosions or during turbulent but nonexplosive mixing of magma with water, steam, or water sprays. The occurrence of turbulent shedding and the resulting fragment sizes depend on the timescale for rind growth and the timescale between disturbances that remove or disintegrate glassy rinds. Turbulent shedding is directly observable in some small littoral jets at Kilauea. Calculations suggest that, in the presence of liquid water or water sprays, glassy rinds having a thickness of microns to millimeters should form in milliseconds to seconds. This is similar to the timescale between turbulent velocity fluctuations that can shred lava globules and remove such rinds. The fraction of a deposit consisting of fine ash should increase with the duration of this process: Large‐scale Surtseyan jets generate hundreds or thousands of shedding events; bubble bursts or tephra jets at Kilauea's coast may produce only a few.

Journal of Geophysical Research B: Solid Earth↗

USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗

USGS2012

Analysis of ERS 1 synthetic aperture radar data of frozen lakes in northern Montana and implications for climate studies

Lakes that freeze each winter are good indicators of regional climate change if key parameters, such as freeze-up and breakup date and maximum ice thickness, are measured over a decade-scale time frame. Synthetic aperture radar (SAR) satellite data have proven to be especially useful for measurement of climatologically significant parameters characteristic of frozen lakes. In this paper, five lakes in Glacier National Park, Montana, have been studied both in the field and using Earth Remote-Sensing Satellite (ERS) 1 SAR data during the 1992-1993 winter. The lakes are characterized by clear ice, sometimes with tubular or rounded bubbles, and often with a layer of snow ice on top of the clear ice. They are also often snow covered. Freeze-up is detected quite easily using ERS 1 SAR data as soon as a thin layer of ice forms. The effect of snow ice on the backscatter is thought to be significant but is, as yet, undetermined. On the five lakes studied, relative backscatter was found to increase with ice thickness until a maximum was reached in February. Breakup, an often ill-defined occurrence, is difficult to detect because surface water causes the SAR signal to be absorbed, thus masking the ice below. Comparison of the bubble structure of thaw lakes in northern Alaska with lakes in northern Montana has shown that the ice structure is quite different, and this difference may contribute to differential SAR signature evolution in the lakes of the two areas.

Journal of Geophysical Research↗

USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗