Hostname: page-component-7b9c58cd5d-dkgms Total loading time: 0 Render date: 2025-03-16T07:38:08.194Z Has data issue: false hasContentIssue false

Higher population densities of the sea urchin Diadema antillarum linked to wave sheltered areas in north Puerto Rico Archipelago

Published online by Cambridge University Press:  20 June 2014

Ruber Rodríguez-Barreras*
Affiliation:
Department of Biology, University of Puerto Rico, Río Piedras, PO Box 23360 San Juan, Puerto Rico 00931-3360
María E. Pérez
Affiliation:
Department of Mathematics, and Center of Applied Tropical and Conservation, University of Puerto Rico, Río Piedras, PO Box 23355 San Juan, Puerto Rico 00931-335
Alex E. Mercado-Molina
Affiliation:
Department of Biology, University of Puerto Rico, Río Piedras, PO Box 23360 San Juan, Puerto Rico 00931-3360
Stacey M. Williams
Affiliation:
Institute for Social-Ecological Research, PO Box 3151, Lajas, Puerto Rico 00667
Alberto M. Sabat
Affiliation:
Department of Biology, University of Puerto Rico, Río Piedras, PO Box 23360 San Juan, Puerto Rico 00931-3360
*
Correspondence should be addressed to: R. Rodríguez-Barreras, Department of Biology, University of Puerto Rico, Rio Piedras, PO Box 23360 San Juan, Puerto Rico 00931-3360 email: ruber.rodriguez@outlook.com
Rights & Permissions [Opens in a new window]

Abstract

The long-spined sea urchin Diadema antillarum has been the focus of multiple studies since the mass mortality event in the 1980s. The recovery of this key herbivore in the wider Caribbean is essential for the well-being of coral reefs. This study examined the population density and structure of D. antillarum at seven northern fringing reefs of Puerto Rico between 2011 and 2013. The total mean density of the sea urchins in northern Puerto Rico was 0.9 ±0.3 ind m−2. Densities of D. antillarum significantly differed among sites, but not temporally. Differences in mean sizes were significant among sites and seasons. Areas with higher densities of D. antillarum showed lower cover of non-calcareous algae. Wave exposure was correlated with the abundance of the sea urchin. This study indicates that the observed abundance of D. antillarum has not yet returned to pre-mortality levels. However, densities showed some degree of recovery when compared with previous studies, enabling at least some degree of control on fleshy macroalgae communities. No significant changes in density occurred between 2011 and 2013, and sites with higher densities were generally located in leeward areas. The low relative abundance of small size individuals points towards recruitment limitation as an explanation for the limited recovery of D. antillarum.

Type
Research Article
Copyright
Copyright © Marine Biological Association of the United Kingdom 2014 

INTRODUCTION

Over the past five decades coral reefs in the Caribbean have changed dramatically (Gardner et al., Reference Gardner, Côté, Gill, Grant and Watkinson2003; Hughes et al., Reference Hughes, Graham, Jackson, Mumby and Steneck2010; De'ath et al., Reference De'ath, Fabricius, Sweatman and Puotinen2012). Regional declines in coral cover have been asynchronous and due to cumulative factors, such as, hurricanes, disease outbreaks, bleaching, pollution and overfishing (Bythell et al., Reference Bythell, Gladfelter and Bythell1993; Hughes, Reference Hughes1994; Bythell et al., Reference Bythell, Hillis-Starr and Rogers2000; Kramer et al., Reference Kramer, Kramer and Ginsburg2003; Hughes et al., Reference Hughes, Graham, Jackson, Mumby and Steneck2010). One of the most dramatic shifts in community structure occurred after the massive die-off of Diadema antillarum (Echinodermata, Echinoidea). This species is considered a keystone herbivore in the western Atlantic, exerting top-down control on algal abundance (Carpenter & Edmunds, Reference Carpenter and Edmunds2006). The long-spined sea urchin has been the focus of many studies after its epizootic event between 1983 and 1984 (Lessios et al., Reference Lessios, Cubit, Robertson, Shulman, Parker, Garrity and Levings1984), when their mortality reached up to 99.9% at some Caribbean locations (Hughes et al., Reference Hughes, Keller, Jackson and Boyle1985; Hunte et al., Reference Hunte, Côté and Tomascik1986).

Herbivores maintain the stability and diversity of the highly-productive coral-dominated state of tropical reefs (Sammarco et al., Reference Sammarco, Levinton and Ogden1974; Sammarco, Reference Sammarco1982; Littler et al., Reference Littler, Littler and Titlyanov1991; Hughes, Reference Hughes1994; Carpenter & Edmunds, Reference Carpenter and Edmunds2006). Diadema antillarum and parrot-fish have the greatest impact on algal abundance and the widest feeding range of algal types (Steneck Reference Steneck1988). Profound changes took place on reefs after the collapse of D. antillarum populations (Hughes et al., Reference Hughes, Keller, Jackson and Boyle1985), shifting the structure and function of Caribbean coral reef ecosystems (Liddell & Ohlhorst, Reference Liddell and Ohlhorst1986). In some Caribbean reef locations, the benthic algal cover increased between 100% and 250% (Phinney et al., Reference Phinney, Muller-Karger, Dustan and Sobel2001). Non-calcareous algae of the genera Gelidium, Coelothrix and Amphiroa became more abundant in Puerto Rico (Vicente, Reference Vicente1987). The potential role of the sea urchin as the agent facilitating the switch between the current algae-dominated state in many Caribbean reefs to the original coral dominated one is considered a ‘new opportunity’ for coral reefs (Knowlton, Reference Knowlton2001), given the demise of herbivorous fish due to overfishing (Idjadi et al., Reference Idjadi, Haring and Precht2010).

Pre-mortality data available for Puerto Rico are useful to understand the dynamics in the recovery process locally (Craft, Reference Craft1975; Rivera & Vicente, Reference Rivera and Vicente1976; Bauer, Reference Bauer1980; Vicente & Goenaga, Reference Vicente and Goenaga1984). Densities dropped from 12.7–13.8 ind m−2 (Craft, Reference Craft1975; Bauer, Reference Bauer1980) to less than 1.0 ind m−2 (Ruiz-Ramos et al., Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011; Soto-Santiago & Irizarry-Soto, Reference Soto-Santiago and Irizarry-Soto2013). Presently, D. antillarum is going through a slow and non-uniform recovery process in many Caribbean reefs where populations still exhibit low densities (Bak et al., Reference Bak, Carpay and Ruyter van Steveninck1984; Vicente & Goenaga, Reference Vicente and Goenaga1984; Hughes et al., Reference Hughes, Keller, Jackson and Boyle1985; Hunte et al., Reference Hunte, Côté and Tomascik1986; Miller et al., Reference Miller, Adams, Ogden, Ogden and Ebersole2003; Lessios, Reference Lessios2005; Weil et al., Reference Weil, Torres and Ashton2005; Debrot & Nagelkerken, Reference Debrot and Nagelkerken2006; Noriega et al., Reference Noriega, Pauls and Mónaco2006; Steiner & Williams. Reference Steiner and Williams2006; Sellers et al., Reference Sellers, Casey, Burge and Koepfler2009; Idjadi et al., Reference Idjadi, Haring and Precht2010; Martín-Blanco et al., Reference Martín-Blanco, González-Sansón, Pina-Amargós and Clero-Alonso2010, Reference Martín-Blanco, Clero-Alonso, González Sansón and Pina-Amargós2011; Ruiz-Ramos et al., Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011; Levitan et al., Reference Levitan, Edmunds and Levitan2014). Differences in recovery may be explained by physical factors such as wave exposure (Chollett et al., Reference Chollett, Mumby, Müller-Karger and Hu2012). Higher densities of the species have been linked to wave sheltered habitats (Clemente & Hernández, Reference Clemente and Hernández2008; Hernández et al., Reference Hernández, Clemente, Sangil and Brito2008). Leeward areas seems to favour higher densities of Diadema and seems to be an important oceanographic factor modeling population dynamics (Debrot & Nagelkerken, Reference Debrot and Nagelkerken2006). The effect of unidirectional hydrodynamic forces on D. antillarum could be explained because its morphological features do not allow a large adhesive area to attach, being considered weakly resistant to wave exposure (Tuya et al., Reference Tuya, Cisneros-Aguirre, Ortega-Borges and Haroun2007).

On the other hand, the recovery of D. antillarum populations at some Caribbean sites has been related to reductions in macroalgae cover and increase of the recruitment for reef-building corals (Bellwood et al., Reference Bellwood, Hughes, Folke and Nystrom2004; Carpenter & Edmunds, Reference Carpenter and Edmunds2006; Idjadi et al., Reference Idjadi, Haring and Precht2010; Soto-Santiago and Irizarry-Soto, Reference Soto-Santiago and Irizarry-Soto2013). It is crucial to assess the current status of D. antillarum populations a decade after the last available data in order to monitor the population dynamics of this important herbivore. The main objectives of this study were to: (1) determine the current density of D. antillarum and compare our results with previous studies to determine whether D. antillarum populations have recovered; (2) analyse spatial and temporal changes in density and size structure; and (3) explore the relationship between the abundance of D. antillarum and wave exposure and in relation to some of the main reef functional groups (calcareous algae, fleshy macroalgae, turf and corals).

MATERIALS AND METHODS

Study area

Surveys were conducted between August 2011 and August 2013 at seven shallow-water fringing reefs (1–3 m deep). The sites were Cerro Gordo A (CGA—18°16′51.40″N 65°17′12.21″W), Cerro Gordo B (CGB—18°28′53.78″N 66°20′25.12″W), Isla Verde (IVD—18°26′49.79″N 66°0′49.88″W), Luquillo (LQY—18°23′18.46″N 65°43′5.52″W), Tamarindo 1 (TM1—18°18′55.29″N 65°19′5.82″W), Tamarindo 2 (TM2—18°19′0.78″N 65°19′2.60″W) and Melones (MLN—18°18′15.70″N 65°18′42.27″W) (Figure 1). The first four sites are characterized by a very narrow shelf and high energy sandy beaches, due to the effect of north-east trade winds and North Atlantic winter storms. Due to the high annual precipitation levels and the discharge of rivers (Williams et al., Reference Williams, Benavides-Serrato, García-Arrarás, Hernández-Delgado, Rodríguez-Barreras, Alvarado-Barrientos and Solis-Marin2013), high sediment loads are common in CGA, CGB, IVD and LQY. On the other hand, TM1, TM2, and MLN are located in Culebra Island. This island has no rivers, precipitation is low and it is characterized by volcanic and intrusive rocks, with some limestone deposits, whereas sites along the north coast (CGA, CGB, IVD and LQY) are made up by carbonate rocks (USGS, Reference Veve and Taggart1996).

Fig. 1. Study area and sampling sites in Puerto Rico. Cerro Gordo A (CGA), Cerro Gordo B (CGB), Isla Verde (IVD), Luquillo (LQY), Tamarindo 1 (TM1), Tamarindo 2 (TM2) and Melones (MLN).

Survey methodology

Mean densities were estimated using a belt-transect methodology (Sellers et al., Reference Sellers, Casey, Burge and Koepfler2009). At each site, eight fixed transects of 10 m2 (5 m × 2 m) were placed parallel to the coast and separated 10 m from each other. All individuals within transects were counted. All crevices and small holes were carefully inspected to avoid missing any individual. Field sampling was carried out every 6 months. We selected August and February based on differences in local precipitation and temperature (NOAA Coral Reef Watch, 2000), because both abiotic factors have some influence on echinoid population dynamics (Lugo-Ascorbe, Reference Lugo-Ascorbe2004; Williams et al., Reference Williams, García-Sais and Capella2009). With the use of a caliper, we measured the test diameter of fifty sea urchins, which were collected randomly every season by site (caliper error ±0.05 mm). We divided populations in three size-categories: small, <40.0 mm; medium, 41.0–60.0 mm; and large, >60.0 mm following the method of Miller et al. (Reference Miller, Adams, Ogden, Ogden and Ebersole2003).

Substrate complexity, cover and wind analysis

The rope-and-chain method was used to measure the substrate complexity at all sites. It was expressed as the ratio between the total length of a chain and the length of the same chain when molded to the reef surface (complexity index = total length of the chain/total distance when moulding the bottom). The range varied from 1 (flat surface) to higher values indicating more complex topography (Alvarez-Filip et al., Reference Alvarez-Filip, Dulvy, Gill, Cote and Watkinson2009). Five PVC tube quadrats of 25 cm2 per transect were photographed and analysed using CPCe-4.1 (Kohler & Gill, Reference Kohler and Gill2006). We selected 25 points randomly per picture and estimated percentage cover of the following categories: fleshy macroalgae (NCAL), calcareous algae (CAL), turf and live coral (LCOR), using only data from August 2012. Furthermore, we analysed the relationship between the abundance of D. antillarum and the wave exposure based on the chronic stress map given by Chollett et al. (Reference Chollett, Mumby, Müller-Karger and Hu2012). Wave exposure values varied on a scale from 1 (low) to 9 (high). Sites in Culebra (TM1, TM2 and MLN) were under moderate wave action (5.5), LQY occupied an intermediate position with 7.0, whereas IVD, CGA, and CGB were the most affected sites with 7.9. We re-evaluated CGB because Chollet et al. (Reference Chollett, Mumby, Müller-Karger and Hu2012) did not reflect the existence of a natural protection against wave exposure. Culebra sites have a similar protection against high wave action due to the presence of Luis Peña Key in front of these sites. Because of this, we decreased the category of this site from 7.9 to 5.5.

Statistical analysis

We fitted a two-way repeated measures analysis of variance (ANOVA) model to detect differences in mean density between sites. The model incorporates time as fixed factor, and transect (fixed factor) nested inside site (random factor). Mean size was analysed using a two-way ANOVA model with time as fixed factor, and site as random one. Fixed effects were tested using F tests; models involving different random effects were compared using likelihood ratio tests. Mean density data followed a normal distribution. A power transformation was applied for reaching normality and homogeneity of variance (Box & Cox, Reference Box and Cox1964). Result suggested the square (not the square root) for mean size (suggesting a normal distribution for the area of the transversal section). We ran a post-hoc Bonferroni test for multiple comparisons of means for density and size (Zar, Reference Zar2010). Pearson's correlation test was used to determine the relationship between the mean density of D. antillarum and cover non-calcareous algae (fleshy macroalgae), coral, and substrate complexity, but not turf. Potential differences among small, medium and large categories for each site and season were analysed by fitting log-linear models. Furthermore, we ran a non-parametric Spearman correlation test between wave exposure and the mean density of D. antillarum. All statistical analyses were performed in the free-license statistical software R v.3.0.1, with a pval = 0.05 (R Core Team, 2013). Packages MASS and nlme were used for estimation of power transformations, and for complex repeated measured analysis of variance respectively (Venables & Ripley, Reference Venables and Ripley2002; Pinheiro & Bates, Reference Pinheiro, Bates and DebRoy2013).

RESULTS

Spatial and temporal abundance

Overall mean density (± 95% confidence interval) of D. antillarum in Puerto Rico was 1.148 ±0.345 ind m−2 at the beginning of the study, and 0.946 ±0.295 ind m−2 at the end. The highest density recorded was 1.398 ±0.342 ind m−2 in February 2012, whereas a minimum of 0.946 ±0.295 ind m−2 was in August 2013. Sites were generally stable temporally in mean density, except in CGA and LQY, where densities increased 0.85 and 0.90 ind m−2 from August of 2011 to February of 2012 respectively. However, temporal differences were not detected among sites (Table 1).

Table 1. Two-way repeated measure analysis of variance on the abundance (a) and the two-way analysis of variance test on the size of the sea urchin Diadema antillarum (b). *this value comes from a χ2 test that compares models with and without site as a random effect.

Spatial differences in density were found (Table 1). These differences among sites were related with the low densities reported in IVD and CGB (Figure 2). Differences were found between CGB and CGA (Bonferroni, P = 0.025), and TM2 (Bonferroni, P = 0.004); and also between IVD and CGB (Bonferroni, P = 0.033) and TM2 (Bonferroni, P = 0.006). During this study mean densities at CGB, TM2, and MLN were above 1.0 ind m−2 (Figure 2). At TM2 the maximum peak of abundance was 2.402 ±0.624 ind m−2 during February of 2012, whereas at IVD and CGA reported densities were below 0.50 ind m−2.

Fig. 2. Spatial and temporal density of Diadema antillarum in Puerto Rico from August 2011 to August 2013 at CGA, CGB, IVD, LQY, TM1, TM2 and MLN (see Methods for acronyms). The X axis represents time with A-11 (August 2011), F-12 (February 2012), etc. The Y axis is the mean density of D. antillarum (ind m−2). Bars represent 95% confidence interval of the mean and red line represents a polynomial fit. Last graph shows the overall temporal change of mean density using pooled data.

Size structure

The mean size of D. antillarum fluctuated throughout the study, with maximum values during summer (6.51 ±0.14 cm) and minimums ones during winter (6.23 ±0.15 cm; Figure 3). However, seasonal differences were not significant (Table 1). The mean size peaked at IVD during February of 2013 with 7.91 ±0.77 cm; while TM2 reached the lowest mean size of 4.04 ±1.08 cm during February of 2012. We found spatial and site–time interaction differences on size (Table 1). Spatial differences were related mainly with TM2, the site with the lowest mean size, due to the high abundance of medium and small sizes found in this site (Bonferroni, P < 0.00001), and also among IVD and the other stations (Bonferroni, P < 0.0001; Figure 3).

Fig. 3. Seasonal size-classes distribution of Diadema antillarum at CGA, CGB, IVD, LQY, TM1, TM2 and MLN in Puerto Rico (see Materials and Methods for acronyms). Black bar: small size (<40 mm), white: medium size (41–60 mm) and dark grey: large size (>60 mm). Last graph shows the temporal variation of mean size of D. antillarum of all sites.

Small individuals (<40 mm) were the less common size class, contributing with 14% of 1534 echinoids measured during 2 yr at all sites. The exception was TM2, where we found an unusual population size structure, dominated by medium and small size-classes (Figure 3). Small sea urchins represented up to 44% during August 2011, whereas large sea urchins did not reached 10% of abundance at any time in TM2. When we considered new recruits (<25 mm) as a separate size-class, they were only 1.10% of the overall individuals measured during the whole study at all sites; and the majority of them were reported during August 2011. No juveniles were detected at IVD, CGB and TM1, whereas the two sites with the greatest number of juveniles were TM2 (8%) and CGB (4%) during August of 2012.

The large size-class (>60 mm) was the most abundant in space and time, followed by the medium size (41–60 mm); while the small category (<40 mm) was the scarcest one. The cohort distribution analysis by fitting log-linear models showed that population structures changed according to site, and these changes were affected by the season (P = 2.59 × 10−8). Sites located near each other such as CGA/CGB and TM1/TM2 showed different size structure. Medium and small cohorts were more abundant at CGB and TM2 while large cohort was dominant at CGA and TM1.

Substrate complexity and cover

Substrate complexity varied from 1.06 ±0.02 units at LQY to 1.37 ±0.16 units at MLN. Fleshy macroalgae (NCAL) varied from a minimum of 1.8 ±2.75% in TM2, to a maximum of 43.6 ±13.21% in CGA (Figure 4). Turf showed a similar pattern, with a minimum of 1.45 ±3.69% at TM2 and a maximum of 18.18 ±7.01% at CGA. Calcareous algae were more abundant at sites with both less fleshy macroalgae and turf (Figure 4). The lowest percent of calcareous algae was found at CGA with 1.09 ±1.87%, while TM1 displayed the highest percent with 10.18 ±10.93%. Coral cover was low at CGA and IVD (<10%), but greater than 15% at the other sites, peaking at TM2 (29.09 ±12.15%, Figure 4). We found a significant inverse relationship between the abundance of D. antillarum and the cover of fleshy macroalgae (Pearson, r 2 = 0.909, r = −0.954, P = 0.0009). We found low correlation between the abundance of D. antillarum and coral cover (Pearson, r 2 = 0.305, r = 0.416, P = 0.0093), and no significant relationship between sea urchin abundance and substrate complexity (Pearson, r 2 = 0.541, r = 0.295, P = 0.179). In addition, mean density of D. antillarum was correlated with wave exposure (r2 = 0.519, r = –0.580, P = 0.0003).

Fig. 4. Spatial variation of the percentage cover of fleshy macroalgae (NCAL), turf, calcareous algae (CALG) and live coral (LCOR) at CGA, CGB, IVD, LQY, TM1, TM2 and MLN in Puerto Rico (see Materials and Methods for acronyms). Last graph shows the relationship between density of Diadema antillarum and percentage cover of fleshy macroalgae (NCAL). Box represents 95% confidence interval of the mean.

Mean density of D. antillarum was correlated with wave exposure followed the wave action map on Challet et al. (2012) (r =−0.58, P = 0.0003). However, the wave exposure map did not take into account the existence of small key close to CGB that protect this site from extreme wave action. We adjusted the value of CGB (from 7.9 to 5.5) and obtained a higher correlation (r = −0.80, P = 0.00001).

DISCUSSION

There has been a slow and patchy recovery of D. antillarum populations throughout the Caribbean region (Miller et al., Reference Miller, Adams, Ogden, Ogden and Ebersole2003; Weil et al., Reference Weil, Torres and Ashton2005; Steiner & Williams, Reference Steiner and Williams2006; Debrot & Nagelkerken, Reference Debrot and Nagelkerken2006; Brown-Saracino et al., Reference Brown-Saracino, Peckol, Curran and Robbart2007; Miller et al., Reference Miller, Adams, Ebersole and Ruiz2007; Harborne et al., Reference Harborne, Renaud and Tyler2009; Seller et al., Reference Sellers, Casey, Burge and Koepfler2009; Idjadi et al., Reference Idjadi, Haring and Precht2010; Martín-Blanco et al., Reference Martín-Blanco, González-Sansón, Pina-Amargós and Clero-Alonso2010, Reference Martín-Blanco, Clero-Alonso, González Sansón and Pina-Amargós2011; Ruiz-Ramos et al., Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011; Soto-Santiago and Irizarry-Soto, Reference Soto-Santiago and Irizarry-Soto2013; Levitan et al. Reference Levitan, Edmunds and Levitan2014). Recently, relative high densities have been reported in Jamaica, with up to 4.77 ind m−2 (Idjadi et al., Reference Idjadi, Haring and Precht2010), whereas Florida reported a maximum of 0.33 ind m−2 after more than two decades of the die-off (Chiappone et al., Reference Chiappone, Rutten, Swanson and Miller2008).

The available historical data related to D. antillarum populations is fragmented and scarce in Puerto Rico, limiting our understanding of the recovery dynamics of this important herbivore. Populations of the sea urchin dramatically collapsed from 10 ind m−2 to close to zero after the die-off occurred in 1984 in Puerto Rico (Bak et al., Reference Bak, Carpay and Ruyter van Steveninck1984). Bauer (Reference Bauer1980) reported densities of 13.8 ind m−2 in San Juan a few years before the mass mortality. Current densities remain far from pre-mortality levels (Weil et al. Reference Weil, Torres and Ashton2005). Twenty years after the mass mortality, Ruiz-Ramos et al. (Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011) reported D. antillarum densities of less than 1.0 ind m−2 at 26 localities around Puerto Rico. We surveyed one of their sites and found a small increase in D. antillarum density, from 0.04 ind m−2 (Ruiz-Ramos et al., Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011) to 0.36 0.213 ind m−2. A similar pattern was seen at TM1 where density increased from 0.67 in d m−2 reported in 2004 (Ruiz-Ramos et al., Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011) to 1.09 ±0.411 ind m−2 nine years later. Average densities close to 1 ind m−2 are not far from desired densities between 2 and 3 ind m−2 to control fleshy macroalgae (Steiner & Williams, Reference Steiner and Williams2006). However, the pre-mortality densities for D. antillarum in the north coast of Puerto Rico about 10–14 ind m−2 (Craft, Reference Craft1975; Bak et al., Reference Bak, Carpay and Ruyter van Steveninck1984) was excessive and could have been causing significant reef erosion (Hunter, Reference Hunter1977).

Differences in recovery may be explained by local factors such as predation, overfishing, sedimentation and wave exposure (Harborne et al., Reference Harborne, Renaud and Tyler2009; Seller et al., Reference Sellers, Casey, Burge and Koepfler2009; Ruiz-Ramos et al., Reference Ruiz-Ramos, Hernández-Delgado and Schizas2011; Chollett et al., Reference Chollett, Mumby, Müller-Karger and Hu2012). Wave exposure may model population dynamic at local scale (Chollett et al., Reference Chollett, Mumby, Müller-Karger and Hu2012), and it seems to be controling population abundances of D. antillarum in north-eastern Puerto Rico Archipelago. Diadema sp. is considered weakly resistant to unidirectional hydrodynamic forces, because their morphological features does not allow a large adhesive area to attach to the substrate (Tuya et al., Reference Tuya, Cisneros-Aguirre, Ortega-Borges and Haroun2007). Our analysis based on the wave exposure map of Chollett et al. (Reference Chollett, Mumby, Müller-Karger and Hu2012) supports that leeward sites (TM1, TM2, CGB and MLN) tend to have greater densities of D. antillarum (r = −0.80, P = 0.00001). These four sites are protected from high wave action because of the existence of small keys in front of them. On the contrary, IVD and CGA do not have any physical protection against wind/wave action, and are exposed to more wave action, especially during the winter season. In consonance, lower densities of D. antillarum were reported in these two sites. This result agrees with previous studies where higher densities of Diadema sp. have been linked with wave sheltered areas (Debrot & Nagelkerken, Reference Debrot and Nagelkerken2006; Clemente & Hernández, Reference Clemente and Hernández2008; Hernández et al., Reference Hernández, Clemente, Sangil and Brito2008).

Size structure of D. antillarum changed seasonally with a reduction of mean size observed during the winter season (Figure 3). Previous studies in Puerto Rico have found temporal variability in settlement and recruitment that have an influence in mean sizes (Lugo, Reference Lugo-Ascorbe2004; Williams et al., Reference Williams, García-Sais and Capella2009). The species usually has a maximum spawning peak between October and November (Lessios, Reference Lessios1981), with permanence in the water column of between 35 and 75 d (Leber et al., Reference Leber, Lorenzen, Main and Vaughan2008). A reduction of mean size on February could be a consequence of recruitment. Small individuals are usually more abundant during winter in the Caribbean according to Hunte & Younglao (Reference Hunte and Younglao1988); however, we did not observe a clear seasonal pattern related with the abundance of small individuals among sites (Figure 3). Furthermore, the presence of new recruits (<25 mm) was limited to 1.7% of the 1534 sea urchins measured during the whole study. This suggests a serious deficiency of recruits on Puerto Rico, and points towards recruitment limitation as a possible explanation for the lack of recovery of the species in Puerto Rico. Current low densities of D. antillarum may be having a negative indirect effect on their recovery. A study conducted on the sister species Diadema africanum, former Diadema aff. antillarum (Rodríguez et al., Reference Rodríguez, Hernández, Clemente and Coppard2013), showed a positive relationship between the abundance of adults and the abundance of juveniles (Tuya et al., Reference Tuya, Ortega-Borges, Del Rosario-Pinilla and Haroun2006).

Coral cover and calcareous algae tend to increase, whereas non-calcareous become less abundant when the echinoid is present (Edmunds & Carpenter, Reference Edmunds and Carpenter2001; Carpenter & Edmunds, Reference Carpenter and Edmunds2006). Coral coverage has declined worldwide while non-calcareous seaweeds have increased (Fabricius, Reference Fabricius, Dubinsky and Stambleredts2011). Mean live coral cover reported in this study was less than 25%, but it was even lower at sites with low densities of D. antillarum (Figures 2 and 4). Our results are consistent with others studies conducted in the Caribbean (Williams & Polunin, Reference Williams and Polunin2001; Myhre & Acevedo-Gutiérrez, Reference Myhre and Acevedo-Gutiérrez2007; Martín-Blanco et al., Reference Martín-Blanco, Clero-Alonso, González Sansón and Pina-Amargós2011) where the abundance of the black sea urchin correlates negatively with cover of fleshy macroalgae. However, a recent study found an increment of algal dominance in coral reef systems despite of the presence of D. antillarum (Lacey et al., Reference Lacey, Fourqurean and Collado-Videz2013). This apparent lack of relationship could be a consequence of the current low densities, insufficient to maintain a successful control of seaweeds.

Sea urchin adults not only remove macroalgae while grazing; they create adequate conditions for larval settlement (Butman, Reference Butman1987; Karlson & Levitan, Reference Karlson and Levitan1990; Miller et al., Reference Miller, Adams, Ebersole and Ruiz2007; Rogers & Lorenzen, Reference Rogers and Lorenzen2008). Diadema larvae require a hardground free of non-calcareous algae. The presence of secondary metabolites produced by algae, such as tannins and phenols, can affect the fertilization success, larvae survival, and the detection of the sea urchin adult chemical signals (Pennington, Reference Pennington1985; Butman, Reference Butman1987; Hay & Fenical, Reference Hay and Fenical1988). The observed negative correlation between sea urchin density and fleshy macroalgae (Figure 4), which supports the importance of the species as seaweed controllers (Edmunds & Carpenter, Reference Edmunds and Carpenter2001). The maintenance of relative high abundance of fleshy macroalgae in CGA and IVD (Figure 4) could preclude D. antillarum settlement, and consequently, the recovery of this important herbivore.

The necessity of refuge against predators offers a congruent explanation for the existence of a positive relationship between the abundance of sea urchins and substrate complexity (Tuya et al., Reference Tuya, Martin and Luque2004). It is known that higher environmental complexity provides better refuge for D. antillarum (Valdez & Villalobos, Reference Valdez and Villalobos1978; Weil et al., Reference Weil, Torres and Ashton2005). However, the lack of correlation between mean density and substrate complexity found here, and in other recent study (Martín-Blanco et al., Reference Martín-Blanco, González-Sansón, Pina-Amargós and Clero-Alonso2010) suggests a nonlinear response between rugosity and urchin abundance. At a certain threshold of abundance, habitat features should determine site-to-site variability of population densities, but at low densities there would be available living space for all individual even on less complex hardgrounds.

To conclude, the observed abundance of D. antillarum has not yet returned to pre-mortality levels. However, population densities of the species showed some degree of recovery when we compare with previous studies, enabling at least some degree of control on fleshy macroalgae communities (Figure 4). No significant changes in density occurred between 2011 and 2013, and sites with higher densities were generally located in leeward areas. Further studies are required to understand how this factor may affect the population dynamics of this important herbivore across a wave gradient. Southern coast of the island should be considered given it is sheltered against winter storms. The low relative abundance of small size individuals is worrisome and points towards recruitment limitation as an explanation for the limited recovery of D. antillarum in north Puerto Rico Archipelago three decades after the die-off event.

ACKNOWLEDGEMENTS

We would like to thank the Center for Applied Tropical Ecology and Conservation (CATEC) of the University of Puerto Rico. In addition, the first author wants to express his gratitude to Pedro M. Alcolado, PhD, for his valuable comments and suggestions.

FINANCIAL SUPPORT

The Sea Grant Program PD-295 and R92-110 provided economic and logistic support.

References

REFERENCES

Alvarez-Filip, L., Dulvy, N.K., Gill, J.A., Cote, I.M. and Watkinson, A.R. (2009) Flattening of Caribbean coral reefs: region-wide declines in architectural complexity. Proceedings of the Royal Society, B 276, 30193025.Google Scholar
Bak, R.P.M., Carpay, M.J.E. and Ruyter van Steveninck, E.D. (1984) Densities of the sea urchin Diadema antillarum before and after mass mortalities on the coral reefs of Curacao. Marine Ecology 17, 105108.CrossRefGoogle Scholar
Bauer, J.C. (1980) Observations of geographical variations in population density of the echinoid Diadema antillarum within the Western North Atlantic. Bulletin of Marine Science 30, 509515.Google Scholar
Bellwood, D.R., Hughes, T.P., Folke, C. and Nystrom, M. (2004) Confronting the coral reef crisis. Nature 429, 827833.Google Scholar
Bythell, J.C., Gladfelter, E.H. and Bythell, M. (1993) Chronic and catastrophic natural mortality of three common Caribbean reef corals. Coral Reefs 12, 143152.CrossRefGoogle Scholar
Bythell, J.C., Hillis-Starr, Z.M. and Rogers, C.S. (2000) Local variability but landscape stability in coral reef communities following repeated hurricane impacts. Marine Ecology Progress Series 204, 93100.Google Scholar
Box, G.E.P. and Cox, D.R. (1964) An analysis of transformations (with discussion). Journal of the Royal Statistical Society, B 26, 211252.Google Scholar
Brown-Saracino, J., Peckol, P., Curran, H.A. and Robbart, M.L. (2007) Spatial variation in sea urchins, fish predators, and bioerosion rates on coral reefs of Belize. Coral Reefs 26, 7178.Google Scholar
Butman, C.A. (1987) Larval settlement of soft sediment in vertebrates: the spatial scales of pattern explained by active habitat selection and the emerging role of hydrodynamical processes. Oceanography and Marine Biology: an Annual Review 25, 113165.Google Scholar
Carpenter, R.C. and Edmunds, P.J. (2006) Local and regional scale recovery of Diadema promotes recruitment of scleractinian corals. Ecology Letters 9, 271280.Google Scholar
Chiappone, M., Rutten, L.M., Swanson, D.W. and Miller, S.L. (2008) Population status of the urchin Diadema antillarum in the Florida Keys 25 years after the Caribbean mass mortality. In Proceedings of the 11th International Coral Reef Symposium, Fort Lauderdale, Florida, 7–11 July, pp. 706–710.Google Scholar
Chollett, I., Mumby, P.J., Müller-Karger, F. E. and Hu, C. (2012) Physical environments of the Caribbean Sea. Limnology and Oceanography 57, 12331244.CrossRefGoogle Scholar
Clemente, S. and Hernández, J.C. (2008) Influence of wave exposure and habitat complexity in determining spatial variation of the sea urchin Diadema aff. antillarum Echinoidea: Diadematidae) populations and macroalgal cover (Canary Islands–Eastern Atlantic Ocean). Biología Tropical 56, 229254.Google Scholar
Craft, L.L. (1975) Aspects of the biology of the crab Percnon gibbesi (Milne Edwards) and its commensal association with the sea urchin Diadema antillarum Phillipi. MSc thesis. University of Puerto Rico.Google Scholar
De'ath, G., Fabricius, K.E., Sweatman, H. and Puotinen, M. (2012) The 27-year decline of coral cover on the Great Barrier Reef and its causes. Proceedings of the National Academy of Sciences of the United States of America 109, 1799517999.CrossRefGoogle ScholarPubMed
Debrot, A.O. and Nagelkerken, I. (2006) Recovery of the long-spined sea urchin Diadema antillarum in Curacao (Netherlands Antilles) linked to lagoonal and wave sheltered shallow rocky habitats. Bulletin of Marine Science 79, 415424.Google Scholar
Edmunds, P.J. and Carpenter, R.C. (2001) Recovery of Diadema antillarum reduces macroalgal cover and increases abundance of juvenile corals on a Caribbean Reef. Proceedings of the National Academy of Sciences of the United States of America 98, 50675071.Google Scholar
Fabricius, K.F. (2011) Factors determining the resilience of coral reefs to eutrophication: a review and conceptual model. In Dubinsky, Z. and Stambleredts, N. (eds) Ecosystems in transition. Coral reefs. Dordrecht, The Netherlands: Springer, pp. 493505.CrossRefGoogle Scholar
Gardner, T.A., Côté, I.M., Gill, J.A., Grant, A. and Watkinson, A.R. (2003) Long-term region-wide declines in Caribbean corals. Science 301, 958960.Google Scholar
Harborne, A.R., Renaud, P.G. and Tyler, E.H. (2009) Reduced density of the herbivorous urchin Diadema antillarum inside a Caribbean marine reserve linked to increased predation pressure by fishes. Coral Reefs 28, 783791.Google Scholar
Hay, M.E. and Fenical, W. (1988) Marine plant–herbivore interactions: the ecology of chemical defense. Annual Review of Ecology and Systematics 19, 111145.Google Scholar
Hernández, J.C., Clemente, S., Sangil, C. and Brito, A. (2008) The key role of the sea urchin Diadema aff. antillarum in controlling macroalgae assemblages throughout the Canary Islands (eastern subtropical Atlantic): a spatio-temporal approach. Marine Environmental Research 66, 259270.Google Scholar
Hughes, T.P., Keller, B.D., Jackson, J.B. and Boyle, M.J. (1985) Mass mortality of the Echinoid Diadema antillarum Philippi in Jamaica. Bulletin of Marine Science 36, 377–84.Google Scholar
Hughes, T.P. (1994) Catastrophes, phase shifts, and large-scale degradation of a Caribbean coral reefs. Science 9, 15471551.Google Scholar
Hughes, T.P., Graham, N.A.J., Jackson, J.B.C., Mumby, P.J. and Steneck, R.S. (2010) Rising to the challenge of sustaining coral reef resilience. Trends in Ecology and Evolution 25, 633642.Google Scholar
Hunte, W., Côté, I. and Tomascik, T. (1986) On the dynamics of the mass mortality of Diadema antillarum in Barbados. Coral Reefs 4, 135139.Google Scholar
Hunte, W. and Younglao, D. (1988) Recruitment and population recovery of Diadema antillarum (Echinodermata; Echinoidea) in Barbados. Marine Ecology Progress Series 45, 109119.Google Scholar
Hunter, I.G. (1977) Sediment production by Diadema antillarum on a Barbados fringing reef. In Proceedings of the Third International Coral Reef Symposium, Miami, Florida: Rosenstiel School, University of Miami, May, Volume 2, pp. 105–110.Google Scholar
Idjadi, J.A., Haring, R.N. and Precht, W.F. (2010) Recovery of the sea urchin Diadema antillarum promotes scleractinian coral growth and survivorship on shallow Jamaican reefs. Marine Ecology Progress Series 403, 91100.Google Scholar
Karlson, R.H. and Levitan, D.R. (1990) Recruitment-limitation in open populations of Diadema antillarum: an evaluation. Oecologia 82, 4044.CrossRefGoogle ScholarPubMed
Knowlton, N. (2001) Sea urchin recovery from mass mortality: new hope for Caribbean coral reefs? Proceedings of the National Academy of Sciences of the United States of America 98, 48224824.Google Scholar
Kohler, K.E., and Gill, S.M. (2006) Coral Point Count with Excel extensions (CPCe): a Visual Basic program for the determination of coral and substrate coverage using random point count methodology. Computers and Geosciences 32, 12591269.Google Scholar
Kramer, P.A., Kramer, P.R. and Ginsburg, R.N. (2003) Assessment of the Andros island reef system, Bahamas (Part 1: Stony corals and algae). Atoll Research Bulletin 496, 77100.Google Scholar
Lacey, E., Fourqurean, J.W. and Collado-Videz, L. (2013) Increased algal dominance despite presence of Diadema antillarum populations on a Caribbean coral reef. Bulletin of Marine Science 89, 603620.Google Scholar
Leber, K., Lorenzen, K., Main, K.M.M. and Vaughan, D. (2008) Developing restoration methods to aid in recovery of a key herbivore Diadema antillarum on Florida coral reefs. Mote Technical Report 1295, 26 pp.Google Scholar
Lessios, H.A. (1981) Reproductive periodicity of the echinoids Diadema and Echinometra on the two coasts of Panama. Journal of Experimental Marine Biology and Ecology 50, 4761.Google Scholar
Lessios, H.A, Cubit, J.D., Robertson, D.R., Shulman, M.J., Parker, M.R., Garrity, S.D. and Levings, S.C. (1984) Mass mortality of Diadema antillarum on the Caribbean coast of Panama. Coral Reefs 3, 173182.Google Scholar
Lessios, H.A. (2005) Diadema antillarum populations in Panama twenty years following mass mortality. Coral Reefs 24, 125127.Google Scholar
Levitan, D., Edmunds, P.J. and Levitan, K.E. (2014) What makes a species common? No evidence of density-dependent recruitment or mortality of the sea urchin Diadema antillarum after the 1983–1984 mass mortality. Oecologia 175, 117128.CrossRefGoogle ScholarPubMed
Liddell, W.D. and Ohlhorst, S.L. (1986) Changes in benthic community composition following the mass mortality of Diadema at Jamaica. Journal of Experimental Marine Biology and Ecology 95, 271278.Google Scholar
Littler, M.M., Littler, D.S. and Titlyanov, E.A. (1991) Comparisons of N-limited and P-limited productivity between high granitic islands versus low carbonate atolls in the Seychelles archipelago. A test of the relative-dominance paradigm. Coral Reefs 10, 199209.CrossRefGoogle Scholar
Lugo-Ascorbe, M.A. (2004) Population status of the black sea urchin Diadema antillarum (Philippi) in La Parguera, Puerto Rico, 20 years after the mass mortality event. MSc thesis. Department of Marine Sciences, University of Puerto Rico.Google Scholar
Martín-Blanco, F., González-Sansón, G., Pina-Amargós, F. and Clero-Alonso, L. (2010) Abundance, distribution and size structure of Diadema antillarum (Echinodermata: Diadematidae) in South Eastern Cuban coral reefs. Biología Tropical 58, 663676.Google Scholar
Martín-Blanco, F., Clero-Alonso, F., González Sansón, G. and Pina-Amargós, F. (2011). Influence of Diadema antillarum populations (Echinodermata: Diadematidae) on algal community structure in Jardines de la Reina, Cuba. Biología Tropical 59, 11491163.Google Scholar
Miller, R.J., Adams, A.J., Ogden, N.B., Ogden, J.C. and Ebersole, J.P. (2003) Diadema antillarum 17 years after mass mortality: is recovery beginning on St Croix? Coral Reefs 22, 181187.Google Scholar
Miller, R.J., Adams, A.J., Ebersole, J.P. and Ruiz, E. (2007) Evidence for positive density-dependent effects in recovering Diadema antillarum populations. Journal of Experimental Marine Biology and Ecology 349, 215222.Google Scholar
Myhre, S. and Acevedo-Gutiérrez, A. (2007) Recovery of sea urchin Diadema antillarum populations is correlated to increased coral and reduced macroalgal cover. Marine Ecology and Progress Series 329, 205210.CrossRefGoogle Scholar
NOAA Coral Reef Watch (2000) NOAA Coral Reef Watch Operational 50-km Satellite Coral Bleaching Degree Heating Weeks Product, Jan. 1 2001–Dec. 31 2011. Silver Spring, MD: NOAA Coral Reef Watch. Available at: http://coralreefwatch.noaa.gov/satellite/hdf/index.php (accessed 22 April 2014).Google Scholar
Noriega, N., Pauls, S. and Mónaco, C. (2006) Abundancia de Diadema antillarum (Echinodermata: Echinoidea) en las costas de Venezuela. Biología Tropical 54, 793802.Google Scholar
Pennington, J.T. (1985) The ecology of fertilization of echinoid eggs: the consequence of sperm dilution, adult aggregation, and synchronous spawning. Biological Bulletin. Marine Biological Laboratory, Woods Hole 169, 417430.Google Scholar
Pinheiro, J., Bates, D. and DebRoy, S. (2013) nlme: linear and nonlinear mixed effects models. R package version 3.1-111. Vienna: R Foundation for Statistical Computing.Google Scholar
Phinney, J.T., Muller-Karger, F., Dustan, P. and Sobel, J. (2001) Using remote sensing to reassess the mass mortality of Diadema antillarum 1983–1984. Conservation Biology 15, 885–881.Google Scholar
R Core Team (2013) R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.Google Scholar
Rivera, J.A. and Vicente, V.P. (1976) Diadema antillarum (Philippi): enhancement of coral species in Thalassia beds. Proceedings of the Association of Island Marine Laboratories of the Caribbean 12, 18.Google Scholar
Rodríguez, A., Hernández, J.C., Clemente, S. and Coppard, S.E. (2013) A new species of Diadema (Echinodermata: Echinoidea: Diadematidae) from the eastern Atlantic Ocean and a neotype designation of Diadema antillarum (Philippi, 1845). Zootaxa 3636, 144170.Google Scholar
Rogers, A. and Lorenzen, K. (2008) Recovery of Diadema antillarum and the potential for active rebuilding measures: modelling population dynamics. In Proceedings of the 11th International Coral Reef Symposium. Fort Lauderdale, FL, 7–11 July, pp. 956–960.Google Scholar
Ruiz-Ramos, D.V., Hernández-Delgado, E.A. and Schizas, N.V. (2011) Population status of the long-spined urchin Diadema antillarum in Puerto Rico 20 years after a mass mortality event. Bulletin of Marine Science 87, 113127.Google Scholar
Sammarco, P.W., Levinton, J.S. and Ogden, J.C. (1974) Grazing and control of coral reef community structure by Diadea antillarum. Philippi (Echinodermata: Echinoidea): a preliminary study. Journal of Marine Research 32, 4753.Google Scholar
Sammarco, P.W. (1982) Effects of grazing by Diadema antillarum Philippi (Echinodermata: Echinoidea) on algal diversity and community structure. Journal of Experimental Marine Biology and Ecology 65, 83105.Google Scholar
Sellers, J.A., Casey, L.O., Burge, E.J. and Koepfler, E.T. (2009) Population growth and distribution of Diadema antillarum at Discovery Bay, Jamaica. The Open Marine Biology Journal 3, 105111.Google Scholar
Soto-Santiago, F.J. and Irizarry-Soto, E. (2013) The sea urchin Diadema antillarum (Echinodermata, Equinoidea), algal cover and juvenile coral densities in La Parguera, Puerto Rico. Cuadernos de Investigación, UNED 5, 3339.Google Scholar
Steiner, S.C. and Williams, S.M. (2006) A recent increase in the abundance of the echinoid Diadema antillarum in Dominica (Lesser Antilles): 2001–2005. Biología Tropical 54, 97103.Google Scholar
Steneck, R.S. (1988) Herbivory on coral reefs: a synthesis. In Proceedings of the sixth International Coral Reef Symposium, Australia, Volume 1, 8–12 August, pp. 37–49.Google Scholar
Tuya, F., Martin, J.A., and Luque, A. (2004) Patterns of nocturnal movement of the long-spined sea urchin Diadema antillarum (Philippi) in Gran Canaria. Helgoland Marine Research 58, 2631.Google Scholar
Tuya, F., Ortega-Borges, L., Del Rosario-Pinilla, A.B. and Haroun, R.J. (2006) Spatio-temporal variability in a key herbivore, the long-spined black sea urchin (Diadema antillarum, Echinodermata: Echinoidea) in the Canary Islands. Journal of the Marine Biological Association of the United Kingdom 86, 791797.Google Scholar
Tuya, F., Cisneros-Aguirre, J., Ortega-Borges, L. and Haroun, R.J. (2007) Bathymetric segregation of sea urchins on reefs of the Canarian Archipelago: role of flow-induced forces. Estuarine, Coastal and Shelf Science 73, 481488.Google Scholar
USGS (United States Geological Survey) (1996) Atlas of ground-water resources in Puerto Rico and the U.S. Virgin Islands. In Veve, T.D. and Taggart, B.E. (eds) Water resources investigation report, 94–4198. Reston, VA: USGS, 151 pp.Google Scholar
Valdez, M.F. and Villalobos, C.R. (1978) Distribución espacial, correlación con el substrato y grado de agregación de Diadema antillarum Philippi (Echinodermata: Echinoidea). Biología Tropical 26, 237245.Google Scholar
Venables, W.N. and Ripley, B.D. (2002) Modern applied statistics with S. 4th editionNew York: Springer, 495 pp.Google Scholar
Vicente, V.P. (1987) The ecology of the encrusting demosponge Chondrilla nucula (Schmidt) in a coral reef community in Puerto Rico. PhD thesis. Department of Marine Sciences, University of Puerto Rico.Google Scholar
Vicente, V.P. and Goenaga, C. (1984) Mortandades masivas del erizo de mar Diadema antillarum (Philippi) en Puerto Rico. Centro para Estudios Energéticos y Ambientales, Div. Ecología Marina, University of Puerto Rico. Reporte CEER-M-195, 30 pp.Google Scholar
Weil, E., Torres, J.L. and Ashton, M. (2005) Population characteristics of the sea urchin Diadema antillarum in La Parguera, Puerto Rico, 17 years after the mass mortality event. Biología Tropical 53, 219231.Google Scholar
Williams, I.D. and Polunin, N.V.C. (2001) Large-scale associations between macroalgal cover and grazer biomass on mid-depth reefs in the Caribbean. Coral Reefs 19, 358366.CrossRefGoogle Scholar
Williams, S.M., García-Sais, J.R. and Capella, J. (2009) Temporal variation of early larval stages of the long-spined sea urchin Diadema antillarum in La Parguera, Puerto Rico. Caribbean Journal of Science 45, 110117.Google Scholar
Williams, S.M., Benavides-Serrato, M., García-Arrarás, J.E., Hernández-Delgado, E.A. and Rodríguez-Barreras, R. (2013) Review of echinoderm research in Puerto Rico, with the focus on biological and ecological aspects. In Alvarado-Barrientos, J.J. and Solis-Marin, F.A. (eds) Echinoderm research and diversity in Latin America. New York: Springer, pp. 437469.Google Scholar
Zar, J.H. (2010) Biostatistical analysis. 5th edition. Upper Saddle River, NJ: Prentice-Hall, 944 pp.Google Scholar
Figure 0

Fig. 1. Study area and sampling sites in Puerto Rico. Cerro Gordo A (CGA), Cerro Gordo B (CGB), Isla Verde (IVD), Luquillo (LQY), Tamarindo 1 (TM1), Tamarindo 2 (TM2) and Melones (MLN).

Figure 1

Table 1. Two-way repeated measure analysis of variance on the abundance (a) and the two-way analysis of variance test on the size of the sea urchin Diadema antillarum (b). *this value comes from a χ2 test that compares models with and without site as a random effect.

Figure 2

Fig. 2. Spatial and temporal density of Diadema antillarum in Puerto Rico from August 2011 to August 2013 at CGA, CGB, IVD, LQY, TM1, TM2 and MLN (see Methods for acronyms). The X axis represents time with A-11 (August 2011), F-12 (February 2012), etc. The Y axis is the mean density of D. antillarum (ind m−2). Bars represent 95% confidence interval of the mean and red line represents a polynomial fit. Last graph shows the overall temporal change of mean density using pooled data.

Figure 3

Fig. 3. Seasonal size-classes distribution of Diadema antillarum at CGA, CGB, IVD, LQY, TM1, TM2 and MLN in Puerto Rico (see Materials and Methods for acronyms). Black bar: small size (<40 mm), white: medium size (41–60 mm) and dark grey: large size (>60 mm). Last graph shows the temporal variation of mean size of D. antillarum of all sites.

Figure 4

Fig. 4. Spatial variation of the percentage cover of fleshy macroalgae (NCAL), turf, calcareous algae (CALG) and live coral (LCOR) at CGA, CGB, IVD, LQY, TM1, TM2 and MLN in Puerto Rico (see Materials and Methods for acronyms). Last graph shows the relationship between density of Diadema antillarum and percentage cover of fleshy macroalgae (NCAL). Box represents 95% confidence interval of the mean.