SEARCH · Search USGS

Results for “Physics of Fluids”

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 181 records · Page 10Linked to original sources

USGS2018

Seismicity induced by wastewater injection in Washington County, Ohio: Influence of preexisting structure, regional stress regime, and well operations

Recent seismicity in Washington County, Ohio, has been suggested to be induced by wastewater disposal operations despite injection ~2 km above the Precambrian basement. We investigated the relationships between disposal well locations and operational histories, spatiotemporal patterns of seismicity enhanced by waveform correlation, and mapped subsurface structures. We also analyzed proxies for in situ stress conditions to evaluate existing principal stress orientations and magnitudes. Double‐difference relocations of small ( M L ≤ 2.1) seismicity revealed linear clusters defining a NE‐SW orientation that parallels the trend of small‐amplitude folds in Upper Paleozoic rocks within 1–2 km of the injection well. Hypocenter depths of ~3.1 to 4.3 km indicate seismicity occurred along faults below the Silurian injection interval in lower Paleozoic strata and Precambrian crystalline basement. Similar folds/faults imaged in nearby seismic reflection lines were associated with basement‐involved fault systems that cut the injection interval targeted by the injection well, providing permeability pathways for fluid pressure increases that initiated slip. Interpreted faults were also optimally oriented in the regional stress field for reactivation due to injection‐related decreases in effective stress. Well completion records suggest that during the recorded seismicity, the Long Run‐1 well was injecting disposal fluids close to or above the fracture pressure of the targeted Silurian injection interval. A comparison of observed seismicity rates and monthly injection records from the disposal well also indicate a close correspondence between earthquake frequency and injection rate. Taken together, these results provide a physical framework for understanding why wastewater disposal operations induced the Washington County seismic sequence.

Ohio↗

USGS1993

Geohydrologic systems in Kansas — Geohydrology of the lower aquifer unit in the western interior plains aquifer system

The purpose of the investigation is to provide a description of the principal geohydrologic systems in Upper Cambrian through Lower Cretaceous rocks in Kansas. This investigation was made as part of the Central Midwest Regional Aquifer-System Analysis (CMRASA). The CMRASA is one of several major investigations by the U.S. Geological Survey of regional aquifer systems in the United States. These regional investigations are designed to increase knowledge of the flow regime and hydrologic properties of major aquifer systems and to provide quantitative information for the assessment, development, and management of water supplies. The CMRASA study area includes all or parts of 10 Central Midwestern States (Jorgensen and Signor, 1981), as shown of the envelope cover. This Hydrologic Investigations Atlas, which consists of a series of chapters, presents a description of the physical framework and geohydrology of principal aquifers and confining systems in Kansas. Chapter H presents the geohydrology of the upper aquifer unit in the Western Interior Plains aquifer system. The physical framework of the aquifer system in relation to other systems is described by maps and sections showing areal extent and the thickness of rocks that compose the unit. The physical framework of the upper aquifer unit is described in detail in chapter D of the atlas (Hansen and others, in press). The hydrology of the system in relation to that of other systems is described in this chapter by maps showing the altitude of fluid levels and the direction of water movement within the unit. The chemical composition of water in the system is described by maps that show the distribution of dissolved-solids concentrations and the differences in water types on the basis of principal chemical constituents. Chapter A of this atlas series (Wolf and others, 1990) describes the relation of principal geohydrologic systems in Kansas and presents a more detailed discussion of the methods and data used to prepare and ensure consistency among the sets of maps.

Kansas↗

USGS2004

Diagenetic overprinting of the sphaerosiderite palaeoclimate proxy: are records of pedogenic groundwater δ 18 O values preserved?

Meteoric sphaerosiderite lines (MSLs), defined by invariant ??18O and variable ??13C values, are obtained from ancient wetland palaeosol sphaerosiderites (millimetre-scale FeCO3 nodules), and are a stable isotope proxy record of terrestrial meteoric isotopic compositions. The palaeoclimatic utility of sphaerosiderite has been well tested; however, diagenetically altered horizons that do not yield simple MSLs have been encountered. Well-preserved sphaerosiderites typically exhibit smooth exteriors, spherulitic crystalline microstructures and relatively pure (> 95 mol% FeCO3) compositions. Diagenetically altered sphaerosiderites typically exhibit corroded margins, replacement textures and increased crystal lattice substitution of Ca2+, Mg2+ and Mn2+ for Fe2+. Examples of diagenetically altered Cretaceous sphaerosiderite-bearing palaeosols from the Dakota Formation (Kansas), the Swan River Formation (Saskatchewan) and the Success S2 Formation (Saskatchewan) were examined in this study to determine the extent to which original, early diagenetic ??18O and ??13C values are preserved. All three units contain poikilotopic calcite cements with significantly different ??18O and ??13C values from the co-occurring sphaerosiderites. The complete isolation of all carbonate phases is necessary to ensure that inadvertent physical mixing does not affect the isotopic analyses. The Dakota and Swan River samples ultimately yield distinct MSLs for the sphaerosiderites, and MCLs (meteoric calcite lines) for the calcite cements. The Success S2 sample yields a covariant ??18O vs. ??13C trend resulting from precipitation in pore fluids that were mixtures between meteoric and modified marine phreatic waters. The calcite cements in the Success S2 Formation yield meteoric ??18O and ??13C values. A stable isotope mass balance model was used to produce hyperbolic fluid mixing trends between meteoric and modified marine end-member compositions. Modelled hyperbolic fluid mixing curves for the Success S2 Formation suggest precipitation from fluids that were < 25% sea water. ?? 2004 International Association of Sedimentologists.

Sedimentology↗

USGS2016

Wave attenuation in the shallows of San Francisco Bay

Waves propagating over broad, gently-sloped shallows decrease in height due to frictional dissipation at the bed. We quantified wave-height evolution across 7 km of mudflat in San Pablo Bay (northern San Francisco Bay), an environment where tidal mixing prevents the formation of fluid mud. Wave height was measured along a cross shore transect (elevation range−2mto+0.45mMLLW) in winter 2011 and summer 2012. Wave height decreased more than 50% across the transect. The exponential decay coefficient λ was inversely related to depth squared (λ=6×10 −4 h −2 ). The physical roughness length scale k b , estimated from near-bed turbulence measurements, was 3.5×10 −3 m in winter and 1.1×10 −2 m in summer. Estimated wave friction factor f w determined from wave-height data suggests that bottom friction dominates dissipation at high Re w but not at low Re w . Predictions of near-shore wave height based on offshore wave height and a rough formulation for f w were quite accurate, with errors about half as great as those based on the smooth formulation for f w . Researchers often assume that the wave boundary layer is smooth for settings with fine-grained sediments. At this site, use of a smooth f w results in an underestimate of wave shear stress by a factor of 2 for typical waves and as much as 5 for more energetic waves. It also inadequately captures the effectiveness of the mudflats in protecting the shoreline through wave attenuation.

California↗

USGS2007

BAT3 Analyzer: Real-time data display and interpretation software for the multifunction bedrock-aquifer transportable testing tool (BAT3)

The BAT3 Analyzer provides real-time display and interpretation of fluid pressure responses and flow rates measured during geochemical sampling, hydraulic testing, or tracer testing conducted with the Multifunction Bedrock-Aquifer Transportable Testing Tool (BAT3) (Shapiro, 2007). Real-time display of the data collected with the Multifunction BAT3 allows the user to ensure that the downhole apparatus is operating properly, and that test procedures can be modified to correct for unanticipated hydraulic responses during testing. The BAT3 Analyzer can apply calibrations to the pressure transducer and flow meter data to display physically meaningful values. Plots of the time-varying data can be formatted for a specified time interval, and either saved to files, or printed. Libraries of calibrations for the pressure transducers and flow meters can be created, updated and reloaded to facilitate the rapid set up of the software to display data collected during testing with the Multifunction BAT3. The BAT3 Analyzer also has the functionality to estimate calibrations for pressure transducers and flow meters using data collected with the Multifunction BAT3 in conjunction with corroborating check measurements. During testing with the Multifunction BAT3, and also after testing has been completed, hydraulic properties of the test interval can be estimated by comparing fluid pressure responses with model results; a variety of hydrogeologic conceptual models of the formation are available for interpreting fluid-withdrawal, fluid-injection, and slug tests.

Open-File Report↗

USGS2011

Direction of unsaturated flow in a homogeneous and isotropic hillslope

The distribution of soil moisture in a homogeneous and isotropic hillslope is a transient, variably saturated physical process controlled by rainfall characteristics, hillslope geometry, and the hydrological properties of the hillslope materials. The major driving mechanisms for moisture movement are gravity and gradients in matric potential. The latter is solely controlled by gradients of moisture content. In a homogeneous and isotropic saturated hillslope, absent a gradient in moisture content and under the driving force of gravity with a constant pressure boundary at the slope surface, flow is always in the lateral downslope direction, under either transient or steady state conditions. However, under variably saturated conditions, both gravity and moisture content gradients drive fluid motion, leading to complex flow patterns. In general, the flow field near the ground surface is variably saturated and transient, and the direction of flow could be laterally downslope, laterally upslope, or vertically downward. Previous work has suggested that prevailing rainfall conditions are sufficient to completely control these flow regimes. This work, however, shows that under time-varying rainfall conditions, vertical, downslope, and upslope lateral flow can concurrently occur at different depths and locations within the hillslope. More importantly, we show that the state of wetting or drying in a hillslope defines the temporal and spatial regimes of flow and when and where laterally downslope and/or laterally upslope flow occurs.

Water Resources Research↗

USGS1971

A paleohydrologic model for mineralization of the White Pine copper deposit, northern Michigan

Pertinent physical properties of the upper Keweenawan rocks can be measured or inferred within a sufficiently narrow range to make the quantitative evaluation of various paleohydrologic models for the origin of the White Pine copper deposit feasible. The approach is illustrated here by calculations for models that involve lateral migration of fluids through the subjacent Copper Harbor Conglomerate to the site of the deposit and stripping of copper from these solutions where they percolated upward through the Nonesuch Shale. The calculations reveal limitations to theories of origin that would not be evident from purely qualitative consideration; some of these limitations could be useful to exploration. For example, if the water was yielded by compaction of the Copper Harbor Conglomerate and contained 50 ppm Cu, there must have been significant convergence of solution paths toward White Pine. Surface water entering the Copper Harbor Conglomerate on the north limb of the Lake Superior syncline is an adequate source if it could be shown that the point of entry was once significantly higher in altitude than the water table at White Pine; this model implies a major copper deposit at great depth north of the axis of the syncline. Ground water entering the Nonesuch Shale up dip from White Pine is not a possible source of mineralizing solutions.

Michigan↗

USGS2018

CVPM 1.1: A flexible heat-transfer modeling system for permafrost

The Control Volume Permafrost Model (CVPM) is a modular heat-transfer modeling system designed for scientific and engineering studies in permafrost terrain, and as an educational tool. CVPM implements the nonlinear heat-transfer equations in 1-D, 2-D, and 3-D Cartesian coordinates, as well as in 1-D radial and 2-D cylindrical coordinates. To accommodate a diversity of geologic settings, a variety of materials can be specified within the model domain, including organic-rich materials, sedimentary rocks and soils, igneous and metamorphic rocks, ice bodies, borehole fluids, and other engineering materials. Porous materials are treated as a matrix of mineral and organic particles with pore spaces filled with liquid water, ice, and air. Liquid water concentrations at temperatures below 0 ∘ C due to interfacial, grain-boundary, and curvature effects are found using relationships from condensed matter physics; pressure and pore-water solute effects are included. A radiogenic heat-production term allows simulations to extend into deep permafrost and underlying bedrock. CVPM can be used over a broad range of depth, temperature, porosity, water saturation, and solute conditions on either the Earth or Mars. The model is suitable for applications at spatial scales ranging from centimeters to hundreds of kilometers and at timescales ranging from seconds to thousands of years. CVPM can act as a stand-alone model or the physics package of a geophysical inverse scheme, or serve as a component within a larger Earth modeling system that may include vegetation, surface water, snowpack, atmospheric, or other modules of varying complexity.

Geoscientific Model Development↗

USGS2014

Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model

Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

Hydrology and Earth System Sciences↗

USGS1990

Textural development of clayey and quartzofeldspathic fault gouges relative to their sliding behavior

Many of the secondary fault structures developed during triaxial friction experiments have been generally correlated with the structures of natural fault zones. Therefore, any physical differences that can be found between laboratory samples that slide stably and those that show stick-slip motion may help to identify the cause of earthquakes. We have examined petrographically the run products of many triaxial friction experiments using clayey and quartzofeldspathic gouges, which comprise the principal types of natural fault gouge material. The examined samples were tested under a wide range of temperature, confining and fluid pressure, and velocity conditions. The clayey and quartzofeldspathic gouges show some textural differences, owing to their different mineral contents and grain sizes and shapes. In the clayey gouges, for example, a clay mineral fabric and kink band sets are commonly developed, whereas in the quartzofeldspathic gouges fracturing and crushing of the predominately quartz and feldspar grains are important processes. For both types of gouge, however, and whatever the pressure-temperature-velocity conditions of the experiments, the transition from stable sliding to stick-slip motion is correlated with: (i) a change from pervasive deformation of the gouge layer to localized slip in subsidiary shears; and (ii) an increase in the angle betweem the shears that crosscut the gouge layer (Riedel shears) and ones that form along the gouge-rock cylinder boundaries (boundary shears). This suggests that the localization of shear within a fault zone combined with relatively high Riedel-shear angles are somehow connected with earthquakes. Secondary fracture sets similar to Riedel shears have been identified at various scales in major strike-slip faults such as the San Andreas of the western United States (Wallace, 1973) and the Luhuo and Fuyun earthquake faults of China (Deng and Zhang, 1984; Deng et al. , 1986). The San Andreas also contains locked and creeping sections that correspond to the stick-slip and stably sliding experimental samples, respectively. We plan to study the physical structure of the San Andreas fault, to see if the experimentally observed differences related to sliding behavior can also be distinguished in the field.

Physics and Chemistry of the Earth↗

USGS2014

Simultaneous Gaussian and exponential inversion for improved analysis of shales by NMR relaxometry

Nuclear magnetic resonance (NMR) relaxometry is commonly used to provide lithology-independent porosity and pore-size estimates for petroleum resource evaluation based on fluid-phase signals. However in shales, substantial hydrogen content is associated with solid and fluid signals and both may be detected. Depending on the motional regime, the signal from the solids may be best described using either exponential or Gaussian decay functions. When the inverse Laplace transform, the standard method for analysis of NMR relaxometry results, is applied to data containing Gaussian decays, this can lead to physically unrealistic responses such as signal or porosity overcall and relaxation times that are too short to be determined using the applied instrument settings. We apply a new simultaneous Gaussian-Exponential (SGE) inversion method to simulated data and measured results obtained on a variety of oil shale samples. The SGE inversion produces more physically realistic results than the inverse Laplace transform and displays more consistent relaxation behavior at high magnetic field strengths. Residuals for the SGE inversion are consistently lower than for the inverse Laplace method and signal overcall at short T2 times is mitigated. Beyond geological samples, the method can also be applied in other fields where the sample relaxation consists of both Gaussian and exponential decays, for example in material, medical and food sciences.

Journal of Magnetic Resonance↗

USGS1997

Review of electric and magnetic fields accompanying seismic and volcanic activity

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.

Surveys in Geophysics↗

USGS2007

Dynamic triggering

A growing body of evidence demonstrates that dynamic stresses propagating as seismic waves from large earthquakes are capable of triggering additional earthquakes ranging from aftershocks in the near-field (within one or two source dimensions of the mainshock epicenter) to remotely triggered earthquakes at distances exceeding 10 000 km. Most of the triggered earthquakes are small (generally M ≤ 3) except within the near field, where dynamic stresses may trigger slip on subadjacent fault segments leading to complex rupture comprising of several large earthquakes of comparable magnitude. Crustal surface waves with periods of 15–30 s and peak dynamic stresses greater than ∼0.01 MPa seem to be most efficient in triggering remote seismicity. Current models for dynamic triggering fall under two broad groups: one appealing to Coulomb failure with various friction laws, and the other appealing to the activation of crustal fluids either hydrous or magmatic. No single model appears capable of accounting for the wide variation observed in the nature of triggered activity. Spatial sampling of dynamic triggering on a global scale is still woefully inadequate because of the limited distribution of adequate seismic networks. From the limited data currently available, it appears that extensional stress regimes hosting geothermal and volcanic activity are more susceptible to remote dynamic triggering than compressional stress regimes, although remote triggering is not limited to extensional regimes. Instances of remote triggering in the few areas with continuous, high-resolution deformation instrumentation (all volcanic or geothermal areas) include distinctive deformation transients, suggesting that the locally triggered seismicity in these areas may be a secondary response to a more fundamental aseismic process that likely involves some form of fluid transport or phase change. Recent evidence for triggering by solid Earth tides and ocean loading in convergent plate margins provides a low-frequency, low-amplitude reference point for the spectrum of stresses capable of dynamic triggering. Remaining challenges include establishing better sampling of the distribution of triggered seismicity and better constraints on physical models for the triggering process.

Book chapter↗

USGS2014

Appendix C—Report on research in the field of ground water being conducted by oil companies

In view of the shortness of time since the appointment of the writer to the Committee on Ground Water this report is confined to the technology and problems in the Gulf Coast Oil Province. Of course, many of the methods and practices would apply to most parts of the country however, some would differ materially from one region to another. The writer wishes to acknowledge the suggestions and comments by F. H. LAHEE and PAUL WEAVER. Having been stationed in Houston, Texas, in the heart of the Gulf Coast Area for four and a half years, the writer has had an opportunity to view the great similarity of the problems confronting the petroleum geologist and engineer and the ground‐water hydrologist. Both groups deal with the accumulation, movement, and withdrawal of fluid from underground strata, yet each group is content to study its own literature and use its own terminology without much concern for the other. The petroleum and ground‐water engineer, independently of one another, have developed mathematical formulas for the determination of permeability from field‐data. These formulas use the same basic principles of physics and the initial papers on the subject by both groups were published within two years of one another. Because of the similarity in the technology and problems of the petroleum engineer and the ground‐water hydrologist there is a definite need for closer cooperation. Some of the problems are so closely related that their solution rests in cooperative studies.

Eos, Transactions, American Geophysical Union↗

USGS2016

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface↗

USGS2019

Valid debris-flow models must avoid hot starts

Debris-flow experiments and models commonly use “hot-start” initial conditions in which downslope motion begins when a large force imbalance is abruptly imposed. By contrast, initiation of natural debris flows almost invariably results from small perturbations of static force balances that apply to debris masses poised in steep channels or on steep slopes. Models that neglect these static balances may violate physical law. Here we assess how the effects of hot starts are manifested in physical experiments, analytical dam-break models, and numerical models in which frictional resistance is too small to satisfy static force balances in debris-flow source areas. We then outline a numerical modeling framework that avoids use of hot starts. In this framework an initial static force balance is gradually perturbed by increasing pore-fluid pressure that may trigger the onset of debris motion. Subsequent increases in pore-fluid pressure, driven by debris motion, may then reduce the debris frictional strength, leading to high flow mobility.

Conference Paper↗

USGS1986

Shallow subsurface temperature surveys in the Basin and Range province, U.S.A.-I. Review and evaluation

Temperature surveys at depths of 1–2 m have had varying success in geothermal exploration in the Basin and Range province. The most successful surveys have identified patterns of near-surface thermal-fluid flow within areas of less than 2 km 2 . Results have been less consistent in larger areas where zones of hydrothermal upflow are less well known, nongeothermal perturbing factors are significant and lateral variations in shallow subsurface temperature are small. Nongeothermal perturbations can be minimized by use of mean annual temperatures instead of synoptic temperatures, by physically based simulation of ground temperatures or by statistical modeling.

Geothermics↗

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↗