SEARCH · Search USGS
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
An analytical solution to the response of a long trapezoidal‐section dam on a foundation consisting of an elastic half‐space and subjected to simulated earthquake morion is developed. An optimum seismic design is achieved when the cross section of the dam is triangular. The effect of soil structure interaction is to lower the strain occurring in the dam.
A concise review of the fheory of adaptive modeling, identification, and control of dynamic structural systems based on discrete‐time recordings is presented. Adaptive methods have four major advantages over the classical methods: (1) Removal of the noise from the signal is done over the whole frequency band; (2) time‐varying characteristics of systems can be tracked; (3) systems with unknown characteristics can be controlled; and (4) a small segment of the data is needed during the computations. Included in the paper are the discrete‐time representation of single‐input single‐output (SISO) systems, models for SISO systems with noise, the concept of stochastic approximation, recursive prediction error method (RPEM) for system identification, and the adaptive control. Guidelines for model selection and model validation and the computational aspects of the method are also discussed in the paper. The present paper is the first of two companion papers. The theory given in the paper is limited to that which is necessary to follow the examples for applications in structural dynamics presented in the second paper.
This paper examines the rationale behind the semiempirical formulation of a generalized viscoplastic fluid (GVF) model in the light of the Reiner-Rivlin constitutive theory and the viscoplastic theory, thereby identifying the parameters that control the rheology of granular flow. The shear-rate number ( N ) proves to be among the most significant parameters identified from the GVF model. As N → 0 and N →∞, the GVF model can reduce asymptotically to the theoretical stress versus shear-rate relations in the macroviscous and grain-inertia regimes, respectively, where the grain concentration ( C ) also plays a major role in the rheology of granular flow. Using available data obtained from the rotating-cylinder experiments of neutrally buoyant solid spheres dispersing in an interstitial fluid, the shear stress for granular flow in transition between the two regimes proves dependent on N and C in addition to some material constants, such as the coefficient of restitution. The insufficiency of data on rotating-cylinder experiments cannot presently allow the GVF model to predict how a granular flow may behave in the entire range of N ; however, the analyzed data provide an insight on the interrelation among the relevant dimensionless parameters.
This paper assesses the validity of the generalized viscoplastic fluid (GVF) model in light of the established constitutive relations in two asymptotic flow regimes, namely, the macroviscous and grain-inertia regimes. A comprehensive review of the literature on constitutive relations in both regimes reveals that except for some material constants, such as the coefficient of restitution, the normalized shear stress in both regimes varies only with the grain concentration, C . It is found that Krieger-Dougherty's relative viscosity, μ * ( C ), is sufficiently coherent among the monotonically nondecreasing functions of C used in describing the variation of the shear stress with C in both regimes. It not only accurately represents the C -dependent relative viscosity of a suspension in the macroviscous regime, but also plays a role of the radial distribution function that describes the statistics of particle collisions in the grain-inertia regime. Use of μ * ( C ) alone, however, cannot link the two regimes. Another parameter, the shear-rate number, N , is needed in modeling the rheology of neutrally buoyant granular flows in transition between the two asymptotic regimes. The GVF model proves compatible with most established relations in both regimes.
This study examines the rationale of Hilbert-Huang transform (HHT) for analyzing dynamic and earthquake motion recordings in studies of seismology and engineering. In particular, this paper first provides the fundamentals of the HHT method, which consist of the empirical mode decomposition (EMD) and the Hilbert spectral analysis. It then uses the HHT to analyze recordings of hypothetical and real wave motion, the results of which are compared with the results obtained by the Fourier data processing technique. The analysis of the two recordings indicates that the HHT method is able to extract some motion characteristics useful in studies of seismology and engineering, which might not be exposed effectively and efficiently by Fourier data processing technique. Specifically, the study indicates that the decomposed components in EMD of HHT, namely, the intrinsic mode function (IMF) components, contain observable, physical information inherent to the original data. It also shows that the grouped IMF components, namely, the EMD-based low- and high-frequency components, can faithfully capture low-frequency pulse-like as well as high-frequency wave signals. Finally, the study illustrates that the HHT-based Hilbert spectra are able to reveal the temporal-frequency energy distribution for motion recordings precisely and clearly.
To experimentally validate a recently developed theory for predicting the stability of cut slopes under unsaturated conditions, the authors measured increasing strain localization in unsaturated slope cuts prior to abrupt failure. Cut slope width and moisture content were controlled and varied in a laboratory, and a sliding door that extended the height of the free face of the slope was lowered until the cut slope failed. A particle image velocimetry tool was used to quantify soil displacement in the x - y "> x - y (horizontal) and x - z "> x -z (vertical) planes, and strain was calculated from the displacement. Areas of maximum strain localization prior to failure were shown to coincide with the location of the eventual failure plane. Experimental failure heights agreed with the recently developed stability theory for unsaturated cut slopes (within 14.3% relative error) for a range of saturation and cut slope widths. A theoretical threshold for sidewall influence on cut slope failures was also proposed to quantify the relationship between normalized sidewall width and critical height. The proposed relationship was consistent with the cut slope experiment results, and is intended for consideration in future geotechnical experiment design. The experimental data of evolution of strain localization presented herein provide a physical basis from which future numerical models of strain localization can be validated.
For unsteady open-channel flows having N interacting unknown variables, a system of N mutually independent, partial differential equations can be used to describe the flow-field. The system generally belongs to marching-type problems and permits transformation into characteristic equations that are associated with N distinct characteristics directions. Because characteristics can be considered 'wave' or 'disturbance' propagation, a fluvial system so described can be viewed as adequately definable using these N component waves. A numerical algorithm to solve the N families of characteristics can then be introduced for formulation of an N-component flow-simulation model. The multimode method of characteristics (MMOC), a new numerical scheme that has a combined capacity of several specified-time-interval (STI) schemes of the method of characteristics, makes numerical modeling of such N-component riverine flows feasible and attainable. Merging different STI schemes yields different kinds of MMOC schemes, for which two kinds are displayed herein. With the MMOC, each characteristics is dynamically treated by an appropriate numerical mode, which should lead to an effective and suitable global simulation, covering various types of unsteady flow. The scheme is always linearly stable and its numerical accuracy can be systematically analyzed. By increasing the N value, one can develop a progressively sophisticated model that addresses increasingly complex river-mechanics problems.
Besides providing an exact solution for steady-state heat conduction processes (Laplace-Poisson equations), the CVBEM (complex variable boundary element method) can be used for the numerical error analysis of domain model solutions. For problems where soil-water phase change latent heat effects dominate the thermal regime, heat transport can be approximately modeled as a time-stepped steady-state condition in the thawed and frozen regions, respectively. The CVBEM provides an exact solution of the two-dimensional steady-state heat transport problem, and also provides the error in matching the prescribed boundary conditions by the development of a modeling error distribution or an approximate boundary generation. Consequently, this error evaluation can be used to develop highly accurate CVBEM models of the heat transport process, and the resulting model can be used as a test case for evaluating the precision of domain models based on finite elements or finite differences.
An infiltrometer was used at multiple locations at a site in Lee County, Fl. to define the spatial variability in infiltration parameters. Water-level data from a well at this site were collected hourly and used to determine the temporal variability in recharge. These results were used to define recharge in a representative stochastic numerical model of the aquifer. Model results without recharge compare well with existing analytical solutions for spatial head variability. Simulations with representative recharge events indicate that recharge produces a significant to dominant effect on head variability, which creates dispersion of contaminants, and that small-scale spatial and temporal recharge variations are the predominant mechanism causing the head variations.
A method of achieving a dimensionless collapse of erosion-rate data for cohesive sediments is proposed and shown to work well for data collected in flume-erosion tests on mixtures of sand and mud (silt plus clay sized particles) for a wide range of mud fraction. The data collapse corresponds to a dimensional erosion law of the form E ∼ ( τ − τ c ) m "> E ∼ ( τ − τ c ) m E∼(τ−τc)m , where E "> E E is erosion rate, τ "> τ τ is shear stress, τ c "> τ c τc is the threshold shear stress for erosion to occur, and m ≈ 7 / 4 "> m ≈ 7 / 4 m≈7/4 . This result contrasts with the commonly assumed linear erosion law E = k d ( τ − τ c ) "> E = k d ( τ − τ c ) E=kd(τ−τc) , where k d "> k d kd is a measure of how easily sediment is eroded. The data collapse prompts a re-examination of the way that results of the hole-erosion test (HET) and jet-erosion test (JET) are customarily analyzed, and also calls into question the meaningfulness not only of proposed empirical relationships between k d "> k d kd and τ c "> τ c τc , but also of the erodibility parameter k d "> k d kd itself. Fuller comparison of flume-erosion data with hole-erosion and jet-erosion data will require revised analyses of the HET and JET that drop the assumption m = 1 "> m = 1 m=1 and, in the case of the JET, certain simplifying assumptions about the mechanics of jet scour.
To better understand controls on debris-flow entrainment and travel distance, we examined topographic and drainage network characteristics of initiation locations in two separate debris-flow prone areas located 700 km apart along the west coast of the U.S. One area was located in northern California, the other in southern Oregon. In both areas, debris flows mobilized from slides during large storms, but, when stratified by number of contributing initiation locations, median debris-flow travel distances in Oregon were 5 to 8 times longer than median distances in California. Debris flows in Oregon readily entrained channel material; entrainment in California was minimal. To elucidate this difference, we registered initiation locations to high-resolution airborne LiDAR, and then examined travel distances with respect to values of slope, upslope contributing area, planform curvature, distance from initiation locations to the drainage network, and number of initiation areas that contributed to flows. Results show distinct differences in the topographic and drainage network characteristics of debris-flow initiation locations between the two study areas. Slope and planform curvature of initiation locations (landslide headscarps), commonly used to predict landslide-prone areas, were not useful for predicting debris-flow travel distances. However, a positive, power-law relation exists between median debris-flow travel distance and the number of contributing debris-flow initiation locations. Moreover, contributing area and the proximity of the initiation locations to the drainage network both influenced travel distances, but proximity to the drainage network was the better predictor of travel distance. In both study areas, flows that interacted with the drainage network flowed significantly farther than those that did not. In California, initiation sites within 60 m of the network were likely to reach the network and generate longtraveled flows; in Oregon, the threshold was 80 m.
The U.S. Geological Survey (USGS) principally has used Price AA and Price pygmy mechanical current meters for measurement of discharge. New technologies have resulted in the introduction of alternatives to the Price meters. One alternative, the FlowTracker acoustic Doppler velocimeter, was designed by SonTek/YSI to make streamflow measurements in wadeable conditions. The device measures a point velocity and can be used with standard midsection method algorithms to compute streamflow. The USGS collected 55 quality-assurance measurements with the FlowTracker at 43 different USGS streamflow-gaging stations across the United States, with mean depths from 0.05to0.67m, mean velocities from 13 to 60 cm/s, and discharges from 0.02 to 12.4m3/s. These measurements were compared with Price mechanical current meter measurements. Analysis of the comparisons shows that the FlowTracker discharges were not statistically different from the Price meter discharges at a 95% confidence level. ?? 2007 ASCE.
This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.
A mathematical model is developed and tested for the production of natural gas with water encroachment and gas entrapment. The model is built on the material and volumetric balance relations, the Schilthuis water drive model, and a gas entrapment mechanism which assumes that the rate of gas entrapment is proportional to the volumetric rate of water influx. This model represents an alternative to the large grid models because of its low computer, maintenance, and manpower costs.
A predictive, coastal erosion/shoreline change model has been developed for a small coastal segment near Drew Point, Beaufort Sea, Alaska. This coastal setting has experienced a dramatic increase in erosion since the early 2000’s. The bluffs at this site are 3-4 m tall and consist of ice-wedge bounded blocks of fine-grained sediments cemented by ice-rich permafrost and capped with a thin organic layer. The bluffs are typically fronted by a narrow ( ∼ 5 m wide) beach or none at all. During a storm surge, the sea contacts the base of the bluff and a niche is formed through thermal and mechanical erosion. The niche grows both vertically and laterally and eventually undermines the bluff, leading to block failure or collapse. The fallen block is then eroded both thermally and mechanically by waves and currents, which must occur before a new niche forming episode may begin. The erosion model explicitly accounts for and integrates a number of these processes including: (1) storm surge generation resulting from wind and atmospheric forcing, (2) erosional niche growth resulting from wave-induced turbulent heat transfer and sediment transport (using the Kobayashi niche erosion model), and (3) thermal and mechanical erosion of the fallen block. The model was calibrated with historic shoreline change data for one time period (1979-2002), and validated with a later time period (2002-2007).
A mass balance approach was used to determine the most important nonpoint source of volatile organic compounds (VOCs) in storm water from an asphalt parking lot without obvious point sources (e.g., gasoline stations). The parking lot surface and atmosphere are important nonpoint sources of VOCs, with each being important for different VOCs. The atmosphere is an important source of soluble, oxygenated VOCs (e.g., acetone), and the parking lot surface is an important source for the more hydrophobic VOCs (e.g., benzene). VOCs on the parking lot surface appear to be concentrated in oil and grease and organic material in urban particles (e.g., vehicle soot). Except in the case of spills, asphalt does not appear to be an important source of VOCs. The uptake isotherm of gaseous methyl tert -butyl ether on urban particles indicates a mechanism for dry deposition of VOCs from the atmosphere. This study demonstrated that a mass balance approach is a useful means of understanding non-point-source pollution, even for compounds such as VOCs, which are difficult to sample.
A 13-year study of tritium transport through a field-scale earthen liner was conducted by the Illinois State Geological Survey to determine the long-term performance of compacted soil liners in limiting chemical transport. Two field-sampling procedures (pressure-vacuum lysimeter and core sampling) were used to determine the vertical tritium concentration profiles at different times and locations within the liner. Profiles determined by the two methods were similar and consistent. Analyses of the concentration profiles showed that the tritium concentration was relatively uniformly distributed horizontally at each sampling depth within the liner and thus there was no apparent preferential transport. A simple one-dimensional analytical solution to the advective-dispersive solute transport equation was used to model tritium transport through the liner. Modeling results showed that diffusion was the dominant contaminant transport mechanism. The measured tritium concentration profiles were accurately modeled with an effective diffusion coefficient of 6 ?? 10-4 mm2/s, which is in the middle of the range of values reported in the literature.
The original paper discusses factors that may have contributed to the occurrence and long runout of a disastrous landslide near the community of Oso, Washington, on March 22, 2014. The paper reinforces a prior finding that the long runout likely resulted from liquefaction of wet colluvium that was rapidly loaded by landslide debris impinging from upslope ( Iverson et al. 2015 ). However, the original paper proposes a landslide geometry and chronology that differ from those presented in previous reports, and these differences have implications for interpreting the landslide’s mechanics. The paper also includes factual errors as well as misrepresentations that result from selective inclusion and exclusion of findings previously reported by others. The purposes of this discussion are to place the original paper in context and to rectify some of its mistakes.