Search NASA⌕ Search

SEARCH · Search NASA

Results for “large rivers”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Yana ring structure, North-Eastern Siberia: A possible counterpart of coronae on Venus

An arch of Verkhoyansky Range in NE Siberia, Russia is considered as W and S parts of the rim of the large ring structure, referred hereafter as Yana Ring Structure (YRS). It has general appearance similar to coronae on Venus and could be considered as the terrestrial counter part of coronae. The structure is located between 63 and 70 deg N, 125 and 140 deg E. The outer diameter of YRS (between the rim feet) is 680-700 km. The rim crest diameter is 500-520 km. The inner diameters (between the inner feet of the rim) is 320-350 km. The rim of YRS consist of the following topographic features: W segment is Orulgan Range (N part of Verkhoyansky Range), SW and S segments consist of S part of Verkhoyansky Range, SE segment is Elginskoye Ploskogorye (Tableland), E segment is Chersky Range (its part W of Indigirka River), NE segment is Selennyakhsky Range, and NW segment is Kular Range. Inside this ring of mountains Yanskoye Ploskogorye (Tableland) is located. The circular outlines of YRS are marked with the large river valleys located outside: Lena River (lower of Aldan River mouth) is located W from YRS, Aldan River (lower of Amga River mouth) is from S, Indigirka River (between 64 and 67 deg N lat) is from E, and Selennyakh River (the part running N to S) is from NE. The ring patterns inside the mountaineous rim of YRS are marked with the tributaries of Yana River. The largest ring here, 300 km in diameter, consists of Bytangay (with Bellyakh) River at W, Nelgese and Derbeke Rivers at S/SE, and Tuostakh River at E. The typical altitudes of YRS mountain rim are 1000-2000 m. The maximal ones are: 2400 m in W part of the rim, 2300 m in S, 1600 m in SE, 2200-2500 m in E, and 1300-1500 m in N. The central part of YRS have typical altitudes of 400-600 m. There are lowlands with altitudes 50-200 m to the S, W, N, and NW outside YRS Oymyakonskoye Nagorye (Highland) is located SE of YRS at 500-1500 m altitudes. It should be recognized from this description YRS have a general appearance of circular highland with its central plateau located 1000-2000 m lower than surrounding mountain rim, but 300-400 m higher than plains of lowlands located outside the rim. Such topographic patterns are typical for coronae on Venus. The diameter of YRS is equal to the largest coronae on Venus. There are only 3 or 4 structures of such or larger size on the whole Venus. This calls to compare morphology and geologic evolution of YRS with morphology and models of coronae evolutions especially the largest ones, to search the similarities and/or differences in geologic histories of the structures on the Earth and Venus.

Burba, G. A.↗
DOE OSTI2024

Zirconium Sludge Criticality Calculations in Large Process Tanks

The Savannah River Site’s H-Canyon facility has been tasked with accelerating the disposition of used research reactor fuels. • Future missions will involve dissolution of “Non-Aluminum” clad fuels – not easily dissolved – Nitric acid is not capable of dissolving stainless steel, Zirconium, and Hastelloy, clad fuels by itself • The electrolytic dissolver has been (re) selected as a disposition path for these fuels. 2

Devine, Nathan P.↗
DOE OSTI2024

Zirconium Sludge Criticality Calculations in Large Process Tanks

The Savannah River Site’s (SRS) H-Canyon recently re-activated its 6.3D electrolytic dissolver to dissolve stainless steel clad research reactor fuels from Japan. A potential new mission was identified to dissolve other long-cooled fuel stored at SRS that is not aluminum clad. In the basic process of dissolution, the spent fuel dissolves into nitric acid (bulk solution). The eight ft. high and four ft. radius 6.3D dissolver is unique in that it has a platinum-coated niobium insert to resist corrosion, and the niobium basket forms an anode and cathode arrangement within the charge chute (Fig 1,2). DC current flows from anode to cathode, at up to 10,000 amps at 45 volts, through the potential gradient in the nitric acid solution, which allows the stainless steel (or other non-aluminum cladding) of the fuel to dissolve.

Wade, Brindley↗
NASA NTRS1994

Long river profiles, tectonism, and eustasy: A guide to interpreting fluvial terraces

Along three rivers at the Mendocino triple junction, northern California, strath, cut, and fill terraces have formed in response to tectonic and eustatic processes. Detailed surveying and radiometric dating at multiple sites indicate that lower reaches of the rivers are dominated by the effects of oscillating sea level, primarily aggradation and formation of fill terraces during sea level high stands, alternating with deep incision during low stands. A eustasy-driven depositional wedge extends tens of kilometers upstream on all rivers (tapering to zero thickness). This distance is greater than expected from studies of the effects of check dams on much smaller streams elsewhere, due in part to the large size of these rivers. However, the change in gradient is nearly identical to other base level rise studies: the depositional gradient is about half that of the original channel. Middle to upper reaches of each river are dominated by the effects of long-term uplift, primarily lateral and vertical erosion and formation of steep, unpaired strath terraces exposed only upstream of the depositional wedge. Vertical incision at a rate similar to that of uplift has occurred even during the present sea level high stand along rivers with highest uplift rates. Strath terraces have steeper gradients than the modern channel bed and do not merge with marine terraces at the river mouth; consequently, they cannot be used to determine altitudes of sea level high stands. Strath formation is a continuous process of response to long-term uplift, and its occurrence varies spatially along a river depending on stream power, and hence position, upstream. Strath terraces are found only along certain parts of a coastal stream: upstream of the aggradational effects of oscillating sea level, and far enough downstream that stream power is in excess of that needed to transport the prevailing sediment load. For a given size river, the greater the uplift rate, the greater the rate of vertical incision and, consequently, the less the likelihood of strath terrace formation and preservation.

Merritts, Dorothy J.↗
DOE OSTI2020

Characterizing the diverse hydrogeology underlying rivers and estuaries using new floating transient electromagnetic methodology

We report the hydrogeology below large surface water features such as rivers and estuaries is universally under-informed at the long reach to basin scales (tens of km+). This challenge inhibits the accurate modeling of fresh/saline groundwater interfaces and groundwater/surface water exchange patterns at management-relevant spatial extents. Here we introduce a towed, floating transient electromagnetic (TEM) system (i.e. FloaTEM) for rapid (up to 15 km/h) high resolution electrical mapping of the subsurface below large water bodies to depths often a factor of 10 greater than other towed instruments. The novel FloaTEM system is demonstrated at a range of diverse 4th through 6th-order riverine settings across the United States including 1) the Farmington River, near Hartford, Connecticut; 2) the Upper Delaware River near Barryville, New York; 3) the Tallahatchie River near Shellmound, Mississippi; and, 4) the Eel River estuary, on Cape Cod, near Falmouth, Massachusetts. Airborne frequency-domain electromagnetic and land-based towed TEM data are also compared at the Tallahatchie River site, and streambed geologic scenarios are explored with forward modeling. A range of geologic structures and pore water salinity interfaces were identified. Process-based interpretation of the case study data indicated FloaTEM can resolve varied sediment-water interface materials, such as the accumulation of fines at the bottom of a reservoir and permeable sand/gravel riverbed sediments that focus groundwater discharge. Bedrock layers were mapped at several sites, and aquifer confining units were defined at comparable resolution to airborne methods. Terrestrial fresh groundwater discharge with flowpaths extending hundreds of meters from shore was also imaged below the Eel River estuary, improving on previous hydrogeological characterizations of that nutrient-rich coastal exchange zone. In summary, the novel FloaTEM system fills a critical gap in our ability to characterize the hydrogeology below surface water features and will support more accurate prediction of groundwater/surface water exchange dynamics and fresh-saline groundwater interfaces.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Ventifact Formation in the Mojave Desert: Field Analogs for Martian Processes

Two field studies in the Mojave Desert, California, shed light on processes of ventifact formation. The field sites are located on a ridge at Little Cowhole Mountain, which lies approximately 12 km south of Baker, and on an unnamed ridge situated along the northern boundary of the Mojave River Sink (Rasor Road site). The rocks at Little Cowhole Mountain are a blue-grey marble/dolomite, whereas those at Rasor Road are Miocene volcanic rocks (basalt). At both sites the abrasive agent is a fine-grained aeolian sand which was probably derived largely from the Mojave River. There are minimal modem inputs of sand to either site: abrasion occurs as a result of unique climatic and topographic conditions which allow pre-existing sand to be recycled from one aspect of the ridge to the other. Climatic conditions are well suited for ventifact formation. Owing, to the dry climate (marked by low average relative humidity, infrequent dew, and low annual rainfall), rates of chemical weathering are low. Where resurfacing of the rocks by sand abrasion proceeds at a rate greater than weathering, the ventifacts are considered "active." Active ventifacts are found atop and straddling the ridge crests, in the zone of maximum wind velocity and sediment supply. Inactive ventifacts occur where modem weathering Processes exceed abrasion rates; principally on the basal two-thirds of the hillslope, where wind velocity and sediment supply are lower. At intermediate locations between the slope base and crest, ventifacts are either active or inactive, depending on local conditions. The presence of relict ventifacts at the study sites, as well as elsewhere in the eastern Mojave Desert, suggests that the conditions for venti fact formation must have been more intense and extensive in the past. Together, the studies illustrate that the processes that interact to form ventifacts are highly complex, and must be studied at many scales. Small-scale effects, such as local topography, plant cover, or even the spatial distribution of boulders, strongly influence the formation of each individual ventifact. Mesoscale effects (hillslope form, total availability of sediment, seasonality of winds, etc.) determine the distribution of overall erosional energy and the location of active and relict forms.

Laity, J. E.↗
NASA NTRS2011

Changes in the Carbon Cycle of Amazon Ecosystems During the 2010 Drought

Satellite remote sensing was combined with the NASA-CASA carbon cycle simulation model to evaluate the impact of the 2010 drought (July through September) throughout tropical South America. Results indicated that net primary production (NPP) in Amazon forest areas declined by an average of 7% in 2010 compared to 2008. This represented a loss of vegetation CO2 uptake and potential Amazon rainforest growth of nearly 0.5 Pg C in 2010. The largest overall decline in ecosystem carbon gains by land cover type was predicted for closed broadleaf forest areas of the Amazon River basin, including a large fraction of regularly flooded forest areas. Model results support the hypothesis that soil and dead wood carbon decomposition fluxes of CO2 to the atmosphere were elevated during the drought period of 2010 in periodically flooded forest areas, compared to forests outside the main river floodplains.

Potter, Christophera↗
DOE OSTI2025

A Comparative Study of Physics‐Informed and Data‐Driven Neural Networks for Compound Flood Simulation at River‐Ocean Interfaces: A Case Study of Hurricane Irene

Simulating compound flooding (CF) at the river-ocean interface within large-scale Earth System Models (ESMs) presents significant challenges due to complex interactions between river discharge, storm surge, and tides. This study assesses the comparative advantages of physics-informed and data-driven machine learning (ML) approaches for enhancing local ESM performance. We systematically compare data-driven neural network models (i.e., CNNs, U-Net, Long Short-Term Memory (LSTM), Gated Recurrent Unit), and physics-informed neural network (PINN) models, including vanilla PINN and a finite-difference-based PINN (FD-PINN). Specifically, FD-PINN is introduced to enhance computational efficiency, accelerating vanilla PINNs by ∼6.5 times while improving accuracy. To enhance data-driven model training, a new data-generation approach is developed to sample historical fluvial and coastal flood events, which ensures a robust data set for extreme event prediction. The models are evaluated using a realistic one-dimensional river domain extracted from an ESM's river mesh and the Hurricane Irene event as an independent test case. Results show that FD-PINN achieves accurate predictions with significantly reduced computational costs relative to vanilla PINNs. Among data-driven models, the best overall performance is achieved by a CNN-LSTM hybrid, which balances accuracy and efficiency. While a fully connected CNN (CNN-FC) provides the best accuracy, it incurs high computational cost. Architectures lacking strong temporal modeling tend to underperform on unseen events. These findings highlight the importance of sequence-aware designs for robust generalization. This study reveals the trade-offs between physics-informed and data-driven models and proposes an adaptive hybrid framework for integrating ML into ESMs to enhance local flood simulations.

Earth Systems Modeling↗
DOE OSTI2024

Increasing phosphorus loss despite widespread concentration decline in US rivers

The loss of phosphorous (P) from the land to aquatic systems has polluted waters and threatened food production worldwide. Systematic trend analysis of P, a nonrenewable resource, has been challenging, primarily due to sparse and inconsistent historical data. Here, we leveraged intensive hydrometeorological data and the recent renaissance of deep learning approaches to fill data gaps and reconstruct temporal trends. We trained a multitask long short-term memory model for total P (TP) using data from 430 rivers across the contiguous United States (CONUS). Trend analysis of reconstructed daily records (1980–2019) shows widespread decline in concentrations, with declining, increasing, and insignificantly changing trends in 60%, 28%, and 12% of the rivers, respectively. Concentrations in urban rivers have declined the most despite rising urban population in the past decades; concentrations in agricultural rivers however have mostly increased, suggesting not-as-effective controls of nonpoint sources in agriculture lands compared to point sources in cities. TP loss, calculated as fluxes by multiplying concentration and discharge, however exhibited an overall increasing rate of 6.5% per decade at the CONUS scale over the past 40 y, largely due to increasing river discharge. Results highlight the challenge of reducing TP loss that is complicated by changing river discharge in a warming climate.

Science & Technology - Other Topics↗
DOE OSTI2021

A novel construct for scaling groundwater–river interactions based on machine-guided hydromorphic classification

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. Predictive numerical models are needed to understand responses of river corridors to environmental change and to support sustainable watershed management. We posit that systematic hydromorphic classification provides a scaling construct that facilitates extrapolation of outputs from local-scale mechanistic models to reduced-order models applicable at reach and watershed scales. This in turn offers the potential to improve large-scale predictions of river corridor hydrobiogeochemical processes. Here we present a new machine-guided hydromorphic classification methodology that addresses the key requirements of this objective, and we demonstrate its application to a segment of the Columbia River in the northwestern United States. The resulting hydromorphic classes form spatially coherent and physically interpretable hydromorphic units that exhibit distinct behaviors in terms of distributions of subsurface residence times (a primary control on critical biogeochemical reactions). This approach forms the basis of ongoing research that is evaluating the formulation of reduced-order models and transferability of results to other river reaches and larger scales.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Optical Characterisation of Suspended Particles in the Mackenzie River Plume (Canadian Arctic Ocean) and Implications for Ocean Colour Remote Sensing

Climate change significantly impacts Arctic shelf regions in terms of air temperature, ultraviolet radiation, melting of sea ice, precipitation, thawing of permafrost and coastal erosion. Direct consequences have been observed on the increasing Arctic river flow and a large amount of organic carbon sequestered in soils at high latitudes since the last glacial maximum can be expected to be delivered to the Arctic Ocean during the coming decade. Monitoring the fluxes and fate of this terrigenous organic carbon is problematic in such sparsely populated regions unless remote sensing techniques can be developed and proved to be operational. The main objective of this study is to develop an ocean colour algorithm to operationally monitor dynamics of suspended particulate matter (SPM) on the Mackenzie River continental shelf (Canadian Arctic Ocean) using satellite imagery. The water optical properties are documented across the study area and related to concentrations of SPM and particulate organic carbon (POC). Robust SPM and POC : SPM proxies are identified, such as the light backscattering and attenuation coefficients, and relationships are established between these optical and biogeochemical parameters. Following a semi-analytical approach, a regional SPM quantification relationship is obtained for the inversion of the water reflectance signal into SPM concentration. This relationship is reproduced based on independent field optical measurements. It is successfully applied to a selection of MODIS satellite data which allow estimating fluxes at the river mouth and monitoring the extension and dynamics of the Mackenzie River surface plume in 2009, 2010 and 2011. Good agreement is obtained with field observations representative of the whole water column in the river delta zone where terrigenous SPM is mainly constrained (out of short periods of maximum river outflow). Most of the seaward export of SPM is observed to occur within the west side of the river mouth. Future work will require the validation of the developed SPM regional algorithm based on match-ups with field measurements, then the routine application to ocean colour satellite data in order to better estimate the fluxes and fate of SPM and POC delivered by the Mackenzie River to the Arctic Ocean.

Mackenzie River plume↗
DOE OSTI2025

Rescue of Unique Liquid IP-2 Portable Tanks at Savannah River Site

Introduction • Savanah River Site procured six large IP-2 portable tanks capable of transporting liquids in 2003 to support layup of F Canyon at SRS. • Tanks are being rescued and recertified to support eventual disposition of up to 171,000 L of contaminated solvent solutions from H Canyon at SRS and to potentially transport Depleted Uranium solutions to support the ABD Program.

Fee, Nathaniel C. [Savannah River Nuclear Solution↗
DOE OSTI2021

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗
NASA NTRS1986

Paleodrainages of the Eastern Sahara - The radar rivers revisited (SIR-A/B implications for a mid-tertiary Trans-African drainage system)

The images obtained by the Shuttle Imaging Radar (SIR)-A and -B systems over the southwestern Egypt and northwestern Sudan were coregistered with the Landsat images and the existing maps to aid in extrapolations of the buried paleodrainages ('radar rivers'), first discovered by SIR-A. Field observations explain the radar responses of three types of radar rivers, RR-1 (broad, aggraded valleys filled with alluvium), RR-2 (braided channels inset in the RR-1 valleys), and RR-3 (narrow, long, bedrock-incised channels). A generalized model of the radar rivers, based on field studies and regional geologic relations, shows inferred changes in river regimen since the large valleys were established during the later Paleogene-early Neogene. It is suggested that a former Trans-African master stream system may have flowed from headwaters in the Red Sea Hills southwestward across North Africa, discharging into the Atlantic at the Paleo-Niger delta, prior to the Neogene domal uplifts and building of volcanic edifices across the paths of these ancient watercourses.

Mccauley, J. F.↗
NASA NTRS2009

"The Waters of Meridiani" - Further Support for a Fluvial Interpretation of the Ridged, Layered Units

A relatively unknown terrestrial fluvial environment, the mesoscale megafan, provides analogs for various Martian landscapes, including the etched unit (etched unit, Unite E of Arvidson et al., 2003; ridge-forming unit R of Edgett, 2005) of the Sinus Meridiani region on Mars. A global survey of Earth shows that megafans are very large partial cones of dominantly fluvial sediment with radii on the order of hundreds of km, and very low slopes. Responsible fluvial processes are sufficiently different from those of classical arid alluvial fans and deltas that it is useful to class megafans as separate features. The megafan model calls into question two commonly held ideas. 1. Earth examples prove that topographic basins per se are unnecessary for the accumulation of large sedimentary bodies. 2. River channels are by no means restricted to valleys (Meridiani sediments are termed a "valley-ed volume" of Edgett). These perspectives reveal unexpected parallels with features at Meridiani-several channel-like features that are widespread, mostly as ridges inverted by eolian erosion; channel networks covering thousands of sq km, especially on intercrater plains; and regional relationships of sediment bodies situated immediately downstream of highland masses. These all suggest that fluvial explanations are at least part of the Meridiani story.

Wilkinson, Justin↗
DOE OSTI2020

Assessment of the Corrosion Risks to a Double Shell Tank During the Selective Removal of a Cesium-rich Supernate Layer - 20227

The Direct Feed Low Activity Waste (DF-LAW) activity will provide the initial feed to the Waste Treatment and Immobilization Plant at the Hanford Tank Farm facility. Feed preparation for DF-LAW will be carried out in multiple double shell waste tanks (DST) utilizing the Tank Side Cesium Removal (TSCR) process. Tank AP-106 will be an integral part of the TSCR process and hence key to the DFLAW activity. Prior to being utilized for this purpose, the Cs-137 concentration levels in the AP-106 supernate must be reduced significantly. The AP-106 re-purposing process manages normal waste tank addition procedures such that waste acceptance criteria objectives for Cs-137 removal are achieved. As this process was being developed, several concerns regarding tank integrity during and after repurposing of Tank AP-106, and the associated tanks (i.e., AP-102, AN-101, and AP-108), were addressed. The potential corrosion degradation mechanisms identified based on process knowledge included: general and pitting corrosion, caustic stress corrosion cracking, concentration cell corrosion, and microbiologically induced corrosion (MIC). A systematic approach was undertaken to assess the risk of corrosion during the AP-106 re-purposing process. The assessment considered the initial and present condition of the tank, the service history of the tank, the environmental conditions created in each tank during each stage of the process, process timing, and corrosion testing and modeling that approximated tank conditions. Based on this information, recommendations were provided to minimize the exposure of the tank to corrosive conditions. These recommendations included time limits for each stage of the process, staging of corrosion inhibitor solutions should process delays occur, and waste sampling with potential chemistry adjustments at defined steps in the process. Additionally, recommendations for post-process supernate chemistry adjustments and tank inspections were put forth. The greatest risk for this process occurs during the manipulation of large quantities of Columbia River water (i.e., greater than 450,000 gallons) within AP-106, the associated tanks and the supporting infrastructure (i.e., transfer line piping). This risk becomes more pronounced if there is a process upset and the tank walls remain exposed to a corrosive condition longer than intended. As long as the river water is blended into the supernate in a timely fashion, and the waste is adjusted appropriately with corrosion inhibitors, corrosion may be minimized and the final condition of the supernate in AP-106 and the associated tanks is expected to be in compliance with corrosion chemistry requirements. Contingency plans will be made to minimize the impact of a process upset. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Improving Elemental Mercury Recovery in DWPF by pH and Ionic Strength Studies Using Laboratory Scale Replica – 26141

In order to safely disposition nuclear material at large scales, the Savannah River Site (SRS) constructed the Salt Waste Processing Facility (SWPF) for the removal of actinides, the Saltstone facility for preparation of a low-level cementitious grout, and the Defense Waste Processing Facility (DWPF) for the vitrification of high-level waste (HLW). A few years after processing began, high concentrations of mercury were discovered in cementitious waste at the Saltstone facility which initiated a search for the source of mercury during processing at these SRS facilities. Mercury serves as a catalyst in the dissolution of spent nuclear fuel, aluminum-actinide alloys, for actinide recovery, but not enough is known about the behavior and properties of mercury within these process streams.

Pina, Jeanette [Savannah River National Laboratory↗
DOE OSTI2023

Linking Stream Chemistry to Subsurface Redox Architecture

As drinking-water scarcity grows worldwide, we need to improve predictions of the quantity and quality of our water resources. An overarching problem for model improvement is that we do not know the geological structure of aquifers in sufficient detail. In this work, we demonstrate that mineral-water reactions imprint structure in the subsurface that impacts the flow and transport of some chemical species. Specifically, pyrite, a ubiquitous mineral, commonly oxidizes and depletes in the upper layers of the weathering profile in most humid watersheds, only remaining at depths of meters. We hypothesize that variations in concentrations (C) of pyrite-derived sulfate released into rivers as a function of discharge (q) reflect the rate-limiting step and depth of pyrite-oxidizing layers. We found that log C - log q behaviors thus differ in small and large watersheds in the Susquehanna River Basin as well as in selected watersheds in the Western United States. Although coal mining changes pyrite oxidation from closed to open system with respect to O 2 , patterns in stream chemistry as a function of discharge are consistent with deep and shallow pyrite oxidation zones in small and large watersheds respectively. Therefore, understanding the subsurface patterns of mineral reactions and how they affect the architecture of aquifers will elucidate patterns of changing river chemistry and our ability to manage water resources in the future under accelerated land use and climate change.

54 ENVIRONMENTAL SCIENCES↗