SEARCH · USGS Science

Results for “Energy”

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

USGS2003

Stratigraphic framework of sediment-starved sand ridges on a mixed siliciclastic/carbonate inner shelf; west-central Florida

Seismic reflection profiles and vibracores have revealed that an inner shelf, sand-ridge field has developed over the past few thousand years situated on an elevated, broad bedrock terrace. This terrace extends seaward of a major headland associated with the modern barrier-island coastline of west-central Florida. The overall geologic setting is a low-energy, sediment-starved, mixed siliciclastic/carbonate inner continental shelf supporting a thin sedimentary veneer. This veneer is arranged in a series of subparallel, shore-oblique, and to a minor extent, shore-parallel sand ridges. Seven major facies are present beneath the ridges, including a basal Neogene limestone gravel facies and a blue-green clay facies indicative of dominantly authigenic sedimentation. A major sequence boundary separates these older units from Holocene age, organic-rich mud facies (marsh), which grades upward into a muddy sand facies (lagoon or shallow open shelf/seagrass meadows). Cores reveal that the muddy shelf facies is either in sharp contact or grades upward into a shelly sand facies (ravinement or sudden termination of seagrass meadows). The shelly sand facies grades upward to a mixed siliciclastic/carbonate facies, which forms the sand ridges themselves. This mixed siliciclastic/carbonate facies differs from the sediment on the beach and shoreface, suggesting insignificant sediment exchange between the offshore ridges and the modern coastline. Additionally, the lack of early Holocene, pre-ridge facies in the troughs between the ridges suggests that the ridges themselves do not migrate laterally extensively. Radiocarbon dating has indicated that these sand ridges can form relatively quickly (???1.3 ka) on relatively low-energy inner shelves once open-marine conditions are available, and that frequent, high-energy, storm-dominated conditions are not necessarily required. We suggest that the two inner shelf depositional models presented (open-shelf vs. migrating barrier-island) may have co-existed spatially and/or temporally to explain the distribution of facies and vertical facies contacts. ?? 2003 Elsevier B.V. All rights reserved.

Marine Geology↗

USGS2002

Identification of lithofacies using Kohonen self-organizing maps

Lithofacies identification is a primary task in reservoir characterization. Traditional techniques of lithofacies identification from core data are costly, and it is difficult to extrapolate to non-cored wells. We present a low-cost automated technique using Kohonen self-organizing maps (SOMs) to identify systematically and objectively lithofacies from well log data. SOMs are unsupervised artificial neural networks that map the input space into clusters in a topological form whose organization is related to trends in the input data. A case study used five wells located in Appleton Field, Escambia County, Alabama (Smackover Formation, limestone and dolomite, Oxfordian, Jurassic). A five-input, one-dimensional output approach is employed, assuming the lithofacies are in ascending/descending order with respect to paleoenvironmental energy levels. To consider the possible appearance of new logfacies not seen in training mode, which may potentially appear in test wells, the maximum number of outputs is set to 20 instead of four, the designated number of lithosfacies in the study area. This study found eleven major clusters. The clusters were compared to depositional lithofacies identified by manual core examination. The clusters were ordered by the SOM in a pattern consistent with environmental gradients inferred from core examination: bind/boundstone, grainstone, packstone, and wackestone. This new approach predicted lithofacies identity from well log data with 78.8% accuracy which is more accurate than using a backpropagation neural network (57.3%). The clusters produced by the SOM are ordered with respect to paleoenvironmental energy levels. This energy-related clustering provides geologists and petroleum engineers with valuable geologic information about the logfacies and their interrelationships. This advantage is not obtained in backpropagation neural networks and adaptive resonance theory neural networks. ?? 2002 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗

USGS2013

Groundwater ages and mixing in the Piceance Basin natural gas province, Colorado

Reliably identifying the effects of energy development on groundwater quality can be difficult because baseline assessments of water quality completed before the onset of energy development are rare and because interactions between hydrocarbon reservoirs and aquifers can be complex, involving both natural and human processes. Groundwater age and mixing data can strengthen interpretations of monitoring data from those areas by providing better understanding of the groundwater flow systems. Chemical, isotopic, and age tracers were used to characterize groundwater ages and mixing with deeper saline water in three areas of the Piceance Basin natural gas province. The data revealed a complex array of groundwater ages (<10 to >50,000 years) and mixing patterns in the basin that helped explain concentrations and sources of methane in groundwater. Age and mixing data also can strengthen the design of monitoring programs by providing information on time scales at which water quality changes in aquifers might be expected to occur. This information could be used to establish maximum allowable distances of monitoring wells from energy development activity and the appropriate duration of monitoring.

Colorado↗

USGS2002

Ontogenetic behavior and migration of Chinese sturgeon, Acipenser sinensis

The Chinese sturgeon, Acipenser sinensis , is an anadromous protected species that presently only spawns in the Yangtze River. Using laboratory experiments, we examined the behavioral preference of young Chinese sturgeon to physical habitat (water depth, illumination intensity, substrate color, and cover) and monitored their downstream migration. Hatchling free embryos were photopositive, preferred open habitat, and immediately upon hatching, swam far above the bottom using swim-up and drift. Downstream migration peaked on days 0–1, decreased about 50% or more during days 2–7, and ceased by day 8. Days 0–1 migrants were active both day and night, but days 2–7 migrants were most active during the day. After ceasing migration, days 8–11 embryos were photonegative, preferred dark substrate and sought cover. Free embryos developed into larvae and began feeding on day 12, when another shift in behavior occurred–larvae returned to photopositive behavior and preferred white substrate. The selective factor favoring migration of free embryos upon hatching and swimming far above the bottom may be avoidance of benthic predatory fishes. Free embryos, which must rely on yolk energy for activity and growth, only used 19 cumulative temperature degree-days for peak migration compared to 234 degree-days for growth to first feeding larvae, a 1 : 12 ratio of cumulative temperature units. This ratio suggests that sturgeon species with large migratory embryos, like Chinese sturgeon, which require a high level of energy to swim during migration, may migrate only a short time to conserve most yolk energy for growth.

Jingzhou↗

USGS1999

Air blasts generated by rockfall impacts: Analysis of the 1996 Happy Isles event in Yosemite National Park

The July 10, 1996, Happy Isles rockfall in Yosemite National Park, California, released 23,000 to 38,000 m 3 of granite in four separate events. The impacts of the first two events which involved a 550-m free fall, generated seismic waves and atmospheric pressure waves (air blasts). We focus on the dynamic behavior of the second air blast that downed over 1000 trees, destroyed a bridge, demolished a snack bar, and caused one fatality and several injuries. Calculated velocities for the air blast from a two-phase, finite difference model are compared to velocities estimated from tree damage. From tornadic studies of tree damage, the air blast is estimated to have traveled <108–120 m/s within 50 m from the impact and decreased to <10–20 m/s within 500 m from the impact. The numerical model simulates the two-dimensional propagation of an air blast through a dusty atmosphere with initial conditions defined by the impact velocity and pressure. The impact velocity (105–107 m/s) is estimated from the Colorado Rockfall Simulation Program that simulates rockfall trajectories. The impact pressure (0.5 MPa) is constrained by the kinetic energy of the impact (10 10 –10 12 J) estimated from the seismic energy generated by the impact. Results from the air blast simulations indicate that the second Happy Isles air blast (weak shock wave) traveled with an initial velocity above the local sound speed. The size and location of the first impact are thought to have injected <50 wt% dust into the atmosphere. This amount of dust lowered the local atmospheric sound speed to ∼220 m/s. The discrepancy between calculated velocity data and field estimated velocity data (∼220 m/s versus ∼110 m/s) is attributed to energy dissipated by the downing of trees and additional entrainment of debris into the atmosphere not included in the calculations.

Journal of Geophysical Research B: Solid Earth↗

USGS2008

Constraining fault constitutive behavior with slip and stress heterogeneity

We study how enforcing self-consistency in the statistical properties of the preshear and postshear stress on a fault can be used to constrain fault constitutive behavior beyond that required to produce a desired spatial and temporal evolution of slip in a single event. We explore features of rupture dynamics that (1) lead to slip heterogeneity in earthquake ruptures and (2) maintain these conditions following rupture, so that the stress field is compatible with the generation of aftershocks and facilitates heterogeneous slip in subsequent events. Our three-dimensional fmite element simulations of magnitude 7 events on a vertical, planar strike-slip fault show that the conditions that lead to slip heterogeneity remain in place after large events when the dynamic stress drop (initial shear stress) and breakdown work (fracture energy) are spatially heterogeneous. In these models the breakdown work is on the order of MJ/m2, which is comparable to the radiated energy. These conditions producing slip heterogeneity also tend to produce narrower slip pulses independent of a slip rate dependence in the fault constitutive model. An alternative mechanism for generating these confined slip pulses appears to be fault constitutive models that have a stronger rate dependence, which also makes them difficult to implement in numerical models. We hypothesize that self-consistent ruptures could also be produced by very narrow slip pulses propagating in a self-sustaining heterogeneous stress field with breakdown work comparable to fracture energy estimates of kJ/M2. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2008

Constitutive relationships and physical basis of fault strength due to flash heating

We develop a model of fault strength loss resulting from phase change at asperity contacts due to flash heating that considers a distribution of contact sizes and nonsteady state evolution of fault strength with displacement. Laboratory faulting experiments conducted at high sliding velocities, which show dramatic strength reduction below the threshold for bulk melting, are well fit by the model. The predicted slip speed for the onset of weakening is in the range of 0.05 to 2 m/s, qualitatively consistent with the limited published observations. For this model, earthquake stress drops and effective shear fracture energy should be linearly pressure-dependent, whereas the onset speed may be pressure-independent or weakly pressure-dependent. On the basis of the theory, flash weakening is expected to produce large dynamic stress drops, small effective shear fracture energy, and undershoot. Estimates of the threshold slip speed, stress drop, and fracture energy are uncertain due to poor knowledge of the average ontact dimension, shear zone thickness and gouge particle size at seismogenic depths. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Composition and variation of noise recorded at the Yellowknife Seismic Array, 1991-2007

We analyze seismic noise recorded on the 18 short-period, vertical component seismometers of the Yellowknife Seismic Array (YKA). YKA has an aperture of 23 km and is sited on cratonic lithosphere in an area with low cultural noise. These properties make it ideal for studying natural seismic noise at periods of 1-3 s. We calculated frequency-wave number spectra in this band for over 6,000 time windows that were extracted once per day for 17 years (1991-2007). Slowness analysis reveals a rich variety of seismic phases originating from distinct source regions: R g waves from the Great Slave Lake; L g waves from the Atlantic, Pacific, and Arctic Oceans; and teleseismic P waves from the north Pacific and equatorial mid-Atlantic regions. The surface wave energy is generated along coastlines, while the body wave energy is generated at least in part in deep-water, pelagic regions. Surface waves tend to dominate at the longer periods and, just as in earthquake seismograms, L g is the most prominent arrival. Although the periods we study are slightly shorter than the classic double-frequency microseismic band of 4-10 s, the noise at YKA has clear seasonal behavior that is consistent with the ocean wave climate in the Northern Hemisphere. The temporal variation of most of the noise sources can be well fit using just two Fourier components: yearly and biyearly terms that combine to give a fast rise in microseismic power from mid-June through mid-October, followed by a gradual decline. The exception is the R g energy from the Great Slave Lake, which shows a sharp drop in noise power over a 2-week period in November as the lake freezes. The L g noise from the east has a small but statistically significant positive slope, perhaps implying increased ocean wave activity in the North Atlantic over the last 17 years. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

The impact of the 2009-10 El Niño Modoki on U.S. West Coast beaches

High-resolution beach morphology data collected along much of the U.S. West Coast are synthesized to evaluate the coastal impacts of the 2009–10 El Niño. Coastal change observations were collected as part of five beach monitoring programs that span between 5 and 13 years in duration. In California, regional wave and water level data show that the environmental forcing during the 2009–10 winter was similar to the last significant El Niño of 1997–98, producing the largest seasonal shoreline retreat and/or most landward shoreline position since monitoring began. In contrast, the 2009–10 El Niño did not produce anomalously high mean winter-wave energy in the Pacific Northwest (Oregon and Washington), although the highest 5% of the winter wave-energy measurements were comparable to 1997–98 and two significant non-El Niño winters. The increase in extreme waves in the 2009–10 winter was coupled with elevated water levels and a more southerly wave approach than the long-term mean, resulting in greater shoreline retreat than during 1997–98, including anomalously high shoreline retreat immediately north of jetties, tidal inlets, and rocky headlands. The morphodynamic response observed throughout the U.S. West Coast during the 2009–10 El Niño is principally linked to the El Niño Modoki phenomena, where the warm sea surface temperature (SST) anomaly is focused in the central equatorial Pacific (as opposed to the eastern Pacific during a classic El Niño), featuring a more temporally persistent SST anomaly that results in longer periods of elevated wave energy but lower coastal water levels.

Geophysical Research Letters↗

USGS2012

Uncovering the nonadiabatic response of geosynchronous electrons to geomagnetic disturbance

We describe an energy spectrum method for scaling electron integral flux, which is measured at a constant energy, to phase space density at a constant value of the first adiabatic invariant which removes much of the variation due to reversible adiabatic effects. Applying this method to nearly a solar cycle (1995 - 2006) of geosynchronous electron integral flux (E>2.0MeV) from the GOES satellites, we see that much of the diurnal variation in electron phase space density at constant energy can be removed by the transformation to phase space density at constant μ (4000 MeV/G). This allows us a clearer picture of underlying non-adiabatic electron population changes due to geomagnetic activity. Using scaled phase space density, we calculate the percentage of geomagnetic storms resulting in an increase, decrease or no change in geosynchronous electrons as 38%, 7%, and 55%, respectively. We also show examples of changes in the electron population that may be different than the unscaled fluxes alone suggest. These examples include sudden electron enhancements during storms which appear during the peak of negative Dst for μ-scaled phase space density, contrary to the slow increase seen during the recovery phase for unscaled phase space density for the same event.

Journal of Geophysical Research↗

USGS2020

Seasonal variation in sediment delivery across the bay-marsh interface of an estuarine salt marsh

Sediment transport across bay–marsh interfaces depends on wave energy, vegetation, and marsh-edge morphology, and varies over a range of timescales. We investigated these dynamics in a tidal salt marsh with a gently-sloped, vegetated edge adjacent to northern San Francisco Bay. Spartina foliosa (cordgrass) inhabits the lower marsh and Salicornia paciﬁca (pickleweed) predominates on the marsh plain. We measured suspended-sediment concentration (SSC) and hydrodynamics in bay shallows and along a 100-m cross-shore transect in the marsh, during winter and summer. Four-year averaged accretion measured with marker-horizon plots was twice as great along the marsh transect as adjacent to a tidal creek, 50 m from the bay. We estimated deposition and trapping eﬃciency from the time-series data to assess its variation with season and wave energy. At high tide the transition zone (between cordgrass and pickleweed) was usually erosional, the pickleweed zone was depositional, and both erosion and deposition increased with wave energy, as did the landward position of maximum deposition. Erosion from the transition zone accounted for approximately one-third of the sediment ﬂux into the pickleweed. In the pickleweed zone, SSC, the diﬀerence between ﬂood- and ebb-tide SSC and trapping eﬃciency were greater in summer than winter for comparable wave conditions, which we attribute to increased sediment trapping by dense summer cordgrass. Moderate waves in summer (46%) accounted for more annual accretion in the pickleweed zone than larger waves in winter (28%), although the contribution of winter storms was diminished by the dry winter during the study.

California↗

USGS2021

Using turbulence to identify preferential areas for grass carp (Ctenopharyngodon idella) larvae in streams: A laboratory study

In this experimental series, we studied the swimming capabilities and response of grass carp ( Ctenopharyngodon idella ) larvae to flow turbulence in a laboratory flume. We compared three different experimental configurations, representing in‐stream obstructions commonly found in natural streams (e.g., a gravel bump, a single vertical cylinder, and patches of submerged rigid vegetation). Grass carp larvae (postgas bladder emergence) were introduced to each experimental configuration and subjected to a variety of hydrodynamic forces of different magnitudes and scales. We varied the flow velocities and water depths and found ranges of turbulent kinetic energy and Reynolds stresses that triggered a response in larval trajectories, identified by measured horizontal and vertical swimming speeds for each flow condition. Larvae apparently actively avoided areas with increased levels of turbulence by swimming away, moving faster in short bursts, and expending more energy. In addition to the magnitude of turbulent kinetic energy, the length scale and time scale of turbulent eddies also influenced the larvae response. These findings support the development of new strategies for controlling the spread of grass carp larvae in rivers, as well as the development of numerical tools incorporating active swimming capabilities to predict larval transport in streams.

Water Resources Research↗

USGS2022

Physical controls on the hydrology of perennially ice-covered lakes, Taylor Valley, Antarctica (1996-2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. We apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996-2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2-17 m higher, depending on the lake, relative to the 2020 level.

JGR Earth Surface↗

USGS2023

Development and application of an Infragravity Wave (InWave) driver to simulate nearshore processes

Infragravity waves are key components of the hydro-sedimentary processes in coastal areas, especially during extreme storms. Accurate modeling of coastal erosion and breaching requires consideration of the effects of infragravity waves. Here, we present InWave, a new infragravity wave driver of the Coupled Ocean-Atmopshere-Waves-Sediment Transport (COAWST) modeling system. InWave computes the spatial and temporal variation of wave energy at the wave group scale and the associated incoming bound infragravity wave. Wave group-varying forces drive free infragravity wave growth and propagation within the hydrodynamic model of the coupled modeling system, which is the Regional Ocean Modeling System (ROMS) in this work. Since ROMS is a three-dimensional model, this coupling allows for the combined formation of undertow currents and infragravity waves. We verified the coupled InWave-ROMS with one idealized test case, one laboratory experiment, and one field experiment. The coupled modeling system correctly reproduced the propagation of gravity wave energy with acceptable numerical dissipation. It also captured the transfer of energy from the gravity band to the infragravity band, and within the different infragravity bands in the surf zone, the measured three-dimensional flow structure, and dune morphological evolution satisfactorily. The idealized case demonstrated that the infragravity wave variance depends on the directional resolution and horizontal grid resolution, which are known challenges with the approach taken here. The addition of InWave to COAWST enables novel investigation of nearshore hydro-sedimentary dynamics driven by infragravity waves using the strengths of the other modeling components, namely the three-dimensional nature of ROMS and the sediment transport routines.

Journal of Advances in Modeling Earth Systems↗

USGS2023

Subsurface porewater flow accelerates talik development under the Alaska Highway, Yukon: A prelude to road collapse and permafrost thaw?

The presence of taliks (perennially unfrozen zones in permafrost areas) adversely affects the thermal stability of infrastructure in cold regions, including roads. The role of heat advection on talik development and feedback on permafrost degradation has not been quantified methodically in this context. We incorporate a surface energy balance model into a coupled groundwater flow and energy transport numerical model (SUTRA-ice). The model, calibrated with long-term observations (1997–2018 on the Alaska Highway), is used to investigate and quantify the role of heat advection on talik initiation and development under a road embankment. Over the 25-year simulation period, the new model is driven by reconstructed meteorological data and has a good agreement with near surface soil temperatures. The model successfully reproduces the increasing depth to the permafrost table (mean absolute error <0.2 m), and talik development. The results demonstrate that heat advection provides an additional energy source that expedites the rate of permafrost thaw and roughly doubles the rate of permafrost table deepening, compared to purely conductive thawing. Talik initially formed and grew over time under the combined effect of water flow, snow insulation, road construction and climate warming. Talik formation creates a new thermal state under the road embankment, resulting in acceleration of underlying permafrost degradation, due to the positive feedback of heat accumulation created by trapped unfrozen water. In a changing climate, mobile water flow will play a more important role in permafrost thaw and talik development under road embankments, and is likely to significantly increase maintenance costs and reduce the long-term stability of the infrastructure.

Yukon↗

USGS2025

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗

USGS2010

Evaluation of 11 equations for determining evaporation for a small lake in the North Central United States

Eleven equations for calculating evaporation were compared with evaporation determined by the energy budget method for Williams Lake, Minnesota. Data were obtained from instruments on a raft, on land near the lake, and at a weather station 60 km south of the lake. The comparisons were based on monthly values for the open-water periods of 5 years, a total of 22 months. A modified DeBruin-Keijman, Priestley-Taylor, and a modified Penman equation resulted in monthly evaporation values that agreed most closely with energy budget values. To use these equations, net radiation, air temperature, wind speed, and relative humidity need to be measured near the lake. In addition, thermal surveys need to be made to determine change in heat stored in the lake. If data from distant climate stations are the only data available, and they include solar radiation, the Jensen-Haise and Makkink equations resulted in monthly evaporation values that agreed reasonably well with energy budget values.

Minnesota↗

USGS1997

The physics of debris flows

Recent advances in theory and experimentation motivate a thorough reassessment of the physics of debris flows. Analyses of flows of dry, granular solids and solid-fluid mixtures provide a foundation for a comprehensive debris flow theory, and experiments provide data that reveal the strengths and limitations of theoretical models. Both debris flow materials and dry granular materials can sustain shear stresses while remaining static; both can deform in a slow, tranquil mode characterized by enduring, frictional grain contacts; and both can flow in a more rapid, agitated mode characterized by brief, inelastic grain collisions. In debris flows, however, pore fluid that is highly viscous and nearly incompressible, composed of water with suspended silt and clay, can strongly mediate intergranular friction and collisions. Grain friction, grain collisions, and viscous fluid flow may transfer significant momentum simultaneously. Both the vibrational kinetic energy of solid grains (measured by a quantity termed the granular temperature) and the pressure of the intervening pore fluid facilitate motion of grains past one another, thereby enhancing debris flow mobility. Granular temperature arises from conversion of flow translational energy to grain vibrational energy, a process that depends on shear rates, grain properties, boundary conditions, and the ambient fluid viscosity and pressure. Pore fluid pressures that exceed static equilibrium pressures result from local or global debris contraction. Like larger, natural debris flows, experimental debris flows of ???10 m3 of poorly sorted, water-saturated sediment invariably move as an unsteady surge or series of surges. Measurements at the base of experimental flows show that coarse-grained surge fronts have little or no pore fluid pressure. In contrast, finer-grained, thoroughly saturated debris behind surge fronts is nearly liquefied by high pore pressure, which persists owing to the great compressibility and moderate permeability of the debris. Realistic models of debris flows therefore require equations that simulate inertial motion of surges in which high-resistance fronts dominated by solid forces impede the motion of low-resistance tails more strongly influenced by fluid forces. Furthermore, because debris flows characteristically originate as nearly rigid sediment masses, transform at least partly to liquefied flows, and then transform again to nearly rigid deposits, acceptable models must simulate an evolution of material behavior without invoking preternatural changes in material properties. A simple model that satisfies most of these criteria uses depth-averaged equations of motion patterned after those of the Savage-Hutter theory for gravity-driven flow of dry granular masses but generalized to include the effects of viscous pore fluid with varying pressure. These equations can describe a spectrum of debris flow behaviors intermediate between those of wet rock avalanches and sediment-laden water floods. With appropriate pore pressure distributions the equations yield numerical solutions that successfully predict unsteady, nonuniform motion of experimental debris flows.

Reviews of Geophysics↗