Model shows abrupt loss of soil organic carbon following disturbance in seagrass ecosystems

Code for running the bifurcation diagrams and the sensitivity analysis of seagrass-soil model Contact: [e-mail] This repository contains the code used to conduct the figures of: Model shows abrupt loss of soil organic carbon following disturbance in seagrass ecosystems. Each figure from the related study has its own folder, which includes the necessary scripts to rerun simulations, the output of those simulations, and the code to plot the results. By navigating to any figure’s folder, you can reproduce the simulations and visualise the results. The repository is organised to facilitate reproducibility and further exploration of the ecosystem model and its behavior under various scenarios. We performed our analysis using R version 3.6.3. Do not forget to add your working directory if you want to save the figures. Sensitivity analysis (Figure 4, Figure S11, Figure S15 & Figure S16) The “sensitivity” folder contains subfolders with the scripts required to run the global sensitivity analysis using the Sobol method for each scenario/case, along with the resulting outputs. In this analysis, higher numbers in folder names indicate a more deteriorated meadow, meaning it’s closer to the point of collapse. The analysis was conducted across different scenarios for different cases: “f” denotes the feedback case, while “no_f” represents the no feedback case. To recreate the figures, you can plot the pie charts for each scenario/case by running the sensitivity_plot.R script after setting the working directory to the appropriate subfolder.海草-土壤模型分岔图与敏感性分析运行代码 联系方式：[e-mail] 本仓库包含用于复现研究《扰动后海草生态系统土壤有机碳突降模型》相关图表的代码。该研究的每张图表均对应独立文件夹，内含重新运行模拟所需的脚本、模拟输出结果以及绘图代码。进入任意图表对应的文件夹，即可复现模拟过程并可视化结果。本仓库的组织形式旨在提升研究可重复性，并便于进一步探索该生态系统模型在不同场景下的行为表现。 本次分析基于R 3.6.3版本完成。 若需保存生成的图表，请务必设置工作目录。 敏感性分析（对应图4、补充材料图S11、S15及S16） "sensitivity"文件夹下包含多个子文件夹，每个子文件夹均配有针对特定场景/案例的索伯尔（Sobol）全局敏感性分析运行脚本，以及对应的分析输出结果。本分析中，文件夹名称的数值越高，代表海草床退化程度越严重，即越接近崩溃临界点。本次分析针对不同场景与案例展开："f"代表存在反馈的场景，"no_f"则代表无反馈的场景。若需复现相关图表，可将工作目录切换至对应子文件夹后，运行sensitivity_plot.R脚本以生成各场景/案例的饼图。