A positive trajectory for corals at Little Cayman Island
By Carrie Manfrino From PLOS ONE
Abstract
Coral reefs are damaged by natural disturbances and local and global anthropogenic stresses. As stresses intensify, so do debates about whether reefs will recover after significant damage. True headway in this debate requires documented temporal trajectories for coral assemblages subjected to various combinations of stresses; therefore, we report relevant changes in coral assemblages at Little Cayman Island. Between 1999 and 2012, spatiotemporal patterns in cover, densities of juveniles and size structure of assemblages were documented inside and outside marine protected areas using transects, quadrats and measurements of maximum diameters. Over five years, bleaching and disease caused live cover to decrease from 26% to 14%, with full recovery seven years later. Juvenile densities varied, reaching a maximum in 2010. Both
Introduction
Coral reefs occupy less than 0.01% of the marine environment, yet they harbor up to 25% of marine biodiversity, yield approximately 25% of the fish catch in developing nations, and generate up to 30% of export earnings in 100 countries that promote reef-related tourism [1]. These benefits may disappear because coral reefs around the world are being degraded by local and global anthropogenic stressors that cause damage beyond that due to natural disturbances, such as large storms, hurricanes, exposure to unusually low tides, and freshwater inputs [2]. Local anthropogenic
Such statistics have engendered debates about the sustainable management of coral reefs, including the value of marine protected areas with or without no-take zones [7], [8]. Managed areas of any type, even those with stringent enforcement, provide no direct protection from natural cataclysms and global perturbations [4], [7], [9], [10]. Nevertheless, managed areas of sufficient size and connectivity may promote resilience through bottom-up effects, e.g., by providing refuges for reproducing corals and enhancing regional recruitment through larval exchange [4]. Moreover, protection from overfishing in no-take zones and protected areas that connect multiple habitats may generate beneficial top-down effects by fostering robust fish and invertebrate assemblages that include herbivores grazing on algae that could otherwise usurp the open space needed by coral larvae or regenerating fragments [4], [11]. Regardless of the mechanism, understanding the conditions that promote recovery of damaged coral reefs remains critical to the formulation of effective management actions.
Studies documenting recovery of corals have yielded mixed conclusions. Fine-scale surveys and experimental studies indicate diverse responses or the lack of a short-term response, which may be due to i) local oceanographic, meteorological and ecological conditions that ameliorate or exacerbate stress from disturbances; ii) variable growth and regeneration rates among coral taxa; iii) variation in larval supply; iv) differences in type and extent of the most recent disturbance; and v) unique interactions among coincident and sequential disturbances [12]–[21]. Data analyzed in broad-scale meta-analyses also suffer from the effects of these influences, along with biases from the non-random placement of managed areas and significant variation in enforcement [22]–[24]. Although two meta-analyses report little or slow recovery of corals in managed areas [22], [23], one study describes increased coral cover within protected areas and decreased cover on unprotected reefs [24]. In contrast, a recent study ascribes recovery from major disturbance to minimal local anthropogenic stresses and long distance recruitment rather than the presence of a protected area [25]. Given these inconsistencies, there remains a need for long-term studies employing sampling designs that document conditions within and outside of managed areas [14], [23], [24]. Eventually, a synthesis of relevant results will deliver insights required for sustainable management of coral reefs and the values they deliver.
As a contribution to this goal, this long-term study examines temporal trajectories for coral reefs in protected and unprotected areas off Little Cayman Island. All data were derived from nondestructive observations that were in full compliance with conditions set by the Cayman Islands Marine Conservation Board. Data collected during 10 of the last 14 years provide the basis for testing null hypotheses related to the interactions between marine protected areas and the effects of bleaching and disease, i.e., bleaching and disease led to no significant temporal changes or spatial differences in i) cover of live coral, ii) abundance of juvenile corals, and iii) taxonomic composition and size structure of live coral assemblages within and outside of marine protected areas.
Materials and Methods
Little Cayman, Cayman Islands is a 17 × 2 km, flat, low-lying island that lies 120 km northeast of Grand Cayman, 145 km south of Cuba and 10 km southwest of Cayman Brac (Fig. 1). Live coral cover around the island has been documented to be equivalent to or higher than other locations in the Caribbean [26]. Two key factors contribute to this status, i.e., the lack of freshwater rivers on the island and minimal anthropogenic stress due to nutrient inputs, commercial fishing, groundings or other activities undertaken by 5% of Agaricia spp., Diploria spp., Montastraea spp., Porites spp. or Siderastrea spp. showed signs of bleaching in 2009 [28].
Figure 1. Maps showing the Caribbean region (inset) and Little Cayman Island.
Marine Protected Areas shaded green and buoys used in sampling indicated by •. Coordinates for buoys: 1 = 19°39′18″N, 80°01′24″W; 2 = 19°40′17″N, 80°06′12″W; 3 = 19°41′06″N, 80°04′41″W; 4 = 19°41′33″N, 80°04′07″W; 5 = 19°41′40″N, 80°04′10″W; 6 = 19°42′03″N, 80°03′25″W; 7 = 19°42′25″N, 80°00′44″W; 8 = 19°40′51″N, 80°01′24″W; 9 = 19°39′26″N, 80°05′22″W.
doi:10.1371/journal.pone.0075432.g001
In July and early August, shallow (9–13 m) spur and groove reefs around Little Cayman were surveyed in 1999, annually from 2002 to 2007, and annually from 2009 to 2012 [9]. Surveys were based on the Atlantic and Gulf Rapid Reef Assessment techniques that are designed to document broad-scale and long-term status and trends [29]. Surveys employed haphazard, 10-m line intercept transects along spurs. To ensure dispersion and more rigorous representation of island-wide conditions, transects were distributed around multiple, permanent mooring buoys on either the leeward or windward side of the island. Across all years of sampling, coral assemblages were surveyed along 8–37 transects inside marine protected areas and 18–61 transects outside marine protected areas as dictated by weather and logistics.
Cover of live corals was analyzed for 9 years when ≥10 transects were surveyed at ≥2 mooring buoys both inside and outside marine protected areas (Table 1). The distances beneath each transect line occupied by live coral were measured to the nearest 1 cm. The total distance along each replicate transect occupied by live coral was converted to proportional cover for the transect by dividing it by the total length of the transect minus the distance occupied by sandy substrate unsuitable for corals.
Coral recruitment was estimated by counting colonies that were ≤2 cm in diameter. These colonies were counted in 625-cm2 quadrats placed every 2 m along the 10-m transects. For 10 years when ≥85 quadrats were sampled at ≥2 mooring buoys both within and outside marine protected areas, counts were summed to the transect level (Table 1). The resulting, replicate sums were scaled to numbers m−2.
To document the contribution of coral skeletons to the structure of reefs, coral colonies under each transect line were identified to the lowest possible taxonomic level, and the maximum diameters of their skeletons were measured. In 10 years when ≥50 colonies were sampled at ≥2 mooring buoys both within and outside marine protected areas, size frequency distributions were created by pooling data from transects surveyed around each mooring buoy (Table 1). The analyses focused on 26 species that represented the common corals in the region. Size classes ranged to ≥190 cm in 10-cm intervals, with the upper size class selected so that at least one species had 10 or more colonies in this bin. Counts per m for all size classes were scaled to a 100-m, linear search area to account for differing numbers of 10-m transects.
Live coral cover and numbers of juveniles m−2 were analyzed with a general linear model that treated year of sampling and degree of protection treated as orthogonal, fixed factors. Proportions for live coral cover were arcsin transformed prior to analysis, Type III sums of squares were used to account for unbalanced replication, and Cochran’s and Anderson–Darling tests were applied to residuals to evaluate compliance with the assumptions of homoscedasticity and normality, respectively.
Scaled counts for combinations of species and size classes were analyzed with a multivariate permutation analysis of variance (PERMANOVA) with years of sampling and levels of protection as orthogonal, fixed factors [30]. Thus, the PERMANOVA assessed significant differences among years of sampling and between levels of protection for both taxonomic composition and size frequency distributions of live coral assemblages, with the latter metric providing insights into differences or changes in the structure of the reef. A permutation analysis of multivariate dispersion (PERMDISP; [30]) was used to determine if any significant difference was due to increased variation among samples rather than a difference in the mean numbers of colonies in combinations of size classes and taxa. Finally, a similarities percentage routine (SIMPER; [30]) identified the primary combinations of size classes and taxa contributing to any significant difference detected by the PERMANOVA.
Results
The proportional cover of live coral along line intercept transects exhibited statistically similar temporal trajectories in protected and unprotected areas (Table 2). Transformation did not yield homoscedasticity, but the results appeared reliable. The lack of a significant interaction and significant difference between protected and unprotected areas certainly represented robust results due to the unidirectional bias introduced by heteroscedasticity. The significant variation among years also appeared informative due to the level of significance and results of Dunnett’s multiple comparisons with 1999 treated as a control. Relative to 1999 levels, mean coral cover decreased significantly from 2002 to 2004 (26% to 14%), remained significantly lower through 2010, returned to a level that was not significantly different in 2011, and remained at that level in 2012 (Fig. 2). In both 2010 and 2011, mean cover of live coral increased by 5%. Thus, coral cover increased despite a localized bleaching event recorded in 2009 [28].
Densities of coral colonies ≤2 cm in diameter also exhibited statistically similar temporal trajectories in protected and unprotected areas (Table 2). According to Tukey’s tests, numbers of juveniles m−2 in 2010 were significantly higher than densities in all other years, and densities recorded in 2011 also were relatively high, albeit not statistically different from densities recorded in 2004 and 2005 (Fig. 3). Although both live coral cover and densities of juveniles increased contemporaneously, densities were not correlated significantly with coral cover in the following year (Pearson correlation coefficient = 0.689, p = 0.06).
In total, size classes for 26 species yielded 186 categories. The numbers of coral colonies in these categories differed significantly among years and between protected and unprotected areas (Table 3). Only the differences among years reflected increased variability among replicate samples as shown by the results for the permutation analyses of dispersion (PERMDISP; Table 4). Notably, coral assemblages displayed consistent temporal trajectories within and outside protected areas and consistent spatial patterns across years as evidenced by the nonsignificant interaction term (Table 3).
According to the SIMPER analysis for the 45 possible pairs of years, 41 to 79 combinations of size class and coral species accounted for 90% of the differences, with results for two size classes of Montastraea annularis and Porites astreoides illustrating key patterns (Fig. 4). Larger standard errors for means from some years illustrated the increased variability detected by the permutation analysis of dispersion (Table 4). Nevertheless, clear patterns were discernible. In 2010, relatively large numbers of
Out of the 186 combinations of size classes and coral taxa, a SIMPER analysis indicated that 75 accounted for 90% of the difference in the single comparison of assemblages within and outside marine protected areas. Numbers of colonies in different size classes for M. annularis, M. cavernosa, M. faveolata, M. franksi, P. astreoides, and S. siderea illustrated key patterns (Fig. 5). In this case, the consistency of standard errors reinforced the lack of a significant result for the permutation analysis of dispersion (Table 4). Means did vary between protected and unprotected areas, but the differences were
Discussion
Surveys spanning 14 years documented several ecologically significant results for coral assemblages off Little Cayman Island. Live coral cover decreased from 26% to 14% following thermal stress, bleaching and disease [9], [27], [28], but cover recovered to previous levels from 2010 onward. In 2010, significantly higher densities of colonies ≤2 cm in diameter provided evidence of a relatively large recruitment event. Minor changes in assemblage structure were detected via analyses of size classes of 26 species of coral, with decreases in some size classes documented after 2002 for species of Montastraea and M. meandrites, but not for other corals. In addition, only minor differences in assemblages within and outside marine protected areas were detected, with some size classes of some species more common within protected areas and others more common outside protected areas. Notably, temporal trajectories for all metrics did not differ significantly between locations within and outside marine protected areas, i.e., there were no significant interactions in the statistical analyses.
During recovery of live coral cover off Little Cayman Island, the increase from 14% to 20% between 2009 and 2010 and the increase from 20% to 25% between 2010 and 2011 were slightly less than the median increase in cover of 8% y−1 derived from previous reports of colony growth or colony recovery after a variety of disturbances on 67 reefs with <1% to 74% coral cover [13], [16], [17], [21], [31]–[33]. Thus, recovery of live coral cover at Little Cayman Island was on par with recoveries reported from other locations, and it occurred at the same rate within and outside protected areas.
In combination, an increase in densities of juvenile colonies and increases in small colonies of 13 coral species highlighted a relatively strong pulse of recruitment. Throughout the 14 years, reefs at Little Cayman Island received reasonable numbers of recruits (2–12 m−2). In fact, these densities of coral colonies ≤2 cm in diameter were i) higher than records from the Great Barrier Reef where colonies ≤2 cm in diameter represented juveniles (0.1–0.8 m−2; [12]); ii) similar to values from the Northern Line Islands where juveniles were designated as 1–5 cm colonies (1–10 m−2; [19]); iii) lower than maximum densities observed in Jamaica where juveniles were classified as 2–4 cm colonies (1–212 m−2; [34], [35]); and iv) lower than maximum densities observed in Curaçao, Bonaire, St. John in the U.S. Virgin Islands, the Florida Keys and Belize where colonies ≤5 cm were considered juveniles (1–44 m−2; [36]–[38]). Numbers of juveniles were not significantly correlated with increases in live coral cover during the following year, which suggested the need for further investigation to determine how variation in prevailing oceanographic conditions affects recruitment and how recruitment relates to recovery of live coral cover [17].
Overall, the size frequency distributions of hard coral assemblages in protected and unprotected areas at Little Cayman Island remained relatively consistent, with 87% of the maximum differences among paired years being ≤2 colonies per 100-m, linear search area. Numbers of larger colonies of M. annularis, M. cavernosa, M. faveolata and M. franksi and M. meandrites decreased in years following bleaching events and disease outbreaks [9], but numbers did not continue to decline as reported for these species at other Caribbean locations [39], [40]. In addition, numbers of P. astreoides, A. palmata, C. natans, M. alcicornis, and S. siderea did not decline consistently during the years surveyed, which suggested resistance to bleaching and disease [32], [41], [42]. At Little Cayman Island, P. astreoides and S. siderea were reported to suffer the lowest rates of bleaching and mortality [9], [28], but the same species, congeners and A. palmata have been reported to suffer bleaching and mortality elsewhere [42]–[49]. Perhaps, local oceanographic conditions, including adjacent deep water, ameliorated the effects of widespread thermal stress [9], [27], [28]. Colony size, physiological condition and genetic variation also influence resistance to stress, with greater tolerance documented for some coral species, larger colonies and populations exposed to repeated stress [41], [42]. In fact, some larger colonies of massive species, e.g., Montastraea species, persisted at Little Cayman Island. In contrast, the relative absence of local stressors at Little Cayman Island would not enhance resistance. Nevertheless, local adaptation by corals can be promoted if phylogeographical barriers reduce gene flow across a species’ range, as reported for Caribbean populations of A. palmata [50]. Further work is needed to disentangle the influences of genetics, physiology and environmental conditions on the persistence of coral colonies.
This study of coral reefs off Little Cayman Island demonstrated recovery after disease combined with regional thermal stress to cause a decrease in live coral cover from 26% to 14% over a five-year period. The temporal trajectories of the decline and subsequent recovery were similar within and outside protected areas, and the assemblage composition also remained similar across protected and unprotected sites throughout the decline and recovery. In addition, juvenile corals achieved similar densities within and outside protected areas. Key factors shaping recovery of corals off Little Cayman probably included the isolated geographic setting; stringent protection of a significant portion of the reefs resulting in healthy populations of herbivorous fishes and preservation of key trophic links; and minimal stress from local human activities.
Although any documented recovery of coral is encouraging, it is unlikely that such positive effects will spread throughout the Caribbean unless protection from local stresses is improved [51]–[53]. Without such improvements, recovery from natural cataclysms, including those exacerbated by global change, remains unlikely. Ultimately, management of local and global stresses will be required to sustain coral reefs and ensure their capacity to recover from disturbance [1].
Acknowledgments
We thank V. Coelho, M. Brandt, C. Dell, K. Wilson, J. Moye, A. Little, S. Barry; K. Lohr and numerous interns and staff from the Little Cayman Research Centre (LCRC) for field assistance and logistical support. This project was part of the Long-Term Assessment and Monitoring Program at LCRC. Metadata from this study will be made available through the Global Coral Reef Monitoring Network of the International Union for the Conservation of Nature and the Atlantic and Gulf Rapid Reef Assessment, and the corresponding author will respond to requests for data. This paper is contribution 2013:2 of the Central Caribbean Marine Institute.
Citation: Manfrino C, Jacoby CA, Camp E, Frazer TK (2013) A Positive Trajectory for Corals at Little Cayman Island. PLoS ONE 8(10): e75432. doi:10.1371/journal.pone.0075432
Editor: Senjie Lin, University of Connecticut, United States of America
Received: February 14, 2013; Accepted: August 14, 2013; Published: October 9, 2013
Copyright: © 2013 Manfrino et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: These authors have no support or funding to report.
Competing interests: The authors have declared that no competing interests exist.
For more on this article with PP presentation material and additional images go to:
http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0075432