2.1 Study sites

Starting in February 2016, we studied the semi-enclosed coral reef lagoon of Bouraké (South Province, Grande Terre, New Caledonia) and adjacent reference reefs (Fig. 1). The mangrove forest of Bouraké is located in a semi-arid geographical area, and it lacks river input typical of mangrove estuaries.

Figure 1Map of the study site (top panel) showing the semi-enclosed mangrove lagoon of Bouraké and reference reefs (R1 and R2). Photos (a) and (b) are aerial pictures (taken at 130 m above the Bouraké lagoon) of the inner (bottom) and middle reefs (a) and the outer reef (b). White dots and yellow squares indicate the sampling stations used for physical and chemical environmental monitoring, respectively. Georep New Caledonian database (https://georep.nc/, last access: January 2021) and QGis software were used to build the figure (top panel).

In Bouraké, a channel of more than 80 m wide and 0.5 to 6 m deep penetrates a dense mangrove forest made up of Avicennia marina and Rhizophora stylosa. It expands into side pools and a large reservoir in the inner part of the mangrove system. A preliminary calculation of the surface and seawater volume of the main area (i.e., without the large seawater reservoir covered by the mangrove forest) using the 3ds Max Model gives underestimated values of 192 100 m² and 246 060 m³, respectively. Terraces extend from the mangrove forest on both sides of the channel and form diverse and compact reefs at their ends. Coral species are less abundant and diverse near the mangrove roots, where the bottom substrate is composed of fine sediment and mud. There, reefs are exposed to air only at low tide during the spring tides (1.1 m maximum tidal range). We subdivided the lagoon into three main areas: an external, an intermediate, and an inner reef (hereafter called outer, middle, and inner reefs; Fig. 1a and b). A series of sampling stations were selected, and some of the main seawater chemical and physical environmental parameters were measured at the study site and two adjacent reference reefs, namely R1, a typical fringing reef, and R2, a semi-enclosed shallow and relatively large bay. Both reefs are located 4.8 and 3.2 km from the entrance of the Bouraké lagoon, respectively.

2.2 Environmental monitoring

2.3 Benthic community characterization and distribution

The benthic community and bottom substrate of the Bouraké lagoon, referred hereafter as biotic and abiotic descriptors, respectively, were assessed in April 2018. A total of 24 30 m long georeferenced transects (T1–T24) were laid in the lagoon along the terraces' edge at similar depths (i.e., ∼1 m), targeting coral-dominated benthic assemblages. On each transect, a 0.5×0.5 m PVC quadrat was placed every meter, and a picture was taken with a waterproof photo camera (Nikon AW130) parallel to the substrate. We made a general description of the bottom (i.e., the various substrates) and a list of the most common and identifiable sessile species for each transect. For each of the 835 pictures collected, we estimated the cover of abiotic (i.e., mud, sand, rock, rubble, dead corals, and unreadable) and biotic descriptors (i.e., branching, massive and soft corals, sponges, macroalgae, and “others”) with photoQuad software both by automatic multi-scale image segmentation regions and manual grid cell counts when necessary.

We used the photos of quadrats, the many other pictures collected during fieldwork, and laboratory morphological observations on collected samples to produce a non-exhaustive species list of corals, macroalgae, and sponges. Corals were comprehensively sampled throughout the Bouraké lagoon and at the reference reef R2, while dominant macroalgal and sponge species were collected in the Bouraké lagoon alone and only if they were encountered at least three times along a transect, likely leading to an underestimation of their diversity. Coral diversity was assessed through photographic and sampling during time-based open-search swims and scuba dives (Hill and Wilkinson, 2004). Whenever possible, scleractinian corals were identified to species level in situ by photographic sampling only. When identification was doubtful in vivo, or when taxa were characterized by small corallite size (<1 mm in diameter) and required additional morphological examination of the skeleton to confirm identification (e.g., genera Acropora, Montipora, and Porites), a fragment of the colony was collected, tagged, cleaned in sodium hypochlorite overnight, rinsed in freshwater, and dried. The reference collection is housed at the Institut de Recherche pour le Développement (IRD), Nouméa. Microscopic examination of the skeletal features allowed species-level identification following the reference literature (Veron and Wallace, 1984; Wallace, 1999; Veron, 2000). Sponges were identified based on their spicules' morphological characteristics (i.e., shape, length, and width) or using a series of morphological descriptors (e.g., shape, size, color, texture, surface ornamentations, fibers) for species without spicules. In the lab, a subsample of the collected sponges was immediately digested using HNO₃, and spicules were measured with an optical stereomicroscope. For species without spicules, hand-cut sections of the choanosome and ectosome were observed under a stereomicroscope. Species were identified using the taxonomical keys of Rützler (1978), Hooper and Van Soest (2002), and Pons et al. (2017). Macroalgae and coral identification was based on morphological and anatomical observations following the dedicated literature and referring to specimens housed at IRD Nouméa.

Abiotic and biotic cover percentages averaged per transect and species richness, calculated as the number of species in the transect, were plotted using non-metric multidimensional scaling (nMDS) based on Bray–Curtis dissimilarities (“vegan” package in R) of square-root-transformed data. Finally, the best number of clusters for the whole Bouraké lagoon was determined using the gap statistic method (“cluster” and “factorextra” packages in R) and used for the hierarchical clustering representation (Ward, 1963). The cluster separation was verified with a two-way ANOSIM (ANalysis Of Similarities). Within each cluster, the benthic community and bottom substrate were averaged between transects, and the dominant biotic and abiotic descriptors were selected and used to define the cluster.

Principal component analysis (PCA) was used to visualize the correlation between environmental parameters and the benthic descriptors of transects inside the Bouraké lagoon. The analysis was performed in R (package “FactorMineR”, version 3.2.4; R Core Team, 2018) using data of biotic descriptors (averaged per transect) and the seawater parameters averaged for each sampling area (i.e., outer, middle, and inner reefs).

3.1 Variability in physical environmental parameters

Tidal phases. The predicted tide in the study area (SHOM data) is close to what we measured at R1 and R2 (Fig. 2). It is semidiurnal with a diurnal inequality (M₂=0.44 m, S₂=0.16 m) and has a small diurnal component (S₁=0.15 m). It varies between 0.4 and 1.7 m, depending on the phase of the moon, with a mean of 1.1 m. The tidal range inside the Bouraké lagoon is lower and varies between 0.4 and 1.2 m, with a mean of 0.9 m. At B1 and B2, the tidal signal has an average lag in the predicted and measured tides at the reference stations of 1.5 h for low tide and up to 45 min for high tide.

Figure 2Temporal shift between the averaged tides recorded at the reference (Sts R1 and R2) and at the Bouraké (Sts B1 and B2) reefs and the predicted tide calculated for the study area (see SHOM: http://data.shom.fr, last access: February 2021, for the Nouméa harbor).

Download

Temperature. Our results indicate that the water temperature has an annual cycle with lower values during winter (May–October) and higher values in summer (November–April). Daily averaged temperatures at R2, B1, and B2 were not significantly different, while the temperature at R1 was cooler in summer and warmer in winter (Fig. 3a; Table 1) compared to the other stations. On a weekly basis, the averaged diel temperature variation was significantly different between stations: 1.34±0.39 ^∘C at R1 and 3.73±0.74 ^∘C at B2 (Fig. 3b; Table 1). In a single day, we recorded temperature fluctuations of up to 3.4 ^∘C at R1 and 6.5 ^∘C at B2.

Figure 3Long-term temperature monitoring from October 2018 to April 2020 at the reference (Sts R1 and R2) and Bouraké (Sts B1 and B2) reefs. Data are plotted using weekly averaged temperature (a) and weekly averaged diel (b) changes.

Download

Seawater pH. During the entire study period, pH was measured during 22, 72, 31, and 72 semidiurnal tidal cycles at R1, R2, B1, and B2, respectively. We overlaid all data at a single tidal phase of a 12 h (Fig. 4a) and a 24 h cycle (Fig. 4b).

Figure 4Seawater pH_T (a, b) and dissolved oxygen (DO) (c, d) variations recorded at the reference (Sts R1 and R2) and Bouraké (Sts B1 and B2) reefs. Data were overlaid at a single tidal phase (12 h) (a, c) and at a 24 h cycle (b, d). Data are 22, 72, 31, and 72 semidiurnal tidal cycles for pH and 36, 79, 34, and 42 semidiurnal tidal cycles for DO, for R1, R2, B1, and B2, respectively.

Download

The pH differed significantly between stations. During the studied period, we recorded means of 8.01±0.04 and 7.89±0.08 pH_T units at R1 and R2 and 7.80±0.22 and 7.67±0.23 pH_T units at B1 and B2, respectively. At both B1 and B2, pH was strongly correlated with the tidal cycle (Fig. 4a; 82 % and 73 % of the total variance were explained by the tidal harmonic analysis, respectively). In contrast, it was only marginally correlated with the 24 h cycle (Fig. 4b). During each tidal phase, the pH changed on average by about 0.6 units and reached a minimum of 7.23 and a maximum of 8.06 at B2 at low and high tides, respectively (data not shown). The pH oscillations were mainly semidiurnal ($M_{\mathrm{2}}=\mathrm{0.27}\pm \mathrm{0.01}$ and $S_{\mathrm{2}}=\mathrm{0.05}\pm \mathrm{0.01}$ pH_T units for B1 and $M_{\mathrm{2}}=\mathrm{0.24}\pm \mathrm{0.01}$ and $S_{\mathrm{2}}=\mathrm{0.09}\pm \mathrm{0.01}$ pH_T units for B2). At stations R1 and R2, pH changed on average by about 0.1 pH_T units and was mostly dependent on the 24 h cycle (Fig. 4b; $S_{\mathrm{1}}=\mathrm{0.058}\pm \mathrm{0.004}$ and $S_{\mathrm{1}}=\mathrm{0.049}\pm \mathrm{0.007}$ pH_T units for R1 and R2, respectively). Simultaneous short-term pH measurements showed significant spatial differences (Fig. 5; Table 1) between (i) stations A and C, (ii) stations B1 and B2, and (iii) the outer reef in the Bouraké lagoon (St U) and the station outside the system (St T). There were also significant depth-related differences between shallow reefs and bottom water at stations R and S. In both spatial and depth-related analyses, differences were approximately 0.05 to 0.1 pH_T units, and we found the lowest values at stations C and B2 and the bottom of the channel (Sts R and S).

Figure 5Spatial and depth-related short-term changes (24 h) in pH_T at stations in the Bouraké lagoon (see Fig. 1 for the station locations). Spatial variations were assessed (a) within the inner reef (Sts A–C); (b) between the inner and the outer reefs (Sts B1 and B2); (c) between the outer reef (but still inside the Bouraké lagoon) and the reef outside the semi-enclosed lagoon (St U vs. T). Depth-related variations were assessed between shallow reefs and the bottom of the channel (d) at the inner reef (St R) and (e) at the outer reef (St S). Tide (black line) refers to St B2.

Download

Dissolved oxygen (DO). DO was measured during 36, 79, 34, and 42 semidiurnal tidal cycles at R1, R2, B1, and B2, respectively. We overlaid all data at a single tidal phase of a 12 h (Fig. 4c) and a 24 h cycle (Fig. 4d). As with pH, the mean diel DO was higher at the reference reefs than in the Bouraké lagoon. During the studied period, mean DO values were 4.89±1.18 and 5.23±0.89 mg L⁻¹ at B1 and B2, respectively, and 6.45±0.95 and 6.48±1.05 mg L⁻¹ at R1 and R2, respectively.

At stations B1 and B2, DO was strongly correlated with the tidal cycle (Fig. 4c; 82 % and 72 % of the total DO variance were explained by the tidal harmonic analysis, respectively) and only marginally with the 24 h cycle (Fig. 4d). DO oscillations were mainly semidiurnal ($M_{\mathrm{2}}=\mathrm{1.17}\pm \mathrm{0.08}$ mg L⁻¹ and $M_{\mathrm{2}}=\mathrm{0.83}\pm \mathrm{0.09}$ mg L⁻¹ for B1 and B2, respectively) with a substantial diurnal component ($S_{\mathrm{1}}=\mathrm{1.12}\pm \mathrm{0.04}$ mg L⁻¹ and $S_{\mathrm{1}}=\mathrm{0.681}\pm \mathrm{0.04}$ mg L⁻¹ for B1 and B2, respectively). During a semidiurnal tidal cycle, DO was lower at low tide (Fig. 4c; 3.7 and 4.6 mg L⁻¹ at B1 and B2, respectively) and higher at high tide (Fig. 4c; 6.0 and 6.1 mg L⁻¹ at B1 and B2, respectively). The minimum (1.89 mg L⁻¹) and the maximum (7.24 mg L⁻¹) DO values were both measured at B1 during low and high tides, respectively (data not shown). During a 24 h cycle, DO was lower in the early morning (Fig. 4d; 4.0 and 4.3 mg L⁻¹ at B1 and B2, respectively) and higher in the middle of the day (Fig. 4d; 5.8 and 5.4 mg L⁻¹ at B1 and B2, respectively). In a single day, we recorded DO fluctuations of up to 6.37 mg L⁻¹ at R1 and 4.91 mg L⁻¹ at B2. At stations R1 and R2, DO was mostly dependent on the 24 h cycle ($S_{\mathrm{1}}=\mathrm{1.34}\pm \mathrm{0.06}$ mg L⁻¹ and $S_{\mathrm{1}}=\mathrm{1.53}\pm \mathrm{0.09}$ mg L⁻¹ for R1 and R2, respectively) with lower values during the night (Fig. 4d; 5.5 mg L⁻¹ at both stations) and higher values in the middle of the day (Fig. 4d; 7.8 mg L⁻¹ at both stations).

Table 1Summary of the statistical analyses applied to the seawater physical and chemical parameters collected between 2016 and 2020 at the Bouraké lagoon (Sts B1, B2, A–R, outer, middle, and inner) and at reference reefs (Sts R1 and R2). Differences in temperature, salinity, and chemical parameters between the Bouraké lagoon and reference reefs and spatial and vertical differences in pH and DO inside the Bouraké lagoon. K–W: Kruskal–Wallis's test; W: Wilcoxon test; GLM: general linear model; D: Dunn's test; T: Tukey's test; C: Conover's test.

Download Print Version | Download XLSX

Figure 6Spatial and depth-related short-term measurements (24 h) of dissolved oxygen (DO) at stations in the Bouraké lagoon (see Fig. 1 for the location of the stations). Spatial variations were assessed at (a) the inner reef (Sts A–G), (b) the middle reef (Sts H–N), and (c) the outer reef (Sts O–Q). Depth-related variations between shallow reefs and bottom water within the Bouraké lagoon were assessed (d) at the inner reef (St R) and (e, f) at the middle reef (Sts I and N). Tide (black line) refers to St B2.

Download

Simultaneous short-term DO measurements (Fig. 6; Table 1) did not show significant differences over a 24 h cycle at the inner (Sts A–G; except for some specific deviations, e.g., St A), the middle (Sts H–N), or the outer (Sts O–Q) reef. Dissolved oxygen did not change on a vertical gradient at stations R and I, but the bottom DO was significantly lower at station N.

Salinity (S). Salinity was measured only during two short-term periods in winter 2019 and summer 2020 (Fig. 7).

Figure 7Short-term salinity monitoring during (a) the winter of 2019 (from 15 to 18 July 2019) at the reference (Sts R1 and R2) and Bouraké (Sts B1 and B2) reefs and (b) the summer of 2020 (from 29 November to 4 December 2020) at the reference (Sts R1 and R2) and Bouraké (Sts T and B2) reefs. Tide (black line) refers to St B2.

Download

Salinity variations in the Bouraké lagoon were strongly correlated with the tidal cycle, and the highest mean values were measured during both winter and summer. Salinity was lower at reference reefs during both the winter and summer seasons (Fig. 7a and b). The mean salinity during the winter of 2019 was 35.49±0.04 at R1, 35.67±0.10 at R2, 36.00±0.35 at B1, and 36.17±0.47 at B2 with significant differences between stations (Table 1). Salinity in the Bouraké lagoon peaked at 37.22 at B2 during low tide, while it was 35.65 at R1. During the summer of 2020, the mean salinity increased significantly in the Bouraké lagoon, at 37.22±0.53 at B2 reaching 39.37, while it remained lower at R1 at 35.42±0.10 and R2 at 35.71±0.16. During summer 2020, we also measured short-term variations in salinity at station T, which is outside the Bouraké system. Despite its distance from the lagoon entrance, the water mass discharged during falling tide increased seawater salinity also at this station where we measured a maximum value of 37.33 during low tide.

3.2 Diel cycles of carbonate chemistry and chemical parameters

We monitored diel cycles in June 2017 and July 2019 at two reference reefs and three reefs within the semi-enclosed lagoon of Bouraké (i.e., the outer, the middle, and the inner reefs). The seawater carbonate chemistry differed significantly between the reference and the Bouraké lagoon reefs (Figs. 8a–c and S1 in the Supplement; Tables 1 and 2).

Figure 8Seawater carbonate chemistry measured (pH_T) and calculated (DIC and Ω_ara) during diel cycles in 2017 and 2019 (pooled data; see also Fig. S1) at R1, R2, and the Bouraké reefs (outer, middle, and inner) (a–c). Boxes (n=6–14) represent the interquartile range (25th and 75th percentiles); the horizontal line is the median, and the whiskers represent the data range (i.e., minimum and maximum). Changes are illustrated across a 12 h tidal phase in the Bouraké lagoon (d–f). Linear regression lines are plotted for each falling (HT to LT) and rising (LT to HT) tide. Shaded boxes have only a graphical significance and only suggest the tide changes during measurements.

Download

The reference reefs R1 and R2 had higher pH, ${\mathrm{CO}}_{\mathrm{3}}^{\mathrm{2}-}$, and Ω_ara and lower DIC and pCO₂ than the outer, middle, or inner reefs in the Bouraké lagoon. The range of values in the lagoon was similar to our previous measurements (see above), with levels of pCO₂>2000 µatm and levels of Ω_ara<1. By assigning a tidal phase (from 0 to 12 h) to each sample taken in the Bouraké lagoon (all reefs combined), our diel measurements showed significant regressions, either positive or negative, depending on the seawater carbonate parameter and the tidal phase (Figs. 8d–f and S1; Table S3). For example, the DIC reached the highest value of 2635 µmol L⁻¹ during low tide at the inner reef with an average of 2315.1±168.1 µmol L⁻¹, while ${\mathrm{CO}}_{\mathrm{3}}^{\mathrm{2}-}$ reached the lowest value of 50.10 µmol kg⁻¹ during low tide at the inner reef with an average of 109.02±45.59 µmol kg⁻¹.

Most of the chemical parameters were, in general, more concentrated (up to 5 times) in the Bouraké lagoon than at the reference reefs R1 and R2, and they increased from the outer to the inner reef (Figs. 9a–d and S2; Tables 1 and 2).

Figure 9Orthosilicic acid, organic carbon, and nitrogen parameters measured during diel cycles in 2017 and 2019 (pooled data; see also Fig. S2) at R1, R2, and the Bouraké reefs (outer, middle, and inner) (a–d). Boxes (n=6–14) represent the interquartile range (25th and 75th percentiles); the horizontal line is the median, and the whiskers represent the data range (i.e., minimum and maximum). Changes are illustrated across a 12 h tidal phase in Bouraké (e–h). Linear regression lines are plotted for each falling (HT to LT) and rising (LT to HT) tide. Shaded boxes have only a graphical significance and only suggest the tide changes during measurements.

Download

In addition, despite the Bouraké lagoon receiving new seawater during the rising tide, only Si(OH)₄ concentrations related to the change in tide showed either a positive or negative correlation with the falling and the rising tide (Figs. 9e–h and S2; Table S3). It reached the highest value of 16.74 µmol L⁻¹ during low tide at the inner reef with an average of 11.93±4.27 µmol L⁻¹ and the lowest value of 3.38 µmol L⁻¹ during high tide at the outer reef with an average of 8.35±3.14 µmol L⁻¹, 3 to 6 times higher than at the reference stations R2 and R1. The lack of a clear effect of the tide for the other chemicals causes their accumulation inside the lagoon system. Ratios of organic and inorganic carbon, nitrogen, and nutrients showed significant differences between the Bouraké lagoon and the reference reefs (Tables 1 and 2) with higher average POC:PON, POC:DOC, Si(OH)₄:DIN, Si(OH)₄:DIP (DIN and DIP signify dissolved inorganic nitrogen and phosphorus), and DIN:DIP (${\mathrm{PO}}_{\mathrm{4}}^{\mathrm{3}-}$) and lower average DIC:DIN and DOC:TOC in the Bouraké lagoon compared to reference reefs.

3.3 Benthic community distribution and species identification

Benthic community distribution. The cluster analysis and nMDS allowed the 24 transects in the Bouraké lagoon to be grouped in six clusters corresponding to habitat descriptors, namely clusters A, B1, B2, C1, C2, and C3 (Fig. 10a; see Fig. S3 for details of the community descriptors and Table S4 for transect averaged data). The nMDS produced an exhaustive representation of the Bouraké lagoon's benthic communities and abiotic features (2-dimensional stress=0.136) confirmed by ANOSIM (global R=0.948, p=0.001). For the benthic community distribution, a clear separation (95 % confidence interval represented by the ellipse in Fig. 10a) was found for cluster A and B2, including most transects located at the middle and inner reefs (Fig. 10b).

Figure 10Benthic community and bottom substrate characterization in the Bouraké lagoon. Hierarchical Ward's method cluster dendrogram (a) and non-metric multidimensional scaling (nMDS) of abiotic and biotic percent cover and species richness per transect (T1–T24) (a). Colored circles represent clusters of transects identified in the dendrogram; circle size corresponds to species richness. Ellipses represent 95 % confidence intervals. Dominant biotic and abiotic descriptors were used to describe the benthic community and bottom substrate for each cluster in (b). The satellite image is from Georep New Caledonian database (https://georep.nc/), and the QGis software was used for transect georeferencing.

For both, branching corals were the dominant biotic descriptor with an average of 81 % and 31 % coverage, respectively (Fig. S3), with a maximum of 96 % in cluster A (T23; Table S4). Among branching corals, the two most abundant genera were Acropora and Montipora (data not shown). For the abiotic features of the substrate, dead corals and rubble characterized transects in clusters A and B2, respectively. Species richness (Fig. 10a) was high in all transects of the Bouraké lagoon except T13 and T24, both found at the middle and inner reefs' limits. The abiotic substrate of cluster C1 was characterized by 66 % mud, 12 % rocks, and 10 % sand (see Table S4 for detailed cover data per transect). Only a few branching corals (<10 %) were found, but species richness was relatively high (6 at T9 and 10 at T10). Cluster C2 was distinct, which is not surprising due to its location in a relatively shallow convergence zone that divides the lagoon into two parts. There, the substrate is made of coarse sand (13 %) and rocks (12 %) and is mainly colonized by macroalgae and sponges (31 % and 32 %, respectively). Species richness in the area was heterogeneous and ranged from 4 to 12. Dictyota spp. and Halimeda discoidea were the main macroalgal species, while Rhabdastrella globostellata was the dominant sponge species in the area. Cluster B1, located at the outer reef, is characterized by an abundance of soft corals (48 %) and rubble (21 %), as well as high biological richness (Fig. 10a). Cluster C3 is characterized by coarse sand (49 %), rocks (17 %), and a few benthic organisms such as macroalgae (10 %) and soft corals (8 %).

Principal component analysis (PCA) reduced the multicollinearity problem and the first two principal components accounted for a cumulative 71.2 % of the dataset variance. In particular. The PCA plot (Fig. 11) allowed us to assign the transects to three distinct groups depending on their position in the study area.

Figure 11Principal component analysis (PCA) between benthic community (macroalgae, sponges, corals, branching corals, and massive corals) data and environmental parameters (PON, POC, DOC, Si(OH)₄, ${\mathrm{PO}}_{\mathrm{4}}^{\mathrm{3}-}$, NO_x, and ${\mathrm{NH}}_{\mathrm{4}}^{+}$) measured for each transect (n=24). The length of the vectors indicates the contribution of each parameter. The first dimension explains 48.1 % of the variance, and the second dimension explains 23.1 % of the variance. Colors refer to the cluster subdivision (Fig. 10).

Download

The distinction among transect groups in the biplot is clear, and they are grouped based on their location at the outer, middle, or inner lagoon reefs. The outer reef sites are characterized by sponges, macroalgae, and soft corals, which appear to correlate with NO_x. Conversely, branching corals are the shared and characteristic feature of the middle and inner reef sites. However, these two reef groups differ in their major correlation to ${\mathrm{NH}}_{\mathrm{4}}^{+}$ for middle reef sites and POC, DOC, Si(OH)₄, and ${\mathrm{PO}}_{\mathrm{4}}^{\mathrm{3}-}$ for the inner sites.

Species identification. A total of 66 coral species were identified from the Bouraké lagoon (Table S5). Most of the species belong to the family Acroporidae (17 Acropora species, 2 Anacropora, and 4 Montipora) and Merulinidae (15 species). In total, 58 species were identified from the reference reef R2, 29 of which were also found in the Bouraké lagoon. Overall, the hard coral fauna within the lagoon was characterized by species commonly found around Grande Terre, with some of them typically found in turbid coastal environments such as Pseudosiderastrea tayamai and Heliofungia actiniformis. They were not encountered at the reference site. Remarkably, two New Caledonian endemic coral species, Cantharellus noumeae and Polycyathus fulvus, both described from coastal habitats characterized by terrigenous sediment inputs (Nouméa lagoon and Prony Bay, respectively), were common at the innermost reef of the Bouraké lagoon. The Bouraké lagoon also had 28 species of macroalgae (14 Phaeophyceae, 9 Chlorophyta, and 5 Rhodophyta), including the crustose coralline algae Lithothamnium sp., Lithophyllum sp., and Hydrolithon reinboldii, and 11 sponge species belonging to eight families of the class Demospongiae (Tables S5 and S6), which were found mainly in the coral matrix and sediment.

4.1 Physical and chemical characteristics of the Bouraké lagoon

The Bouraké lagoon covers an estimated area of ca. 20 ha that penetrates a mangrove forest, which is large enough to assess the combined effects of extreme environmental conditions at an established coral reef ecosystem. Our multi-scale approach confirmed previous findings (Camp et al., 2017), showing that the Bouraké lagoon is hottest (ranging from 17.5–33.8 ^∘C), more deoxygenated (ranging from 1.87–7.24 mg L⁻¹), and more acidic (ranging from 7.23–8.10 pH_T units) when compared to neighboring reefs. Besides, we found that salinity was significantly higher than at the reference reefs during both winter and summer (with maxima of 37.22 and 39.37, respectively).

We found several marked differences in the environmental conditions between the Bouraké lagoon and the reference reefs, both in the absolute range and in the variability in the measured environmental parameters.

First, the seawater temperature is higher in summer in the Bouraké lagoon (Fig. 3), but it is also colder during winter, resulting in an annual temperature range of 17.5–33.8 ^∘C. We compared temperatures recorded at Bouraké lagoon to those of the reference St R2, which showed the most typical temperature range for shallow water temperatures in the south of New Caledonia (i.e., 22–28 ^∘C; Varillon et al., 2021). We notice that in Bouraké, temperatures were 40 % of the time above 28 ^∘C during the summer of 2020, while winter temperatures were on average 46.5 % of the time lower than 22 ^∘C. While warming is considered the main threat for coral reefs, low temperatures (<20 ^∘C) can cause coral bleaching by inducing responses similar to high temperatures, including a reduction in the Symbiodiniaceae cell density and chlorophyll a content (e.g., Saxby et al., 2003; Hoegh-Guldberg and Fine, 2004; Hoegh-Guldberg et al., 2005; Kemp et al., 2011; Bellworthy and Fine, 2021). The negative effect of cold temperatures is even more substantial during neap tides when colonies on the reef crest are exposed to air for hours at low temperatures during cold winters. For example, Porter et al. (1982) and Davis (1982) reported >90 % coral mortality in shallow (<2 m) reefs of the Dry Tortugas following the winter of 1976–1977 when temperatures reached 14 ^∘C. This is consistent with our observations (data not shown) during the cold winter of 2019 when we found that the upper 12–20 cm of several massive and branching corals had died.

Second, in the Bouraké lagoon, benthic assemblages are continuously exposed to large fluctuations in the main environmental parameters toward suboptimal values. Some of these environmental fluctuations are entirely predictable. For instance, marine organisms are exposed to a temperature fluctuation of about twice the reference reefs' amplitude (up to 6.5 vs. 3.5 ^∘C) in a single day. Dissolved oxygen fluctuations were similar between stations but in a significantly different range: 3.7 to 6.8 mg L⁻¹ at B2 and 5.4 to 7.8 mg L⁻¹ at the reference stations. According to the organisms tolerance to DO fluctuation, which is quite unknown in corals, low DO concentrations can change fish tidal migration in the mangrove (Dubuc et al., 2019). Besides, our tidal modeling revealed that, at the reference reefs, pH, DO (Fig. 4), and temperatures (data not shown) slightly increased in the afternoon and decreased during the night. This finding agrees with what should be expected from reef metabolic activities and daily cycles, but, in the Bouraké lagoon, these parameters, including salinity (Fig. 7), are entirely driven by tides. Here, seawater pH and DO varied between extremely low values at low tide and close-to-normal values during high tide (see also Fig. 8a–c for pH). Finally, we found that the timing of the tide was out of phase between sites, with a delay of about 45 min at high tide and 1.5 h at low tide in the Bouraké lagoon (Fig. 2).

The unique environmental conditions measured in the Bouraké lagoon are linked to its unique topographical and geomorphological characteristics, the resulting water circulation, and the direction of the tide. New water from the lagoon enters through the channel at each rising tide and flows into the semi-enclosed lagoon towards the large mangrove area behind it. This water initially had ambient values of pH, temperature, and dissolved oxygen, but, during the trip, it mixes with the acidic, warm, and deoxygenated water in the system and the mangrove area, therefore gradually changing from its original values. Inside the mangrove forest, we hypothesize that the water chemistry further changes due to the metabolic reactions in the sediments and mangrove roots (e.g., Alongi et al., 2004; Bouillon et al., 2007; Gleeson et al., 2013; Call et al., 2015). Conversely, for a falling tide, the seawater becomes gradually more acidic, hot, and oxygen-depleted because the water that resided in the mangrove area gradually drains out of the system. This takes about 6 h, during which the vast reservoir of shallow mangrove water continues to be chemically altered, becoming increasingly acidic, oxygen-depleted, and hot. As a result, we measured significant spatial differences in pH between the outer reef (the entry of the lagoon) and the inner reef (near the mangrove forest), as well as a considerable delay in the synchronization of the tidal shift (Fig. 5b). Interestingly, because the volume of seawater discharged in 6 h is so large, it affects also the area outside the system where we measured similar seawater conditions as inside, even if it mixes with the main lagoon's water (see Fig. 5c, St U vs. St T). It means that the area (and the organisms) affected by the suboptimal parameters is larger than previously thought. The species living in this area have likely developed specific mechanisms to withstand the drastically fluctuating environmental conditions, and as such, they warrant further attention.

Since the fluctuations are linked to tidal phases, it could be argued that organisms living in the Bouraké lagoon may benefit from periods of normal conditions at high tide, during which they can recover from the stress they have experienced at low tide (e.g., Rivest et al., 2017). While this could be partially the case for species living on the outer reef, close to the main lagoon, the environmental conditions inside the Bouraké lagoon rarely reach normal values (Fig. 8a–c) and persist longer as the low tide is delayed by 1.5 h compared to the reference reef (Fig. 2).

Preliminary results from a hydrodynamic model of the study site suggest that tide-associated water mass movements are spatially heterogeneous and likely to play an essential role in shaping coral resilience to extreme conditions (see discussion below). Indeed, one can imagine a single water mass moving with the same physical characteristics from the mangrove area towards the outer reefs or in the opposite direction depending on the tide. However, the complex geomorphology of the Bouraké lagoon, its bottom topography, and the complex web of coral reefs and mangrove trees on the edges deviate and probably change the seawater's physical and chemical properties. We measured significant spatial differences in pH within each reef area (inner, middle, and outer reefs; Fig. 5a and b), as well as throughout the water column (i.e., between the surface and the bottom; Fig. 5d and e). In general, bottom seawater was 0.1–0.2 pH_T units lower than the surface probably due to a pumping mechanism by the water mass of more acidic pore water from the sediments. The pH also differed spatially within the inner reef by up to 0.3 pH_T units (for instance St A vs. C; Fig. 5a), perhaps due to the influence of stagnant water on the mangrove edges or a more intense metabolic activity by the local benthic community. Dissolved oxygen concentrations between the surface and the bottom were only significantly different at St N at the middle reef (Fig. 6f). This station is mainly characterized by mud as the current is lower than anywhere else in the Bouraké lagoon. Here, the sediment's biological activity possibly traps the oxygen making it less available to the water column.

Table 2Summary of the seawater physical and chemical data (mean±SD) measured in July 2017 and June 2019 and calculated using the CO₂SYS package. Data were pooled and averaged per station; nd=not determined. Reference reefs: Sts R1 and R2; Bouraké lagoon stations: outer, middle, and inner.

Download Print Version | Download XLSX

We can assume that, throughout the Bouraké lagoon, organisms are exposed to extreme and fluctuating suboptimal physical parameters, such as pH, and DO which are more pronounced on the bottom and last longer – and with more extreme values – at the inner reefs. It is also true for some of the seawater chemical parameters we measured, which show higher concentrations in the Bouraké lagoon than at the reference reefs (see Table 2). For instance, we found that orthosilicic acid, phosphate, dissolved and particulate organic carbon, and particulate organic nitrogen were 1.4- to 4.2-fold more concentrated at the inner reef than at station R2. Within the Bouraké lagoon in general, chemicals increased from the outer to the inner reef, and except for [NO_x], [NH₄]⁺, and [PO₄]³⁻, they never return to “normal” values despite the Bouraké system receiving new seawater during the rising tide. The high concentrations in orthosilicic acid and organic compounds, both dissolved and particulate, are due to the combined effect of more acidic and organic-rich water coming out of the mangrove forest during a falling tide (Fig. 9) and the system's particular morphology, limiting the complete renewal of the seawater, especially at the inner reefs.

4.2 Effects of physical and chemical conditions on species distribution

Extreme environmental conditions, such as those measured in the Bouraké lagoon, are known to strongly affect the metabolism, growth, and even survival of several marine organisms, particularly those tolerating only a narrow range of environmental changes such as scleractinian corals (e.g., Coles and Jokiel, 1977; Hoegh-Guldberg and Smith, 1989; Hoegh-Guldberg, 1999; Fitt et al., 2001). For instance, in situ studies at volcanic CO₂ seeps have shown that chronic exposure to ocean acidification (near-future pCO₂ levels) can cause a reduction in coral diversity and lower the recruitment and abundances of structurally complex hermatypic corals. Moreover, shifts in competitive interactions between taxa and a decrease in cover and richness of soft corals and sponges were also observed (e.g., Fabricius et al., 2011; Enochs et al., 2015; Sunday et al., 2017; Agostini et al., 2018). However, the Bouraké lagoon features high coral, sponge, macroalgae, and crustose coralline algae (CCA) cover and species richness, adding to the checklist proposed by Camp et al. (2017) an additional 46 coral species (totaling 66 species and 33 genera), 28 species of macroalgae, and 11 species of sponges (Table S5). Such a high number of species has never been reported for mangrove systems characterized by marginal conditions, at least concerning the three major benthic group taxa investigated here (e.g., Yates et al., 2014; Camp et al., 2019). In general, given the extensive literature demonstrating the negative effects of suboptimal abiotic conditions on corals, our findings are unexpected and suggest that corals in the Bouraké lagoon may have developed unique survival and adaptive mechanisms. Some studies have reported similar findings, but none have ever reported such high coral species diversity (e.g., Yates et al., 2014; Schoepf et al., 2015; Shamberger et al., 2018; Camp et al., 2019). For instance, in Palau (Micronesia), the highest coral cover (>60 %) and species diversity (21 scleractinian genera) were found at the lowest pH study site of the Nikko Bay station (Barkley et al., 2015). Studies in the Virgin Islands' mangrove system have reported over 30 coral species growing in marginal conditions (Yates et al., 2014; Rogers, 2017). Similarly, in a recent study of two Australian mangrove lagoons, characterized by relatively extreme environmental conditions (low pH, low oxygen, and high temperature), Camp et al. (2019) identified 12 and 29 scleractinian coral species in the Woody Isles and at Howick Island, respectively. Among the 66 species we identified, 29 of the most abundant coral species in New Caledonia were found both at the reference reef and in the Bouraké lagoon, indicating that species living at our study site are not very different from a typical bay-sheltered fringing reef. There is not an apparent selection in the coral form since massive (e.g., Porites spp.), corymbose (e.g., Pocillopora digitata), phaceloid (Galaxea fascicularis), and branching (Acropora spp.) species are equally abundant, as well as in their thermal and pH tolerance. Further experiments are underway at both the molecular and phenotypic level to better understand the mechanisms of resilience used by the corals from Bouraké.

Remarkably, the two coral species, currently considered New Caledonian endemics, thrive in the innermost benthic assemblages of the Bouraké lagoon, making it not only a unique natural laboratory but a potential conservation priority site. In the Bouraké lagoon, benthos species richness was very high throughout the system, and the distribution of branching and massive corals was spatially heterogeneous (Figs. 10b and S3). Coral cover decreased near the mangrove forest, where the accumulation of fine sediments and exposure to air during low tide do not allow corals' survival (Fig. S3). Branching corals were particularly abundant at the inner and middle reefs, with the highest cover of 96 % at the inner reef (T23; Table S4). They became sparse at the outer reef, even absent at the system's entrance where soft and massive corals dominated. Macroalgae and sponges, including CCA (data not shown), were quite cryptic throughout the Bouraké lagoon but found almost everywhere in the coral matrix and buried in the sediment. They are particularly abundant in a shallow area that divides the lagoon into two parts (transects T6–T8). There, a Venturi effect generates a strong current, and the substrate consists of coarse sand and rocks. The high flow and the bottom characteristics are both likely to promote macroalgae and sponge occurrences (averaged 30 % and 32 % cover, respectively). The sponge Rhabdastrella globostellata is particularly abundant (up to 40 % cover) and forms massive banks embedding adjacent coral colonies. In general, macroalgae and sponge diversity was not particularly rich, with only 28 and 11 species, respectively, likely due to the sampling method. However, in the Indo-Pacific mangrove systems, sponge diversity is, in general, not extremely high, especially if compared to the Caribbean (Nagelkerken et al., 2008). Few studies are available for Indo-Pacific mangrove sponges. The highest diversity was reported in a study in the mangrove forest of Bangka Island (North Sulawesi, Indonesia), where 19 species were found (Calcinai et al., 2017). Our findings on mangrove sponges in such a unique mangrove area in New Caledonia add to the limited knowledge of sponge species diversity in Pacific mangrove systems.

Water flow may play a critical role in the response of organisms to acidification (Cornwall et al., 2014; Hurd, 2015; Comeau et al., 2019), warming (Schoepf et al., 2018), and deoxygenation (Hughes et al., 2020). It has been reported to affect the diffusion boundary layer (i.e., a thin layer of stagnant water located around aquatic organisms) of corals, CCA, and other calcareous macroalgae, altering their ability to calcify at low pH and to excrete metabolic wastes (reviewed in Nelson and Altieri, 2019). Knowing that flow speed could play a critical role which deserves further investigations and measurements at the study site, we found two hypotheses to explain the high diversity of the Bouraké species and their resilience to suboptimal parameters: (i) the species may benefit from the constant fluctuations of the physical and chemical parameters, and (ii) the species may benefit from heterotrophic inputs of the mangrove (in term of nutrients and organic matter).

Regarding environmental fluctuations, we measured averaged daily variations in temperature, DO, and pH of about 3.79 ^∘C, 3.80 mg L⁻¹, and 0.63 pH_T units, respectively. The temperature fluctuated by up to 6.55 ^∘C in a single day. Frequent exposure to stressful temperatures can induce acclimatization or adaptation in corals if the period of exposure is short enough to avoid mortality (Oliver and Palumbi, 2011; Palumbi et al., 2014; Schoepf et al., 2015; Rivest et al., 2017; Safaie et al., 2018). For example, in the Kimberley region (Australia), corals regularly exposed to temperatures up to 37 ^∘C, with daily variations of up to 7 ^∘C, appear less prone to bleaching and grow and calcify at rates comparable to corals in more thermally stable environments (Dandan et al., 2015; Schoepf et al., 2015, 2020). However, they are not immune to extreme heat stress events (Le Nohaïc et al., 2017). Another example is the corals at the back reef and shallow pools around the island of Ofu (American Samoa), which despite occasional daily fluctuations of up to 6 ^∘C and an average daily temperature of 29 ^∘C in summer (Piniak and Brown, 2009; Oliver and Palumbi, 2011) sustain reasonable levels of coral cover (25 %–26 %) and high diversity (Craig et al., 2001). Some studies tend to indicate that corals depended on the a priori “frontloading” of genes involved in heat resistance in the coral host and/or the host's ability to adjust its physiology during short-term (weeks) environmental changes (Barshis et al., 2013; Palumbi et al., 2014). However, physiological adjustments have an energetic cost that corals cannot sustain in the long term, affecting other metabolic functions such as calcification. It could explain in part the lower calcification rates observed in some coral species, as seems to be the case for corals in the Bouraké lagoon (e.g., Camp et al., 2017). There is also evidence that exposure to short-term oscillations in low pH, alone or in combination with high temperature, can mitigate the impact of extreme environmental stressors on corals (Warner et al., 1996; Oliver and Palumbi, 2011; Dufault et al., 2012; Schoepf et al., 2015; Safaie et al., 2018). Species in the Bouraké lagoon appear to be thriving despite the large seawater pH variability (up to 0.6 pH_T units). In general, species exposed to pCO₂ fluctuations above 500 µatm in their natural environment display enhanced plasticity to future ocean acidification scenarios (Vargas et al., 2017; Carstensen and Duarte, 2019). For example, fluctuations in pH could increase CCA tolerance to ocean acidification by providing respite periods at ambient pH (Rivest et al., 2017), during which organisms can calcify at a standard rate, compensating for decreased calcification during periods of lower pH (Comeau et al., 2013; Cornwall et al., 2018). Diurnal variability in pH, typical of more variable pH habitats, may confer tolerance to resident calcifying species via the selection of individuals better adapted to survive in these environments (Rivest et al., 2017; Kapsenberg and Cyronak, 2019), although trans-generational acclimatization is required (Cornwall et al., 2020).

In addition to changes in temperature and pH, which are the two most commonly tested environmental parameters, oxygen must also be considered. Despite the impact of reduced oxygen levels on and natural fluctuations in coral reefs having yet to be established, low DO (2–4 mg L⁻¹) seems to increase the sensitivity of branching corals, resulting in a decline in coral health, bleaching, and tissue loss (Haas et al., 2014; Hughes et al., 2020; Alderdice et al., 2021). Our study shows that although the Bouraké system can reach conditions close to hypoxia for several coral species (<3 mg L⁻¹; Fig. 4), these latter seem to have promoted compensation mechanisms that allow them to survive in these conditions. The natural laboratory of Bouraké, where DO fluctuates with the tide, in combination with other environmental stressors, offers a perfect setting to test the practically unknown effects of deoxygenation in reef-building corals exposed to acid and hot conditions (Nelson and Altieri, 2019; Hughes et al., 2020).

Although corals appear to possess cellular mechanisms to counteract short-term osmotic changes (Mayfield and Gates, 2007), high and fluctuating salinity is a possible additional stress that corals living in the Bouraké lagoon have to face daily, adding up to the already long list of suboptimal environmental parameters. In situ studies have demonstrated that acute and prolonged decrease in salinity can affect the coral photosynthetic efficiency, resulting in a reduction of the amount of energy transferred to corals (Muthiga and Szmant, 1987; Manzello and Lirman, 2003), and induce coral death (e.g., Jokiel et al., 1993). The response of corals to a change in salinity is related to the strength and duration of the hypo-hypersaline exposure and the species tolerance. For example, Stylophora pistillata seems able to acclimate more effectively to hypo- rather than to hyper-saline conditions (Ferrier-Pagès et al., 1999). Further experiments are needed to assess the effect of high and fluctuating salinity on the physiology of corals in Bouraké.

Besides the hypothesis that environmental variability improves the metabolism of organisms, particularly their resilience to extreme conditions, a series of other physical (e.g., current flow) and chemical parameters (e.g., organic matter) in the Bouraké lagoon may work in combination to offset or enhance these effects. Mangrove habitats are highly productive ecosystems and are sites of intense carbon processing, with a high potential impact on the global carbon budget (e.g., Borges et al., 2003; Dittmar et al., 2006; Bouillon et al., 2014). In the Bouraké lagoon, benthic communities might have access to a range of heterotrophic inputs, nutrients, carbon, and nitrogen sources. These sources can be metabolized by the species to increase their energy budget and cope with the suboptimal parameters, but they can also become toxic if too concentrated or depleted, leading to functional limitations. We measured exceptionally high concentrations of organic and inorganic carbon and nitrogen but also of some nutrients, notably silicates and phosphorus, and we confirmed the potential contribution of the mangrove in those inputs, especially during the falling tide (Figs. 8, 9, S1, and S2; Table 2). We found that dissolved organic carbon contributes significantly to the TOC pool (POC+DOC), with a concentration increase between the reference sites and the Bouraké lagoon, reaching maximum values at the inner reef. High organic matter content can increase DOC availability to corals, providing the sustainable energy to withstand extreme environmental conditions (Levas et al., 2015). Some studies showed that the high organic matter of turbid reefs can support elevated coral heterotrophy that can facilitate energy maintenance during periods of stress (Anthony and Fabricius, 2000).

We also found that nutrients could partially explain the distribution of organisms throughout the Bouraké lagoon (Fig. 11). Indeed, NO_x concentrations were higher at the outer reef, dominated by sponges, macroalgae, and soft corals, while ${\mathrm{NH}}_{\mathrm{4}}^{+}$ was high at the middle reef, and POC, DOC, Si(OH)₄, and ${\mathrm{PO}}_{\mathrm{4}}^{\mathrm{3}-}$ were higher at the inner reef. Both the middle and inner reefs are characterized by the highest branching coral cover. Nitrates can accumulate inside the host cells, possibly favoring sponges and macroalgae, while ammonium is the preferred source used by coral symbionts (Raven et al., 1992). However, the nutrient negative or positive effect on corals' physiology is difficult to demonstrate experimentally (Atkinson et al., 1995; Szmant, 2002; Bongiorni et al., 2003). Van De Waal et al. (2009) have shown that the performance of organisms depends on concentrations of, and ratios between, different nutrients, as well as between organic and inorganic matter/components and the possible imbalance due to environmental changes. We found that the Bouraké lagoon is mostly N-limited ($\text{DIN}:\text{DIP}<\mathrm{10}$ and $\mathrm{Si}(\mathrm{OH}{)}_{\mathrm{4}}:\text{DIN}>\mathrm{1}$), which confirms the findings of Justić et al. (1995) and is similar to the conditions observed in other New Caledonian mangroves (Jacquet et al., 2006). Nutrient limitation has been demonstrated to lower the temperature effect at which coral bleaching occurs (Wiedenmann et al., 2013; Ezzat et al., 2016a, 2019), which contrasts with the resilience of Bouraké corals to the warming in the summer of 2016 (10 %–20 % bleaching only) compared to other reefs in New Caledonia (up to 90 % bleaching) (Benzoni et al., 2017). Coral symbionts recycle their host's metabolic wastes and take up dissolved inorganic nitrogen (DIN) and phosphorus (DIP) from seawater (Grover et al., 2003; Pernice et al., 2012; Rosset et al., 2015), both of which are used to produce vital organic molecules. Nutrient starvation can occur when the availability of one type of essential nutrient decreases, resulting in an imbalanced N:P ratio of inorganic nutrients in seawater (Wiedenmann et al., 2013; D'Angelo and Wiedenmann, 2014). Based on our measurements, the N:P ratio in the Bouraké lagoon was 2.8:1 and 4.6:1 at the middle and the reference reefs, respectively, which is lower than the range calculated for average reef waters of 4.3:1 to 7.2:1 (Crossland et al., 1984; Furnas et al., 1995). Although these data should be considered cautiously, they suggest that the Bouraké lagoon seawater is not limited in phosphorus, an essential nutrient in coral resilience to bleaching (Ezzat et al., 2016b; Rosset et al., 2017).

Seawater in the Bouraké lagoon was not limited in Si ($\mathrm{Si}(\mathrm{OH}{)}_{\mathrm{4}}:\text{DIP}>\mathrm{10}$ and $\mathrm{Si}(\mathrm{OH}{)}_{\mathrm{4}}:\text{DIN}>\mathrm{1}$; Justić et al., 1995), and Si was 6 times more abundant at the inner reef than elsewhere in the system (Fig. 9; Table 2). Silicates are an important source of nutrition and skeletal construction for primary producers such as sponges. Orthosilicic acid is a biologically available form of silicon that is poorly soluble at a very low pH. In the Bouraké lagoon, pH fluctuates with the tide, and one can expect that orthosilicic acid would rise during the falling tide. This elevated concentration of orthosilicic acid could be involved in the growth and persistence of the Bouraké large banks of Rhabdastrella globostellata, as previously reported for other sponges in the shallow hydrothermal vents in the Pacific Ocean (Maldonado et al., 1999; Cárdenas and Rapp, 2013; Bertolino et al., 2017).

We are still missing information about light, turbidity, current, bacteria biomass, phytoplankton sources, and other biological communities to fully understand this complex and dynamic system and the functioning of this complex and dynamic coral ecosystem. However, our data already demonstrate the Bouraké lagoon's uniqueness as a natural laboratory for studying the adaptive responses of corals and other reef species to the combination of multiple suboptimal environmental parameters, which are, to some extent, worse than those projected for the future. Our investigations indicate that the geomorphology of this unique site has not changed for the last 80 to 100 years and certainly longer. With this in mind, we assume that, at least in the previous century, the environmental conditions remained unchanged, and the corals of the Bouraké lagoon have experienced the current extreme conditions for several generations. Most importantly, the species found in the Bouraké lagoon are not unique to the mangrove habitat. They are common throughout New Caledonia and the southwest Pacific region, suggesting that they have used specific strategies to cope with the suboptimal environmental conditions.

We are also unsure how the different parameters will fluctuate at bay-sheltered reefs over time and under future climate change conditions and if the variability measured in the Bouraké lagoon is representative of the natural fluctuations expected for coastal habitats in the future. Indeed, earth system model simulations, which are mainly based on open-ocean-system models, project that the seasonal amplitude of pH and pCO₂ will increase by 81 % in the future climate (Kwiatkowski and Orr, 2018). Projections also indicate that ocean acidity extremes will be more frequent (Burger et al., 2020), which could mean that future diel pH variability will increase even more at sites with the most significant variability today. The situation could be similar for seawater temperature, a parameter that is already high in Bouraké and likely close to the coral's thermal tolerance threshold. If this were to happen in the following decades, the coral reef of Bouraké would face incredibly harsh environmental conditions since physical and biogeochemical parameters measured during this study already exceed future climatic simulations.