Modèle conceptuel du simulateur

Dans cette partie, les différents modules de Syst'N sont présentés de manière détaillée

Module Absorption et fixation

Formalisme général

Les formalismes de ce module sont largement repris du modèle Azodyn (Jeuffroy et al., 1999).
L’absorption globale des cultures est estimée par comparaison entre les besoins azotés de la plante et un potentiel maximal d’absorption spécifique. Dans le cas des cultures fixatrices (i.e. légumineuses), une partie de l’assimilation totale d’azote est issue de la fixation symbiotique.
Les étapes de calcul de ce module sont les suivantes :

Calcul des besoins en azote de la plante (NBESOIN)

Le calcul des besoins en azote du couvert se fait à partir de la courbe de dilution maximale de l’azote dans la plante :
NBESOIN(j) = COEF_RAC_N x [(MSaer(j) x %Nmax(j)) – (MSaer(j-1) x %Nmax(j-1))]
Avec :
COEF_RAC_N : coefficient de prise en compte des racines dans les besoins en azote de la plante
MSaer(j) : quantité de matière sèche aérienne au jour j (T.ha-1)
%Nmax : teneur en azote maximale dans les parties aériennes (voir le calcul ici)

Courbe de dilution maximale de l’azote dans les parties aériennes

La teneur en azote maximale se calcule comme suit :
Si MSaer(j) < MSaer_seuil alors %Nmax = amax x MSaer_seuil-bmax
Si MSaer(j)MSaer_seuil alors %Nmax = amax x MSaer(j)-bmax
Avec
MSaer(j) : quantité de matière sèche aérienne (T.ha-1)
MSaer_seuil : quantité de matière sèche en-deçà de laquelle la courbe maximale décrit un plateau (quantités de matière sèche faibles)(T.ha-1)
%Nmax : teneur en azote maximale dans les parties aériennes (%)
amax, bmax : paramètres de la courbe de dilution maximale de l’azote dans les parties aériennes, connus par espèce

panneau information Le courbe maximale a été validée jusqu’à floraison, en général, pour les différentes espèces. Après ce stade, plusieurs travaux montrent, soit que la courbe initiale peut être maintenue, soit que la teneur en azote maximale des parties aériennes de la culture reste stable. Ici, on propose de garder la même courbe maximale

Calcul de l’absorption potentielle maximale (NABS_POT)

L’absorption potentielle maximale d’azote par la plante au jour j est calculée avec l’équation suivante :
NABS_POT (j) = Vabs_max x ∆SOMTEMP(j)
Avec
Vabs_max : vitesse maximale d’absorption d’azote propre à chaque espèce (kgN.d°j-1 )
∆SOMTEMP(j) : somme dde degrés-jours du jour (degrés-jours) (voir dans Interventions)

panneau information Quand la température moyenne de l’air dépasse la température maximale de croissance de la culture alors ∆SOMTEMP(j) = 0 et l’absorption potentielle maximale d’azote est elle aussi nulle. D’un point de vue physiologique cela traduit le fait qu’au-delà de cette température maximale de croissance, la plante ferme ses stomates, pour limiter son évapotranspiration, arrêtant ainsi la circulation d’eau dans la plante (pas d’entrée ni de sortie) et donc l’absorption d’azote.

Calcul de la demande en azote de la plante (NA_D)

La demande en azote de la plante est définie comme :
NA_D(j) = min(NBESOIN(j); NABS_POT (j)) x COEFF_DENSITE
NBESOIN(j) : besoin en azote du couvert (kgN.ha-1)
NABS_POT (j) : quantité d’azote maximale potentiellement absorbable par le couvert (kgN.ha-1)
COEFF_DENSITE : coefficient de prise en compte de la densité de peuplement UNIQUEMENT pour les cultures intermédiaires et les repousses. Il est égal à 1 si le peuplement est défini comme « dense » par l’utilisateur dans l’IHM d’entrée, et à 0.5 s’il est défini comme « clairsemé ». Si pas de renseignement sur la densité alors il est par défaut égal à 1. De même, il est systématiquement égal à 1 pour les cultures principales

Calcul de la fixation symbiotique (NFIX)

La quantité d’azote fixé est calculée en fonction de la demande en azote de la plante, de la quantité d’azote minéral présente dans le sol, du stade végétatif de la plante, et du stress hydrique :
NFIX(j) = min(NA_D(j) ; NfixN_MIN ; NfixSTADE ; NfixEAU)
Avec
NFIX(j): quantité d’azote fixé (kgN.ha-1)
NA_D: besoin en azote de la culture (kgN.ha-1)
NfixN_MIN : fonction liée à la présence d’azote minéral dans le sol (kgN.ha-1)
NfixSTADE : quantité d’azote que peut fixer la plante en fonction de son stade végétatif(kgN.ha-1)
NfixEAU : fonction liée au stress hydrique (kgN.ha-1)

Fonction liée à la présence d’azote minéral dans le sol

La quantité d’azote fixable en fonction de la quantité d’azote minéral présente dans le sol est calculée comme suit :
Avant début floraison :

NfixN_MIN(j) = NA_D(j) x -2.96 x Nmin(j) + 142,77
100
Après début floraison :
NfixN_MIN(j) = NA_D(j) x -3,27 x Nmin(j) + 110.1901
100
Avec
NA_D(j) (j) : besoin en azote de la culture (kgN.ha-1)
Nmin(j) : quantité d’azote minéral dans la couche de sol colonisée par les racines (kgN.ha-1)

Fonction liée au stade végétatif de la culture

La quantité d’azote fixable selon le stade de la culture est calculée comme suit :
Si SOMTEMP(j) < SOMTEMP_DEBUTFIX alors NfixSTADE(j) = 0
Si SOMTEMP_DEBUTFIXSOMTEMP(j) < SOMTEMP_DRG alors

NfixSTADE(j) = 0,0000357 x [ 0,0119 + SOMTEMP_DRG - SOMTEMP(j)] x MSaer(j)
0,0000357
Si SOMTEMP_DRGSOMTEMP(j)SOMTEMP_FSLA alors NfixSTADE(j) = 0,0119 x MSaer(j)
Si SOMTEMP(j) > SOMTEMP_FSLA alors NfixSTADE(j) = 0
Avec
SOMTEMP(j) : somme de températures depuis le semis (degrés-jour)
SOMTEMP_DEBUTFIX : somme de températures entre le semis et le début de fixation (degrés-jour)
SOMTEMP_DRG : somme de températures entre le semis et le stade DRG (degrés-jour)
SOMTEMP_FSLA : somme de températures entre le semis et la fin de stade limite d’avortement (degrés-jour)
MSaer(j) = quantité de matière sèche aérienne au jour j (kgMS/ha) (voir module Croissance aérienne)

Fonction liée au stress hydrique

La quantité d’azote fixable en fonction de l’humidité du sol (i.e. d’un potentiel stress hydrique) est calculée comme suit :

IndStressH = HHpF
HCCHpF
Après floraison, si IndStressH < 1/3 pendant 7 jours consécutifs alors NfixEAU = 0 (i.e. la fixation symbiotique est alors inhibée de façon irréversible)
Avec
H : humidité volumique de la couche de sol colonisée par les racines (% volumique)
HCC : humidité volumique à la capacité au champ de la couche de sol colonisée par les racines (% volumique)
HpF : humidité volumique au point de flétrissement de la couche de sol colonisée par les racines (% volumique)

Calcul de l’azote à absorber dans le sol par la culture (NABS_SOL)

La quantité d’azote (restant) à absorber dans le sol par la culture est calculée comme suit :
NABS_SOL = NA_D - NFIX
Pour chaque couche élémentaire de sol colonisée par les racines (1 cm par 1 cm), on calcule la quantité d’azote absorbée par la plante dans cette couche en fonction du besoin d’absorption de la plante et de la quantité d’azote minéral dans cette couche par rapport à la quantité d’azote minéral totale sur la profondeur racinaire :

NABS [i] = NABS_SOL x NMIN [i]
Nmin_tot_rac<
NABS [i] : quantité d’azote minéral absorbée par la culture dans la couche de sol i (kgN.ha-1)
NMIN [i] : quantité d’azote minéral disponible dans la couche de sol i (kgN.ha-1)
Nmin_tot_rac : quantité d’azote minéral totale disponible dans la couche de sol colonisée par les racines (kgN.ha-1)

Enfin, la quantité totale d’azote réellement absorbé (NABS_SOL) dans le sol est égale à :
NABS = ∑NABS[i]

Calcul de l’azote total contenu dans la plante (NPLANTE)

La quantité d’azote totale contenue dans la plante le jour j est calculée comme suit :
NPLANTE (j) = NPLANTE (j-1) + NFIX (j) + NABS (j)
Avec
NPLANTE (j-1) : quantité d’azote totale contenue dans la plante au jour précédent (kgN.ha-1)
NFIX (j) : quantité d’azote fixé par la plante le jour j (kgN.ha-1)
NABS (j) : quantité d’azote absorbé dans le sol par la plante le jour j (kgN.ha-1)

panneau attention Dans Syst’N l’absorption d’azote est complétement déconnectée du statut hydrique du sol. Un stress hydrique n’impacte donc pas l’absorption d’azote : ne correspond pas à la réalité agronomique