Compromisator_ID · Module_4Carbone

Module carbone du Compromisator

Le module compare, LandUnit par LandUnit, l’état initial et l’état final d’un scénario d’usage des sols. Il calcule aujourd’hui la variation de stock de carbone du sol et de la biomasse aérienne ; les émissions des effluents, du transport, de la fertilisation et des intrants sont à intégrer depuis les prototypes Ocelet.

État au 6 novembre 2026. Les anciennes tables de stocks ont été supprimées le 3 novembre, et le module ne tourne pas avec CARBONE_LUC=true tant que les tables GIEC ne sont pas lues.

codé dans Ocelet prévu, ou calculé dans un prototype (CarboneGenerique, CarboneRealiste) et absent de Module_4Carbone table saisie, pas encore lue

Entrées, sous-modules et sorties

MODULES PRÉCÉDENTS DU RUNRotation des usagesModule_1LandUseProductionrendements, résidusAppariement5 passes, flux de MOTERRITOIREun million de LandUnits de 50 m, 249 995 haLandUnitfractions d’usage en m²Usages50 codes cartographiésSol et climat43 profils, SOCrefZonezones 1, 2 et 3Itinéraire techniqueapports, travail du solHaies4 composantessurfaces m²résidusapports de MOMODULE CARBONEStocks du sol et aériensΔC, CO₂ du solItinéraire techniqueF_MG, F_IÉlevage et effluentsCH₄, N₂O, NH₃TraitementCH₄, N₂O, NH₃TransportCO₂ fossileFertilisation et solN₂O, NH₃, NO₃⁻, CO₂AmontCO₂e importéBilan GES et azoteCO₂e, bilan Nfractions, profil,zone, ITKtonnes de matière, trajets, apports par fractionTables de paramètresGIEC, procédés, véhiculesSORTIES carbone_luc/landunits.csvΔC sol, CO₂, statutsurfaces_usages.csvha retirés, ajoutéstrajectoire.csvannée par annéebiomasse_aerienne_*ΔC aérien, sommeemissions_postes.csvgaz, poste, support, datebilan_ges.csvCO₂e et bilan N par scénariokg
Le module reçoit les fractions d’usage du territoire, les tonnes de matière organique (MO) et les trajets de l’appariement. Il se divise en huit sous-modules, dont un seul est codé.

Toutes les émissions suivent le même contrat. Le flux d’azote ou de carbone restant après une perte devient l’activité de l’étape suivante, ce qui ferme le bilan de matière du bâtiment jusqu’au sol. Chaque valeur exportée porte son niveau de référence (3 mesure locale, 2 référence régionale, 1 GIEC, 0 inconnu) et sa méthode, facteur ou modèle dynamique. Les formules en N × EF donnent de l’azote, converti en N₂O par 44/28.

Variable d’activitéquantité transmise par le modèlekg N, kg TAN, kg MS volatile,km, t, ha, kWhsupport spatial et dateMécanismefacteur par défaut du GIECfacteur local ou régionalou modèle dynamique(ALFAM2, RothC, modèle N₂O)Table de paramètresune ligne par processus, matière,système et climat, avec sa sourceÉmissiongaz en kg, par postepar support et par dateniveau de référence exportéBilan× PRG, N₂O 267, CH₄ 27CO₂ biogénique à partsomme par scénarioFlux restantN et C moins la perte,transmis à l’étape suivantedevient la variable d’activité de l’étape suivante
Code actuelRotationCarboneProductionfractionsPrévu sans appariement (I7)RotationProductionCarbonerendements, résidusPrévu avec appariement (I8)RotationProductionAppariementCarboneflux de MO, ITK
Position de l’appel du module dans un run, aujourd’hui et après les étapes I7 et I8 du plan du 3 novembre 2026.

Échelle du territoire

Les émissions amont sont celles des intrants importés, rattachées à la quantité achetée par un facteur d’analyse de cycle de vie. L’électricité est produite sur l’île, hors du système agricole. Le transport routier relie les offres de matière aux parcelles.

ATMOSPHÈREHORS ÎLEEngrais minéraux· CO₂ synthèse de NH₃· N₂O acide nitrique· CO₂ fret maritimeAliments du bétail· N₂O et CO₂ des cultures· CO₂ changement d’usage· CO₂ fret maritimeCarburant· CO₂ extraction, raffinage· CO₂ fret maritimeN, kgtLLA RÉUNIONCentrales électriqueshors système agricoleCO₂ par kWh du réseauSYSTÈME AGRICOLE SIMULÉÉlevages et installationsbâtiments, fosses, fumières,méthaniseurs, plateformesCH₄ NH₃ N₂OTransport routiertrajets de l’appariement,aller chargé, retour à videCO₂Parcelles (LandUnits)apports, résidus, cultures,stock de carbone du solNH₃ N₂O NO₃⁻ CO₂MOMOkWh des procédésengrais vers les parcellesaliments vers les élevagesdieselproduits exportés(N et C qui sortent)émis hors île

Sous-module transport

Le sous-module part des paires acceptées par l’appariement. Chaque paire donne la matière, la quantité en tonnes de matière brute et la distance routière d’un aller entre l’offreur et le receveur. La quantité divisée par la charge d’un véhicule donne le nombre de rotations ; chaque rotation parcourt l’aller en charge et, le plus souvent, le retour à vide. La consommation de gazole dépend du véhicule et de sa charge, et sa combustion donne le CO₂.

Paire de l’appariementtype_biomassematièreq_transpt de matière brutedist_kmroute, un allerid_s, id_roffreur, receveurRotationsquantité / chargepar véhicule, arrondicharge utile,masse volumiqueKilomètresen charge etretour à videpart du retourà videGazoleconsommation selonla charge, manutentionL/100 km à videet à pleine chargeCO₂combustiondu gazolefacteur decombustionTables de paramètres, par classe de véhicule et par matièreamont du gazole, vers le sous-module Amontpar paire,passe,matière,débouchéFond orange, émission calculée. Bordure en tirets, étape à coder. Bordure en pointillés, table à créer. Le transport ne retire ni azote ni carbone au flux de matière.

Trois postes voisins restent hors du sous-module. L’extraction et le raffinage du gazole, ainsi que le fret maritime des imports, relèvent de l’amont ; le gazole de l’épandage au champ relève de la fertilisation. Une paire dont l’installation est sur le site de l’offreur n’a qu’un coût de manutention. Dans la passe 5, la matière des imports débarqués au port n’est pas renseignée, ce qui empêche aujourd’hui de choisir le véhicule.

Champs de l’appariement lus par le sous-module

Chaque passe RP_<passe>_<graine>/ écrit un shp_paths.csv, une ligne par paire acceptée. Le sous-module lit les paires en mémoire une fois l’étape I8 faite, ou ce fichier en attendant.

ChampContenuUsage dans le sous-module
id_pair, id_s, id_rPaire, offreur, receveurClé de rattachement des émissions
type_biomasseMatière transportée (lisier porcin, compost, canne…)Choix du véhicule et de sa charge utile
type_usageDestination (SugarMill, Methaniseur, fertilisation40…)Rattachement de l’émission à un débouché
q_transpQuantité livrée, en t de matière brute (MB)Nombre de rotations
dist_kmDistance routière d’un aller, calculée par Dijkstra sur le réseauKilomètres parcourus
x_s, y_s, x_r, y_rCoordonnées de l’offreur et du receveurZone ERIS, carte des trajets

La quantité est en MB, comme dans le staging de la méthanisation qui somme q en MB (Staging/ST_Supplyer_manager.oclt:485). Le champ cout n’est pas lu.

Tables de paramètres du transport

transport_vehicules.csv porte une ligne par classe de véhicule (tonne à lisier, benne agricole, camion benne, semi-remorque, citerne). Ses colonnes sont la charge utile (t), le volume (m³), la consommation à vide et à pleine charge (L/100 km), la part du retour faite à vide (0 à 1), le gazole de chargement et de déchargement (L/t) et la source de chaque valeur. transport_matieres.csv associe chaque type_biomasse à une classe de véhicule et à une masse volumique (t/m³). Pour une matière peu dense (déchets verts, bagasse), le volume limite la charge avant la masse.

Les facteurs du gazole vont dans la table du bilan, en deux lignes séparées, la combustion pour ce sous-module et l’extraction avec le raffinage pour l’amont. Le prototype CarboneGenerique utilise 3,2 kg CO₂e/L amont inclus ; les deux parts sont à prendre dans la Base Empreinte de l’ADEME [à vérifier].

Calcul par paire

charge_tmin(charge utile, volume × masse volumique)
→
n = ⌈q / charge_t⌉rotations, arrondi par paire
→
kmn × dist en charge, n × dist × retour_vide à vide
→
L gazoleconso interpolée selon le taux de charge, plus manutention × q
→
CO₂L × facteur de combustion

L’arrondi des rotations se fait par paire, un véhicule ne mélangeant pas deux receveurs. Le taux de charge vaut q / (n × charge_t), et la consommation en charge s’interpole entre la valeur à vide et la valeur pleine selon ce taux. La quantité de gazole L passe aussi au sous-module amont pour sa part d’extraction et de raffinage. Chaque ligne garde la passe, la matière, le débouché, l’offreur, le receveur et la zone ERIS de l’offreur. Le transport ne retire rien au flux d’azote ou de carbone, les pertes en citerne fermée étant négligées.

NiveauMécanismeDonnées nécessaires
1Consommation fixe par km, retour à vide complet, comme le prototype (0,35 L/km)Distance, quantité, charge utile
2Consommation interpolée selon le taux de charge, retour à vide par filière, manutentionTable des véhicules complète
3Consommation selon la pente et le type de route, tournées qui enchaînent plusieurs receveursAltitude des cellules du chemin RasterDjikstra, calendrier des épandages

Le niveau 3 compte la surconsommation d’un tracteur chargé qui monte vers les Hauts. Le chemin raster est déjà calculé par RasterDjikstra, et l’altitude de ses cellules suffit. Les sorties prévues sont transport_paires.csv (rotations, km en charge et à vide, gazole, CO₂, niveau par paire) et transport_synthese.csv (somme par passe, matière et débouché, en kg CO₂ par t MB et par kg N livré, ce second ratio se comparant à l’engrais minéral évité).

Ordre de grandeur sur le run F01 du 5 novembre

Calcul hors modèle sur output/jcp_full_effluent_20261105/F01/runs/RB/results/RP_*_1/shp_paths.csv, pour situer les passes entre elles. Les hypothèses, non sourcées, sont 15 t par rotation pour les liquides (lisier, digestat liquide, vinasse) et 25 t pour les autres matières, avec retour à vide complet, 0,35 L/km et 3,2 kg CO₂e/L.

PasseClassePaireskt MBMt·kmRotationskm parcourusm³ gazolet CO₂e
1, canne vers sucreriessolide2251 172,425,5547 0062 048 962717,12 295
2, méthanisationliquide33860,20,164 08021 2357,424
2, méthanisationsolide82776,10,363 59837 86613,342
3, compostageliquide301106,50,347 17945 95516,151
3, compostagesolide1 736209,53,059 630266 76493,4299
4, épandage et autres débouchésliquide2 203430,71,0929 305148 47552,0166
4, épandage et autres débouchéssolide6 8941 015,35,4745 038469 368164,3526
5, imports depuis le portsolide4 6684 682,3236,00190 43619 222 5196 727,921 529
Total17 1927 753,0272,02336 27222 261 1447 791,424 932

La passe 5 fait 86 % du total, dont 3,76 Mt vers « Bruleur » à 50,8 km en moyenne et 0,50 Mt vers « Bruleur_Biogaz ». La nature de ces flux reste à vérifier avant toute interprétation. Les passes 2 à 4, qui portent les effluents, les composts et une partie des produits (fourrage, fruits et légumes), donnent 1 108 t CO₂e. Restent à vérifier aussi les charges et consommations par classe de véhicule avec leur source, et le sens de dist_km sur d’éventuelles routes à sens unique.

GazProcessus et momentVariable d’activitéNiveau 1, paramètreMécanisme plus précisStatut
Transport. Trajets de l’appariement
CO₂e du gazolecombustion du gazole, extraction et raffinage inclus
à chaque trajet, aller chargé et retour à vide
km par trajet × nombre de rotations (t transportées / charge utile)km × 0,35 L/km × 3,2 kg CO₂e/L, amont du carburant inclus, prototypeconsommation selon la charge, la pente et le véhicule (tracteur et tonne, camion)prototype
CO₂ des opérations au champdiesel de l’épandage
à l’épandage
t épandues0,7 L/t dans le prototypeconsommation par matériel et doseprototype
Amont. Intrants importés, hors île, et électricité de l’île, rattachés aux quantités consommées
CO₂e de fabrication des engraissynthèse de l’ammoniac, N₂O des usines d’acide nitrique, fret maritime
à la fabrication, hors île
N minéral complémentaire (besoin moins service organique)5,5 kg CO₂e/kg N dans le prototypefacteur par forme d’engrais issu d’une base d’ACV [à vérifier]prototype
CO₂e de l’électricitécentrales de l’île, hors système agricole
à la consommation des procédés
kWh consommés0,7 kg CO₂e/kWh, pas un facteur réunionnaisfacteur du réseau réunionnais ; le crédit d’export suit le moyen de production évitéprototype
CO₂e des aliments importésculture et transport des aliments du bétail
à l’importation
t d’alimentsnon paramétréchantier 15absent

Échelle de l’installation

Chaque contenant d’effluent émet selon son pool de matière et ses conditions. Le CH₄ vient de la matière sèche volatile sans oxygène, le NH₃ de l’azote ammoniacal (TAN) au contact de l’air, le N₂O des zones aérées où l’azote se nitrifie.

ATMOSPHÈREBâtiment d’élevagecaillebotis ou litièrepool TAN, MS volatileNH₃CH₄Lisierfosse, avec ou sans croûtepool TAN, MS volatile, N org.CH₄NH₃N₂O si croûteFumierfumière ou tas au champpool N org., C, TANN₂ONH₃CH₄Méthaniseurdigesteur fermé, gazomètrepool C → CH₄, N org. → TANCH₄ fuiteCH₄ imbrûléPlateforme de compostandain actif puis maturationpool TAN, N org., CNH₃N₂OCH₄Digestatpost-digesteur ou fossepool TAN élevé, C résiduelNH₃CH₄Compost mûrstock couvert ou au champpool N org. stableNH₃ faibleco-substratépandageépandagesur les LandUnitssur les LandUnitsLe pool est la partie de la matière qui produit le gaz. TAN, azote ammoniacal. MS volatile, matière sèche organique. N org., azote organique.Le N₂O des effluents se forme dans les zones aérées (croûte, fumier, andain). Le facteur du GIEC est nul ou faible pour un lisier sans croûte [à vérifier].
GazProcessus et momentVariable d’activitéNiveau 1, paramètreMécanisme plus précisStatut
Élevage et effluents. Cheptel, bâtiment, fosse à lisier, fumière
CH₄ entériquefermentation dans le rumen
en continu, sur l’île
effectif par espècenon paramétréchantier 15absent
CH₄méthanogenèse anaérobie de la matière organique volatile
en continu au bâtiment puis pendant tout l’entreposage en fosse ou en fumière, plus intense par temps chaud
MS volatile (VS) par matière et systèmeVS × B0 × 0,67 × MCF/100, MCF en % par système et climat, GIEC 2019 éq. 10.23 [à vérifier] ; le prototype prend une part fixe du CMCF mensuel selon la température, la durée d’entreposage et les vidanges de fosseprototype
NH₃volatilisation de l’azote ammoniacal (TAN)
au bâtiment puis en fosse ou en fumière, selon la surface exposée
TAN collectépart du N total (Frac_GasMS, GIEC) ; 5 % du TAN au bâtiment dans le prototypeflux selon pH, température, vent, couverture ou croûteprototype
N₂O directnitrification et dénitrification dans les zones aérées (croûte, fumier solide)
pendant l’entreposage en fumière ou en fosse avec croûte
N excrété par systèmeN excrété entrant dans le système × EF3 × 44/28, GIEC 2019 éq. 10.25 [à vérifier]effet de l’aération, des retournements, de la duréeprototype
NO₃⁻, N ruisseléécoulements des fumières et des fosses
pendant l’entreposage, à chaque pluie
N total entreposéN × fraction lessivée par systèmebilan hydrique de la fumière ou de l’aire non couverteabsent
Traitement. Méthaniseur, plateforme de compostage, tonnes admises par l’appariement
CH₄ de fuitepertes du digesteur, du gazomètre et des soupapes
en continu pendant l’exploitation
CH₄ produit (25 % du C entrant dans le prototype)CH₄ produit × taux de fuite, 3 % dans le prototypetaux par équipement, torchère, couverture du stockage de digestatprototype
CH₄, NH₃ du digestatméthanogenèse résiduelle, TAN accru par la minéralisation (50 % du N organique dans le prototype)
pendant le stockage du digestat
digestat brut, solide et liquidefacteurs de stockage du lisier appliqués au digestatMCF du digestat, couverture de la fosseprototype
NH₃, N₂O, CH₄ du compostvolatilisation en phase chaude, nitrification et dénitrification, poches anaérobies
pendant le compostage et la maturation
TAN, N et C entrantsNH₃-N 25 % du TAN, N₂O-N 1 % du N, CH₄-C 0,5 % du C entrant, prototyperetournements, humidité, rapport C:N, structurantprototype
CO₂ biogéniquerespiration du compost, CO₂ du biogaz et de sa combustion
pendant le traitement
C entrant40 % du C entrant au compostage, prototypecompté à part, hors bilan GESprototype

Échelle de la parcelle

Une LandUnit est une maille de 50 m qui porte plusieurs usages. Le module garde les fractions produites par la rotation ; seules les surfaces changées entrent dans le calcul, et le profil sol et climat reste identique entre les deux états.

État initial (id_ocs_shp)canne 3000,15 haprairie 4000,10 hascénarioÉtat final (id_ocs_postransi_merged)canne 3000,15 halégumes 5200,10 hamaille de 50 m, 0,25 ha nominal, clé x_yProfil sol et climat 15_5, identique dans les deux étatssol volcanique (VOL), tempéré chaud humide, SOCref 138 t C/ha sur 0 à 30 cmComposition nettecanne 300, maintenue 0,15 has’annuleprairie 400, retirée 0,10 ha− A × c prairielégumes 520, ajoutée 0,10 ha+ A × c légumesSupport réalloué 0,10 haSeules les surfaces changées entrent dans ΔC.L’usage dominant n’est pas utilisé.

Après l’étape I3, le stock d’une fraction vaudra SOCref × F_LU × F_MG × F_I ; le code actuel lit un stock par profil et usage dans une table supprimée. SOCref dépend du profil sol et climat, F_LU de la catégorie GIEC du code d’usage, F_MG et F_I de l’itinéraire technique (ITK), à 1 pour l’ITK nominal. Dans l’exemple, le passage de 0,10 ha de prairie (138 t C/ha) à des légumes (138 × 0,69 = 95,2 t C/ha) donne ΔC = −4,28 t C, soit 15,7 t de CO₂. Le facteur 0,69 est celui des cultures annuelles en climat tempéré chaud humide (GIEC 2019, vol. 4, tableau 5.5). L’écart se décomposera en effet de l’usage, effet de l’ITK et interaction.

Classe GIECTypes de sol SeqCOITropical secTropical humideTropical très humideTempéré chaud secTempéré chaud humide
VOL, sols volcaniques (andosols)Aa, A, Ag, Ga, P, défaut50707784138
HAC, argiles à forte activitéHi2140602464
LAC, argiles à faible activitéGh, H, Hh1938521955

SOCref en t C/ha sur 0 à 30 cm (GIEC 2019, vol. 4, tableau 2.3). Les neuf types de sols d’Allo (2019) se réduisent à trois classes. Pour les andosols de La Réunion, C@RUN (Todoroff et al. 2019) donne 145 à 176 t C/ha sur 0 à 30 cm, au-dessus des 50 à 138 t C/ha du GIEC.

Parcelle de 200 m × 200 m à l’échelledemi-haie, 4,94 %0,20 ha, sol de forêt et biomasse ligneusebande non cultivée, 5,76 %0,23 ha, prairie non travailléeculture sous influence, 21,24 %0,85 ha, sol de la culture, facteur 1culture hors influence, 68,06 %2,72 ha, sol de la culturePart de surface de chaque composanteCoupe depuis l’axe de la haie (m)02,55,517,530axeinfluence de 15 m depuis le bord de la haieHypothèse du 4 novembre 2026, haie de 5 m partagée entre parcelles voisines, bande non cultivée de 3 m, influence sur le sol jusqu’à 15 m du bord de la haie.400 m de haie par parcelle de 4 ha, soit 100 m de haie par hectare, appliqués à toute la surface des 27 codes avec haies.
Composantes des 27 codes avec haies, à l’échelle. Hypothèse de 100 m de haie par hectare appliquée à toute leur surface, soit 4 890 km de haies dans G3 (48 904 ha). La biomasse de la haie reste à trancher entre 9,4 t C/km (GIEC 2019, tableau 5.1) et une plantation de 100 t MS/ha, qui donne 2,5 fois plus.
ATMOSPHÈREVÉGÉTATIONBiomasse aériennestock C, Module_4CarboneRacineshors niveau 1CultureN absorbé, récolte exportéeSURFACEApport épandu (TAN, N org., C)Résidus restitués, pailleSOL 0 À 30 CMC organique du solSOCref × facteursN organiquerapport C:NNH₄⁺ammoniumNO₃⁻nitrateSOUS-SOL ET NAPPEminéralisationminéralisationnitrificationhumificationdécompositionNH₃volatilisationCO₂respirationN₂ON₂O, N₂, NOnitrificationdénitrificationlixiviationNO₃⁻ vers l’eauabsorptionexportrestitution
Coupe d’une fraction d’usage. Seuls le carbone organique du sol et la biomasse aérienne sont calculés aujourd’hui.
GazProcessus et momentVariable d’activitéNiveau 1, paramètreMécanisme plus précisStatut
Fertilisation et sol. Fraction d’usage de la LandUnit
NH₃ des apports organiquesvolatilisation du TAN après épandage
de 0 à 7 jours après l’apport, surtout dans les premières heures
N total et TAN épandus par fraction et dateN total × FracGASM, GIEC 2019 tableau 11.3 [à vérifier] ; part du TAN par produit dans le prototypeALFAM2 4.2.14, selon MS, pH, température, vent, technique et délai d’enfouissementprototype
NH₃ des engrais minérauxhydrolyse de l’urée et volatilisation
dans les jours qui suivent l’apport
N minéral par forme d’engrais10 % dans le prototype ; FracGASF, GIEC 2019 tableau 11.3 [à vérifier]fraction par forme (urée, ammonitrate), pH du solprototype
N₂O directnitrification et dénitrification dans le sol
dans les semaines qui suivent l’apport, en pics après les pluies
N appliqué, N des résidus, N minéraliséN × EF1 × 44/28, EF1 0,6 % organique, 1,6 % minéral en climat humide, 0,5 % en climat sec, GIEC 2019 tableau 11.1modèle selon l’humidité du sol, la température et le C labileprototype
NOₓnitrification
après l’apport
N appliqué0,5 % dans le prototypeprototype
NO₃⁻ lixiviédrainage de l’eau en excès
pendant la saison des pluies
toutes les entrées de N (engrais, apports, résidus, N minéralisé)FracLEACH, 24 % en climat humide, 0 sans excès d’eau, GIEC 2019 tableau 11.3bilan hydrique mensuel (pluie moins évapotranspiration), réserve utile, enracinementprototype
N₂O indirectredépôt du NH₃ et dénitrification de l’azote lixivié
après la perte
N volatilisé, N lixivié× EF4 (1,4 % humide, 0,5 % sec) ou EF5 (1,1 %), × 44/28, GIEC 2019 tableau 11.3prototype
CO₂ de l’urée et du chaulagedissolution de l’urée et des carbonates
à l’apport
urée, CaCO₃ apportésurée × 0,20 × 44/12, GIEC 2006 ch. 11 ; chaulage non configuréprototype
N₂O des résidusdécomposition des résidus restitués (paille de canne)
pendant la saison qui suit la récolte
N des résidus aériens et racinaires laissés au solN × EF1 × 44/28biomasse et teneur en N par culture, part exportéeprototype
CH₄, N₂O du brûlagecombustion des résidus avant récolte
au brûlage
masse de résidus brûléeA × M_B × C_f × G_ef × 10⁻³, C_f 0,80 pour la canne, GIEC 2019 vol. 4 éq. 2.27 ; N brûlé retiré du N des résiduspart de la canne brûlée par zoneabsent
CO₂ du solvariation du stock de carbone organique
sur la durée d’ajustement D
ΔC par fractionmoins ΔC × 44/12RothC mensuel (prototype CarboneSolRothC)codé
N₂O de la minéralisationperte de carbone du sol qui libère de l’azote
sur la durée d’ajustement D
C perdu par anC perdu / R × EF1 × 44/28, R 15 après conversion, 10 pour une culture qui reste culture, GIEC 2019 éq. 11.8 [à vérifier]rapport C:N mesuré par type de solabsent
N₂O, NH₃ du pâturageurine et bouses déposées au champ
pendant le pâturage
N excrété au pâturageN × EF3PRPtemps de présence au champ par espèceprototype

Échelles de temps

Le pas de temps d’un sous-module se règle sur la durée de son processus. Le stock du sol s’ajuste linéairement sur D ans à partir de l’année de conversion t0, f(t) = min(1, max(0, (t − t0) / D)), avec D = 20 ans et un horizon de sortie de 20 ans par défaut. Le ΔC aérien est aujourd’hui compté en une fois, alors que les ligneux mettent environ 20 ans à atteindre leur stock.

1 h1 j1 sem.1 mois1 an10 ans20 ansdurée du processus après l’événement déclencheur (échelle logarithmique)NH₃ après épandageFertilisation et solNH₃, CH₄ au bâtimentEffluentsN₂O après apportFertilisation et solCompostageTraitementCH₄, N₂O en fosse ou fumièreEffluentsLixiviation du nitrateFertilisation et solAjustement du stock du solStocksCroissance des haiesStocksLe transport est instantané à cette échelle.
Ordres de grandeur, à préciser avec les données locales.

Plan de fin de développement

Le plan du 6 novembre 2026 couvre les huit sous-modules et reprend les étapes I et P du plan du 3 novembre. Tout le calcul se fait en Ocelet, dans Module_4Carbone ; les scripts Python ou R se limitent à préparer les tables et à lire les sorties. Un sous-module se compose d’une structure Ocelet, de tables CSV sourcées et de contrôles natifs. Chaque valeur exportée porte son poste, son gaz, son support, sa date, son niveau de référence et sa source, et une référence absente donne une valeur inconnue.

contrainte entre F_I et le N₂O des mêmes apports0Socle communEmission, lecteur, bilans N et C1Stocks rebranchéstables GIEC, haies, niveauxarbitrages des haies3Effluents, traitementport de CarboneGeneriqueGIEC ch. 10 et 11 absents4Transportpaires de l’appariementmatière « ? » en passe 55Fertilisation et solALFAM2, N₂O, lixiviationGIEC ch. 10 et 11 absents2Position et ITKF_MG, F_I, contrastes S06Amont et bilanbilan_ges.csv, bilan N7Runs et sensibilités8 scénarios, 2 jeuxBordure en tirets, phase à coder. Flèche pleine, dépendance d’ordre. Flèche en tirets, contrainte entre deux phases. Point orange, blocage au 6 novembre 2026.
Ordre des phases 0 à 7. Une flèche relie une phase à celle qui en dépend ; un repère orange signale un blocage présent au 6 novembre 2026.

État de départ au 6 novembre 2026

Sous-moduleCode dans Module_4CarboneParamètresTests
Stocks du sol et aériensCodé, inexécutable depuis la suppression des tables SeqCOI le 3 novembreTables GIEC saisies le 4 novembre, non lues49 contrôles natifs à réécrire
Itinéraire techniqueAbsentITK nominal décidé (facteurs à 1)Aucun
Élevage et effluentsAbsent, présent dans CarboneGeneriqueFacteurs génériques du prototype47 contrôles du prototype
TraitementAbsent, présent dans CarboneGenerique ; le staging de l’appariement calcule déjà MB, MS, N et CH₄ par méthaniseurFacteurs génériquesContrôles du prototype
TransportAbsent, décrit plus hautAucunAucun
Fertilisation et solAbsent, présent dans CarboneGenerique et CarboneRealiste (ALFAM2)GIEC 2019 tableaux 11.1 et 11.3, ALFAM2 4.2.1437 contrôles de CarboneRealiste
AmontAbsent, présent dans CarboneGeneriqueFacteurs génériquesAucun propre
Bilan GES et azoteAbsentPRG du prototypeBilans C et N fermés dans le prototype

L’appariement multi-passes est porté dans Module_3Allocation (étape I8, 175 fichiers conformes le 4 novembre). Les huit scénarios ont été relancés le 4 novembre dans output/runs/carbone_relance_20261104/.

Phases, tests et critères de fin

PhaseImplémentationParamétrageTestsCritère de fin
0 Socle communStructure Emission (poste, gaz, kg, support, date, niveau, source). Lecteur de table qui refuse une ligne sans source, une valeur négative ou une clé en double. Écrivain unique de emissions_postes.csv. Contrôle de bilan N et C à chaque étape (entrée = pertes + sortie).prg.csv, avec le choix 267 ou 273 pour le N₂OCas synthétiques à une et trois étapes ; table invalide rejetée ; bilan non fermé rejetéContrôles réussis sur données synthétiques, sans sous-module branché
1 Stocks rebranchés
I2, I3, P3, P4, P7, P8
Lecture de correspondance_carbone.csv (code, zone, compartiment), qui remplace usage_reference et family. Calcul SOCref × F_LU × F_MG × F_I depuis les tables GIEC. Résolveur des niveaux 3 à 0. Décomposition des 27 systèmes avec haies et des 8 systèmes mixtes. Variation nulle pour 100 et 200.correspondance_carbone.csv rempli depuis univers_codes.csv ; haie, vergers, sylvopastoral et riz selon les arbitrages du 4 novembreCas analytique du profil 15_5 (ΔC = −4,28 t C) ; surfaces conservées par LandUnit ; composantes à 100 % par code ; aucun niveau 0 ; 49 contrôles réécritsHuit scénarios avec CARBONE_LUC=true et niveau 0 vide
2 Position dans le run et ITK
I4, I5, I7
Appel après la production et après la passe 5. Structure Carbone_itk qui déduit F_MG et F_I de l’ITK de chaque fraction (table de scénario sans appariement, apports reçus avec). Décomposition usage, ITK, interaction. Contrastes Si − S0. Comptage des parcelles à usage initial unique.itk_scenario.csv par code et scénario ; maraîchage à 1, 0,92 en sensibilitéITK nominal identique à la phase 1 ; décomposition fermée sur un cas connu ; S0 − S0 nulMêmes stocks avec et sans appariement à ITK nominal
3 Élevage, effluents, traitementFormules de CarboneGenerique portées contenant par contenant, à partir des activités de l’appariement (inventaire des effluents, tonnes admises, staging). Chaque perte réduit le flux transmis. CH₄ entérique rattaché ici.facteurs_effluents.csv (B0, MCF, EF3, fractions volatilisées) ; procedes_traitement.csvBilans N et C fermés par contenant et sur la chaîne ; égalité avec CarboneGenerique ; quantité nulle, matière sans débouché, digestat sans séparationÉmissions des passes 2 et 3 d’un run réel, bilan d’azote fermé
4 TransportCalcul par paire décrit plus haut ; part amont du gazole transmise à l’amonttransport_vehicules.csv, transport_matieres.csv ; matière de la passe 5 à renseignerPaire analytique ; distance nulle ; somme des paires égale à la synthèseÉmission par paire sur les cinq passes, sans matière inconnue
5 Fertilisation et solNH₃ des apports (ALFAM2 pour les lisiers, facteur pour le reste) et des engrais ; N₂O direct, NOₓ, lixiviation, N₂O indirect (EF4, EF5), CO₂ de l’urée ; résidus, brûlage (éq. 2.27), minéralisation (éq. 11.8). Régime humide ou sec par profil. Règle contre le double compte avec F_I.facteurs_sol_azote.csv, alfam2_parametres.csv, part de canne brûlée par zoneALFAM2 contre le paquet R sur 14 356 points ; cas analytiques avec 44/28 ; bilan d’azote par fractionRun S0 avec appariement donnant toutes les pertes au champ par fraction
6 Amont et bilanAmont par intrant (engrais minéral, gazole, fret maritime, électricité, aliments importés). bilan_ges.csv et bilan d’azote du territoire. CO₂ biogénique à part.amont_intrants.csv, facteur du réseau réunionnais, production évitée par l’électricitéSomme des postes égale au total ; aucun double compte ; bilan d’azote du territoire ferméBilan complet de S0 avec la part de chaque niveau
7 Runs et sensibilités
P9
Huit scénarios, tout au niveau 1 puis meilleure référence, avec et sans appariementSensibilités sur D, biomasse et part des haies, sylvopastoral à 10 % d’arbres, maraîchage à 0,92, rizière inondée, PRG du N₂O, retour à videAnalyses dans 1B2_Analyses_Modèles/Une note d’analyse par série de runs

Les runs suivent la règle de charge du dossier 1B0_Modèles, avec un seul orchestrateur et deux cas au plus en même temps.

Arbitrages et risques ouverts

Arbitrages attendus

  • Biomasse de la haie, 9,4 t C/km (tableau 5.1) ou plantation de 100 t MS/ha
  • Part des parcelles des 27 codes qui porte une haie, 100 % aujourd’hui
  • Catégorie de l’ananas et des friches 220
  • Temps de croissance des ligneux dans le ΔC aérien
  • PRG du N₂O, 267 ou 273
  • Nature des 3,76 Mt de la passe 5 vers « Bruleur »
  • Calcul par paires de conversion, selon le comptage de I4

Risques

  • Mémoire de la JVM, 6 à 10 Go pour le Compromisator et 10 Go pour RasterDjikstra
  • Chapitres 10 et 11 du volume 4 du GIEC 2019 absents du vault, requis par les phases 3 et 5
  • Double compte entre facteurs de gestion du sol et émissions des apports
  • Matière inconnue dans la passe 5, qui bloque le transport et l’amont des imports

Fichiers de sortie

FichierContenuStatut
landunits.csvPar LandUnit, statut, ΔC sol et CO₂codé
surfaces_usages.csvPar LandUnit et usage, surfaces maintenue, retirée et ajoutéecodé
trajectoire.csvPar année, ΔC sol, CO₂ annuel et cumulé, couverturecodé
biomasse_aerienne_*.csvΔC aérien et total sol plus aérien sur les LandUnits renseignées dans les deux compartimentscodé
emissions_postes.csvGaz en kg par poste, support, date et niveau de méthodeprévu
bilan_ges.csvCO₂e par an (N₂O 267, valeur du prototype, 273 dans l’AR6 de 2021 ; CH₄ biogénique 27) et bilan d’azote par scénarioprévu

Une référence absente rend inconnu le support changé de la LandUnit, sans valeur zéro.