1 / 71

Isotope Hydrology Shortcourse

Isotope Hydrology Shortcourse. Residence Time Approaches using Isotope Tracers. Prof. Jeff McDonnell Dept. of Forest Engineering Oregon State University. Outline. Day 1 Morning: Introduction, Isotope Geochemistry Basics Afternoon: Isotope Geochemistry Basics ‘cont, Examples Day 2

calais
Download Presentation

Isotope Hydrology Shortcourse

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Isotope Hydrology Shortcourse Residence Time Approaches using Isotope Tracers Prof. Jeff McDonnell Dept. of Forest Engineering Oregon State University

  2. Outline • Day 1 • Morning: Introduction, Isotope Geochemistry Basics • Afternoon: Isotope Geochemistry Basics ‘cont, Examples • Day 2 • Morning: Groundwater Surface Water Interaction, Hydrograph separation basics, time source separations, geographic source separations, practical issues • Afternoon: Processes explaining isotope evidence, groundwater ridging, transmissivity feedback, subsurface stormflow, saturation overland flow • Day 3 • Morning: Mean residence time computation • Afternoon: Stable isotopes in watershed models, mean residence time and model strcutures, two-box models with isotope time series, 3-box models and use of isotope tracers as soft data • Day 4 • Field Trip to Hydrohill or nearby research site

  3. How these time and space scales relate to what we have discussed so far Bloschel et al., 1995

  4. This section will examine how we make use of isotopic variability

  5. Outline • What is residence time? • How is it determined? modeling background • Subsurface transport basics • Stable isotope dating (18O and 2H) • Models: transfer functions • Tritium (3H) • CFCs, 3H/3He, and 85Kr

  6. Residence time distribution Residence Time • Mean Water Residence Time (aka: turnover time, age of water leaving a system, exit age, mean transit time, travel time, hydraulic age, flushing time, or kinematic age) • tw=Vm/Q • For 1D flow pattern: tw=x/vpw where vpw =q/f • Mean Tracer Residence Time

  7. Why is Residence Time of Interest? • It tells us something fundamental about the hydrology of a watershed • Because chemical weathering, denitrification, and many biogeochemical processes are kinetically controlled, residence time can be a basis for comparisons of water chemistry Vitvar & Burns, 2001

  8. Tracers and Age Ranges • Environmental tracers: • added (injected) by natural processes, typically conservative(no losses, e.g., decay, sorption), or ideal(behaves exactly like traced material)

  9. Modeling Approach • Lumped-parameter models (black-box models): • System is treated as a whole & flow pattern is assumed constant over modeling period • Used to interpret tracer observations in system outflow (e.g. GW well, stream, lysimeter) • Inverse procedure; Mathematical tool: • The convolution integral

  10. Convolution • A convolution is an integral which expresses the amount of overlap of one function h as it is shifted over another function x. It therefore "blends" one function with another • It’s frequency filter, i.e., it attenuates specific frequencies of the input to produce the result • Calculation methods: • Fourier transformations, power spectra • Numerical Integration

  11. Y(w)=F(w)G(w) and • |Y(w)|2=|F(w)| 2 |G(w)| 2 The Convolution Theorem Proof: Trebino, 2002 We will not go through this!!

  12. x(t) g(t) = e -at t g(-t) e -(-at) t e -a(t-t) g(t-t) t x(t)g(t-t) t y(t) Shaded area t t t Multiplication Displacement Folding Integration Convolution: Illustration of how it works Step 1 2 3 4

  13. Example: Delta Function Convolution with a delta function simply centers the function on the delta-function. This convolution does not smear out f(t). Thus, it can physically represent piston-flow processes. Modified from Trebino, 2002

  14. Matrix Set-up for Convolution = [length(x)+length(h)]-1 = length(x) =S y(t) = x(t)*h = 0

  15. Similar to the Unit Hydrograph Precipitation Excess Precipitation Infiltration Capacity Excess Precipitation Time Tarboton

  16. Instantaneous Response Function Unit Response Function U(t) Excess Precipitation P(t) Event Response Q(t) Tarboton

  17. Subsurface Transport Basics

  18. Subsurface Transport Processes • Advection • Dispersion • Sorption • Transformations Modified from Neupauer & Wilson, 2001

  19. Advection Solute movement with bulk water flow t=t1 t2>t1 t3>t2 FLOW Modified from Neupauer & Wilson, 2001

  20. Subsurface Transport Processes • Advection • Dispersion • Sorption • Transformations Modified from Neupauer & Wilson, 2001

  21. Dispersion Solute spreading due to flowpath heterogeneity FLOW Modified from Neupauer & Wilson, 2001

  22. Subsurface Transport Processes • Advection • Dispersion • Sorption • Transformations Modified from Neupauer & Wilson, 2001

  23. Sorption Solute interactions with rock matrix FLOW t2>t1 t=t1 Modified from Neupauer & Wilson, 2001

  24. Subsurface Transport Processes • Advection • Dispersion • Sorption • Transformations Modified from Neupauer & Wilson, 2001

  25. Transformations Solute decay due to chemical and biological reactions MICROBE CO2 t2>t1 t=t1 Modified from Neupauer & Wilson, 2001

  26. Stable Isotope Methods

  27. Stable Isotope Methods • Seasonal variation of 18O and 2H in precipitation at temperate latitudes • Variation becomes progressively more muted as residence time increases • These variations generally fit a model that incorporates assumptions about subsurface water flow Vitvar & Burns, 2001

  28. Seasonal Variation in 18O of Precipitation Vitvar, 2000

  29. Deines et al. 1990 Seasonality in Stream Water

  30. Cin(t)=A sin(wt) Cout(t)=B sin(wt+j) Example: Sine-wave T=w-1[(B/A)2 –1)1/2

  31. Convolution Movie

  32. Transfer Functions Used for Residence Time Distributions

  33. Common Residence Time Models

  34. Piston Flow (PFM) • Assumes all flow paths have transit time • All water moves with advection • Represented by a Dirac delta function:

  35. Exponential (EM) • Assumes contribution from all flow paths lengths and heavy weighting of young portion. • Similar to the concept of a “well-mixed” system in a linear reservoir model

  36. Dispersion (DM) • Assumes that flow paths are effected by hydrodynamic dispersion or geomorphological dispersion • Arises from a solution of the 1-D advection-dispersion equation:

  37. Piston flow = Exponential-piston Flow (EPM) • Combination of exponential and piston flow to allow for a delay of shortest flow paths for tT (1-h-1), and g(t)=0 for t< T (1-h-1)

  38. Heavy-tailed Models • Gamma • Exponentials in series

  39. DM DM Exit-age distribution (system response function) Confined aquifer PFM: g(t’) = (t'-T) Unconfined aquifer EM: g(t’) = 1/T exp(-t‘/T) EM EM EPM PFM PFM EM Maloszewski and Zuber Kendall, 2001

  40. DM Exit-age distribution (system response function) cont… • Partly Confined Aquifer: EPM: g(t’) = /T exp(-t'/T + -1) for t‘≥T (1 - 1/) g(t’) = 0 for t'< T (1-1/ ) Kendall, 2001 Maloszewski and Zuber

  41. Dispersion Model Examples

  42. Residence Time Distributions can be Similar

  43. Uncertainty

  44. Identifiable Parameters?

  45. Review: Calculation of Residence Time • Simulation of the isotope input – output relation: • Calibrate the function g(t) by assuming various distributions of the residence time: • Exponential Model • Piston Flow Model • Dispersion Model

  46. Input Functions • Must represent tracer flux in recharge • Weighting functions are used to “amount-weight” the tracer values according recharge: mass balance!! • Methods: • Winter/summer weighting: • Lysimeter outflow • General equation: where w(t) = recharge weighting function

  47. Model 1 Cin Cout Model 3 g Cin 1- g 1- b Cout Upper reservoir g Model 2 b Direct runoff Cin 1- g Lower reservoir Cout Models of Hydrologic Systems Maloszewski et al., 1983

  48. Soil Water Residence Time Stewart & McDonnell, 2001

  49. Example from Rietholzbach Vitvar, 1998

  50. Model 3… Stable deep signal Uhlenbrook et al., 2002

More Related