SEARCH · Search NASA

Results for “subsurface flow”

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.

SourceResearch connectionFrom yearThrough year

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

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2013

Preliminary Analysis of Grande Ronde Basalt Formation Flow Top Transmissivity as it Relates to Assessment and Site Selection Applications for Fluid/Energy Storage and Sequestration Projects

The following summary comments pertain to a preliminary assessment of transmissivity conditions within the Grande Ronde Basalt Formation flow tops (i.e., interflow zones), as they may relate to a variety of fluid/energy storage or sequestration projects within deep Columbia River Basalt formations. Transmissivity is the principal hydraulic parameter controlling well/reservoir injectivity and/or extraction of fluids within subsurface formations. In this regard, the preliminary findings and statistical relationships developed herein regarding Grande Ronde Basalt flow-top transmissivity are considered relevant for feasibility assessments and site selection applications within deep basalts.

Carbon sequestration↗

DOE OSTI2022

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗

DOE OSTI2019

Biofilm Feedbacks Alter Hydrological Characteristics of Fractured Rock Impacting Sulfidogenesis and Treatment

Flow-through fractures dominate the movement of fluids in a variety of natural as well as engineered subsurface systems. Microbial activities in fractured rock impact subsurface energy recovery, storage, and waste disposal. It has been recognized that understanding how the contrasting permeability between fracture and matrix interacts with microbial metabolism under thermal and hydrological gradients is key to effective utilization of the subsurface, yet such studies are sparse. Microorganisms mediate the production of hydrogen sulfide (also known as souring) in oil-bearing geological formations. We conducted a comprehensive experimental study of a novel 2D fractured rock system to understand these complex interactions and demonstrated how biofilm development can impact fracture flow, which subsequently feedbacks to moderate sulfidogenesis. Elevated temperature relevant to reservoir conditions interacted with the injection of cold fluid and formed a thermal gradient away from fractures, creating thermal niches for microbial activities in the fractured rock. Finally, results showed that while fracture flows were dominant in the beginning, with time, growth of the biofilm in the fractures reduced permeability, effectively moderating the initial fracture-matrix contrast, and limited microbial accessibility to nutrients and subsequent reactions rates.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2012

Conceptualization of the Predevelopment Groundwater Flow System and Transient Water-Level Responses in Yucca Flat, Nevada National Security Site, Nevada

Contaminants introduced into the subsurface of Yucca Flat, Nevada National Security Site, by underground nuclear testing are of concern to the U.S. Department of Energy and regulators responsible for protecting human health and safety. The potential for contaminant movement away from the underground test areas and into the accessible environment is greatest by groundwater transport. The primary hydrologic control on this transport is evaluated and examined through a set of contour maps developed to represent the hydraulic-head distribution within the two major aquifer systems underlying the area. Aquifers and confining units within these systems were identified and their extents delineated by merging and analyzing hydrostratigraphic framework models developed by other investigators from existing geologic information. Maps of the hydraulic-head distributions in the major aquifer systems were developed from a detailed evaluation and assessment of available water-level measurements. The maps, in conjunction with regional and detailed hydrogeologic cross sections, were used to conceptualize flow within and between aquifer systems.

58 GEOSCIENCES↗

DOE OSTI2021

Permafrost Promotes Shallow Groundwater Flow and Warmer Headwater Streams

Abstract The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from 11 headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7% of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5% moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Directly imaging steeply-dipping fault zones in geothermal fields with multicomponent seismic data

For characterizing geothermal systems, it is important to have clear images of steeply-dipping fault zones because they may confine the boundaries of geothermal reservoirs and influence hydrothermal flow. Elastic reverse-time migration (ERTM) is the most promising tool for subsurface imaging with multicomponent seismic data. However, conventional ERTM usually generates significant artifacts caused by the cross correlation of undesired wavefields and the polarity reversal of shear waves. In addition, it is difficult for conventional ERTM to directly image steeply-dipping fault zones. We develop a new ERTM imaging method in this paper to reduce these artifacts and directly image steeply-dipping fault zones. In our new ERTM method, forward-propagated source wavefields and backward-propagated receiver wavefields are decomposed into compressional (P) and shear (S) components. Furthermore, each component of these wavefields is separated into left- and right-going, or downgoing and upgoing waves. The cross correlation imaging condition is applied to the separated wavefields along opposite propagation directions. For converted waves (P-to-S or S-to-P), the polarity correction is applied to the separated wavefields based on the analysis of Poynting vectors. Numerical imaging examples of synthetic seismic data demonstrate that our new ERTM method produces high-resolution images of steeply-dipping fault zones.

58 GEOSCIENCES↗

DOE OSTI2017

Washington Geothermal Play Fairway Analysis Data From Potential Field Studies

A recent study which adapts play fairway analysis (PFA) methodology to assess geothermal potential was conducted at three locations (Mount Baker, Mount St. Helens seismic zone, and Wind River valley) along the Washington Cascade Range (Forson et al. 2017). Potential field (gravity and magnetic) methods which can detect subsurface contrasts in physical properties, provides a means for mapping and modeling subsurface geology and structure. As part of the WA-Cascade PFA project, we performed potential field studies by collecting high-resolution gravity and ground-magnetic data, and rock property measurements to (1) identify and constrain fault geometries (2) constrain subsurface lithologic distribution (3) study fault interactions (4) identify areas favorable to hydrothermal flow, and ultimately (5) guide future geothermal exploration at each location.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Performance, Cost, and Financial Parameters of Geothermal District Heating Systems for Market Penetration Modeling under Various Scenarios

Geothermal district heating (GDH) systems have limited penetration in the U.S., with an estimated installed capacity of only 100 MWth for a total of 21 sites. We see higher deployment in other regions, for example, in Europe with an installed capacity of more than 4,700 MWth for 257 GDH sites. The U.S. Department of Energy Geothermal Vision (GeoVision) Study is currently looking at the potential to increase the deployment in the U.S. and to understand the impact of this increased deployment. This paper reviews 31 performance, cost, and financial parameters as input for numerical simulations describing GDH system deployment in support of the GeoVision effort. The focus is on GDH systems using hydrothermal and Enhanced Geothermal System resources in the U.S.; ground-source heat pumps and heat-to-electricity conversion technology were excluded. Parameters investigated include 1) capital and operation and maintenance costs for both subsurface and surface equipment; 2) performance factors such as resource recovery factors, well flow rates, and system efficiencies; and 3) financial parameters such as inflation, interest, and tax rates. Current values as well as potential future improved values under various scenarios are presented. Sources of data considered include academic and popular literature, software tools such as GETEM and GEOPHIRES, industry interviews, and analysis conducted by other task forces for the GeoVision Study, e.g., on the drilling costs and reservoir performance.

15 GEOTHERMAL ENERGY↗

DOE OSTI2015

Adsorption kinetics dictate monolayer self-assembly for both lipid-in and lipid-out approaches to droplet interface bilayer formation

The droplet interface bilayer (DIB) is a method to assemble planar lipid bilayer membranes between lipid-coated aqueous droplets and has gained popularity among researchers in many fields. Well-packed lipid monolayer on aqueous droplet–oil interfaces is a prerequisite for successfully assembling DIBs. Such monolayers can be achieved by two different techniques: “lipid-in”, in which phospholipids in the form of liposomes are placed in water, and “lipid-out”, in which phospholipids are placed in oil as inverse micelles. While both approaches are capable of monolayer assembly needed for bilayer formation, droplet pairs assembled with these two techniques require significantly different incubation periods and exhibit different success rates for bilayer formation. In our study, we combine experimental interfacial tension measurements with molecular dynamics simulations of phospholipids (DPhPC and DOPC) assembled from water and oil origins to understand the differences in kinetics of monolayer formation. With the results from simulations and by using a simplified model to analyze dynamic interfacial tensions, we conclude that, at high lipid concentrations common to DIBs, monolayer formation is simple adsorption controlled for lipid-in technique, whereas it is predominantly adsorption-barrier controlled for the lipid-out technique due to the interaction of interface-bound lipids with lipid structures in the subsurface. The adsorption barrier established in lipid-out technique leads to a prolonged incubation time and lower bilayer formation success rate, proving a good correlation between interfacial tension measurements and bilayer formation. Also, we clarify that advective flow expedites monolayer formation and improves bilayer formation success rate by disrupting lipid structures, rather than enhancing diffusion, in the subsurface and at the interface for lipid-out technique. Additionally, electrical properties of DIBs formed with varying lipid placement and type are characterized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

A general peridynamics model for multiphase transport of non-Newtonian compressible fluids in porous media

Herein a general state-based peridynamics model is developed to simulate transport of fluids in an arbitrary heterogeneous porous medium. The generality encompasses modeling of multiphase, multi-component flow of non-Newtonian and compressible fluids, which is often encountered in but not limited to subsurface reservoirs. Peridynamic model is especially useful for solving non-local problems, such as crack propagation, since it does not assume spatial continuity of field variables. Thus, the formulation presented here, combined with peridynamics-based damage model, can be used to simulate hydraulic fracturing with complex fluids. To demonstrate its capability to simulate multi-phase flow in porous media, the derived model is verified against the analytical Buckley-Leverett solution for immiscible Newtonian two-phase flow. Further, the non-Newtonian two-phase fluid flow in porous media is verified by simulating the polymer flood process involving immiscible displacement of a Newtonian fluid by a non-Newtonian fluid against a generalized solution obtained by Wu et al. The non-local solutions are shown to be consistent with the corresponding local solutions in limiting cases. Moreover, mass conservation of all the phases is satisfied, irrespective of discretization and extent of non-locality.

42 ENGINEERING↗

DOE OSTI2015

EOS7CA Version 1.0: TOUGH2 Module for Gas Migration in Shallow Subsurface Porous Media Systems

EOS7CA is a TOUGH2 module for mixtures of a non-condensible gas (NCG) and air (with or without a gas tracer), an aqueous phase, and water vapor. The user can select the NCG as being CO2, N2, or CH4. EOS7CA uses a cubic equation of state with a multiphase version of Darcy's Law to model flow and transport of gas and aqueous phase mixtures over a range of pressures and temperatures appropriate to shallow subsurface porous media systems. The limitation to shallow systems arises from the use of Henry's Law for gas solubility which is appropriate for low pressures but begins to over-predict solubility starting at pressures greater than approximately 1 MPa (10 bar). The components modeled in EOS7CA are water, brine, NCG, gas tracer, air, and optional heat.

Oldenburg, CurtisM↗

DOE OSTI2015

EOS7CA Version 1.0: TOUGH2 Module for Gas Migration in Shallow Subsurface Porous Media Systems

Oldenburg, Curtis M.↗

DOE OSTI2013

UBIQUITOUS SOLAR ERUPTIONS DRIVEN BY MAGNETIZED VORTEX TUBES

The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruption events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by the Sun's turbulent convection in subsurface layers. The swirling vortex tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and stretch background magnetic field, and push the surrounding material up, generating shocks. Our simulations reveal complicated high-speed flow patterns and thermodynamic and magnetic structure in the erupting vortex tubes. The main new results are: (1) the eruptions are initiated in the subsurface layers and are driven by high-pressure gradients in the subphotosphere and photosphere and by the Lorentz force in the higher atmosphere layers; (2) the fluctuations in the vortex tubes penetrating into the chromosphere are quasi-periodic with a characteristic period of 2-5 minutes; and (3) the eruptions are highly non-uniform: the flows are predominantly downward in the vortex tube cores and upward in their surroundings; the plasma density and temperature vary significantly across the eruptions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

The Role of Chemistry in Fracture Pattern Development and Opportunities to Advance Interpretations of Geological Materials

Abstract Fracture pattern development has been a challenging area of research in the Earth sciences for more than 100 years. Much has been learned about the spatial and temporal complexity inherent to these systems, but severe challenges remain. Future advances will require new approaches. Chemical processes play a larger role in opening‐mode fracture pattern development than has hitherto been appreciated. This review examines relationships between mechanical and geochemical processes that influence the fracture patterns recorded in natural settings. For fractures formed in diagenetic settings (~50 to 200 °C), we review evidence of chemical reactions in fractures and show how a chemical perspective helps solve problems in fracture analysis. We also outline impediments to subsurface pattern measurement and interpretation, assess implications of discoveries in fracture history reconstruction for process‐based models, review models of fracture cementation and chemically assisted fracture growth, and discuss promising paths for future work. To accurately predict the mechanical and fluid flow properties of fracture systems, a processes‐based approach is needed. Progress is possible using observational, experimental, and modeling approaches that view fracture patterns and properties as the result of coupled mechanical and chemical processes. A critical area is reconstructing patterns through time. Such data sets are essential for developing and testing predictive models. Other topics that need work include models of crystal growth and dissolution rates under geological conditions, cement mechanical effects, and subcritical crack propagation. Advances in machine learning and 3‐D imaging present opportunities for a mechanistic understanding of fracture formation and development, enabling prediction of spatial and temporal complexity over geologic timescales. Geophysical research with a chemical perspective is needed to correctly identify and interpret fractures from geophysical measurements during site characterization and monitoring of subsurface engineering activities.

58 GEOSCIENCES↗

DOE OSTI2015

EOS7CA Version 1.0: TOUGH2 Module for Gas Migration in Shallow Subsurface Porous Media Systems

EOS7CA is a TOUGH2 module for mixtures of a non-condensible gas (NCG) and air with or without a gas tracer, an aqueous phase, and water vapor. The user can select the NCG as being CO 2 , N 2 , or CH 4 . EOS7CA uses a cubic equation of state with a multiphase version of Darcy’s Law to model flow and transport of gas and aqueous phase mixtures over a range of pressures and temperatures appropriate to shallow subsurface porous media systems. The limitation to shallow systems arises from the use of Henry’s Law for gas solubility which is appropriate for low pressures but begins to over-predict solubility starting at pressures greater than approximately 1 MPa (10 bar). The components modeled in EOS7CA are water, brine, NCG, gas tracer, air, and optional heat. The real gas properties module (ZEVCA) has options for Peng-Robinson, Redlich-Kwong, or Soave-Redlich-Kwong equations of state to calculate gas mixture density, enthalpy departure, and viscosity. Transport of the gaseous and dissolved components is by advection and Fickian molecular diffusion. This user guide provides instructions for use and two sample problems as verification and demonstration of EOS7CA.

03 NATURAL GAS↗

DOE OSTI2025

Identifying preferential flow from soil moisture time series: Review of methodologies

Abstract Identifying and quantifying preferential flow (PF) through soil—the rapid movement of water through spatially distinct pathways in the subsurface—is vital to understanding how the hydrologic cycle responds to climate, land cover, and anthropogenic changes. In recent decades, methods have been developed that use measured soil moisture time series to identify PF. Because they allow for continuous monitoring and are relatively easy to implement, these methods have become an important tool for recognizing when, where, and under what conditions PF occurs. The methods seek to identify a pattern or quantification that indicates the occurrence of PF. Most commonly, the chosen signature is either (1) a nonsequential response to infiltrated water, in which soil moisture responses do not occur in order of shallowest to deepest, or (2) a velocity criterion, in which newly infiltrated water is detected at depth earlier than is possible by nonpreferential flow processes. Alternative signatures have also been developed that have certain advantages but are less commonly utilized. Choosing among these possible signatures requires attention to their pertinent characteristics, including susceptibility to errors, possible bias toward false negatives or false positives, reliance on subjective judgments, and possible requirements for additional types of data. We review 77 studies that have applied such methods to highlight important information for readers who want to identify PF from soil moisture data and to inform those who aim to develop new methods or improve existing ones. Core Ideas Soil moisture data can be used to identify the occurrence of preferential flow (PF) and its initiating conditions. Various data‐analysis methods to identify PF differ in susceptibility to error, bias, and subjectivity. These methods can utilize vast amounts of data from soil moisture monitoring networks to develop understanding of when, where, and under what conditions PF occurs. Newly developed methods may lead to better accuracy and reliability, and reduce the need for subjective judgments. Plain Language Summary Preferential flow through soil occurs when a large amount of water is suddenly available, as during an intense storm. This type of flow moves rapidly through the soil in distinct narrow pathways rather than moving evenly throughout the body of soil, with major consequences for groundwater resources, ecosystems, spreading of contaminants, and other vital concerns. Methods of detecting preferential flow have been developed that utilize measurements of soil water content made by sensors installed at various depths. This measurement technology has been widely implemented, many locations now having datasets years in length, and various methods have been developed for using these to identify preferential flow. The various methods are based on different features in the soil moisture records and vary in their advantages and shortcomings. In this review, we explain and evaluate these methods, highlighting important information for their implementation to identify preferential flow from soil moisture data and for efforts to develop new methods or improve existing ones.

Nimmo, John R↗

DOE OSTI2020

Catchment storage and residence time in a periodically irrigated watershed

Understanding anthropogenic impacts on water storage and water flow pathways in catchments is an ongoing challenge in hydrology. Here, we study the dynamics of subsurface storage and residence time of water in a catchment in Berkeley, California, that is within a regional park but contains diverse land use within its perimeter, including a periodically irrigated golf course. Our study combines several isotopic tracers with water budget data to examine sources of water in a stream draining the site. Irrigation water, applied to a small area of the watershed, is a minor component of the water budget. However, geochemical tracers reveal that irrigation water is a significant fraction of stream flow downstream of the golf course during baseflow and during precipitation events. Isotopic tracers indicate that the watershed has a preference to release young water for stream flow generation, resulting in contrasting tritium ages for stream water and groundwater of 1.3 ± 0.5 year and 8.2 ± 1.7 year, respectively. We determined that the older water is a very small component (0.7%) of the stream water in the tail of an assumed exponential distribution. We used the seasonal variation of stable water isotopes in precipitation and stream water over two water years to explain the damping of the isotopic signature of stream water, which yields information about the catchment's response to the input signal. In conclusion, the methods described here may be applicable to other urban or suburban headwater catchments in areas with a component of non-natural recharge from, for example, leaky infrastructure, storm water routing or dry season irrigation.

54 ENVIRONMENTAL SCIENCES↗