Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Soil Research”

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 865 records · Page 48Linked to original sources

Geochemical evidence for a magmatic CO 2 degassing event at Mammoth Mountain, California, September-December 1997

Recent time series soil CO 2 concentration data from monitoring stations in the vicinity of Mammoth Mountain, California, reveal strong evidence for a magmatic degassing event during the fall of 1997 lasting more than 2 months. Two sensors at Horseshoe Lake first recorded the episode on September 23, 1997, followed 10 days later by a sensor on the north flank of Mammoth Mountain. Direct degassing from shallow intruding magma seems an implausible cause of the degassing event, since the gas released at Horseshoe Lake continued to be cold and barren of other magmatic gases, except for He. We suggest that an increase in compressional strain on the area south of Mammoth Mountain driven by movement of major fault blocks in Long Valley caldera may have triggered an episode of increased degassing by squeezing additional accumulated CO 2 from a shallow gas reservoir to the surface along faults and other structures where it could be detected by the CO 2 monitoring network. Recharge of the gas reservoir by CO 2 emanating from the deep intrusions that probably triggered deep long-period earthquakes may also have contributed to the degassing event. The nature of CO 2 discharge at the soil-air interface is influenced by the porous character of High Sierra soils and by meteorological processes. Solar insolation is the primary source of energy for the Earth atmosphere and plays a significant role in most diurnal processes at the Earth surface. Data from this study suggest that external forcing due largely to local orographic winds influences the fine structure of the recorded CO 2 signals.

California↗

Spatiotemporal analysis of black spruce forest soils and implications for the fate of C

Post-fire storage of carbon (C) in organic-soil horizons was measured in one Canadian and three Alaskan chronosequences in black spruce forests, together spanning stand ages of nearly 200 yrs. We used a simple mass balance model to derive estimates of inputs, losses, and accumulation rates of C on timescales of years to centuries. The model performed well for the surface and total organic soil layers and presented questions for resolving the dynamics of deeper organic soils. C accumulation in all study areas is on the order of 20–40 gC/m 2 / yr for stand ages up to ∼200 yrs. Much larger fluxes, both positive and negative, are detected using incremental changes in soil C stocks and by other studies using eddy covariance methods for CO 2 . This difference suggests that over the course of stand replacement, about 80% of all net primary production (NPP) is returned to the atmosphere within a fire cycle, while about 20% of NPP enters the organic soil layers and becomes available for stabilization or loss via decomposition, leaching, or combustion. Shifts toward more frequent and more severe burning and degradation of deep organic horizons would likely result in an acceleration of the carbon cycle, with greater CO 2 emissions from these systems overall.

Journal of Geophysical Research↗

Biological soil crusts: Diminutive communities of potential global importance

Biological soil crusts (biocrusts) are widespread, diverse communities of cyanobacteria, fungi, lichens, and mosses living on soil surfaces, primarily in drylands. Biocrusts can locally govern primary production, soil fertility, hydrology, and surface energy balance, with considerable variation in these functions across alternate community states. Further, these communities have been implicated in Earth system functioning via potential influences on global biogeochemistry and climate. Biocrusts are easily destroyed by disturbances and appear to be exceptionally vulnerable to warming temperatures and altered precipitation inputs, signaling possible losses of dryland functions with global change. Despite these concerns, we lack sufficient spatiotemporal data on biocrust function, cover, and community structure to confidently assess their ecological roles across the extensive dryland biome. Here, we present the case for cross-scale research and restoration efforts coupled with remote-sensing and modeling approaches that improve our collective understanding of biocrust responses to global change and the ecological roles of these diminutive communities at global scales.

Frontiers in Ecology and the Environment↗

Methods for installation, removal, and downloading data from the temperature profiling probe (TROD)

This document provides a brief method overview on the deployment and removal of the temperature profiling probe developed by the U.S. Geological Survey (USGS) in 2015 and referred to as SensorRod or temperature rod (TROD). The TROD is suitable for short- to long-term deployments (days to years) for evaluating thermal gradients in soils and sediments beneath surface water. Applications include evaluating exchange between groundwater and surface water, seepage losses in agricultural canals, infiltration on hillslopes, and timing of snowmelt and runoff. The temperature sensors inside the TROD are iButtons and are meant to be replaceable when needed. The housing is water tight and is made of schedule 80 PVC pipe. Sensor intervals are flexible and can be built to specific applications. A joint patent (U.S 10/180,360.) and license agreement between the USGS and Alpha Mach for the TROD design was obtained in 2019, and a paper in Water Resource Research describes the design and testing (Naranjo and Turcotte, 2015).

Open-File Report↗

A cellular automata downscaling based 1 km global land use datasets (2010–2100)

Global climate and environmental change studies require detailed land-use and land-cover(LULC) information about the past, present, and future. In this paper, we discuss a methodology for downscaling coarse-resolution (i.e., half-degree) future land use scenarios to finer (i.e., 1 km) resolutions at the global scale using a grid-based spatially explicit cellular automata (CA) model. We account for spatial heterogeneity from topography, climate, soils, and socioeconomic variables. The model uses a global 30 m land cover map (2010) as the base input, a variety of biogeographic and socioeconomic variables, and an empirical analysis to downscale coarse-resolution land use information (specifically urban, crop and pasture). The output of this model offers the most current and finest-scale future LULC dynamics from 2010 to 2100 (with four representative concentration pathway (RCP) scenarios—RCP 2.6, RCP 4.5, RCP 6.0, and RCP 8.5) at a 1 km resolution within a globally consistent framework. The data are freely available for download, and will enable researchers to study the impacts of LULC change at the local scale.

Science Bulletin↗

Monitoring gas and heat emissions at Norris Geyser Basin, Yellowstone National Park, USA based on a combined eddy covariance and Multi-GAS approach

We quantified gas and heat emissions in an acid-sulfate, vapor-dominated area (0.04-km 2 ) of Norris Geyser Basin, located just north of the 0.63 Ma Yellowstone Caldera and near an area of anomalous uplift. From 14 May to 3 October 2016, an eddy covariance system measured half-hourly CO 2 , H 2 O and sensible ( H ) and latent ( LE ) heat fluxes and a Multi-GAS instrument measured (1 Hz frequency) atmospheric H 2 O, CO 2 and H 2 S volumetric mixing ratios. We also measured soil CO 2 fluxes using the accumulation chamber method and temperature profiles on a grid and collected fumarole gas samples for geochemical analysis. Eddy covariance CO 2 fluxes ranged from − 56 to 885 g m − 2 d − 1 . Using wavelet analysis, average daily eddy covariance CO 2 fluxes were locally correlated with average daily environmental parameters on several-day to monthly time scales. Estimates of CO 2 emission rate from the study area ranged from 8.6 t d − 1 based on eddy covariance measurements to 9.8 t d − 1 based on accumulation chamber measurements. Eddy covariance water vapor fluxes ranged from 1178 to 24,600 g m − 2 d − 1 . Nighttime H and LE were considered representative of hydrothermal heat fluxes and ranged from 4 to 183 and 38 to 504 W m − 2 , respectively. The total hydrothermal heat emission rate ( H + LE + radiant) estimated for the study area was 11.6 MW and LE contributed 69% of the output. The mean ± standard deviation of H 2 O, CO 2 and H 2 S mixing ratios measured by the Multi-GAS system were 9.3 ± 3.1 parts per thousand, 467 ± 61 ppmv, and 0.5 ± 0.6 ppmv, respectively, and variations in the gas compositions were strongly correlated with diurnal variations in environmental parameters (wind speed and direction, atmospheric temperature). After removing ambient H 2 O and CO 2 , the observed variations in the Multi-GAS data could be explained by the mixing of relatively H 2 O-CO 2 -H 2 S-rich fumarole gases with CO 2 -rich and H 2 O-H 2 S-poor soil gases. The fumarole H 2 O/CO 2 and CO 2 /H 2 S end member ratios (101.7 and 27.1, respectively, on average) were invariant during the measurement period and fell within the range of values measured in direct fumarole gas samples. The soil gas H 2 O/CO 2 end member ratios (~ 15–30) were variable and low relative to the fumarole end member, likely resulting from water vapor loss during cooling and condensation in the shallow subsurface, whereas the CO 2 /H 2 S end member ratio was high (~ 160), presumably related to transport of CO 2 -dominated soil gas emissions mixed with trace fumarolic emissions to the Multi-GAS station. Nighttime eddy covariance ratios of H 2 O to CO 2 flux were typically between the soil gas and fumarole end member H 2 O/CO 2 ratios defined by Multi-GAS measurements. Overall, the combined eddy covariance and Multi-GAS approach provides a powerful tool for quasi-continuous measurements of gas and heat emissions for improved volcano-hydrothermal monitoring.

Wyoming↗

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗

Assessment of calcium status in Maine forests: Review and future projection

Forest harvesting and acidic deposition can cause substantial decreases in the calcium (Ca) inventory of forest soils if such losses are not replenished through mineral weathering, atmospheric deposition, or fertilization. The net balance between losses and gains defines the forest Ca status. Site-specific studies have measured Ca pools and fluxes in Maine forests, but no synthesis has been published. In this paper, I review the literature on forest Ca and assess the current status and potential future trends. Forest soils in Maine are currently at lesser risk of Ca depletion compared with many forest soils in the central and southeastern United States, because levels of acidic deposition and rates of Ca accumulation in trees are lower in Maine. The rate of Ca accumulation in trees is reduced in Maine as a result of lower growth rates and a higher proportion of conifer trees that require less Ca than hardwoods. However, field-scale biogeochemical studies in Maine and New Hampshire, and regional estimates of harvest removals and soil inventories coupled with low weathering estimates, indicate that Ca depletion is a realistic concern in Maine. The synthesis of site-specific and regional data for Maine in conjunction with the depletion measured directly in surrounding areas indicates that the Ca status of many forest soils in Maine is likely declining. Ca status could decrease further in the future if forest growth rates increase in response to climate trends and recovery from insect-induced mortality and excessive harvesting in recent years. Proposed climate change induced reductions in spruce and fir and increases in hardwoods would also increase the risk of Ca depletion. ?? 2005 NRC.

Canadian Journal of Forest Research↗

Initial soil respiration response to biomass harvesting and green-tree retention in aspen-dominated forests of the Great Lakes region

Contemporary forest management practices are increasingly designed to optimize novel objectives, such as maximizing biomass feedstocks and/or maintaining ecological legacies, but many uncertainties exist regarding how these practices influence forest carbon (C) cycling. We examined the responses of soil respiration (R s ) to biomass harvesting and green-tree retention in an effort to empirically assess their impacts on C cycling. We measured R s and soil microclimatic variables over four growing seasons following implementation of these management practices using a fully replicated, operational-scale experiment in aspen-dominated forests in northern Minnesota. Treatments included three levels of biomass removal within harvested areas: whole-tree harvest (no slash deliberately retained), 20% slash retained, and stem-only harvest (all slash retained), and two levels of green-tree retention: 0.1 ha aggregate or none. The relative amount of biomass removed had a negligible effect on R s in harvested areas, but treatment effects were probably obscured by heterogeneous slash configurations and rapid post-harvest regeneration of aspen in all of the treatments. Discrete measurements of R s and soil temperature within green-tree aggregates were not discernible from surrounding harvested areas or unharvested control stands until the fourth year following harvest, when R s was higher in unharvested controls than in aggregates and harvested stands. Growing season estimates of Rs showed that unharvested control stands had higher R s than both harvested stands and aggregates in the first and third years following harvest. Our results suggest that retention of larger forest aggregates may be necessary to maintain ecosystem-level responses similar to those in unharvested stands. Moreover, they highlight the innate complexity of operational-scale research and suggest that the initial impacts of biomass harvest on R s may be indiscernible from traditional harvest in systems where incidental breakage is high.

Great Lakes↗

Permafrost stores a globally significant amount of mercury

Changing climate in northern regions is causing permafrost to thaw with major implications for the global mercury (Hg) cycle. We estimated Hg in permafrost regions based on in situ measurements of sediment total mercury (STHg), soil organic carbon (SOC), and the Hg to carbon ratio (R HgC ) combined with maps of soil carbon. We measured a median STHg of 43 ± 30 ng Hg g soil −1 and a median R HgC of 1.6 ± 0.9 μg Hg g C −1 , consistent with published results of STHg for tundra soils and 11,000 measurements from 4,926 temperate, nonpermafrost sites in North America and Eurasia. We estimate that the Northern Hemisphere permafrost regions contain 1,656 ± 962 Gg Hg, of which 793 ± 461 Gg Hg is frozen in permafrost. Permafrost soils store nearly twice as much Hg as all other soils, the ocean, and the atmosphere combined, and this Hg is vulnerable to release as permafrost thaws over the next century. Existing estimates greatly underestimate Hg in permafrost soils, indicating a need to reevaluate the role of the Arctic regions in the global Hg cycle.

Geophysical Research Letters↗

Aeolian sediments in paleowetland deposits of the Las Vegas Formation

The Las Vegas Formation (LVF) is a well-characterized sequence of groundwater discharge (GWD) deposits exposed in and around the Las Vegas Valley in southern Nevada. Nearly monolithologic bedrock surrounds the valley, which provides an excellent opportunity to test the hypothesis that GWD deposits include an aeolian component. Mineralogical data indicate that the LVF sediments are dominated by carbonate minerals, similar to the local bedrock, but silicate minerals are also present. The median particle size is ~35 μm, consistent with modern dust in the region, and magnetic properties contrast strongly with local bedrock, implying an extralocal origin. By combining geochemical data from the LVF sediments and modern dust, we found that an average of ~25% of the LVF deposits were introduced by aeolian processes. The remainder consists primarily of authigenic groundwater carbonate as well as minor amounts of alluvial material and soil carbonate. Our data also show that the aeolian sediments accumulated in spring ecosystems in the Las Vegas Valley in a manner that was independent of both time and the specific hydrologic environment. These results have broad implications for investigations of GWD deposits located elsewhere in the southwestern U.S. and worldwide.

Nevada↗

Environmental controls, emergent scaling, and predictions of greenhouse gas (GHG) fluxes in coastal salt marshes

Coastal salt marshes play an important role in mitigating global warming by removing atmospheric carbon at a high rate. We investigated the environmental controls and emergent scaling of major greenhouse gas (GHG) fluxes such as carbon dioxide (CO 2 ) and methane (CH 4 ) in coastal salt marshes by conducting data analytics and empirical modeling. The underlying hypothesis is that the salt marsh GHG fluxes follow emergent scaling relationships with their environmental drivers, leading to parsimonious predictive models. CO 2 and CH 4 fluxes, photosynthetically active radiation (PAR), air and soil temperatures, well water level, soil moisture, and porewater pH and salinity were measured during May–October 2013 from four marshes in Waquoit Bay and adjacent estuaries, MA, USA. The salt marshes exhibited high CO 2 uptake and low CH 4 emission, which did not significantly vary with the nitrogen loading gradient (5–126 kg · ha −1 · year −1 ) among the salt marshes. Soil temperature was the strongest driver of both fluxes, representing 2 and 4–5 times higher influence than PAR and salinity, respectively. Well water level, soil moisture, and pH did not have a predictive control on the GHG fluxes, although both fluxes were significantly higher during high tides than low tides. The results were leveraged to develop emergent power law‐based parsimonious scaling models to accurately predict the salt marsh GHG fluxes from PAR, soil temperature, and salinity (Nash‐Sutcliffe Efficiency = 0.80–0.91). The scaling models are available as a user‐friendly Excel spreadsheet named Coastal Wetland GHG Model to explore scenarios of GHG fluxes in tidal marshes under a changing climate and environment.

Massachusetts↗

The association of geochemical anomalies with a negative gravity anomaly in the Chief Mountain-Soda Creek area, Clear Creek County, Colorado

Geochemical studies in the Chief Mountain Soda Creek area, Clear Creek County, Colo., show that anomalously high amounts of Au, Ag, Pb, Zn, Cd, and Bi in mull ash, and Cu and Hg in soil, correspond to a negative gravity anomaly in the Front Range mineral belt. The correspondence of the geochemical anomalies to the negative gravity anomaly suggests the presence of altered and mineralized rock and a Tertiary intrusive body at depth. This gravity anomaly and other similar gravity anomalies in the Front Range mineral belt merit further investigation for possible ore deposits.

Colorado↗

Select techniques for detecting and quantifying seepage from unlined canals

Canal seepage losses affect the ability of water conveyance structures to maximize efficiency and can be a precursor to canal failure. Identification and quantification of canal seepage out of unlined canals is a complex interaction affected by geology, canal stage, operations, embankment geometry, siltation, animal burrows, structures, and other physical characteristics. Seepage out of unlined canals can be coarsely estimated using a mass balance-type approach (water in minus water out with the difference assumed to be a combination of seepage and evapotranspiration). More sophisticated methods are used in some instances but are typically limited efforts aimed at quantifying seepage in a specific location. Seepage is generally broken out into two categories: diffuse and concentrated (or focused) seepage. Diffuse seepage is where the seepage discharges relatively constant over a given area, whereas concentrated (point discharge source) seepage discharges along preferentially focused areas. Diffuse seepage typically occurs in homogeneous conditions where the amount of water flowing into the subsurface is controlled by soil permeability and canal stage. Conversely, concentrated seepage occurs in areas of heterogeneous conditions where water flows into bedrock fractures, rodent burrows or other pre-existing discrete flow-paths. Concentrated seepage can also develop in the advent of sudden or excessive increases in hydraulic gradient which can lead to heaving, cracking, and development of backward erosion piping flow-paths. Concentrated and diffuse seepage can lead to seeps, in this case, a surface expression of water fed by irrigation water on canal embankment or at distal regions away from the canal. This report focuses on work funded by the Research and Development Office from Fiscal Year 2016 through 2021 and the references provided pertain primarily to those efforts. This report also provides a generalized framework for how and when to investigate seepage out of an unlined canal based on the type of seepage, level of understanding about the seepage locations, geology, and knowledge of the subsurface conditions. The various methods used to locate seeps and quantify canal seepage are discussed in further detail, with references provided for the reader. The following seepage investigation scenarios are discussed within the report: 1. Idealized workflow insensitive to time with highest quality data required 2. General workflow sensitive to time with highest quality data required 3. General workflow insensitive to time with lowest cost items preceding more costly techniques 4. Newly developed concentrated seep(s), concern about consequences (time sensitive) 5. Newly developed or rapidly increasing diffuse seepage, concern about consequences (time sensitive) 6. Existing concentrated seep(s), limited concern about consequences, poor geologic understanding 7. Existing concentrated seep(s), limited concern about consequences, good geologic understanding 8. Existing diffuse seepage, limited concern about consequences, poor geologic understanding 9. Existing diffuse seepage, limited concern about consequences, good geologic understanding A workflow is given for each scenario which details recommended steps and the order in which those steps should be taken to maximize efficiency and data quality. The various seepage investigation techniques and estimated costs are discussed in more detail later in this report. The next step is to take the data collected from the various methods and incorporate them into canal operations models to optimize deliveries. This step could also include the development of 3D seepage models to better understand the larger-scale groundwater-surface water interactions and how they are affected by the water delivery system.

Final Report↗

Reconstructed historical land cover and biophysical parameters for studies of land-atmosphere interactions within the eastern United States

Over the past 350 years, the eastern half of the United States experienced extensive land cover changes. These began with land clearing in the 1600s, continued with widespread deforestation, wetland drainage, and intensive land use by 1920, and then evolved to the present-day landscape of forest regrowth, intensive agriculture, urban expansion, and landscape fragmentation. Such changes alter biophysical properties that are key determinants of land-atmosphere interactions (water, energy, and carbon exchanges). To understand the potential implications of these land use transformations, we developed and analyzed 20-km land cover and biophysical parameter data sets for the eastern United States at 1650, 1850, 1920, and 1992 time slices. Our approach combined potential vegetation, county-level census data, soils data, resource statistics, a Landsat-derived land cover classification, and published historical information on land cover and land use. We reconstructed land use intensity maps for each time slice and characterized the land cover condition. We combined these land use data with a mutually consistent set of biophysical parameter classes, to characterize the historical diversity and distribution of land surface properties. Time series maps of land surface albedo, leaf area index, a deciduousness index, canopy height, surface roughness, and potential saturated soils in 1650, 1850, 1920, and 1992 illustrate the profound effects of land use change on biophysical properties of the land surface. Although much of the eastern forest has returned, the average biophysical parameters for recent landscapes remain markedly different from those of earlier periods. Understanding the consequences of these historical changes will require land-atmosphere interactions modeling experiments.

Journal of Geophysical Research D: Atmospheres↗

Scenarios of climate adaptation potential on protected working lands from management of soils

Management of protected lands may enhance ecosystem services that conservation programs were designed to protect. Practices that build soil organic matter (SOM) on agricultural lands also increase soil water holding capacity, potentially reducing climatic water deficit (CWD), increasing actual evapotranspiration (AET) and increasing groundwater recharge (RCH). We developed nine spatially-explicit land use and conservation scenarios (2001 - 2100) in the LUCAS land use change model to address two questions for California working lands (cropland and rangeland): How does land use change limit opportunities to manage soils for hydrologic climate adaptation benefits? To what extent and where can soil management practices increase climate adaptation on protected working lands? Hydrologic benefits [∑(∆CWD, ∆AET, ∆RCH)] due to soil management were simulated in the Basin Characterization Model (a state-wide water balance model) for two Representative Concentration Pathway 8.5 climate models. LUCAS simulated land conversion and new conservation easements with potential for maximum hydrologic benefits. Climate drove differences in lost potential for water benefits due to urbanization (33.9 - 87.6 m3 x 106) in 2050. Conflict between development pressure and potential hydrologic benefits occurred most in Santa Clara County in the San Francisco Bay Area and Shasta County in Northern Sacramento Valley. Hydrologic benefits on easements were similar in magnitude to losses from development. Water savings from management of California Land Conservation (a.k.a. Williamson) Act contract lands were an order of magnitude greater, totaling over 460 m3 x106 annually in a drier climate by 2050. Few counties provide most benefits because of soil properties, climate and land area protected. The increase in hydrologic benefits varies by agricultural practice and adoption rate, land use type and configuration, and terms of conservation agreements. The effectiveness of programs designed to improve climate adaptation at county to state scales will likely increase by taking this variability into consideration.

California↗

Helium in soil gases of the Roosevelt Hot Springs Known Geothermal Resource Ares, Beaver County, Utah

Soil samples were collected in two parallel traverses across the Dome fault zone of the Roosevelt Hot Springs Known Geothermal Resource Area. The samples were sealed in air-tight aluminum cans, and the soil gas was allowed to equilibrate with the atmospheric air in the cans. Gas from the cans was analyzed by mass spectrometry. Samples collected over faults contained anomalously high concentrations of helium. Samples collected close to a geothermal well 884 m deep contained more helium than samples collected near another geothermal well 1370 m deep.

Utah↗

A unified assessment of hydrological and biogeochemical responses in research watersheds in Eastern Puerto Rico using runoff-concentration relations

An examination of the relation between runoff rate, R , and concentration, C , of twelve major constituents in four small watersheds in eastern Puerto Rico demonstrates a consistent pattern of responses. For solutes that are not substantially bioactive (alkalinity, silica, calcium, magnesium, sodium, and chloride), the log( R )–log( C ) relation is almost linear and can be described as a weighted average of two sources, bedrock weathering and atmospheric deposition. The slope of the relation for each solute depends on the respective source contributions to the total river load. If a solute were strictly derived from bedrock weathering, the slope would be −0.3 to −0.4, whereas if strictly derived from atmospheric deposition, the slope would be approximately −0.1. The bioactive constituents (dissolved organic carbon, nitrate, sulfate, and potassium), which are recycled by plants and concentrated in shallow soil, demonstrate nearly flat or downward-arched log( R )–log( C ) relations. The peak of the arch represents a transition from dominantly soil-matrix flow to near-surface macropore flow, and finally to overland flow. At highest observed R (80 to >90 mm/h), essentially all reactive surfaces have become wetted, and the input rate of C becomes independent of R (log( R )–log( C ) slope of –1). The highest R are tenfold greater than any previous study. Slight clockwise hysteresis for many solutes in the rivers with riparian zones or substantial hyporheic flows indicates that these settings may act as mixing end-members. Particulate constituents (suspended sediment and particulate organic carbon) show slight clockwise hysteresis, indicating mobilization of stored sediment during rising stage.

Aquatic Geochemistry↗