Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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 433 records · Page 24Linked to original sources
USGS1996

Transverse dispersion of contaminants in fractured permeable formations

Our interest in understanding some of the mechanisms involved in the mineralization of the Great Bend Prairie aquifer of Kansas by salt water originating from Permian bedrock formations, which are fractured sandstones, has prompted this basic study. The fractured permeable formatiodis represented by a simplified conceptual model incorporating two sets of oblique and parallel fractures embedded in permeable blocks. The domain is initially divided into a completely freshwater zone overlying a completely saline water zone. However, the sharp interface originally existing between the saline and fresh water is subject to dispersion because of mixing in fracture intersections and between the fracture flow and the permeable block flow. Simulations based on the use of an appropriate numerical model developed in the present study have helped us to characterize dispersion of the sharp interface and creation of the transition zone. Relationships of transverse and longitudinal dispersion in the domain are also determined.

Journal of Hydrology↗
USGS1995

Seasonal cycles of dissolved constituents in streamwater in two forested catchments in the mid-Atlantic region of the eastern U.S.A.

Streamwater discharge and chemistry of two small catchments on Catoctin Mountain in north-central Maryland have been monitored since 1982. Repetitive seasonal cycles in stream-water chemistry have been observed each year, along with seasonal cycles in the volume of stream discharge and in groundwater levels. The hypothesis that the observed streamwater chemical cycles are related to seasonal changes in the hydrological flow paths that contribute to streamflow is examined using a combination of data on groundwater levels, shallow and deep groundwater chemistry, streamwater discharge, streamwater chemistry, soil-water chemistry, and estimates of water residence times. The concentrations of constituents derived from rock weathering, particularly bicarbonate and silica, increase in streamwater during the summer when the water table is below the regolith-bedrock interface and stream discharge consists primarily of deep groundwater from the fractured-bedrock aquifer. Conversely, the concentrations in streamwater of atmospherically derived components, particularly sulfate, increase in winter when the water table is above the regolith-bedrock interface and stream discharge consists primarily of shallow groundwater from the regolith. Tritium and chlorofluorocarbon (CFC) measurements suggest that the groundwater in these systems is young, with a residence time of less than several years. The results of this study have implications for the design of large-scale water-quality monitoring programs.

Maryland↗
USGS1991

Ore microscopy of the Paoli silver-copper deposit, Oklahoma

The Paoli silver-copper deposit is located in south-central Oklahoma, 56 km south-southeast from Norman, Oklahoma. It was mined for high-grade silver-copper near the beginning of this century, and intensive exploratory drilling during the early 1970's delineated unmined portions of the deposit. A collaborative study between the U.S.G.S., the Kansas Geological Survey, and the University of Missouri-Rolla was undertaken to provide new information on the character of red bed copper deposits of the Midcontinent region. The Paoli deposit has been interpreted to occur as a roll-front type of deposit. The silver and copper mineralization occurs within paleochannels in the Permian Wellington Formation. The silver-copper interfaces appear to be controlled by oxidation-reduction interfaces that are marked by grey to red color changes in the host sandstone. Ore microscopic examinations of polished thin sections show that unoxidized ore consists of chalcocite, digenite, chalcopyrite, covellite and pyrite; and oxidized ores are characterized by covellite, bornite, hematite and goethite. In sandstone-hosted ores, chalcocite and digenite replace dolomite and border clastic quartz grains. In siltstone-hosted ores, the copper sulfide grains have varied shapes; most are irregular in shape and 5–25 μm across, others have euhedral shapes suggestive of pyrite crystal replacements, and some are crudely spherical and are 120–200 μm across. Chalcopyrite is the predominant copper sulfide at depth. Covellite and malachite replace chalcocite and digenite near the surface. Silver only occurs as native silver; most as irregularly shaped grains 40–80 μm across, but some as cruciform crystals that are up to 3.5 mm across. The native silver has been deposited after copper sulfides, and locally replaces chalcocite. Surficial nodules of pyrite, malachite and hematite locally are present in outcrops at the oxidation-reduction fronts. Polished sections of the nodules show that malachite forms a cement around quartz sand grains, and brecciated pyrite grains are surrounded by rims of hematite and goethite. Dolomite is the principal sandstone cement. Cathodoluminescence microscopic study of the mineral has shown that it was deposited during seven periods before the copper sulfide mineralization.

Oklahoma↗
USGS1992

Distribution and effects of shallow gas on bulk estuarine sediment properties

Gas bubble are present in sediments covering approximately 30% of the main stem of Chesapeake Bay, with bubbles occurring at the sediment-water interface in 18% of the main stem sediments. This biogenic gas is found either in the sediments in the lower salinity reaches of the Bay, or confined to sediments which overline infilled palaeodrainage channels formed during the Wisconsinan low sea level stand (approximately 18 ka). Gas associated with the old drainage network does not correlate with present bathymetry or sedimentological patterns. Some differences between the gas-charged and gas-free sediments are: (1) gas-charged sediments have water contents 10-20% higher than comparable gas-free cores; (2) organic matter is better presented with depth in the gas-charged sediments (upwards of 60% more at one depth); (3 monosulphides are dominant sulphide mineral phase within the gas-charged sediments, comprising over 40% of the total sulphur. Within the gas-free sediments monosulphides are significant only near the sediment-water interface and rapidly become negligible with depth, and; (4) cores of gas-charged sediments are highly colour-banded due to preservation of sulphide mineral variations, while gas-free cores are diagenetically altered to pyrite.

Maryland↗
USGS1994

Use of geoelectrical methods in groundwater pollution surveys in a coastal environment

The pollution of coastal aquifers by old landfills can contaminate valuable and scarce water resources in the freshwater lens utilized seasonably by overcrowded communities. The pollutants will ultimately flow into the sea where they may also cause a coastal water pollution problem. We have detected pollution in the freshwater lens from a sanitary landfill near Provincetown, Cape Cod, using the geoelectrical resistivity method. This survey included Schlumberger geoelectrical depth soundings and a horizontal geoelectrical profile using the Wenner configuration. The geoelectrical survey was conducted at a site along Highway 6 where it passes the coastal town of Provincetown and a sanitary landfill that has been in operation since 1954. The depth soundings suggest the characteristic decrease in resistivity vs. depth from the high resistivity of the unsaturated zone to the low resistivity of the saltwater saturated zone. The freshwater lens is clearly identified by the change in slope of the steeply dipping curve of resistivity versus electrode spacing. Interpretations made using a multilayer program, Geomate, resulted in layer resistivities between 460 and 95 ohm?m for the freshwater lens. A comparison with well water resistivities suggests that a layer resistivity of 230 ohm?m or lower is indicative of pollution in the freshwater lens. The results of the geoelectrical depth soundings were confirmed in the Wenner horizontal profile. Both measurements suggest that the pollutants do not spread evenly as one would expect for a homogeneous and isotropic medium. Instead, a preferred channel for the flow of the pollutants is observed along a path from the landfill toward the shoreline. The depth to the saltwater/freshwater interface or, more specifically, to the low resistivity-high resistivity interface appears to be shallow where the freshwater lens is polluted. This was confirmed by pore water well samples that were highly mineralized. The equilibrium postulated by the Ghyben-Herzberg relation appears to be disturbed in the area of aquifer pollution. This rise in the conductivity boundary is caused by the highly mineralized bottom of the contaminant plume that submerges into the saltwater saturated zone. In the area of high freshwater pollution the groundwater can be subdivided into three layers that show a decrease in resistivity with depth. The formation factor, F, defined as the ratio of bulk aquifer resistivity to pore water resistivity, shows unusually high values between 10 and 12. These high values are unexpected for an unconsolidated sand. Pollution residues are suspected to clog the pores and thus to increase the resistivity. It is possible that iron-oxidizing bacteria and the precipitation of dissolved iron or organic pollutants are the cause of the high values of F. If proven correct, these interesting possibilities could lead to future new applications of the geoelectrical resistivity method in contaminant hydroloy.

Journal of Applied Geophysics↗
USGS1992

Holocene coseismic and aseismic uplift of Isla Mocha, south-central Chile

During the past 6000 years Isla Mocha, a 12 km-long island 30 km off the coast of south-central Chile, experienced a 38 m fall of relative sea level caused primarily by rapid tectonic uplift of the island. As many as 18 raised shorelines (strandlines) record this uplift. Historic accounts of uplift during the great earthquakes (M > 8) of 1835 and 1960 suggest some of the more prominent prehistoric strandlines also emerged during great earthquakes on the interface between the Nazca and South America plates. But the close elevational spacing of strandlines, subdued morphology of strandline beaches, scarcity of exposed bedrock wave-cut platforms, and the extremely high rates of aseismic uplift (ca. 70 mm/yr) of the island since the last great earthquake suggest that many strandlines were raised by aseismic rather than coseismic uplift. Strandline heights and 14 new radiocarbon ages on marine shells show that the present-day uplift rate is more than three times the net rate (ca. 20 mm/yr) of the past 1000 years. The recent high rate probably reflects increased aseismic slip on an inferred thrust fault in the overriding South America plate. Isla Mocha overlies an area of high stress concentration between two major segments of the Chilean subduction zone. The inferred high rate of slip on the thrust fault may be a response to stress changes on the plate interface near the boundary between the segments. ?? 1992.

Quaternary International↗
USGS2022

Wetlands under global change

Wetlands are among the ecosystem types most threatened by global change, including both climate change and other anthropogenic factors such as sea level rise, urban development, deforestation, agricultural land use, drainage, levees, tidal flow restrictions, pollution, eutrophication, and fires. Wetlands not only store disproportionate amounts of carbon compared to other terrestrial ecosystems, but they lie at the terrestrial-aquatic interface crucial to understanding landscape and global scale biogeochemical cycles. In this chapter, we focus on the major global change factors affecting wetlands and the responses of different wetland types to those global change factors. Special attention is given to direct responses to increasing atmospheric carbon dioxide levels. Because of their hydrological connections and placement at the terrestrial-aquatic interface, the conservation of wetlands involves accounting for uncertainties related to interacting stressors. While the past decades have seen many important experimental and observational studies of wetland responses to global change factors, large uncertainties remain, especially within tropical regions where even the basic extent of wetland ecosystems is not well documented.

Book chapter↗
USGS2015

PhreeqcRM: A reaction module for transport simulators based on the geochemical model PHREEQC

PhreeqcRM is a geochemical reaction module designed specifically to perform equilibrium and kinetic reaction calculations for reactive transport simulators that use an operator-splitting approach. The basic function of the reaction module is to take component concentrations from the model cells of the transport simulator, run geochemical reactions, and return updated component concentrations to the transport simulator. If multicomponent diffusion is modeled (e.g., Nernst–Planck equation), then aqueous species concentrations can be used instead of component concentrations. The reaction capabilities are a complete implementation of the reaction capabilities of PHREEQC. In each cell, the reaction module maintains the composition of all of the reactants, which may include minerals, exchangers, surface complexers, gas phases, solid solutions, and user-defined kinetic reactants. PhreeqcRM assigns initial and boundary conditions for model cells based on standard PHREEQC input definitions (files or strings) of chemical compositions of solutions and reactants. Additional PhreeqcRM capabilities include methods to eliminate reaction calculations for inactive parts of a model domain, transfer concentrations and other model properties, and retrieve selected results. The module demonstrates good scalability for parallel processing by using multiprocessing with MPI (message passing interface) on distributed memory systems, and limited scalability using multithreading with OpenMP on shared memory systems. PhreeqcRM is written in C++, but interfaces allow methods to be called from C or Fortran. By using the PhreeqcRM reaction module, an existing multicomponent transport simulator can be extended to simulate a wide range of geochemical reactions. Results of the implementation of PhreeqcRM as the reaction engine for transport simulators PHAST and FEFLOW are shown by using an analytical solution and the reactive transport benchmark of MoMaS.

Advances in Water Resources↗
USGS2009

Modeling Carbon Dioxide, pH and Un-Ionized Ammonia Relationships in Serial Reuse Systems

In serial reuse systems, excretion of metabolic carbon dioxide has a significant impact on ambient pH, carbon dioxide, and un-ionized ammonia concentrations. This impact depends strongly on alkalinity, water flow rate, feeding rate, and loss of carbon dioxide to the atmosphere. A reduction in pH from metabolic carbon dioxide can significantly reduce the un-ionized ammonia concentration and increase the carbon dioxide concentrations compared to those parameters computed from influent pH. The ability to accurately predict pH in serial reuse systems is critical to their design and effective operation. A trial and error solution to the alkalinity–pH system was used to estimate important water quality parameters in serial reuse systems. Transfer of oxygen and carbon dioxide across the air–water interface, at overflow weirs, and impacts of substrate-attached algae and suspended bacteria were modeled. Gas transfer at the weirs was much greater than transfer across the air–water boundary. This simulation model can rapidly estimate influent and effluent concentrations of dissolved oxygen, carbon dioxide, and un-ionized ammonia as a function of water temperature, elevation, water flow, and weir type. The accuracy of the estimates strongly depends on assumed pollutional loading rates and gas transfer at the weirs. The current simulation model is based on mean daily loading rates; the impacts of daily variation loading rates are discussed. Copies of the source code and executable program are available free of charge.

Aquacultural Engineering↗
USGS2015

The data quality analyzer: a quality control program for seismic data

The U.S. Geological Survey's Albuquerque Seismological Laboratory (ASL) has several initiatives underway to enhance and track the quality of data produced from ASL seismic stations and to improve communication about data problems to the user community. The Data Quality Analyzer (DQA) is one such development and is designed to characterize seismic station data quality in a quantitative and automated manner. The DQA consists of a metric calculator, a PostgreSQL database, and a Web interface: The metric calculator, SEEDscan, is a Java application that reads and processes miniSEED data and generates metrics based on a configuration file. SEEDscan compares hashes of metadata and data to detect changes in either and performs subsequent recalculations as needed. This ensures that the metric values are up to date and accurate. SEEDscan can be run as a scheduled task or on demand. The PostgreSQL database acts as a central hub where metric values and limited station descriptions are stored at the channel level with one-day granularity. The Web interface dynamically loads station data from the database and allows the user to make requests for time periods of interest, review specific networks and stations, plot metrics as a function of time, and adjust the contribution of various metrics to the overall quality grade of the station. The quantification of data quality is based on the evaluation of various metrics (e.g., timing quality, daily noise levels relative to long-term noise models, and comparisons between broadband data and event synthetics). Users may select which metrics contribute to the assessment and those metrics are aggregated into a “grade” for each station. The DQA is being actively used for station diagnostics and evaluation based on the completed metrics (availability, gap count, timing quality, deviation from a global noise model, deviation from a station noise model, coherence between co-located sensors, and comparison between broadband data and synthetics for earthquakes) on stations in the Global Seismographic Network and Advanced National Seismic System.

Computers & Geosciences↗
USGS2004

Dynamic modelling of an adsorption storage tank using a hybrid approach combining computational fluid dynamics and process simulation

A computational fluid dynamics (CFD) software package has been coupled with the dynamic process simulator of an adsorption storage tank for methane fuelled vehicles. The two solvers run as independent processes and handle non-overlapping portions of the computational domain. The codes exchange data on the boundary interface of the two domains to ensure continuity of the solution and of its gradient. A software interface was developed to dynamically suspend and activate each process as necessary, and be responsible for data exchange and process synchronization. This hybrid computational tool has been successfully employed to accurately simulate the discharge of a new tank design and evaluate its performance. The case study presented here shows that CFD and process simulation are highly complementary computational tools, and that there are clear benefits to be gained from a close integration of the two. ?? 2004 Elsevier Ltd. All rights reserved.

Computers and Chemical Engineering↗
USGS2007

Barite-forming environments along a rifted continental margin, Southern California Borderland

The Southern California Continental Borderland (SCCB) is part of the broad San Andreas transform-fault plate boundary that consists of a series of fault-bounded, petroleum-generating basins. The SCCB has high heat flow and geothermal gradients produced by thinned continental crust and Neogene volcanism. Barite deposits in the SCCB occur along faults. Barite samples from two sea-cliff sites and four offshore sites in the SCCB were analyzed for mineralogy, chemical (54 elements) and isotopic (S, Sr) compositions, and petrography. Barite from Palos Verdes (PV) Peninsula sea-cliff outcrops is hosted by the Miocene Monterey Formation and underlying basalt; carbonate rocks from those outcrops were analyzed for C, O, and Sr isotopes and the basalt for S isotopes. Cold-seep barite from Monterey Bay, California was analyzed for comparison. SCCB offshore samples occur at water depths from about 500 to 1800 m. Those barites vary significantly in texture and occurrence, from friable, highly porous actively growing seafloor mounds to dense, brecciated, vein barite. This latter type of barite contrasts with cold-seep barite in being much more coarse grained, forms thick veins in places, and completely replaced rock clasts in breccia. The barite samples range from 94 to 99 wt% BaSO4, with low trace-element contents, except for high Sr, Zr, Br, U, and Hg concentrations compared to their crustal abundances. δ34S for SCCB offshore barites range from 21.6‰ to 67.4‰, and for PV barite from 62‰ to 70‰. Pyrite from PV sea-cliff basalt and sedimentary rocks that host the barites averages 7.8‰ and 2.2‰, respectively. Two offshore barite samples have δ34S values (21.6‰, 22.1‰) close to that of modern seawater sulfate, whereas all other samples are enriched to strongly enriched in 34S. 87Sr/86Sr ratios for the barites vary over a narrow range of 0.70830–0.70856 and are much lower than that of modern seawater and also lower than the middle Miocene seawater ratio, the time of deposition of the host rocks for the PV samples. This indicates that the mineralizing fluids were not unaltered seawater. We develop a model in which the barites precipitated both below the sediment–water interface and at the seafloor from low-temperature fluids that circulated along faults. The isotopic, chemical, and textural data indicate that the barites formed by several processes. Mesozoic and Cenozoic basement rocks (crystalline and overlying sedimentary rocks), Quaternary basin fill, turbidite fans, and seawater provided various elements for the barites in different environments. The fluids had a deep-seated source and were hydrothermal in the deeper parts of the system for all the barite types, including the seafloor cold-seep deposit, based on Sr isotopes and regional geothermal gradients. These deep-seated fluids mixed with other fluids as they ascended, including seawater at and near the seafloor. The high δ34S values may have resulted from extreme Rayleigh fractionation during bacterially mediated (semi)closed-system sulfate reduction, probably driven by the flux of methane- and hydrocarbon-bearing fluids through basement rocks and basin sediments. Early diagenetic dolomite and pyrite in the host Monterey Formation in the PV Headland also formed in a zone of sulfate reduction, but within sediment buried only centimeters to a few meters and with a predominantly seawater source for the sulfur. Dolomite occurring with vein barite in the PV Headland formed at temperatures in the range of 40–90 °C. The cold-seep barites have δ34S values near that of modern seawater, although still somewhat fractionated. The barites that precipitated below the sediment–water interface have higher δ34S values, suggesting that the fluids were relatively reduced with molar dissolved barium in excess of dissolved sulfate. Those samples were exposed at the seafloor by uplift along faults and are composed predominantly of massive, brecciated, and vein barite.

California↗
USGS2020

History and evolution of seepage meters for quantifying flow between groundwater and surface water: Part 1 – Freshwater settings

More than 75 years after its introduction, the seepage meter remains the only device for directly quantifying exchange across the sediment-water interface between groundwater and surface water. This device, first presented in the literature in the 1940s, has been in a state of near-constant improvement and design change, necessitating a review of the history and evolution of the device and a description of current best-measurement practices. Part 1 of this two-part review documents the evolution of seepage meters deployed in freshwater settings, including a listing of suggestions for best-measurement and deployment practices. Part 2 covers the same scope for seepage meters deployed in marine settings. Traditional seepage meters isolate a portion of the sediment bed; seepage commonly is determined by routing the volume of flow across that isolated interface to or from a submerged measurement bag over a known time interval. The time-integrated volume is then divided by the bed area covered by the meter to obtain a seepage flux expressed in distance per time. Both the instrument and the measurement are deceptively simple, leading some early users to question the viability of the measurement. Numerous sources of error have been identified and addressed over the decades, resulting in large improvements in measurement consistency and accuracy. Duration of each measurement depends on the seepage rate and can vary from minutes to days, leading to the erroneous and yet common assumption that seepage is relatively stable over time. Designs that replace the measurement bag with a flowmeter eliminate bag-related errors and provide much finer temporal resolution. Resulting data indicate seepage is highly variable in many settings and responds to numerous sub-daily processes, including evapotranspiration, rainfall, seiches and waves. Combining direct measurements from seepage meters with other measurements, such as vertical hydraulic gradients and vertical temperature profiles, provides far better understanding of the processes that control exchange between groundwater and surface water.

Earth-Science Reviews↗
USGS2014

Triggered aseismic slip adjacent to the 6 February 2013 Mw 8.0 Santa Cruz Islands megathrust earthquake

Aseismic or slow slip events have been observed in many subduction zones, but whether they affect the occurrence of earthquakes or result from stress changes caused by nearby events is unclear. In an area lacking direct geodetic observations, inferences can be made from seismological studies of co-seismic slip, associated stress changes and the spatiotemporal nature of aftershocks. These observations indicate that the February 2013 Mw 8.0 Santa Cruz Islands earthquake may have triggered slow or aseismic slip on an adjacent section of the subduction thrust over the following hours to days. This aseismic event was equivalent to Mw 7.6, significantly larger than any earthquakes in the aftershock sequence. The aseismic slip was situated within the seismogenic portion of the subduction interface, and must have occurred to the south of the main seismic slip and most aftershocks in order to promote right-lateral faulting in the upper plate, the dominant deformation style of the aftershock sequence. This plate boundary segment can support either stable sliding (aseismic) or stick-slip (seismic) deformation in response to different driving conditions. The complete lack of aftershocks on the thrust interface implies this pair of megathrust slip episodes (seismic and aseismic) released a substantial portion of the stored strain on the northernmost section of the Vanuatu subduction zone.

Santa Cruz Islands↗
USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗
USGS2008

Weathering of the Rio Blanco quartz diorite, Luquillo Mountains, Puerto Rico: Coupling oxidation, dissolution, and fracturing

In the mountainous Rio Icacos watershed in northeastern Puerto Rico, quartz diorite bedrock weathers spheroidally, producing a 0.2-2 m thick zone of partially weathered rock layers (???2.5 cm thickness each) called rindlets, which form concentric layers around corestones. Spheroidal fracturing has been modeled to occur when a weathering reaction with a positive ??V of reaction builds up elastic strain energy. The rates of spheroidal fracturing and saprolite formation are therefore controlled by the rate of the weathering reaction. Chemical, petrographic, and spectroscopic evidence demonstrates that biotite oxidation is the most likely fracture-inducing reaction. This reaction occurs with an expansion in d (0 0 1) from 10.0 to 10.5 A??, forming 'altered biotite'. Progressive biotite oxidation across the rindlet zone was inferred from thin sections and gradients in K and Fe(II). Using the gradient in Fe(II) and constraints based on cosmogenic age dates, we calculated a biotite oxidation reaction rate of 8.2 ?? 10-14 mol biotite m-2 s-1. Biotite oxidation was documented within the bedrock corestone by synchrotron X-ray microprobe fluorescence imaging and XANES. X-ray microprobe images of Fe(II) and Fe(III) at 2 ??m resolution revealed that oxidized zones within individual biotite crystals are the first evidence of alteration of the otherwise unaltered corestone. Fluids entering along fractures lead to the dissolution of plagioclase within the rindlet zone. Within 7 cm surrounding the rindlet-saprolite interface, hornblende dissolves to completion at a rate of 6.3 ?? 10-13 mol hornblende m-2 s-1: the fastest reported rate of hornblende weathering in the field. This rate is consistent with laboratory-derived hornblende dissolution rates. By revealing the coupling of these mineral weathering reactions to fracturing and porosity formation we are able to describe the process by which the quartz diorite bedrock disaggregates and forms saprolite. In the corestone, biotite oxidation induces spheroidal fracturing, facilitating the influx of fluids that react with other minerals, dissolving plagioclase and chlorite, creating additional porosity, and eventually dissolving hornblende and precipitating secondary minerals. The thickness of the resultant saprolite is maintained at steady state by a positive feedback between the denudation rate and the weathering advance rate driven by the concentration of pore water O2 at the bedrock-saprolite interface. ?? 2008 Elsevier Ltd. All rights reserved.

Geochimica et Cosmochimica Acta↗
USGS2017

Lithological influences on contemporary and long-term regolith weathering at the Luquillo Critical Zone Observatory

Lithologic differences give rise to the differential weatherability of the Earth’s surface and globally variable silicate weathering fluxes, which provide an important negative feedback on climate over geologic timescales. To isolate the influence of lithology on weathering rates and mechanisms, we compare two nearby catchments in the Luquillo Critical Zone Observatory in Puerto Rico, which have similar climate history, relief and vegetation, but differ in bedrock lithology. Regolith and pore water samples with depth were collected from two ridgetops and at three sites along a slope transect in the volcaniclastic Bisley catchment and compared to existing data from the granitic Río Icacos catchment. The depth variations of solid-state and pore water chemistry and quantitative mineralogy were used to calculate mass transfer (tau) and weathering solute profiles, which in turn were used to determine weathering mechanisms and to estimate weathering rates. Regolith formed on both lithologies is highly leached of most labile elements, although Mg and K are less depleted in the granitic than in the volcaniclastic profiles, reflecting residual biotite in the granitic regolith not present in the volcaniclastics. Profiles of both lithologies that terminate at bedrock corestones are less weathered at depth, near the rock-regolith interfaces. Mg fluxes in the volcaniclastics derive primarily from dissolution of chlorite near the rock-regolith interface and from dissolution of illite and secondary phases in the upper regolith, whereas in the granitic profile, Mg and K fluxes derive from biotite dissolution. Long-term mineral dissolution rates and weathering fluxes were determined by integrating mass losses over the thickness of solid-state weathering fronts, and are therefore averages over the timescale of regolith development. Resulting long-term dissolution rates for minerals in the volcaniclastic regolith include chlorite: 8.9 × 10 −14 mol m −2 s −1 , illite: 2.1 × 10 −14 mol m −2 s −1 and kaolinite: 4.0 × 10 −14 mol m −2 s −1 . Long-term weathering fluxes are several orders of magnitude lower in the granitic regolith than in the volcaniclastic, despite higher abundances of several elements in the granitic regolith. Contemporary weathering fluxes were determined from net (rain-corrected) solute profiles and thus represent rates over the residence time of water in the regolith. Contemporary weathering fluxes within the granitic regolith are similar to the long-term fluxes. In contrast, the long-term fluxes are faster than the contemporary fluxes in the volcaniclastic regolith. Contemporary fluxes in the granitic regolith are generally also slightly faster than in the volcaniclastic. The differences in weathering fluxes over space and time between these two watersheds indicate significant lithologic control of chemical weathering mechanisms and rates.

Bisley watersheds↗
USGS2006

Ge/Si and 87Sr/86Sr tracers of weathering reactions and hydrologic pathways in a tropical granitoid system

Ge/Si and 87Sr/86Sr data from primary and secondary minerals, soil waters, and stream waters in a tropical granitoid catchment quantitatively reflect mineral alteration reactions that occur at different levels within the bedrock-saprolite-soil zone. Near the bedrock-saprolite interface, plagioclase to kaolinite reaction yields low Ge/Si and 87Sr/86Sr. Higher in the regolith column, biotite weathering and kaolinite dissolution drive Ge/Si and 87Sr/86Sr to high values. Data from streams at base flow sample the bedrock-saprolite interface zone, while at high discharge solutes are derived from upper saprolite-soil zone. Coupled Ge/Si and 87Sr/86Sr can be effective tools for quantifying the importance of specific weathering reactions, and for geochemical hydrograph separation. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Geochemical Exploration↗