High Coral Recruitment Despite Coralline Algal Loss Under Extreme Environmental Conditions

The crucial role of crustose coralline algae (CCA) in inducing hard coral larval settlement and ensuring the replenishment of coral reefs is widely accepted, and so are the negative effects of anthropogenic CO2 emissions on both CCA abundance and coral development. However, diversified and well-developed coral reef communities have been recently discovered in natural conditions where CCA and corals would not be expected to thrive. Back-reef pools, volcanic CO2 vents, mangrove estuaries, and semi-enclosed lagoons systems can present seawater pH, temperature, and dissolved oxygen values reaching or even exceeding the conditions currently predicted by the Inter Panel on Climate Change (IPCC) for 2100. In the semi-enclosed lagoon of Bouraké (New Caledonia, southwest Pacific Ocean), seawater pHT, dissolved oxygen, and temperatures regularly fluctuate with the tide reaching respectively minimum values of 7.23 pHT units, 2.28 mg O2 L-1, and maximum of 33.85°C. This study reports the effect of such extreme environmental conditions on hard coral recruitment and CCA originally settled at a forereef on artificial substrates that were transplanted over two years in two fringing reef and at the Bouraké lagoon. Our data emphasize the negative effects of the extreme conditions in our study sites on the CCA, which decreased in cover by ca. 80% and lost in the competition with turf algae, which, in turn, increased up to 162% at the end of the two years. Conversely, hard coral recruitment remained high at Bouraké throughout the study, three-fold higher than at two sites located outside Bouraké where environmental conditions were typical for coastal fringing reefs. Our findings show that while such extreme, climate change like-conditions have a direct and adverse effect on CCA abundance, and despite a certain persistence, coral larvae settlement was not affected. Based on previous findings from Bouraké, and the present observations, both coral recruits and adults seem to be unaffected despite the extreme environmental conditions. This study supports previous research illustrating how extreme natural and variable environments may reveal unexpected and positive insights on the processes underlying coral acclimatization and adaptation to global change.


INTRODUCTION
Scleractinian corals begin their life as pelagic larvae eventually settling to the reef, growing and establishing as adult, longlasting colonies (Harrison and Wallace, 1990). Healthy and functional coral reef communities directly depend on coral recruitment success that, in turn, depends on the abundance of their preferred settling substrate i.e., crustose coralline algae (CCA, family Corallinaceae, Lithop hyllaceae, and Porolithaceae) (Heyward and Negri, 1999). Indeed, CCA are known as framework binders and act as a preferential settlement substrate for coral larvae (Morse et al., 1988;Heyward and Negri, 1999;Harrington et al., 2004;Nelson, 2009;Price, 2010). Whilst, CCA thallus surface topography differences imply biotic interaction changes (Figueiredo et al., 1996), this preferred substrate is most likely induced by chemical cues from the CCA wall-associated compounds and to a lesser extent from their bacterial biofilm (Tebben et al., 2015;Goḿez-Lemos et al., 2018). The early-life stage is a harsh period for coral survival since a multitude of biotic and environmental factors such as light (Bessell-Browne et al., 2017), smothering from sedimentation (Fabricius and Wolanski, 2000), and several complex trophic interactions (Fabricius et al., 2015) including competition with turf algae (Swierts and Vermeij, 2016;Jorissen et al., 2020), suspension-feeders (Elmer et al., 2018), grazers and predators (Penin et al., 2011;O'Leary et al., 2013), cause high mortality rates (67-99%) in coral recruits (Wilson and Harrison, 2005). In particular, competitive interactions play a significant role in the survival of both coral recruits and adult colonies (Penin et al., 2011;Swierts and Vermeij, 2016). In this context, the complex and only partially studied interactions between CCA and coral larvae and their consequences on larval settlement, have been shown to play a fundamental role in the replenishment of the reef (O'Leary et al., 2012). Likely, this complex ecological framework could be furthermore compromised in the future when considering climate changes, in particular rising temperatures and ocean acidification as it has been widely described across several studies (e.g., Mccoy and Kamenos, 2015;Cornwall et al., 2021a;Page et al., 2021).
Ocean acidification (OA) is generally described as the rapid and global decline of pH and calcium carbonate saturation in surface seawater (Orr et al., 2005). Several lab-based experiments addressing the effects of OA on CCA have demonstrated reduced calcification and growth rates Comeau et al., 2013), increased dissolution (Anthony et al., 2008;Diaz-Pulido et al., 2012), and reduced recruitment Bradassi et al., 2013). Field studies have shown as well reduced recruitment (Fabricius et al., 2011;Fabricius et al., 2015) or impaired capacity to attract coral larval settlement (Doropoulos et al., 2012). However, it has been recently shown that lab-based versus field experiments tend to conclude that OA has a stronger impact on CCA (Page et al., 2022). Indeed, CCA and corals persist at CO 2 seep sites (Kamenos et al., 2016;Enochs et al., 2020; but see also Fabricius et al., 2011 for an opposite finding) and at low pH environments in Palau (Barkley et al., 2015). For coral recruits, OA has been shown to cause skeletal malformations (Foster et al., 2016), and decrease in availability of settlement cues (Albright et al., 2010;Albright and Langdon, 2011). Thus, in general, OA could affect coral recruitment by dissolving CCA or changing their microbial biofilms compounds, hence decreasing the availability of a supposedly ideal substrate for coral recruits (Webster et al., 2013) and finally affecting the recruit survival (Doropoulos et al., 2012;Doropoulos and Diaz-Pulido, 2013). However, OA is only one of the ongoing threats to coral reefs, and its combination with both ocean warming (OW) and deoxygenation (OD) will likely cause CCA bleaching and mortality Cornwall et al., 2019;Hughes et al., 2020), and impairment of coral larvae settlement (Jorissen and Nugues, 2021). Combining all those environmental parameters and testing their effects using short-term lab-based experiments is challenging. Natural systems that are exposed to natural in-situ gradients can be used to study their integrated effect (Page et al., 2022). Despite several limitations, few study sites have been recognized as natural laboratories to study the response of species and ecosystem levels to future extreme environmental conditions (e.g., Hall-Spencer et al., 2008;Fabricius et al., 2011;Kroeker et al., 2013;Barkley et al., 2015;Agostini et al., 2018). The semienclosed lagoon of Bourakéin New Caledonia is one of these natural laboratories; it hosts well-diversified and complex coral assemblages despite fluctuating acidified, warm, and deoxygenated environments (Camp et al., 2017;Maggioni et al., 2021). In this study, the Bourakésite provided the opportunity to test the long-term responses of CCA and coral recruitment to multiple stressors in a natural setting, which would be hard to reproduce in tank experiments artificially mimicking future climatic conditions.
Here, we assessed the effects of extreme and fluctuating environmental conditions (i.e., low pH, O 2 , and high temperature) on the abundance of CCA, algal turf, and coral recruits. To address this objective, we obtained CCA-covered tiles from a forereef and moved them to three fringing coastal reef sites with distinct environmental conditions. We investigated the effect of the environmental conditions typical of coastal fringing reefs on the CCA cover twice during the two years experiment and we quantified and identified the coral family recruited on the tiles at the end. Based on the reported effect of OA on CCA cover and coral recruitment rates, we hypothesized that in Bourakéboth will be heavily affected by the extreme conditions. that is enclosed between two islets and has an abundant and scattered reef. Station M1 is a typical fringing reef that is open to the lagoon.
In Bourakéextreme conditions are recorded in seawater temperature (range 17.5 -33.8°C), dissolved oxygen (range 1.87 -7.24 mg L -1 ), and pH (range 7.23 -7.92 pH T units) (Maggioni et al., 2021). Those environmental parameters fluctuate regularly according to the tide (ca. 1 m high) and only return close to coastal values during the high tide of spring tides when seawater of the semi-enclosed lagoon is partially renewed with seawater entering from the external lagoon. As a result, daily fluctuations of dissolved oxygen (DO), temperature, and pH fluctuate can reach 4.91 mg O 2 L -1 , 6.50°C, and 0.69 pH T units, respectively. Despite these harsh conditions, at least 66 coral species form compact and wealthy reefs all along the lagoon (Maggioni et al., 2021). In the first assessment of the Bourakélagoon (Camp et al., 2017), the authors found that these coral species spent 44% of the time at pH T 7.7-7.8, thus reaching the worst scenarios RCP7.0 and RCP8.5 predicted in global ocean by the Inter Panel on Climate Change for the end of the century (IPCC Report, 2021). During particularly warm years, such as the El Niño event in 2016 (Camp et al., 2017), they experienced temperatures predicted for the end of this century 71% of the time, while facing now the best case scenarios RCP2.6 and RCP1.9 (IPCC Report, 2021).
The composition of the reef communities was assessed at each station visually, by taking photographs of the corals found along three 30-m line transects, and when necessary by collecting fragments of corals and further observing their skeletal structure under a stereomicroscope for final identification. Most of the corals, which are common fringing reef species, were recognized at the genus and species level, with reference to a species list previously established for B3 and R2 (Maggioni et al., 2021). Because the reef community composition might affect coral recruitment rates and the comparison between sites, we verified that the reef diversity was similar among the three stations using the Shannon's diversity index (H) on identified species trough the stations (Supplementary Tables S1, S4). The three stations were chosen because in addition to have similar reef composition, they allow to compare the effect of a gradient of seawater carbonate chemistry on the larvae and CCA settlement. In fact, as previously found (Maggioni et al., 2021) and also during our measurements (see results section), R2 could be considered as an intermediate station between B3 and M1, regarding pH absolute value. Therefore, this continuum in the pH experienced by the organisms allows us to evaluate the impact of three different pH scenarios: the end of the century, future 30 years, and the actual conditions (i.e., B3, R2 and M1, respectively).
For stations, B3 and R2, total scale seawater pH T , dissolved oxygen, and temperature were periodically undertaken since 2016 (Maggioni et al., 2021) using Seabird SeaFET pH loggers, YSI 600 OMS-M, and Hobo water temperature Pro V2, all settled at 10-min logging intervals. For station M1, seawater pH T and dissolved oxygen content (DO, mg L -1 ) were measured during three 7-14 day periods (August 2017, June 2018, and January 2019). Unfortunately, the temperature probe was lost twice, therefore the temperature variations could not be assessed for this site.

Terracotta Tiles Set-Up and Settlement Analyses
A total of 60 red terracotta tiles (11 x 11 x 1 cm) were deployed after being soaked for 7 months at a forereef (22°17.152'S, 166°1 0.992'E) at 12 meters depth. In September 2017, they were transferred and randomly deployed at the three stations B3, R2, and M1. This experimental approach was necessary to obtain 60 tiles having both the same CCA community and a similar percentage cover that could be difficult to obtain if the tiles were pre-soaked directly at the stations. This same basis cover allowed us to consider in particular the effects of environmental changes. At each station, 20 tiles were secured to the reef with metal rods at 1.5-2 m depth. Tiles were arranged with a 30degree inclination at approximately 10 cm from the substrate. They were positioned ca. 3 m distant along a transect of ca. 60 m. The tiles were therefore photographed underwater before moving them to the three stations (T 0 ), two months after their deployment (T 2 ), and two years after deployment (T 24/26 ). The top sides of 39 (off 60) randomly chosen tiles were photographed in situ at T 0 , using an underwater Canon G16 camera. The top sides of all the tiles were imaged at T 2 . Tiles from B3 and R2 were retrieved in September 2019 (T 24 ). Tiles from M1 were collected in October 2019 (T 26 ). Of the 20 tiles initially installed at St B3 and M1, only 15 and 19 were found, respectively. All the deployed tiles were recovered at R2. In the laboratory, both top and bottom sides were photographed under a dissecting stereomicroscope (tiles were immersed in seawater).
All photographs were analyzed using the photoQuad software (version 1.4) as in Trygonis and Sini (2012), and the percent cover of sessile organisms was assessed. Each photograph was calibrated for size and pixels-segmented to the highest level possible (segments ignored < 10 pixels). PhotoQuad's image segmentation is based on the SRM algorithm developed by Nock and Nielsen (2004). The photograph's size was usually comprised between 3x10 6 pixels to 9x10 6 pixels for a surface area typically around 125 cm 2 . For top sides, three main categories describing the settled organisms were identified: (i) Algal turf, (ii) CCA, and (iii) empty (i.e., tile surface devoid of any recognizable benthic organisms). For bottom tile sides, one more category was assessed, namely: (iv) suspension-feeders (i.e., ascidians, sponges, serpulid worms, and bryozoans). At the end of the experiment, tiles were cleaned for 24 h in a 10% NaCl solution to remove all organic material, rinsed, and dried. The coral recruits settled on the top, bottom, and the four lateral sides (all sides pooled) of the tiles were counted under a stereomicroscope and identified to the lowest possible taxonomic level. We were able to categorize the coral recruits into three families: Acroporidae, Poritidae, and Pocilloporidae.

Statistical Analyses and Data Presentation
The database for the map (Figure 1) was collected from map tiles at www.georep.nc ( © Georep contributors) and customized in QGIS (version 3.4.14). Statistical analyses were conducted and figures were produced using RStudio (R Development Core Team, version 4.1.0, 2021), including the packages: "ggplot2", "stats", "dunn.test", "FSA". Homogeneity of variance was tested using the Levene test and the normality of variance distribution was tested using the Shapiro-Wilk test and checked graphically with a Q-Q plot. When data did not conform to normality or homogeneity of variance, we used the non-parametric Wilcoxon test or the Kruskal-Wallis followed by the Dunn's multiple comparisons test (Bonferroni-adjusted). First, for each site separately we used two factors Wilcoxon test (two groups) for each of the three categories of organisms (CCA, Turf, and Empty) settled on the top sides of tiles to test for statistical differences first after two months of incubation (i.e., T 0 vs T 2 ), and then after two years (i.e., T 2 vs T 24/26 ). To compare significant differences between stations (i.e., B3, R2, and M1) at the start of the study (T 2 ), we used the non-parametric Kruskal-Wallis followed by the Dunn's multiple comparisons test (Bonferroni-adjusted). The latter was also used for each of the four categories identified on the bottom sides of tiles (i.e., Turf, CCA, Suspension-feeders, and Empty) and the three categories identified on the top sides of tiles (i.e., Turf, CCA and Empty) at the end of the experiment (T 24/26 ) to test for statistical differences between stations (B3, R2, M1). The same tests were also used for statistical differences in coral recruits between stations (B3, R2, and M1) from all tiles sides pooled at T 24/26 . Data were graphically reported as median values ± 25 th and 75 th percentiles (box), minimum and maximum values (whiskers), otherwise specified.

Environmental Parameters
Our results showed differences in the main environmental parameters between stations with lower values at B3, intermediates at R2, and normal at M1 (Table 1). At B3, pH was strongly correlated with the tidal cycle decreasing and increasing with the ebb and rising tide, respectively ( Figure 2). Within a tidal cycle, the pH reached a minimum of 7.23 and a maximum of 8.06 pH T units at B3 (Table 1) and it changed on average by about 0.15 units per hour. At stations R2 and M1, pH fluctuated between 7.69 and 8.12 pH T units for R2, and from 7.93 to 8.18 pH T units for M1. Similar to pH values, the mean dissolved oxygen concentration (DO) was lower in the Bourakélagoon than at the other stations. DO concentrations were lower in the early morning at R2 and M1 and higher at mid-day, while they were tidally driven at B3 with low and high values coincident with low and high tide, respectively. Over a 48-hour cycle, DO fluctuates up to 4.25 mg L -1 at station R2 and 3.75 mg L -1 at stations M1 and B3.

Changes in Crustose Coralline Algae and Turf Covers
After 7 months of conditioning at a forereef, i.e., just before they were moved and randomly shared between the three study stations (T 0 ), the top sides of the tiles were predominantly covered in CCA (median percent cover of 69.21%, all tiles pooled), while algal turf was rare (median percent cover of 11.95%). Two main CCA genera were recognized on the tiles on a simple visual check regarding the conceptacles size and disposal: Porolithon and Neogoniolithon. After two months (T 2 ) at the three stations, the CCA median percent cover was significantly different among stations and significantly decreased at B3 from 69.21% to 41.4% (Figure 3; Table 2). CCA cover significantly increased at the two stations R2 and M1, respectively 64.4% and 81.3%. In contrast, turf cover was significantly different between stations, reaching the highest values in B3 (20.4%), and only 6.4% and 9.7% at R2 and M1, respectively. After two years of deployment (T 24/26 ), CCA median percent cover on the top sides of the tiles was lower than the one at T 2 at all stations: this decrease was dramatic at B3 (from 41.4% to 8.4%) and M1 (from 81.3% to 1.5%), but to a much lower extent at R2 (from 64.4% to 49.3%).
Turf cover on the top significantly increased at all stations, dramatically at both B3 and M1 (from 20.4% to 53.5%, and from 9.7% to 91.2%, respectively), to a much lower extent at R2 (from 6.4% to 18.7%). Similarly, Turf cover on the bottom was significantly different between stations (Table S3; Figure S3), low at B3 (31.89%) and high at R2 and M1 (71.79% and 71.56%, respectively). Crustose coralline algae median percent cover on the bottom was similar at B3 (16.12%) and R2 (21.21%), and significantly lower at M1 (3.42%) after two years of immersion Suspension-feeder percent covers on the bottom of tiles in B3 and R2 were as low as 2.74% and 9.57% for B3 and M1, respectively, and rare for R2 (Table S3). Tile surface areas devoid of settlers significantly differed between stations. They were similar at M1 and R2 (6.78% and 6.26%, respectively) and higher at B3 (39.47%).

Spatial Variability in Coral Recruitment
Coral recruitment significantly differed among stations and sides after two years of tile immersion (Table S3). It was found to be ca. 3-fold higher at B3 than R2 and M1 with 18.4, 4.6, and 5.5 recruits per tile (all tiles pooled per station), respectively. Regarding where the larvae preferentially settled on the tile, the number of recruits was always higher on the bottom than on the top, and scarce on the lateral sides ( Figure 4).
Acroporidae was the dominant family of coral recruits settled in our experiment ( Figure S4), with similar relative abundances between stations (65.58%, 73.91%, and 67.33%, at B3, R2, and M1, respectively). Pocilloporidae was the second most abundant Mean (SD), minimum (min) and maximum (max) values of temperature, pH T (total scale), and dissolved oxygen (DO). Values were calculated from data collected during August 2017, June 2018 and January 2019 at station M1, and during 3-year long-term monitoring at B3 and R2 considering both diel and seasonal fluctuations, and compiling several thousands of measurements (Camp et al., 2017;Maggioni et al., 2021). The temperature at station M1 was only occasionally recorded during the 2-year study period, therefore not presented.
family, with recruits who equally settled at all three stations in terms of relative abundance (14.85%, 17.39%, and 16.83% at B3, R2, and M1, respectively). Only a few Poritidae were found at B3 and M1 (17.39%, and 6.52%), and none at M1 which exhibit more unknown genera.

DISCUSSION
Coral larvae are tied to preferred settlement substrates such as CCA (Harrington et al., 2004;Price, 2010), which in turn can be severely affected by environmental conditions, such as OA. In fact, reducing both the availability of preferred substrate for coral settlement, such as CCA and biofilms and increasing turf algae cover, may negatively reduce recruitment rates and coral population replenishment and/or recovery (Doropoulos et al., 2012;O'Leary et al., 2012;Fabricius et al., 2017). Although the CCA we transplanted at study sites suffered from moderate bleaching, partially confounding the effect of environmental conditions on their survival, the present study corroborates previous findings reporting a decline in CCA cover and a shift towards an increase in turf algae under extreme conditions Martin et al., 2008;Barott et al., 2012;Kroeker et al., 2013;Porzio et al., 2013;Enochs et al., 2015). While the CCA abundance decline was expected, surprisingly, coral recruitment rates in Bourakéwere higher than outside, at both the two fringing reefs. Our study sites differed in several ways, one is tidally driven and exhibits variable and extreme pH, temperature and hypoxia while the two others are far less variable and close to ambient values. While our findings appear to contradict results from several OA laboratory-based experiments (e.g., Doropoulos et al., 2012;Webster et al., 2013;Jorissen et al., 2020), they are in accordance with some field experiments (Kamenos et al., 2016;Enochs et al., 2020) showing  Table 2). in general, a lesser impact of OA (Page et al., 2022). All these observations together reinforce the idea that coral recruitment is possible even under extreme conditions and that, particularly for the Bourakésite, corals may have adapted to such variable and extreme environmental conditions.

Shift Between CCA and Turf Algae
CCA cover on the top sides of the tiles remarkably declined at Bourakéafter only 2 months of immersion, while remaining stable or slightly higher at the two other reef stations. This suggests that ca. 40% of CCA settled during the tiles predeployment, dissolved during a short-term exposure, therefore further indicating a strong negative effect of the Bourakénatural environmental conditions on CCA. At the same time, turf algae cover was significantly higher at Bourakéthan at other reef stations. It could be argued that our findings were mostly due to the effect of an abrupt change of environmental conditions, including light exposure, from those characterizing the forereef where the tiles were soaked before deployment to the three study sites. The abrupt deployment could have impacted the benthic organisms present on tiles and indeed, moderate bleaching was evident after the first two months at the study sites (see Figure 5A-C), potentially accelerating a decline in the CCA abundance. However, such deployment effect, should be common across the three stations, and this has not been observed. Indeed, while CCA cover decreased in Bourake, it was stable or slightly increased at the two stations R2 and M1, after two months ( Table 2), which we believe to be a long enough time for any consequence of an abrupt displacement to occur. Furthermore, the initial phase-shift between CCA and turf seen after the first two months at Bourake, was still confirmed after our two-years in situ experiment: a significant loss in CCA cover (ca. 80%) and a significant increase of turf algae (up to 162%) were measured ( Figure 5D-F). While we cannot fully dismiss (A and B) Non-parametric Wilcoxon tests between the beginning (T 0 ) and after two-month incubation (T 2 ); and between the T 2 and 2-year (T 24/26 ) at the three stations. C) Non-parametric Kruskal-Wallis test, followed by the Dunn's multiple comparisons, between stations for CCA, Turf and Empty measured at the beginning (T 2 ) and the end (T 24/26 ). Significant values are in bold (p<0.05).
the hypothesis that the initial deployment has affected the tiles, our 2-years observation did not show a substantial recovery in the color of the settled CCA and we are unable to discern the effect of bleaching from that environmental conditions. A potential caveat of our experimental design is that forereef where the CCA community (mostly Porolithon and Neogoniolithon) had settled on the tiles was characterized by different environmental condition compared to the three study sites where the tiles were kept during the 2-year incubation. Although both CCA families are also commonly found on coastal fringing reefs (Claude Payri, personal communication), and the same families were also found at our study sites, we are unable to verify whether the CCAs we transplanted were the same as those that would be found naturally at the study sites.
Inter-species interactions, in combination with a direct effect of acidification and the extreme environmental conditions in the Bourakésystem on CCA skeletons, could partially explain this cover drop-off despite a certain persistence obtained by altering skeletal mineralogy, although it remains difficult to assess with certainty (Goffredo et al., 2014;Kamenos et al., 2016). In Bourake, seawater pH decreased down to 7.23 and exhibited a mean value of 7.67, thus reducing carbonate ion concentrations ([CO 2− 3 ]) available and the saturation state of seawater with respect to aragonite and calcite (W arag ). Low values of seawater pH have been experimentally shown to cause a decrease in CCA growth and an increase in skeletal dissolution Comeau et al., 2013;Fabricius et al., 2015;Foster et al., 2016), and are known to reduce their skeletal density (Chan et al., 2020;Williams et al., 2021). Our findings are in agreement with several reports of overall CCA cover reduction ranging from 43% to 92% under low pH Kuffner et al., 2008;Doropoulos et al., 2012;Fabricius et al., 2015), and from 60% to 100% in combination with high temperature (Martin and Gattuso, 2009;Diaz-Pulido et al., 2012;Vogel et al., 2016). Dissolved oxygen concentration in Bourakéexhibits a mean value of 5.23 mg L -1 with a minimum of 2.28 mg L -1 ( Table 1) which is equivalent to depletion by ca. 30% compared to both stations R2 and M1. Despite the lack of research on hypoxia effects, we could hypothesize an additional negative effect from low dissolved oxygen on CCA growth as it has been previously reported (Altieri et al., 2017;Nelson and Altieri, 2019).
The decline of CCA cover could also be due to the combined effect of decreased light exposure and high sedimentation. Indeed, under OA conditions, CCA growth significantly decreases under the combination of low light and low flow conditions , and faces dissolution in shade environments (Martin et al., 2013;Bessell-Browne et al., 2017;McNicholl and Koch, 2021). Although light irradiance was not measured and sedimentation rates in the Bourakélagoon have yet to be comprehensively measured, it has been suggested that high turbidity in Bourakécould attenuate solar radiation (Jacquemont et al., 2022). Indeed, the Bourakélagoon is inside a mangrove forest, and when the sediment interstitial water is sucked by the current that empties the mangrove and arrives in the lagoon during ebb tide (spring ebb tide in particular), turbidity temporally increases because of the organic compounds released from the mangrove. This, in turn, and either alone or combined with the high pCO 2 and the warm temperatures in Bourakémay enhance turf algae productivity by alleviating carbon limitation (Johnson et al., 2017). Although such potential negative mechanisms, in Bourakémore than 66 species of corals, and also some CCA on coral rubbles were found (Maggioni et al., 2021). Further measurements of light availability to CCA along diel and seasonal cycles, as well as specific experiments will be necessary to disentangle the contribution of turbidity in combination with high pCO 2 on the CCA growth and dissolution. However, CCA and turf algae cover on the bottom sides of the tiles i.e., the CCA cover median was low FIGURE 4 | Number of coral recruits (per tile) found on the different sides of the tiles immersed during ca. two years at Bouraké(B3) and stations R2 and M1. Data are means ± SD (number of tiles n = 15, 20, and 19, for B3, R2, and M1, respectively). Data from the lateral sides were pooled. Stars represent statistical significance (see Table S3).
in B3 than in R2, and quasi absent in M1, partially confirmed data found on the top sides ( Figure 5G-I). While the quasi absence of CCA at M1 could be due to both the high turf algae cover (ca. 75%) and the presence of suspension feeders such as bryozoans, known to be efficient competitors for space (Glassom et al., 2004;Elmer et al., 2018), the presence of CCA at both B3 and R2 at bottom side support the hypothesis that the decline of CCA found at the top side is due to the combination of the high pCO 2 and other factors such as light and sedimentation.
Our results emphasized the extent to which CCA growth is limited at Bourakéwhen compared to the station R2 where CCA cover decreased by 24% during the two years of deployment. This result was a bit unexpected because although previous measurements suggested that R2 could show averaged pH values lower than for a typical reef (Maggioni et al., 2021), the study stations were selected having in mind that R2 and M1 were two reef stations having the main environmental parameters (e.g., seawater pH) close to ambient values typical for a coastal fringing reef. However, in terms of seawater pH, R2 turned out to be intermediate between Bourakéand M1 according to recent data on the same site (Maggioni et al., 2021). With regard to R2, the low seawater pH (mean of 7.89 and a minimum of 7.69 pH T units) might explain why we measured some decrease in CCA cover, since similar pH values have led experimentally to a significant loss in CCA cover Doropoulos et al., 2012;Comeau et al., 2013;Fabricius et al., 2015). In fact, this presupposed reference station actually exhibited a pH close to values expected in global ocean from the RCP7.0 scenario in 30-40 years (IPCC Report, 2021). These unexpected environmental conditions widened our opportunities to study the effect of OA on CCA and coral recruitment. The low pH was likely due to the site topography. Indeed, R2 is located in a "culde-sac" between two islets with sparse but abundant reefs. It is characterized by a fine sediment bottom, and a low seawater renewal rate, likely such as for the Rock Islands in Palau (Micronesia) (Barkley et al., 2015;Shamberger et al., 2018;Kurihara et al., 2021). The combined effect of biological  activities and carbonate disequilibrium likely caused such low pH values which are frequent and typical for many coastal fringing reefs with limited seawater exchanges with the lagoon. In contrast to Bourake, turbid waters were never observed at R2 that benefits from greater light intensity; despite low pH, environmental conditions would still allow CCA persistence. At station M1, CCA was completely dissolved after two years and tiles were recovered by turf algae. While the effect of pH, in this case, can be discarded (ambient pH of ca. 8.07 pH T units), the combined effect of lower grazing pressure and wellestablished wave exposure may partially explain this evolution. It has been shown that grazer presence depends directly on factors such as overfishing, tidal height and wave exposure (Doropoulos et al., 2017). In particular, turf algae seems more abundant where wave action is greater and grazers are less present. Conversely, turf algae are less abundant in sheltered areas where grazers are more active (Underwood and Jernakoff, 1984;Doropoulos et al., 2017). While Bourakéand R2 are sheltered from wave action and quite isolated sites, M1 is an open bay much more subject to wave action and a more intense fishing pressure due to its proximity to Noumeá. Our findings are in agreement with previous studies reporting both limited CCA growing on the top of the tiles (O'Leary et al., 2012;O'Leary et al., 2013) and an increase in the turf algal cover (Penin et al., 2011). In addition, we observed a high sedimentation rate in M1 with fine particles being trapped by the dense algal cover. This accumulation of sediment certainly contributed to smothering the underlying CCA.
Coral Recruitment Seems to be Maintained in BourakeÓ cean acidification, warming and deoxygenation have been shown to negatively affect coral larvae settlement by reducing the growth of favorite substrate (e.g., Diaz-Pulido et al., 2012;Doropoulos et al., 2012) and altering chemical biofilm and microbiological cues for settlement onto CCA (Albright et al., 2010;Albright and Langdon, 2011;Webster et al., 2013;Fabricius et al., 2017;. Our findings seem to provide an alternative scenario using a natural extreme environment where coral recruits survived. In fact, although the extreme conditions reduced the CCA cover in Bourake, after two years of exposure we observed a 3-fold increase in coral recruits when compared with the stations R2 and M1. Recruitment rates recorded in this study are difficult to compare with the literature because values are quite variable. For instance, in French Polynesia, it varies from 0 to 10.8 recruits tile -1 (Gleason, 1996;Adjeroud et al., 2007;Bramanti and Edmunds, 2016), and it was less in the Virgin Islands (average of 3 recruits tile -1 ; Edmunds, 2021) and in Palau (0.3 to 3.9 recruits tile -1 ; Victor, 2008), while no data is available for New Caledonia. There might be several explanations for our uncommon and counterintuitive findings. Firstly, because the whole set of tiles was maintained in situ for two years, potentially acted as substrates during two consecutive spawning events. The latter follows the full moon in November, coincident with the mass spawning period on the Great Barrier Reef in Australia (Baird et al., 2010). During the first spawning event, which happened 1-2 months after deployment, the CCA cover was still quite high (41.4%) and still may have allowed tiles to attract coral larvae. If coral recruitment intensity on tiles was the result of a single spawning event, it was quite massive and unusual. According to the disruptive effect of OA on the coral skeleton, especially during early life stages (Cohen et al., 2009;Foster et al., 2016), coral recruits would have been dissolved after two years or at least would have shown important malformation which was not the case (see below some considerations about the site M1). Most likely, our tiles received recruits also during the second spawning event, in spite of the low CCA cover. The role of biofilm acting as an ideal chemical cue for coral larvae settlement (Tebben et al., 2015) could be crucial in this case. Alternatively, coral spawning was either low during two consecutive years at R2 and M1 or recruits did not survive and were covered by competitors, such as bryozoans. The unexpected higher recruitment in B3 could also be due to high concentrations of organic carbon and high levels of nutrients such as [Si(OH) 4 ], [NO x ] and [PO 4 ] 3which were 8.88, 0.54 and 0.19 µmol L -1 , respectively, in Bouraké (Maggioni et al., 2021) and were, in general, more concentrated (up to 5 times) than at station R2. These concentrations coupled together would have increased coral heterotrophic rates and helped recruits to cope with the extreme conditions (Anthony, 2006;Houlbrèque and Ferrier-Pagès, 2009;Ragazzola et al., 2013;Goldberg, 2018). Increasing the energy budget through heterotrophic nutrition could enhance juvenile growth and facilitate coral recruitment under high temperatures (Huffmyer et al., 2021). This energy budget enhancement could also benefit coral larvae under deoxygenation which has recently been shown to increase the coral larvae energy needs at a dissolved oxygen concentration of around 2 mg L -1 (Alderdice et al., 2021).
The role of water circulation in the semi-enclosed lagoon of Bourakéhas also to be considered since a close water circulation could enhance coral larvae settlement via a decrease in larval dispersal (O'Connor et al., 2007). The particular morphology and hydrodynamics of Bourakécould make this lagoon a real recruitment hotspot as recently observed in St. John, US Virgin Islands (Edmunds, 2021). Recruitment was especially noticeable on the bottom sides of the tiles for each station, likely because tiles act as a shelter against grazing and predation and enhance coral recruitment compared with the top or lateral sides (Penin et al., 2011). Indeed, the top and lateral sides are prone to a more intense grazing activity. In addition, recruits face high levels of sedimentation on top sides which induces smothering (Fabricius and Wolanski, 2000), settlement limitation (Ricardo et al., 2017) and decreases in fertilization, embryogenesis and larval development (Jones et al., 2015). This response was similar among the three stations despite the different environmental conditions.
Our recruit identification revealed the same proportions of corals taxa settled between stations with the Acroporidae as the dominant family (70% of all recruits) followed by the Pocilloporidae (ca. 16%) and only a few Poritidae. This is in agreement with previous studies showing that Acroporidae and other broadcasting corals have high recruitment rates, especially in low watercirculation systems like Bouraké (Oprandi et al., 2019). The dominance of branching morphologies (Acroporidae and Pocilloporidae) may be due to their ability to compete with vertical CCA growth for space in light-limited environments (Jorissen et al., 2020). On the contrary, recruits from massive species such Poritidae, which mainly have a planar growth form, were rare, likely because they are unable to face competitive interactions with other coral species, CCA or suspension feeders. In agreement with this hypothesis, we found only a few recruits of Poritidae at R2 (6.52%) and zero at M1, the latter exhibiting a high density of bryozoans. At site M1 some recruit malformations were found and most of the recruits were classified as unknown genera. Such skeleton malformations are likely to be due to the competition with other organisms recovering the recruits, than to the effect of dissolving seawater carbonate values.
In this study, we used one of the few natural sites exhibiting extreme and highly variable environmental conditions including high temperature, low DO and low pH. Therefore, our experimental setting benefits from the natural conditions, but in a unique setting where CCA, adult corals and their recruits have been exposed to extreme conditions for several generations over likely hundreds of years. In particular, high pCO 2 and low dissolved oxygen in Bourakeẃ ere even worse than what is expected for 2100. Surprisingly, healthy and diversified coral populations still persist as likely because of unexpected high coral recruitment (this study) that gives a glimpse of hope for the future of coral reefs facing climate change. However, our findings need to be used with caution as the Bourakésystem is not a typical reef ecosystem. In fact, although our results showed that recruitment is possible even under extreme conditions, it does not mean that it might be the case for other reefs in the future. For instance, it has been shown that OA can cause morphological changes to the coral skeleton and thus, reduce their bulk density and wave resistance (Fantazzini et al., 2015;Tambutteé t al., 2015). Further measurements on the skeleton of coral recruits and adults (i.e., SEM observations) from Bourakéare underway to shed light on this possibility. Importantly, habitats displaying natural variability like Bourakéhave been hypothesized to be populated by resilient coral reef organisms (Padilla-Gamiño et al., 2016;Vargas et al., 2017;Kapsenberg and Cyronak, 2019;Enochs et al., 2020). Coral recruits could as well benefit from such high natural variability (Jiang et al., 2017;Jiang et al., 2020), although some studies described small influences from variability on some corals and CCA (Cornwall et al., 2018), increased negative effect of OA on CCA calcification (Johnson et al., 2019), or no systematic coral tolerance to warming (Klepac and Barshis, 2020). Defining how variability in multiple co-varying parameters could influence tolerance in biological responses of coral reef organisms remains a significant challenge for the scientific community (Kapsenberg et al., 2015;Boyd et al., 2016;Rivest et al., 2017). Even though the effect of environmental parameters has not been fully deciphered yet, we could assume that Bourake's corals could have acquired specific mechanisms providing recruits with the ability to face extreme environmental conditions (Putnam et al., 2016;Putnam et al., 2017). Which mechanism is used and to what extent acclimation vs adaptation is involved was not the objective of this study and would need further experiments using genomic approaches and reciprocal transplantations between Bourakeá nd reference stations. However, the present study provided evidence that survival and specific high coral recruitment are possible in natural and extreme habitats such as Bourakéleading to lifelong adapted coral reefs.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
RR-M conceived and designed the project. CT, RR-M and FB collected the data. FB identified the corals. CP counted coral recruits and identified them. CT performed the data analysis and drafted the manuscript in collaboration with RR-M. All coauthors read and edited the final version of the manuscript.

FUNDING
This study was supported by the French grant scheme Fonds Pacifique (no. 1976;project SuperCoraux). CT received a PhD fellowship from the University of Western Brittany (UBO) and the Brittany region (France). This study has also received financial support from the CNRS through the MITI interdisciplinary programs and was supported by ISblue project, Interdisciplinary graduate school for the blue planet (ANR-17-EURE-0015) and cofunded by a grant from the French government under the program "Investissements d'Avenir". Funding for FB was provided by KAUST baseline research funding and award FCC/1/1973-49-01.

ACKNOWLEDGMENTS
We wish to express our thanks to the captains of the IRD vessels and the diving technical staff for their fieldwork assistance. We thank Dr. Mehdi Adjeroud for providing us with the tiles used and for scientific advice on coral recruitment. We thank Dr. Claude Payri and Dr. Christopher Cornwall for their help and expertise regarding the assessment of organisms percentage cover and for CCA taxonomy and health review. We also thank the Province Sud of New-Caledonia for sample collection permits (no. 3413-2019). We are indebted to Federica Maggioni for environmental data sharing. We thank Yuri Artioli for his help throughout the reviewing and editing process of this article.