USGS Science⌕ Search

Geology topics

R.P. Denlinger

Publications and source records attributed to R.P. Denlinger.

16 recordsLinked to original sources
USGS2012

Geodetic analysis of reservoir depletion at the Geyser steam field in northern California

Reservoir depletion at the Geysers from 1974 to 1977 is evident in measured changes in gravity, surface strain, and pore pressure drainage. The drainage area increased about 20%, the maximum gravity decrease was about −120 μGal, and the maximum elevation change was about 6 cm during this period. Since the net mass withdrawal is known, it may be combined with the gravity change to estimate a drainage volume. The maximum drainage volume is 25 km 3 . Because the depth of the caprock and extent of the pore pressure drainage is known, this volume limit implies that no significant mass withdrawal occurs below a depth of 4 km. The ratios of surface elevation changes to horizontal contraction coinciding with the drainage area imply an equant drainage geometry, assuming that fluid production produces negative dilatation. Using the same cylindrical geometry used to model the gravity, negative dilatational strain rates of 4 to 5 × 10 −5 /yr were found to produce the observed surface displacements. The likelihood of boiling in the system to produce steam combined with apparent large bulk moduli (from seismic velocities) and small pore pressure declines suggest that most of the reservoir contraction is due to cooling.

Journal of Geophysical Research Solid Earth↗
USGS2012

Thermal stresses due to cooling of a viscoelastic oceanic lithosphere

Theories based upon thermal contraction of cooling oceanic lithosphere provide a successful basis for correlating seafloor bathymetry and heat flow. The horizontal components of the contraction of the lithosphere as it cools potentially give rise to large thermal stresses. Current methods to calculate these stresses assume that on the time scales of cooling, the lithosphere initially behaves as an inviscid fluid and instantly freezes into an elastic solid at some critical temperature. These instant-freezing methods inaccurately predict transient thermal stresses in rapidly cooling silicate glass plates because of the temperature dependent rheology of the material. The temperature dependent rheology of the lithosphere may affect the transient thermal stress distribution in a similar way, and for this reason we use a thermoviscoelastic model to estimate thermal stresses in young oceanic lithosphere. This theory is formulated here for linear creep processes that have an Arrhenius rate dependence on temperature. Our results show that the stress differences between instant freezing and linear thermoviscoelastic theory are most pronounced at early times (0–20 m.y.) when the instant freezing stresses may be twice as large. The solutions for the two methods asymptotically approach the same solution with time. A comparison with intraplate seismicity shows that both methods underestimate the depth of compressional stresses inferred from the seismicity in a systematic way.

Journal of Geophysical Research Solid Earth↗
USGS2012

Cooling rate and thermal structure determined from progressive magnetization of the dacite dome at Mount St. Helens, Washington

Our study of a magnetic anomaly associated with the recently active dacite dome at Mount St. Helens suggests that the dome consists of a hot, nonmagnetized core surrounded by a cool, magnetized carapace and flanking talus. The talus does not contribute to the anomaly because its constituent blocks are randomly oriented. Temporal changes in the magnetic anomaly indicate that the magnetized carapace thickened at an average rate of 0.03±0.01 m/d from 1984 to 1986. Petrographic and rock magnetic properties of dome samples indicate that the dominant process responsible for these changes is magnetization of extensively oxidized rock at progressively deeper levels within the dome as the rock cools through its blocking temperature, rather than subsequent changes in magnetization caused by further oxidation. Newly extruded material cools rapidly for a short period as heat is conducted outward in response to convective heat loss from its surface. The cooling rate gradually declines for several weeks, and thereafter the material cools at a relatively constant rate by convective heat loss from its interior along fractures that propagate inward. The rate of internal convective heat loss through fractures varies with rainfall, snowmelt, and large-scale fracturing during subsequent eruptive episodes. In accordance with a model for solidification of the 1959 lava lake at Kilauea Iki, Hawaii, we picture the dome's magnetized carapace as being a two-phase, porous, convective zone separated from the nonmagnetized core of the dome by a thin, single-phase conductive zone. As a consequence of the heat balance between the conductive and convective zones, the blocking-temperature isotherm migrates inward at a relatively constant rate. If the dome remains inactive, the time scale for its complete magnetization is estimated to be 18–36 years, a forecast which can be refined by shallow drilling into the dome and by continuing studies of its growing magnetic anomaly.

Journal of Geophysical Research Solid Earth↗
USGS2008

Computing nonhydrostatic shallow-water flow over steep terrain

Flood and dambreak hazards are not limited to moderate terrain, yet most shallow-water models assume that flow occurs over gentle slopes. Shallow-water flow over rugged or steep terrain often generates significant nonhydrostatic pressures, violating the assumption of hydrostatic pressure made in most shallow-water codes. In this paper, we adapt a previously published nonhydrostatic granular flow model to simulate shallow-water flow, and we solve conservation equations using a finite volume approach and an Harten, Lax, Van Leer, and Einfeldt approximate Riemann solver that is modified for a sloping bed and transient wetting and drying conditions. To simulate bed friction, we use the law of the wall. We test the model by comparison with an analytical solution and with results of experiments in flumes that have steep (31°) or shallow (0.3°) slopes. The law of the wall provides an accurate prediction of the effect of bed roughness on mean flow velocity over two orders of magnitude of bed roughness. Our nonhydrostatic, law-of-the-wall flow simulation accurately reproduces flume measurements of front propagation speed, flow depth, and bed-shear stress for conditions of large bed roughness.

Journal of Hydraulic Engineering↗
USGS2006

Usoi dam wave overtopping and flood routing in the Bartang and Panj Rivers, Tajikistan

The Usoi dam was created in the winter of 1911 after an enormous seismogenic rock slide completely blocked the valley of the Bartang River in the Pamir Mountains of southeastern Tajikistan. At present the dam impounds 17 million cubic meters of water in Lake Sarez. Flood volume and discharge estimates were made for several landslide generated floods that could overtop the dam. For landslide volumes of 200, 500, and 1,000 million cubic meters, estimated overtopping flood volumes were 2, 22, and 87 million cubic meters of water, respectively. Estimated peak discharge at the dam for these three flood scenarios were 57,000, 490,000, and 1,580,000 m3/s, based on triangular hydrographs of 70-, 90-, and 110-s durations, respectively. Flood-routing simulations were made for the three landslide-induced overtopping floods over a 530-km reach of the Bartang and Panj Rivers below the Usoi dam. A one-dimensional flow model using a Riemann numerical solution technique was selected for the analysis. For the 87 million cubic meter volume overtopping flood scenario, the peak flows were approximately 1, 100, 800, and 550 m3/s at locations 50, 100, and 150 km downstream of the dam respectively. ?? Springer 2006.

Natural Hazards↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain 2. Numerical predictions and experimental tests

Numerical solutions of the equations describing flow of variably fluidized Coulomb mixtures predict key features of dry granular avalanches and water-saturated debris flows measured in physical experiments. These features include time-dependent speeds, depths, and widths of flows as well as the geometry of resulting deposits. Threedimensional (3-D) boundary surfaces strongly influence flow dynamics because transverse shearing and cross-stream momentum transport occur where topography obstructs or redirects motion. Consequent energy dissipation can cause local deceleration and deposition, even on steep slopes. Velocities of surge fronts and other discontinuities that develop as flows cross 3-D terrain are predicted accurately by using a Riemann solution algorithm. The algorithm employs a gravity wave speed that accounts for different intensities of lateral stress transfer in regions of extending and compressing flow and in regions with different degrees of fluidization. Field observations and experiments indicate that flows in which fluid plays a significant role typically have high-friction margins with weaker interiors partly fluidized by pore pressure. Interaction of the strong perimeter and weak interior produces relatively steep-sided, flat-topped deposits. To simulate these effects, we compute pore pressure distributions using an advection-diffusion model with enhanced diffusivity near flow margins. Although challenges remain in evaluating pore pressure distributions in diverse geophysical flows, Riemann solutions of the depthaveraged 3-D Coulomb mixture equations provide a powerful tool for interpreting and predicting flow behavior. They provide a means of modeling debris flows, rock avalanches, pyroclastic flows, and related phenomena without invoking and calibrating Theological parameters that have questionable physical significance.

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

A dynamic balance between magma supply and eruption rate at Kilauea volcano, Hawaii

The dynamic balance between magma supply and vent output at Kilauea volcano is used to estimate both the volume of magma stored within Kilauea volcano and its magma supply rate. Throughout most of 1991 a linear decline in volume flux from the Kupaianaha vent on Kilauea's east rift zone was associated with a parabolic variation in the elevation of Kilauea's summit as vent output initially exceeded then lagged behind the magma supply to the volcano. The correspondence between summit elevation and tilt established with over 30 years of data provided daily estimates of summit elevation in terms of summit tilt. The minimum in the parabolic variation in summit tilt and elevation (or zero elevation change) occurs when the magma supply to the reservoir from below the volcano equals the magma output from the reservoir to the surface, so that the magma supply rate is given by vent flux on that day. The measurements of vent flux and tilt establish that the magma supply rate to Kilauea volcano on June 19, 1991, was 217,000±10,000 m 3 /d (or 0.079±0.004 km 3 /yr). This is close to the average eruptive rate of 0.08 km 3 /yr between 1958 and 1984. In addition, the predictable response of summit elevation and tilt to each east rift zone eruption near Puu Oo since 1983 shows that summit deformation is also a measure of magma reservoir pressure. Given this, the correlation between the elevation of the Puu Oo lava lake (4 km uprift of Kupaianaha and 18 km from the summit) and summit tilt provides an estimate for magma pressure changes corresponding to summit tilt changes. The ratio of the change in volume to the change in reservoir pressure ( dV / dP ) during vent activity may be determined by dividing the ratio of volume erupted to change in summit tilt ( dV / d tilt ) by the ratio of pressure change to change in summit tilt ( dP / d tilt ). This measure of dV / dP , when combined with laboratory measurements of the bulk modulus of tholeitic melt, provides an estimate of 240±50 km 3 for the volume of Kilauea's magma reservoir. This estimate is much larger than traditional estimates but consistent with seismic tomographic imaging and geophysical modeling of Kilauea's magma system.

Journal of Geophysical Research B: Solid Earth↗
USGS1994

Role of olivine cumulates in destabilizing the flanks of Hawaiian volcanoes

The south flank of Kilauea Volcano is unstable and has the structure of a huge landslide; it is one of at least 17 enormous catastrophic landslides shed from the Hawaiian Islands. Mechanisms previously proposed for movement of the south flank invoke slip of the volcanic pile over seafloor sediments. Slip on a low friction de??collement alone cannot explain why the thickest and widest sector of the flank moves more rapidly than the rest, or why this section contains a 300 km3 aseismic volume above the seismically defined de??collement. It is proposed that this aseismic volume, adjacent to the caldera in the direction of flank slip, consists of olivine cumulates that creep outward, pushing the south flank seawards. Average primary Kilauea tholeiitic magma contains about 16.5 wt.% MgO compared with an average 10 wt.% MgO for erupted subaerial and submarine basalts. This difference requires fractionation of 17 wt.% (14 vol.%) olivine phenocrysts that accumulate near the base of the magma reservoir where they form cumulates. Submarine-erupted Kilauea lavas contain abundant deformed olivine xenocrysts derived from these cumulates. Deformed dunite formed during the tholeiitic shield stage is also erupted as xenoliths in subsequent alkalic lavas. The deformation structures in olivine xenocrysts suggest that the cumulus olivine was densely packed, probably with as little as 5-10 vol.% intercumulus liquid, before entrainment of the xenocrysts. The olivine cumulates were at magmatic temperatures (>1100??C) when the xenocrysts were entrained. Olivine at 1100??C has a rheology similar to ice, and the olivine cumulates should flow down and away from the summit of the volcano. Flow of the olivine cumulates places constant pressure on the unbuttressed seaward flank, leading to an extensional region that localizes deep intrusions behind the flank; these intrusions add to the seaward push. This mechanism ties the source of gravitational instability to the caldera complex and deep rift systems and, therefore, limits catastrophic sector failure of Hawaiian volcanoes to their active growth phase, when the core of olivine cumulates is still hot enough to flow. ?? 1994 Springer-Verlag.

Bulletin of Volcanology↗
USGS1990

Limiting equilibrium and liquefaction potential in infinite submarine slopes

Stability evaluation of submarine slopes is hampered by the difficulty of making field measurements. Owing to the scarcity of detailed field data, stability is commonly assessed by assuming homogenous infinite slopes with steady seepage. For these conditions, it is necessary to measure only the slope angle, friction angle, cohesion, and pore pressure at some distance into the sediment to evaluate stability. Examination of available data shows that conditions close to those required for liquefaction are necessary for Coulomb failure in many continental shelf areas. This favors long landslide runouts and flow of sediment subsequent to failure. -from Authors

Marine Geotechnology↗