The effect of global warming on Western Mediterranean seagrasses: towards an agent-based modelling approach

Eva Llabrés,Aina Blanco-Magadán,Marta Sales,Tomàs Sintes

DOI: https://doi.org/10.3354/meps14298

2023-05-12

Abstract:Seagrasses are key structural elements in coastal ecosystems, and studying how temperature affects these species is crucial to anticipate the implications of global warming. In this work, we use an empirically-based numerical model to study the combined dynamics of {\it Posidonia oceanica} and {\it Cymodocea nodosa} and their resilience to sea warming. The model is parametrised using seagrass growth rates measured at the Western Mediterranean Sea. Under favorable growth conditions, our simulations predict the emergence of a coexistence region at the front between mono-specific meadows. This region can be characterised by its width and local shoot densities, which are found to depend on the coupling parameter between {\it Posidonia oceanica} and {\it Cymodocea nodosa} species. Such regions have been empirically observed in Ses Olles de Son Saura (Balearic Islands, Western Mediterranean Sea). A comparison between the field measurements at the study site with the model predictions has been used to fit the value of the coupling parameter. Field data also relates the width of the coexistence region to the average length of {\it Posidonia oceanica} leaves at the front. Remarkably, a linear relationship is found between the coupling parameter and the leaf length. In the presence of sea warming, the model predicts an exponential decay in the population of {\it Posidonia oceanica}, which is highly sensitive to temperature. This behaviour is a direct consequence of the clonal nature of the plant and can be characterised by the model parameters. Considering a scenario of high greenhouse emissions, our model forecasts that {\it Posidonia oceanica} meadows will lose 70$\%$ of their population by the year $2050$. {\it Cymodocea nodosa}, with higher thermal resilience, acts as an opportunistic species conquering the space left by the degraded {\it Posidonia oceanica}.

Populations and Evolution

What problem does this paper attempt to address?

The problems that this paper attempts to solve are: **To study the impact of global warming on seagrasses in the Western Mediterranean (especially Posidonia oceanica and Cymodocea nodosa), and how these two seagrasses interact and change in the case of rising temperatures**. Specifically, the paper focuses on:

1. **Response of seagrass ecosystems to global warming**: Predict the changes of seagrass populations under different greenhouse gas emission scenarios through numerical models, especially the decline of Posidonia oceanica when the temperature rises.
2. **Interaction between the two seagrasses**: Study how Posidonia oceanica and Cymodocea nodosa compete for space in the case of rising temperatures, and explore the coupling parameters between them (i.e., the degree of influence of one seagrass on the other).
3. **Characteristics of the co - existence area**: Through model simulation and field measurement, study the width and density of the co - existence area formed by the two seagrasses at the junction, and analyze the relationship between these characteristics and temperature changes.

### Summary of the main problems in the paper

- **Threats of climate change to seagrass ecosystems**: Since the warming rate of the Mediterranean Sea is three times faster than the global average, this poses a serious threat to the local seagrass ecosystems, especially for the temperature - sensitive species such as Posidonia oceanica.
- **Modeling of multi - species interactions**: Use agent - based numerical models to simulate the growth dynamics and interactions of Posidonia oceanica and Cymodocea nodosa under different temperature conditions.
- **Predicting future changes**: Through different greenhouse gas emission scenarios (RCP2.6, RCP4.5, RCP8.5), predict the decline rate of Posidonia oceanica in the coming decades and its impact on the ecosystem.

### Formula presentation

The formulas involved in the paper include:

- **Change of the mortality rate of P. oceanica over time**:
\[
\mu(t)=\alpha t+\mu_0
\]
where \(\alpha = 0.028\lambda\), \(\lambda\) is the annual growth rate of seawater temperature (unit: °C/yr), and \(\mu_0\) is the intrinsic mortality rate (unit: yr⁻¹).

- **Population dynamic equation of a single P. oceanica patch**:
\[
\frac{dN(t)}{dt}=(\nu_0-\mu(t))N(t)
\]
where \(N(t)\) is the number of plants at time \(t\), and \(\nu_0\) is the branching rate.

- **Exponential decay of the P. oceanica population**:
\[
N(t)=N_0\exp\left[-\frac{\alpha}{2}t^2+(\nu_0 - \mu_0)t\right]
\]
where \(N_0\) is the initial number of plants.

Through these models and formulas, researchers can quantitatively evaluate the impact of global warming on seagrass ecosystems and provide a scientific basis for future conservation measures.