# Recovery of coral assemblages despite acute and recurrent disturbances on a South Central Pacific reef

Mehdi Adjeroud, Mohsen Kayal, Claudie Iborra-Cantonnet, Julie Vercelloni, Pauline Bosserelle, Vetea Liao, Yannick Chancerelle, Joachim Claudet, Lucie Penin

► To cite this version:

Mehdi Adjeroud, Mohsen Kayal, Claudie Iborra-Cantonnet, Julie Vercelloni, Pauline Bosserelle, et al. Recovery of coral assemblages despite acute and recurrent disturbances on a South Central Pacific reef. Scientific Reports, 2018, 8 (1), pp.9680. 10.1038/s41598-018-27891-3. hal-01945698

HAL Id: hal-01945698

https://univ-perp.hal.science/hal-01945698

Submitted on 5 Dec 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# Recovery of coral assemblages despite acute and recurrent disturbances on a South Central Pacific reef

Received: 1 March 2018

Accepted: 1 June 2018

Published online: 26 June 2018

Mehdi Adjeroud<fnref n="1" />,², Mohsen Kayal³, Claudie Iborra-Cantonnet<fnref n="1" />, Julie Vercelloni⁴,⁵, Pauline Bosserelle⁶, Vetea Liao³, Yannick Chancerelle³, Joachim Claudet²,⁷ & Lucie Penin²,⁸

Coral reefs are increasingly threatened by various types of disturbances, and their recovery is challenged by accelerating, human-induced environmental changes. Recurrent disturbances reduce the pool of mature adult colonies of reef-building corals and undermine post-disturbance recovery from newly settled recruits. Using a long-term interannual data set, we show that coral assemblages on the reef slope of Moorea, French Polynesia, have maintained a high capacity to recover despite a unique frequency of large-scale disturbances which, since the 1990s, have caused catastrophic declines in coral cover and abundance. In 2014, only four years after one of the most extreme cases of coral decline documented, abundance of juvenile and adult colonies had regained or exceeded pre-disturbance levels, and no phase-shift to macroalgal dominance was recorded. This rapid recovery has been achieved despite constantly low coral recruitment rates, suggesting a high post-disturbance survivorship of recruits. However, taxonomic differences in coral susceptibility to disturbances and contrasting recovery trajectories have resulted in changes in the relative composition of species. In the present context of global coral reef decline, our study establishes a new benchmark for the capacity of certain benthic reef communities to sustain and recover their coral cover from repeated, intense disturbances.

Coral reefs are threatened by various types of local and large-scale anthropogenic and natural disturbances that cause widespread mortalities of scleractinian corals, the primary framework builders and key components of reef health and biodiversity<fnref n="1" />⁻³. Many reefs have undergone striking phase shifts in their benthic communities, which classically involves the replacement of corals by fleshy macroalgae or other non reef-building organisms, an undesirable state which provides less ecosystem goods and services⁴. Given the likely continuation of such disturbances in the next decades, coral reefs are expected to be highly vulnerable to future environmental changes⁵⁻⁸. In this context, understanding the determinants of coral recovery and identifying disturbance regimes that cause shifts in species composition is crucial for the preservation of reef ecosystems.

While coral reefs from various regions have undergone severe and likely irreversible degradation⁹,¹⁰, some reefs have shown the capacity to recover from severe disturbances¹¹⁻¹³. Recovery of coral populations following disturbances relies on the arrival of newly settling larval recruits, as well as on the growth and propagation of surviving coral colonies¹³,¹⁴. Given the slow growth and long life span of most coral species, understanding mechanisms of coral recovery requires long-term surveys of population dynamics and demography¹²,¹⁵,¹⁶. However, due

<footnote n="1">¹Institut de Recherche pour le Développement (IRD), UMR 9220 ENTROPIE, UPVD 52 avenue Paul Alduy, 66860, Perpignan, France. ²Laboratoire d'Excellence CORAIL, Perpignan, France. ³USR 3278 CNRS-EPHE-UPVD, Centre de Recherches Insulaires et Observatoire de l'Environnement, UPVD 52 avenue Paul Alduy, 66860, Perpignan, France. ⁴Global Change Institute, The University of Queensland, Brisbane, QLD 4072, Queensland, Australia. ⁵ARC Centre of Excellence for Coral Reef Studies, The University of Queensland, Brisbane, QLD 4072, Queensland, Australia. ⁶Pacific Community (SPC), Fisheries, Aquaculture and Marine Ecosystem division, BP D5, 98848, Noumea, New Caledonia. ⁷National Center for Scientific Research, PSL Université Paris, CRIOBE, USR 3278 CNRS-EPHE-UPVD, Maison des Océans, 195 rue Saint-Jacques, 75005, Paris, France. ⁸Université de la Réunion, UMR 9220 ENTROPIE, 15 avenue René Cassin CS 92003, 97744, Saint Denis, Cédex 9, La Réunion, France. Correspondence and requests for materials should be addressed to M.A. (email: mehdi.adjeroud@ird.fr)</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location of the nine stations sampled on the outer reef slope around Moorea, representing three locations (Vaipahu, Tiahura, Haapiti), and three depths (6, 12, and 18 m). Stations codes are abbreviated as follows: the first letter indicates the location (V: Vaipahu; T: Tiahura; H: Haapiti), and the associated number represents the depth (6, 12, and 18 m). Dashed lines represent the approximate extent of the reef front surrounding Moorea. Note that the position of the stations is provided schematically and does not follow the spatial scale of the map. This map was created using Adobe Illustrator CS5 (http://www.adobe.com/fr/products/illustrator.html).</figcaption></figure>

to the difficulty in collecting such data over a time-period that allows for coral decline-recovery cycles to occur, few studies have examined the interplay between early recruitment, dynamics of juvenile and adult populations, and disturbances at a relevant time-scale¹⁵. Consequently, there is still a critical knowledge gap in the ecological drivers of coral community trajectory, particularly during recovery processes that follow disturbances¹⁶,¹⁷.

Here, we analyzed a unique long-term dataset to examine the interannual variability in cover (1991–2014; at a reference site) and abundance (2003–2014; across different locations and depths) of the three benthic stages of corals (recruits, juveniles, and adults)¹⁸ on the outer reef slope of Moorea Island (17°30′S, 149°50′W) in French Polynesia (Fig. 1). Since the 1980s, this reef has been affected by a remarkably high frequency and intensity of the three most destructive large-scale disturbances that affect Indo-Pacific reefs¹⁰,¹⁹: coral bleaching events, cyclones, and outbreaks of the coral-killing crown-of-thorns seastar *Acanthaster* spp. (COTS). Moorea is thus a unique study site to examine the temporal dynamics and recovery trajectories of coral populations affected by acute and recurrent disturbances, and the data set analyzed here is one of the most comprehensive in the Indo-Pacific.

## Results and Discussion

Prior to the establishment of our survey in 1991, Moorea's outer reefs had already recovered from two mass bleaching events (1984 and 1987), and a COTS outbreak (1980–1982) that reduced coral cover from ~45% in 1979 to ~12% in 1982²⁰. Between 1991 and 2014, the reefs were impacted by two cyclones, five bleaching episodes (associated with high sea-surface temperatures), and a catastrophic COTS outbreak that caused two additional episodes of major coral decline (Fig. 2a). The first major decline resulted from the combined impacts of a cyclone and a bleaching event that reduced live coral cover from ~51% in 1991 to ~22% in 1993. Coral cover increased rapidly after 1993 and recovered pre-disturbance values of 45–50% that persisted from 2001 to 2006, despite subsequent bleachings in 1994, 2002, and 2003 resulting from thermal anomalies that were similar in extent and prevalence to the 1991 event, but did not cause mass mortality²¹. The second major episode of coral decline resulted from the successive occurrence of a COTS outbreak in 2006–2010 followed by a cyclone in February 2010. This constituted one of the most extreme cases of coral reef decline documented, with live coral cover decreasing from ~49% in 2005 to <1% in 2010. Coral cover increased gradually thereafter, reaching ~17% in 2014 (Fig. 2a; Supplementary Table S1), and has since recovered to pre-disturbance levels (~50%) at several reef slope locations in 2016²².

Both major coral decline events were accompanied by a rapid increase in turf algae (from ~16% to ~50% cover; Fig. 2a; Supplementary Table S2), which is generally among the first to colonize vacant reef space⁴. However, in contrast with the 'typical' phase shifts that are increasingly reported on coral reef ecosystems⁴, a consecutive proliferation of fleshy macroalgae was not observed at Moorea (Fig. 2a; Supplementary Table S3). Resistance of Moorea's reefs to a transition to macroalgal dominance probably results from the high grazing
<figure><img src="image_3.png" /></figure>

**Figure 2.** Impact of eight large-scale disturbances between 1991 and 2014 on (a) cover of corals (all 18 genera), macroalgae and turf algae, and (b) cover of the four dominant coral genera at 12 m depth at Tiahura reef, Moorea island. Dotted lines represent standard deviation. A portion of the reef is shown through time (photos A–D). (A) Coral dominate the healthy reef (coral cover >45%). (B) Turf algae have colonized dead coral colonies following outbreaks of the coral-killing crown-of-thorns seastar *Acanthaster* spp. (COTS) (<5% coral cover). (C) Mostly dead and weakened coral colonies were swept away by the cyclone Oli in February 2010 (<1% coral cover). (D) Juvenile colonies of *Pocillopora* spp. have recolonized the substrate (~17% coral cover). © Photos Mohsen Kayal (A–C), Yannick Chancerelle (D).

pressure by herbivorous fish, most notably the parrotfishes (Scarinae), whose densities remain stable, or even increase, following disturbances at Moorea, and which are considered well above those needed to prevent proliferation of macroalgae<sup>23,24</sup>.

The response of coral populations to disturbances and their recovery trajectories differed among the four dominant genera (Fig. 2b), which represent >80% of the overall coral cover. *Acropora*, *Montipora* and *Pocillopora* cover dropped following the disturbances of 1991 (~ two-fold decrease), but between 1992 and 2005, *Acropora* showed a relatively continuous recovery trend, *Montipora* declined further as a result of the bleaching events of 1994 and 2002, and *Pocillopora* recovered partially within four years and remained relatively stable (13–17%) thereafter. In contrast, the cover of *Porites* was unaffected by the first five disturbances, and its cover increased gradually from ~4% to ~12% between 1991 and 2005 to become the second most abundant taxa. Since the second major coral decline (2006–2010), which affected all of the four genera, recovery has been mostly driven by *Pocillopora*, whose percent cover had increased rapidly from 2010 to 2014 (from ~0% to ~9%).

Similar to changes in percent live coral cover, a general decreasing trend in the abundance of juvenile and adult coral colonies was observed between 2003 and 2010 (~ ten-fold decrease for adults), quickly followed by a six-fold increase between 2010 and 2014 (Fig. 3a,b). This pattern was recorded at all locations and depths, even if the trajectories varied among coral taxa and life stages (see generalized linear mixed-effect model; Supplementary Tables S4 to S6). Adult coral trajectories varied with depth for all genera except *Acropora*, which differed among locations (see Supplementary Figs S1 to S5). Juvenile coral trajectories were consistent across depths and locations for *Pocillopora* and *Montipora* but varied with depths for *Porites* and with locations for *Acropora* (see Supplementary Figs S6 to S10). The post-disturbance increase in the overall juvenile and adult abundance was mainly driven by *Pocillopora* and, to a lesser extent, *Porites*. In contrast, no signs of recovery were found for *Acropora* juveniles and adults at any locations or depths. The temporal dynamics of coral recruits were different to those described for adults and juveniles (Fig. 3c). For recruits, no obvious trend was recorded during the monitoring period (2001–2014), except an overall decreasing trajectory of *Acroporidae*, and a slight increase of *Pocilloporidae* at 6 m depth at Vaipahu after the occurrence of cyclone Oli in 2010 (see Supplementary Figs S11 to S15). In fact, recruitment rates at Moorea were always low (<140 recruits.m⁻²) compared to other Pacific reefs characterized by a dominance of *Pocilloporidae* and a low contribution of *Acroporidae* recruits, similar to those of high-latitude or sub-tropical Indo-Pacific reefs¹⁷,²⁵.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Temporal dynamics in (a) adult, (b) juvenile, and (c) recruit coral abundance impacted by recent disturbances. All three locations and three depths sampled around Moorea were pooled. Dotted lines represent standard deviation.</figcaption></figure>

Moorea is among the few Indo-Pacific reefs that have shown the capacity to recover from severe and recur-
rent large-scale disturbances<sup>12,26</sup>. However, the post-disturbance dynamics of coral assemblages reported here
are marked by a number of striking differences compared to those dominating contemporary literature. Notably,
despite the frequency and severity of the disturbances that have affected Moorea in the last three decades, and the
two major declines in coral cover and abundance that rival some of the biggest losses reported elsewhere<sup>9</sup>, this
reef has shown, on two occasions, a remarkable capacity to quickly recover. Furthermore, we found successful
recovery at all habitats surveyed, across various depths and locations. Successful recovery is generally associated
with significant inputs of coral recruits<sup>27,28</sup>, but at Moorea, recovery was achieved despite constantly low recruit-
ment rates. Larval input from other less impacted islands may explain this successful recovery, as suggested by
the strong genetic connectivity among populations of *Pocillopora meandrina*, the dominant species on the outer
reef slope within the Society Archipelago<sup>29</sup>. Alternatively, coral populations in habitats less impacted by distur-
bances, such as deeper reef slopes (≥30 m depth) and lagoons, may produce enough larvae to support recoloni-
zation through self-seeding<sup>30,31</sup>. Moreover, the recruits that colonized the reefs after the 2006–2010 decline have
probably benefited from a higher availability of suitable substrate, with a low level of macroalgal cover due to the
high grazing pressure by herbivorous fish, and a lower post-settlement mortality due to reduced predation and
competition with other benthic organisms<sup>18,19,32</sup>.

Our results highlight the contrasted susceptibility among coral taxa to disturbances and recovery trajectories
that led to significant changes in the structure of the assemblages between 1991 and 2014. This time-interval was
characterized by short periods of intense reorganization in response to the disturbances of 1991 and 2006–2007,
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Changes in the structure of the coral assemblage from 1991 to 2014 at 12 m depth at Tiahura reef. Non-metric multidimensional scaling (nMDS) using the Bray-Curtis dissimilarity index of the cover of coral genera recorded annually. The scale of the circle for each year is proportional to the overall coral cover. Vectors link points by chronology.</figcaption></figure>

followed by longer periods of slower changes (Fig. 4). This outcome is in accordance with the process of 'recovery without resilience'²⁰, highlighting that successful recovery in terms of overall coral cover and abundance is not necessarily accompanied by a reassembly of the relative abundances of component taxa, thus preventing a 'true' and complete recovery. These discrepancies in the capacity of coral taxa for population maintenance and recovery are most likely linked to differences in life history traits such as reproduction strategies, stress tolerance, growth capacities, and competitive abilities³³,³⁴. Indeed, *Pocillopora* clearly exhibits an opportunistic life strategy in French Polynesia, with a high population turnover resulting from comparatively higher recruitment rates, a high potential for larval dispersal and colonization, fast growth, and low resistance to perturbations¹⁸,³⁴. On the other hand, massive *Porites* shows a resistant life strategy, with slow turnover among individuals resulting from the production of fewer off-spring endowed with a high capacity for survival, slow growth to large colony sizes, as well as lower rates of decline in presence of harsh conditions and catastrophic events¹⁸,³⁴. *Acropora* and *Montipora* are competitive species characterized by high aptitudes to preempt space and resources in optimal environments through faster individual growth and by a high susceptibility to disturbances¹⁹,³⁴. A high potential for recovery following physical disturbance, such as wave damage during cyclones, has also been observed for *Acropora*¹¹, which has been related to the capacity of this taxon for asexual propagation through high survival and reattachment of colony fragments³⁴. However, recurrent disturbances have probably overwhelmed the capacity of competitive coral taxa for recovery, and have progressively pushed communities towards the predominance of opportunistic and stress-tolerant species. At Moorea, *Pocillopora* and *Porites* are thus clearly the contemporary and, most probably, future ecological 'winners', whereas *Acropora* and *Montipora* appear to be the 'losers'³²,³⁴.

The outcomes of this study support the hypothesis that, regardless of their ultimate fate, some reefs will undergo gradual changes in the structure of their coral communities in response to major stress, rather than collapse abruptly and irreversibly (see the concept of 'novel coral reef ecosystems'³⁵). However, successful recovery is conditioned by a low frequency of large-scale disturbances that should allow juvenile corals to reach maturity and contribute to the larval pool, but this scenario is challenged by some recent projections⁶,⁸. Another critical point now is to determine if changing coral assemblages will maintain the ecological and socio-economic goods and services that they have provided until recently⁵.

## Methods

**Study site.** This study was conducted at Moorea (17°30′S, 149°50′W), located 17 km to the northwest of Tahiti in the Society Archipelago, French Polynesia. Approximately 17000 people live on Moorea. It comprises 134 km² of land, 49 km² of reefs and lagoon, and has a circumference of 61 km (Fig. 1). The island is surrounded by a narrow coral reef of 2 km maximum width, which can be divided into three major reef habitats: fringing reef, barrier reef flat (separated from the fringing reef by a narrow sandy channel), and outer reef slope (separated from the barrier reef flat by the reef front). The reef slope is largely free of direct anthropogenic disturbances and nutrient loads are low³⁶,³⁷. The tides are semi-diurnal with an amplitude rarely exceeding 40 cm. Surveys of coral abundance were undertaken on the outer reef slope in three water depths (6, 12, and 18 m) at each of three locations (Vaipahu, Tiahura, Haapiti), giving a total of nine stations (Fig. 1). Vaipahu and Tiahura are situated on the northern shore subjected to the north swell that occurs exclusively during the austral summer (November to April), whereas Haapiti is located on the western shore and is directly exposed to the prevailing southwest swell, which is of higher amplitude and occurs throughout the year.
**Sampling design.** The percent cover of scleractinian corals (at genus level; 18 genera were recorded during the study period; see Supplementary Table S7), turf algae and macroalgae was recorded at 12 m depth at Tiahura, along four permanent transects of 25 m length, oriented parallel to each other and to the reef front (this station was established on a neighboring spur adjacent to the one selected for coral abundance surveys – see below). For the purpose of the present analyses, each transect was treated as a statistical replicate in all temporal contrasts. We used the Point Intercept Transect Method, with points placed every 0.25 m, to estimate cover of corals and algae. Data were collected once a year from 1991 to 2014, in March-April¹¹. The significance of the interannual variability of percent cover was evaluated using the nonparametric test of Friedman, because of the absence of normality in the data set, and because these data were not independent among years. When a significant inter-annual variability was detected, the Wilcoxon test was used a *posteriori* to compare values of two different years.

The non-reproductive benthic phase of corals can be separated into a recruit stage, i.e. newly settled corals invisible to the naked eye on natural substrata, generally aged less than one year, and a juvenile stage, observable in situ (>1 cm) and typically aged at least one year¹⁶,¹⁸. In the present study, recruits were less than three months old; juveniles were defined as colonies with maximum diameters between 1 and 5 cm, whereas colonies >5 cm were regarded as adults¹⁶,¹⁸. Juvenile and adult corals were sampled along 3 replicate permanent 10 × 1 m belt-transects at each station, once a year between 2003 and 2011 and in 2014, at each of the nine stations. Transects were parallel to the depth contour and separated by ~2 m. All the colonies at least partially enclosed in the belt transects were counted and identified to the genus level (18 and 19 genera were recorded during the study period for juvenile and adult corals, respectively; see Supplementary Table S7)¹⁸. Our analysis of coral cover, and abundance of juvenile and adult colonies focused on the four major coral genera (*Pocillopora*, *Acropora*, *Porites*, and *Montipora*) that together constitute >80% of the overall coral assemblages in this reef system.

To assess variation in coral settlement, 20 unglazed terracotta tiles (11 × 11 × 1 cm) were deployed at each station following the direct attachment method³⁸. Two consecutive sets of tiles were deployed for three months (September to December and December to March)²⁵ annually over a period of 13 years (2001–2002 to 2013–2014). This time frame was chosen because the main recruitment period in Moorea is December to March, with comparatively fewer recruits in September–December²⁵. For the purpose of this survey which focused on interannual changes rather than seasonal variability, we summed the number of coral recruits over these two consecutive periods (Sept. –Dec. and Dec. –March) to estimate ‘annual’ recruitment rates. Once collected from the reefs, settlement tiles were bleached and sun-dried, and then coral recruits present on each surface of the tiles (upper-side, under-side, and edges) were counted and identified using a dissecting microscope. At this stage of development, the micro-architecture of the corallum is not sufficiently developed to allow high resolution identification, and only three families (Acroporidae, Pocilloporidae, Poritidae) of coral recruits can be reliably distinguished³⁹. All other families that could not be consistently distinguished were compiled into a category of “others”. At Moorea, the family Acroporidae only encompasses three genera (*Acropora*, *Montipora*, and the uncommon *Astreopora*), whereas Pocilloporidae and Poritidae are only represented by one genus (*Pocillopora* and *Porites*, respectively).

**Statistical analyses.** We used generalized linear mixed effects models (GLMMs) to evaluate if the trajectories (covariable time) of coral taxa at the different life stages varied with reef locations (Tiahura, Vaipahu and Haapiti), depths (6, 12, and 18 m) and stations (location × depth). Coral abundances were modelled using negative binomial distribution, as a number of coral colonies was recorded each year in an additive fashion.

For each coral taxa and stages, four models (see equations below) were used to fit the observed data. The model 0, called the null-model, fit the abundance of coral recruits, juveniles or adults without explanatory variables. In this case, the variability in counts of coral colonies was purely random. In Model 1, Time was added as an explanatory variable, in order to take into account the temporal variability. In Models 2 and 3, Location and Depth were added to Model 1, respectively. In Model 4, Station was added as an explanatory variable to Model 1. Random slopes and intercepts and measurements error terms were also formulated for each model. Each explanatory variable was considered as categorical with Time centered on the year 2007 and Haapiti 6 m used as an arbitrary Location and Depth reference.

The abundance (Ab) for each coral stages (f) at a station i and at time j was modelled as follows:

Model 0: no covariables (no significant spatio-temporal pattern)

Level 1: $Ab_{fjj} = \beta_{0i} + \varepsilon_{ij}$

Level 2: $\beta_{0i} = \beta_0 + u_{0i}$

Model 1: Time (temporal variability)

Level 1: $Ab_{fjj} = \beta_{0i} + \beta_{1i}(Time)_j + \varepsilon_{ij}$

Level 2: $\beta_{0i} = \beta_0 + u_{0i}$

$\beta_{1i} = \beta_1 + u_{1i}$

Model 2: Location (different trajectories among locations)

Level 1: $Ab_{fjj} = \beta_{0i} + \beta_{1i}(Time)_j + \varepsilon_{ij}$

Level 2: $\beta_{0i} = \beta_0 + \beta_2(Location)_j + u_{0i}$

$\beta_{1i} = \beta_1 + \beta_2(Location)_{ij} + u_{1i}$

Model 3: Depth (different trajectories among depths)

Level 1: $Ab_{fjj} = \beta_{0i} + \beta_{1i}(Time)_j + \varepsilon_{ij}$

Level 2: $\beta_{0i} = \beta_0 + \beta_2(Depth)_{ij} + u_{0i}$

$\beta_{1i} = \beta_1 + \beta_2(Depth)_{ij} + u_{1i}$

Model 4: Station (different trajectories among stations, i.e. per depths within locations)

Level 1: $Ab_{fjj} = \beta_{0i} + \beta_{1i}(Time)_j + \varepsilon_{ij}$

Level 2: $\beta_{0i} = \beta_0 + \beta_2(Location)_{ij} + \beta_3(Depth(Location))_{ij} + u_{0i}$

$\beta_{1i} = \beta_1 + \beta_2(Location)_{ij} + \beta_3(Depth(Location))_{ij} + u_{1i}$
The best model formulations were identified using Akaike Information Criterion (AIC) and Chi² deviance tests⁴⁰. Model validity was controlled by plotting the estimated data vs. observed data and model residuals. All statistics were performed in R version 3.1.0⁴¹ complemented by the package glmmADMB⁴².

Model outputs were then used to determine coral recovery using pairwise comparisons where the first year of survey was considered as reference. Post-hoc analyses were conducted at the level of stations, depths (all locations pooled), locations (all depths pooled), and for all stations pooled. The multcomp package⁴³ was used to perform these tests.

Finally, we used non-metric multidimensional scaling (nMDS) to examine interannual variation in the composition of the coral assemblage based on the Bray-Curtis dissimilarity index of the cover of coral genera recorded annually. This analysis was encoded in R version 3.1.0⁴¹.

**Data availability.** Data used in this paper can be accessed at http://www.criobe.pf.

**References**

1. Bellwood, D. R., Hughes, T. P., Folke, C. & Nystrom, M. Confronting the coral reef crisis. *Nature* **429**, 827–833 (2004).

2. Hoeg-Guldberg, O. & Bruno, J. F. The impact of climate change on the world's marine ecosystems. *Science* **328**, 1523–1528 (2010).

3. Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. *Nature* **543**, 373–377 (2017).

4. McManus, J. W. & Polsenberg, J. F. Coral-algal phase shifts on coral reefs: ecological and environmental aspects. *Prog. Oceanogr.* **60**, 263–279 (2004).

5. Hughes, T. P. et al. Coral reefs in the Anthropocene. *Nature* **546**, 82–90 (2017).

6. Hughes, T. P. et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. *Science* **359**, 80–83 (2018).

7. Hoeg-Guldberg, O. et al. Coral reefs under rapid climate change and ocean acidification. *Science* **318**, 1737–1742 (2007).

8. van Hooidonk, R. et al. Local-scale projections of coral reef futures and implications of the Paris Agreement. *Sci. Rep.* **6**, 39666 (2016).

9. Bruno, J. F. & Selig, E. R. Regional decline of coral cover in the Indo-Pacific: timing, extent, and subregional comparisons. *PLoS ONE* **2**(8), e711 (2007).

10. De'ath, G., Fabricius, K. E., Sweatman, H. & Puotine, M. The 27-year decline of coral cover on the Great Barrier Reef and its causes. *Proc. Natl. Acad. Sci. USA* **109**, 17995–17999 (2012).

11. Adjeroud, M. et al. Recurrent disturbances, recovery trajectories, and resilience of coral assemblages on a south Central Pacific reef. *Coral Reefs* **28**, 775–780 (2009).

12. Graham, N. A. J., Nash, K. L. & Kool, J. T. Coral reef recovery dynamics in a changing world. *Coral Reefs* **30**, 283–294 (2011).

13. Gilmour, J. P., Smith, L. D., Heyward, A. J., Baird, A. H. & Pratchett, M. S. Recovery of an isolated coral reef system following severe disturbance. *Science* **340**, 69–71 (2013).

14. Mumby, P. J., Hastings, A. & Edwards, H. J. Thresholds and the resilience of Caribbean coral reefs. *Nature* **450**, 98–101 (2007).

15. Connell, J. H., Hughes, T. P. & Wallace, C. C. A 30-year study of coral abundance, recruitment, and disturbance at several scales in space and time. *Ecol. Monogr.* **67**, 461–488 (1997).

16. Adjeroud, M., Kayal, M. & Penin, L. The importance of recruitment processes in the dynamics and resilience of reef coral assemblages. In *Marine Animal Forests*, S. Rossi, L. Bramanti, A. Gori, C. Orejas, Eds, 549–569 (Springer, 2017).

17. Graham, N. A. J., Jennings, S., MacNeil, M. A., Mouillot, D. & Wilson, S. K. Predicting climate-driven regime shifts versus rebound potential in coral reefs. *Nature* **518**, 94–97 (2015).

18. Penin, L. et al. Early post-settlement mortality and the structure of coral assemblages. *Mar. Ecol. Prog. Ser.* **408**, 55–64 (2010).

19. Kayal, M. et al. Predator Crown-of-Thorns Starfish (*Acanthaster planci*) outbreak, mass mortality of corals, and cascading effects on reef fish and benthic communities. *PLoS ONE* **7**(10), e47363 (2012).

20. Berumen, M. L. & Pratchett, M. S. Recovery without resilience: persistent disturbance and long-term shifts in the structure of fish and coral communities at Tiahura Reef, Moorea. *Coral Reefs* **25**, 647–653 (2006).

21. Carroll, A. G., Harrisson, P. L. & Adjeroud, M. Susceptibility of coral assemblages to successive bleaching events at Moorea, French Polynesia. *Mar. Freshw. Res.* **68**, 760–771 (2017).

22. Edmunds, P. J. Unusually high coral recruitment during the 2016 El Niño in Moorea, French Polynesia. *PLoS ONE* **12**(10), e0185167 (2017).

23. Holbrook, S. J., Schmitt, R. J., Adam, T. C. & Brooks, A. J. Coral reef resilience, tipping points and the strength of herbivory. *Sci. Rep.* **6**, 35817 (2016).

24. Lamy, T., Galzin, R., Kulbicki, M., Lison de Loma, T. & Claudet, J. Three decades of recurrent declines and recoveries in corals belie ongoing change in fish assemblages. *Coral Reefs* **35**, 293–302 (2016).

25. Adjeroud, M., Penin, L. & Carroll, A. Spatio-temporal heterogeneity in coral recruitment around Moorea, French Polynesia: implications for population maintenance. *J. Exp. Mar. Biol. Ecol.* **341**, 204–218 (2007).

26. Connell, J. H. Disturbance and recovery of coral assemblages. *Coral Reefs* **16**(suppl), S101–S113 (1997).

27. Doropoulos, C., Ward, S., Roff, G., Gonzales-Rivero, M. & Mumby, P. J. Linking demographic processes of juvenile corals to benthic recovery trajectories in two common reef habitats. *PLoS ONE* **10**(5), e0128535 (2015).

28. Lukoschek, V., Cross, P., Torda, G., Zimmerman, R. & Willis, B. L. The importance of coral larval recruitment for the recovery of reefs impacted by cyclone Yasi in the Central Great Barrier Reef. *PLoS ONE* **8**(6), e65363 (2013).

29. Magalon, H., Adjeroud, M. & Veille, M. Patterns of genetic variation do not correlate with geographical distance in the reef building coral *Pocillopora meandrina* in the South Pacific. *Mol. Ecol.* **14**, 1861–1868 (2005).

30. Tsounis, G. & Edmunds, P. J. The potential for self-seeding by the coral *Pocillopora* spp. in Moorea, French Polynesia. *PeerJ* **4**, e2544 (2016).

31. Bongaerts, P. et al. Deep reefs are universal refuges: Reseeding potential varies among coral species. *Sci. Adv.* **3**, e1602373 (2017).

32. Bramanti, L. & Edmunds, P. J. Density-associated recruitment mediates coral population dynamics on a coral reef. *Coral Reefs* **35**, 543–553 (2016).

33. Darling, E. S., Alvarez-Filip, L., Oliver, T. A., McClanahan, T. R. & Côté, I. M. Evaluating life-history strategies of reef corals from species traits. *Ecol. Lett.* **15**, 1378–1386 (2012).

34. Kayal, M., Vercelloni, J., Wand, M. P. & Adjeroud, M. Searching for the best bet in life-strategy: A quantitative approach to individual performance and population dynamics in reef-building corals. *Ecol. Complex.* **23**, 73–84 (2015).

35. Graham, N. A. J., Cinner, J. E., Norström, A. V. & Nystrom, M. Coral reefs as novel ecosystems: embracing new futures. *Curr. Opin. Env. Sust.* **7**, 9–14 (2014).

36. Schrimm, M., Buscail, R. & Adjeroud, M. Spatial variability of the biogeochemical composition of surface sediments in an insular coral reef ecosystem: Moorea, French Polynesia. *Estuar. Coast. Shelf Sci.* **60**, 515–528 (2004).

37. Gil, M. A., Goldenberg, S. U., Thai Bach, A. L., Mills, S. C. & Claudet, J. Interactive effects of three pervasive marine stressors in a post-disturbance coral reef. *Coral Reefs* **35**, 1281–1293 (2016).

38. Mundy, C. N. An appraisal of methods used in coral recruitment studies. *Coral Reefs* **19**, 124–131 (2000).
39. Babcock, R. C., Baird, A. H., Piromvaragorn, S., Thomson, D. P. & Willis, B. L. Identification of Scleractinian coral recruits from Indo-Pacific reefs. *Zool. Stud.* **42**, 211–226 (2003).

40. Akaike, H. A Bayesian analysis of the minimum AIC procedure. *Ann. Inst. Stat. Math.* **30**, 9–14 (1978).

41. R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2014).

42. Bolker, B., Skaug, H., Magnusson, A. & Nielsen A. Getting started with the glmmADMB package, http://glmmadmb.r-forge.r-project.org/glmmADMB.html (2012).

43. Hothorn, T., Bretz, F., Westfall, P., Heibergeer, R. M. & Schuetzenmeister A. Simultaneous Inference in General Parametric Models, package "Multcomp", http://cran.r-project.org/web/packages/multcomp/multcomp.pdf (2014).

**Acknowledgements**

We thank Jane Ballard for constructive comments and editorial help. This research was supported by grants from the Total Foundation, CRISP, BioAqua, and GOPS.

**Author Contributions**

Conception and design: M.A., L.P. and M.K. Data collections: M.A., L.P., M.K., P.B., V.L. and Y.C. Analysis and interpretation of data: M.A., L.P., P.B., M.K., J.V., J.C. and C.I.C. Drafting manuscript: M.A., L.P. and M.K. All authors reviewed the manuscript.

**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-018-27891-3.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018