Partitioning plant transpiration and soil evaporation with eddy covariance and stable isotope method in North China Plain Prof. Mei Xurong, Theme Leader Scientist Director General, IEDA, CAAS Director, Key Laboratory of Agro-environment, MOA Director, Consultant Group for Agricultural Disaster Mitigation, MOA meixr@ieda.org.cn
Introduction Water & Food in the North China Plain Total Area: 303,585 km 2 Population: 300Million Water resource: 13.5 BCM(5%) Rainfall: 550-650 mm Cropping System: Winter wheat- Summer Maize Cereal Production: 70% of wheat 30% of Maize
Annual P & ET 年蒸散量 年降水量 (mm) 水量 (mm ) Water Scarcity in the North China Plain 180 160 140 120 100 80 60 40 20 0 亏缺水量 Water Deficit 作物需水量 ET 0 降水量 Rainfall Oct Nov Dec Jan Feb Mar Apr May June July August Sept Month Winter Wheat Summer Maize 1000 800 600 400 200 Precipitation 降水 蒸散 Evapotranspiration -30% ET c! 0 1982 1986 1990 1994 1998 2002 2006 year 年代
Water Scarcity in the North China Plain Ground Water Table 2001.12
Improve Water Productivity in NCP Crop Yield(kg) Water Productivi ty (WP) 3 Water Consumption (m ) Biomass Harvest Index Evaporatio n Transpiration Option1: Increase crop yield by using same quantity of water Option2: Reduce water consumption while maintain the yield
Improve Water Management in NCP ET a = T r + E s = P + I D + dw Soil evaporation (E s ) is a physically controlled flux, but plant transpiration (T r ) is strongly influenced by plant physiology and can also be affected by abiotic environmental condition E s can be a major flux in sparse vegetation (low LAI) or wetting soil surface after irrigation/precipitation E s is non productive water use and accounts 1/3 ~ 1/2 of actual evapotranspiration (ET a ) which depend upon the water management practices Need to understand the E s regime and means to control
Conventional approaches for partitioning ET a Combination of soil lysimeters for evaporation (E s ) and sap flow sensors/chambers for plant transpiration (T r ) - poor spatial representation Eddy covariance system /Bowen ratio energy balance system / weighting macro-lysimeter for evapotranspiration (ET a ) and soil lysimeters for evaporation (E s ) or sap flow sensors/ chamber system for plant transpiration (T r ) - scale transformation, fetch length Theatrical methods such as Shuttleworth-Wallace model, dual crop coefficients method and time series analysis method parameters uncertainty
Incorporating measurements of isotopic concentration of water in soil and plant and air vapor as a tracer can overcome the limitations of those conventional methods ET partitioning from isotopes of canopy vapor δ ET using the Keeling plot approach δ E estimated by Craig-Gordon model
ET partitioning from isotopes of canopy vapor ET a δ ET ET a = T r + E s 1=T r /ET a +E s /ET a T r δ T δ T = δ T T r /ET a + δ T E s /ET a δ ET ET a = δ T T r + δ E E s δ ET = δ T T r /ET a + δ E E s /ET a E s δ E F s = δ ET δ T δ E δ T =E s /ET a δ = (R sample /R standard 1) * 1000
δ ET using the Keeling plot approach Keeling relationship for water vapor Slope (m) Intercept (b)
δ E estimated by Craig-Gordon model E L V S h V L V S h V L V 1 h /1000 1 h Δξ is an isotopic diffuse coefficient; δ S δ V is the isotopic composition of liquid water at the evaporating front; the isotopic composition of the background atmospheric water vapor; ε L-V the temperature dependent equilibrium fraction factor α L-V is (1- ε L-V )X1000; h the relative humidity normalized to the temperature of the soil surface
Shortcomings Traditional cold-trap method is time consuming and labor-intensive, and has limited most studies to short period (several days), small scale (chamber scale), and low time precision (daily) There were some differences between the measurements and real values because of sampling pollution, isotope fraction from condensation, and assumption for isotopic steady state in soil plantatmosphere
Water isotope analyzer with liquid water injector Field deployable water vapor isotope analyzer Real-time and continuous measurements of 18 O and D in air vapor and liquid water by tunable diode laser absorption spectroscopy provide an opportunity to perform in situ and continuous evapotranspiration partitioning on diurnal timescale
Objectives: Assess the accuracy of isotopic method in partitioning ET over irrigated wheat field in North China Pain to guide the E s management Partitioning E and ET by using mini-lysimeters and eddy covariance system Partitioning E and T using in situ measurements by air vapor isotopic analyzer Comparing the results estimated by conventional methods (EC + mini-lysimeter) with the ones done by isotopic methods
Materials and Methods 昌平
Sprinkler irrigation Rainfall 75 mm 45mm Soil moisture dynamics and irrigation
EC System Picarro System BR System EC/BR T a R n R a V w P RH T s G
Soil evaporation-mls Soil water content profile
160cm 100cm 80cm 30cm 5cm 2 H / 18 O air vapor isotope analyzer
Isotopic composition analysis of water in soil and plant samples Soil & plant sampling Extracting liquid water Purifying liquid water Isotopic composition analysis
Results Canopy cover variation and fitted curve during experimental period
18 O/ 16 O hour scale dynamics with rainfall and irrigation and weather parameters P+I The results showed that δ 18 O composition of air vapor at two heights correlated significantly with VPD and R n with mean correlation coefficients about 0.696 (n=1250,α<0.001) and 0.704 (n=1250,α<0.001).
The relationship between δ 18 O and 1/vapor H 2 O content Jointing Booting Filling Maturing
Sprinkler Irrigation 18 O/ 16 O D/H
O18/O16 discrimination ( ) O18/O16 discrimination ( ) O18/O16 discrimination ( ) O18/O16 discrimination ( ) δ ET -The Keeling plot 1/Vapor H2O 0.00E+00 4.00E-05 8.00E-05 1.20E-04 1.60E-04 0.0 1/Vapor H 2 O 0.00E+00 2.00E-05 4.00E-05 6.00E-05 8.00E-05 1.00E-04 1.20E-04-6.0-5.0-10.0 y = -52591x - 12.273 R² = 0.418 n=75 <0.001 9-May -10.0-15.0 y = -48006x - 13.22 R² = 0.405 n=75 <0.001-14.0-20.0-18.0-25.0-22.0 1/ Vapor H2O 2.00E-05 6.00E-05 1.00E-04 1.40E-04 1.80E-04 0.0 1/ Vapor H2O 3.00E-05 4.00E-05 5.00E-05 6.00E-05 7.00E-05 8.00E-05 9.00E-05-5.0-10.0 y = -42308x - 15.72 R² = 0.653 n=85 <0.001 3-Jun -10.0-14.0 y = -87843x - 10.7 R² = 0.813 n=80 <0.001 12-Jun -15.0-18.0-20.0-22.0-25.0-26.0-30.0 2012/7/26 25-30.0
ET_EC + Es_Micro-lysimeter
Fs-SI The relative contribution of soil evaporation to evapotranspiration (Fs) estimated by eddy covariance micro-lysimeters method and its comparison with isotopic method during the experimental period 1.00 1.0 E T / ET F T 0.90 0.80 0.70 0.60 ISS 法 MLS 法 0.8 0.6 0.4 0.2 y = 1.1967x R 2 = 0.8468 0.50 17/Apr 27/Apr 7/May 17/May 27/May 6/Jun 16/Jun Date 0.0 0.0 0.2 0.4 0.6 0.8 1.0 Fs-EC_MLS
The relationship between Es/ET and canopy cover (Cc)
The ratio of plant transpiration to evapotranspiation (Tr/ET) Crop sowingtillering tilleringjointing tilleringheading floweringfilling Fillingmature Total Wheat 0.58 0.32 0.76 0.85 0.78 0.71 Maize 0.33 0.72 0.82 0.80 0.68 30% of Evaporation!
Conclusions There was a good consistent between the estimated Es/ET by the stable isotopic method and ones by the conventional method, indicating that combination of Keeling plot method with in situ continuous measurements of water vapor stable isotope composition can accurately partition evapotranspiration in wheat field of NCP Managing evaporation in the low LAI period is the effective solutions to reduce the non-productive water use while improve water productivity, eg., mulching, insufficient irrigation, etc.
Thank You! 2012/7/26 31