Search USGSSearch

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources

Performance of four turbulence closure models implemented using a generic length scale method

A two-equation turbulence model (one equation for turbulence kinetic energy and a second for a generic turbulence length-scale quantity) proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235] is implemented in a three-dimensional oceanographic model (Regional Oceanographic Modeling System; ROMS v2.0). These two equations, along with several stability functions, can represent many popular turbulence closures, including the k-kl (Mellor-Yamada Level 2.5), k-??, and k-?? schemes. The implementation adds flexibility to the model by providing an unprecedented range of turbulence closure selections in a single 3D oceanographic model and allows comparison and evaluation of turbulence models in an otherwise identical numerical environment. This also allows evaluation of the effect of turbulence models on other processes such as suspended-sediment distribution or ecological processes. Performance of the turbulence models and sediment-transport schemes is investigated with three test cases for (1) steady barotropic flow in a rectangular channel, (2) wind-induced surface mixed-layer deepening in a stratified fluid, and (3) oscillatory stratified pressure-gradient driven flow (estuarine circulation) in a rectangular channel. Results from k-??, k-??, and gen (a new closure proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235]) are very similar for these cases, but the k-kl closure results depend on a wall-proximity function that must be chosen to suit the flow. Greater variations appear in simulations of suspended-sediment concentrations than in salinity simulations because the transport of suspended-sediment amplifies minor variations in the methods. The amplification is caused by the added physics of a vertical settling rate, bottom stress dependent resuspension, and diffusive transport of sediment in regions of well mixed salt and temperature. Despite the amplified sensitivity of sediment to turbulence models in the estuary test case, the four closures investigated here all generated estuarine turbidity maxima that were similar in their shape, location, and concentrations.

Ocean Modelling

Comment on the treatment of residual water content in “A consistent set of parametric models for the two-phase flow of immiscible fluids in the subsurface” by L. Luckner et al.

Luckner et al. [1989] (hereinafter LVN) present a clear summary and generalization of popular formulations used for convenient representation of porous media fluid flow characteristics, including water content ( θ ) related to suction ( h ) and hydraulic conductivity ( K ) related to θ or h. One essential but problematic element in the LVN models is the concept of residual water content ( θ r ; in LVN, θ w , r ). Most studies using θ r determine its value as a fitted parameter and make the assumption that liquid flow processes are negligible at θ values less than θ r . While the LVN paper contributes a valuable discussion of the nature of θ r , it leaves several problems unresolved, including fundamental difficulties in associating a definite physical condition with θ r , practical inadequacies of the models at low θ values, and difficulties in designating a main wetting curve.

Water Resources Research

Pallid sturgeon in the Lower Mississippi Region: Hematology and genome information

This project (Project 1448-43270-2M-002) has been coordinated through the Natchitoches National Fish Hatchery (NNFH) and the U.S. Geological Survey’s National Wetlands Research Center (NWRC). From November 2001 to April 2002, over 280 sturgeon of the genus Scaphirhynchus (including pallid sturgeon, shovelnose, and their hybrids) were sampled from the outflow channel of the Old River Control Structure Complex (ORCC) in Concordia Parish, La. In the overall project, several datasets were collected (see Appendix), including species identification by using microsatellites and morphometric characters, food habits, physical anomalies, information on blood cells, and pathologic evidence of iridovirus – the first indication in the lower Mississippi population of pallid sturgeon. In this study, data on blood cells were obtained from the sturgeon collected monthly from approximately 20 different animals at each sampling time. This report presents preliminary information on differential blood cell identifications in sturgeon, data on comparative genomic DNA content and DNA degradation, and summaries and interpretations of data collected in light of available scientific literature addressing blood parameters of fish and sturgeon, in particular. Results obtained from collection and examination of blood and body fluids are often essential in establishing the health of fish (Blaxhall, 1972; Fange, 1992). Blood cells and sperm cells can be obtained nondestructively from fishes, even from small specimens that weigh less than 100 g (Stoskopf, 1992a). For flow cytometry assays, whereby cells are analyzed individually in a fluid stream, less than 1 :L of blood is needed. Examinations of blood by microscopy and flow cytometry were performed at NWRC in assisting in the efforts directed at recovery of the pallid sturgeon population in the Lower Mississippi River Basin.

Open-File Report

Induced earthquake families reveal distinctive evolutionary patterns near disposal wells

The timing of events in seismic sequences can provide insights into the physical processes controlling fault slip. In southern Kansas, the rate of earthquakes rose rapidly starting in 2013 following expansion of energy production into the area, demanding the disposal of large volumes of wastewater into deep wells. Seismicity catalogs that are complete to low magnitudes can provide insights into the physical processes that induce seismicity near wastewater disposal. We develop a catalog of over 130,000 earthquakes recorded in southern Kansas from mid‐March 2014 through December 2017 by applying a matched filter algorithm to an original catalog of 5,831 template earthquakes. Detections have nearly identical waveforms to their associated template event and represent slip on nearly co‐located sections of a fault. We select template events with at least 100 associated detections and examine the characteristics of these prolific families of earthquakes. We find that families located close (<10 km) to areas with significant volumes of injected fluids have near‐Poissonian interevent times and the families remain active over longer durations. Families farther from high‐volume injection wells show strong clustering of interevent times and shorter sequence durations. We conclude that increasing pore fluid pressures from nearby disposal of large volumes of wastewater is the primary driver of these long duration episodes, with earthquake‐earthquake interactions driving sequences at greater distance from the wells.

Kansas

Lattice-Boltzmann simulation of coalescence-driven island coarsening

The first-order phase separation in a thin fluid film was simulated using a two-dimensional lattice-Boltzman model (LBM) with fluid-fluid interactions. The effects of the domain size on the intermediate asymptotic island size distribution were also discussed. It was observed that the overall process is dominated by coalescence which is independent of island mass. The results show that the combined effects of growth, coalescence, and Ostwald ripening control the phase transition process in the LBM simulations.

Journal of Chemical Physics

An effective medium inversion algorithm for gas hydrate quantification and its application to laboratory and borehole measurements of gas hydrate-bearing sediments

The presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments.

Geophysical Journal International

Debris flow monitoring in the Acquabona watershed on the Dolomites (Italian Alps)

In 1997 a field monitoring system was installed in Acquabona Creek in the Dolomites (Eastern Italian Alps) to observe the hydrologic conditions for debris flow occurrence and some dynamic properties of debris flow. The monitoring system consists of three remote stations: an upper one located at the head of a deeply-incised channel and two others located downstream. The system is equipped with sensors for measuring rainfall, pore pressures in the mobile channel bottom, ground vibrations, debris flow depth, total normal stress and fluid pore-pressure at the base of the flow. Two video cameras record events at the upper channel station and one video is installed at the lowermost station. During summer 1998, three debris flows (volumes from less than 1000 m 3 up to 9000 m 3 ) occurred at Acquabona. The following results were obtained from a preliminary analysis of the data: 1) All of the flows were triggered by rainfalls of less than 1 hour duration, with peak rainfall intensities ranging from 4.8 to 14.7 mm / 10 minute. 2) Debris flows initiated in several reaches of the channel, including the head of the talus slope. 3) The initial surges of the mature flows had a higher solid concentration and a lower velocity (up to 4 m/s) than succeeding, more dilute surges (more than 7 m/s). 4) Total normal stress and pore fluid pressures measured at the base of the flow. (mean depth about 1.1 m) were similar (about 15 kPa), indicating a completely liquefied flow. 5) Peak flows entrained debris at a rate of about 6 m 3/m of channel length and channel bed scouring was proportional to the local slope gradient and was still evident in the lower channel where the slope was 7°.

Acquabona Creek, Alps, Boite River, Dolomites

Seasonal uranium distributions in the coastal waters off the Amazon and Mississippi Rivers

The chemical reactivity of uranium was investigated across estuarine gradients from two of the world’s largest river systems: the Amazon and Mississippi. Concentrations of dissolved (<0.45 μm) uranium (U) were measured in surface waters of the Amazon shelf during rising (March 1990), flood (June 1990) and low (November 1991) discharge regimes. The dissolved U content was also examined in surface waters collected across estuarine gradients of the Mississippi outflow region during April 1992, August 1993, and November (1993). All water samples were analyzed for U by isotope dilution inductively coupled plasma mass spectrometry (ICP-MS). In Amazon shelf surface waters uranium increased nonconservatively from about 0.01 μg I −1 at the river’s mouth to over 3 μg I −1 at the distal site, irrespective of river discharge stage. Observed large-scale U removal at salinities generally less than 15 implies a) that riverine dissolved U was extensively adsorbed by freshly-precipitated hydrous metal oxides (e.g., FeOOH, MnO 2 ) as a result of flocculation and aggregation, and b) that energetic resuspension and reworking of shelf sediments and fluid muds on the Amazon shelf released a chemically reactive particle/colloid to the water column which can further scavenge dissolved U across much of the estuarine gradient. In contrast, the estuarine chemistry of U is inconclusive within surface waters of the Mississippi shelf-break region. U behavior is most likely controlled less by traditional sorption and/or desorption reactions involving metal oxides or colloids than by the river’s variable discharge regime (e.g., water parcel residence time during estuarine mixing, nature of particulates, sediment storage and resuspension in, the confined lower river), and plume dispersal. Mixing of the thin freshwater lens into ambient seawater is largely defined by wind-driven rather than physical processes. As a consequence, in the Mississippi outflow region uranium predominantly displays conservative behavior; removal is evident only during anomalous river discharge regimes. ‘Products-approach’ mixing experiments conducted during the Flood of 1993 suggest the importance of small particles and/or colloids in defining a depleted U versus salinity distribution.

Louisiana

Natural occurrence and significance of fluids indicating high pressure and temperature

Most natural minerals have formed from a fluid phase such as a silicate melt or a saline aqueous solution. Fluid inclusions are tiny volumes of such fluids that were trapped within the growing crystals. These inclusions can provide valuable but sometimes ambiguous data on the temperature, pressure, and composition of these fluids, many of which are not available from any other source. They also provide “visual autoclaves” in which it is possible to watch, through the microscope, the actual phase changes take place as the inclusions are heated. This paper reviews the methods of study and the results obtained, mainly on inclusions formed from highly concentrated solutions, at temperatures ≥500°C. Many such fluids have formed as a result of immiscibility with silicate melt in igneous or high-temperature metamorphic rocks. These include fluids consisting of CO 2 , H 2 O, or hydrosaline melts that were <50% H 2 O. From the fluid inclusion evidence it is clear that a boiling, very hot, very saline fluid was present during the formation of most of the porphyry copper deposits in the world. Similarly, from the inclusion evidence it is clear that early (common) pegmatites formed from essentially silicate melts and that the late, rare-element-bearing and chamber-type pegmatites formed from a hydrosaline melt or a more dilute water solution. The evidence on whether this change in composition from early to late solutions was generally continuous or involved immiscibility is not as clear.

Physics and Chemistry of the Earth

Hematite-bearing materials surrounding Candor Mensa in Candor Chasma, Mars: Implications for hematite origin and post-emplacement modification

The Valles Marineris canyon system on Mars is of enduring scientific interest in part due to the presence of interior mounds that contain extensive layering and water-altered minerals, such as crystalline gray hematite and hydrated sulfates. The presence of hematite and hydrated sulfate minerals is important because their host rock lithologies provide information about past environments that may have supported liquid water and may have been habitable. This work further defines the association and relationship between hematite-bearing materials and low albedo (presumably aeolian) deposits and layered materials, identifies physical characteristics that are strongly correlated with the presence of hematite, and refines hypotheses for the origin and post-emplacement modification (including transport) of these hematite-bearing and associated materials. There are only three regions surrounding Candor Mensa where hematite has been identified, even though morphologic properties are similar throughout the entire mensa. Three possible explanations for why hematite is only exposed in these regions include: (1) the topographic structure of the mensa walls concentrates hematite at the base of the layered deposits, influencing the ability to detect hematite from orbit; (2) the presence of differing amounts of &ldquo;dark mantling material&rdquo; and hematite-free erosional sediment; (3) the potential fracturing of the mensa and the influence of these structures on fluid flow and subsequent digenesis. The observations of hematite-bearing materials in this work support the hypothesis that hematite is eroding from a unit in the Candor Mensa interior layered deposits (ILD) and is being concentrated as a lag deposit adjacent to the lower layers of Candor Mensa and at the base in the form of dark aeolian material. Due to the similar geologic context associated with hematite-bearing and ILD materials throughout the Valles Marineris canyon system, the insight gained from studying these materials surrounding Candor Mensa can likely be applicable to similar layered deposits throughout Valles Marineris.

Icarus

Intensity of quartz cathodoluminescence and trace-element content in quartz from the porphyry copper deposit at Butte, Montana

Textures of hydrothermal quartz revealed by cathodoluminescence using a scanning electron microscope (SEM-CL) reflect the physical and chemical environment of quartz formation. Variations in intensity of SEM-CL can be used to distinguish among quartz from superimposed mineralization events in a single vein. In this study, we present a technique to quantify the cathodoluminescent intensity of quartz within individual and among multiple samples to relate luminescence intensity to specific mineralizing events. This technique has been applied to plutonic quartz and three generations of hydrothermal veins at the porphyry copper deposit in Butte, Montana. Analyzed veins include early quartz-molybdenite veins with potassic alteration, pyrite-quartz veins with sericitic alteration, and Main Stage veins with intense sericitic alteration. CL intensity of quartz is diagnostic of each mineralizing event and can be used to fingerprint quartz and its fluid inclusions, isotopes, trace elements, etc., from specific mineralizing episodes. Furthermore, CL intensity increases proportional to temperature of quartz formation, such that plutonic quartz from the Butte quartz monzonite (BQM) that crystallized at temperatures near 750 ??C luminesces with the highest intensity, whereas quartz that precipitated at ???250 ??C in Main Stage veins luminesces with the least intensity. Trace-element analyses via electron microprobe and laser ablation-ICP-MS indicate that plutonic quartz and each generation of hydrothermal quartz from Butte is dominated by characteristic trace amounts of Al, P, Ti, and Fe. Thus, in addition to CL intensity, each generation of quartz can be distinguished based on its unique trace-element content. Aluminum is generally the most abundant element in all generations of quartz, typically between 50 and 200 ppm, but low-temperature, Main Stage quartz containing 400 to 3600 ppm Al is enriched by an order of magnitude relative to all other quartz generations. Phosphorous is present in abundances between 25 and 75 ppm, and P concentrations in quartz show little variation among quartz generations. Iron is the least abundant of these elements in most quartz types and is slightly enriched in CL-dark quartz in pyrite-quartz veins with sericitic alteration. Titanium is directly correlated with both temperature of quartz precipitation, and intensity of quartz luminescence, such that BQM quartz contains hundreds of ppm Ti, whereas Main Stage quartz contains less than 10 ppm Ti. Our results suggest that Ti concentration in quartz is controlled by temperature of quartz precipitation and that increased Ti concentrations in quartz may be responsible for increased CL intensities.

American Mineralogist

Fracture-mesh faulting in the swarm-like 2020 Maacama sequence revealed by high-precision earthquake detection, location, and focal mechanisms

In August of 2020, an earthquake sequence initiated within the Maacama fault zone in northern California, raising questions about its relationship with the larger-scale fault. To investigate the faulting geometry and its implications for physical processes driving seismicity, we applied an integrated, multi-faceted seismic analysis including waveform-correlation-based event detection, relative relocation, and high-precision focal mechanisms. To determine mechanisms for a large population of very small earthquakes (predominantly M < 1), we combined correlation-derived polarity analysis with a recently developed technique to determine S/P amplitude ratios from single seismic components. Finally, we applied an iterative stress inversion to distinguish between the likely fault and auxiliary planes. This analysis reveals that although the sequence initiated on a right-lateral fault, plausibly a strand of the Maacama Fault, it also activated numerous en echelon left-lateral conjugate faults. Together, these interlocking faults form a fracture mesh, consistent with structures associated elsewhere with fluid-induced faulting and earthquake swarms.

California

Relationships between sliding behavior and internal geometry of laboratory fault zones and some creeping and locked strike-slip faults of California

In order to relate fault geometries to sliding behavior, maps of recently active breaks within the Hayward fault of central California, which is characterized by fault creep, have been examined and compared to maps of the San Andreas fault. The patterns of recent breaks of the Hayward fault are consistent with those found within the creeping section of the San Andreas, and they appear to have plausible physical explanations in the findings of laboratory experiments. The distinguishing geometric features of the examined locked and creeping faults are: (1) P-type second-order traces predominate over R(Riedel)-type traces in creeping sections; and (2) R-type second-order traces make smaller angles to the local fault strike in creeping sections than they do in locked sections. Two different maps of the Hayward fault gave similar results, supporting the inference that the patterns identified are basic characteristics of the fault rather than artifacts of a particular mapping procedure. P shears predominate over R shears under laboratory conditions that allow dilation within the fault zone. In our own experiments, P-shear development was favored by the generation of excess pore-fluid pressures. We propose that creep in California faults also is the result of fluid overpressures that are maintained in a low-permeability gouge zone and that significantly lower effective stresses, thus helping to stabilize slip and producing high values of the ratio P/R. Small R-trace angles may also be an indicator of low effective stresses, but the evidence for this is not conclusive because other factors can also affect the size of the angles.

California

Geology of the Marble exploration hole 4, Nevada Test Site, Nye County, Nevada

The Marble exploration hole 4 was drilled to determine the chemical and physical properties of marble from the surface to depth of about 1,200 feet. The drill hole is in the north-central part of the Tippipah Spring quadrangle, approximately 40 miles north of Mercury, Nev. The collar of the hole is at an elevation of 5,490.28 feet and has the Nevada State grid coordinates of N. 904,480.13 and E. 637,758.38. The hole is accessible by a graded dirt road easily traversed by heavy construction equipment, that connects to the highway between Mercury and Groom Lake. The hole was core drilled by Minerals Engineering Company, Grand Junction, Colo., with a 2,000-foot capacity Portadrill drilling rig and conventional NX coring equipment. Drilling was started July 13, and completed August 30, 1959. The drilling fluid was water for most of the hole. At some places in the hole, where loss of drilling fluid was high, the water was thickened with bentonitic mud. At 1,103 feet, presumed caving of the hole caused the drill stem to stick so that it could not be withdrawn; subsequently about 920 feet of the drill stem was recovered. The hole was then reamed to 800 feet and cemented from 785 to 855 feet. A whipstock was emplaced at depth of 805 feet so that the hole, which was vertical to this depth, would be deflected 4&deg; from the vertical. Core was recovered for this part of the hole from 1,088 to 1,187 feet, the total depth. Seven feet of 8-inch outside diameter casing was set from the collar. This report summarizes the information obtained during preparation of the lithologic log of the core and presents results of chemical analyses of marble samples collected from surface near the drill hole. The report was prepared by the U.S. Geological Survey on behalf of the Albuquerque Operations Office, U.S. Atomic Energy Commission. The writers acknowledge the assistance of Mr. John Foster, drilling foreman for Minerals Engineering Company and Mr. Walter A. Johnson, field engineer for Holmes and Narver, Inc., the engineering-contracting firm.

Nevada

Downhole physical property-based description of a gas hydrate petroleum system in NGHP-02 Area C: A channel, levee, fan complex in the Krishna-Godavari Basin offshore eastern India

India’s second National Gas Hydrate Program expedition, NGHP-02, collected logging while drilling and sediment core data in Area C offshore eastern India, to investigate controls on the distribution and peak saturations of methane gas hydrate occurrences in buried channel, levee and fan deposits. Physical property results are presented here for the four Area C coring sites: NGHP-02-07, targeting an upper continental-slope channel deposit; NGHP-02-08 and -09, targeting levee deposits on either side of a channel further downslope, and NGHP-02-05, targeting a sequence of fan deposits extending out from the slope base. Coarse-grained sediment exists at each site, but site-specific differences in clay distribution provide significant controls on the gas hydrate distribution and saturation. At NGHP-02-07, only the upper ~4 m of a ~42 m-thick, relatively low clay-content, coarse-grained interval is inferred to be gas hydrate-bearing. NGHP-02-07 has a relatively thin, high-permeability overburden seal, and methane-rich fluid likely leaks from the primary reservoir. NGHP-02-08’s levee deposit seal is similarly permeable near the reservoir, but becomes less permeable toward the seafloor. Relative to NGHP-02-07, methane is retained more effectively in the NGHP-02-08 reservoir, but that reservoir is interbedded with layers of high clay-content, low gas hydrate saturation sediment, limiting the maximum gas hydrate content for NGHP-02-08. NGHP-02-09, the second levee deposit site, has a thicker, less permeable overburden than NGHP-02-08, combined with >50 m-thick, low clay-content reservoir sediments. Correspondingly, NGHP-02-09 has a thicker gas hydrate-bearing reservoir with consistently higher gas hydrate saturations than NGHP-02-08. NGHP-02-05 has abundant coarse-grained material spread over nearly the entire drilling interval, but the sediment is poorly sorted. Gas hydrate is distributed among several primarily coarse-grained layers, but gas hydrate saturations are limited by relatively high clay contents, and an overlying seal that is too thin and permeable to effectively retain methane in the reservoir.

Marine and Petroleum Geology

Detection of crystalline hematite mineralization on Mars by the Thermal Emission Spectrometer: evidence for near-surface water

The Thermal Emission Spectrometer (TES) instrument on the Mars Global Surveyor (MGS) mission has discovered a remarkable accumulation of crystalline hematite (α-Fe 2 O 3 ) that covers an area with very sharp boundaries approximately 350 by 350–750 km in size centered near 2°S latitude between 0° and 5°W longitude (Sinus Meridiani). Crystalline hematite is uniquely identified by the presence of fundamental vibrational absorption features centered near 300, 450, and >525 cm −1 and by the absence of silicate fundamentals in the 1000 cm −1 region. Spectral features resulting from atmospheric CO 2 , dust, and water ice were removed using a radiative transfer model. The spectral properties unique to Sinus Meridiani were emphasized by removing the average spectrum of the surrounding region. The depth and shape of the hematite fundamental bands show that the hematite is crystalline and relatively coarse grained (>5–10 μm). Diameters up to and greater than hundreds of micrometers are permitted within the instrumental noise and natural variability of hematite spectra. Hematite particles <5–10 μm in diameter (as either unpacked or hard-packed powders) fail to match the TES spectra. The spectrally derived areal abundance of hematite varies with particle size from ∼10% (>30 μm diameter) to 40–60% (10 μm diameter). The hematite in Sinus Meridiani is thus distinct from the fine-grained (diameter <5–10 μm), red, crystalline hematite considered, on the basis of visible, near-IR data, to be a minor spectral component in Martian bright regions like Olympus-Amazonis. Sinus Meridiani hematite is closely associated with a smooth, layered, friable surface that is interpreted to be sedimentary in origin. This material may be the uppermost surface in the region, indicating that it might be a late stage sedimentary unit or a layered portion of the heavily cratered plains units. We consider five possible mechanisms for the formation of coarse-grained, crystalline hematite. These processes fall into two classes depending on whether they require a significant amount of near-surface water: the first is chemical precipitation that includes origin by (1) precipitation from standing, oxygenated, Fe-rich water (oxide iron formations), (2) precipitation from Fe-rich hydrothermal fluids, (3) low-temperature dissolution and precipitation through mobile ground water leaching, and (4) formation of surface coatings, and the second is thermal oxidation of magnetite-rich lavas. Weathering and alteration processes, which produce nanophase and red hematite, are not consistent with the coarse, crystalline hematite observed in Sinus Meridiani. We prefer chemical precipitation models and favor precipitation from Fe-rich water on the basis of the probable association with sedimentary materials, large geographic size, distance from a regional heat source, and lack of evidence for extensive groundwater processes elsewhere on Mars. The TES results thus provide mineralogic evidence for probable large-scale water interactions. The Sinus Meridiani region may be an ideal candidate for future landed missions searching for biotic and prebiotic environments, and the physical characteristics of this site satisfy all of the engineering requirements for the missions currently planned.

Journal of Geophysical Research E: Planets

Elements of an improved model of debris-flow motion

A new depth-averaged model of debris-flow motion describes simultaneous evolution of flow velocity and depth, solid and fluid volume fractions, and pore-fluid pressure. Non-hydrostatic pore-fluid pressure is produced by dilatancy, a state-dependent property that links the depth-averaged shear rate and volumetric strain rate of the granular phase. Pore-pressure changes caused by shearing allow the model to exhibit rate-dependent flow resistance, despite the fact that the basal shear traction involves only rate-independent Coulomb friction. An analytical solution of simplified model equations shows that the onset of downslope motion can be accelerated or retarded by pore-pressure change, contingent on whether dilatancy is positive or negative. A different analytical solution shows that such effects will likely be muted if downslope motion continues long enough, because dilatancy then evolves toward zero, and volume fractions and pore pressure concurrently evolve toward steady states. ?? 2009 American Institute of Physics.

Conference Paper

The stable isotope geochemistry of acid sulfate alteration

Acid sulfate wall-rock alteration, characterized by the assemblage alunite + kaolinite + quartz + or - pyrite, results from base leaching by fluids concentrated in H 2 SO 4 . Requisite amounts of H 2 SO 4 can be generated by different mechanisms in three principal geologic environments: (1) by atmospheric oxidation of sulfides in the supergene environment, (2) by atmospheric oxidation at the water table in the steam-heated environment of H 2 S released by deeper, boiling fluids, and (3) by the disproportionation of magmatic SO 2 to H 2 S and H 2 SO 4 during condensation of a magmatic vapor plume at intermediate depths in magmatic hydrothermal environments in silicic and andesitic volcanic terranes. In addition, coarse vein alunite may form in a magmatic steam environment from rapid release of an SO 2 -rich magmatic vapor phase at high temperature and low pressure or from the oxidation of a more reduced magmatic vapor by entrained atmospheric oxygen in the carapace of a volcanic edifice.Alunite [KAl 3 (SO 4 ) 2 (OH) 6 ] contains four stable isotope sites and complete analyses (delta D, delta 18 O (sub SO 4 ) , delta 18 O OH , and delta 34 S) are now possible. Except for delta 18 O OH in magmatic hydrothermal alunites, primary values are usually retained. In cooperation with many colleagues, over 500 measurements have been made on nearly 200 samples of alunite and associated minerals from 23 localities, and 55 additional analyses have been taken from the literature. This survey confirms that kinetic factors play an important role in the stable isotope systematics of alunite and acid sulfate alteration. To a very large extent they form the isotopic basis for distinguishing between environments of acid sulfate alteration, and they provide important insights into attendant processes. Stable isotope analyses of alunite, often in combination with those on associated sulfides and kaolinite, permit recognition of environments of formation and provide information on origins of components, processes (including rates), physical-chemical environments, and temperatures of formation.Supergene acid sulfate alteration may form over any sulfide zone when it is raised above the water table by tectonics or exposed by erosion. It may overprint earlier acid sulfate assemblages, particularly the magmatic hydrothermal assemblages which are pyrite rich such as at El Salvador, Chile; Rodalquilar, Spain; and Goldfield, Nevada. Supergene alunite normally has delta 34 S values virtually identical to precursor sulfides unless bacteriogenic reduction of aqueous sulfate takes place in standing pools of water. delta D values are close to that of local meteoric water unless extensive evaporation occurs. delta D and delta 18 O OH values of supergene alunites from a range of latitudes fall in a broad zone parallel to the meteoric water line much the way delta D and delta 18 O values of associated halloysite-kaolinite fall near the kaolinite line of Savin and Epstein (1970). delta 18 O (sub SO 4 ) values are kinetically controlled and will reflect the hydro-geochemistry of the environment. delta 18 O (sub SO 4 -OH) alues are grossly out of equilibrium and large negative values are definitive of a supergene origin.In steam-heated environments, such as those at the Tolfa district, Italy, and Marysvale, Utah, and numerous modern geothermal systems, acid sulfate alteration zones are characterized by pronounced vertical zoning. Such acid sulfate alteration may occur over adularia-sericite-type base and precious metal ore deposits such as the one at Buckskin, Nevada. Initial delta 18 O (sub SO 4 ) and delta 34 S values are kinetically controlled, but delta 18 O (sub SO 4 ) values usually reach equilibrium with fluids, and even delta 34 S values may reflect partial exchange with H 2 S where the residence time of aqueous sulfate is sufficient. Most alunites of steam-heated origin have delta 34 S values the same as those of precursor H 2 S (and as related sulfides, if present) and delta D values similar to that of local meteoric water. In the samples analyzed, most delta 18 O (sub SO 4 -OH) values give reasonable temperatures of 90 degrees to 160 degrees C, indicating that delta 18 O (sub SO 4 ) and delta 18 O OH values reflect a close approach to equilibrium with the fluid. The delta 18 O (sub SO 4 ) and delta 18 O OH values also reflect the degree of exchange of the meteoric fluids with wall rock. Coeval kaolinites typically have delta 18 O and delta D values to the left of the kaolinitc line.Magmatic hydrothermal, acid sulfate alteration zones in near-surface epithermal deposits such as Summitville, Colorado. Julcani, Peru, and Red Mountain and Lake City, Colorado, are characterized by vertical orientation and horizontal zoning, the presence of coeval pyrite, PO 4 analogues of alunite, zunyite, and later gold, pyrite and enargite, and often other Cu-As-Sb-S minerals. Acid sulfate alteration assemblages also occur as late stages in the porphyry-copper deposit at E1 Salvador, Chile. In the examples studied, magmatic hydrothermal alunites have delta D values close to those for magmatic water. delta 34 S values are 16 to 28 per mil larger than those for associated pyrite, reflecting equilibrium between aqueous H 2 S and SO 4 formed by the disproportionation of magmatically derived SO 2 . delta 18 O (sub SO 4 ) values are usually 8 to 18 per mil and vary systematically with delta 34 S values, reflecting variations in temperature and/or H 2 S/SO 4 fluid ratios. Further variation in delta 18 O (sub SO 4 ) values may result if SO 2 condenses in mixed magmatic meteoric water fluids. delta 18 O (sub SO (sub 4-) OH) values of magmatic hydrothermal alunites are generally unsuitable for temperature determinations because of retrograde exchange in the OH site, but delta 34 S (sub alunite-pyrite) values provide reliable temperature estimates.Magmatic steam environments appear to occur over a range of depths and are characterized by monomineralic veins of coarse alunite in variably alunitized and kaolinized wall rocks containing minor pyrite. Alunite formed in the magmatic steam environment can usually be recognized by delta 34 S near delta 34 S (sub Sigma S) values and delta D and delta 18 O (sub SO 4 ) values near magmatic values. Magmatic steam alunite differs from magmatic hydrothermal alunite by having delta 34 S close to delta 34 S (sub Sigma S) values of the system. delta 18 O (sub SO (sub 4-) OH) values of most magmatic steam alunite give temperatures ranging from 90 degrees to 210 degrees C but, for reasons which are not understood, some temperatures as well as calculated delta 18 O (sub H 2 O) values are too low for presumed precipitation from a magmatic vapor phase. Magmatic steam environments may occur over porphyry-type mineralization as at Red Mountain, Colorado, and Alunite Ridge, Utah, and over or adjacent to adularia-sericite-type deposits in volcanic domes as at Cactus, California.

Economic Geology