pith. sign in

arxiv: 2604.18818 · v1 · submitted 2026-04-20 · 🧮 math.DS

Qualitative Analysis of a Three-Stage Anaerobic Digestion Model with Microbial Decay

Radhouane Fekih-Salem This is my paper

Pith reviewed 2026-05-10 03:31 UTC · model grok-4.3

classification 🧮 math.DS
keywords anaerobic digestionmicrobial mortalitythree-stage modelhydrolysis mechanismsequilibria stabilityAM2 modelsubstrate inhibitiondynamical systems
0
0 comments X

The pith

Incorporating microbial mortality into the hydrolysis stage of a three-stage anaerobic digestion model creates a more general framework with enriched qualitative dynamics compared to the classical AM2 model.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends the classical AM2 model by adding microbial mortality terms specifically to the hydrolysis of particulate organic matter, while keeping distinct dilution and decay rates across the three stages of hydrolysis, acidogenesis, and methanogenesis. Two hydrolysis formulations are studied: a simple first-order rate and a biomass-dependent rate that separates hydraulic washout from intrinsic cell death. For each version the authors prove that solutions exist, stay positive, and remain bounded, then locate every equilibrium and determine its local stability as functions of the operating parameters. The mortality terms produce different persistence thresholds between the two hydrolysis cases and allow the model to reflect biological cell death, which the authors argue makes the dynamics richer and more realistic than in the original AM2 system.

Core claim

We extend the AM2 model by incorporating microbial mortality into the hydrolysis stage, considering both first-order and biomass-dependent hydrolysis mechanisms along with non-monotonic growth functions for substrate inhibition. We establish the well-posedness of the model, prove positivity and boundedness of solutions, characterize all equilibria, and analyze their local stability with respect to key operating parameters. The resulting framework is more general and biologically relevant than the classical model, thereby enriching the qualitative dynamics.

What carries the argument

The three-stage ODE system with microbial mortality added only to the hydrolysis equations, using distinct dilution and decay rates per stage, and non-monotonic growth rates to capture substrate inhibition.

If this is right

  • The biomass-dependent hydrolysis mechanism produces persistence conditions that differ from those of the first-order mechanism.
  • Every equilibrium of the system is explicitly located and its local stability is determined by the dilution rates and decay rates.
  • Solutions remain positive and bounded for all time, preserving biological interpretability.
  • The extended model admits a broader set of long-term behaviors than the classical AM2 model because mortality is treated separately from washout.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The stability conditions could be used to set dilution rates that avoid washout in industrial biogas reactors.
  • Adding mortality to the other stages or to multi-species versions of the model might reveal further changes in community structure.
  • Global stability or bifurcation analysis would be a natural next step to map how parameters move the system between different stable states.
  • Direct comparison with time-series data from real digesters would test whether the predicted stability boundaries appear when decay rates are measured independently.

Load-bearing premise

Microbial mortality can be modeled by adding distinct decay rates specifically to the hydrolysis stage and that non-monotonic growth functions adequately represent substrate inhibition effects.

What would settle it

An experiment that varies microbial decay rates in a lab-scale digester and finds no corresponding shift in the observed persistence thresholds or stability of steady states would challenge the claim that the mortality terms enrich the dynamics in the predicted way.

read the original abstract

We extend a three-stage anaerobic digestion model by incorporating microbial mortality into the hydrolysis of particulate organic matter. The model describes hydrolysis, acidogenesis, and methanogenesis, each with distinct dilution and decay rates, and accounts for non-monotonic growth in order to capture substrate inhibition. Two hydrolysis mechanisms are considered: a first-order formulation and a biomass-dependent one. The latter distinguishes hydraulic washout from intrinsic mortality and leads to different persistence conditions. For both models, we establish well-posedness and prove the positivity and boundedness of solutions. We then characterize all equilibria and analyze their local stability with respect to key operating parameters. The inclusion of microbial mortality provides a more general and biologically relevant framework, thereby enriching the qualitative dynamics compared to the classical AM2 model.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript extends the classical three-stage anaerobic digestion (AM2) model by incorporating microbial mortality specifically into the hydrolysis stage, with distinct dilution and decay rates per stage (hydrolysis, acidogenesis, methanogenesis) and non-monotonic growth functions to capture substrate inhibition. Two hydrolysis formulations are analyzed: standard first-order and biomass-dependent. For both, the authors establish well-posedness, positivity and boundedness of solutions, fully characterize equilibria, and perform local stability analysis with respect to key operating parameters such as dilution rates. The central claim is that including microbial mortality yields a more general, biologically relevant framework that enriches the qualitative dynamics relative to the classical AM2 model, including altered persistence conditions under the biomass-dependent hydrolysis.

Significance. If the local stability and equilibrium results hold and the enrichment is demonstrated, the work would provide a mathematically rigorous extension of standard chemostat-type models for anaerobic digestion, potentially useful for predicting washout thresholds and stability in bioprocess control. The stage-specific rates and non-monotonic kinetics address realistic biological features like differential mortality and inhibition, which could inform parameter estimation in applications. However, the significance hinges on whether the new terms produce dynamics (e.g., modified bifurcation sets or attractors) that cannot be recovered simply by setting decay rates to zero in the same equations.

major comments (2)
  1. [Abstract and §1] Abstract and §1 (Introduction): The assertion that the biomass-dependent hydrolysis 'leads to different persistence conditions' and that the model 'enriches the qualitative dynamics compared to the classical AM2 model' is not supported by any explicit comparison. No bifurcation diagrams, persistence thresholds, or global attractor descriptions are contrasted with the zero-mortality case (i.e., all decay rates set to zero while retaining the same non-monotonic growth functions). Local stability analysis alone does not establish enrichment; a direct side-by-side analysis of equilibria or stability regions is required to substantiate the central claim.
  2. [§2] §2 (Model formulation), equations for hydrolysis terms: The biomass-dependent hydrolysis is presented as distinguishing hydraulic washout from intrinsic mortality, but the manuscript does not verify that the resulting persistence conditions differ qualitatively from those obtained by uniform decay rates across stages or by monotonic growth. Without this, the 'more general framework' remains a parameter extension rather than a demonstrated source of new phenomena.
minor comments (2)
  1. [§2] Notation for growth functions (e.g., μ_i(S)) should be clarified to distinguish the non-monotonic form from standard Monod kinetics, perhaps with an explicit plot or parameter range where inhibition occurs.
  2. [Stability sections] In the stability analysis sections, the Jacobian matrices at equilibria are computed but the characteristic equations could be simplified or presented more explicitly to aid verification of the Routh-Hurwitz criteria.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and constructive suggestions. We address each major comment below and will revise the manuscript accordingly to strengthen the claims with explicit comparisons.

read point-by-point responses

  1. Referee: [Abstract and §1] Abstract and §1 (Introduction): The assertion that the biomass-dependent hydrolysis 'leads to different persistence conditions' and that the model 'enriches the qualitative dynamics compared to the classical AM2 model' is not supported by any explicit comparison. No bifurcation diagrams, persistence thresholds, or global attractor descriptions are contrasted with the zero-mortality case (i.e., all decay rates set to zero while retaining the same non-monotonic growth functions). Local stability analysis alone does not establish enrichment; a direct side-by-side analysis of equilibria or stability regions is required to substantiate the central claim.

    Authors: We agree that an explicit side-by-side comparison is necessary to fully substantiate the enrichment claim. The manuscript derives distinct persistence conditions for the biomass-dependent hydrolysis (Theorems 3.3 and 4.2) that arise specifically from separating hydraulic washout and intrinsic mortality, conditions which do not appear when decay rates are uniformly zero. However, we acknowledge the absence of direct contrasts such as stability region diagrams or threshold comparisons. In the revised version, we will add a dedicated remark in §1 and a short comparative analysis in §5, analytically contrasting the washout thresholds and local stability boundaries with the zero-decay case while retaining the non-monotonic growth functions. This will demonstrate that the new terms produce persistence conditions not recoverable by parameter tuning alone. revision: yes

  2. Referee: [§2] §2 (Model formulation), equations for hydrolysis terms: The biomass-dependent hydrolysis is presented as distinguishing hydraulic washout from intrinsic mortality, but the manuscript does not verify that the resulting persistence conditions differ qualitatively from those obtained by uniform decay rates across stages or by monotonic growth. Without this, the 'more general framework' remains a parameter extension rather than a demonstrated source of new phenomena.

    Authors: The referee correctly notes that we have not provided an explicit verification against uniform decay rates or monotonic growth functions. The stage-specific dilution and decay rates, combined with the biomass-dependent hydrolysis term, are introduced to reflect biological distinctions, and the resulting equilibria conditions (Proposition 2.1) depend on the biomass term in a manner that uniform decay does not replicate. To address this, we will insert a clarifying remark in §2 and include a brief analytical comparison in the results section showing how uniform decay alters the stability boundaries relative to the stage-specific case. We maintain that the framework is more general, but agree to supply the requested verification to demonstrate the new phenomena. revision: partial

Circularity Check

0 steps flagged

No circularity; standard extension and analysis of explicit ODE model

full rationale

The paper defines an extended three-stage model by adding microbial mortality terms (with stage-specific dilution and decay rates) and non-monotonic growth functions to the classical AM2 equations. It then applies textbook dynamical-systems arguments: well-posedness via standard existence theorems, positivity/boundedness via comparison or invariant-set methods, exhaustive equilibrium enumeration, and local stability via Jacobian eigenvalues or Routh-Hurwitz criteria. None of these steps reduce to a fitted parameter renamed as prediction, a self-definitional loop, or a load-bearing self-citation whose content is merely the present claim. The statement that the biomass-dependent hydrolysis “leads to different persistence conditions” is a direct algebraic consequence of the new decay term appearing in the biomass equation, not a tautology. The enrichment claim is therefore an assertion about the mathematical consequences of the added terms rather than a circular re-derivation of the inputs.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The model relies on several biological rate parameters treated as inputs rather than derived, along with standard mathematical assumptions for ODE systems.

free parameters (3)
  • dilution rates
    Distinct dilution rates for each of the three stages are model parameters.
  • decay rates
    Microbial decay rates added for each stage.
  • growth function parameters
    Parameters defining the non-monotonic growth rates.
axioms (2)
  • standard math The system is governed by a system of ordinary differential equations with unique solutions.
    Invoked for well-posedness proofs.
  • domain assumption Growth rates are non-monotonic to capture substrate inhibition.
    Biological modeling choice stated in the abstract.

pith-pipeline@v0.9.0 · 5422 in / 1445 out tokens · 75348 ms · 2026-05-10T03:31:22.492150+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

  1. [1]

    Benyahia, T

    B. Benyahia, T. Sari, B. Cherki and J. Harmand, Bifurcation and stability analysis of a two step model for monitoring anaerobic digestion processes, J. Proc. Control, 221008–1019 (2012)

    work page 2012

  2. [2]

    Bernard, Z

    O. Bernard, Z. Hadj-Sadok, D. Dochain, A. Genovesi and J.P. Steyer, Dynami- cal model development and parameter identification for an anaerobic wastewater treatment process,Biotechnol. Bioeng.,75, 424–438 (2001)

    work page 2001

  3. [3]

    Fekih-Salem, N

    R. Fekih-Salem, N. Abdellatif, T. Sari and J. Harmand, Analyse mathématique d’un modèle de digestion anaérobie à trois étapes, Arima Journal,17, 53–71 (2014)

    work page 2014

  4. [4]

    Fekih-Salem, Y

    R. Fekih-Salem, Y. Daoud, N. Abdellatif and T. Sari, A mathematical model of anaerobic digestion with syntrophic relationship, substrate inhibition and distinct removal rates, SIAM J. Appl. Dyn. Syst20, 1625–1654 (2021)

    work page 2021

  5. [5]

    Nouaoura, N

    S. Nouaoura, N. Abdellatif, R. Fekih-Salem and Tewfik Sari., Mathematical analysis of a three-tiered model of anaerobic digestion, Siam Journal on Applied Mathemat- ics,81, 1264–1286 (2021)

    work page 2021

  6. [6]

    Nouaoura, R

    S. Nouaoura, R. Fekih-Salem, N. Abdellatif and Tewfik Sari., Mathematical analysis of a three-tiered food-web in the chemostat, Discrete and Continuous Dynamical Systems - Series B,26, 5601-–5625 (2021)

    work page 2021

  7. [7]

    T.Sari, BestOperatingConditionsforBiogasProductioninSomeSimpleAnaerobic Digestion Models, Processes,2, 258 (2022)

    work page 2022

  8. [8]

    Sari and J

    T. Sari and J. Harmand, A model of a syntrophic relationship between two microbial species in a chemostat including maintenance, Math. Biosci.,275, 1–9 (2016)

    work page 2016

  9. [9]

    Sari and B

    T. Sari and B. Benyahia, The operating diagram for a two-step anaerobic digestion model, Nonlinear Dyn.,105, 2711–2737 (2021)

    work page 2021

  10. [10]

    Sari and M.J

    T. Sari and M.J. Wade, Generalised approach to modelling a three-tiered microbial food-web, Math. Biosci.,291, 21–37 (2017)

    work page 2017