CR journee Capsis 12 juin 2003

5ème journée Capsis – 12 juin 2003 – Montpellier

Compte-rendu

F. de Coligny – 23 juin 2003

Participants (20) :

Philippe Ancelin (LRBB), François de Coligny (AMAP), Thomas Cordonnier (Cemagref), Guillaume Cornu (Cirad Forêt), Benoît Courbaud (Cemagref), François Courbet (INRA), Christine Deleuze (AFOCEL), Philippe Dreyfus (INRA), Pierre Dubus (Cirad), André Gavaland (INRA), François Goreaud (Cemagref), Sylvie Gourlet-Fleury (Cirad Forêt), Hélène Joly (CEFE), Isabelle Lecomte (INRA), Benoît Loussier (Cemagref), Céline Meredieu (INRA), Olivier Pain (AFOCEL), Sandrine Perret (Cemagref), Thomas Perrier (AMAP), Laurent Saint-André (Cirad Forêt), Patrick Vallet (LERFoB)

Excusés (grèves des transports) : Ariane Angelier (ONF), Jean-Claude Chopart (ONF)

Compte-rendu :

La 5ème journée Capsis s'est tenue au Cirad à Montpellier le 12 juin 2003. Cette réunion annuelle est un lieu d'échange et de discussion autour de la plate-forme de modélisation de la dynamique forestière Capsis. La réunion organisée par l'UMR AMAP s'organise autour de présentation des projets en cours ou prévus et laisse une large part aux discussions.

1. Capsis 4.1.2 – Bilan, avancées, perspectives – F. de Coligny

La première présentation a été l'occasion d'évoquer les avancées concernant la plate-forme dans son ensemble et qui concernent donc tous les chantiers d'intégration de modèles. Les points principaux :

Le bilan d'activité au mois de juin 2003 fait état d'environ 21 modules, 4 bibliotèques, 9 familles regroupant une centaine d'extensions, 2 pilotes (interactif, script), 4 sessions de formation (une quarantaine de personnes), la participation à 6 projets nationaux ou européens.

2. Pilotage par l'azote du modèle de croissance Eucalypt-dendro – L. Saint-André.

Le projet concerne plusieurs disciplines en interaction...

...et met en oeuvre plusieurs types de modèles :

  1. Ecophysiologique (lumière -> biomasse)

  2. Architecture (probabiliste)

  3. Dendrométrique (variables synthétiques, évolution)

  4. Cycles biogéochimiques (compartiments / flux)

Connexions :

Description du module Eucalypt-dendro : modèle piloté par la croissance en hauteur dominante (Hdom), fonction de (site, âge, densité). "Habillage" ultérieur : profil de tige, biomasse (tronc, houppier, racines, éléments nutritifs), épaisseur de l'écorce. Le modèle est intégré dans Capsis pour des plantations clonales correspondant à des ensembles de modèles.

Les forces du modèle : simplicité, entrées compatibles avec la gestion forestière, prise en compte de la compétition entre arbres. Faiblesses : ne simule pas les changements à long terme, pas d'influence des saisons, pas de rétroaction – bilan entrées / sorties en éléments minéraux.

L'intégration des cycles : le pas de temps va de l'heure à quelques jours, des éssais très instrumentés et suivis sont disponibles, adapté pour les modèles à base écophysiologiques.

Etude du retour des éléments nutritifs au sol (chute de litière, turn-over racines), décomposition de la litière, conséquences sur la croissance.

Site d'étude : plusieurs zones : rémanents balayés, site en gestion courante et site avec deux fois plus de rémanents. En première analyse, la gestion des rémanents a un effet sur la croissance en hauteur moyenne.

Pour le couplage dans Capsis, la machinerie est déjà en place, pas de difficultés majeures attendues.

3. Genetics : bibliothèque dédiée aux informations génétiques et aux processus de reproduction. Intérêt pour la modélisation et la simulation conjointe des structures démographiques, dendrométriques et génétiques – C. Pichot, P. Dreyfus

La présentation a été illustrée sur le cas du module Ventoux de P. Dreyfus : retour du Hêtre et du Sapin pectiné dans des peuplements de Pins (noir, sylvestre, à crochets) de l'arrière pays méditerranéen. Il s'agit d'une situation assez fréquente dans l'arrière pays méditerranéen.

Les éléments abordés :

Pour plus de détails, consultez l'url du projet : couplagebrg.free.fr

4. Selva sous Capsis - G. Cornu

Le module Selva implémente le modèle Selva de forêt tropicale de Sylvie Gourlet-Fleury, qui a été concu grâce au dispositif de Paracou en Guyane.

Discussion : proposition pour une bibliothèque mathématique dans Capsis pour accueillir des outils de génération de nombres aléatoire ou d'ajustement non linéaire.

5. Une démarche d'utilisation du mode Script - T. Cordonnier

Contexte de l'étude : rôle de protection des forêts de montagne contre les risques naturels (avalanches, chutes de blocs). Recherche de structures stables des peuplements résistantes et pérennes dans le temps. Existent-elles ? Peut-on les maintenir ? Peut-on optimiser simultannément plusieurs fonctions pour la forêt ?

L'exposé s'appuie sur le module "mountain" (pessières de montagne) de Benoît Courbaud pour étudier les conséquences d'un régime de perturbations (prélèvement de tiges) sur la stabilité de la surface terrière, des structures et des critères de protection (élancement moyen / vent, somme des diamètres / chutes de blocs, surface des trouées / avalanches).

La recherche de tendances se fait par répétitions (pour la stochasticité), par exemple sur une surface d'un hectare, à partir d'un sol nu, sur 1000 ans, étude de n modalités sur k répétitions.

On utilise le mode script de Capsis pour piloter les répetitions et extraire des variables de sortie dans un fichier texte repris dans Excel pour traitement et représentation graphique. Le script est organisé en modalités > répétitions et on obtient un fichier par modalité. Des graphiques sont produits pour montrer l'évolution des critères surveillés.

6. Utilisation du mode script pour estimer la variablilité des simulations - F. Goreaud, B. Courbaud

Les raisons d'étudier la variablilité : Capsis fournit des prédictions en fonction d'un état initial et de paramètres. Quelles confiance avoir en ces résultats ? Pour répondre, on peut :

Les modèles de Capsis sont stochastiques, au moins généralement pour la mortalité et la régénération.

L'état initial peut être aléatoire (ex: à partir de données moyennes) et il est possible de simuler un peuplement virtuel, correspondant à un état stochastique. Le modèle donne des prédictions et les résultats diffèrent à chaque simulation, il faut donc évaluer la variabilité.

Le problème n'est pas simple car les résultats de simulation sont complexes (plusieurs variables, évolution - non ponctuel, stochasticité).

Comme outil, on peut choisir une approche intuitive et étudier les courbes d'évolution moyenne / écrat type pour un pas de temps et une variable ou bien utiliser d'autres pistes (intervalles de confiance, périodicité de déphasage, attracteurs...).

Le mode script a été utilisé pour jouer un grand nombre de simulations (au moins 100 répétitions pour étudier la moyenne et l'écart type avec précision), ce qui est difficile avec le mode interactif de Capsis, non adapté.

Un script boucle sur n répétitions : initialisation du peuplement, évolutions prévues pour la durée souhaitée, récupération des données objectifs pour obtenir une courbe d'évolution pour chaque simulation et chaque variable.

Résultats sur 1/4 d'hectare :

Sur 100 simulations, peu de variabilité rencontrée pour N : écart type maximum de 14, ce qui est peu, avec un écart type relatif maximum de 3,5%, ce qui est faible.

Une seule simulation suffit pour tirer des conclusions.

La variabilité de N est beaucoup plus forte : écart type maximum de 35, 35% en relatif.

En conclusion, on insiste sur l'importance d'estimer la variabilité. En projet : calcul en réseau sur un ferme de calcul et vers un calcul intégré dans Capsis ?

Note : ferme de calcul : réalisé par F. Goreaud depuis la réunion.

7. Projet Safe : modélisation de plantations agroforestières - I. Lecomte.

Rappel du contexte : Safe est un projet Européen concernant l'agroforesterie. Ses objectifs sont de mesurer les potentiels des pratiques, prédire l'impact sur l'agriculture européenne, fournir un cadre réglementaire.

Utilisation d'un modèle au pas de temps journalier :

Au pas de temps annuel :

Si possible, construire un modèle générique paramétrable pour plusieurs espèces, pour comparer des scénarios (élagage, éclaircie, rotations...), des géométries de plantes différentes et évaluer en aval les impacts économiques.

La parcelle type comporte des bandes d'arbres, des bandes non cultivées et des bandes de cultures intercalaires. La scène utile dans Capsis est découpée en cellules de terrain de 1 m2 et en couches souterraines correspondant à des horizons pédologiques différents. Elle est plantée d'arbres moyens dont on considère les voisins identiques.

Processus mis en oeuvre :

Actuellement, les points traités :

Discussion sur l'exportation de données, il serait possible d'utiliser un langage d'exploration de structures de données (XPATH, XML), Le Cemagref utilise un explorateur Sim Explorer qui propose des formats, des traducteurs, des parsers. Certains sont partisans d'en rester à l'exportation de fichiers ascii avec séparateurs (virgule, tabulation...).

8. Outils logiciels de l'AFOCEL : OASIS, CUBESTIM, AFOSYLV - O. Pain, C. Deleuze

L'intervention concerne l'approche modélisation de l'AFOCEL et les outils mis en oeuvre auprès des utilisateurs gestionnaires.

D'après les résultats de l'IFN, les ressources forestières en production (proche récolte) sont pour les futaies régulières feuillus de 38% et pour les futaies régulières résineux de 86%. Les 3 espèces importantes sont l'Epicea, le Douglas et le Pin maritime.

Les modèles sont construits à partir du réseau d'éssais AFOCEL, répartis sur toute la France avec de nombreuses essences et des mesures à la fréquence des éclaircies qui nourissent une base de données.

Les principes observés sont : des variables d'entrée accessibles au gestionnaire (hauteur dominante...), des formes contraintes et une orientation vers des modèles génériques.

CUBESTIM s'adresse à l'exploitant, il s'appuie sur des fonctions de forme pour estimer les volumes produits. Il donne des estimations avec intervalles de confiance. Il a été testé par l'ONF qui l'a jugé assez satisfaisant.

OASIS s'adresse au gestionnaire, il modéliseles accroissements en hauteur dominante et individuelle, l'auto-éclaircie, la sylviculture et jusqu'au bilan économique. Le moteur d'OASIS est un modèle Arbre Indépendant des Distances qui prédit à partir d'une distribution de circonférences à t0 l'accroissement en surface terrière et les hauteurs individuelles. Il utilise des scénarios prédéfinis ou construits pas à pas (âge, type d'éclaircie) et propose des sorties dendrométriques (hauteur dominante, effectif...).

La partie économie gère des échéanciers, des jeux de produits et permet la comparaison entre la vente en bloc avec prix au m3 avec la vente à l'unité de produit avec un prix par produit. Un module FOREXPER permet l'expertise d'une parcelle à un certain âge.

AFOSYLV (en cours de développement) vise à traiter une propriété (plusieurs parcelles). Il propose des comparaisons prédit / réel sur des mêmes graphiques pour aider à faire correspondre le théorique avec le réel. La démarche est importante pour le gestionnaire qui apréhende mieux ce qu'il fait en s'appuyant sur ce qu'il connaît. Le logiciel contient des tarifs de cubages français, il est possible d'en créer de nouveaux personnalisés. Une grande liberté est donnée concernant les tarifs, la fertilité et les courbes de prix. Les résultats de simulation sont obtenus au niveau parcelle, propriété ou pour plusieurs propriétés.

En conclusion, les modèles proposés concernent les essences majeures de production et sont des modèles génériques. Ils sont implémentés dans des logiciels concrets, ouverts, personnalisables et à disposition des acteurs de la filière bois.

9. Fagacées 2.0 et nouvel algorithme d'éclaircie utilisant le RDI - P. Vallet

Fagacées est un modèle de croissance (Jean-François Dhôte) pour le Chêne Sessile et le Hêtre. Il étudie des peuplements réguliers, non spatialisés, le moteur de la croissance est la hauteur dominante du peuplement, il s'applique à une surface de 1 hectare et utilise un pas de simulation de 3 ans.

L'évolution se calcule comme suit :

Algorithme d'éclaircie sur RDI (Relative Density Index = indice de densité relatif).

1. définir la population suceptible d'éclaircie

2. le taux d'éclaircie (strictement par la bas, standard, par le haut)

Fagacées 2.0 sous Capsis utilise une gestion individuelle des arbres et peut calculer jusqu'à 150000 arbres en régénération naturelle.

A l'initialisation, choix de l'espèce (chêne ou hêtre), régénération ou plantation, fertilité : (hauteur, âge, date) pour prendre en compte des modèles de changement de productivité (Jean-Christophe Hervé) et date de naissance (pour la même raison). Le modèle peut fournir des courbes d'évolution du volume en fonction de différentes dates de naissance du peuplement (1500, 1900, 2000 : 3 courbes) afin d'appréhender les changements de productivité.

Eclaircies automatiques : on considère une bande passante autour du RDI objectif. Quand le peuplement atteint la limite supérieure de la bande, l'éclaircie le ramène à la limite inférieure.

Le module Fagacées propose également un outil de profil de tige.

Note : depuis la réunion, le module Fagacées a été fusionné avec capsis4.1.2 et il est passé sous CVS. P. Vallet devrait pouvoir préparer sa diffusion pour évaluation à l'ONF dès cet été.

10. Le module Pin Maritime (PP3), l'extension Economics, le module Pin Laricio (PNL) - T. Labbé, C. Meredieu, C. Orazio, S. Perret

10.1 Le module Pin Maritime (PP3) - T. Labbé, C. Meredieu

Le module PP3 s'appuie sur les relations de Bernard Lemoine, il propose le mécanisme d'éclaircie par Densité et KG migré par Philippe Dreyfus depuis Capsis2 et dispose d'un visualiseur texte spécifique.

Dans le cadre du projet CapForêt, une connexion est en cours avec le logiciel ForestGales de la Forestry Commision (GB). Ce logiciel calcule des vitesses et probabilités de retour pour les chablis et volis sur le peuplement considéré. La communication est gérée depuis Capsis (par une extension ForestGales) et repose sur un échange de fichiers. Le fichier aller (.fgi) contient une ligne par étape du projet Capsis-PP3, il est passé en entrée de l'exécutable ForestGales (une "dll" tournant sous windows) et le fichier résultat est interprété et "tartiné" dans le projet de départ pour établir ensuite des graphiques.

On prévoit d'ores et déjà de passer non plus les peuplements mais les diamètres d'arbres individuels, de prendre en compte la profondeur de l'alios et les différences d'ancrages des arbres de bordure ou à l'intérieur du peuplement.

Pour ce qui est de la diffusion, l'objectif est de fournir une version de démonstration de PP3 à des coopérateurs du GIS (ONF, CPFA) à l'automne 2003, d'avoir une version avec licence disponible sur internet début 2004 et de mettre au point une formation Capsis et PP3 pour le premier trimestre 2004.

En ce qui concerne les améliorations à apporter, on prévoit :

10.2 L'extension économique dans Capsis - T. Labbé, C. Meredieu, C. Orazio

Etudier les réactions du propriétaire forestier par rapport à la tempête et évaluer les traductions économiques.

Améliorations envisagées :

Problèmes en suspend :

10.3 Le module Pin Laricio (PNL) - T. Labbé, C. Meredieu, S. Perret

Le modèle Pin Laricio implémenté dans Capsis est basé classiquement sur la croissance en hauteur dominante. Le module Laricio de Capsis2 (vesions précédente) était très complet : croissance, mortalité, branchaison, profil de tige, tarif de cubage, effet de la densité sur Hdom et liaison avec WinEpifn (évaluation de la qualité de la ressource sur pied). Dans Capsis4, le module n'est pas encore aussi complet.

A l'initialisation, le module propose la création de peuplement virtuel comme dans Capsis2 (nombre de tiges, surface, fertilité). On a rajouté la possibilité de définir la hauteur dominante en fonction de la station (T. Labbé). Il est également possible de charger un peuplement réel par un procédé semblable au module PP3.

Reste à faire :

Le module Laricio est de constitution simple et robuste, il est le fruit de la collaboration INRA / Cemagref qu'il est important de poursuivre.

11. Le module Cèdre dans Capsis - F. Courbet

C'est un modèle arbre indépendant des distances (MAID) concernant des peuplements purs réguliers de Cèdre de l'Atlas du quart sud-est de la France.

Pour la croissance, le module Capsis4-Cèdre comporte :

Il manque encore la croissance individuelle en hauteur.

Pour la branchaison, les propriétés suivantes ne sont pas encore incluses dans le module Capsis, amis sont connues :

Concernant les tiges : propriétés non encore incluses mais connues :

En entrée, on peut choisir de charger un peuplement virtuel ou un fichier inventaire.

En perspectives :

L'objectif est de fournir le module au gestionnaire pour aider à mettre au point des modèles de sylviculture.

Discussion : quels avantages offre le passage d'un MAID (arbres types = classes, avec un effectif) à un MADD (modèle arbre individuel) ?

Un argument est d'éviter les problèmes rencontrés avec les limites de classes dans les MAID. Les MADD sont plus lourds parce qu'ils font gérer plus d'objets, mais ils sont plus simples dans la démarche.

Il est noté qu'on peut éviter le problème de limites de classes avec des classes plus fines et que les MAID proposent des présentations plus synthétiques, adaptées au niveau forestier. Dans l'exemple du projet Ventoux (2800 ha étudiés), il aurait été bien trop lourd d'utiliser des MADD.

L'étude de la branchaison et des cernes se marie mieux avec une approche individuelle.

12. Fiesta : modèle de croissance et de branchaison du Pin d'Alep au service de la gestion des coupures de combustible - P. Dreyfus

Ce projet est en rapport avec le projet européen FireStar (coord. J.-C. Valette - INRA URFM Avignon) qui concerne l'interface forêt-habitat. Il s'agit de décrire et modéliser le combustible (arbustes, litières, arbres...) et sa répartition, d'utiliser des modèles physiques de comportement du feu et enfin des méthodes d'estimation du risque pour des "cibles" humaines ou bâtiments. Pour réduire ces risques, on crée, au niveau des coupures de combustible (pare-feux) des discontinuités horizontales (espacer les arbres) et verticales (débroussaillage et élagage pour éloigner le houppier du sol).

L'implémentation dans Capsis4 concerne la croissance et branchaison du Pin d'Alep dans le but d'étudier la gestion des coupures de combustible.

Le modèle de croissance utilisé pour le Pin d'Alep est un MAID avec branchaison (relations allométriques) et profil de tige.

Il existe deux modules Pin d'Alep dans Capsis : NRG et Fiesta. Des lacunes subsistent :

Démonstration : les modules NRG et Fiesta, avec un mécanisme d'éclaircie par espacement et un visualiseur montrant la branchaison. Les branches (diamètre à l'insertion) sont calculées à la volée lors de la visualisation, elles ne sont pas stockées dans l'arbre.

13. Elaboration d'un guide territorial Douglas du Grand Massif Central - A. Angelier

Ariane Angelier n'ayant pas pu venir à cause des grèves de transport, elle a envoyé sa présentation qui a été montrée à l'assemblée. Les grandes lignes sont présentées ici.

Il y a un intérêt fort d'Auvergne-Limousin pour un guide Douglas. Une enquête a été menée auprès des DT de l'ONF et un groupe de travail a été créé (Auvergne-Limousin, Bourgogne-Champagne-Ardennes, Méditerranée).

La première réunion (11 avril 2003) a traité de la notion de "guide biogéologique" s'appuyant sur la notion d'entité bio-géographique. Des objectifs communs ont été définis : établir un guide de sylviculture "complet" pour répondre aux attentes suivantes :

Une deuxième réunion s'est tenue le 19 mai pour illustrer l'intérêt d'une mutualisation, autour de la régénération naturelle, la présentation des résultats de la recherche menée par les STIR sur 8 sites et des discussions avec les représentants des Agences et UT.

Les thèmes déjà en partie traités :

Les thèmes restant à traiter :

désignation et élagage,

éclaircies (normes, rattrapage).

Rappel de l'objectif : plan d'action et délais précis pour 3 DT fin 2003 / début 2004 pour début de la diffusion et de la formation.

Conclusions

Depuis 1999, date de reprise du projet par AMAP, l'accent a été mis sur les relations avec les modélisateurs. Fin 2002, une réunion à Paris a permis d'aborder avec des gestionnaires la relance de l'utilisation de Capsis pour leurs besoins (compte-rendu sur le site Capsis https://membres.lycos.fr/coligny/).

Des groupes de travail ont été proposés pour certains modules et des réflexions sont en cours concernant leur diffusion pour évaluation. C'est notamment le cas pour le groupe Pin (diffusion du module PP3 Pin Maritime – C. Meredieu, T. Sardin (ONF)) et le groupe Chêne-Hêtre (diffusion du module Fagacées – J.-F. Dhôte, P. Vallet, M. Legay (ONF)).

Note : dans cette optique, les deux projets PP3 et Fagacées sont passés prioritairement sous le système de gestion de versions CVS depuis la réunion du 12 juin. Il est à présent possible pour C. Meredieu et P. Vallet de finaliser et d'assembler des versions d'évaluation pour les gestionnaires en supprimant les modules non encore diffusables.

F. Courbet a souligné que les modules Cèdre et Ventoux doivent être prochainement évalués par J.-M. Courdier et B. Courbaud a ajouté que le module Mountain est déjà chez Xavier Gauquelin avec lequel il co-encadre un stagiaire.

L'idée de formations couplées à l'utilisation de Capsis4 et de modules semble faire son chemin. Les gestionnaires semblent demandeurs. Elles pourraient avoir lieu à l'occasion de la diffusion des versions d'évaluation aux gestionnaires.

Des discussions sur la gestion des priorités dans le projet, qu'est ce qui doit passer avant ? Jusqu'à présent, les chantiers impliquant plus d'acteurs ont été priorisés en travail de fond (formations, développement concernant la communauté - noyau, interface graphique, extensions communes, mode script, gestion des graphiques, des groupes, etc - interventions pour les projets Démo-génétique, Safe, AMAPsim, CapForêt...). Les demandes ponctuelles (courtes) sont traitées rapidement (par interruption) par téléphone, par mail ou rendez vous à très court terme (problème technique, demande de conseil, demande d'évolution, assistance stagiaire). Les nouveaux projets sont intercalés dans le planning au plus tôt (sessions de travail en collaboration de 1 à 3 jours, le modélisateur repart avec son module / sa bibliothèque en autonomie). Les tâches les moins prioritaires sont les demandes individuelles nécessitant de la programmation (ajout de fonctionnalités dans un module spécifique ex: plusieurs espèces, correction d'erreurs non bloquantes avec solution de contournement...).

Il semble intéressant de viser à montrer des données réelles au milieu de données simulées (domaine de calibration...) pour que l'utilisateur puisse évaluer les résultats de simulation par rapport aux mesures réelles (détection visuelle d'un résultat hors calibration). L'option domaine de calibration, quand disponible, pourrait être associée à une légende (données qui on servi à calibrer le modèle, données de la propriété X...).

Le principe de la réunion annuelle apparaît satisfaisant et il est reconduit.