Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrologic models”

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 847 records · Page 47

A Distinct Atmospheric Mode for California Precipitation

The hydrologic cycle in California is strongly influenced by wet-season (November–April) precipitation. Here, we demonstrate the existence of an influential mode of North Pacific atmospheric pressure variability that regulates wet-season precipitation variability over both northern and southern California. This mode, named as the “California precipitation mode” (CPM), is statistically distinct from other well-known modes of pressure variability such as the Pacific-North American pattern. In addition to controlling wet-season mean precipitation, positive days of the CPM coincide with up to 90% of the extreme (>99th percentile) precipitation days and 76% of detected atmospheric rivers (ARs) days, while the negative days correspond with 60% of the dry days. CMIP6 models capture the CPM remarkably well, including its statistical separation from the other well-known modes of pressure variability. The models also reproduce the CPM's strong association with California wet-season precipitation, giving confidence in the models’ dynamics relating to regional hydrologic extremes. However, the models also exhibit biases in regional hydrologic extremes. The CPM is a useful way to understand the origins of those biases and select the more credible models for further analysis: Models with unrealistically strong gradients in the CPM pressure pattern generally oversimulate larger wet extremes and produce excessively long dry intervals in the historical period. Thus the hydrologic biases can be traced to the particular aspects of North Pacific atmospheric dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertainty characterization in a coupled human-natural system: Modeling agricultural adaptation in the Great Lakes Region

The Great Lakes Region's water quality and ecological health are threatened by the export of nutrients from agricultural lands, which causes eutrophication, hypoxia, and destructive algal blooms. The intensification of hydrologic cycles brought about by climate change is expected to exacerbate nutrient loading in the region, and, at the same time, agricultural adaptation to changing conditions is also expected to affect loading through shifting amounts and timing of fertilization. Quantifying these future effects and their interactions necessitates modeling both the human and natural processes as a coupled system, by pairing land use and agricultural management with hydrologic modeling. At the same time, compounding uncertainties arising from the complex interactions in both systems significantly limit our predictive understanding of the region's impacts. This study utilizes the Soil and Water Assessment Tool (SWAT), developed for simulating the impact of various farmer decisions on watershed functions in Western Lake Erie watersheds, and an under-development agent-based model (ABM) for agricultural management decisions. The aim of this study is to use global sensitivity analysis on the coupled ABM and SWAT models to quantify how uncertainty in both models interactively affects nutrient loading. To do so, we will conduct Sobol sensitivity analysis experiments at different levels of coupling assumptions to quantify how various uncertain factors (e.g., soil moisture and crop choice) and their interactions affect our estimates of nutrient loading. The results of this analysis will allow us to quantify how complex interactions and dependencies between both systems amplify the effect of uncertainties. Insights gained from this study will have broader implications for modeling the adaptive co-evolution of human and natural systems under climate change and can inform effective management of nutrient loading in the Great Lakes Region.

Climate Change↗

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗

A Regional, Integrated Monitoring System for the Hydrology of the Pan-Arctic Land Mass

Work under this NASA contract developed a system for monitoring and historical analysis of the major components of the pan-Arctic terrestrial water cycle. It is known as Arctic-RIMS (Regional Integrated Hydrological Monitoring System for the Pan-Arctic Landmass). The system uses products from EOS-era satellites, numerical weather prediction models, station records and other data sets in conjunction with an atmosphere-land surface water budgeting scheme. The intent was to compile operational (at 1-2 month time lags) gridded fields of precipitation (P), evapotranspiration (ET), P-ET, soil moisture, soil freeze/thaw state, active layer thickness, snow extent and its water equivalent, soil water storage, runoff and simulated discharge along with estimates of non-closure in the water budget. Using "baseline" water budgeting schemes in conjunction with atmospheric reanalyses and pre-EOS satellite data, water budget fields were conjunction with atmospheric reanalyses and pre-EOS satellite data, water budget fields were compiled to provide historical time series. The goals as outlined in the original proposal can be summarized as follows: 1) Use EOS data to compile hydrologic products for the pan-Arctic terrestrial regions including snowcover/snow water equivalent (SSM/A MODIS, AMSR) and near-surface freeze/thaw dynamics (Sea Winds on QuikSCAT and ADEOS I4 SSMI and AMSR). 2) Implement Arctic-RIMS to use EOS data streams, allied fields and hydrologic models to produce allied outputs that fully characterize pan-Arctic terrestrial and aerological water budgets. 3) Compile hydrologically-based historical products providing a long-term baseline of spatial and temporal variability in the water cycle.

Serreze, Mark↗

How Moist and Dry Intrusions Control the Local Hydrologic Cycle in Present and Future Climates

Models disagree on the magnitude of the intensification of hydrologic cycle extremes under climate change. Changes in these extremes are expected to scale with the Clausius-Clapeyron relation but may locally diverge due to the uncertain response of the atmospheric general circulation. Recent work has identified that the circulation contribution (the “dynamic” component) is the largest contributor to uncertainty in changes in hydrologic cycle extremes. In this work, we develop an exact budget for changes in wet and dry extremes based on changes in the associated moist and dry intrusions behind those extremes. We decompose changes into four contributors: thermodynamic, low-level convergence, anomalous meridional moisture advection, and zonal wave activity convergence. In analyzing these changes, we show that low-level convergence is the main dynamic contributor to the tropical and subtropical annual response, consistent with weakening of the Hadley circulation, and we show that anomalous meridional moisture advection is the primary dynamic contributor annually for the extratropics, consistent with storm track changes. We achieve a clean separation between dynamic and thermodynamic contributions through an empirical scaling between extremes and their associated wave activities, showing again that shifting storm tracks and expanding subtropics are critical to the dynamic response of the transient hydrologic cycle. We infer that changes in the large-scale circulation are predominantly responsible for the annual, dynamic response of hydrologic cycle extremes and therefore an important contributor to uncertainty in future projections.

Hydrological cycle, hydrological extremes, dynamic↗

Machine Learning Calibration of Groundwater Table Depth in ELM: Impact on Land Surface Hydrology and Land‐Atmosphere Fluxes

Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.

Fang, Yilin [Pacific Northwest National Laboratory↗

Some aspects of hydrology in the GLAS GCM

The response of hydrology in the Laboratory for Atmospheric Sciences General Circulation Model is evaluated. The results show that the distribution of precipitation agrees fairly well with observations of precipitation and that the model and that the model tends to maintain the hydrological balance.

Godbole, R.↗

Evaluation of a Regional Crop Model Implementation for Sub-National Yield Assessments in Kenya

CONTEXT: Cropping system models can be used to both assess regional food security and to monitor and predict agricultural drought. Agriculture in Kenya is extremely important to both the economy and food security of the country. OBJECTIVE: This study evaluated a regional implementation of a widely used crop model, the Decision Support System for Agrotechnology Transfer (DSSAT), within a coupled modeling framework, the Regional Hydrologic Extremes Assessment System (RHEAS), over Kenya. The goal of this study was to assess the ability of RHEAS to simulate the annual variability of maize yields at the county level and evaluate the uncertainty inherent in the model and inputs. METHODS: The RHEAS system implements a stochastic ensemble approach to account for field scale variabilities in crop management practices and underlying soil and weather conditions. Satellite-derived datasets were used to evaluate the land surface component of the system and seasonally disaggregated yield for 5 years was used to assess the performance of the cropping system model. RESULTS AND CONCLUSIONS: The median correlation between RHEAS and satellite-derived soil moisture and evapotranspiration estimates were 0.78, and 0.51, respectively, indicating that the model is able to capture the key drivers of the hydrological budget. Overall, RHEAS simulated yearly yield variations with a median correlation of 0.7 with reported yields, with the best performance in the short rains season. However, across both seasons, the RHEAS model was positively biased on the order of ~1.6 MT/ha. The overall median unbiased RMSE was 0.66 MT/ha. The RHEAS system shows skill at simulating extreme departures in anomalies, and a majority of the time (62.5%) the reported yields fall within the interquartile range of the simulations. SIGNIFICANCE: One of the most important areas of improvement for the next generation of agricultural data and models is to better understand and communicate the inherent uncertainties. This is especially critical in data-limited regions. Here we present a modeling system and its implementation that begins to address these concerns. We demonstrate the ability to simulate broad trends in yields at the county level for sub-annual yields with skills that commensurate previous national/annual level studies.

Crop model↗

Spatial Scaling of Snow Observations and Microwave Emission Modeling During CLPX and Appropriate Satellite Sensor Resolution

Accurate estimates of snow water equivalent and other properties play an important role in weather, natural hazard, and hydrological forecasting and climate modeling over a range of scales in space and time. Remote sensing-derived estimates have traditionally been of the "snapshot" type, but techniques involving models with assimilation are also being explored. In both cases, forward emission models are useful to understand the observed passive microwave signatures and developing retrieval algorithms. However, mismatches between passive microwave sensor resolutions and the scales of processes controlling subpixel heterogeneity can affect the accuracy of the estimates. Improving the spatial resolution of new passive microwave satellite sensors is a major desire in order to (literally) resolve such subpixel heterogeneity, but limited spacecraft and mission resources impose severe constraints and tradeoffs. In order to maximize science return while mitigating risk for a satellite concept, it is essential to understand the scaling behavior of snow in terms of what the sensor sees (brightness temperature) as well as in terms of the actual variability of snow. NASA's Cold Land Processes Experiment-1 (CLPX-1: Colorado, 2002 and 2003) was designed to provide data to measure these scaling behaviors for varying snow conditions in areas with forested, alpine, and meadow/pasture land cover. We will use observations from CLPX-1 ground, airborne, and satellite passive microwave sensors to examine and evaluate the scaling behavior of observed and modeled brightness temperatures and observed and retrieved snow parameters across scales from meters to 10's of kilometers. The conclusions will provide direct examples of the appropriate spatial sampling scales of new sensors for snow remote sensing. The analyses will also illustrate the effects and spatial scales of the underlying phenomena (e.g., land cover) that control subpixel heterogeneity.

Kim, Edward J.↗

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Scientific goals of the Cooperative Multiscale Experiment (CME)

Mesoscale Convective Systems (MCS) form the focus of CME. Recent developments in global climate models, the urgent need to improve the representation of the physics of convection, radiation, the boundary layer, and orography, and the surge of interest in coupling hydrologic, chemistry, and atmospheric models of various scales, have emphasized the need for a broad interdisciplinary and multi-scale approach to understanding and predicting MCS's and their interactions with processes at other scales. The role of mesoscale systems in the large-scale atmospheric circulation, the representation of organized convection and other mesoscale flux sources in terms of bulk properties, and the mutually consistent treatment of water vapor, clouds, radiation, and precipitation, are all key scientific issues concerning which CME will seek to increase understanding. The manner in which convective, mesoscale, and larger scale processes interact to produce and organize MCS's, the moisture cycling properties of MCS's, and the use of coupled cloud/mesoscale models to better understand these processes, are also major objectives of CME. Particular emphasis will be placed on the multi-scale role of MCS's in the hydrological cycle and in the production and transport of chemical trace constituents. The scientific goals of the CME consist of the following: understand how the large and small scales of motion influence the location, structure, intensity, and life cycles of MCS's; understand processes and conditions that determine the relative roles of balanced (slow manifold) and unbalanced (fast manifold) circulations in the dynamics of MCS's throughout their life cycles; assess the predictability of MCS's and improve the quantitative forecasting of precipitation and severe weather events; quantify the upscale feedback of MCS's to the large-scale environment and determine interrelationships between MCS occurrence and variations in the large-scale flow and surface forcing; provide a data base for initialization and verification of coupled regional, mesoscale/hydrologic, mesoscale/chemistry, and prototype mesoscale/cloud-resolving models for prediction of severe weather, ceilings, and visibility; provide a data base for initialization and validation of cloud-resolving models, and for assisting in the fabrication, calibration, and testing of cloud and MCS parameterization schemes; and provide a data base for validation of four dimensional data assimilation schemes and algorithms for retrieving cloud and state parameters from remote sensing instrumentation.

Cotton, William↗

Influence of Vegetation on Simulation of the Water Balance and Hydrological Response to El Niño–Southern Oscillation in Western Tropical South America

Land surface models (LSMs) rely on vegetation parameters for use in hydrological and energy balance analysis, monitoring, and forecasting. This study examines the influence that vegetation representation in the Noah-Multiparameterization (Noah-MP) LSM has on hydrological simulations across the diverse climate zones of western tropical South America (WTSA), with specific consideration of hydrological variability associated with El Niño–Southern Oscillation (ENSO). The influence of model representation of vegetation on simulated hydrology is evaluated through three simulation experiments that use 1) satellite-derived constant MODIS; 2) satellite-derived time-varying MODIS; and 3) the Noah-MP dynamic leaf model. We find substantial differences in vegetation fields between these simulations, with the Noah-MP dynamic leaf model diverging significantly from satellite-derived vegetation fields in many ecoregions. Impacts on simulated hydrology were, however, found to be modest across climate zones, except for select extreme events. Also, although impacts on hydrology under ENSO-induced variability were small, we find that the Noah-MP dynamic leaf model simulates a positive relationship between rainfall and vegetation in humid ecoregions of WTSA, where satellite observations may indicate the opposite. The relatively small sensitivity of simulated hydrology to vegetation scheme suggests that the performance of hydrological monitoring and forecasting in WTSA that uses Noah-MP is largely unaffected by the choice of vegetation scheme, such that using a simple climatological default is generally no worse than adopting more complicated options. The presence of some differences between the time-varying and constant MODIS simulations for hydrologic extremes, however, indicates that time-varying MODIS configuration might be more suitable for hydrological hazards applications.

South America↗

The SMAP Level-4 ECO Project: Improving Terrestrial Flux Estimates Through Coupled Hydrology-Vegetation Data Assimilation

Simulations of hydrologic and vegetation states as well as water, energy and carbon fluxes from the land surface to the atmosphere are crucial for a wide range of applications, including agricultural advisories, forecasts of (short-term) atmospheric behavior and seasonal weather predictions including forecasts of extreme events, such as heatwaves or droughts. The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Eco-Hydrology (L4-ECO) project aims to improve modeled estimates of the terrestrial water, energy and carbon fluxes and states by developing a fully-coupled hydrology-vegetation data assimilation system. This system is developed around the NASA Goddard Earth Observing System (GEOS) Catchment-CN land surface model, which combines land hydrology and energy balance components of the GEOS Catchment model with dynamic vegetation components of the Community Land Model version 4. Catchment-CN fully couples the terrestrial water, energy and carbon cycles, allowing feedbacks from the land hydrology to the biosphere and vice versa.Here, we implement a calibration of the Catchment-CN vegetation parameterization against observations of the fraction of absorbed photosynthetically active radiation (FPAR) from the Moderate Resolution Imaging Spectroradiometer (MODIS) to improve the model's standalone skill. Later, the DA algorithm used to produce the SMAP L4 soil moisture product will be adapted to Catchment-CN to assimilate SMAP brightness temperatures and inform the model's land hydrology component. Finally, the DA system will be further extended to assimilate MODIS FPAR observations in order to constrain the model's dynamic vegetation component.In this presentation, we demonstrate that the Catchment-CN parameter calibration leads to more realistic vegetation simulations and reduces the root mean squared error between modeled and observed vegetation states across the model's various plant functional types. We also show that the assimilation of SMAP observations is able to improve the average correlation, bias and unbiased RMSE between the modeled surface and root zone soil moisture estimates, and ground observations from the SMAP core validation sites.

Kolassa, Jana↗

Evaluation of Methods for Causal Discovery in Hydrometeorological Systems

Understanding causal relations is of utmost importance in hydrology and climate research for systems identification, prediction, and understanding systems behavior in a changing climate. Traditionally, researchers in hydrometeorology attempted to study causal questions by conducting controlled experiments using numerical models. This approach, however, in most cases of interest provides uncertain results because the models are approximate representation of the natural system. An alternative approach that has recently drawn significant attention in several fields is to infer causal relations from purely observational data. It possesses several traits to its utility particularly in hydrometeorology due to the rapid accumulation of in situ and remotely sensed data records. The first objective of this study is to present a brief description of four causal discovery methods (Granger causality, Transfer Entropy, graph-based algorithms, and Convergent Cross Mapping) with special emphasis on the assumptions on which they are built. Second, using synthetic data generated from a hydrological model, we assess their performance in retrieving causal information taking into account sensitivity to sample size and presence of noise. Last, we use causal analysis to examine and formulate hypotheses on causal drivers of evapotranspiration in a shrubland region during summer and winter seasons. An interpretation of the hypotheses based on canopy seasonal dynamics and evapotranspiration processes is presented. It is hoped that the results presented here can be useful in guiding researchers studying hydrometeorological systems as to which causal method is most appropriate to the characteristics of the system under study.

54 ENVIRONMENTAL SCIENCES↗

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗