Development of multi-functional measurement devices for vadose zone characterization

Size: px

Start display at page:

Download "Development of multi-functional measurement devices for vadose zone characterization"

Stanley Garrett
5 years ago
Views:

1 Development of multi-functional measurement devices for vadose zone characterization Jan Hopmans University of California, Davis, CA, USA Yasushi Mori Shimane University, Japan Annette Pia Mortensen Geological Institute, Copenhagen University, Denmark Gerard Kluitenberg Kansas State University Carlos Vaz EMBRAPA, Sao Carlos, Brazil Keith Bristow CSIRO, Townsville, Australia

2 SOIL PROPERTIES ARE NOTORIOUSLY HETEROGENEOUS, IN BOTH SPACE AND TIME QUESTIONS: Measurement Scale???? Measurement Types??? Measurement Instruments?? Measurement Locations??? Measurement Times??????

3 MULTI-FUNCTIONAL instruments, ensuring identical measurement volumes Rubber Septum Acrylic Tube Vent Tube Plug PVC Pipe Copper Tube Rubber Stopper Sample Bottle Sampling Tube PVC Dowel Porous Ceramic Cup

4 Combined tensiometer-solution sampling probe Rubber Acrylic Tube Vent Tube Plug PVC Pipe Copper Tube Rubber Tensiometer Sample Bottle Sampling Tube PVC Dowel Porous Ceramic Cup Soil Solution Sampling

Coiled Cone Penetrometer-TDR Probe o o o o o Two parallel copper wires are wrapped around inner core (as double helix); Wires are connected to conductor and ground of

5 Coiled Cone Penetrometer-TDR Probe o o o o o Two parallel copper wires are wrapped around inner core (as double helix); Wires are connected to conductor and ground of coaxial cable; Signal is analyzed by cable tester; Long wire length (about 30 cm) ensures accurate travel time measurement,and Narrow wire spacing ensure high depth resolution

6 Hammer penetrometer parallel wires (steel) tip (steel) 2.4 cm

7 Calibration ε coil Columbia Yolo Sand Fit Columbia Fit Yolo Fit Sand Polynomial fit Yolo θ (cm 3 cm -3 ) Penetration Resistance (MPa) to 1.3 g cm g cm 1.3 to 1.4 g cm r 2 = 0.72 RMSE = to 1.5 g cm to 1.6 g cm -3 equation [3] Water Content (cm 3 cm -3 )

8 Field development, Brazil

Combined Tensiometer-TDR Reflection Coefficient 1.0 0.9 0.8 0.7 0.6 0.5 0.

9 Combined Tensiometer-TDR Reflection Coefficient dried glass beads sat. glass beads water Time (ns)

10 Laboratory Calibration 24 Oso Flaco 22 Ottawa SRI Columbia 20 Lincoln 18 ε coil Soil Water Content (cm 3 cm -3 )

11 A heat pulse probe? Instrumentation in the vadose zone TDR Time domain reflectometry Tensiometers Tracer experiments Heat pulse probe multifunctional Apply heat as tracer Heat, water and solute transport

12 Multifunctional heat pulse probe Temperature, T Thermal properties Heat capacity, C Heat conductivity, λ 0 Thermal diffusivity, κ Heat dispersion, D Hydraulic properties same time + same place + same scale Water flux, q w Water content, θ Electrical conductivity, EC a

13 Heat pulse probe design Heat pulse heater thermistor 6 needles 1 mm diameter 6 mm spacing 28 mm long 25 mm wide Electrical conductivity Wenner array

14 I. Multi-step outflow experiment Multiplexer AM416 MF-HPP Multi-plexer AM416 Heating Sensitometer Tensiometer and Outflow Four-electrode Datalogger CR10 Thermistor Datalogger CR10

15 I. Analytical solutions of heat transport De Vries (1952)- Thermal Conduction ( ) ( ) ; ', t t t r Ei t t r Ei C q t r T > κ κ κ π = Ren et al. (2000) Thermal Convection ( ) t t ; ds s s V r exp s C q T t t t h d d > = κ κ π ( ) t t ; ds s s V r exp s C q T t t t h u u > + = κ κ π

16 II. Experimental flow column pump rainmaking device multiplexer AM416 tensiometer and pressure transducer multiplexer AM416 HPP datalogger CR10

17 II. Numerical solution heat transport Heat conduction and convection Homogenoues and isotropic media Thermal equilibrium between phases x C V 2 2 T T T w T = κ + V 2 2 θ w t x z Cb z h Cq θcv = = C C w w w w bulk b z T U H Water flow r r Cbulk = Cs(1 φ) + C w θ D

18 Fitting of temperature response curve 0.7 Fitting 0.6 I. Analytical solution Temperature difference [C] Downstream Upstream Transverse II. Numerical solution (HYDRUS 2D) Parameters Thermal conductivity λ Water content θ 0 Water flux q w Time [s]

19 I and II: Calibration of MF-HPP Needle distance, r Calibration in agar solution heat capacity for water no convective heat transport Calibration in porous media heat capacity for material saturated conditions Wenner array Calibration in porous media varying EC and water content Temperature difference [K] ECb (ds m -1 ) M 0.06 M 0.10 M Fitted Time [s] θ by MFHPP (m 3 m -3 )

20 I.Unsaturated hydraulic functions Cumulative outflow (cm 3 ) and -matric head (cm) Cumulative outflow Matric head Time (hours) -Matric head (cm) M 0.06 M 0.10 M Volumetric water content (m 3 m -3 ) Log 10 Hydraulic conductivity (m s -1 ) Inverse Modeling of Multi-step Outflow

21 I and II. Estimation of thermal conductivity, λ ( heat conduction only) 2.0 Thermal conductivity λ [W/mK] Mortensen et al. [2003] Mori et al [2003] Hopmans and Dane [1986] Thermal conductivity λ 0 = b + b θ + b θ 0 Thermal diffusivity κ = λ 0 1 /C bulk Water content [cm3/cm3]

22 II. Solute transport Electrical conductivity EC = θτ(θ) EC + bulk water EC solid v= cm/s 2.5 θ= θ=0.22 θ=0.32 EC [ms/cm] θ= Ψ=20 cm Time [min]

23 I and II: Estimation of water content, θ (Heat conduction only) Real water content Estimated water content Real water content (m3m-3) Static Transient Estimated water content (m3m-3)

24 Water flux effect on temperature signature Temperature ( o C) No Flow 1 m/d, downstream needle 1 m/d, upstream needle 10 m/d, downstream needle 10 m/d, upstream needle Time (s) (B)

25 I and II: Estimation of water flux, q w (heat conduction and convection) Accurate range: to 0.01 cm/s or 10 to 1000 cm/day 8.E-5 7.E-5 6.E-5 Real flux [m/s] 5.E-5 4.E-5 3.E-5 2.E-5 1.E-5 Unsaturated Saturated 0.E+0 0.E+0 1.E-5 2.E-5 3.E-5 4.E-5 5.E-5 6.E-5 7.E-5 8.E-5 Estimated flux [m/s]

26 Water content error for high water fluxes However, these high water fluxes only occur under saturated conditions

27 Thermal dispersion effect on temperature signature (q w = 1 m d -1 ) Temperature ( o C) Dispersivity = 0 Downstream Upstream Transverse Temperature ( o C) Downstream Upstream Transverse Dispersivity = Time (s) Time (s) Temperature ( o C) Downstream Upstream Transverse Dispersivity = 0.01 Temperature ( o C) Downstream Upstream Transverse Dispersivity = Time (s) Time (s)

28 Conclusions Advantages Simultaneous measurement of water flow, solute and heat transport properties within the same sample volume Limitations In situ calibration of probe Sensibility to needle spacing Accurate water flux range is limited to cm/day Future work Compare heat and solute dispersivities Improve for water flux density < 10 cm/day Field applications

29 Multi-Functional Sensors: Are being developed, that Measure approximately same volumes, BUT Issues of probe-soil contact Measurement interferences What is real measurement volume, contributing to the measurement???

30 REV (Bear) or Relativist Concept (Baveye and Sposito, 1984)???? SENSOR Spatial weighting function, ω: ω = V ( xyzdv,, ) 1 θ ω θ mac ( xyz,, ) = ( xyz,, ) ( xyzdv,, ) mic V

Indirect estimation of soil thermal properties and water flux using heat pulse probe measurements: Geometry and dispersion effects

WATER RESOURCES RESEARCH, VOL. 38, NO. 1, 1006, 10.1029/2000WR000071, 2002 Indirect estimation of soil thermal properties and water flux using heat pulse probe measurements: Geometry and dispersion effects