Search USGS⌕ Search

SEARCH · Search USGS

Results for “Acoustics”

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,045 records · Page 58Linked to original sources
USGS2019

Surface imaging functions for elastic reverse time migration

Reverse time migration is often used to interpret acoustic or three‐component seismic recordings by creating an image of subsurface seismic reflectors. Here I describe elastic reverse time migration imaging functions that are cast as waveform misfit sensitivity kernels of contrasts in material parameters across hypothetical seismic discontinuities, that is, specular reflectors. The proposed “surface” imaging functions are theoretically applicable to either reflected or converted waves in order to estimate the location and reflectivity of these discontinuities. The surface imaging functions, as well as volumetric sensitivity kernels that target point diffractors, are tested on sets of synthetic surface array recordings that sample the 3‐D seismic wavefield on simple 2‐D structures generated using a 2.5‐D spectral element method. These tests illustrate that in contrast with the volumetric sensitivity kernels, the reflectivity is generally dominated by a high‐amplitude peak that coincides with input locations of discontinuities. Passive recordings of microseismicity, shot gathers, or a combination thereof can be potentially interpreted with the new surface imaging functions to yield useful reflectivity images.

Journal of Geophysical Research Solid Earth↗
USGS2018

Submarine deposition of a subaerial landslide in Taan Fiord, Alaska

A large subaerial landslide entered Taan Fiord, Alaska, on 17 October 2015 producing a tsunami with runup to 193 m. We use LiDAR data to show the slide volume to be 76 + 3/−4 million cubic meters and that 51,000,000 m 3 entered Taan Fiord. In 2016, we mapped the fjord with multibeam bathymetry and high‐resolution seismic data. Landslide and postlandslide deposits extend 6 km downfjord, are up to 70 ± 11 m thick, and have a total volume of ~147,000,000 m 3 . Seismic data image a blocky landslide unit and two units deposited immediately after the landslide. The blocky landslide unit is ~65,000,000 m 3 . We infer it consists dominantly of subaerially derived material and secondarily of fjord floor sediment. The overlying units are likely megaturbidites presumably deposited within minutes to days after the landslide. We infer that these deposits dominantly consist of fjord floor material mobilized and suspended as the slide entered and traveled downfjord. The lower postlandslide unit is up to 35 ± 6 m thick, and the upper unit is up to 12 ± 3 m thick. These deposits are distinctive and will leave a lasting record of the event. This subaerial‐to‐submarine landslide deposit is distinct from other submarine landslide deposits studied in Alaskan fjords because it has a much greater thickness, larger and more angular blocks, distinctive postlandslide megaturbidites, and a higher‐amplitude acoustic signature of the blocky deposit. The tight constraints on the landslide source and deposit volumes, topography, bathymetry, and tsunami runup heights and flow directions should make this a benchmark site for landslide‐tsunami models.

Alaska↗
USGS2023

Dynamics of the December 2020 ash-poor plume formed by lava-water interaction at the summit of Kilauea Volcano, Hawaii

On 20 December 2020, after more than 2 years of quiescence at Kīlauea Volcano, Hawaiʻi, renewed volcanic activity in the summit crater caused boiling of the water lake over a period of ∼90 min. The resulting water-rich, electrified plume rose to 11–13 km above sea level, which is among the highest plumes on record for Kīlauea. Although conventional models would infer a high mass flux from explosive magma-water interaction, the plume was not associated with an infrasound signal indicative of “explosive” activity, nor did it produce a measurable ash-fall deposit. We use multisensor data to characterize lava-water interaction and plume generation during this opening phase of the 2020–21 eruption. Satellite, weather radar, and eyewitness observations revealed that the plume was rich in water vapor and hydrometeors but transported less ash than expected from its maximum height. Volcanic lightning flashes detected by ground-based cameras were confined to freezing altitudes of the upper cloud, suggesting that the ice formation drove the electrification of this plume. The low acoustic energy from lava-water interaction points to a weakly explosive style of hydrovolcanism. Heat transfer calculations show that the lava to water heat flux was sufficient to boil the lake within 90 min. Limited mixing of lava and water inhibited major steam explosions and fine fragmentation. Results from one-dimensional plume modeling suggest that the models may underpredict plume height due to overestimation of crosswind air-entrainment. Our findings shed light on an unusual style of volcanism in which weakly explosive lava-water interaction generated an outsized plume.

Hawaii↗
USGS2023

Time-lapse seafloor surveys reveal how turbidity currents and internal tides in Monterey Canyon interact with the seabed at centimeter-scale

Here we show how ultra-high resolution seabed mapping using new technology can help to understand processes that sculpt submarine canyons. Time-lapse seafloor surveys were conducted in the axis of Monterey Canyon, ∼50 km from the canyon head (∼1,840 m water depth) over an 18-month period. These surveys comprised 5-cm resolution multibeam bathymetry, 1-cm resolution lidar bathymetry, and 2-mm resolution stereophotographic imagery. Bathymetry data reveal centimeter-scale textures that would be undetectable by more traditional survey methods. Upward-looking Acoustic Doppler Current Profilers at the site recorded the flow character of internal tides and the passage of three turbidity currents, while sediment cores collected from the site record flow deposits. Combined with flow and core data, the bathymetry shows how turbidity currents and internal tides modify the seabed. The turbidity currents drape sediment across the site, infilling bedform troughs and smoothing erosional features carved by the internal tides (e.g., rippled scours). Turbidity currents with speeds of 0.9–3.3 m/s failed to cause notable bedform movement, which is surprising given that flows with similar speeds produced rapid bedform migration elsewhere, including the upper Monterey Canyon. The lack of migration may be related to the character of the underlying substrate or indicate that turbidity currents at the site lack dense, near-bed layers. The scale of scours produced by the internal tides (≤0.7 m/s) approaches the scale of features recorded in the ancient rock record. Thus, these results illustrate how the scale gap between seabed mapping technology and the rock record may eventually be bridged.

California↗
USGS2024

Laboratory hydrofractures as analogs to tectonic tremors

The fracture of Earth materials occurs over a wide range of time and length scales. Physical conditions, particularly the stress field and Earth material properties, may condition rupture in a specific fracture regime. In nature, fast and slow fractures occur concurrently: tectonic tremor events are fast enough to emit seismic waves and frequently accompany slow earthquakes, which are too slow to emit seismic waves and are referred to as aseismic slip events. In this study, we generate simultaneous seismic and aseismic processes in a laboratory setting by driving a penny-shaped crack in a transparent sample with pressurized fluid. We leverage synchronized high-speed imaging and high-frequency acoustic emission (AE) sensing to visualize and listen to the various sequences of propagation (breaks) and arrest (sticks) of a fracture undergoing stick-break instabilities. Slow radial crack propagation is facilitated by fast tangential fractures. Fluid viscosity and pressure regulate the fracture dynamics of slow and fast events, and control the inter-event time and the energy released during individual fast events. These AE signals share behaviors with observations of episodic tremors in Cascadia, United States; these include: (a) bursty or intermittent slow propagation, and (b) nearly linear scaling of radiated energy with area. Our laboratory experiments provide a plausible model of tectonic tremor as an indicative of hydraulic fracturing facilitating shear slip during slow earthquakes.

AGU Advances↗
USGS2023

UAS-Based Observations of Infrasound Directionality at Stromboli Volcano, Italy

Infrasound (low frequency sound waves) can be used to monitor and characterize volcanic eruptions. However, infrasound sensors are usually placed on the ground, thus providing a limited sampling of the acoustic radiation pattern that can bias source size estimates. We present observations of explosive eruptions from a novel uncrewed aircraft system (UAS)-based infrasound sensor platform that was strategically hovered near the active vents of Stromboli volcano, Italy. We captured eruption infrasound from short-duration explosions and jetting events. While potential vertical directionality was inconclusive for the short-duration explosion, we find that jetting events exhibit vertical sound directionality that was observed with a UAS close to vertical. This directionality would not have been observed using only traditional deployments of ground-based infrasound sensors, but is consistent with jet noise theory. This proof-of-concept study provides unique information that can improve our ability to characterize and quantify the directionality of volcanic eruptions and their associated hazards.

Stromboli Volcano↗
USGS2024

Siting considerations for satellite observation of river discharge

With growing global capability for satellite measurement of river discharge (flow) comes a need to understand and reduce error in satellite-based discharge measurements. Satellite-based discharge estimates are based on measurements of water surface width, elevation, slope, and potentially velocity. Site selection is important for reducing error and uncertainty in both conventional and satellite-based discharge measurements because geomorphic river characteristics have strong control over the relationships between discharge and width, water surface elevation (or depth), slope, and velocity. A large ground-truth data set of 8,445 conventional hydraulic measurements, collected by acoustic Doppler current profilers at 503 stations in the United States, was developed and quality assured to examine correlation between river discharge and water surface width, depth, velocity, and cross-sectional area. A separate database of river surface slope and discharge time-series was developed from paired continuous monitoring stations to examine slope-discharge correlations. Results show that discharge correlates most strongly with velocity, cross-sectional area, depth, width, and slope, in that order. Uncertainty of satellite discharge estimates is affected by observed hydraulic variable and reach-specific variability in observed variable(s) characteristics including range of variability, georegistration accuracy, and stability over time of relationships between discharge and observed hydraulic variable.

Water Resources Research↗
USGS2025

The first instrumentally detected hydrothermal explosion in Yellowstone National Park

Hydrothermal explosions are one of the geological hazards most likely to impact people in Yellowstone National Park, but their frequency is poorly known. Infrasound and seismic sensors identified an explosion in Norris Geyser Basin on 15 April 2024, at 14:56 MDT (20:56 UTC)—the first instrumentally detected hydrothermal explosion in the Yellowstone region. The event affected an area tens of meters across, resulting in fractured ground, a shallow explosion crater, and a field of ejecta. There were no immediate geophysical precursors, but in the preceding years elevated discharge of thermal water altered the color, temperature, and level of a nearby small lake. Expanded seismo-acoustic monitoring in Yellowstone National Park could be useful for detecting small hydrothermal explosions and constraining their frequency, magnitude, energy release, and locations—information that could be used to better assess and mitigate hazards for the millions of people that visit the park each year.

Wyoming↗
USGS2026

Interplay between tectonics and submarine mass transport deposits in Cortes Basin: New high-resolution geophysics in the Outer California Borderland

The Outer California Borderland (OCB) is an active transform plate boundary offshore Southern California, where the relationship between faulting and submarine mass transport deposits (MTDs) remains poorly understood. Onshore paleoseismic data provide high-resolution earthquake records, whereas marine geophysical data capture longer-term histories. Offshore fault systems pose hazards to infrastructure and dense coastal populations, particularly when linked to submarine landslides. We present new high-resolution geophysical data set (cruise SR2303), including bathymetric and CHIRP sub-bottom data integrated with legacy seismic reflection data and chronostratigraphic constraints from ODP Site 1012 to examine Quaternary MTD recurrence and tectonic controls in the Cortes Basin, OCB. Bathymetry shows deformational features, including slide scarps and previously unmapped fault segments with evidence of Holocene activity. CHIRP profiles reveal 10 stacked MTDs in the East Cortes Basin and 8 in the West Cortes Basin, spanning ∼752 ka with an average recurrence of ∼83.6 ± 1 ka. Acoustic imaging shows 7 MTD intervals coinciding with fault offset increments and fault growth suggesting earthquake-triggered mass wasting. A strong association between MTD occurrences and sea-level extremes also supports glacio-eustatic contribution to slope failure. Stratigraphic correlations suggest quasi-synchronous MTDs across the eastern and western areas, likely triggered by larger eathquakes in the Quaternary. Although the identified MTDs occur relatively far from the Southern California coast, they still pose a potential tsunamigenic hazard requiring further assessment. Moreover, if linked to earthquakes along major strike-slip faults, for example, the Ferrelo fault, the MTDs may provide valuable proxies to constrain rupture scenarios and fault connectivity within the understudied OCB.

California↗
USGS2026

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2012

Fractal properties of tremor and gas piston events observed at Kilauea Volcano, Hawaii

Studies the fractal properties of shallow volcanic tremor and gas piston events associated with magma degassing at Kilauea Volcano, Hawaii, using data from two dense short-baseline arrays of seismographs deployed near the active crater of Puu Oo on the east rift of the volcano. The existence of a categorically stable attractor characterizing both types of activities strongly suggests that the excitation mechanism of tremor is similar to that of gas piston events, which in turn are correlated with visual observations at the volcanic vent. Low values derived for the dimension of the attractor in phase space points to significant self-organization in the process of generation of tremor and offers general constraints on the dimensionality of attractors derived from models of acoustic emission associated with magma flow, vesiculation, and degassing. -from Authors

Hawaii↗
USGS2012

Modeling tidal exchange and dispersion in Boston Harbor

Tidal dispersion and the horizontal exchange of water between Boston Harbor and the surrounding ocean are examined with a high-resolution (200 m) depth-averaged numerical model. The strongly varying bathymetry and coastline geometry of the harbor generate complex spatial patterns in the modeled tidal currents which are verified by shipboard acoustic Doppler surveys. Lagrangian exchange experiments demonstrate that tidal currents rapidly exchange and mix material near the inlets of the harbor due to asymmetry in the ebb/flood response. This tidal mixing zone extends roughly a tidal excursion from the inlets and plays an important role in the overall flushing of the harbor. Because the tides can only efficiently mix material in this limited region, however, harbor flushing must be considered a two step process: rapid exchange in the tidal mixing zone, followed by flushing of the tidal mixing zone by nontidal residual currents. Estimates of embayment flushing based on tidal calculations alone therefore can significantly overestimate the flushing time that would be expected under typical environmental conditions. Particle-release simulations from point sources also demonstrate that while the tides efficiently exchange material in the vicinity of the inlets, the exact nature of dispersion from point sources is extremely sensitive to the timing and location of the release, and the distribution of particles is streaky and patchlike. This suggests that high-resolution modeling of dispersion from point sources in these regions must be performed explicitly and cannot be parameterized as a plume with Gaussian-spreading in a larger scale flow field.

Massachusetts↗
USGS1996

Analysis of fractures intersecting Kahi Puka Well 1 and its relation to the growth of the island of Hawaii

As part of the Hawaii Scientific Drilling Project, Kahi Puka Well 1 penetrated about 275 m of Mauna Loa basalts overlying a sequence of Mauna Kea flow units as it was drilled and cored to a total depth of 1053 m below land surface. A borehole televiewer (BHTV) was run in most of the well in successive stages prior to casing in order to obtain magnetically oriented acoustic images of the borehole wall. A total of 283 individual fractures were identified from this log and characterized in terms of strike and dip. These data are divided into three vertical sections based upon age and volcanic source, and lower hemisphere stereographic plots identify two predominant, subparallel fracture subsets common to each section. Assuming that most of the steeply dipping fractures observed in the BHTV log are tensile features generated within basalt flows during deposition and cooling, this fracture information can be combined with models of the evolution of the island of Hawaii to investigate the depositional history of these Mauna Loa and Mauna Kea basalts over the past 400 kyr. The directions of high-angle fractures appear to be generally parallel to topography or to the coastline at the time of deposition, as is supported by surface mapping of modern flows. Consequently, an overall counterclockwise rotation of about 75?? in the strike of these fractures from the bottom to the top of the well represents a systematic change in depositional slope direction over time. We attribute the observed rotation in the orientations of the two predominant fracture subsets over the past 400 kyr to changes in the configurations of volcanic sources during shield building and to the structural interference of adjacent volcanoes that produces shifts in topographic patterns.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Submarine landslides

Landslides are common on inclined areas of the seafloor, particularly in environments where weak geologic materials such as rapidly deposited, finegrained sediment or fractured rock are subjected to strong environmental stresses such as earthquakes, large storm waves, and high internal pore pressures. Submarine landslides can involve huge amounts of material and can move great distances: slide volumes as large as 20,000 km3 and runout distances in excess of 140 km have been reported. They occur at locations where the downslope component of stress exceeds the resisting stress, causing movement along one or several concave to planar rupture surfaces. Some recent slides that originated nearshore and retrogressed back across the shoreline were conspicuous by their direct impact on human life and activities. Most known slides, however, occurred far from land in prehistoric time and were discovered by noting distinct to subtle characteristics, such as headwall scarps and displaced sediment or rock masses, on acoustic-reflection profiles and side-scan sonar images. Submarine landslides can be analyzed using the same mechanics principles as are used for occurrences on land. However, some loading mechanisms are unique, for example, storm waves, and some, such as earthquakes, can have greater impact. The potential for limited-deformation landslides to transform into sediment flows that can travel exceedingly long distances is related to the density of the slope-forming material and the amount of shear strength that is lost when the slope fails.

Reviews of Geophysics↗
USGS1998

The excitation and characteristic frequency of the long-period volcanic event: An approach based on an inhomogeneous autoregressive model of a linear dynamic system

We present a method to quantify the source excitation function and characteristic frequencies of long-period volcanic events. The method is based on an inhomogeneous autoregressive (AR) model of a linear dynamic system, in which the excitation is assumed to be a time-localized function applied at the beginning of the event. The tail of an exponentially decaying harmonic waveform is used to determine the characteristic complex frequencies of the event by the Sompi method. The excitation function is then derived by operating an AR filter constructed from the characteristic frequencies to the entire seismogram of the event, including the inhomogeneous part of the signal. We apply this method to three long-period events at Kusatsu-Shirane Volcano, central Japan, whose waveforms display simple decaying monochromatic oscillations except for the beginning of the events. We recover time-localized excitation functions lasting roughly 1 s at the start of each event and find that the estimated functions are very similar to each other at all the stations of the seismic network for each event. The phases of the characteristic oscillations referred to the estimated excitation function fall within a narrow range for almost all the stations. These results strongly suggest that the excitation and mode of oscillation are both dominated by volumetric change components. Each excitation function starts with a pronounced dilatation consistent with a sudden deflation of the volumetric source which may be interpreted in terms of a choked-flow transport mechanism. The frequency and Q of the characteristic oscillation both display a temporal evolution from event to event. Assuming a crack filled with bubbly water as seismic source for these events, we apply the Van Wijngaarden-Papanicolaou model to estimate the acoustic properties of the bubbly liquid and find that the observed changes in the frequencies and Q are consistently explained by a temporal change in the radii of the bubbles characterizing the bubbly water in the crack.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Information on stress conditions in the oceanic crust from oval fractures in a deep borehole

Oval images etched into the vail of a deep borehole were detected in DSDP (Deep Sea Drilling Project) Hole 504B, eastern equatorial Pacific Ocean, from analysis of an acoustic televiever log. Initial speculation as to the origin of these remarkable features considered the intersection of pillow basalts with the wellbore. However, a systematic inspection of these ovals has identified intriguing consistencies in appearance that cannot be explained satisfactorily by a random, coincidental distribution of pillow lavas. As an alternative hypothesis, Mohr-Coulomb failure criterion is used to account for the generation and orientation of similarly curved, stress-induced fractures. Consequently, these oval features can be interpreted as fractures and related directly to stress conditions in the oceanic crust at this site. The azimuth of the oval center corresponds to the orientation of maximum horizontal principal stress (S H ), and the oval width, which spans approximately 180° of the borehole, is aligned vith the azimuth of minimum horizontal principal stress (S h ). The oval height is controlled by the fracture angle and thus is a function of the coefficient of internal friction of the rock.

Geophysical Research Letters↗
USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗
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↗