Post on 06-Jan-2016
description
Mortalité progressive des talles :
0.0 0.5 1.0 1.5 2.0
02
46
8
stress hydrique moyen
dis
trib
utio
n
0.0 0.2 0.4 0.6 0.8 1.0
0.00.2
0.40.6
0.81.0
x/10
pgam
ma(x,
shap
e = 4)
distribution des transpiration des apex
nous avons choisi une loi gamma :
k : paramètre de forme θ : paramètre d’échelle x : seuil de mortalité
T/Tp =
espérance = k θ variance = k θ²
la fréquence de mortalité dépend du seuil la valeur du seuil est observée par les écophysiologistes, les valeurs diffèrent selon la sensibilité à la sécheresse
par conséquent, la forme de la distribution dépend de la variance de la population nous devons donner une relation entre la variance et la moyenne (kθ
nous avons testé comme constante la variance ou le coefficient de variation et choisi la variance
Propriétés
Définitions
0
1
2
3
4
0 50 100 150 200 250 300 350 400 450 500Days from 1/1/2003
0
1000
2000
3000
4000
Den
sity
(nb
tille
rs /
m²)
0
1
2
3
4
0 50 100 150 200 250 300 350 400 450 500Days from 1/1/2003
T/Tp
, LAI
0
1000
2000
3000
4000
0
1
2
3
4
0 50 100 150 200 250 300 350 400 450 500Days from 1/1/2003
0
1000
2000
3000
4000
Sensibilité de la mortalité de talles aux paramètres et au milieuRalentissement de la consommation d’eau quand l’humidité moyenne du sol (teta) est inférieure à tetstomate (Brisson et al 1998)
ZRACCUMLRACZ
RAYONPSISTOCUMLRACZ
EOPTETSTOMATE
1ln
102ln
80
14
INTRODUCTION DE LA DYNAMIQUE DE TALLES POUR MODÉLISER LA MORTALITÉ ESTIVALE DES PRAIRIES EN CAS DE SÉCHERESSE FORTE
Françoise Ruget, Sylvain Satger, Florence Volaire, François Lelièvre, Amirouche Abdessemed, S.amuel Buis
Durée de sécheresse et seuil de mortalité : évolution du nombre de talles
TETA <
0
1
2
3
4
0 50 1001502002503003504004505000
1000
2000
3000
4000
0
1
2
3
4
0 50 100 150 200 250 300 350 400 450 500Days from 1/1/2003
0
1000
2000
3000
4000
Den
sity
(nb
tille
rs /
m²)
0
1
2
3
4
0 50 100 150 200 250 300 350 400 450 500Days from 1/1/2003
0
1000
2000
3000
4000
début sécheresse 01/05 01/05 01/05
fin sécheresse 15/08 08/09 15/09
ETP correspondante
690 820 870
seuil de mortalité T/Tp
0.02 0.03 0.04
0 1000 2000 3000 40000
5000
10000
Montp
ellier
2003
sol 102
0
5000
10000
sol 109 0 1000 2000 3000 40000
5000
10000
sol 113
0
5000
10000
sol 126
0 1000 2000 3000 40000
5000
10000
Montp
ellier
2006
sol 102
0
5000
10000
sol 1090
5000
10000
sol 1130
2000
4000
6000
sol 126
0
1000
2000
3000
Buga
c 200
3
sol 1020
5000
10000
sol 109 0 1000 2000 3000 40000
5000
10000
sol 1130 1000 2000 3000 4000
0
5000
10000
sol 126
0
1000
2000
Buga
c 200
6
sol 1020
5000
10000
sol 109 0 1000 2000 3000 40000
5000
10000
sol 113
0
5000
10000
sol 126
Densité
fermeture stomatique psistoenracinement : densité racinaire initiale densinitial production de longueur racinaire draclong déclenchement de sénescence racinaire debsenrac vitesse de progression des racines en profondeur croirac action racines : rayon
variable d’état paramètreclimat : ETP
variable de forçage
longueur racinaire totale cumlracz
profondeur maximale des racines zrac
profondeur maximale d’évaporation du sol
évaporation du sol
Distribution de la densité finale de Distribution de la densité finale de talles pour 4 lieux-années et 4 solstalles pour 4 lieux-années et 4 sols
1 2 3 40
20
40
60
psisto0.02 0.04 0.06 0.08
0
10
20
croirac10 20 30 40
0
10
20
draclong0.2 0.4 0.6 0.8
0
10
20
lvfront
350 400 450 500 5500
5
10
15
debsenrac0.02 0.03 0.04
0
10
20
30
rayon20 30 40 50
0
10
20
zesx0.5 1 1.5
0
10
20
densinitial(1)
0.5 1 1.50
10
20
densinitial(2)0.5 1 1.5
0
50
100
densinitial(3)0.5 1 1.5
0
20
40
60
densinitial(4)0.5 1 1.5
0
10
20
densinitial(5)
30 40 50 600
5
10
15
zrac0
usm9 Bugac 03 S102
Distribution des paramètres nécessaires pour obtenir une reprise entre le 21 septembre et le 31 décembreDistribution des paramètres nécessaires pour obtenir une reprise entre le 21 septembre et le 31 décembre
2 4 6 8 10 12 14
2
4
6
8
10
12
14
usm9 Bugac 03 S102
Distribution finale de densités : les paramètres font que toutes les valeurs de densités sont possibles, le milieu n’est pas déterminant
psisto, croirac, draclong, lvfront, debsenrac, rayon, zesx, densinitial(1),densinitial(2), densinitial(3), densinitial(4), densinitial(5), zrac0
conditions initiales
zrac0
densinitial (1) à (5)
Paramètres :
Corrélations entre paramètres = combinaisons conduisant à ces situations