Approaches in modelling tritium uptake by crops

Size: px

Start display at page:

Download "Approaches in modelling tritium uptake by crops"

Edwin Flowers
6 years ago
Views:

1 Approaches in modelling tritium uptake by crops EMRAS II Approaches for Assessing Emergency Situations Working Group 7 Tritium Accidents Vienna January 2010 D. Galeriu, A Melintescu

2 History Different models and equations have been proposed to express the uptake kinetics of tritiated water.the first is C TFWT :HTO concentration in the plant at the considered time t (Bq L -1 ) C : steady-state TFWT concentration (Bq L -1 ) k : rate constant for HTO uptake (h -1 ) t : time after the beginning of exposure (h) But C =1.1*ρ a / ρ s C ah ρ s is water vapor density in leaf stomatal pore (g /m3), ρ a is the water vapor density in atmosphere (g /m3), C ah is the air water HTO concentration (Bq/L) k = ρ s /(1.1*W*r) W water content of leaf (g /m2), r leaf resistance to water transport (h/m) The above relationships were used to interpet experimental dat aon various plants and environmental conditions. Many results will follow

3 Atarashi 1997

4 From Ichimasa

5 From Ichimasa Other values in Cecile Boyer thesis and paper

6 Mesures dans l eau l tissulaire : conditions d éd éclairement C HTO laitues (Bq L 1 ) / C HTO air (Bq L 1 ) 1,20 1,00 0,80 0,60 0,40 0,20 C t HTO laitues 1 / 2 = HTO = C α (1 e ln(2) k air k.t t1 / 2 = 2,9 h, α = t1 / 2 = 22,4 h, α = ) t1 / 2 = 1,5 h, α = 0,42 matures 0,21 0,43 témoins prémontaison jeunes 0, Durée de l'exposition (h) 7 e Congrès National de la SFRP juin Angers 6

7 Mesures dans l eau l tissulaire : conditions d éd éclairement C HTO laitues (Bq L 1 ) / C HTO air (Bq L 1 ) 1,20 1,00 0,80 0,60 0,40 0,20 C t HTO laitues 1 / 2 = HTO = C α (1 e ln(2) k air k.t t1 / 2 = 2,9 h, α = t1 / 2 = 22,4 h, α = ) t1 / 2 = 1,5 h, α = 0,42 matures 0,21 0,43 témoins prémontaison jeunes 0, Durée de l'exposition (h) 7 e Congrès National de la SFRP juin Angers 7

8 Rate constant k shows a large variability between plants and environmental conditions. Clearly depends on light, temperature, humidity and development stage of plants We must asses the uptake by the vegetation canopy, not for a single leaf Keum use a single value for morning, all plants, Gazaxi (2002) use single values for day and night ETMOD (1994) use seasonal value of leaf resistance by macro plants categories (binome) UFOTRI scale leaf resistance to canopy by dividing leaf resistance to leaf area index In land atmosphere interaction, exchange velocity is used (inverse of resistance) due to atmospheric resistance, boundary layer resistance and canopy resistance Follows excerpts form a lecture last year (A Melintescu)

9 Resistance Approaches to Deposition and Exchange Total Surface, R c Atmospheric source Aerodynamic, R a Boundary, R b Stomatal, R s Cuticular, R ct Ground, R g for various surfaces Similitude between water vapour transport and electric circuits, because in both cases the transport is due to specific gradients: - specific humidity for water - electric potential for electricity Resistance to environmental transport is defined by analogy with resistance in electric circuits, both of them being the ratio between potential difference and flux Aerodynamic resistance R a depends on turbulence and wind speed Boundary layer resistance R b depends on turbulence, wind speed and surface properties Total surface resistance R c can be split up into canopy and ground related resistance Canopy resistance depends on surface properties, temperature, photosynthetically active radiation (PAR), humidity, water content in soil For HT deposition, ground resistance depends on the rates of diffusion and oxidation in soil, and is much lower than the canopy resistance Deposition velocity=1/(r a +R b +R c ) This is also an exchange velocity at air to plant (soil) interface

10 Turbulent eddies are responsible for transporting material through the surface boundary layer Transport processes associated with the transfer of heat, mass and momentum modify the properties of the the atmosphere. A distinct aspect of the boundary layer is its turbulent nature. Momentum must be transferred downward. A force is needed to change momentum transfer from one level to another. This drag force or shear stress is also equivalent to the momentum flux density Logarithmic wind profile Visualization of momentum transfer u* - friction velocity K von Karmann s constant (=0.40) z - height above the ground z 0 roughness parameter. It defines the effectiveness of a canopy to absorb momentum; valid only for very short vegetation and for a neutrally stratified atmosphere d - Zero-Plane Displacement Height. It represents the level at which surface drag acts on the roughness elements or level which would be obtained by flattening out all the roughness elements into a smooth surface.

11 Turbulent eddies are responsible for transporting material through the surface boundary layer. The aerodynamic resistance determines the rate that momentum, and other scalars, are transported between a given level in the atmosphere and the vegetation s effective surface sink. The aerodynamic resistance is expressed as: ψc - adiabatic correction function Surrounding the leaf and covering the surface of the soil is a thin skin of unperturbed air - the boundary layer Heat and water vapor must be transferred through this layer through molecular diffusion (conduction). The long timescale involved can be represented by a large resistance - the boundary layer resistance. The magnitude of this resistance depends mainly on the depth of the boundary layer and is proportional to leaf size/wind speed. z c - scalar roughness length, Sc - Schmidt number Pr Prandtl number. constant is often assumed to equal 2 over closed canopies, but can be much greater over rough incomplete canopies

12 R a, R b - affected by wind speed, crop height, leaf size, and atmospheric stability; - decrease with increasing wind speed and crop height Smaller resistances are expected over tall forests than over short grass and under unstable atmospheric thermal stratification, than under neutral and stable stratification. When wind speeds are 4 m s -1 theoretical boundary layer resistances over a 0.1 m tall grass, a 1.0 m crop and a 10 m conifer forest are about 60, 20 and 10 s m -1, respectively Experimental measurements show that both R a and R b are less than 20 s m -1 during the day over a temperate deciduous forest. Greater R a values (up to 150 s m -1 ) occur at night when turbulent mixing is reduced. Canopy resistance is predominant FOREST R a and R b vary between 4-18 s/m Surface resistance, mainly canopy, varies between s/m

13 Pojanie Khummongkol

14 Pojanie Khummongkol

16 Canopy resistance physiological models The canopy resistance (R c ) is a function of the canopy stomatal resistance (R stom ), the canopy cuticle resistance (R cuticle ), and the soil resistance (R soil ). These resistances are affected by leaf area, stomatal physiology, soil ph, and the presence and chemistry of liquid drops and films. The stomatal, leaf surface (cuticle) and soil resistances act in parallel, causing R c to be formulated as: Big-Leaf resistance models have electrical analogy - current flow (mass or energy flux density) is equal to the ratio between a potential and the sum of the resistances to the flow: C a concentration of a scalar in the atmosphere over the vegetation C 0 internal concentration

temperature q air air vapour in atmosphere Scalling from leaf to canopy -classic: R

17 Stomatal cavity common pathway for water and CO 2 Leaf = Σ stomata E qin = ρa r a q + r air c E evaporation ρ a air density q in saturated air vapour at leaf temperature q air air vapour in atmosphere Scalling from leaf to canopy -classic: R c = R leaf /LAI -big leaf: integral over all canopy as a single leaf -physiological approach

Jarvis approach light, temperature, water vapour deficit, and soil water deficit behave independently as modifying factors (0, 1) - minimal leaf resistance R c- min is plant characteristic

18 Jarvis approach light, temperature, water vapour deficit, and soil water deficit behave independently as modifying factors (0, 1) - minimal leaf resistance R c- min is plant characteristic Physiological approach link between water and CO 2 pathway to photosynthesis (A n ), taking into account different diffusion coefficients Ball-Berry scheme uses m and b as semi-empirical coefficients inconvenience

20 Leuning, improvement of Ball Berry Cs - the CO2 concentration at the leaf surface Ci - the CO2 concentration in the plant interior An - the net assimilation rate- leaf MOSES g l,c and g l,w are leaf conductance for CO2 and water vapor

Jacobs-Calvet-Ronda (preferred and tested) - assumes that C conductance is determined by ratio between photosynthetic rate and the concentration difference of CO 2 for leaf surface and leaf

CO 2 concentration in the plant interior f 0 - the maximum value of (C i - Γ )/(C s - Γ) D 0 - the value of Ds at which the stomata close Γ CO 2 compensation point g l,c leaf C conductance; g

21 Jacobs-Calvet-Ronda (preferred and tested) - assumes that C conductance is determined by ratio between photosynthetic rate and the concentration difference of CO 2 for leaf surface and leaf interior g min,c - the cuticular conductance A g - the gross assimilation rate- leaf D s - the vapour pressure deficit at plant level C -the CO s 2 concentration at the leaf surface C i -the CO 2 concentration in the plant interior f 0 - the maximum value of (C i - Γ )/(C s - Γ) D 0 - the value of Ds at which the stomata close Γ CO 2 compensation point g l,c leaf C conductance; g l,w leaf water conductance; g c,c C canopy conductance; g c,w - water canopy conductance For canopy - integrate on LAI We use gross canopy photosynthesis rate from WOFOST; Data base exist advantage

22 Water vapor deficit and soil water deficit Vegetation type Low vegetation C3 Low vegetation C4 Lobos Rice and phalaris grass Forest temperate Boreal forest f o a d (kpa -1 ) relative conductance 1.2 C3 teo 1 Do= Do=1 Do= VPD [kpa] stomatal conductance m/s stomatal conductance and humidity defficit -C3 and C4 grass humidity deficit g/kg g_c3 g_c4

23 Photosynthesis, at canopy level Many approaches in literature Need to considers sun and shaded leaves, nitrogen influence on photosynthesis rate, leaf orientation, leaf area profile etc. Scaling from leaf to canopy We simplify using WOFOST In land atmosphere interaction they use Ball Berry and Farquar models

24 NASA/LBA-ECO (CD36)

27 Canopy resistence 1.00E E+03 fgri fmai fpot fbet fwba fsba fwwh s/m 1.00E E time after 4 am, unit=0.5 h

28 atm+bound resist s/m hour of day Sum of atmospheric and boundary layer resistances

EXCHANGE VELOCITY APPROACH AND OBT FORMATION IN PLANTS DURING THE DAYTIME

EXCHANGE VELOCITY APPROACH AND OBT FORMATION IN PLANTS DURING THE DAYTIME Anca Melintescu PhD Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest - Magurele, ROMANIA ancameli@ifin.nipne.ro,