Conseils de recherche
Les critères de recherche 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2017; 12(3): e0173513.
Published online 2017 March 8. doi:  10.1371/journal.pone.0173513
PMCID: PMC5342305

Population structure of the hydrocoral Millepora platyphylla in habitats experiencing different flow regimes in Moorea, French Polynesia

Christina A. Kellogg, Editor


While the fire coral Millepora platyphylla is an important component of Indo-Pacific reefs, where it thrives in a wide range of environments, the ecological and biological processes driving its distribution and population structure are not well understood. Here, we quantified this species’ population structure in five habitats with contrasting hydrodynamic regimes in Moorea, French Polynesia; two in the fore reef: mid and upper slopes, and three in the lagoon: back, fringing and patch reefs. A total of 3651 colonies of fire corals were mapped and measured over 45,000 m2 of surveyed reef. Due to the species’ sensitivity to fragmentation in response to strong water movement, hydrodynamic conditions (e.g. waves, pass and lagoonal circulation) corresponded to marked differences in colony size distributions, morphology and recruitment dynamics among habitats. The size structure varied among reef habitats with higher proportions of larger colonies in calm nearshore reefs (fringing and patch reefs), while populations were dominated by smaller colonies in the exposed fore reefs. The highest densities of fire corals were recorded in fore reef habitats (0.12–0.20 n.m-2) where the proportion of recruits and juveniles was higher at mid slope populations (49.3%) than on the upper slope near where waves break (29.0%). In the latter habitat, most colonies grew as vertical sheets on encrusting bases making them more vulnerable to colony fragmentation, whereas fire corals were encrusting or massive in all other habitats. The lowest densities of M. platyphylla occurred in lagoonal habitats (0.02–0.04 n.m-2) characterized by a combination of low water movement and other physical and biological stressors. This study reports the first evidence of population structure of fire corals in two common reef environments and illustrates the importance of water flow in driving population dynamic processes of these reef-building species.


Coral reefs exhibit a remarkable diversity of organisms that reside within highly variable environments resulting in strong spatial variability in species’ distribution patterns [1]. For scleractinian corals, spatial differences in temperature, light, water flow and water quality conditions can influence their distribution and population dynamics [25]. Millepora hydrocorals, also called fire corals, are an important component of reefs communities worldwide where they, similar to scleractinian corals, contribute to reef accretion and community dynamics [6,7]. Fire corals can colonize a wide range of reef environments through sexual reproduction [7,8] and colony fragmentation [7,9]. Fire corals have been reported to grow faster than scleractinian corals [7,10] and often grow into large colonies that preempt space and compete with scleractinian corals [11,12]. On the other hand Millepora species also contribute to the survival of corals during Acanthaster outbreaks as this corallivorous predator tends to avoid Millepora species [7,13].

Hydrodynamic forces in the form of water-displacement, velocity and acceleration have been recognized as a key factor in determining the shape and occurrence of many reef-building organisms [1416]. In coral reef ecosystems, the magnitude of water flow is mostly related to the wave energy dispersal [5]. On barrier reefs, the amount of wave energy is highest on the reef crest, where waves break, and subsequently attenuates towards fore reef and lagoonal environments [5,17]. Inside lagoons, internal waves and flows drive circulation and water exchange with the surrounding ocean [18,19]. Such variation in hydrodynamic regimes, combined with other physical (e.g. light, nutrients and disturbances) and biological factors (e.g. colony size and shape), differently affect the performance of reef benthic organisms [2023] resulting in corresponding changes in population structure and community composition.

Water flow can drive the spatial distribution in adult populations through the distribution and dilution of larval settlement cues [24] and dispersal of reproductive propagules [25,26]. Many studies have related the contribution of recruitment to colony size variation in scleractinian corals (e.g., [27,28]) and the size structure of a population often reflects other species specific responses to environmental conditions and disturbances as well [2931]. The size-frequency distributions of fire coral populations could therefore provide insights on which biotic (e.g., recruitment of larvae and asexually produced fragments) and abiotic (e.g., wave energy) factors influence their population structure and dynamics.

Water flow also influences colony growth and morphology [32,33]. Under the increasing influence of hydrodynamics, delicate branching corals transform into growth forms able to withstand strong water movement such as compact, robust plating or thick branching morphologies [34,35]. Such inter- and intraspecific variation resulting in different coral morphologies affects not only their mechanical strength but also their ability to compete for space [3638] and capture light and food [23,39]. Branching and plating scleractinian corals, such as Pocillopora damicornis and Acropora hyacinthus, often grow into large and delicate arborescent colonies in areas of relatively high water flow [23,39,40], but this growth strategy also renders them extremely vulnerable to breakage when large waves and storm events occur, often resulting in fragmentation with some mortality [41]. Asexual reproduction through colony fragmentation can be a successful reproductive strategy to sustain local population growth in some species of scleractinian corals [42,43]. Fire corals are also known for their extensive intra- and interspecific morphological variability across hydrodynamic gradients with consequences due to their vulnerability to wave-induced breakage [7,9,36,44]. Determining to what degree the population structure of fire corals depends on the differences in water flow among common reef habitats has so far not been determined.

In this study, we investigated whether and how different reef habitats with contrasting water regimes affect the population structure of Millepora platyphylla, Hemprich & Ehrenberg 1834, the only species of fire coral found in French Polynesia [45,46]. Surveys of M. platyphylla were conducted in five habitats on the north shore of Moorea (Society Archipelago, French Polynesia) with differing amounts of water flow: fore reef habitats with high water movement, especially on the upper slope and decreasing with depth to the mid slope. Lagoonal habitats (back reefs, fringing reefs and patch reefs) are sheltered from waves and oceanic swell, except during storms, and water movement in these habitats is less than on the fore reef [5,17]. We examined colony size distribution, morphological variability and recruitment dynamics to assess to what degree the variability in the population structure of M. platyphylla among reef habitats can be attributed to different flow regimes.

Materials and methods

Model species

Millepora platyphylla is a gonochoric broadcast spawner that reproduces sexually by producing medusoids and planula larvae [7]. The medusoids are released into the water column and the gametes are released in one hour post-spawning during the medusoid’s swim. Then, external fertilization and embryogenesis occur after which the larvae sink and move epibenthically (i.e. crawling) on the reef substratum and metamorphose into calcifying polyps within one day after spawning [8]. M. platyphylla can also reproduce asexually through fragmentation [7,9].

Study sites and field surveys

Between April and December 2013, a series of surveys were conducted on the north shore of Moorea, French Polynesia, in the South Pacific Ocean (17,5267 S, 149,8348 W), at four different locations (Tiahura, Papetoai, Cook’s Bay and Temae). Five habitats with contrasting water flow regimes were selected; two in a fore reef environment: the mid slope (13 m depth) and the upper slope (6 m depth), and three in the lagoon (< 1 m depth): the back reef, fringing reef and patch reef (Fig 1 and S1 Table). These habitats greatly differ in terms of water flow. The fore reef experiences strong wave action from incoming waves that break on the reef crest with gradual swell wave attenuation towards deeper waters [5,47]. Because of this strong linear relationship between wave forcing and water depth [47], the colonies of M. platyphylla growing within fore reef habitats are exposed to lower wave energy on the mid slope compared to the those growing on the upper slope, near where the waves break. In the lagoon, the wave energy disperses from the reef crest towards nearshore reefs [5]. A recent study on wave energy across reef environments revealed that the reef crest dissipated 70% of the incident swell wave energy with gradual wave attenuation from the back reef to the shore [17]. Consequently, we assumed that wave energy is higher on the back reef, near the reef crest compared to the fringing reef, a nearshore reef. Although the patch reef is located in a nearshore narrow channel, the wave energy there is also higher than on the fringing reef due to its proximity to the reef crest and to the currents that run on either side of the channel (i.e. pass circulation). Variations in other physical constraints exist between the fore reef and lagoonal habitats in terms of e.g., temperature, water clarity, nutrient and disturbances [48,49], which make them highly contrasting reef environments.

Fig 1
Aerial views of the locations of each transect in the five surveyed habitats in Moorea, French Polynesia.

Within each habitat, three 300 m long by 10 m wide belt transects were laid over the reef parallel to shore, at least 30 m apart, resulting in a total of 45,000 m2 reef area being surveyed. The north shore of Moorea is ~16 km long with a fore reef area of ~3.15 km2 and a back reef with a hard-bottom area of ~4.58 km2 [50]. We performed six belt transects of 0.003 km2 on the fore reef, which is ~0.1% of the total fore reef area on the north shore of Moorea. We also performed six belt transects of 0.003 km2 on the back reef area (i.e. both back and patch reef habitats), which is ~0.07% of the total back reef area. All colonies of M. platyphylla that were at least 50% within the transect borders were measured, photographed and georeferenced using SCUBA. No specific permit was required at the time of field work for sampling protocols described herein and our surveys did not involve endangered/protected species and did not require animal tissue/skeleton collection.

Spatial distribution patterns

All M. platyphylla colonies were georeferenced by determining their position along the transect-line (0 to 300 m) and straight-line distance from both sides of the transect (0 to 10 m). From these measures, each colony was mapped with x and y coordinates, from which the distribution index (DI) and mean neighborhood distance (ND) were calculated using the spdep package [51] in R [52]. The DI is based on Ripley’s method [53] and calculated for each transect to determine whether colonies were having a contagious (DI > 1), random (DI ≈ 1) or homogenous (DI < 1) pattern of distribution [54]. The mean distance to each colony’s 10 nearest neighbors was estimated and the mean ND was calculated for each transect. The mean colony density (n.m-2) and cover (%) were also calculated for each transect (i.e. 3000 m2). Using these variables, variability in the spatial distribution among habitats was quantified by one-way PERMANOVA tests in PRIMER 6 software [55], since assumptions of parametric testing could not be met. Pair-wise tests followed the PERMANOVA to assess the degree of similarity among habitats. In order to determine how different habitats with contrasting water regimes affect the spatial distribution of M. platyphylla, we assumed that swell wave energy exposure decreases with habitat depth and its proximity to the coastline, as demonstrated in previous studies [17,47]. Consequently, the density, cover, DI and ND were regressed against the mean depth and mean distance from shore estimated from the three transects within each of the five surveyed habitats and Pearson’s r coefficient was used to test for significant correlations.

Colony size distribution

The size-frequency distributions of M. platyphylla populations were generated from estimates of colony sizes computed from 2D photographs. Photographs were taken from above the colony and included a plate of known dimensions positioned next to each colony. For bigger colonies, pictures were taken from a larger distance, and for 5 colonies (out of the 3561) 2 photographs were required to photograph the whole colony. Each colony size, standardized as the projected surface, was then measured (in cm2) using ImageJ 1.4f software [56]. The size-frequency distribution for each transect was given as percentages of all colonies belonging to 10 size classes on a logarithmic scale. Data were then analyzed using basic statistical measures of size hierarchies [57]: the coefficient of variation (CV) and skewness (g1), indicative of the relative abundance of small and large colonies within a population. CV and g1 were computed for each habitat per transect together with standard descriptive statistics, such as 95% percentile of the mean (describes the maximum colony size reached within a population, see [58]) and the probability that the data are normally distributed (Kolmogorov-Smirnov test, Pnorm). Differences in size-frequency distributions among habitats were quantified using one-way PERMANOVA based on normalized abundances. Spearman’s rank coefficient and pair-wise tests followed the PERMANOVA to assess the degree of similarity among habitats.

Recruitment dynamics

The mean abundance and proportion of recruits, juveniles and adults were estimated for each transect whereby the three life stages were defined based on colony size. Colonies with a total size (surface) below 1 cm2 were considered as recruits and were most likely the result of sexual reproduction. Larger colonies with a size between 1 and 20 cm2 were classified as juveniles based on previous studies on coral recruitment using both settlement plate experiments [59] and field surveys [60]. While the origin of each colony (sexual or asexual) could not be confirmed from field surveys, both size classes were considered as non-reproductive in contrast to colonies above 20 cm2 based on previous studies on other reef-building taxa [61]. Differences in abundances and proportions (i.e. the fraction in the entire population) of early life stages (i.e. both recruits and juveniles) among habitats were quantified using one-way PERMANOVA, followed by a pair-wise test. Pearson’s correlation coefficient was used to determine whether the abundance of early life stages increased with the abundance and cover of adults, and whether differences in their proportions among habitats correlate with water movement, i.e. with depth and distance from shore used as proxies.

Colony morphology

For each colony’s morphology, the maximum height, from the base to the highest part of the colony (rounded to the nearest half cm), was recorded and linked to the colony size data previously described. Colonies below 20 cm2 were removed from this analysis to only retain the mean height and size of adults for each transect. Adult colonies were assigned to one of these three morphologies: 1) encrusting: thin colonies growing against the substratum (Fig 2A); 2) “sheet tree”: encrusting bases with platelike outgrowths facing wave energy (see [36]) (Fig 2B) and 3) massive: solid colonies, roughly hemispherical in shape (Fig 2C). Differences in proportions of each of the three morphologies among habitats were quantified with one-way PERMANOVA and pair-wise tests.

Fig 2
Morphologies of Millepora platyphylla colonies in habitats experiencing contrasting hydrodynamic regimes.

Population structure assessment

Similarities in population structure based on the following parameters: density, cover, DI, ND, mean adult colony size and height, and proportion of recruits, juveniles and adults were calculated and visualized using a hierarchical complete-linkage agglomerative clustering (CLUSTER) method and a non-parametric multidimensional scaling (MDS) ordination on normalized data in PRIMER 6 software. Multivariate PERMANOVA on aforementioned characteristics was used to determine differences in population structure of M. platyphylla among the five surveyed habitats, i.e., those on the fore reef (mid and upper slopes) and those in lagoonal habitats (back, fringing and patch reefs).


Spatial distribution of Millepora platyphylla

M. platyphylla was found in all habitats, but its population composition differed among habitats. A total of 3651 colonies of M. platyphylla were counted in the five surveyed habitats. Most colonies (48.2%) occurred on the upper slope, whereas M. platyphylla colonies on patch reefs accounted only for 5.2% of all colonies (S2 Table). Colony density differed among habitats (PERMANOVA test, P < 0.01) and was higher on the upper slope (0.20 ± 0.03 n.m-2, N = 1761) and mid slope (0.12 ± 0.05 n.m-2, N = 1075), i.e. fore reef habitats, compared to lagoonal habitats (back reef: 0.03 ± 0.01 n.m-2, N = 324, fringing reef: 0.04 ± 0.03 n.m-2, N = 302 and patch reef: 0.02 ± 0.00 n.m-2, N = 189) (Fig 3A). M. platyphylla’s cover also differed among habitats (PERMANOVA test, P < 0.01) and was again highest on the upper slope (3.2 ± 0.4%, Fig 3B). Colonies on the fringing reef, mid slope and upper slope occurred in a contagious pattern of distribution (DI: 2.74–4.18), while colonies in the back and patch reefs were more evenly distributed (≤ 1.93) (Fig 3C, PERMANOVA test, P < 0.05). Colonies occurred closer together on the mid slope (6.64 ± 1.86 m) and upper slope (4.09 ± 0.34 m) where the average distance among neighboring fire coral colonies was 4.3 times smaller compared to lagoonal habitats (back reef: 18.39 ± 1.10 m, fringing reef: 14.31 ± 4.41 m and patch reef: 36.51 ± 2.95 m) (Fig 3D, PERMANOVA test, P < 0.01).

Fig 3
Index describing the spatial distribution of Millepora platyphylla colonies across the five surveyed habitats.

Size structure of Millepora platyphylla

Across all habitats, 85% of the surveyed colonies were smaller than 1000 cm2 and approximately one third (30%) of aforementioned colonies fell in recruit and juvenile size classes, i.e. were smaller than 20 cm2. The size-frequency distributions of M. platyphylla populations differed among certain habitats, but were similar among the lagoonal fringing and patch reefs with populations dominated by a combination of small (≤ 32 cm2, comprising both small recruits and juveniles) and large colonies (≥ 2050 cm2) resulting in bimodal size-frequency distributions (Fig 4, Spearman’ rank coefficient 87.9%, P < 0.05). All fire coral populations were characterized by relatively symmetrical size distributions (g1: –0.01–0.71), but the degree of skewness was again lower on the fringing and patch reefs (Table 1). The maximum colony size differed among habitats and was smallest on the mid slope (95%: 1295 cm2) and back reef (2512 cm2) compared to other populations (upper slope: 8514 cm2, fringing reef: 7107 cm2 and patch reef: 9890 cm2) (Table 1). With 64% of all colonies falling in a few medium size classes (32–512 cm2, Fig 4), colonies comprising back reef populations were very similar relative to each other as indicated by the lowest coefficient of variation (CV: 0.33, Table 1) of all habitats. Overall, the composition of M. platyphylla populations in terms of colony density and size differed among the five habitats, except between the two lagoonal habitats, the fringing and patch reefs, located closest to shore.

Fig 4
Size-frequency distributions of Millepora platyphylla across the five surveyed habitats.
Table 1
Index describing the population size structure and recruitment for Millepora platyphylla across the five habitats surveyed.

Recruitment of Millepora platyphylla

In total, 71 recruits (2%) and 1094 juveniles (30%) were observed within the five surveyed habitats (S2 Table). The abundance of recruits and juveniles differed among habitats (PERMANOVA test, P < 0.05 and P < 0.01, respectively) with 96% of all recruits and juveniles occurring in fore reef habitats (48% for both the mid and upper slopes) and only small numbers were observed in lagoonal habitats (Table 1). The fraction of the entire population consisting of recruits and juveniles differed among habitats (PERMANOVA test, P < 0.01). Recruits were found in low proportions in most habitats with the highest value recorded on the mid slope (3.4 ± 0.4%), while no recruit was observed on the back reef (Table 1; Fig 5). The mid slope habitat sheltered the highest proportion of juvenile colonies (45.9 ± 5.1%), while lower values were found in all other reef habitats (11.8–27.3%) with no significant difference between the fringing reef (24.0 ± 0.2%) and upper slope (27.3 ± 3.5%, Table 1; Fig 5). Only in lagoonal habitats did the abundance of adults, not their total cover, and the abundance of both recruits and juveniles increase simultaneously suggesting the presence of a stock-recruitment relationship (Fig 6A), that was not observed in fore reef habitats (Fig 6B). No significant stock recruitment relationship was found in both lagoonal and fore reef habitats when only considering recruit and adult’s abundances.

Fig 5
Recruitment dynamics across the five surveyed habitats.
Fig 6
Stock-recruitment relationship between the abundance of adults and coral new recruits and juveniles.

Morphology of Millepora platyphylla

M. platyphylla colonies ranged in size from 0.18 cm2 to 189 062 cm2 (projected surface) and 0.1 cm up to 130 cm in height. Mean colony size and height of adults (i.e., all colonies > 20cm2) differed among habitats (PERMANOVA tests, P < 0.01). Fire corals were approximately 4 times larger on average in the upper reef slope (2308 ± 115 cm2), fringing (1590 ± 329 cm2) and patch reef (3090 ± 1294 cm2), compared to colonies growing in the back reef and mid slope (510 ± 97 cm2 and 819 ± 73 cm2, respectively, Table 1). The average height of fire coral colonies was highest in fringing (24 ± 5 cm) and patch reefs (25 ± 8 cm), i.e. nearshore habitats (S2 Table). Morphologies of adult colonies differed among habitats (PERMANOVA test, P < 0.01). Massive morphologies dominated nearshore reefs (fringing reef: 79.7 ± 8.3% and patch reef: 59.0 ± 9.9%) whereas colonies were mostly encrusting on the mid slope (79.9 ± 1.1%) and back reef (74.5 ± 5.2%) (Fig 7 and S3 Table). On the fringing reef, no colony with the sheet tree morphology was observed. On the upper slope, 69.5% (± 3.2) of the colonies displayed the sheet tree morphology, while the remaining colonies were only encrusting.

Fig 7
Morphology of Millepora platyphylla adult colonies across the five surveyed habitats.

Population structure assessment

Combining all variables into a single multivariate analysis, the population structure of M. platyphylla varied significantly among reef habitats (PERMAVOVAs, P < 0.01). Based on MDS, two main clusters can be distinguished: one with populations from wave exposed fore reef habitats, i.e. mid and upper slopes, and a second cluster consisting of populations from lagoonal habitats, i.e. back, fringing and patch reefs (Fig 8). The main differences between these two clusters are that fore reef populations are characterized by a high relative abundance of recruits and juveniles (mid slope), or a higher density and cover (upper slope). Populations from lagoonal habitats are characterized by large colony size and height (both fringing and patch reefs) and widely spaced colonies. Back reef populations are characterized by the dominance of adult colonies.

Fig 8
Non-metric multidimensional scaling (MDS) plot of Millepora platyphylla population structure across the five surveyed habitats.

Using all variables of each transect surveyed within the fore reef habitat (i.e. six replicates), we found that adult colonies became smaller (i.e. colony size decreases) with increasing depth (r = –1.00, P < 0.001; N = 6). Fractions of recruits and juveniles increased with increasing depth (r = 0.91, P < 0.05 and r = –0.92, P < 0.05; N = 6), while total cover decreased (r = –0.97, P < 0.01; N = 6). Colonies grew in an encrusting morphology at mid depth and in the sheet tree morphology in shallow waters (Fig 7). Among shallow lagoonal habitats, we found that adult colonies became smaller towards the back reef, far from shore (r = –0.76, P < 0.01; N = 9), where wave energy was higher and colonies mostly occurred in the encrusting morphology (Fig 7). Total cover increased with increasing distance from the coast (r = –0.69, P < 0.05; N = 9).


Distribution patterns: Fore reef versus lagoonal habitats

In Moorea, M. platyphylla colonized a wide range of habitats reflecting its ability to adapt and survive in a large variety of environmental settings. This study is, to our knowledge, the most extensive sampling ever conducted to assess local patterns in population structure of Milleporid corals. Reef habitats where M. platyphylla colonies were found were selected because of their difference in water regimes according to their depth and proximity to the coastline (see Materials and Methods for details). Due to M. platyphylla’s sensitivity, especially of larger colonies to fragmentation induced by wave action and/or water movement (i.e. currents) (Fig 2), we sought for possible relationships between hydrodynamic conditions and the population structure of fire corals on Moorea. Differences in population size structure, recruitment and morphology existed among habitats and confirmed expected relationships between such characteristics and the amount of water flow in several of the five surveyed habitats (i.e. mid slope, upper slope, back reef, fringing reef, and patch reef). The highest densities of fire corals, including that of recruits and juveniles, occurred on the exposed fore reef (i.e. mid and upper slopes) whereby colonies were often observed growing in contagious pattern of distribution. In calm lagoonal environments (i.e. back, fringing and patch reefs) fire coral colonies occurred in low densities, where the number of recruits and juvenile was low and colonies grew in a random pattern of distribution. Variability in density among fore reef and lagoonal habitats has been described for numerous other sessile organisms and related to a large number of environmental factors such as water flow, solar irradiance, sedimentation and/or species’ life history traits (e.g., reproductive mode, competitive ability, morphological plasticity; [20,6265].

Size structure and morphological variations

Although differences in size-frequency distributions among habitats were found, e.g. few larger colonies in the calm fringing and patch reefs, smaller colonies in the mid slope and medium size colonies in the back reef, the degree of skweness was similar among all habitats with populations consisting of both small and large colonies. These distributions likely reflect low mortality in small size classes, as well as the persistence of the larger ones [66]. Our results showed that the proportion of recruits and juveniles was highest on the mid slope, an exposed reef where wave energy is reduced due to increased depth [47]. Earlier reports have also shown the influence of depth and water flow on the recruitment dynamics in some scleractinian coral species in many reef locations [6769]. These studies revealed an increase in the occurrence of recruits and juveniles with increasing depth. Another study compiling juvenile data of all coral species surveyed in Palmyra Atoll (Central Pacific) has shown that most juveniles were growing at middle depth (i.e. 14 m) in a fore reef habitat [63], as for M. platyphylla. Water flow is also considered as an important factor influencing a colony’s morphology [70,71], generally showing a transition from easily fragmented morphologies towards more robust morphologies with increasing water movement [34,35]. This study shows a similar trend whereby large and high colonies were more common in protected nearshore habitats (i.e. fringing and patch reef) and small and encrusting in exposed mid slope and back reef habitats. On the upper slope, near where the waves break, fire corals are large, but largely encrusting, and of the unusual sheet tree morphology of Millepora that was only observed in low proportions in all other habitats (0–9%).

Population structure and dynamics

Fire corals, like many other reef-building organisms, reproduce through both asexual and sexual reproduction with a dimorphic life cycle, with a pelagic dispersive phase (i.e. medusoids and larvae), followed by a sessile adult phase [7,8]. If dispersal distances are small due to low water movement or retention, the spatial distribution of adults could influence the distribution of young colonies as previously shown for scleractinian corals [68,69,72,73]. On Moorea the abundance of M. platyphylla recruits and juveniles could not be related to adult population size in fore reef habitats. The proportions of recruits and juveniles were highest at mid-depths (13 m) on the fore reef where wave energy and solar irradiance are lower compared to the shallower depth habitats (as described in [63]). Low wave energy can indeed increase settlement success of both coral larvae and fragments [74]. At shallow depths (6 m) on the fore reef, high wave energy and irradiance can reduce the abundance of settlement cues [74], but also indirectly affect settler survival through high grazing pressure by herbivorous fishes at this depth which constitutes a major source of mortality for juvenile corals on the upper slope in Moorea [59]. The abundance of coral fragments that re-attached to the reefs can also be reduced due to high wave energy and subsequent increased mortality [75]. Such physical and biological constraints in a dynamic environment likely reduce local recruitment rates and prevent high coral cover [76]. However, the highest cover of M. platyphylla (3.2%) occurred on the upper slope where wave breaking first occurs, i.e. wave energy is the highest. Many studies investigating spatial distributions in coral reef communities often find that high energy reef zones restrict species’ distributions and cover [77,78]. M. platyphylla shows the opposite trend: we observed high density and cover in the upper slope, a high energy reef zone, where colonies are growing in a contagious pattern of distribution. Such differences in fire coral distribution patterns in habitats of high energy are mostly related to the wave-vulnerable sheet tree morphology of Millepora. This growth form occurred nearly exclusive on the upper slope, while colonies were massive or encrusting in other habitats. The unusual sheet tree morphology observed in the upper slope has been described as a successful strategy exploited by Millepora to preempt the space and to compete with other coral taxa [12,36]. Waves can easily break the blades and enhance population growth through clonal propagation [9], while the encrusting bases remain intact and grow through horizontal stolonal spreading [12]. The fact that M. platyphylla can rapidly overtake newly available space through clonal propagation and stolonal spreading may explain the increase of fire coral cover on Moorea’s reefs following the massive decline in coral cover from the Acanthaster outbreaks and cyclone Oli in 2010 [79]. Between 2006 and 2010, M. platyphylla’s cover was stable at approximately 1% at 6 m on the fore reef, i.e. more than 3 times lower than in 2013 at the same location. On the other hand, fragmentation usually induces corals to regress in size and increases mortality, especially in small size classes [80]. Here, the sheet tree morphology is more easily fragmented, but the unilateral growth of Millepora allows them to reach larger colony sizes. This study shows that asexual reproduction through fragmentation and stolonal spreading likely plays a key role in structuring M. platyphylla populations where water flow is high and where fire corals face wave-induced breakage.

In the lagoonal environment, the wave energy is reduced by the reef crest [17] likely explaining the positive stock-recruitment relationship found in these habitats. There is evidence showing that the fecundity in populations of sessile marine broadcast spawners, such as Millepora species, is strongly determined by the local density of adults [81,82], and especially where water movement is reduced and local retention occurs. The low abundance of early life stages observed in all lagoonal habitats may result from competition with macroalgae and sediment smothering affecting back, patch and fringing reefs inside the lagoon of Moorea [83] and/or solar irradiance [74,76]. The presence of macroalgae and high sedimentation can additionally reduce adults’ fecundity [8486], larval settlement cues [87], larval survival [88] and settlement space [89]. Poorer water quality compared to fore reef habitats also likely contributes to the low abundance of Milleporid corals inside the lagoon [90,91]. The back and patch reefs are the nearest to the reef crest where waves break resulting in low residence times and high flushing rates from large incoming waves that break on the north shore of Moorea during the austral summer [18]. These dynamics of water flow are known to negatively affect local recruitment rates of sexual propagules in scleractinian corals on the back reef of Moorea [25] and could also apply to M. platyphylla. In the lagoon, fire corals are characterized by wave-tolerant morphologies (i.e. encrusting and massive) suggesting that asexual reproduction through colony fragmentation is less likely of structuring importance compared to fore reef habitats. Colonies on the fringing reef, where wave energy is typically low, were distributed in patches. In Moorea, the fringing reef is exposed to large waves in the austral summer [18], which has the potential to enhance the breakage of the colonies during short periods. Subsequent calm periods can facilitate fragment survival and reattachment resulting in the patches of M. platyphylla observed.

Implications for population maintenance and recovery

It must be noted that abundance of recruits and juveniles was likely underestimated in this study as they are difficult to find due to their small size during field surveys. Still, we identified 32% of colonies below 20 cm2 (recruits and juveniles), a higher fraction than observed for 14 different genera of scleractinians corals in Moorea (13–29%, see [68]). The abundance of fire corals around Moorea is also higher compared to more diverse and healthy reefs, such as the Great Barrier Reef [48] and the shallow fringing reefs in the US Virgin Islands [92]. Our results thus suggest that fire coral populations are relatively resilient in the face of recent and major disturbances that have impacted Moorea’s reefs (Acanthaster outbreaks, cyclones and bleaching events). The maintenance and recovery of fire coral populations are foremost sustained by the growth of remnant colonies, local recruitment through sexual reproduction where wave energy is low and clonal propagation in high wave energy zones.

Supporting information

S1 Table

Locations of each transect surveyed in the five habitats.


S2 Table

Index describing the spatial distribution, recruitment and morphology for Millepora platyphylla across the five surveyed habitats.


S3 Table

Average percentages of adult colonies for Millepora platyphylla with encrusting, sheet tree and massive morphology across surveyed habitats.



We are grateful to Franck Lerouvreur and Marc Besson who contributed to the fieldwork and Emilie Boissin for valuable discussion. Thanks to the CRIOBE staff for technical and logistic support. We also thank CAF Bourmaud and one anonymous reviewer for insightful comments on the draft manuscript.

Funding Statement

CED is supported by a graduate scholarship from the Fonds Québécois de Recherche sur la Nature et les Technologies.

Data Availability

Data Availability

All relevant data are within the paper and its Supporting Information files.


1. Wilkinson C. Status of coral reefs of the world: 2008. Global Coral Reef Monitoring Network and Reef and Rainforest Research Center; Townsville, Australia; 2008.
2. Glynn PW. Coral reef bleaching: facts, hypotheses and implications. Glob Chang Biol. 1996; 2: 495–509.
3. Vermeij MJA, Bak RPM. How are coral populations structured by light? Marine light regimes and the distribution of Madracis. Mar Ecol Prog Ser. 2002; 233: 105–116.
4. Fabricius KE. Effects of terrestrial runoff on the ecology of corals and coral reefs: review and synthesis. Marine Poll Bull. 2005; 50: 125–146. [PubMed]
5. Monismith SG. Hydrodynamics of coral reefs. Annu Rev Fluid Mech. 2007; 39: 37–55.
6. Nagelkerken I, Nagelkerken WP. Loss of coral cover and biodiversity on shallow Acropora and Millepora reefs after 31 years on Curaçao, Netherlands Antilles. B Mar Sci. 2004; 74: 213–223.
7. Lewis JB. Biology and ecology of the hydrocoral Millepora on coral reefs. Adv Mar Biol. 2006; 50: 1–55. doi: 10.1016/S0065-2881(05)50001-4 [PubMed]
8. Bourmaud C, Leung JKL, Bollard S, Gravier-Bonnet N. Mass spawning events, seasonality and reproductive features in Milleporids (Cnidaria, Hydrozoa) from Reunion Island. Mar Ecol. 2013; 34: 14–24.
9. Edmunds PJ. The role of colony morphology and substratum inclination in the success of Millepora alcicornis on shallow coral reefs. Coral Reefs. 1999; 18: 133–140.
10. Attalla TM, Hanafy MH, Aamer MA. Growth rates of the two reef-building species, Acropora humilis and Millepora platyphylla at Hurghada, Red Sea, Egypt. Egypt J Aquat Biol Fish. 2011; 15: 1–5.
11. Andréfouët S, Benzoni F, Payri C. A monospecific Millepora reef in Marquesas Islands, French Polynesia. Coral Reefs. 2014; 33: 463.
12. Dubé CE, Boissin E, Planes S. Overgrowth of living scleractinian corals by the hydrocoral Millepora platyphylla in Moorea, French Polynesia. Mar Biodiv. 2016; 46: 329–330.
13. Kayal M, Kayal E. Colonies of the fire coral Millepora platyphylla constitute scleractinian survival oases during Acanthaster outbreaks in French Polynesia. Mar Biodiv. 2016; 1–4.
14. Madin JS, Black KP, Connolly SR. Scaling water motion on coral reefs: from regional to organismal scales. Coral Reefs. 2006; 25: 635–644.
15. Williams GJ, Smith JE, Conklin EJ, Gove JM, Sala E, Sandin SA. Benthic communities at two remote Pacific coral reefs: effect of reef habitat, depth, and wave energy gradients on spatial patterns. PeerJ. 2013; 1: e81 doi: 10.7717/peerj.81 [PMC free article] [PubMed]
16. Lenihan HS, Hench JL, Holbrook SJ, Schmitt RJ, Potoski M. Hydrodynamics influence coral performance through simultaneous direct and indirect effects. Ecology. 2015; 96: 1540–1549.
17. Ferrario F, Beck MW, Storlazzi CD, Micheli F, Shepard CC, Airoldi L. The effectiveness of coral reefs for coastal hazard risk reduction and adaptation. Nat Commun. 2014; 13: 5. [PMC free article] [PubMed]
18. Hench JL, Leichter JJ, Monismith SG. Episodic circulation and exchange in a wave-driven coral reef and lagoon system. Limnol Oceanogr. 2008; 53: 2681–2694.
19. Leichter JJ, Alldredge AL, Bernardi G, Brooks AJ, Carlson CA, Carpenter RC, et al. Biological and physical interactions on a tropical island coral reef: Transport and retention processes on Moorea, French Polynesia. Oceanography. 2013; 26: 52–63.
20. Gili JM, Hughes RG. The ecology of marine benthic hydroids. Oceanogr Mar Biol. 1995; 33: 351–426.
21. Hoogenboom MO, Connolly SR. Defining fundamental niche dimensions of corals: synergistic effects of colony size, light, and flow. Ecology. 2009; 90: 767–780. [PubMed]
22. Madin JS, Hoogenboom MO, Connolly SR. Integrating physiological and biomechanical drivers of population growth over environmental gradients on coral reefs. J Exp Biol. 2012; 215: 968–976. doi: 10.1242/jeb.061002 [PubMed]
23. Darling ES, Alvarez-Filip L, Oliver TA, McClanahan TR, Coté IM. Evaluating life-history strategies of reef corals from species traits. Ecol Lett. 2012; 15: 1378–1386. doi: 10.1111/j.1461-0248.2012.01861.x [PubMed]
24. Koehl MAR, Hadfield MG. Soluble settlement cue in slowly moving water within coral reefs induces larval adhesion to surfaces. J Mar Syst. 2004; 49: 75–88.
25. Edmunds PJ, Leichter JJ, Adjeroud M. Landscape-scale variation in coral recruitment in Moorea, French Polynesia. Mar Ecol Prog Ser. 2010; 414: 75–89.
26. Gleason DF, Hofmann DK. Coral larvae: From gametes to recruits. J Exp Mar Biol Ecol. 2011; 408: 42–57.
27. Meesters EH, Hilterman M, Kardinaal E, Keetman M, de Vries M, Bak RPM. Colony size-frequency distributions of scleractinian coral populations: spatial and interspecific variation. Mar Ecol Prog Ser. 2001; 209: 43–54.
28. Harris A, Wilson S, Graham N, Sheppard C. Scleractinian coral communities of the inner Seychelles 10 years after the 1998 mortality event. Aquat Conserv. 2014; 24: 667–679.
29. Hughes TP. Population dynamics based on size rather than age: a general model with a reef coral example. Am Nat. 1984; 123: 778–795.
30. Albright R, Langdon C. Ocean acidification impacts multiple early life history processes of the Caribbean coral Porites astreoides. Glob Chang Biol. 2011; 17: 2478–2487.
31. Madin JS, Hughes TP, Connolly SR. Calcification, storm damage and population resilience of tabular corals under climate change. PLoS One. 2012; 7: e46637 doi: 10.1371/journal.pone.0046637 [PMC free article] [PubMed]
32. Storlazzi CD, Brown EK, Field ME, Rodgers K, Jokiel PL. A model for wave control on coral breakage and species distribution in the Hawaiian Islands. Coral Reefs. 2005; 24: 43–55.
33. Madin JS, Connolly SR. Ecological consequences of major hydrodynamic disturbances on coral reefs. Nature. 2006; 444: 477–480. doi: 10.1038/nature05328 [PubMed]
34. Kaandorp JA, Sloot PMA. Morphological models of radiate accretive growth and the influence of hydrodynamics. J Theor Biol. 2001; 209: 257–274. doi: 10.1006/jtbi.2001.2261 [PubMed]
35. Chindapol N, Kaandorp JA, Cronemberger C, Mass T, Genin A. Modelling growth and form of the scleractinian coral Pocillopora verrucosa and the influence of hydrodynamics. PLoS Comput Biol. 2013; 9: e1002849 doi: 10.1371/journal.pcbi.1002849 [PMC free article] [PubMed]
36. Jackson JBC. Morphological strategies of sessile animals In: Larwood G, Rosen BR, editors. Biology and systematics of colonial organisms. Academic Press, London; 1979. pp. 499–555.
37. Denny MW. Biology and the mechanics of the wave-swept environment. Princeton University Press, New Jersey; 2014.
38. Swierts T, Vermeij MJA. Competitive interactions between corals and turf algae depend on coral colony form. PeerJ. 2016; 4: e1984 doi: 10.7717/peerj.1984 [PMC free article] [PubMed]
39. Baird AH, Hughes TP. Competitive dominance by tabular corals: an experimental analysis of recruitment and survival of understorey assemblages. J Exp Mar Biol Ecol. 2000; 251: 117–132. [PubMed]
40. Kaandorp JA. Morphological analysis of growth forms of branching marine sessile organisms along environmental gradients. Mar Biol. 1999; 134: 295–306.
41. Madin JS, Baird AH, Dornelas M, Connolly SR. Mechanical vulnerability explains size-dependent mortality of reef corals. Ecol Lett. 2014; 17: 1008–1015. doi: 10.1111/ele.12306 [PMC free article] [PubMed]
42. Baums IB, Miller MW, Hellberg ME. Geographic variation in clonal structure in a reef-building Caribbean coral, Acropora palmata. Ecol Monogr. 2006; 76: 503–519.
43. Aranceta-Garza F, Balart EF, Reyes-Bonilla H, Cruz-Hernández P. Effect of tropical storms on sexual and asexual reproduction in coral Pocillopora verrucosa subpopulations in the Gulf of California. Coral Reefs. 2012; 31: 1157–1167.
44. Weerdt WH. Transplantation experiments with Caribbean Millepora species (Hydrozoa, Coelenterata), including some ecological observations on growth forms. Bijdr Dierkd. 1981; 51: 1–19.
45. Bosserelle P, Berteaux-Lecellier V, Chancerelle Y, Hédouin L, Nugues M, Wallace C, et al. Guide d’identification des coraux de Moorea. CRIOBE, Moorea; 2014.
46. Dubé CE, Planes S, Zhou Y, Berteaux-Lecellier V, Boissin E. Genetic diversity and differentiation in reef-building Millepora species, as revealed by cross-species amplification of fifteen novel microsatellite loci. PeerJ. 2017; 5: e2936 doi: 10.7717/peerj.2936 [PMC free article] [PubMed]
47. Hearn CJ. Wave-breaking hydrodynamics within coral reef systems and the effect of changing relative sea-level. J Geophys Res. 1999; 104: 30007–30019.
48. Done TJ. Patterns in the distribution of coral communities across the central Great Barrier Reef. Coral Reefs. 1982; 1: 95–107.
49. Witman JD. Physical disturbance and community structure of exposed and protected reefs: A case study from St-John, U.S. Virgin Islands. Am Zool. 1992; 32: 641–654.
50. Tsounis G, Edmunds PJ. The potential for self-seeding by the coral Pocillopora spp. in Moorea, French Polynesia. PeerJ. 2016; 4: e2544 doi: 10.7717/peerj.2544 [PMC free article] [PubMed]
51. Bivand RS, Pebesma E, Gómez-Rubio V. Hello World: Introducing spatial data In: Applied Spatial Data Analysis with R. Springer, New York; 2013. pp. 1–16.
52. R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.
53. Ripley BD. The second-order analysis of stationary point processes. J Appl Probab. 1976; 13: 255–266.
54. Urban DL. Spatial analysis in ecology Duke’ Nicholas School of the Environment. Levine Science Research Center, Duke University, Durham; 2000.
55. Clarke KR, Somerfield PJ, Gorley RN. Testing of null hypotheses in exploratory community analyses: similarity profiles and biota-environment linkage. J Exp Mar Biol Ecol. 2008; 366: 56–69.
56. Abràmoff MD, Magalhães PJ, Ram SJ. Image processing with ImageJ. Biophotonics Int. 2004; 11: 36–43.
57. Bendel RB, Higgins SS, Teberg JE, Pyke DA. Comparison of skewness coefficient, coefficient of variation and Gini coefficient as inequality measures within populations. Oecol. 1989; 78: 394–400. [PubMed]
58. Soong K. Colony size as a species character in massive reef corals. Coral Reefs. 1993; 12: 77–83.
59. Penin L, Michonneau F, Baird AH, Connolly SR, Pratchett MS, Kayal M, et al. Early post-settlement mortality and the structure of coral assemblages. Mar Ecol Prog Ser. 2010; 408: 55–64.
60. Sandin SA, Smith JE, DeMartini EE, Dinsdale EA, Donner SD, Friedlander AM, et al. Baselines and degradation of coral reefs in the Northern Line Islands. PLoS One. 2008; 3: e1548 doi: 10.1371/journal.pone.0001548 [PMC free article] [PubMed]
61. Harrison PL, Wallace CC. Reproduction, dispersal and recruitment of scleractinian corals In: Dubinsky ed. Coral reefs. Elsevier, Amsterdam; 1990. pp. 133–207.
62. Vermeij MJA, Frade PR, Jacinto RIR, Debrot AO, Bak RPM. Effects of reproductive mode on habitat-related differences in the population structure of eight Caribbean coral species. Mar Ecol Prog Ser. 2007; 351: 91–102.
63. Roth MS, Knowlton N. Distribution, abundance, and microhabitat characterization of small juvenile corals at Palmyra Atoll. Mar Ecol Prog Ser. 2009; 376: 133–142.
64. Doropoulos C, Ward S, Roff G, González-Rivero M, Mumby PJ. Linking demographic processes of juvenile corals to benthic recovery trajectories in two common reef habitats. PLoS One. 2015; 10: e0128535 doi: 10.1371/journal.pone.0128535 [PMC free article] [PubMed]
65. Adjeroud M, Mauguit Q, Penin L. The size-structure of corals with contrasting life-histories: A multi-scale analysis across environmental conditions. Mar Environ Res. 2015; 112: 131–139. doi: 10.1016/j.marenvres.2015.10.004 [PubMed]
66. Adjeroud M, Pratchett MS, Kospartov MC, Lejeusne C, Penin L. Small-scale variability in the size structure of scleractinian corals around Moorea, French Polynesia: patterns across depths and locations. Hydrobiologia. 2007; 589: 117–126.
67. Chiappone M, Sullivan KM. Distribution, abundance and species composition of juvenile scleractinian corals in the Florida reef tract. B Mar Sci. 1996; 58: 555–569.
68. Penin L, Adjeroud M, Pratchett MS, Hughes TP. Spatial distribution of juvenile and adult corals around Moorea (French Polynesia): Implications for population regulation. B Mar Sci. 2007; 80: 379–390.
69. Penin L, Adjeroud M. Relative importance of recruitment and post-settlement processes in the maintenance of coral assemblages in an insular, fragmented reef system. Mar Ecol Prog Ser. 2013; 473: 149–162.
70. Veron JEN. Corals of the world. Australian Institute of Marine Sciences; Townsville, Australia; 2000.
71. Todd PA. Morphological plasticity in scleractinian corals. Biol Rev Camb Philos Soc. 2008; 83: 315–337. [PubMed]
72. Edmunds PJ. Patterns in the distribution of juvenile corals and coral reef community structure in St. John, US Virgin Islands. Mar Ecol Prog Ser. 2000; 202: 113–124.
73. Vermeij MJA, Sandin SA. Density-dependent settlement and mortality structure the earliest life phases of a coral population. Ecology. 2008; 89: 1994–2004. [PubMed]
74. Price N. Habitat selection, facilitation, and biotic settlement cues affect distribution and performance of coral recruits in French Polynesia. Oecol. 2010; 163: 747–758. [PMC free article] [PubMed]
75. Highsmith RC. Reproduction by fragmentation in corals. Mar Ecol Prog Ser. 1982; 7: 207–226.
76. Hughes TP, Connell JH. Multiple stressors on coral reefs: a long-term perspective. Limnol Oceanogr. 1999; 44: 932–940.
77. Dollar SJ. Wave stress and coral community structure in Hawaii. Coral Reefs. 1982; 1: 71–81.
78. Storlazzi CD, Field ME, Dykes JD, Jokiel PL, Brown E. Wave control on reef morphology and coral distribution: Molokai, Hawaii. In: WAVES 2001 conference proceedings. 2002; 1: 784–793.
79. Kayal M, Vercelloni J, Lison de Loma T, Bosserelle P, Chancerelle Y, Geoffroy S, 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. 2012; 7: e47363 doi: 10.1371/journal.pone.0047363 [PMC free article] [PubMed]
80. Wallace CC. Reproduction, recruitment and fragmentation in nine sympatric species of the coral genus Acropora. Mar Biol. 1985; 88: 217−233.
81. Coma R, Lasker HR. Effects of spatial distribution and reproductive biology on in situ fertilization rates of a broadcast-spawning invertebrate. Biol Bull. 1997; 193: 20–29. [PubMed]
82. Hughes TP, Baird AH, Dinsdale EA, Moltschaniwskyj NA, Pratchett MS, Tanner JE, et al. Supply-side ecology works both ways: the link between benthic adults, fecundity, and larval recruits. Ecology. 2000; 81: 2241–2249.
83. Galzin R, Pointier JP. Moorea Island, Society Archipelago. In Proceedings of the Fifth International Coral Reef Symposium. 1985; 1: 73–102.
84. Richmond RH. Coral reefs: present problems and future concerns resulting from anthropogenic disturbance. Am Zool. 1993; 33: 524–536.
85. Hughes TP, Rodrigues MJ, Bellwood DR, Ceccarelli D, Hoegh-Guldberg O, McCook L, et al. Phase shifts, herbivory and the resilience of coral reefs to climate change. Curr Biol. 2007; 17: 360–365. doi: 10.1016/j.cub.2006.12.049 [PubMed]
86. Foster NL, Box SJ, Mumby PJ. Competitive effects of macroalgae on the fecundity of the reef-building coral Montastraea annularis. Mar Ecol Prog Ser. 2008; 367: 143–152.
87. Kuffner IB, Walters LJ, Becerro MA, Paul VJ, Ritson-Williams R, Beach KS. Inhibition of coral recruitment by macroalgae and cyanobacteria. Mar Ecol Prog Ser. 2006; 323: 107–117.
88. Gilmour J. Experimental investigation into the effects of suspended sediment on fertilisation, larval survival and settlement in a scleractinian coral. Mar Biol. 1999; 135: 451–462.
89. Box SJ, Mumby PJ. Effect of macroalgal competition on growth and survival of juvenile Caribbean corals. Mar Ecol Prog Ser. 2007; 342: 139–149.
90. Fichez R, Adjeroud M, Bozec Y-M, Breau L, Chancerelle Y, Chevillon C, et al. A review of selected indicators of particle, nutrient and metal inputs in coral reef lagoon systems. Aquat Living Resour. 2005; 18: 125–147.
91. Cooper TF, Uthicke S, Humphrey C, Fabricius KE. Gradients in water column nutrients, sediment parameters, irradiance and coral reef development in the Whitsunday Region, central Great Barrier Reef. Estuar Coast Shelf Sci. 2007; 74: 458–470.
92. Brown D, Edmunds PJ. Long-term changes in the population dynamics of the Caribbean hydrocoral Millepora spp. J Exp Mar Biol Ecol. 2013; 441: 62–70.

Articles from PLoS ONE are provided here courtesy of Public Library of Science