OYSTER : Filtrage benthique

  • Ce module complémetnaire permet de prendre en compte la présence d’huitres sur le fond qui filtrent le matériel phytoplanctonique.
  • Dans cette option, les huitres ne sont pas simulées en tant que variable d’état, mais leur densité par maille agit sur le terme de “filtrage benthique” (Puits d’azote au fond)
  • Ce “filtrage benthique” correspond à la filtration des huitres qui dépend des MES et de la température.
  • Les paramètres qui sont utilisés dans ce module sont donnnés dans la namelist namoysterfiltbent dans parabiolo.txt (Modele Huitre : parametres dans les namelist)

Note

Ce module est activé en utilisant la clé key_oyster_benthos (+ key_biolo )

La distribution des huitres

  • La distribution spatiale des huitres est determinée par le nombre d’huitres nbhuitre(i,j)
  • nbhuitre doit être initialisé, estimé ou lu dans un fichier.
  • Cette définition est réalisée par l’utilisateur dans bioloinit.F90
  • la hauteur des tables y peuvent être également définies (hautable) (en m)

La filtration

  • L’intensité de la filtration dépend de la température et de la quantité de matière en suspension (MES) dans l’eau.

  • Un effet de colmatage est pris en compte :

    \(Colmat_{huitre} = \min (0.,X_{colmat_thr} - MES)\)
    \(Effettemp_{huitre} = X_{Filt-temp1} ∗ [TEMPERATURE−X_{Tempopt-Filt}]^2\)
  • Si MES < p_huitre_thr_mes

    \(Filt= X_{filtmax} - Effettemp_{huitre}\)
  • sinon

    \(Filt= X_{filt1}*MES + X_{filt2} - Effettemp_{huitre}\)
    \(Filt= Filt * (X_{DW}^{X_{expfilt}}\)
    \(Filt= Filt * \exp ( X_{colmat} * Colmat_{huitre} )\)
    Si Filt < 0 ==> Filt=0.
    \(T_{filtbenth} = Filt * nbhuitre(i,j) /surf(i,j)\)
  • avec :

    • \(X_{colmat_thr}\) = p_huitre_thr_colmat
    • \(X_{Filt-temp1}\) = p_huitre_temp1
    • \(X_{Tempopt-Filt}\) = p_huitre_temp_opt
    • \(X_{filtmax}\) = p_huitre_filt_std
    • \(X_{filt1}\) = p_huitre_loi_mes1
    • \(X_{filt2}\) = p_huitre_loi_mes2
    • \(X_{DW}\) = p_huitre_DW
    • \(X_{expfilt}\) = p_huitre_exp_allom
    • \(X_{colmat}\) = p_huitre_loi_colmat