Skip to content
CPGE Futur

Exercice 2 : Impacts des perturbations anthropiques sur les écosystèmes

🌐 Énoncé du problème / Cadre d’étude

Section titled “🌐 Énoncé du problème / Cadre d’étude”

Les écosystèmes terrestres et aquatiques sont soumis à des pressions anthropiques majeures, parmi lesquelles l’eutrophisation des milieux aquatiques et la perturbation du cycle global du carbone occupent une place centrale. L’objectif de ce problème est d’analyser, à l’aide de modèles mathématiques compartimentés, la dynamique de la matière organique et de l’oxygène dissous dans un lac soumis à des apports accrus en nitrates (NO3−\text{NO}_3^-) et phosphates (PO43−\text{PO}_4^{3-}), ainsi que l’impact de ces perturbations sur le puits de carbone planétaire.

On considère un lac de volume V=1,0×107 m3V = 1{,}0 \times 10^7 \, \text{m}^3, considéré comme parfaitement agité. Ce lac est alimenté par une rivière à un débit constant Q=5,0×103 m3.j−1Q = 5{,}0 \times 10^3 \, \text{m}^3.\text{j}^{-1} et se dévide au même rythme, de sorte que le volume VV reste constant.

Soit C(t)C(t) la concentration en phosphore total dans le lac (exprimée en mg.m−3\text{mg}.\text{m}^{-3}). Le phosphore est le facteur limitant de la production primaire phytoplanctonique. Le bilan massique du phosphore dans le lac s’écrit selon l’équation différentielle :

dVdt=0  ⟹  VdCdt=QCin−QC(t)−ksVC(t)\frac{dV}{dt} = 0 \implies V \frac{dC}{dt} = Q C_{\text{in}} - Q C(t) - k_s V C(t)

Où :

  • Cin=50 mg.m−3C_{\text{in}} = 50 \, \text{mg}.\text{m}^{-3} est la concentration en phosphore dans les eaux entrantes (enrichies par les activités agricoles).
  • ks=0,05 j−1k_s = 0{,}05 \, \text{j}^{-1} est le coefficient global de sédimentation du phosphore (précipitation et enfouissement dans le sédiment).

Dans une seconde partie, nous étudierons la répercussion de la biomasse algale sur la concentration en oxygène dissous O2\text{O}_2 et l’apparition potentielle d’une zone d’hypoxie.

❓ Questions progressives (incluant une analyse critique du modèle)

Section titled “❓ Questions progressives (incluant une analyse critique du modèle)”

  1. Analyse préliminaire et mise en équation

    • Rappeler brièvement le rôle biogéochimique du phosphore dans l’eutrophisation des écosystèmes aquatiques.
    • Mettre l’équation différentielle régissant l’évolution de C(t)C(t) sous la forme canonique : dCdt+αC(t)=β\frac{dC}{dt} + \alpha C(t) = \beta Exprimer α\alpha et β\beta en fonction de QQ, VV et ksk_s. Préciser leurs dimensions et unités respectives.

  2. Résolution et régime permanent

    • Résoudre cette équation différentielle en supposant qu’initialement, à t=0t = 0, la concentration dans le lac est C(0)=C0C(0) = C_0.
    • Déterminer l’expression de la concentration à l’équilibre stationnaire CeqC_{\text{eq}} lorsque t→+∞t \to +\infty. Calculer sa valeur numérique en mg.m−3\text{mg}.\text{m}^{-3}.
    • Calculer le temps caractéristique τ\tau du système (exprimé en jours) et interpréter sa signification écologique.

  3. Couplage biogéochimique : Impact sur l’oxygène dissous

    • La production primaire nette génère une biomasse végétale dont la décomposition bactérienne consomme de l’oxygène selon une stœchiométrie globale où 1 mg1 \, \text{mg} de phosphore conduit, via la stœchiométrie de Redfield, à la consommation d’environ 150 mg150 \, \text{mg} de O2\text{O}_2 lors de la minéralisation.
    • Soit O(t)O(t) la concentration en oxygène dissous. On pose l’équation simplifiée de son bilan : dOdt=KL(Osat−O(t))−Rbio\frac{dO}{dt} = K_L (O_{\text{sat}} - O(t)) - R_{\text{bio}} où KL=0,20 j−1K_L = 0{,}20 \, \text{j}^{-1} est le coefficient de réoxygénation par l’atmosphère, Osat=10,0 mg.m−3O_{\text{sat}} = 10{,}0 \, \text{mg}.\text{m}^{-3} la concentration à saturation, et Rbio=γ⋅ks⋅CeqR_{\text{bio}} = \gamma \cdot k_s \cdot C_{\text{eq}} le taux de consommation biologique lié à la dégradation de la matière organique sédimentée/suspendue, avec γ=150\gamma = 150.
    • Calculer RbioR_{\text{bio}} à l’équilibre, puis déterminer la concentration stationnaire en oxygène dissous OeqO_{\text{eq}}.

  4. Analyse critique du modèle et esprit concours

    • Le seuil critique d’hypoxie létale pour la faune ichtyologique est fixé à 2,0 mg.m−32{,}0 \, \text{mg}.\text{m}^{-3}. Le lac étudié est-il en crise dystrophe ?
    • Discuter les limites majeures de ce modèle compartimenté simple face à la réalité d’un écosystème lacustre complexe (on attend au moins deux arguments critiques rigoureux portant sur les hypothèses biologiques ou physiques).

💡 Pistes & Indications méthodologiques

Section titled “💡 Pistes & Indications méthodologiques”
  • Q1 : Pensez à diviser l’équation initiale par le volume VV pour obtenir la forme standard d’une équation différentielle linéaire du premier ordre à coefficients constants. Vérifiez l’homogénéité des termes (flux massique par uniteˊ de volume\text{flux massique par unité de volume}).
  • Q2 : La solution générale d’une équation du type y′+αy=βy' + \alpha y = \beta est la somme de la solution générale de l’équation sans second membre et d’une solution particulière constante.
  • Q3 : Procédez par substitution directe des valeurs trouvées à la question précédente. Attention à bien mener les applications numériques avec leurs unités.
  • Q4 : Pour l’esprit critique, pensez aux variations saisonnières (température et stratification thermique de type thermocline/hypolimnion), ainsi qu’à la non-linéarité des interactions biologiques (boucles de rétroaction, prédation, recyclage interne des nutriments par les sédiments).

📋 Corrigé détaillé pas à pas

Section titled “📋 Corrigé détaillé pas à pas”

1. Analyse préliminaire et mise en équation

Section titled “1. Analyse préliminaire et mise en équation”
  • Rôle du phosphore : Le phosphore sous forme de phosphate (PO43−\text{PO}_4^{3-}) est généralement le facteur limitant de la production primaire dans les écosystèmes dulçaquicoles. Un apport massif anthropique lève cette limitation, entraînant une prolifération algale (bloom phytoplanctonique).
  • Mise en équation : L’équation initiale est : VdCdt=QCin−QC(t)−ksVC(t)V \frac{dC}{dt} = Q C_{\text{in}} - Q C(t) - k_s V C(t) En divisant par le volume VV : dCdt=QVCin−(QV+ks)C(t)\frac{dC}{dt} = \frac{Q}{V} C_{\text{in}} - \left(\frac{Q}{V} + k_s\right) C(t) On identifie la forme canonique dCdt+αC(t)=β\frac{dC}{dt} + \alpha C(t) = \beta avec :
    • α=QV+ks\alpha = \frac{Q}{V} + k_s (dimension de temps−1\text{temps}^{-1}, unité : j−1\text{j}^{-1})
    • β=QVCin\beta = \frac{Q}{V} C_{\text{in}} (dimension d’une concentration par unité de temps, unité : mg.m−3.j−1\text{mg}.\text{m}^{-3}.\text{j}^{-1})

2. Résolution et régime permanent

Section titled “2. Résolution et régime permanent”

  • Résolution de l’équation différentielle : L’équation homogène associée dCdt+αC(t)=0\frac{dC}{dt} + \alpha C(t) = 0 a pour solution générale : CH(t)=Ke−αt(K∈R)C_H(t) = K e^{-\alpha t} \quad (K \in \mathbb{R}) Une solution particulière constante CpC_p s’obtient en annulant la dérivée : 0+αCp=β  ⟹  Cp=βα0 + \alpha C_p = \beta \implies C_p = \frac{\beta}{\alpha} La solution générale est donc : C(t)=Ke−αt+βαC(t) = K e^{-\alpha t} + \frac{\beta}{\alpha} En utilisant la condition initiale C(0)=C0C(0) = C_0, on trouve K=C0−βαK = C_0 - \frac{\beta}{\alpha}. Ainsi : C(t)=(C0−βα)e−αt+βαC(t) = \left(C_0 - \frac{\beta}{\alpha}\right) e^{-\alpha t} + \frac{\beta}{\alpha}

  • Concentration à l’équilibre stationnaire (t→+∞t \to +\infty) : lim⁡t→+∞C(t)=Ceq=βα=QVCinQV+ks=Cin1+ksVQ\lim_{t \to +\infty} C(t) = C_{\text{eq}} = \frac{\beta}{\alpha} = \frac{\frac{Q}{V} C_{\text{in}}}{\frac{Q}{V} + k_s} = \frac{C_{\text{in}}}{1 + \frac{k_s V}{Q}}

    • Application numérique : Calculons d’abord le taux de renouvèlement QV\frac{Q}{V} : QV=5,0×103 m3.j−11,0×107 m3=5,0×10−4 j−1\frac{Q}{V} = \frac{5{,}0 \times 10^3 \, \text{m}^3.\text{j}^{-1}}{1{,}0 \times 10^7 \, \text{m}^3} = 5{,}0 \times 10^{-4} \, \text{j}^{-1} On a ks=5,0×10−2 j−1k_s = 5{,}0 \times 10^{-2} \, \text{j}^{-1}. Ainsi, α=5,0×10−4+5,0×10−2=5,05×10−2 j−1\alpha = 5{,}0 \times 10^{-4} + 5{,}0 \times 10^{-2} = 5{,}05 \times 10^{-2} \, \text{j}^{-1}. La concentration à l’équilibre vaut : Ceq=5,0×10−4×505,05×10−2=0,0250,0505≈0,495 mg.m−3C_{\text{eq}} = \frac{5{,}0 \times 10^{-4} \times 50}{5{,}05 \times 10^{-2}} = \frac{0{,}025}{0{,}0505} \approx 0{,}495 \, \text{mg}.\text{m}^{-3} (Note : Bien que faible en apparence dans la masse d’eau en raison de la sédimentation rapide, ce flux soutient une biomasse algale considérable).

  • Temps caractéristique : Le temps caractéristique du système est τ=1α=15,05×10−2≈19,8 jours\tau = \frac{1}{\alpha} = \frac{1}{5{,}05 \times 10^{-2}} \approx 19{,}8 \, \text{jours}. Ce temps indique la rapidité avec laquelle le lac s’adapte à une modification des apports amont (inertie du système d’environ 20 jours).

3. Couplage biogéochimique : Impact sur l’oxygène dissous

Section titled “3. Couplage biogéochimique : Impact sur l’oxygène dissous”

  • Calcul du taux de consommation biologique RbioR_{\text{bio}} : Le taux de minéralisation de la matière organique liée au phosphore à l’équilibre est proportionnel au flux de phosphore sédimenté/consommé par unité de temps : Rbio=γ⋅(ksCeq)ou assimileˊ au flux globalR_{\text{bio}} = \gamma \cdot (k_s C_{\text{eq}}) \quad \text{ou assimilé au flux global} Plus rigoureusement, la consommation d’oxygène est liée à la dégradation de la biomasse produite, que l’on modélise par Rbio=γ⋅ks⋅CeqR_{\text{bio}} = \gamma \cdot k_s \cdot C_{\text{eq}}.

    • Application numérique : Rbio=150×(5,0×10−2 j−1)×0,495 mg.m−3=150×0,02475=3,7125 mg.m−3.j−1R_{\text{bio}} = 150 \times (5{,}0 \times 10^{-2} \, \text{j}^{-1}) \times 0{,}495 \, \text{mg}.\text{m}^{-3} = 150 \times 0{,}02475 = 3{,}7125 \, \text{mg}.\text{m}^{-3}.\text{j}^{-1}

  • Détermination de la concentration stationnaire en oxygène OeqO_{\text{eq}} : À l’équilibre pour l’oxygène, dOdt=0\frac{dO}{dt} = 0, donc : KL(Osat−Oeq)−Rbio=0  ⟹  Oeq=Osat−RbioKLK_L (O_{\text{sat}} - O_{\text{eq}}) - R_{\text{bio}} = 0 \implies O_{\text{eq}} = O_{\text{sat}} - \frac{R_{\text{bio}}}{K_L}

    • Application numérique : Oeq=10,0−3,71250,20=10,0−18,56=−8,56 mg.m−3O_{\text{eq}} = 10{,}0 - \frac{3{,}7125}{0{,}20} = 10{,}0 - 18{,}56 = -8{,}56 \, \text{mg}.\text{m}^{-3}
    • Résultat mathématique brut aberrant : La valeur obtenue est négative, ce qui traduit une limite physique évidente du modèle simple (la concentration réelle ne peut descendre en dessous de 00). Cela signifie que la demande biologique en oxygène dépasse largement la capacité maximale de réoxygénation par l’atmosphère dans les conditions du modèle, conduisant à une anoxie totale (Oeq=0 mg.m−3O_{\text{eq}} = 0 \, \text{mg}.\text{m}^{-3}).

4. Analyse critique du modèle et esprit concours

Section titled “4. Analyse critique du modèle et esprit concours”

  • Risque de crise dystrophe : Puisque la concentration en oxygène chute à 0 mg.m−30 \, \text{mg}.\text{m}^{-3} (bien en deçà du seuil critique d’hypoxie de 2,0 mg.m−32{,}0 \, \text{mg}.\text{m}^{-3}), le lac est en situation de crise dystrophe majeure (mortalité piscicole massive, apparition de conditions anaérobies favorables aux bactéries sulfato-réductrices produisant du H2S\text{H}_2\text{S}).

  • Limites critiques du modèle compartimenté :

    1. Hypothèse de mélange parfait (agitation homogène) : En réalité, les lacs profonds subissent une stratification thermique saisonnière (épilimnion chaud et oxygéné, hypolimnion froid et isolé). L’oxygène ne se diffuse pas uniformément, et l’hypolimnion s’anoxie beaucoup plus vite que ne le prédit un modèle à un seul compartiment homogène.
    2. Dynamique non linéaire et rétroactions sédimentaires : Le modèle suppose un coefficient de sédimentation ksk_s constant. Or, en condition d’anoxie, le fer présent dans les sédiments change de degré d’oxydation (Fe3+→Fe2+\text{Fe}^{3+} \to \text{Fe}^{2+}), ce qui libère massivement le phosphore piégé (libération benthique ou internal loading). Ce recyclage interne amplifie l’eutrophisation de manière autocatalytique, un phénomène totalement ignoré par l’équation linéaire stationnaire proposée ici.