The Salvatori 2022 model — "A System Dynamics Approach to Model Photosynthesis at Leaf Level Under Fluctuating Light" — is a deliberately coarse-grained, system-dynamics description of leaf photosynthesis under changing illumination. Rather than resolving every carrier in the electron transport chain, it lumps the process into just six state variables: the active fraction of photosystem II (PSII), an electron-transport intermediate (Q), the non-photochemical quenching state (P_NPQ), the NADPH/NADP+ redox pool, and the activation state of RuBisCO (R). Incident light drives PSII, electrons flow through the chain to build up reducing power, cyclic electron flow and the resulting transmembrane ΔpH feed ATP production and activate the NPQ machinery, and the NADPH/NADP+ ratio in turn activates RuBisCO, which assimilates CO₂.
The model's goal is to capture the dynamic, leaf-level response of photosynthesis to fluctuating light with as few components as possible, linking the fast photoprotective regulation (NPQ) to the slower activation of carbon fixation. Because it abstracts the light and carbon reactions into a small, transparent set of coupled fluxes, it is included in GreenSloth as a minimal, high-level alternative to the more mechanistically detailed photosynthesis models — useful for illustrating the overall regulatory logic of the system without the complexity of a full kinetic electron-transport chain.
Scheme
Model definition
Variables
| Symbol | ID | Initial value | Unit |
|---|---|---|---|
| E_PSII | 0 | ||
| Q | 0 | ||
| P_NPQ | 0 | ||
| NADP | 5 | ||
| NADPH | 5 | ||
| R | 0.001 |
Parameters
| Symbol | ID | Value | Unit |
|---|---|---|---|
| PAR | 0 | ||
| alpha | 0.78 | ||
| c_in | 0.23 | ||
| Ec_PSII | 157.56 | ||
| v_ETR | 0.78 | ||
| v_d | 0.08 | ||
| Qc | 0.07 | ||
| v_NPQ | 70.58 | ||
| v_p | 0.07 | ||
| v_C | 11.75 | ||
| eta_NADPH | 5.07 | ||
| eta_NADP | 0.89 | ||
| v_R | 0.00089 | ||
| d | 8.4 | ||
| c_y | -4 |
Derived quantities
| Symbol | ID | Equation | Unit |
|---|---|---|---|
| pH | |||
| ETR | |||
| A |
Reactions
| Symbol | ID | Rate | Stoichiometry | Unit |
|---|---|---|---|---|
| Energy_input | ||||
| ETR_out | ||||
| ETR_in | ||||
| energy_dissipation | ||||
| NPQ | ||||
| NPQ_activation | ||||
| Carbon_assimilation | ||||
| RuBisCO_activation |
Curation
- July 2026: Maintenance by Marvin van Aalst
- June 2026: Initial implementation by Josha Ebeling
- June 2026: Validation by Elouen Corvest
This model was validated by reproducing the following figures of the original publication. Due to copyright restrictions, a side-by-side comparison of the original figures is not shown here.
The publication did not provide any code, but the differential equations are sufficient to reproduce the model and results. The description of the simulation performed is clearly stated, so the figures were replicated with high accuracy.