Carlot et al., (2026) Unravelling marine benthic functioning shifts under ocean acidification – data

Data availability.
The raw observational data underlying this study are subject to a temporary embargo, as they form part of the ongoing Beatriu de Pinós grant (2024BP00106; expected completion: 31 August 2027). To ensure full reproducibility, we provide the complete processed datasets, model outputs, and fully annotated analysis code used to generate all results, tables, and figures in the manuscript. Upon completion of the grant, the raw data will be made publicly available in this repository.

Figure reproducibility.
All figures and figure panels presented in the manuscript were graphically post-processed for publication purposes. As a result, figures generated directly from the provided code may differ slightly in layout or formatting from those shown in the published article. However, the code reproduces all figures in their analytical (raw) form and fully replicates the underlying results.

Repository structure.
📁 Data is the folder where you might find the data needed to reproduce the figures.
• Figure 1 • Figure 2 • Figure 3 • Figure 4 •

📁 R_Script hosts the scripts used for the current analyses.
The analysis scripts are indexed and correspond to the different phases of the project conducted over the past two years. Once the raw data are released, these scripts will serve as the complete processing and analysis pipeline, from raw measurements to final results.
Each script is documented and linked below:\

• 0. Titration_alkalinity.R: Quantifies total alkalinity (AT) from processed titration samples.
• 1. Quality_Check_O2_Sensors.R: Formats and performs quality control on oxygen data obtained from sensors.
• 2. MiniDots.R : Extracts oxygen data to estimate dark respiration (DR) and gross photosynthetic rates (GP) of the samples.
• 3. Alkalinity.R: Converts changes in total alkalinity into calcification rates.
• 4. PI_Photo.R: Derives the relationships between calcification rate (CR), dark respiration (DR), gross photosynthesis (GP), and irradiance (photosynthesis–irradiance curves).
• 5. PAR_profiles.R: Extracts light (PAR) profiles recorded during incubation experiments.
• 6. Vizualisation.R: Visualizes the three pH conditions across the three study sites (Figure 1).
• 7. Nutrients.R: Extracts and processes nutrient uptake data.
• 8. Covers_communities.R: Quantifies community cover on experimental tiles across different sampling times.
• 9. Biomass_communities.R: Estimates community biomass on experimental tiles across sampling times (Figure 2).
• 10. Conditions_T0.R: Defines initial (pre-experimental) conditions (Table 1).
• 11. Changes_communities.R: Quantifies temporal changes in ecosystem functions.
• 12. Long-term: Assesses long-term functional changes on historical tiles (Figure 4).
• 13. Review: Addresses reviewer comments from Round 1, including the generation of the revised Figure 3.

As the raw observational data are not yet publicly available, the following scripts correspond to the figures in the manuscript and currently work with the processed data:
• Figure 1 • Figure 2 • Figure 3 • Figure 4 •

📁 Outputs hosts the main figures.
• Figure 1 • Figure 2 • Figure 3 • Figure 4 •