Hostname: page-component-745bb68f8f-s22k5 Total loading time: 0 Render date: 2025-02-06T14:51:45.004Z Has data issue: false hasContentIssue false

Revisiting Connell: competition but not as we know it

Published online by Cambridge University Press:  25 May 2017

Jenna M. Gordon
Affiliation:
Marine Biology and Ecology Research Centre, School of Biological and Marine Sciences, Plymouth University, Davy Building 622, Drake Circus, Plymouth PL4 8AA, UK
Antony M. Knights*
Affiliation:
Marine Biology and Ecology Research Centre, School of Biological and Marine Sciences, Plymouth University, Davy Building 622, Drake Circus, Plymouth PL4 8AA, UK
*
Correspondence should be addressed to: A. M. Knights, Marine Biology and Ecology Research Centre, School of Biological and Marine Sciences, Plymouth University, Davy Building 622, Drake Circus, Plymouth PL4 8AA, UK email: antony.knights@plymouth.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Space is one of the primary limiting resources for organisms on the intertidal rocky shore. This paper examined the effect of reduced density on key traits (mortality and growth) on the intertidal barnacles, Chthamalus montagui and Semibalanus balanoides, on the mid-shore in Plymouth, UK. Intra- and interspecific treatments comprising of C. montagui and S. balanoides were manipulated to reduce densities at two similar sites. Changes in mortality and operculum growth were assessed over an 8-week period using digital photography. Covariates of growth included nearest neighbour distance, competition between closest pairs and initial size. Conflicting patterns were observed when comparing growth rates between treatments and sites. At Site 1, interspecific treatments had a lower growth rate than intraspecific treatments, whereas at Site 2, interspecific growth rates were higher. ANCOVA showed that nearest neighbour distance had no significant effect on growth, but when comparing differences in growth of closest neighbouring pairs, C. montagui treatment showed evidence of competition whereas S. balanoides did not. ANCOVA analysis indicated no difference in growth between each outcome of pair competition, suggesting winners are initially bigger than losers. Comparisons of mortality between treatments indicated mortality over time with no significant differences observed between treatments, but response surface methodology (RSM) revealed no effects of competition on mortality of S. balanoides, but negative effects of both intra- and interspecific competition on C. montagui survivorship. Examination of natural populations of barnacles in the mid-shore indicated there was strong spatial variation in growth rates, perhaps driven by small-scale differences within sites.

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

INTRODUCTION

The effect of competition on natural populations and community dynamics has long been of interest to ecologists (Connell, Reference Connell1961a, Reference Connellb; Tilman, Reference Tilman1982; Strong et al., Reference Strong, Simberloff, Abele and Thistle1984), especially the role of intraspecific and interspecific resource competition as determining factors of the structure of both terrestrial and marine populations (Hart & Marshall, Reference Hart and Marshall2009; Caro et al., Reference Caro, Guiñez, Ortiz and Castilla2011; Shinen & Navarrete, Reference Shinen and Navarrete2014).

Competition acts as an important feedback loop that controls population density and growth rate of the population and individuals (Begon et al., Reference Begon, Townsend and Harper2006). Intraspecific competition between individuals of the same species commonly leads to mortality when the resources needed to sustain them, such as food and space, become limiting (Moore, Reference Moore1935; Hixon et al., Reference Hixon, Pacala and Sandin2002; Begon et al., Reference Begon, Townsend and Harper2006; Knights & Walters, Reference Knights and Walters2010). This mechanism can be described by a logistic growth model which describes the negative effect of population size on growth rate until the carrying capacity is reached (birth rate equals death rate) and resources are no longer limited (Hixon et al., Reference Hixon, Pacala and Sandin2002; Neal, Reference Neal2004). While a change in population size is the emergent result, intraspecific competition may in fact alter survivorship at the scale of the individual rather than at the population level, for example, as neighbouring individuals compete for a resource (Begon et al., Reference Begon, Townsend and Harper2006) reducing their growth and size (Hixon et al., Reference Hixon, Pacala and Sandin2002).

In addition to intraspecific competition, population size has also been shown to be dependent on interspecific interactions. The Lotka–Volterra model (Lotka, Reference Lotka1925; Volterra, Reference Volterra1926) states that when two species occur together, the growth rates of both species are affected by the presence of each other (sensu ‘interference competition’; Neal, Reference Neal2004) as a result of direct competition for the same resource which affects their growth or survival (Feldhamer et al., Reference Feldhamer, Dricamer, Vessey, Merritt and Krajewski2007; Reece et al., Reference Reece, Urry, Cain, Wasserman, Minorsky and Jackson2011). This can affect the distribution and abundance of different species in natural communities. Gause (Reference Gause1934) used laboratory experiments to observe resource competition between two closely related species, namely Paramecium caudatum and Stylonychia mytilus. When grown separately, each population grew rapidly before reaching asymptote at their carrying capacity (k), yet when both species co-occurred, growth rates were reduced and the carrying capacity was lower. Moreover, S. mytilus appeared to partially outcompete P. caudatum, evident as greater reductions in the growth of P. caudatum than S. mytilus.

A comparison of the relative strengths of intra- and interspecific competition provides an indication of how species may coexist. The Lotka–Volterra model (Lotka, Reference Lotka1925; Volterra, Reference Volterra1926) suggests that species co-exist when intraspecific competition is stronger than interspecific competition (Connell, Reference Connell1983; Ying et al., Reference Ying, Liao, Wang, Hui and Liu2014); the competitively superior species is ‘self-limited’ by competition between individuals below a density threshold that is necessary to eliminate the other species (Connell, Reference Connell1983). If the competitively superior species is not ‘self-limited’, then the weaker competitor may be eliminated unless external sources of mortality, such as predation or disturbance, limit the population of the superior species (Shinen & Navarrete, Reference Shinen and Navarrete2014). While a well-established principle, recent studies have suggested that species are not ‘co-existing’ per se, but instead a species is slowly driving another to extinction (Siepielski & McPeek, Reference Siepielski and McPeek2010; Shinen & Navarrete, Reference Shinen and Navarrete2014). Nonetheless, the relative strength of competition is important to determine if we are to predict the possibility of co-existence between species and understand changes in the density of competing species (Connell, Reference Connell1983).

Barnacles have long been used to test hypotheses of intra- and interspecific competition (e.g. Connell, Reference Connell1961a, Reference Connellb; Wethey, Reference Wethey1983; Jenkins et al., Reference Jenkins, Murua and Burrows2008) and have been shown to demonstrate both intra- and interspecific competition by crushing or overgrowing neighbours (Connell, Reference Connell1961a, Reference Connellb; Wethey, Reference Wethey1983; Jenkins et al., Reference Jenkins, Murua and Burrows2008). In intertidal systems, barnacles are excluded from higher regions of the shore by physical stress (e.g. desiccation) and reduced in number on lower shore heights by biological control (e.g. predation and competition). Their small size, dense concentrations and intertidal location make barnacles an ideal model organism for manipulation in field experiments (Leslie, Reference Leslie2005; López et al., Reference López, Catalán, Barriga and López2014). For instance, the survival of individuals can be determined accurately by simply mapping the position of all the members of a population and then following the same individuals by regular censuses (Connell, Reference Connell1961b).

Space is one of the primary limiting resources for barnacles on the intertidal rocky shore (Connell, Reference Connell1961a, Reference Connellb; Leslie, Reference Leslie2005). At high densities, intraspecific competition for space may negatively affect survival, cause changes in growth rates, and reduce reproductive activity (Barnes & Powell, Reference Barnes and Powell1950; López et al., Reference López, Catalán, Barriga and López2014) and success. In contrast, high densities have been reported to facilitate survival by buffering individuals from interspecific competitive pressures, consumers, physical disturbance and physiological distress (Bertness, Reference Bertness1989; Leslie, Reference Leslie2005). It is suggested that barnacles in dense aggregations grow more slowly than adjacent isolated individuals as food in the water flowing over the surface is shared among more individuals (Moore, Reference Moore1935).

In a seminal study by Connell (Reference Connell1961a), it was found that barnacles within the mid-intertidal zone undergo significant interspecific competition. Undertaken on the Isle of Cumbrae in the Firth of Clyde, Connell demonstrated interspecific competition between two co-existing species, Chthamalus stellatus (now recognized as Chthamalus montagui and referred to as C. montagui herein, see Southward, Reference Southward1976) and Semibalanus balanoides. Chthamalus montagui generally occurred above S. balanoides and was shown to be able to settle lower on the shore, but was unable to survive as a result of being eliminated by S. balanoides over a 1-year period. Connell argued that the short supply of a common resource caused the exclusion of C. montagui as space for attachment and growth was limited, and the poor survival of C. montagui in the lower shore was as a result of being outcompeted by the faster growing species, S. balanoides (Connell, Reference Connell1961a).

This paper revisits Connell's study (Reference Connell1961a), re-testing his assumptions of intra- and interspecific competition by way of manipulated densities and combinations of S. balanoides and C. montagui. The study aims to determine whether (1) growth rates vary between intra- and interspecific treatments, and (2) if survival rate varies between intra- and interspecific treatments over time. Small-scale effects on growth between closely interacting (neighbouring) individuals are also tested.

MATERIALS AND METHODS

The study was carried out at Mount Batten, Plymouth, UK (Figure 1) between September and December 2014. Mount Batten is a headland of limestone protruding into Plymouth Sound (see Knights et al., Reference Knights, Firth, Thompson, Yunnie, Hiscock and Hawkins2016 for description of the area). Two sites on the shore were identified and defined as limestone rock surfaces of similar aspect, gradient, tidal exposure and orientation (~ south-west facing; see Figure 1 and www.EMODnet.eu for more information) where Semibalanus balanoides and Chthamalus montagui are locally abundant and coexist in the mid-shore. Locations were at the same tidal height and separated by >50 m. The study area was intentionally limited in order to reduce variability caused by differences in tidal exposure, salinity, temperature and light (Connell, Reference Connell1961b).

Fig. 1. Geographic location of the study in Plymouth, UK. Inset: the intertidal rocky shore at Mount Batten (50°21′N 4°07′W).

To test for evidence and strength of intraspecific and interspecific competition in CM and SB, three treatments were established in areas characterized by 100% cover of adult barnacles. Intraspecific treatment patches contained either SB or CM, and interspecific (mixed) treatment patches contained both species. Patches were 5 × 5 cm and located randomly within the mid-shore (determined using tide-tables) where both species are roughly equally abundant. The location of each patch was recorded using GPS (Garmin eTrex10, USA) and barnacles were removed from each patch using forceps to manipulate the density and occurrence of species within each patch (Figure 2B, D) to allow a response surface experimental design to be used (see Inouye, Reference Inouye2001) to test hypotheses about competition. The area around each 5cm2 patch was cleared of all barnacles using a paint scraper to reduce the likelihood of edge effects (Volkenborn et al., Reference Volkenborn, Hedtkamp, van Beusekom and Reise2007). Each quadrat (N = 30 per site) was photographed (Panasonic DMC-FS16) prior to the removal of any barnacles, following manipulation at Time 0 and subsequently at 2, 4, 6 and 8 weeks post-manipulation. Photographs were used to calculate density (and therefore mortality) and to estimate individual growth over time.

Fig. 2. Schematics and photograph of 5 cm2 manipulative treatment plots used to test for evidence of intraspecific (A, B – Chthamalus montagui only) and interspecific (C) competition and (D) the response surface experimental design with example species densities in each patch at Time-0 as recommended by Inouye (Reference Inouye2001).

Estimates of growth rate by species and treatment

Growth was measured as the change in the length of the operculum over time (Wethey, Reference Wethey1983; Jenkins et al., Reference Jenkins, Murua and Burrows2008, Burrows et al., Reference Burrows, Jenkins, Robb and Harvey2010), rather than the total length of the barnacle (rostro-carinal) as this metric can be severely affected by crowding, as well as other micro-topographical features of the rock surface. Operculum length of individual barnacles was measured using a photograph (e.g. Figure 2B) that had been scaled in the image analysis programme, ImageJ (Schneider et al., Reference Schneider, Rasband and Eliceiri2012), and individuals geo-referenced allowing their growth to be tracked over time. Average growth rate (per species and per 2-week period) was calculated for each patch and treatment type.

Mortality of barnacles

Barnacle mortality was measured using images from Time-0 and after 8 weeks. Images of patches were overlaid with one-another so that individual barnacles could firstly be identified, speciated (in mixed patches) and a binary code applied to whether they were dead (0) or alive (1) after 8 weeks. Individuals were identified as dead by the absence of opercular plates. These data were then used to calculate proportional mortality of each species in each patch.

Growth and distance from a neighbour

Photographs were also used to calculate the nearest neighbour distances (NND) of all barnacles in a patch. Images were imported in ArcGIS (ArcInfo 10.2.2), georeferenced and a point applied to the centre of each barnacle in an image. The ordinal distance between the centre of all barnacles were calculated using the spatial analysis toolkit (nearest neighbour tool) and NNDs used as a covariate to test the hypothesis that shorter NNDs would lead to a reduction in operculum growth as a result of competition.

Evidence of ‘winners’ and ‘losers’ in intraspecific patches

To determine if small-scale intraspecific competition occurs between barnacles in close proximity to each other, barnacles closest to each other (based on NND) were paired to test if either individual was outcompeting its ‘pair’ (measured as a difference in growth rate). Individuals were classified as the ‘winner’ and ‘loser’ based on growth rate differences; the individual that grew more was the ‘winner’. When no evidence of growth rate differences was seen (i.e. equal growth in paired individuals), the contest was considered a ‘draw’. The operculum length of each barnacle at the start of the experiment was also used as a covariate to account for potential differences in growth rate based on starting size (Moore, Reference Moore1939).

Statistical analysis

Growth rates were compared using an orthogonal two-factor ANOVA with the factors: (1) Site (2 levels, random), and (2) Treatment (3 levels: CM only, SB only, CM + SB (mixed)). The outcome of paired competitions as determined by growth rate was also tested using two-factor ANOVA with the factors: (1) Outcome (3 levels: Win, Lose or Draw), and (2) Treatment (2 levels: CM only, SB only). Significant differences between means were compared using post-hoc pairwise comparisons (Tukey HSD, ‘car’ package).

The effect of NND on operculum growth was testing using Analysis of Covariance (ANCOVA) with the fixed factors, Treatment (3 levels: CM only, SB only, CM + SB (mixed)), Outcome (see above) and continuous covariate (NND in mm). Step-wise model simplification of the maximal ANCOVA model was used to test between slope and intercept parameters for the relationship between growth and initial size (covariate) and categorical factor levels. Akaike Information Criteria (AIC) was used determine the best-fitting model (Burnham & Anderson, Reference Burnham and Anderson2004).

Change in barnacle density and mortality over 8 weeks between treatments was compared using linear mixed effects models (lmer) with an auto-regression (1) correlation structure to account for possible temporal autocorrelation effects. Local regression (loess) was used to describe change in average density over time between treatments. For all analyses above, data were tested for residual normality and homoscedasticity prior to statistical analyses, and in cases of significance (P < 0.05), data were log transformed. Multicollinearity was examined using the variance inflation factor (VIF), values of which were <10, and therefore unlikely to affect regression outcomes (O'Brien, Reference O'Brien2007).

To test for evidence of intraspecific and interspecific competition, mortality within plots was compared using a response surface method (after Box et al., Reference Box, Hunter and Hunter2005; see Inouye, Reference Inouye2001, Lenth, Reference Lenth2009 and Figure 7D for mixture details) and ANOVA tests. The effect of CM and SB density on mortality of CM and SB was modelled using first-order polynomial regression models and ANOVA tests, R 2 for goodness-of-fit and plotted using perspective plots.

All analyses were performed using the software R (R Core Team, 2016) using the R packages, ‘graphics’, ‘lme4’ (Bates et al., Reference Bates, Maechler, Bolker and Walker2015) and ‘rsm’ (Lenth, Reference Lenth2009).

RESULTS

Variation of operculum growth in relation to treatment (Tr) and site (Si)

There was marked spatial variation in operculum growth rates among treatments between sites (Table 1, Si × Tr, P < 0.0001). In intraspecific treatments, CM grew ~16% more at Site 2 than Site 1, whereas for SB, the pattern was reversed with SB growing ~21% more at Site 1 than 2 (Figure 3). In the interspecific treatment containing both barnacle species, average growth was significantly higher (~44%) at Site 2 than Site 1. A comparison of the three treatments across sites indicated little variation in operculum growth at Site 1 irrespective of treatment, whereas at Site 2, growth was highest in the CM + SB and CM only treatments, and lowest in the SB only treatment (Figure 3).

Fig. 3. Mean (±SE; N = 5) log growth of barnacle opercula in three treatments: CM + SB (mixed species plots), CM (Chthamalus montagui only) and SB (Semibalanus balanoides only) at two sites after 8 weeks. Post-hoc pairwise comparison results are shown as letters above columns, where different letters indicates significant differences between ‘treatment × site’ combinations (P < 0.01).

Table 1. ANOVA of operculum growth (mm) in relation to site and treatment at Mount Batten, Plymouth, in 2014.

Variation in operculum growth in relation to nearest neighbour distance

When data were combined within a single analysis without consideration of barnacle ‘pairings’, there was no effect of nearest neighbour distance on the operculum growth of an individual in any of the treatments (Table 2, NND, P > 0.01). Operculum growth was highly variable, especially when NNDs were short (<5 mm) (Figure 4) and individuals located further away from a neighbour (i.e. large NND) showed no significant increase in growth in comparison to those with barnacles closer to them.

Fig. 4. Scatterplot of operculum growth in individual barnacles after 8 weeks in relation to distance (mm) from its nearest neighbour. Different treatment combinations are shown: Semibalanus balanoides + Chthamalus montagui (white); Chthamalus montagui only (grey); and Semibalanus balanoides only (black).

Table 2. Analysis of covariance (ANCOVA) of barnacle operculum growth in relation to nearest neighbour distance and treatment (SB only; CM only; CM + SB).

However, when barnacles were ‘paired’ to those closest to each other and classified as either a ‘winner’, ‘loser’ or ‘draw’ based on their growth after 8 weeks, ANOVA indicated significant differences in growth between treatment (Table 3, Outcome × Treatment, P < 0.05). In general, losers grew 60% less than individuals classified as winners. Post-hoc comparisons revealed clear differences in growth between winners, losers and drawing individuals within CM and SB treatments (Figure 5). In the CM treatment, winners grew significantly more than individuals who lost, but losers also grew more than those paired individuals who drew. In the SB treatment, winners grew significantly more than individuals who lost or drew, but there were no significant differences between losing individuals and drawing pairs (Figure 5).

Fig. 5. Growth of barnacles who ‘win’ (white), ‘lose’ (grey) or ‘draw’ (black) in a contest with their nearest neighbour in intraspecific competition treatments (Chthamalus montagui (CM) and Semibalanus balanoides (SB)). Post-hoc pairwise comparison outcomes are shown as letters above columns, where different letters indicates significant differences between groups (P < 0.05).

Table 3. Analysis of Variance (ANOVA) of barnacle growth in relation to outcome (win, lose, draw) and intraspecific treatments (SB and CM).

Effect of operculum size on growth over time

Growth was dependent on the initial size of the barnacle (at Time-0) and the outcome of the ‘competition’ between paired individuals (Table 4, P < 0.01). The ANCOVA revealed that the slopes were parallel but the intercepts (starting body size) were significantly different, indicating that the initial operculum size affected the outcome of the contest i.e. a larger individual was more likely to ‘win’ a contest (Figure 6). There was a negative relationship between growth and initial operculum size suggesting smaller individuals had a greater scope-for-growth than larger individuals.

Fig. 6. Relationship between initial operculum size (mm) and growth after 8 weeks for barnacles that ‘win’ (white circles, dashed line), ‘lose’ (grey circles, dotted line), and ‘draw’ (black circles, solid line). Significant regressions are shown (P < 0.05).

Table 4. Analysis of Covariance (ANCOVA) of barnacle growth in relation to initial size and outcome (win, lose, draw).

Variation in survival of barnacles in relation to time (Ti) and treatment (Tr)

There was a significant decrease in survival over time (Table 5, Ti, P < 0.001). Overall there was a 57% decrease in total barnacle survival after 8 weeks, although densities approached asymptote after 4 weeks suggesting a period of rapid early mortality of individuals (0 – <4 weeks) followed by relatively little mortality (Figure 7). There was no significant difference in survival between the three treatments over time (Table 5, Figure 7).

Fig. 7. Panel plot showing survival over time (weeks) for each treatment combination in 10 replicate quadrats: Mixed (CM + SB), CM (Chthamalus montagui) and SB (Semibalanus balanoides). A LOESS smoother with a span of 1 was fitted to aid visual interpretation.

Table 5. Linear mixed effects model summary testing the effect of time (Ti) and intraspecific treatment on survival.

Response Surface Methodology (RSM) for testing competition

Fitting a response surface model indicated that in mixed plots, the density of both CM and SB had a significant effect on the mortality of CM (Table 6, P < 0.05). Mortality was highest when (i) CM occurred at intermediate densities (~40 individuals per 25 cm2) in conjunction with high (~60 individuals per 25 cm2) densities of SB, or (ii) when densities of CM were low and SB were absent (Table 6, Figure 8). In contrast, model reduction revealed no significant effect of CM density or SB density on SB mortality (F 1,13 = 0.5983, P = 0.45) suggesting no negative effects of intra- or interspecific competition on SB at the patch scale.

Fig. 8. Response surface plot of Chthamalus montagui proportional mortality in 25 cm2 plots after 8 weeks.

Table 6. Response Surface Model Fitting (ANOVA) of Chthamalus montagui and Semibalanus balanoides barnacle survival in mixed treatments.

R 2 (CM) = 0.38.

DISCUSSION

The results of this study demonstrate that Semibalanus balanoides and Chthamalus montagui show strong spatial variation in operculum growth in areas where their distributions overlap in their intertidal. The composition of the patch – either a single species (CM or SB only) or a mixture of both species (CM + SB) – had a greater influence on growth than density. A comparison of the growth of all individuals in a patch with consideration of distance between individuals (nearest neighbour distance) suggested little or no competition was occurring between individuals, but comparisons at smaller spatial scales focusing on pairs of individuals neighbouring each other revealed clear ‘winners’ and ‘losers’ in terms of individual growth with few ‘draws’. An individual with a larger operculum from the outset of a paired competition tended to result in that individual ‘winning’ the growth competition, although larger individuals tended to grow less than smaller conspecifics. Comparing mortality rates at the patch level using a surface response method suggests that CM are undergoing competition, both with conspecifics and SB, whereas SB show no evidence of mortality related to competition.

Connell (Reference Connell1961a) investigated growth rates of intertidal barnacles, discovering growth rates of S. balanoides were greater than that of Chthamalus montagui. While direct interactions between S. balanoides and C. montagui have not been investigated, based on Connell's findings it was predicted that S. balanoides would have a higher operculum growth rate than C. montagui when co-existing. Treatments containing only a single species (CM or SB) indicated considerable variation in operculum growth depending on location. Differences in growth were not consistent across species, with growth of C. montagui greater at Site 2 over Site 1, and for S. balanoides, greater at Site 1 over Site 2. This differentiation may be due to relatively small-scale spatial variation in resource availability between sites located on Mount Batten shore. In previous studies, spatial heterogeneity of biotic and abiotic conditions that vary among rocks within a shore has been found to influence fluctuations in population growth (Fukaya et al., Reference Fukaya, Okuda, Hori, Yamamoto, Nakaoka and Noda2013). This is especially important in temperate regions, where seasonal fluctuations have been shown to affect growth; in the summer, population growth rates are strongly affected by regional-scale fluctuation, whereas in winter growth is affected more by rock-scale fluctuation (Fukaya et al., Reference Fukaya, Okuda, Nakaoka, Hori and Noda2010). Given Connell (Reference Connell1961a) carried out his study in spring/summer and this study was undertaken in the autumn, rock-scale fluctuations may have influenced the results. López et al. (Reference López, Ramírez, Guaitro and López2010) also showed variation in morphological structures (e.g. filtration and respiration) due to spatial-temporal fluctuations in biotic and abiotic factors including density and wave exposure. Such variation may also account for some of the variation in individual growth rates and explain the differences here.

Shinen & Navarrete (Reference Shinen and Navarrete2014) indicated that in the two barnacle species, Notochthamalus scabrosus and Jehlius cirratus, growth rates were largely consistent and independent of occurring together or separately. Certainly, in this study, growth rates of C. montagui at Site 1 was largely the same irrespective of whether it occurred in isolation or in combination with S. balanoides. In contrast, the growth of S. balanoides was marginally greater when occurring in isolation than when occurring with C. montagui although clear significant differences could not be determined. Connell (Reference Connell1961a) reported differences in growth between intraspecific and interspecific treatments, where he showed the average growth of C. montagui was higher when independent of S. balanoides (Connell, Reference Connell1961a). In contrast, at Site 2 the average growth of individuals in the interspecific treatment (CM + SB) was higher than when S. balanoides occurred on its own and the growth of C. montagui occurring in isolation was marginally lower, contradicting Connell's findings (Reference Connell1961a). Location, even at relatively small spatial scales, appears to play an important role in the effect of interspecific interactions (Sandford & Menge, Reference Sandford and Menge2001). Thus, the experimental area – a unique shore on the Isle of Cumbrae in the Firth of Clyde, Scotland – and level of replication used in Connell's study (Reference Connell1961a) may therefore be less representative of general interaction strength implications than previously thought. The mechanisms that alter growth rates in interspecific treatments remain unclear, but it is suspected that differences in microhabitats, which vary greatly in their degree of physical stress, may be related to differences of growth as this environmental heterogeneity creates distinct selection regimes (Schmidt & Rand, Reference Schmidt and Rand1999; Schmidt et al., Reference Schmidt, Bertness and Rand2015) and this should be explored further to tease apart differences.

Many studies have shown high densities of organisms can negatively affect growth if resources are limited due to increased competition (Barnes & Powell, Reference Barnes and Powell1950; Connell, Reference Connell1961a, Reference Connellb; Leslie, Reference Leslie2005; López et al., Reference López, Catalán, Barriga and López2014). Here, including nearest neighbour distance (NND) between all individuals had no effect on growth (although growth was highly variable) suggesting that resources were not limited in this instance at scale of a patch to the extent that competition was reduced and both species were able to co-exist (Gerwing et al., Reference Gerwing, Drolet, Hamilton and Barbeau2016). Higher densities of organisms has been shown to facilitate growth as the complex structure that is formed (e.g. mussel hummocks) can buffer individuals from physical disturbance, consumers and physiological stress (Bertness, Reference Bertness1989; Leslie, Reference Leslie2005) and elevate individuals exposing them to higher particle fluxes (Bertness et al., Reference Bertness, Gains and Yeh1998). Our results for C. montagui indicate some support for these mechanisms, with individuals on average exhibiting increased growth in areas of higher density.

The relative strength (importance) of intraspecific and interspecific competition has long been debated, with intraspecific competition often implicated as the main driver of negative effects (e.g. reductions in growth) through competition for resources between closely matched individuals (Moore, Reference Moore1935; Hixon et al., Reference Hixon, Pacala and Sandin2002; Begon et al., Reference Begon, Townsend and Harper2006). Here, there were significant differences in growth between neighbouring pairs of C. montagui and S. balanoides indicating clear winners and losers at small-spatial scales. As early as 1939, Moore demonstrated that the initial size of a S. balanoides could influence its growth. Here, initial operculum size was a good predictor of whether an individual would win, lose or draw a contest, with larger individuals tending to win a contest over smaller counterparts. There was also a difference in growth potential, with individuals with a smaller operculum size growing more than larger conspecifics within contest categories.

Differentiation between winners and losers helps to partition much of the variability in growth at the patch level. Interestingly, while there were relatively few ‘drawn’ contests, where they occurred, growth in C. montagui was greatly reduced in drawn contests such that growth was lower than the ‘loser’ in a contest where there was a clear winner and loser. In S. balanoides, growth of the ‘loser’ or those featuring in a ‘drawn’ contest exhibited similar but reduced growth in comparison to a ‘winning’ individual. These results suggest that competition occurs at the spatial scale of individuals, rather than at the patch level, requiring comparisons to be made at the scale of individuals if we are to be able to detect the effect of competition on patch dynamics. The differences in growth of C. montagui and S. balanoides as a result of individual contests suggest that there may be differences, perhaps in the morphology or physiological requirements of the two species, that alter the strength of mechanisms such as interference competition (Shinen & Navarrete, Reference Shinen and Navarrete2010) allowing S. balanoides to be more successful than C. montagui when occurring in close proximity with conspecifics. This explanation is supported by the comparison of mortality rates between species, which show negative effects of competition on CM, but not SB.

Interspecific competition can also have a negative effect on a species' survival (Connell, Reference Connell1961a; Wethey, Reference Wethey1983; Jenkins et al., Reference Jenkins, Murua and Burrows2008). Connell (Reference Connell1961a) showed direct competition between S. balanoides and C. montagui and indicated crowding was an important cause of death of C. montagui. Crowding can lead to an elongation of the calcareous and exoskeleton structures (Bertness et al., Reference Bertness, Gains and Yeh1998; López et al., Reference López, López, Burgosb, Arriagadac and González2007), especially in S. balanoides, which have been shown to grow tall, thin-walled and dependent on neighbours for structural support (Connell, Reference Connell1961a). Aggregations can be extremely fragile (Bertness, Reference Bertness1989) such that water motion (e.g. from wave exposure) can reduce barnacle survival (Connell, Reference Connell1961b; Gaylord, Reference Gaylord1999). Reductions in density, through mechanisms such as disturbance or predation (Knights et al., Reference Knights, Firth and Walters2012), can increase survival and reduce competition between or within species (Bertness, Reference Bertness1989; Bracewell et al., Reference Bracewell, Robinson, Firth and Knights2013). Here, no difference in the rate of survival of SB between intraspecific and interspecific treatments was found, suggesting that density of individuals in interspecific plots was sufficiently low to reduce the impacts of direct competition between the two species, increasing survivorship of both species within the patch.

In summary, this study suggests small-scale spatial distribution of organisms, the density and composition of species within a patch play an important role in determining the strength (or lack of) intraspecific and interspecific competition in intertidal communities. The results show clear variability in the outcome of contests between individuals of the same and different species related to these factors and indicates that an evaluation of the mechanism at the scale of a patch may not capture the effects of those processes effectively. The role of competition in affecting the structure and functioning of intertidal shores is clearly important, however, the mechanisms may not be as generic as previously thought.

ACKNOWLEDGEMENTS

We would like to thank Louis Evans and Tom Anderson for their assistance in the field and the technical support of Plymouth University. We would especially like to thank Dr Sarah Lane for valued discussions regarding animal contests.

References

REFERENCES

Barnes, H. and Powell, H.T. (1950) The development, general morphology and subsequent elimination of barnacle populations, Balanus crenatus and Balanus balanoides, after a heavy initial settlement. Journal of Animal Ecology 19, 175179.Google Scholar
Bates, D., Maechler, M., Bolker, B. and Walker, S. (2015) Fitting linear mixed-effect models using lme4. Journal of Statistical Software 67, 148.Google Scholar
Begon, M., Townsend, C. and Harper, J. (2006) Ecology: from individuals to ecosystems, 4th edition. Oxford: Blackwell Publishing, pp. 132163.Google Scholar
Bertness, M.D. (1989) Intraspecific competition and facilitation in a northern acorn barnacle population. Ecology 70, 257268.Google Scholar
Bertness, M.D., Gains, S.D. and Yeh, S.M. (1998) Making mountains out of barnacles: the dynamics of acorn barnacle hummocking. Ecology 79, 13821398.Google Scholar
Box, G.E.P., Hunter, J.S. and Hunter, W.G. (2005) Statistics for experimenters, 2nd edition. Hoboken, NJ: Wiley-Interscience.Google Scholar
Bracewell, S.A., Robinson, L.A., Firth, L.B. and Knights, A.M. (2013) Predicting free-space occupancy on novel artificial structures by an invasive intertidal barnacle using a removal experiment. PLoS ONE 8, 17.Google Scholar
Burnham, K.P. and Anderson, D.R. (2004) Multimodel inference understanding AIC and BIC in model selection. Sociological Methods and Research 33, 261304.Google Scholar
Burrows, M.T., Jenkins, S.R., Robb, L. and Harvey, R. (2010) Spatial variation in size and density of adult and post-settlement Semibalanus balanoides: effects of oceanographic and local conditions. Marine Ecology Progress Series 398, 207219.Google Scholar
Caro, A.U., Guiñez, R., Ortiz, V. and Castilla, J.C. (2011) Competition between a native mussel and a non-indigenous invader for primary space on intertidal rocky shores in Chile. Marine Ecology Progress Series 428, 177185.Google Scholar
Connell, J.H. (1961a) The influence of interspecific competition and other factors on the distribution of the barnacle Chthamalus stellatus. Ecological Society of America 42, 710723.Google Scholar
Connell, J.H. (1961b) Effects of competition, predation by Thais lapillus, and other factors on natural populations of the barnacle Balanus balanoides. Ecological Society of America 31, 62103.Google Scholar
Connell, J.H. (1983) On the prevalence and relative importance of interspecific competition: evidence from field experiment. American Naturalist 122, 661662.Google Scholar
Feldhamer, G.A., Dricamer, L.C., Vessey, S.H., Merritt, J.F. and Krajewski, C. (2007) Mammalogy: adaptation, diversity, ecology. Baltimore, MD: The Johns Hopkins University Press, pp. 481482.Google Scholar
Fukaya, K., Okuda, T., Hori, M., Yamamoto, T., Nakaoka, M. and Noda, T. (2013) Variable processes that determine population growth and an invariant mean-variance relationship of intertidal barnacles. Ecological Society of America 4, 48.Google Scholar
Fukaya, K., Okuda, T., Nakaoka, M., Hori, M. and Noda, T. (2010) Seasonality in the strength and spatial scale of processes determining intertidal barnacle population growth. Journal of Animal Ecology 79, 12701279.Google Scholar
Gause, G.F. (1934) Experimental analysis of Vito Volterra's mathematical theory of the struggle for existence. Science 79, 1617.Google Scholar
Gaylord, B. (1999) Detailing agents of physical disturbance: wave-induced velocities and accelerations on a rocky shore. Journal of Experimental Marine Biology and Ecology 239, 85124.Google Scholar
Gerwing, T.G., Drolet, D., Hamilton, D.J. and Barbeau, M.A. (2016) Relative importance of biotic and abiotic forces on the composition and dynamics of a soft-sediment intertidal community. PLoS ONE 11, 115.Google Scholar
Hart, S.P. and Marshall, D.J. (2009) Spatial arrangement affects population dynamics and competition independent of community composition. Ecological Society of America 90, 14851491.Google Scholar
Hixon, M.A., Pacala, S.W. and Sandin, S.A. (2002) Population regulation: historical context and contemporary challenges of open vs closed systems. Ecology 83, 14901508.Google Scholar
Inouye, B.D. (2001) Response surface experimental designs for investigating interspecific competition. Ecology 82, 26962706.Google Scholar
Jenkins, S.R., Murua, J. and Burrows, M.T. (2008) Temporal changes in the strength of density-dependent mortality and growth in intertidal barnacles. Journal of Animal Ecology 77, 573584.Google Scholar
Knights, A.M., Firth, L.B., Thompson, R.C., Yunnie, A.L.E., Hiscock, K. and Hawkins, S.J. (2016) Plymouth – a world harbour through the ages. Regional Studies in Marine Science. Available at http://dx.doi.org/10.1016/j.rsma.2016.02.002.Google Scholar
Knights, A.M., Firth, L.B. and Walters, K. (2012) Interactions between multiple recruitment drivers: post-settlement predation mortality overrides flow-mediated recruitment. PLoS ONE 7, e35096.Google Scholar
Knights, A.M. and Walters, K. (2010) Recruit-recruit interactions, density-dependent processes and population persistence in the eastern oyster Crassostrea virginica. Marine Ecology Progress Series 404, 7990.Google Scholar
Lenth, R.V. (2009) Response-surface methods in R, using rsm. Journal of Statistical Software 32, 117.Google Scholar
Leslie, H.M. (2005) Positive intraspecific effects trump negative effects in high-density barnacle aggregations. Ecology 86, 27162725.Google Scholar
López, B.A., Catalán, A.M., Barriga, D.A. and López, D.A. (2014) Morphological response of the exoskeleton in the intertidal barnacle Jehlius cirratus growing at different densities. Journal of Crustacean Biology 1, 16.Google Scholar
López, B.A., Ramírez, R.P., Guaitro, S.Y. and López, D.A. (2010) Interspecific differences in the phenotypic plasticity of intertidal barnacles in response to habitat changes. Journal of Crustacean Biology 30, 357365.Google Scholar
López, D.A., López, B.A., Burgosb, I.C., Arriagadac, S.E. and González, M.L. (2007) Consequences of base modification in hummocks of the barnacle Austromegabalanus psittacus. New Zealand Journal of Marine and Freshwater Research 41, 291298.Google Scholar
Lotka, A.J. (1925) Elements of physical biology. Baltimore, MD: Williams and Wilkins.Google Scholar
Moore, H.B. (1935) The biology of Balanus balanoides. IV. Relation to environmental factors. Journal of the Marine Biological Association of the United Kingdom 20, 279307.Google Scholar
Moore, H.B. (1939) The biology of Balanus balanoides. I. Growth rate and its relation to size, season and tide level. Journal of the Marine Biological Association of the United Kingdom 17, 851868.Google Scholar
Neal, D. (2004) Introduction to population biology. Cambridge: Cambridge University Press, pp. 6870.Google Scholar
O'Brien, R.M. (2007) A caution regarding rules of thumb for Variance Inflation Factors. Quality and Quantity 41, 673690.Google Scholar
R Development Core Team (2016) R: a language and environment for statistical computing. Available at http://www.r-project.org/Google Scholar
Reece, J.B., Urry, L.A., Cain, M.L., Wasserman, S.A., Minorsky, P.V. and Jackson, R.B. (2011) Campbell biology. Boston, MA: Pearson Education, pp. 12001300.Google Scholar
Sandford, E. and Menge, B.A. (2001) Spatial and temporal variation in barnacle growth in a coastal upwelling system. Marine Ecology Progress Series 209, 143157.Google Scholar
Schmidt, P.S., Bertness, M.D. and Rand, D.M. (2015) Environmental heterogeneity and balancing selection in the acorn barnacle Semibalanus balanoides. The Royal Society 282(1804), 379383.Google Scholar
Schmidt, P.S. and Rand, D.M. (1999) Intertidal microhabitat and selection at MPI: interlocus contrasts in the northern acorn barnacle, Semibalanus balanoides. Evolution 53, 135146.Google Scholar
Schneider, C.A., Rasband, W.S. and Eliceiri, K.W. (2012) NIH Image to ImageJ: 25 years of image analysis. Nature Methods 9, 671675.Google Scholar
Shinen, J.L. and Navarrete, S.A. (2010) Coexistence and intertidal zonation of chthamalid barnacles along central Chile: interference competition or a lottery for space? Journal of Experimental Marine Biology and Ecology 392, 176187.Google Scholar
Shinen, J.L. and Navarrete, S.A. (2014) Lottery coexistence on rocky shores: weak niche differentiation or equal competitors engaged in neutral dynamics? American Society of Naturalists 183, 342362.Google Scholar
Siepielski, A.M. and McPeek, M.A. (2010) On the evidence for species coexistence: a critique of the coexistence program. Ecological Society of America 91, 31533164.Google Scholar
Southward, A.J. (1976) Taxonomic status and distribution of Chthamalus stellatus (Cirripedia) in Northeast Atlantic region – with a key to common intertidal barnacles of Britain. Journal of the Marine Biological Association of the United Kingdom 56, 10071028.Google Scholar
Strong, D.R., Simberloff, D., Abele, L. and Thistle, A. (1984) Ecological communities: conceptual issues and evidence. Princeton, NJ: Princeton University Press.Google Scholar
Tilman, D. (1982) Resource competition and community structure. Princeton, NJ: Princeton University Press, pp. 5175.Google Scholar
Volkenborn, N., Hedtkamp, S.I.C., van Beusekom, J.E.E. and Reise, K. (2007) Effects of bioturbation and bioirrigation by lugworms (Arenicola marina) on physical and chemical sediment properties and implications for intertidal habitat succession. Estuarine, Coastal and Shelf Science 74, 331343.Google Scholar
Volterra, V. (1926) Fluctuations in the abundance of a species considered mathematically. Nature 118, 558560.Google Scholar
Wethey, D.S. (1983) Intrapopulation variation in growth of sessile organisms: natural populations of the intertidal barnacle Balanus balanoides. Oikos 40, 1423.Google Scholar
Ying, Z., Liao, J., Wang, S., Hui, L. and Liu, Y. (2014) Species coexistence in a lattice-structured habitat: effects of species dispersal and interactions. Journal of Theoretical Biology 359, 184191.Google Scholar
Figure 0

Fig. 1. Geographic location of the study in Plymouth, UK. Inset: the intertidal rocky shore at Mount Batten (50°21′N 4°07′W).

Figure 1

Fig. 2. Schematics and photograph of 5 cm2 manipulative treatment plots used to test for evidence of intraspecific (A, B – Chthamalus montagui only) and interspecific (C) competition and (D) the response surface experimental design with example species densities in each patch at Time-0 as recommended by Inouye (2001).

Figure 2

Fig. 3. Mean (±SE; N = 5) log growth of barnacle opercula in three treatments: CM + SB (mixed species plots), CM (Chthamalus montagui only) and SB (Semibalanus balanoides only) at two sites after 8 weeks. Post-hoc pairwise comparison results are shown as letters above columns, where different letters indicates significant differences between ‘treatment × site’ combinations (P < 0.01).

Figure 3

Table 1. ANOVA of operculum growth (mm) in relation to site and treatment at Mount Batten, Plymouth, in 2014.

Figure 4

Fig. 4. Scatterplot of operculum growth in individual barnacles after 8 weeks in relation to distance (mm) from its nearest neighbour. Different treatment combinations are shown: Semibalanus balanoides + Chthamalus montagui (white); Chthamalus montagui only (grey); and Semibalanus balanoides only (black).

Figure 5

Table 2. Analysis of covariance (ANCOVA) of barnacle operculum growth in relation to nearest neighbour distance and treatment (SB only; CM only; CM + SB).

Figure 6

Fig. 5. Growth of barnacles who ‘win’ (white), ‘lose’ (grey) or ‘draw’ (black) in a contest with their nearest neighbour in intraspecific competition treatments (Chthamalus montagui (CM) and Semibalanus balanoides (SB)). Post-hoc pairwise comparison outcomes are shown as letters above columns, where different letters indicates significant differences between groups (P < 0.05).

Figure 7

Table 3. Analysis of Variance (ANOVA) of barnacle growth in relation to outcome (win, lose, draw) and intraspecific treatments (SB and CM).

Figure 8

Fig. 6. Relationship between initial operculum size (mm) and growth after 8 weeks for barnacles that ‘win’ (white circles, dashed line), ‘lose’ (grey circles, dotted line), and ‘draw’ (black circles, solid line). Significant regressions are shown (P < 0.05).

Figure 9

Table 4. Analysis of Covariance (ANCOVA) of barnacle growth in relation to initial size and outcome (win, lose, draw).

Figure 10

Fig. 7. Panel plot showing survival over time (weeks) for each treatment combination in 10 replicate quadrats: Mixed (CM + SB), CM (Chthamalus montagui) and SB (Semibalanus balanoides). A LOESS smoother with a span of 1 was fitted to aid visual interpretation.

Figure 11

Table 5. Linear mixed effects model summary testing the effect of time (Ti) and intraspecific treatment on survival.

Figure 12

Fig. 8. Response surface plot of Chthamalus montagui proportional mortality in 25 cm2 plots after 8 weeks.

Figure 13

Table 6. Response Surface Model Fitting (ANOVA) of Chthamalus montagui and Semibalanus balanoides barnacle survival in mixed treatments.