Water Budget II: Evapotranspiration P = Q + ET + G + ΔS

Size: px
Start display at page:

Download "Water Budget II: Evapotranspiration P = Q + ET + G + ΔS"

Transcription

1 Water Budget II: Evapotranspiration P = Q + ET + G + ΔS

2

3 Evaporation Transfer of H 2 O from liquid to vapor phase Diffusive process driven by Saturation (vapor density) gradient ~ (ρ s ρ a ) Aerial resistance ~ f(wind speed, temperature) Energy to provide latent heat of vaporization (radiation) Transpiration is plant mediated evaporation Same result (water movement to atmosphere) Summative process = evapotranspiration (ET) Dominates the fate of rainfall ~ 95% in arid areas ~ 70% for all of North America

4 ET is the sum of Evapotranspiration Evaporation: physical process from free water Soil Plant intercepted water Lakes, wetlands, streams, oceans Transpiration: biophysical process modulated by plants (and animals) Controlled flow through leaf stomata Species, temperature and moisture dependent

5 Global Rates

6 National Rates

7 Spatial Variation in P - ET Controls streamflow (mm/yr) Beck et al. (2015)

8 Globally (P ET in mm/yr) Beck et al. (2015)

9 Four Requirements for ET TP Energy Water Vapor Pressure Gradient NP Wind

10 NASA 3850 zettajoules per year

11 Radiation Budget Energy Inputs R total = Total Solar Radiation Inputs on a horizontal plane at the Earth s Surface R net = R total reflected radiation = R total * (1 albedo) Albedo (α) values Snow 0.9 Hardwoods 0.2 Water 0.05 Flatwoods pine plantation 0.15 Flatwoods clear cut Burn Asphalt 0.05

12 Energy and Temperature The simplest conceptualization of the ET process focuses solely on temperature. Blaney-Criddle Method: ET = p * (0.46*T mean + 8) Where p is the mean daytime hours T mean is the mean daily temp (Max+Min/2) ET (mm/day) is treated as a monthly variable

13 Vapor Deficit Drives the Process Distance between actual conditions and saturation line Greater distances = larger evaporative potential Slope of this line (δ) is an important term for ET models Usually measured in mbar/ C Graph shows mass water per mass air as a function of T

14 Wind Boundary layer saturates under quiescent conditions Inhibits further ET UNLESS air is replaced Turbulence at boundary layer is therefore necessary to ensure a steady supply of undersaturated air

15 Water Availability: PET vs. AET PET (potential ET) is the expected ET if water is not limiting Given conditions of: wind, Temperature, Humidity AET (actual ET) is the amount that is actually abstracted (realizing that water may be limiting) AET = α * PET Where α is a function of soil moisture, species, climate ET:PET is low in arid areas due to water limitation ET ~ PET in humid areas due to energy limitation

16 Budyko Revisited Budyko Curve shows annual catchment hydrology Can neglect G and ΔS, focus on ET (P ~ ET + Q) Aridity index on the x-axis (potential ET:rain) Evaporative index on the y axis (actual ET:rain)

17 Methods of Estimating ET Since ET is the largest flux OUT of the watershed, we need good estimates Techniques have focused mostly on predicting capacity (i.e., PET, where water is not limiting) energy balance methods mass transfer or aerodynamic methods combination of energy and mass transfer (Penman equation) pan evaporation data

18 Evaporation from a Pan Mass balance equation S = I 0 National Weather Service Class A type Installed on a wooden platform in a grassy location Filled with water to within 2.5 inches of the top Evaporation rate is measured by manual readings or with an analog output evaporation gauge H2 E p H1 = P Pans measure more evaporation than natural water bodies because: E = P ( H2 H1) 1) less heat storage capacity (smaller volume) 2) heat transfer 3) wind effects

19 Energy Balance Method Assumes energy supply the limiting factor (Budyko). sensible heat transfer to air net radiation energy used in evaporation H s R n Q e heat stored in system G heat conducted to ground (typically neglected) Consider energy balance on a small lake with no water inputs (or evaporation pan)

20 Energy Budget Energy in = Energy out (conservation law) Energy In = R total Energy Out Albedo Latent Heat Sensible Heat Soil Heat Flux If R total = 800 cal/cm 2 /day and α = 0.2 R net = 800 * (1 0.2) = 640 cal/cm 2 /day

21 Energy Budget Estimates of ET R net = λe + H + G We want to know E E = (R net (H+G))/λ What are evaporative losses if: R total = 800 cal/cm 2 /day Albedo = 0.2 λ = 586 cal/g H = 100 cal/cm 2 /d (convected heat) G = 50 cal/cm 2 /d (soil heating)

22 Static Computation R net = λe + H + G = 800 * (1 0.2) = = (586*E) E = 0.84 cm/d Annual ET = 0.84 * 365 * 1 m/100 cm = 3.07 m R tot = 800 cal/cm/d Albedo = 0.2 λ =586 cal/g H = 100 cal/cm/d lost G = 50 cal/cm/d lost to ground

23 Energy Budget Bowen Ratio β = H/λE G

24 Mass Transfer (Aerodynamic) Method Assumes rate of turbulent mass transfer of water vapor E = B( u)( e from evaporating surface to s e( z)) atmosphere is limiting factor Mass transfer is controlled 0.102u B( u) = by (1) vapor gradient (e s e) 2 and (2) wind velocity (u) ln z2 z o which determines rate at u which vapor is carried away. B( u) = (1 + ) 100

25 Combination Method (Penman) Evaporation can be estimated by aerodynamic method (E a ) when energy supply not limiting and energy method (E r ) when vapor transport not limiting Typically both factors limiting so use combination of above methods E = γ E r + E a + γ + γ Weighting factors sum to 1. = vapor pressure deficit γ = psychrometric constant = 4098 e s ( ) T γ 66.8Pa / 0 C

26 Combination Method (Penman) Penman is most accurate and commonly used method if meteorological information is available. Need: net radiation, air temperature, humidity, wind speed If not available use Priestley-Taylor approximation: E = α + γ Based on observations that second term (advection) in Penman equation typically small in low water stress areas. The α term is crop coefficient that assumes no advection limitation. Usually >1 (1.2 to 1.7), suggesting that actual ET is greater than what is predicted from radiation alone. E r

27 Diurnal Water Level Variation (White, 1932) Diel variation in water level yields ET (during the day) and net groundwater flux (at night) Curiously, not widely used 0.51 h (cm/hr) Water Level ET = S y (S + 24 x h) Exfiltration = S y (24 x h) S :00 12:00 0:00 12:00 0:00

28 Actual Diurnal Data Water Level (m) Water Level (m) :00 12:00 0:00 12:00 0:00 12: h ET/S s y h s ET/S y 0:00 12:00 0:00 12:00 0:00 12:00 Nighttime slope is groundwater flow (inflow is UP, outflow is DOWN) Assuming constant GW flow, daytime slope is ET + GW. Specific yield (S y )

29 What is Specific Yield? How much water (in units of cm) drains out of a soil; also called dynamic drainable porosity

30 Eddy Correlation Measures Moving air is combination of multiple scales of eddies (turbulent gyres) Measure gas concentrations (water vapor, CO 2, CH 4 ) in up-welling vs. down-welling air is ecosystem fluxes (correlation between concentration and flow direction)

31 Time Scales of Variability Controls on ET create variability at scales from seconds to centuries Eddies change ET at the time scale of seconds Diel cycles affect water fluxes over 24 hours Weather patterns affect fluxes at days to weeks Water availability Vapor deficit Wind and energy Climate variability at decadal and beyond

32 High Resolution ET Observations

33

34

35 Total System ET Ordered Process Intercepted Water Transpiration Surface Water Soil Water Why? Implications for: Cloud forests Understory vegetation in wetlands Deep rooted arid ecosystems

36 Evapotranspiration has Multiple Components

37 Interception Surface tension holds water falling on forest vegetation. Leaf Storage Fir 0.25 Pines 0.10 Interception Loss (% of rainfall) Hardwoods 0.05 Hardwoods 10-20% (less LAI) Litter 0.20 Conifers 20-40% SP Plantations Mixed slash and Cypress Florida Flatwoods 20%

38 Transpiration Plant mediated diffusion of soil water to atmosphere Soil-Plant-Atmosphere Continuum (SPAC) Transpiration and productivity are tightly coupled Transpiration is the primary leaf cooling mechanism under high radiation Provides a pathway for nutrient uptake and matrix for chemical reactions Worldwide, water limitations are more important than any other limitation to plant productivity CO 2 H 2 O 1 : 300

39 Transpiration Dominates the Evaporation Process Trees have: Large surface area More turbulent air flow Conduits to deeper moisture sources T/ET Hardwood ~80% White Pine~60% Flatwoods ~75%

40 Cover Evaporation Interception Transpiration Forest 10% 30% 60% Meadow 25% 25% 50% Ag 45% 15% 40% Bare 100%

41 The SPAC (soil-plant-atmosphere continuum) Ψ w (atmosphere) -95 MPa Ψ w (small branch) -0.8 MPa Ψ w (stem) -0.6 MPa Ψ w (soil) -0.1 MPa Ψ w (root) -0.5 MPa

42 The driving force of transpiration is the difference in water vapor concentration, or vapor pressure difference, between the internal spaces in the leaf and the atmosphere around the leaf

43 Transpiration The physics of evaporation from stomata are the same as for open water. The only difference is the conductance term. Conductance is a two step process stomata to leaf surface leaf surface to atmosphere

44 Transpiration

45 How Does Water Get to the Leaf? Water is PULLED, not pumped. Water within the whole plant forms a continuous network of liquid columns from the film of water around soil particles to absorbing surfaces of roots to the evaporating surfaces of leaves. It is hydraulically connected.

46 Even a perfect vacuum can only pump water to a maximum of a little over 30 feet. At this point the weight of the water inside a tube exerts a pressure equal to the weight of the atmosphere pushing down So why doesn t the continuous column of water in trees taller than 34 feet collapse under its own weight? And how does water move UP a tall tree against the forces of gravity? > 100 meters

47 Water is held up by the surface tension of tiny menisci ( menisci is the plural of meniscus) that form in the microfibrils of cell walls, and the adhesion of the water molecules to the cellulose in the microfibrils cell wall microfibrils of carrot

48 Cohesion-Tension Theory: (Böhm, 1893; Dixon and Joly, 1894) The cohesive forces between water molecules keep the water column intact unless a threshold of tension is exceeded (embolism). When a water molecule evaporates from the leaf, it creates tension that pulls on the entire column of water, down to the soil.

49 ? ET = Rain * 0.80 ET = Rain * ,000 mm * 0.80 = 800 mm 1,000 mm * 0.95 = 950 mm G = 200 mm Assume Q & ΔS = 0 G = P - ET G = 50 mm 4x more groundwater recharge from open stands than from highly stocked plantations. NRCS is currently paying for growing more open stands, mainly for wildlife.

50

51 Controls on Stand Water Use More leaves per area = more water use Foresters don t measure LAI Proxies for LAI Basal area Age Height Site Index

52 A Fair Comparison Pine stands are clear-cut every years (low ET) Compare water yield (Rain ET) over entire rotation

53 Ecosystem Service Water Yield Forest management may yield new water Win-win for other forest services Who pays and how much?

54 Trading Environmental Priorities? Water for Carbon Water for Energy Jackson et al (Science)

55 Surface Water Evaporation Air Temp Air relative humidity Water temp Wind Radiation Water Quality Actual surface water evaporation ~ pan evaporation * 0.7

56 Soil Water Evaporation Stage 1. For soils saturated to the surface, the evaporation rate is similar to surface water evaporation. Stage 2. As the surface dries out, evaporation slows to a rate dependent on the capillary conductivity of the soil. Stage 3. Once pore spaces dry, water loss occurs in the form of vapor diffusion. Vapor diffusion requires more energy input than capillary conduction and is much, much, slower. Note that for soils under a forest canopy, R net, vapor pressure deficit, and turbulent transport (wind) are lower than for exposed soils.

57 Soil water loss with different cover

58 Rooting Depth Effects Surface 2 months later

59 Streamflow Next Time