USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cells”

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.

At least 1,405 records · Page 78Linked to original sources
USGS2002

Preliminary cellular-automata forecast of permit activity from 1998 to 2010, Idaho and Western Montana

For public land management in Idaho and western Montana, the U.S. Forest Service (USFS) has requested that the U.S. Geological Survey (USGS) predict where mineral-related activity will occur in the next decade. Cellular automata provide an approach to simulation of this human activity. Cellular automata (CA) are defined by an array of cells , which evolve by a simple transition rule, the automaton . Based on exploration trends, we assume that future exploration will focus in areas of past exploration. Spatial-temporal information about mineral-related activity, that is permits issued by USFS and Bureau of Land Management (BLM) in the last decade, and spatial information about undiscovered resources, provide a basis to calibrate a CA. The CA implemented is a modified annealed voting rule that simulates mineral-related activity with spatial and temporal resolution of 1 mi 2 and 1 year based on activity from 1989 to 1998. For this CA, the state of the economy and exploration technology is assumed constant for the next decade. The calibrated CA reproduces the 1989–1998-permit activity with an agreement of 94%, which increases to 98% within one year. Analysis of the confusion matrix and kappa correlation statistics indicates that the CA underestimates high activity and overestimates low activity. Spatially, the major differences between the actual and calculated activity are that the calculated activity occurs in a slightly larger number of small patches and is slightly more uneven than the actual activity. Using the calibrated CA in a Monte Carlo simulation projecting from 1998 to 2010, an estimate of the probability of mineral activity shows high levels of activity in Boise, Caribou, Elmore, Lincoln, and western Valley counties in Idaho and Beaverhead, Madison, and Stillwater counties in Montana, and generally low activity elsewhere.

Idaho, Montana↗
USGS2004

Terrestrial microorganisms at an altitude of 20,000 m in Earth's atmosphere

A joint effort between the U.S. Geological Survey's (USGS) Global Desert Dust and NASA's Stratospheric and Cosmic Dust Programs identified culturable microbes from an air sample collected at an altitude of 20,000 m. A total of 4 fungal (Penicillium sp.) and 71 bacteria colonyforming units (70 colonies of Bacillus luciferensis believed to have originated from a single cell collected at altitude and one colony of Bacillus sphaericus) were enumerated, isolated and identified using a morphological key and 16S rDNA sequencing respectively. All of the isolates identified were sporeforming pigmented fungi or bacteria of terrestrial origin and demonstrate that the presence of viable microorganisms in Earth's upper atmosphere may not be uncommon.

Aerobiologia↗
USGS2004

Physiological changes in largemouth bass exposed to paper mill effluents under laboratory and field conditions

We report here on studies designed to asses the effects of paper mill effluents on non-reproductive functions of free-ranging and captive Florida largemouth bass (Micropterus salmoides floridanus) This was accomplished by conducting an outdoor tank study, in which fish were exposed to well water or to 10%, 20%, 40%, and 80% full strength effluent for 28 or 56 days, and by sampling largemouth bass from sites within the St. Johns River, Florida, upstream and downstream from a paper mill plant. Blood and plasma samples from fish from the tank study and from fish sampled from the ambient sites were analyzed for over 20 variables. We also determined liver and spleen weights and examined them histologically. The most significant finding from the tank study was an increase in the concentration of albumin and hepatosomatic index for bass exposed to ???20% effluents for 56 days. Spleenosomatic index and number of melanomacrophage centers were decreased in bass from effluent-dominated sites (Palatka and Rice Creek), whereas concentrations of calcium, phosphorous, glucose, and creatinine were elevated in fish from these sites, compared to fish from reference streams. Fish from Rice Creek also had fewer red blood cells, and male bass from Palatka had lower concentrations of cholesterol. Plasma concentrations of albumin and hepatic concentrations of glutathione were elevated in males from Palatka, and both females and males from Rice Creek had higher concentrations of globulin. These results indicate a complex pattern of effects of paper mill effluents on several physiological functions. However, despite the myriad of treatment and site-related effects, most physiological parameters fell within normal ranges when compared to reports on largemouth bass and other freshwater species.

Ecotoxicology↗
USGS2001

Statistical self-similarity of hotspot seamount volumes modeled as self-similar criticality

The processes responsible for hotspot seamount formation are complex, yet the cumulative frequency-volume distribution of hotspot seamounts in the Easter Island/Salas y Gomez Chain (ESC) is found to be well-described by an upper-truncated power law. We develop a model for hotspot seamount formation where uniform energy input produces events initiated on a self-similar distribution of critical cells. We call this model Self-Similar Criticality (SSC). By allowing the spatial distribution of magma migration to be self-similar, the SSC model recreates the observed ESC seamount volume distribution. The SSC model may have broad applicability to other natural systems.

Geophysical Research Letters↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 1. Theory and computation

To establish a theoretical basis for predicting and interpreting the behavior of rapid mass movements on Earth's surface, we develop and test a new computational model for gravity-driven motion of granular avalanches across irregular, three-dimensional (3-D) terrain. The principles embodied in the model are simple and few: continuum mass and momentum conservation and intergranular stress generation governed by Coulomb friction. However, significant challenges result from the necessity of satisfying these principles when deforming avalanches interact with steep and highly variable 3-D terrain. We address these challenges in four ways. (1) We formulate depth-averaged governing equations that are referenced to a rectangular Cartesian coordinate system (with z vertical) and that account explicitly for the effect of nonzero vertical accelerations on depth-averaged mass and momentum fluxes and stress states. (2) We compute fluxes of mass and momentum across vertical cell boundaries using a high-resolution finite volume method and Roe-type Riemann solver. Our algorithm incorporates flux difference splitting, an entropy correction for the flux, and eigenvector decomposition to embed the effects of driving and resisting forces in Riemann solutions. (3) We use a finite element method and avalanche displacements predicted by Riemann solutions to compute Coulomb stresses conjugate to the displacements in 3-D stress space. (4) We test the model output against analytical solutions, a sand cone conceptual experiment, and (in a companion paper) data from detailed laboratory experiments. Model results illustrate a complex interplay of basal traction and internal stress, and they successfully predict not only the gross behavior but also many details of avalanche motion from initiation to deposition.

Journal of Geophysical Research Earth Surface↗
USGS2005

Response of seismicity to Coulomb stress triggers and shadows of the 1999 Mw=7.6 Chi-Chi, Taiwan, earthquake

The correlation between static Coulomb stress increases and aftershocks has thus far provided the strongest evidence that stress changes promote seismicity, a correlation that the Chi-Chi earthquake well exhibits. Several studies have deepened the argument by resolving stress changes on aftershock focal mechanisms, which removes the assumption that the aftershocks are optimally oriented for failure. Here one compares the percentage of planes on which failure is promoted after the main shock relative to the percentage beforehand. For Chi-Chi we find a 28% increase for thrust and an 18% increase for strike-slip mechanisms, commensurate with increases reported for other large main shocks. However, perhaps the chief criticism of static stress triggering is the difficulty in observing predicted seismicity rate decreases in the stress shadows, or sites of Coulomb stress decrease. Detection of sustained drops in seismicity rate demands a long catalog with a low magnitude of completeness and a high seismicity rate, conditions that are met at Chi-Chi. We find four lobes with statistically significant seismicity rate declines of 40-90% for 50 months, and they coincide with the stress shadows calculated for strike-slip faults, the dominant faulting mechanism. The rate drops are evident in uniform cell calculations, 100-month time series, and by visual inspection of the M ??? 3 seismicity. An additional reason why detection of such declines has proven so rare emerges from this study: there is a widespread increase in seismicity rate during the first 3 months after Chi-Chi, and perhaps many other main shocks, that might be associated with a different mechanism. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Forecasting the evolution of seismicity in southern California: Animations built on earthquake stress transfer

We develop a forecast model to reproduce the distibution of main shocks, aftershocks and surrounding seismicity observed during 1986-200 in a 300 ?? 310 km area centered on the 1992 M = 7.3 Landers earthquake. To parse the catalog into frames with equal numbers of aftershocks, we animate seismicity in log time increments that lengthen after each main shock; this reveals aftershock zone migration, expansion, and densification. We implement a rate/state algorithm that incorporates the static stress transferred by each M ??? 6 shock and then evolves. Coulomb stress changes amplify the background seismicity, so small stress changes produce large changes in seismicity rate in areas of high background seismicity. Similarly, seismicity rate declines in the stress shadows are evident only in areas with previously high seismicity rates. Thus a key constituent of the model is the background seismicity rate, which we smooth from 1981 to 1986 seismicity. The mean correlation coefficient between observed and predicted M ??? 1.4 shocks (the minimum magnitude of completeness) is 0.52 for 1986-2003 and 0.63 for 1992-2003; a control standard aftershock model yields 0.54 and 0.52 for the same periods. Four M ??? 6.0 shocks struck during the test period; three are located at sites where the expected seismicity rate falls above the 92 percentile, and one is located above the 75 percentile. The model thus reproduces much, but certainly not all, of the observed spatial and temporal seismicity, from which we infer that the decaying effect of stress transferred by successive main shocks influences seismicity for decades. Finally, we offer a M ??? 5 earthquake forecast for 2005-2015, assigning probabilities to 324 10 ?? 10 km cells.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Role of microbial iron reduction in the dissolution of iron hydroxysulfate minerals

Iron-hydroxysulfate minerals can be important hosts for metals such as lead, mercury, copper, zinc, silver, chromium, arsenic, and selenium and for radionuclides such as 226Ra. These mineral-bound contaminants are considered immobilized under oxic conditions. However, when anoxic conditions develop, the activities of sulfate- or iron-reducing bacteria could result in mineral dissolution, releasing these bound contaminants. Reduction of structural sulfate in the iron-hydroxysulfate mineral jarosite by sulfate-reducing bacteria has previously been demonstrated. The primary objective of this work was to evaluate the potential for anaerobic dissolution of the iron-hydroxysulfate minerals jarosite and schwertmannite at neutral PH by iron-reducing bacteria. Mineral dissolution was tested using a long-term cultivar, Geobacter metallireducens strain GS-15, and a fresh isolate Geobacter sp. strain ENN1, previously undescribed. ENN1 was isolated from the discharge site of Shadle Mine, in the southern anthracite coalfield of Pennsylvania, where schwertmannite was the predominant iron-hydroxysulfate mineral. When jarosite from Elizabeth Mine (Vermont) was provided as the sole terminal electron acceptor, resting cells of both G. metallireducens and ENN1 were able to reduce structural Fe(III), releasing Fe+2, SO4-2, and K+ ions. A lithified jarosite sample from Utah was more resistant to microbial attack, but slow release of Fe+2 was observed. Neither bacterium released Fe+2 from poorly crystalline synthetic schwertmannite. Our results indicate that exposure of jarosite to iron-reducing conditions at neutral pH is likely to promote the mobility of hazardous constituents and should therefore be considered in evaluating waste disposal and/or reclamation options involving jarosite-bearing materials.

Journal of Geophysical Research G: Biogeosciences↗
USGS2006

Vital effects in coral skeletal composition display strict three-dimensional control

Biological control over coral skeletal composition is poorly understood but critically important to paleoenvironmental reconstructions. We present microanalytical measurements of trace-element abundances as well as oxygen and carbon isotopic compositions of individual skeletal components in the zooxanthellate coral Colpophyllia sp. Our data show that centers of calcification (COC) have higher trace element concentrations and distinctly lighter isotopic compositions than the fibrous components of the skeleton. These observations necessitate that COC and the fibrous skeleton are precipitated by different mechanisms, which are controlled by specialized domains of the calicoblastic cell-layer. Biological processes control the composition of the skeleton even at the ultra-structure level. Copyright 2006 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2007

Airborne dust transport to the eastern Pacific Ocean off southern California: Evidence from San Clemente Island

Islands are natural dust traps, and San Clemente Island, California, is a good example. Soils on marine terraces cut into Miocene andesite on this island are clay-rich Vertisols or Alfisols with vertic properties. These soils are overlain by silt-rich mantles, 5-20 cm thick, that contrast sharply with the underlying clay-rich subsoils. The silt mantles have a mineralogy that is distinct from the island bedrock. Silt mantles are rich in quartz, which is rare in the island andesite. The clay fraction of the silt mantles is dominated by mica, also absent from local andesite, and contrasts with the subsoils, dominated by smectite. Ternary plots of immobile trace elements (Sc-Th-La and Ta-Nd-Cr) show that the island andesite has a composition intermediate between average upper continental crust and average oceanic crust. In contrast, the silt and, to a lesser extent, clay fractions of the silt mantles have compositions closer to average upper continental crust. The silt mantles have particle size distributions similar to loess and Mojave Desert dust, but are coarser than long-range-transported Asian dust. We infer from these observations that the silt mantles are derived from airborne dust from the North American mainland, probably river valleys in the coastal mountains of southern California and/or the Mojave Desert. Although average winds are from the northwest in coastal California, easterly winds occur numerous times of the year when "Santa Ana" conditions prevail, caused by a high-pressure cell centered over the Great Basin. Examination of satellite imagery shows that easterly Santa Ana winds carry abundant dust to the eastern Pacific Ocean and the California Channel Islands. Airborne dust from mainland North America may be an important component of the offshore sediment budget in the easternmost Pacific Ocean, a finding of potential biogeochemical and climatic significance.

Journal of Geophysical Research D: Atmospheres↗
USGS2008

In situ Raman spectroscopic investigation of the structure of subduction-zone fluids

In situ Raman spectra of synthetic subduction-zone fluids (KAlSi3O8-H2O system) were measured to 900?? and 2.3 GPa using a hydrothermal diamond-anvil cell. The structures of aqueous fluid and hydrous melt become closer when conditions approach the second critical endpoint. Almost no three-dimensional network was observed in the supercritical fluid above 2 GPa although a large amount of silicate component is dissolved, suggesting that the physical and chemical properties of these phases change drastically at around the second critical endpoint. Our experimental results indicate that the fluids released from a subducting slab change from aqueous fluid to supercritical fluid with increasing depth under the volcanic arcs. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Diffuse flow hydrothermal manganese mineralization along the active Mariana and southern Izu-Bonin arc system, western Pacific

Abundant ferromanganese oxides were collected along 1200 km of the active Izu-Bonin-Mariana arc system. Chemical compositions and mineralogy show that samples were collected from two deposit types: Fe-Mn crusts of mixed hydrogenetic/hydrothermal origin and hydrothermal Mn oxide deposits; this paper addresses only the second type. Mn oxides cement volcaniclastic and biogenic sandstone and breccia layers (Mn sandstone) and form discrete dense stratabound layers along bedding planes and within beds (stratabound Mn). The Mn oxide was deposited within coarse-grained sediments from diffuse flow systems where precipitation occurred below the seafloor. Deposits were exposed at the seabed by faulting, mass wasting, and erosion. Scanning electron microscopy and microprobe analyses indicate the presence of both amorphous and crystalline 10 ?? and 7 ?? manganate minerals, the fundamental chemical difference being high water contents in the amorphous Mn oxides. Alternation of amorphous and crystalline laminae occurs in many samples, which likely resulted from initial rapid precipitation of amorphous Mn oxides from waxing pulses of hydrothermal fluids followed by precipitation of slow forming crystallites during waning stages. The chemical composition is characteristic of a hydrothermal origin including strong fractionation between Fe (mean 0.9 wt %) and Mn (mean 48 wt %) for the stratabound Mn, generally low trace metal contents, and very low rare earth element and platinum group element contents. However, Mo, Cd, Zn, Cu, Ni, and Co occur in high concentrations in some samples and may be good indicator elements for proximity to the heat source or to massive sulfide deposits. For the Mn sandstones, Fe (mean-8.4%) and Mn (12.4%) are not significantly fractionated because of high Fe contents in the volcaniclastic material. However, the proportion of hydrothermal Fe (nondetrital Fe) to total Fe is remarkably constant (49-58%) for all the sample groups, regardless of the degree of Mn mineralization. Factor analyses indicate various mixtures of two dominant components: hydrothermal Mn oxide for the stratabound Mn and detrital aluminosilicate for the Mn-cemented sandstone; and two minor components, hydrothermal Fe oxyhydroxide and biocarbonate/biosilica. Our conceptual model shows that Mn mineralization was produced by hydrothermal convection cells within arc volcanoes and sedimentary prisms that occur along, the flanks and within calderas. The main source of hydrothermal fluid was seawater that penetrated through fractures, faults, and permeable volcanic edifices. The fluids were heated by magma, enriched in metals by leaching of basement rocks and sediments, and mixed with magmatic fluids and gases. Dikes and sills may have been another source of heat that drove small-scale circulation within sedimentary prisms. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Pathogen and chemical transport in the karst limestone of the Biscayne aquifer: 3. Use of microspheres to estimate the transport potential of Cryptosporidium parvum oocysts

The vulnerability of a municipal well in the Northwest well field in southeastern Florida to potential contamination by Cryptosporidium parvum oocysts was assessed in a large‐scale, forced‐gradient (convergent) injection and recovery test. The field study involved a simultaneous pulse introduction of a nonreactive tracer (SF 6 , an inert gas) and oocyst‐sized (1.6, 2.9, and 4.9 μ m diameter) carboxylated polystyrene microspheres into karst limestone of the Biscayne aquifer characterized by a complex triple (matrix, touching‐vug, and conduit) porosity. Fractional recoveries 97 m down gradient were inversely related to diameter and ranged from 2.9% for the 4.9 μ m microspheres to 5.8% for 1.6 μ m microspheres. Their centers of mass arrived at the pumping well approximately threefold earlier than that of the nonreactive tracer SF 6 (gas), underscoring the need for use of colloid tracers and field‐scale tracer tests for these kinds of evaluations. In a modified triaxial cell using near in situ chemical conditions, 2.9 and 4.9 μ m microspheres underestimated by fourfold to sixfold the attachment potential of the less electronegative 2.9–4.1 μ m oocysts in the matrix porosity of limestone core samples. The field and laboratory results collectively suggested that it may take 200–300 m of transport to ensure even a 1‐log unit removal of oocysts, even though the limestone surfaces exhibited a substantive capability for their sorptive removal. The study further demonstrated the utility of microspheres as oocyst surrogates in field‐scale assessments of well vulnerability in limestone, provided that differences in attachment behaviors between oocysts and microspheres are taken into account.

Water Resources Research↗
USGS2008

Effects of turbulence on hydraulic heads and parameter sensitivities in preferential groundwater flow layers

A conduit flow process (CFP) for the Modular Finite Difference Ground‐Water Flow model, MODFLOW‐2005, has been created by the U.S. Geological Survey. An application of the CFP on a carbonate aquifer in southern Florida is described; this application examines (1) the potential for turbulent groundwater flow and (2) the effects of turbulent flow on hydraulic heads and parameter sensitivities. Turbulent flow components were spatially extensive in preferential groundwater flow layers, with horizontal hydraulic conductivities of about 5,000,000 m d −1 , mean void diameters equal to about 3.5 cm, groundwater temperature equal to about 25°C, and critical Reynolds numbers less than or equal to 400. Turbulence either increased or decreased simulated heads from their laminar elevations. Specifically, head differences from laminar elevations ranged from about −18 to +27 cm and were explained by the magnitude of net flow to the finite difference model cell. Turbulence also affected the sensitivities of model parameters. Specifically, the composite‐scaled sensitivities of horizontal hydraulic conductivities decreased by as much as 70% when turbulence was essentially removed. These hydraulic head and sensitivity differences due to turbulent groundwater flow highlight potential errors in models based on the equivalent porous media assumption, which assumes laminar flow in uniformly distributed void spaces.

Water Resources Research↗
USGS2010

Parametric study of the physical properties of hydrate-bearing sand, silt, and clay sediments: 1. Electromagnetic properties

The marked decrease in bulk electrical conductivity of sediments in the presence of gas hydrates has been used to interpret borehole electrical resistivity logs and, to a lesser extent, the results of controlled source electromagnetic surveys to constrain the spatial distribution and predicted concentration of gas hydrate in natural settings. Until now, an exhaustive laboratory data set that could be used to assess the impact of gas hydrate on the electromagnetic properties of different soils (sand, silt, and clay) at different effective stress and with different saturations of hydrate has been lacking. The laboratory results reported here are obtained using a standard geotechnical cell and the hydrate-formed tetrahydrofuran (THF), a liquid that is fully miscible in water and able to produce closely controlled saturations of hydrate from dissolved phase. Both permittivity and electrical conductivity are good indicators of the volume fraction of free water in the sediment, which is in turn dependent on hydrate saturation. Permittivity in the microwave frequency range is particularly predictive of free water content since it is barely affected by ionic concentration, pore structure, and surface conduction. Electrical conductivity (or resistivity) is less reliable for constraining water content or hydrate saturation: In addition to fluid-filled porosity, other factors, such as the ionic concentration of the pore fluid and possibly other conduction effects (e.g., surface conduction in high specific surface soils having low conductivity pore fluid), also influence electrical conductivity.

Journal of Geophysical Research↗
USGS2010

Parametric study of the physical properties of hydrate‐bearing sand, silt, and clay sediments: 2. Small‐strain mechanical properties

The small‐strain mechanical properties (e.g., seismic velocities) of hydrate‐bearing sediments measured under laboratory conditions provide reference values for calibration of logging and seismic exploration results acquired in hydrate‐bearing formations. Instrumented cells were designed for measuring the compressional (P) and shear (S) velocities of sand, silts, and clay with and without hydrate and subject to vertical effective stresses of 0.01 to 2 MPa. Tetrahydrofuran (THF), which is fully miscible in water, was used as the hydrate former to permit close control over the hydrate saturation S hyd and to produce hydrate from dissolved phase, as methane hydrate forms in most natural marine settings. The results demonstrate that laboratory hydrate formation technique controls the pattern of P and S velocity changes with increasing S hyd and that the small‐strain properties of hydrate‐bearing sediments are governed by effective stress, σ ′ v and sediment specific surface. The S velocity increases with hydrate saturation owing to an increase in skeletal shear stiffness, particularly when hydrate saturation exceeds S hyd ≈ 0.4. At very high hydrate saturations, the small strain shear stiffness is determined by the presence of hydrates and becomes insensitive to changes in effective stress. The P velocity increases with hydrate saturation due to the increases in both the shear modulus of the skeleton and the bulk modulus of pore‐filling phases during fluid‐to‐hydrate conversion. Small‐strain Poisson's ratio varies from 0.5 in soft sediments lacking hydrates to 0.25 in stiff sediments (i.e., subject to high vertical effective stress or having high S hyd ). At S hyd ≥ 0.5, hydrate hinders expansion and the loss of sediment stiffness during reduction of vertical effective stress, meaning that hydrate‐rich natural sediments obtained through pressure coring should retain their in situ fabric for some time after core retrieval if the cores are maintained within the hydrate stability field.

Journal of Geophysical Research↗
USGS2010

Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration

Shallow rainfall-induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two-layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure-diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil-water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one-dimensional slope-stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall-induced pore-pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS2011

Spatial connectivity in a highly heterogeneous aquifer: From cores to preferential flow paths

This study investigates connectivity in a small portion of the extremely heterogeneous aquifer at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi. A total of 19 fully penetrating soil cores were collected from a rectangular grid of 4 m by 4 m. Detailed grain size analysis was performed on 5 cm segments of each core, yielding 1740 hydraulic conductivity ( K ) estimates. Three different geostatistical simulation methods were used to generate 3‐D conditional realizations of the K field for the sampled block. Particle tracking calculations showed that the fastest particles, as represented by the first 5% to arrive, converge along preferential flow paths and exit the model domain within preferred areas. These 5% fastest flow paths accounted for about 40% of the flow. The distribution of preferential flow paths and particle exit locations is clearly influenced by the occurrence of clusters formed by interconnected cells with K equal to or greater than the 0.9 decile of the data distribution (10% of the volume). The fraction of particle paths within the high‐ K clusters ranges from 43% to 69%. In variogram‐based K fields, some of the fastest paths are through media with lower K values, suggesting that transport connectivity may not require fully connected zones of relatively homogenous K . The high degree of flow and transport connectivity was confirmed by the values of two groups of connectivity indicators. In particular, the ratio between effective and geometric mean K (on average, about 2) and the ratio between the average arrival time and the arrival time of the fastest particles (on average, about 9) are consistent with flow and advective transport behavior characterized by channeling along preferential flow paths.

Water Resources Research↗