Geology Reports⌕ Search

Geology topics

P.T. Delaney

Publications and source records attributed to P.T. Delaney.

17 recordsLinked to original sources
USGS2012

Rapid intrusion of magma into wet rock: Groundwater flow due to pore pressure increases

Analytical and numerical solutions are developed to simulate the pressurization, expansion, and flow of groundwater contained within saturated, intact host rocks subject to sudden heating from the planar surface of an igneous intrusion. For most rocks, water diffuses more rapidly than heat, assuring that groundwater is not heated along a constant-volume pressure path and that thermal expansion and pressurization adjacent to the intrusion drives a flow that extends well beyond the heated region. The forcing parameter for pressurization and flow is α Δ T , where α is a thermal expansion coefficient reflecting the overall expansion of water heated through the temperature difference Δ T between the initial ambient and intrusive values. Pore pressure increases due to heating are greatest when the intrusion is emplaced rapidly and where the intrusive contact is impervious to groundwater contained in stiff, impermeable rocks with high thermal diffusivities and porosities. The maximum velocity of water flowing in pores decays with the inverse square root of time and is insensitive to hydraulic properties of the host rocks. Pressures are lessened and flow directions are reversed with the onset of hydrothermal convection. This occurs at times ranging from hours to weeks after onset of intrusion. As magma rises into near-surface rocks, steam can be generated. Solutions indicate that pressure increases and velocities are sensitive to the overall amount of expansion rather than the behavior of the water-steam transition. Both the overall thermal expansion coefficient α and the temperature difference Δ T are greater in shallow (<1 km) environments than in deep (∼5 km) ones. Thus, for rocks with similar transport properties, pressure increases due to heating are greatest in shallow environments. Although solutions can be applied to rocks with a wide variety of properties, pressure increases are calculated for compliant quartz-rich sedimentary rocks with a porosities between 1 and 20% and permeabilities between 1 darcy and 1 μdarcy, subject to temperature increases of 500 and 1000 K at depths ranging from 0.1 to 5 km in a region of hydrostatic pressures and normal geothermal gradient. Such rocks, with porosities greater than 5%, permeabilities less than a 0.1 mdarcy, and drained hydrostatic compressibilities of 10 −4 /MPa, undergo pressure increases greater than 10 MPa (100 bars)for conditions typical of water table depths of 2.5 km and heating to 500 K above ambient. Similar rocks, but with permeabilities less than 1 mdarcy, undergo pressure increases of 10 MPa for conditions typical of 1 km water table depth. Rocks commonly considered to be good aquifers undergo pressure increases of less than 1 MPa, primarily because of their high permeability. Although these estimates neglect the effects of fracturing and brecciation that may accompany such pressure increases, calculations indicate that pressure increases due to heating of cool groundwater can lead to failure of host rocks by a phreatic mechanism.

Journal of Geophysical Research Solid Earth↗
USGS2000

Comment on 'Volume of magma accumulation or withdrawal estimated from surface uplift or subsidence, with application to the 1960 collapse of Kilauea volcano' by P.T. Delaney and D.F. McTigue

In volcanoes that store a significant quantity of magma within a subsurface summit reservoir, such as Kilauea, bulk compression of stored magma is an important mode of deformation. Accumulation of magma is also accompanied by crustal deformation, usually manifested at the surface as uplift. These two modes of deformation - bulk compression of resident magma and deformation of the volcanic edifice - act in concert to accommodate the volume of newly added magma. During deflation, the processes reverse and reservoir magma undergoes bulk decompression, the chamber contracts, and the ground surface subsides. Because magma compression plays a role in creating subsurface volume of accommodate magma, magma budget estimates that are derived from surface uplift observations without consideration of magma compression will underestimate actual magma volume changes.

Hawai'i↗
USGS2000

Spreading volcanoes

As volcanoes grow, they become ever heavier. Unlike mountains exhumed by erosion of rocks that generally were lithified at depth, volcanoes typically are built of poorly consolidated rocks that may be further weakened by hydrothermal alteration. The substrates upon which volcanoes rest, moreover, are often sediments lithified by no more than the weight of the volcanic overburden. It is not surprising, therefore, that volcanic deformation includes-and in the long term is often dominated by-spreading motions that translate subsidence near volcanic summits to outward horizontal displacements around the flanks and peripheries. We review examples of volcanic spreading and go on to derive approximate expressions for the time volcanoes require to deform by spreading on weak substrates. We also demonstrate that shear stresses that drive low-angle thrust faulting from beneath volcanic constructs have maxima at volcanic peripheries, just where such faults are seen to emerge. Finally, we establish a theoretical basis for experimentally derived scalings that delineate volcanoes that spread from those that do not.

Annual Review of Earth and Planetary Sciences↗
USGS1999

Stabilization of volcanic flanks by dike intrusion: An example from Kilauea

Dike propagation and dilation increases the compression of adjacent rocks. On volcanoes, especially oceanic shields, dikes are accordingly thought to be structurally destabilizing. As compression is incremented, volcanic flanks are driven outward or downslope and thus increase their susceptibility to destructive earthquakes and giant landslides. We show, however, that the 2-m-thick dike emplaced along the east rift zone of Kilauea in 1983 actually stabilized that volcano's flank. Specifically, production of flank earthquakes dropped more than twofold after 1983 as maximum downslope motion slowed to 6 cm ?? year-1 from approximately 40 cm ?? year-1 during 1980-1982. As much as 65 cm of deflationary subsidence above Kilauea's summit and upper rift zones accompanied the dike intrusion. According to recent estimates, this deflation corresponds to a reduction in magma-reservoir pressure of approximately 4 MPa, probably about as much as the driving pressure of the 1983 dike. The volume of the dike, approximately 0.10-0.15 km3, is orders of magnitude less than the estimated 200- to 250-km3 volume of Kilauea's reservoir of magma and nearby hot, mushy rock. Thus, deflation of that reservoir reduces the compressional load on the flank over a much larger area than intrusion of the dike adds to it, particularly at the dominant depth of seismicity, 8-9 km. A Coulomb block model for flank motion during intervals between major earthquakes requires the low-angle fault beneath Kilauea's flank to exhibit slip weakening, conducive to earthquake instability. Accordingly, the triggering mechanism of destructive earthquakes, several of which have struck Hawaii during the past 150 years, need not require stresses accumulated by dike intrusions.

Bulletin of Volcanology↗
USGS1998

Sources of crustal deformation associated with the Krafla, Iceland, eruption of September 1984

A decade-long plate-boundary rifting episode in northern Iceland ended with the September 1984 fissure eruption of Krafla volcano. We apply a nonlinear inversion method to geodetic data collected before and after the eruption to infer the location, geometry, and strengths of deformation sources associated with the eruption. The net outflow of magma from a 3-km-deep magma chamber beneath the Krafla caldera was 30−120× 10 6 m³. A similar volume of magma, 50−70×10 6 m³, was emplaced in a 1-meter-wide, ∼9-km-long dike extending from the surface to ∼7 km depth. Furthermore, at least 110×10 6 m³ of magma erupted. Accordingly, a surplus of magma must have been expelled from a second reservoir, the location of which, although uncertain, is likely to lie at depths greater than ∼5 km beneath Krafla volcano. It would be difficult to detect this deeper source because of the narrow aperture of the geodetic networks.

Geophysical Research Letters↗
USGS1997

Crustal structure, evolution, and volcanic unrest of the Alban Hills, Central Italy

The Alban Hills, a Quaternary volcanic center lying west of the central Apennines, 15-25 km southeast of Rome, last erupted 19ka and has produced approximately 290 km3 of eruptive deposits since the inception of volcanism at 580 ka. Earthquakes of moderate intensity have been generated there at least since the Roman age. Modern observations show that intermittent periods of swarm activity originate primarily beneath the youngest features, the phreatomagmatic craters on the west side of the volcano. Results from seismic tomography allow identification of a low-velocity region, perhaps still hot or partially molten, more than 6 km beneath the youngest craters and a high-velocity region, probably a solidified magma body, beneath the older central volcanic construct. Thirty centimeters of uplift measured by releveling supports the contention that high levels of seismicity during the 1980s and 1990s resulted from accumulation of magma beneath these craters. The volume of magma accumulation and the amount of maximum uplift was probably at least 40 ?? 106 m3 and 40 cm, respectively. Comparison of newer levelings with those completed in 1891 and 1927 suggests earlier episodes of uplift. The magma chamber beneath the western Alban Hills is probably responsible for much of the past 200 ka of eruptive activity, is still receiving intermittent batches of magma, and is, therefore, continuing to generate modest levels of volcanic unrest. Bending of overburden is the most likely cause of the persistent earthquakes, which generally have hypocenters above the 6-km-deep top of the magma reservoir. In this view, the most recent uplift and seismicity are probably characteristic and not precursors of more intense activity.

Alban Hills↗
USGS1997

Physical processes of shallow mafic dike emplacement near the San Rafael Swell, Utah

Some 200 shonkinite dikes, sills, and breccia bodies on the western Colorado Plateau of south-central Utah were intruded from approximately 3.7 to 4.6 Ma, contemporaneous with mafic volcanism along the nearby plateau margin. Thicknesses of dikes range to about 6 m; the log-normal mean thickness is 85 cm. Despite the excellent exposures of essentially all dikes in strata of the Jurassic San Rafael Group, their number is indeterminate from their outcrop and spacing because they are everywhere greatly segmented. By our grouping of almost 2000 dike segments, most dikes are less than 2 km in outcrop length; the longest is 9 km. Because the San Rafael magmas were primitive and probably ascended directly from the mantle, dike lengths in outcrop are much less than their heights. The present exposures probably lie along the irregular upper peripheries of dikes that lengthen and merge with depth. Orientations of steps on dike contacts record local directions of dike-fracture propagation; about half of the measurements plunge less than 30°, showing that lateral propagation at dike peripheries is as important as the vertical propagation ultimately responsible for ascent. The San Rafael dikes, now exposed after erosion of about 0.5–1.5 km, appear to thicken and shorten upward, probably because near-surface vesiculation enhanced magmatic driving pressures. Propagation likely ceased soon after the first dike segments began to feed nearby sills or vented to initiate small-volume eruptions.

Utah↗
USGS1994

Volume of magma accumulation or withdrawal estimated from surface uplift or subsidence, with application to the 1960 collapse of Kilauea volcano

An elastic point source model proposed by Mogi for magma chamber inflation and deflation has been applied to geodetic data collected at many volcanoes. The volume of ground surface uplift or subsidence estimated from this model is closely related to the volume of magma injection into or withdrawal from the reservoir below. The analytical expressions for these volumes are reviewed for a spherical chamber and it is shown that they differ by the factor 2(1-v), where v is Poisson's ratio of the host rock. For the common estimate v=0.25, as used by Mogi and subsequent workers, the uplift volume is 3/2 the injection volume. For highly fractured rocks, v can be even less and the uplift volume can approach twice the injection volume. Unfortunately, there is no single relation between the inflation of magma reservoirs and the dilation or contraction of host rocks. The inflation of sill-like bodies, for instance, generates no overall change in host rock volume. Inflation of dike-like bodies generates contraction such that, in contrast with Mogi's result, the uplift volume is generally less than the injection volume; for v=0.25, the former is only 3/4 of the latter. Estimates of volumes of magma injection or withdrawal are there-fore greatly dependent on the magma reservoir configuration. Ground surface tilt data collected during the 1960 collapse of Kilauea crater, one of the first events interpreted with Mogi's model and one of the largest collapses measured at Kilauea, is not favored by any one of a variety of deformation models. These models, however, predict substantially different volumes of both magma withdrawal and ground surface subsidence. ?? 1994 Springer-Verlag.

Bulletin of Volcanology↗
USGS1988

FORTRAN 77 programs for conductive cooling of dikes with temperature-dependent thermal properties and heat of crystallization

Temperature histories obtained from transient heat-conduction theory are applicable to most dikes despite potential complicating effects related to magma flow during emplacement, groundwater circulation, and metamorphic reaction during cooling. Here. machine-independent FORTRAN 77 programs are presented to calculate temperatures in and around dikes as they cool conductively. Analytical solutions can treat thermal-property contrasts between the dike and host rocks, but cannot address the release of magmatic heat of crystallization after the early stages of cooling or the appreciable temperature dependence of thermal conductivity and diffusivity displayed by most rock types. Numerical solutions can incorporate these additional factors. The heat of crystallization can raise the initial temperature at the dike contact, ??c1, about 100??C above that which would be estimated if it were neglected, and can decrease the rate at which the front of solidified magma moves to the dike center by a factor of as much as three. Thermal conductivity and diffusivity of rocks increase with decreasing temperature and, at low temperatures, these properties increase more if the rocks are saturated with water. Models that treat these temperature dependencies yield estimates of ??c1 that are as much as 75??C beneath those which would be predicted if they were neglected. ?? 1988.

Computers & Geosciences↗
USGS1986

Deformation of poorly consolidated sediment during shallow emplacement of a basalt sill, Coso Range, California

A 150-m-long, wedge-shaped unit of folded and faulted marly siltstone crops out between undeformed sedimentary rocks on the north flank of the Coso Range, California. The several-meter-thick blunt end of this wedge abuts the north margin of a basaltic sill of comparable thickness. Chaotically deformed siltstone crops out locally at the margin of this sill, and at one locality breccia pipes about one meter in diameter crosscut the sill. The sill extends about 1 km south up the paleoslope, where it merges through continuous outcrop with a lava flow that in turn extends 1.4 km to a vent area marked by more than 100 m of agglutinate and scoria. Apparently, lava extruded at this vent flowed onto unconsolidated sediments, burrowed into them, and fed a sill at about 40 m depth within the sedimentary sequence. The sill initially propagated by wedging between sedimentary beds, but eventually began to push some beds ahead of itself, forming a remarkable train of folds in the process. The sediments apparently were wet at the time of sill emplacement, because hydrothermal alteration is common near the contact between the two rock types and because the breccia pipes that crosscut the sill apparently resulted from phreatic explosions of pore water heated at the base of the cooling sill. Comparison of deformation of the host material at the Coso locality with that reportedly caused by emplacement of sills elsewhere indicates that the character of deformation differs greatly among the various localities. The specific response of host material depends upon such parameters as initial properties of magma and host material, rate of sill growth and attendant rate of strain of host material, and depth of sill emplacement. Some properties may change considerably during an intrusive-deformational episode, thus complicating accurate reconstruction of such an event.

Bulletin of Volcanology↗
USGS1982

Formation and interpretation of dilatant echelon cracks

The relative displacements of the walls of many veins, joints, and dikes demonstrate that these structures are dilatant cracks. We infer that dilatant cracks propagate in a principal stress plane, normal to the maximum tensile or least compressive stress. Arrays of echelon crack segments appear to emerge from the peripheries of some dilatant cracks. Breakdown of a parent crack into an echelon array may be initiated by a spatial or temporal rotation of the remote principal stresses about an axis parallel to the crack propagation direction. Near the parent-crack tip, a rotation of the local principal stresses is induced in the same sense, but not necessarily through the same angle. Incipient echelon cracks form at the parent-crack tip normal to the local maximum tensile stress. Further longitudinal growth along surfaces that twist about axes parallel to the propagation direction realigns each echelon crack into a remote principal stress plane. The walls of these twisted cracks may be idealized as helicoidal surfaces. An array of helicoidal cracks sweeps out less surface area than one parent crack twisting through the same angle. Thus, many echelon cracks grow from a single parent because the work done in creating the array, as measured by its surface area, decreases as the number of cracks increases. In cross sections perpendicular to the propagation direction, echelon cracks grow laterally, each crack overlapping its neighbors, until the mechanical interaction of adjacent cracks limits this growth. Dilation of each crack pinches the tips of adjacent cracks into an asymmetrical form and introduces local stresses that can cause lateral growth along a curving, sigmoidal path. Sigmoidal echelon cracks may link at tip-to-plane intersections, leaving a step in the through-going crack wall. The geometry of dilatant echelon cracks may be used to infer spatial or temporal changes in the orientation of principal stresses in the Earth.

Geological Society of America Bulletin↗