Actuellement en stage dans un Labo de géographie travaillant sur un modèle de la dynamique de la végétation, j'ai besoin de calculer la proportion de racine présente dans un sol pour des herbacées et des arbres en fonction de la profondeur du sol. Toutefois, cela fait appel à des notions mathématiques que je n'ai qu'approcher vaguement et il y a de ça très longtemps. Ainsi, si quelqu'un pouvais m'aider à résoudre ce petit problème cela me ferais beaucoup avancer.
Je vous pose ci-dessous les bases du problème :
- Code: Tout sélectionner
La densite racinaire du modèle decroit exponentiellement avec la profondeur.
On peut ecrire :
D(z) = EXP(-C * Z)
D etant la densite de racines a la profondeur Z.
C varie en fonction du PFT (type de végétation du modèle), du PFT 1 (mais
pour le sol nu cela n'a pas beaucoup de sens, au PFT 13,
C = 5., .8, .8, 1., .8, .8, 1., 1., .8, 4., 4., 4., 4.
L'unite de C est : metre a la puissance -1
Donc en gros C = 0.8 a 1 pour les arbres, et C=4 pour les herbacees.
Quand tu fais l'integrale de EXP(-C * Z) * DZ entre 2 niveaux de
profondeur, cela te donne la proportion de racines dans cette couche.
Entre 0 (la surface) et l'infini tu as le total des racines disponibles.
Tu dois donc pouvoir rapidement calculer la proportion de racines
consideree pour les arbres et les herbacees en fonction de la profondeur
totale de ton sol.
