Original Research ARTICLE
Seaweeds in Two Oceans: Beta-Diversity
- 1Department of Biodiversity and Conservation Biology, University of the Western Cape, Bellville, South Africa
- 2South African Environmental Observation Network, Elwandle Coastal Node, Port Elizabeth, South Africa
- 3Biological Sciences Department and Marine Research Institute, University of Cape Town, Rondebosch, South Africa
- 4Seaweed Unit, Fisheries Branch, Department of Agriculture, Forestry and Fisheries, Roggebaai, South Africa
Several species assembly mechanisms have been proposed to structure ecological communities. We assess the biogeography of seaweeds along 2,900 km of South Africa's coastline in relation to a thermal gradient produced by the Agulhas Current, and contrast this with the environmental structure created by the Benguela Current. We subdivided the coastline into “bioregions” to examine the regional patterning. To investigate the assembly mechanisms, we decomposed Sørensen's β-diversity into “turnover” (βsim) and “nestedness-resultant” (βsne) dissimilarities, and used distance-based redundancy analysis (db-RDA) to relate them to the Euclidean thermal difference, dE, and geographical distance. Moran's eigenvector maps (MEM) were used as an additional set of spatial constraints. Variation partitioning was then used to find the relative strengths of thermal and spatially-structured thermal drivers. Spatial and environmental predictors explained 97.9% of the total variation in βsim and the thermal gradient accounted for 84.2% of this combined pool. βsim was the major component of overall β-diversity in the Agulhas Current region, suggesting niche influences (environmental sorting) as dominant assembly process there. The much weaker thermal gradient in the Benguela Current-influenced region resulted in a high amount of βsne that could indicate neutral assembly processes. The intensification of upwelling during the mid-Pliocene 4.6–3.2 Ma (i.e., historical factors) were likely responsible for setting up the strong disjunction between the species-poor west coast and species-rich south and east coast floras, and this separation continues to maintain two systems of community structuring mechanisms in the Atlantic and Indian Ocean influenced sides of South Africa.
The assembly processes that structure biodiversity across a range of scales form the theme of macroecology (Chave, 2013). This paper deals with the composition of seaweed assemblages along the ~2,900 km South African coastline, including the identification, description and explanation for the spatially-structured patterns at scales from 100s to 1,000s of kilometers. Inshore conditions along this coastline range from cool through warm temperate to fringe tropical (Bolton and Anderson, 2004; Bolton et al., 2004), and are influenced by two major ocean currents in two oceans that set up a strong thermal gradient along the shore (Smit et al., 2013). We investigate how a gradient in seaweed community composition may have arisen in response to this temperature gradient (e.g., Qian and Ricklefs, 2007), and the extent to which deterministic, niche-based species assembly processes could have contributed toward the assembly of the seaweed flora of the region.
Biodiversity may be viewed in terms of α-, γ- and β-diversity, with the former two referring to “local” and “regional” diversity, respectively. β-diversity as defined by Whittaker (1960, 1972) is a measure of variation in species composition from place to place and is comprised of two processes (Baselga, 2012): the replacement of species independently of the difference in species richness (called “turnover,” βsim) and a term that considers this difference in species richness called (“nestedness-resultant,” βsne). The relative contribution of these two processes is not immediately evident from species dissimilarities, and yet such considerations should be implicit in macroecological studies. The ability to decompose β-diversity into turnover and nestedness-resultant components has become possible within the last decade (Baselga, 2010), but the use of this form of β-diversity partitioning has not yet widely permeated the phycological literature (or even in marine studies more broadly; Anderson et al., 2013). Cognisance of spatial scaling should also be deeply rooted in the study of biodiversity more generally, but such questions have only recently begun to be addressed (Barton et al., 2013).
Since studies of β-diversity consider the variation in species assembly processes from place to place, and necessitate an understanding of mechanisms that drive these processes (Davidar et al., 2007), studying β-diversity of adjacent coastal marine bioregions or marine provinces (Spalding et al., 2007) could provide deeper insight into how such processes operate along gradients (Qian and Ricklefs, 2007). Few studies dealing specifically with β-diversity exist for marine biota, especially over spatial areas in excess of 1,000s of kilometers. Recent examples include studies on fish distribution (Zintzen et al., 2011; Anderson et al., 2013), the biogeography of macroalgae along 6,600 km of Australian coastline (Leaper et al., 2011), and a study on deep-sea bivalves (McClain et al., 2011). All four have β-diversity as a primary interest. However, macroecological studies in the terrestrial realm have yielded the most comprehensive insights into the drivers responsible for assembling species into communities (Whittaker, 1960; Davidar et al., 2007; Qian and Ricklefs, 2007; Soininen et al., 2007). What these studies showed is that, perhaps universally, high rates of species turnover are associated with steep environmental gradients, such as those occurring along mountain slopes or across latitudes. In the ocean, the biggest driver latitudinally seems to be temperature (Tittensor et al., 2010; Straub et al., 2016; Stuart-Smith et al., 2017); other drivers such as nutrients, salinity, turbidity, wave action and photoperiod may also structure species composition, but these are likely to depend more on the regional context, and some, such as nutrients (e.g., Waldron and Probyn, 1992), are often strongly correlated with a more distal variable, such as temperature.
Recent decades have seen temperature emerging as a readily available environmental driver that's able to resolve the ocean's heat content across a multitude of temporal and spatial scales. Temperature is also very useful as a predictor of species range limits due to it being one of the key output variables of coupled global ocean-atmosphere circulation models used to project the physical milieu of the future Earth that has experienced anthropogenic climatic change (IPCC, 2013), hence making it the focus of our present study. The controlling effect of seawater temperature on the survival and reproduction of benthic organisms and patterns in the evolution and ecology of biological assemblages at regional scales is well known (van den Hoek, 1982b; Breeman, 1988; Blanchet et al., 2008; Broitman et al., 2008; Byrne et al., 2009; Verbruggen et al., 2009; Wieters et al., 2009; Couce et al., 2012; Potts et al., 2014). It influences the species composition of the biota associated with various thermal zones, which forms patterns on global, regional and local scales (Tittensor et al., 2010; Spalding et al., 2012). At the global scale, ocean currents and well-characterized latitudinal solar heat flux gradients maintain the ocean's thermal regime. At local scales (<10 km) at the land/sea margin, however, additional less well-known physical phenomena contribute to thermal patterns and dynamics that may differ markedly from those at the mesoscale (i.e., spatial scales > 50 km), and the different properties of the temperature regime which set up the biogeographical patterns are less well understood there—this is especially true at regional and local scales near the land where satellite data poorly reflect reality (Smit et al., 2013). Most studies that recognize temperature as a major driver of species distribution take the annual mean temperature (Tittensor et al., 2010); fewer recognize the importance of variability and range (e.g., Couce et al., 2012; Tyberghein et al., 2012). Detailed laboratory culture studies with seaweed species, for example, closely link biogeographical distribution limits with maximum and minimum monthly mean temperatures (van den Hoek, 1982b; Breeman, 1988). The time-integrated “thermal environment” is an amalgam of various statistical properties derived from a multitude of instantaneous temperature recordings, and it is pertinent to question if one or a few of these properties have an overriding imprint on the species assembly.
Understanding the processes influencing species turnover and geographical range is of practical importance to spatial biodiversity planning. This is especially relevant in today's warming world because it will allow us to project future biodiversity responses. Focusing conservation efforts on areas of high β-diversity will ensure the preservation of a diversity of species and their environmental niches. Peaks in β-diversity along ecoclines may indicate boundaries between bioregions, highlighting regions where communities are poised at their environmental limits, thus aligning with long-term monitoring surveys that aim to project ecosystem responses to climatic change.
Using a detailed data set of seaweed presence and absence records coupled with coastal in situ seawater temperature climatologies that are able to resolve the coastal zone, we investigated the thermal properties and species composition of 58 coastal sections spaced around the South African coastline. The coastline is broadly influenced by two major ocean currents, the Benguela Current and the Agulhas Current. One is an eastern boundary upwelling system and defines a cool temperate environment, and the other a western boundary current driving meridional transport of sub-tropical water toward the tip of Africa. We were primarily interested in establishing how these ocean currents, of which the effect at the coast can readily be measured as gradients in thermal properties, may have influenced the species assembly processes operating along an approximately 2,900 km long coastline. To do so, we first link matrices of species dissimilarity to environmental distance matrices and Moran's eigenvector maps (MEM) using distance-based redundancy analysis, and then apply variance partitioning (Peres-Neto and Legendre, 2010) to determine the relative contributions of thermal metrics and other spatially-organized drivers of seaweed community composition (as seen in the β-diversity components, βsim and βsne). Second, we examine how the thermal and spatial structuring agents operate at smaller spatial scales within marine provinces of the region by undertaking a more detailed analysis of the various distance matrices, and also ask which of the thermal properties we examined were most influential in setting up the patterns that emerged. Our findings are congruent with existing knowledge (e.g., Stephenson, 1948; Lombard et al., 2004; Spalding et al., 2007) and add a more nuanced understanding of the processes responsible for structuring the seaweed biodiversity in the sub-region.
2.1. Species Data and Explanatory Variables
We used three sets of data in this analysis. The first comprises distribution records (presence/absence) of 846 macroalgal species belonging to the Divisions Ochrophyta, Rhodophyta and Chlorophyta within each of 58 × 50 km-long sections of the South African coast (mentioned in bold font in the text and listed in Appendix A in Supplementary Materials, Table 1). The seaweed data represent ca. 90% of the known seaweed flora of South Africa in the intertidal and shallow subtidal, but excludes some very small and/or very rare species for which data are insufficient. The data are from verifiable literature sources and our own collections, assembled from information collected by teams of phycologists over three decades (Bolton, 1986, 1996; Stegenga et al., 1997; Bolton and Stegenga, 2002; De Clerck et al., 2005).
The second is a dataset of in situ coastal seawater temperatures (Smit et al., 2013) derived from daily measurements over up to 40 years. The thermal data set was used as the first set of explanatory variables. The following statistical properties were entered into the analysis: the means for the year (annMean), February (febMean, Austral summer) and August (augMean, Austral winter); the annual standard deviation (SD) around the mean for the year, February and August (annSD, febSD and augSD, respectively); and the annual thermal range between the mean temperature of the warmest and coldest months (annRange), and the mean range of February and August temperatures (febRange and augRange).
The third set of explanatory variables was generated to represent the spatial connectivity among coastal sections. Because of the strong environmental gradients along the shore, species community composition should be spatially organized (including the modeled spatially structured environmental variables, exogenous spatial variables and autocorrelation), and this was accounted for in the analysis. To this end we produced a section–section connectivity matrix based on a minimum spanning tree (MST) topology, starting from a geographical distance matrix. This topology focuses on relationships between neighboring sections and discards connections that are further away. In the case of the coastline data, sections are connected only to other sections that are below a certain truncation distance. This largely resulted in a “string-of-beads” series of section–section connections; for example, section 4 is directly connected to only sections 3 and 5; sections at the termini of the coastal “string” of sections are each connected to only one other section. We generated Moran's eigenvector maps (MEM; Dray et al., 2006, 2012) from this connectivity matrix through a principal coordinates analysis (PCoA) using the PCNM function of the PCNM package in R 3.3.3 (R Core Team, 2017), and kept the MEMs with positive spatial correlation. The MEMs are completely orthogonal and represent the spatial structures over the full range of scales from 50 to 2,900 km. Large eigenvectors represent broad spatial scales while smaller ones cover finer features. The spatial data were used as the second set of explanatory variables in multiple regression type analyses (Dray et al., 2012). Details and code are provided in Appendix B (Supplementary Materials) and may be accessed online at https://github.com/ajsmit/Seaweeds_in_Two_Oceans.git.
2.2. Distance-Based Redundancy Analysis and Variance Partitioning
Using the distance-based redundancy analysis (db-RDA; Minchin, 1987) implemented with the capscale function of the R package, vegan (Oksanen et al., 2017), we explored the role of the thermal and spatial descriptors in structuring the seaweed communities across the 58 coastal sections. To represent the biotic data, we decomposed Sørensen's dissimilarity (βsør) into its “nestedness-resultant” (βsne) and “turnover” (βsim) components (Baselga, 2010) using the betapart package (Baselga et al., 2017). This approach was necessary because β-diversity is strongly coupled with α-diversity, and it allowed us to make inferences about the possible drivers of β-diversity. Turnover refers to processes that cause communities to differ due to species being lost and/or gained from section to section, i.e., the species composition changes between sections without corresponding changes in α-diversity. The nestedness-resultant component implies processes that cause species to be gained or lost, and the community with the lowest α-diversity is a subset of the richer community.
The assessment of the biotic ordination within the context of the underlying environmental properties proceeded using the environmental variables' z-scores. To determine the descriptors that best describe the patterns in the seaweed dissimilarity data, we first applied full (global) db-RDAs using the complete sets of thermal and spatial variables, separately for each set. Using the forward selection procedure implemented in the packfor package for R (Blanchet et al., 2008), we reduced the number of variables in each set and retained only those that best fit the biotic data. Forward selection prevents the inflation of the overall type I error and reduces the number of explanatory variables used in the final model, which improves parsimony. The reduced set of thermal variables retained collinear variables (Graham, 2003), which were identified and removed using variance inflation factors (VIF; Dormann et al., 2013). This was not necessary for the MEMs as they are orthogonal by definition. The remaining significant orthogonal thermal and spatial variables were then regressed with βsim and βsne and final db-RDA models produced. The computation of db-RDAs was followed by permutation tests of the adjusted R2 to assess the significance of constraints (thermal and spatial descriptors). Lastly we undertook variance partitioning (Peres-Neto et al., 2006; Peres-Neto and Legendre, 2010) between the environmental (thermal) and spatial predictors using the varpart function in the vegan package. Refer to Appendix B in Supplementary Materials for more information about the methodology and for the R code underlying the analysis.
2.3. Analyses of Dissimilarity and Distance Matrices
We then explored regional patterns in β-diversity. Because connectivity between sections is constrained by their location along the shore and thus direct distances between sections do not apply, the total distance between a pair of arbitrary sections is the cumulative sum of the great circle distances between each consecutive pair of intervening sections along the coast (this is in fact encapsulated by the connectivity matrix used in the PCNM analysis, above). Plots showing the relationship of βsim and βsne with distance are limited because they do not provide a geographical context. To overcome this problem, we used a “network graph” to show spatial relationships in regional species dissimilarity. See Appendix C in Supplementary Materials for details and R code.
The last step of our analysis was applied to the four bioregions recognized for South Africa (Bolton and Anderson, 2004), namely the Benguela Marine Province (BMP; 1–17), the Benguela-Agulhas Transition Zone (B-ATZ; 18–22), the Agulhas Marine Province (AMP; 19–43/44) and the East Coast Transition Zone (ECTZ; 44/45–58). To this end, we calculated an Euclidean distance matrix that encapsulated all pairwise differences between coastal sections for each of the thermal metrics highlighted in the db-RDA (above) using the vegan package; these are called thermal differences (dE) throughout. We then correlated βsim and βsne with geographical distance and dE and provided matching plots in which the four bioregions were color-coded to discern bioregional affiliations and differences. Together these analyses were able to capture β-diversity at two spatial scales: among sections within bioregions, and among all sections for the whole country.
3.1. The Thermal Environment and Species Richness
South Africa's annual mean coastal water temperature ranged from 12.0 ± 0.9°C (mean ± SD) at its north-western limit near the Namibian border (1) to 24.0 ± 1.9°C on the east coast near the Mozambican border (58) (Figure 1). The global latitudinal gradient of diminishing temperature with increasing latitude was seen only along the east coast where the annual mean temperature decreased from ca. 24.5°C near 58 to 17.5°C around 39. The alongshore thermal gradient for this 950 km stretch of coastline was ca. 0.7°C per 100 km, with steeper gradients near 54. The latitudinal gradient largely reversed in direction along the west coast (1–16), i.e., temperatures became slightly cooler further north.
Figure 1. Profile of mean annual temperature along the coast of South Africa. Profiles are indicated as a function of geographical position on a coastal map (top) and as a function of distance away from section 1 (bottom). The latter visualization also indicates the long-term minimum (mean August) and maximum (mean February) temperatures as a gray shaded area around the annual mean temperature. The inset shows the species richness of macroalgae along the coast. Note that although a detailed temperature profile is displayed here, further analyses in this paper proceed with temperature data interpolated to the 58 sections for which seaweed diversity data are available.
On average, these data indicated an increase in inshore annual mean temperatures from west to east (1–58) of 12.1–24.4°C (range: 12.3°C). In February the thermal range was 13.7°C, while in August it was 10.5°C. In August the west–east temperature transition was smooth whereas in February substantial warm fluctuations in the mean monthly temperature were observed in embayments such as 13 and the False Bay sections from 17–18, and 28 and some sections around 35/36 and Algoa Bay from 34–36. In summer the mean monthly temperature gradient steepened between 19–28, and thereafter decreased eastwards along the coast from 33–34.
The number of species within the BMP was low (Figure 1, inset), and many northern sections had fewer than 150 species, and it rose significantly in the warmer section around 12/13, and the southern sections around the Cape Peninsula (16/17). Thereafter, richness increased markedly in the B-ATZ, the AMP and the ECTZ. The highest number of species in any one section was 340 (section 39 near the eastern end of the AMP).
3.2. Environmental Correlates of Seaweed Diversity
db-RDA, forward selection and the assessment of VIF retained augMean, febRange, febSD and augSD as the most parsimonious descriptors of βsim with an adjusted R2 of 0.885, explaining 89.8% of the variation (global permutation test on final model: d.f. = 4, F = 110.16, p = 0.001). The model consisted of two significant canonical axes: CAP1 and CAP2 explained 73.3 and 14.9% of the variation, respectively. The biplot scores (vectors) showed that augMean was heavily loaded along CAP1 and the metrics related to variation around the mean, i.e., febRange and febSD, strongly influenced βsim along CAP2 (Figure 2). Plots of the “lc” scores on geographic axes are given in Figure 3. The scores representing CAP1 increased gradually along the shore from west to east, reflecting the pervasive influence of augMean as coastal sections changed from cool to warm temperate through to sub-tropical thermal regimes. CAP2 site scores were lowest along the southern sections. βsne was only influenced by annMean (Figure 2) along CAP1 that explained 20.3% of the total variation (R2 = −0.140, d.f. = 1, F = −6.018).
Figure 2. db-RDA biplots of lc scores for βsim and βsne constrained by the thermal and spatial MEM variables. The only two canonical axes shown are CAP1 and CAP2 since these capture the bulk of the inertia present in the ordinations. The constraining vectors were selected during the db-RDA steps.
Figure 3. db-RDA site scores for βsim and βsne on geographic coordinates. The site scores indicate the major gradients captured by the environmental and spatial (MEM) constraints. Gray circles indicate negative site scores.
The db-RDA analysis procedure retained 17 significant MEMs that fully encapsulated the spatial dependence within βsim (d.f. = 18, F = 84.055, p = 0.001), resulting in an adjusted R2 of 0.963 and accounting for all of the variation. Fifteen canonical axes were produced of which five were significant. The first one alone explained 77.1% of variance with the remaining axes accounting for ca. 1% or less of the inertia. MEM2, MEM3 and MEM5 were most strongly loaded along CAP1 (Figure 3). These MEMs caused the sections belonging with the BMP to separate out from all other sections, and also for the northern sites of the ECTZ to diverge from the sections and bioregions further south. The scales of spatial dependence captured by these MEMs could all be considered to be broad-scaled. Two canonical axes comprised of two significant MEMs were able to explain 87.3% of the variation in βsne (R2 of 0.437, d.f. = 4, F = 12.06, p = 0.001; Figure 3). CAP1 consumed 79.0% of the inertia due to MEM1 and MEM5.
The partitioning of the variance associated with the seaweed community along the coast (Table 1) was explained jointly by the thermal and spatial variables selected in the preceding db-RDAs. Combining the thermal and spatial predictors (fractions [E+S]) allowed the model to capture 97.9% of the total βsim variance (F = 191.56, p = 0.001), with a residual variance of 2.1%. The thermal variables on their own (i.e., those that are spatially unstructured; fraction [E|S]) were able to account for only 1.8% of the total variation (F = 20.506, p = 0.001), but including some spatially structured thermal properties, [E], raised the proportion of explained variation to 84.2% (F = 110.16, p = 0.001). Pure spatial patterning (i.e., in the absence of temperature influences, perhaps with exogenous environmental influences or autocorrelation; fraction [S|E]) drove 13.7% of the species variation (F = 23.649, p = 0.001), and adding some thermal influences together with spatial descriptors, [S], increased this to 96.1% (F = 80.731, p = 0.001). Turning now to βsne, we see that our explanatory variables were less successful in capturing the variation. The spatial variables, [S], and the spatial plus thermal variables, [E+S], were able to account for 67.8% (F = 12.06, p = 0.001) and 71.4% (F = 9.077, p = 0.001) of the variation, respectively.
3.3. Pairwise Dissimilarities
Network graphs show the spatial relationships of βsim (Figure 4). βsim clearly highlighted the effect of the sharp change in α-diversity between the BMP and the B-ATZ. Sections in the B-ATZ retained similarities with sections as far east as 33 within the AMP. Eastwards from 30 similarities with sections within the ECTZ became apparent, but they generally did not extend past 43. Aside from a very low similarity with sections at the eastern extent of the AMP, the ECTZ sections retained similarities with sections within the same biogeographical province only over very short distances, and again it highlighted the high β-diversity in this region.
Figure 4. Network graphs of pairwise βsim showing sections that are similar to one-another in species composition. Dissimilarity indices range from >0.0 to 0.5 (possible dissimilarities range from 0.0 to 1.0). The 58 sections in Appendix A (Supplementary Materials) appear as dots at the vertices of the network graph, with those belonging with the “active” bioregion shown in white. The pairwise dissimilarities are shown by the colored lines, with blacker lines indicating lower dissimilarity indices (species composition more similar) and redder ones higher dissimilarity indices (species composition more dissimilar).
The overall and regional mean values for the three measures of pairwise β-diversity are presented in Table 2. The overall Sørensen β-diversity (βsør, 0.496 ± 0.287) was larger than that of the bioregions, and only a small fraction of it was comprised of nestedness-resultant β-diversity. Of the four bioregions, the ECTZ had the highest βsør (0.259 ± 0.157). At the bioregional scale, it is important to note that the nestedness-resultant component was about four times larger for the BMP (βsne/βsør = 0.581) than that of the other three (βsne/βsør ranged from 0.097 to 0.170).
Plots of βsim and βsne among all possible section pairs indicated clear differences among the four bioregions in their relationships with geographic and thermal distances (Figure 5). The BMP showed a weak relationship between βsim and distance (r2 = 0.052; regression statistics in Table 3) since much of the compositional variation between sections in this region was due as much to nestedness-resultant β-diversity as it was to turnover (Figure 5) as noted above (see βsne/βsør in Table 2). Within the ECTZ (~900 km long coastline) and the B-ATZ (only 170 km long) the rates were moderate at β = 0.079 (r2 = 0.936) and high at β = 0.109 (r2 = 0.658) per 100 km, respectively. In the AMP it was lower at β = 0.029 per 100 km (r2 = 0.834). Higher rates (β) indicate that communities turned over more rapidly per unit distance of coastline; furthermore, this also provided strong evidence that β-diversity was structured along environmental gradients. β-diversity was less influenced by changes in species numbers between sections in the B-ATZ, AMP and the ECTZ, as these two marine provinces were characterized by a relatively even number of species and hence had low β-diversities attributed to nestedness (Figure 5). βsim expressed with respect to the augMean thermal distance showed similarly steep slopes for three of the bioregions (βs ranging from 0.290 to 0.350, with r2 > 0.605), with that for the BMP about a two-thirds lower (Figure 5; Table 3). With respect to febRange, a significant relationship with βsim existed only for the two transitional areas, the B-ATZ and the ECTZ (r2 of 0.548 and 0.583, respectively; Figure 5). The steepness of the relationship between βsim and the latter thermal metric was lower than that seen with augMean. Concerning febSD, this relationship was steepest for the B-ATZ (β = 0.103, r2 = 0.276) and then the ECTZ (β = 0.082, r2 = 0.310; Figure 5); the same general trend held for augSD, with r2 = 0.310 and r2 = 0.276 for the ECTZ and B-ATZ, respectively. The pattern of βsne with geographic and thermal distance was generally significant but very poor (low r2-values) and with weak gradients (β) (Figure 5). The notable outcome there was that it was the BMP where βsne was strongest even though its r2 values were weak (r2 = 0.205 and 0.164 for the relationship with geographical distance and thermal distance, respectively).
Figure 5. Plots of the turnover (βsim) and nestedness-resultant (βsne) forms of β-diversity as a function of geographical or thermal distance between the coastal sections. The influential dE variables (augMean, febRange, febSD and augSD for βsim, and annMean for βsne) were determined in the db-RDA procedure. βsim, βsne, dE and geographical distance were calculated as differences between coastal section pairs, and data points representing section pairs falling between bioregions are colored yellow and labeled “out.” The strengths and gradients of the regression lines indicated in the panels are presented in Table 3.
Table 3. Linear regressions of pairwise βsim and βsne as a function of the geographical or thermal distance between the section pairs.
In this paper, we considered the drivers of seaweed β-diversity at the scale of bioregions (marine provinces) nested within a ~2,900 km stretch of coastline. We examined β-diversity using a numerical ecological approach that permits assessing its spatial pattern independently of the influence of α-diversity (Baselga et al., 2017), which in our study displays major regional differences. There are fewer species of macroalgae on the Benguela Current dominated west coast, as is also shown by studies on invertebrates and fishes (Thandar, 1989; Turpie et al., 2000). There is a similar species richness per coastal section on the south coast and the east coast. This was not the case up until the early 2000s (Bolton and Stegenga, 2002), when a more restricted data set showed more species per section on the south than east coast. This difference disappeared during the mid-2000s, when major collections on the east coast filled out the species lists available for the region (De Clerck et al., 2005). Continuous sampling over three decades has now resulted in a robust α-diversity data set, and it is still different from the animal data: fish and invertebrates become more abundant per section on the east coast (toward the tropics) but the seaweeds do not. This matches previously published global seaweed diversity patterns, which suggested that biogeographical seaweed diversity does not peak in the tropics. Globally, seaweed species diversity at a biogeographical scale is highest in warm temperate regions (Bolton, 1994; Santelices et al., 2009). Examining spatial trends in α-diversity on its own offers an incomplete understanding of the structuring of biodiversity at regional scales, and more insight may be obtained by looking at species richness (i.e., lists of species, such as which we use in this study) across two scales, sensu Whittaker (1972)—i.e., α-diversity (species richness at the local scale, such as per coastal section) and γ-diversity (total species richness in all coastal sections). In this broad sense, the two scales of species richness are not independent, and the loss or gain in species from section to section will be seen in β-diversity. In this work, we employ Baselga's (2010) framework to represent overall β-diversity as two components, one of which excludes the effect of regional differences in α-diversity.
With the exception of the some Australian studies (Smale et al., 2010, 2011; Waters et al., 2010; Leaper et al., 2011; Wernberg et al., 2013), one in Europe (Tuya et al., 2012) and another two in the northwest Indian Ocean/Red Sea region (Schils and Wilson, 2006; Issa et al., 2016), studies of this scale and nature have so far been infrequently seen for marine macroalgae. Of the above-mentioned studies, only two specifically considered β-diversity (Leaper et al., 2011; Issa et al., 2016). Our data focus on the taxonomic representivity of a region by aggregating species occurrence records within 50 km long sections of coastline, and are blind to the effects of small-scale habitat heterogeneity (e.g., as seen in Smale et al., 2010). For this reason we cannot infer influences of stochastic processes and other aspects of environmental heterogeneity at scales of <50 km; as such we excluded drivers that may be only of local relevance, such as salinity, turbidity (light), and wave action. Rather, our data emphasize broad biogeographic patterns (Lawton, 1999), with induced spatial dependence (Peres-Neto and Legendre, 2010) emerging as the main species assembly process east of the Cape Peninsula. This points to deterministic niche influences that correlate with environmental drivers that function over a broad spatial scale. These environmental and biotic patterns reflect the nature of the two dominant ocean currents of the region, and we see that β-diversity at the scale of the country is moderately high and generally influenced by processes that cause species turnover (βsim). Our selected environmental drivers comprise mostly thermal properties of seawater with a strong spatial dependence across a broad scale (82.4%) and a smaller amount of unknown non-thermal broad-scale spatial influences (13.7%). Only in the species-poor west coast region is there evidence of neutral assembly processes (e.g., dispersal limitation, and stochastic processes as one might expect within kelp forests), as shown by the higher rates of nestedness-resultant β-diversity (βsne) emerging there. We recognize that some environmental drivers such as nutrients can and do co-vary with temperature, but they were not used in our analysis because they are not as readily measured as temperature, they are not always directly produced by coupled global ocean-atmosphere models and thus have less value as predictive (long-term) variables, and because of the issues arising due to collinearity.
The most comprehensive assessment of global marine biogeography (Spalding et al., 2007) places much of South Africa's coast within two realms: Temperate Southern Africa (from southern Angola to around 56) and the western-most edge of the Western Indo-Pacific Realm (across the Western Indian Ocean to Sumatra). The transition into the Western Indian Ocean realm is just visible in section 57 where our data show a steady rise to tropical water temperatures with an annual minimum that exceeds 20°C, between 28.5°S and 29°S. Our study region is too restricted to capture the rates of β-diversity change at the transition of realms. Here we would expect higher rates of turnover than we report for the smaller-scale “Provinces” and transition zones that make up the Temperate Southern African Realm (see below). We would also expect this at higher taxonomic levels, consistent with the definition for a realm (Spalding et al., 2007): this is true for the temperate Southern African seaweed floras, which have high species endemism but very low generic endemism (except for the Fucales, Ochrophyta). We hypothesize that this is because the cool temperate region may be geologically recent (4.6–3.2 Ma, Marlow et al., 2000).
Nested within the realms are marine provinces, areas that according to Spalding et al. (2007) are large, with distinct biotas, some level of endemism (mainly at species level), and distinctive abiotic environments. In this context, the Benguela and Agulhas Currents prescribe the broadest scale hydrographic features whose imprint can be seen, at that scale, on the seaweed flora: they maintain the Benguela Marine Province (BMP) and Agulhas Marine Province (AMP, as per Spalding et al., 2007) of South Africa, separated by an area of transition. The southward flowing Agulhas Current has an overriding effect on the east coast, extending as far as the eastern portion of the Western Cape Province (58–22). Along the east coast it is responsible for setting up a region that transitions from sub-tropical in the north to warm temperate near Cwebe (i.e., the region 58–44), which is termed the East Coast Transition Zone (ECTZ). From 44–22 the Agulhas Current continues to maintain an influence on the coast, but its direct effect is subdued because of the widening of the continental shelf southwards of Cwebe, and the warm-temperate biogeographical region, the AMP, is consequentially formed. The northward flowing Benguela Current from which upwelling is maintained by prevailing south-easterly trade winds influences the remainder of the Western Cape Province (west of 22), particularly from the western side of the Cape Peninsula (16/17) northwards to about 16°S. The influence of the Benguela Current here defines a cool-temperate regime, the BMP, with the range of monthly mean temperatures at most sections intermediate between cold-temperate and warm-temperate, according to the definitions of Lüning et al. (1990).
The oceanographic features of the Benguela and Agulhas Currents set up gradients in seawater temperature that run predominantly from east to west (i.e., August mean temperatures decreasing in this direction), but also from north to south (i.e., with the southern coastline becoming more variable in its thermal regime). Of particular ecological importance is the gradient in August mean temperature. It clearly structures the thermal distance decay curve, which shows a steadily increasing dissimilarity in thermal distance along the shore from 1. The gradient may be smooth or punctuated at irregular intervals by peaks in dissimilarity. Specifically, the mean temperature for Austral winter (August) increases more-or-less smoothly from the west. Within the region, species richness along the coastline reflects the general global trend—for most taxa but ironically not for seaweeds (Bolton, 1994; Santelices et al., 2009)—of diminishing diversity with decreasing temperature, which at that scale is seen as a latitudinal gradient: the cool temperate area of the west coast has a much lower α-diversity than the warm temperate and sub-tropical south and east coasts, with the latter two regions having similar values. This same pattern exists for fish and invertebrates in this region, albeit in data which include coastal waters further from the shore (Griffiths et al., 2010). In the previous presentation of the seaweed dataset (Bolton and Stegenga, 2002), the ECTZ had far fewer species per section than the AMP, but this change in the dataset results from significant research in the period 2000–2005 on the seaweeds of the ECTZ (De Clerck et al., 2005). In contrast, bringing the annual mean temperature or thermal range during Austral summer into the calculation of the thermal distance introduces several large punctuations in dissimilarities between some sections—most notably near Cape Point into False Bay (17), and near 30 and 48. The region of the coastline where low biodiversity community composition changes to highly biodiverse communities is reflected in the β-diversity, which peaks around the Cape Peninsula (16/17). Correlations between β-diversity (resulting from a change in α-diversity) and areas of transition have also been noted in several studies on terrestrial biota (Melo et al., 2009; Tonial et al., 2012), but such investigations are less established for marine biota (Schils and Wilson, 2006; Anderson et al., 2013).
β-diversity is generally lower within individual bioregions than at the scale of the country, because the coastlines are shorter and large differences in diversity cannot approach that which can occur over the entire coastline. Even within bioregions, β-diversity is still generally dominated by species turnover rather than by nestedness-resultant processes. The exception is the BMP where nestedness-resultant (βsne) assembly processes contribute to around 58% of the β-diversity. The generally strong influence of turnover-based processes elsewhere highlights the role of environmental drivers, which here are spatially structured as one would expect of the steep environmental gradients, and hence of niche influences (i.e., environmental filtering) in determining species assembly (Fitzpatrick et al., 2013).
It seems plausible that factors linked with the establishment of oceanographic regimes (historical factors sensu Baselga et al., 2012) that influenced dispersal and speciation along the coast resulted in the current-day gradients in β-diversity within bioregions. This hypothesis was put forward to explain the low endemism and low species richness of Ochrophyta in the upwelling dominated west coast of southern Africa (Bolton, 1986). But is the underlying mechanism the selection of cold tolerant species from amongst a more diverse pool in the wider region (Bolton, 1986), or is it due to a dispersal barrier that prevented elements of the warmer south coast seaweed flora from occupying the west coast? Evidence from seaweed phylogeography postulates different evolutionary origins for seaweeds of the BMP and AMP (Hommersand, 1986; Hommersand and Fredericq, 2003). The mechanism could be the intensification of upwelling during the mid-Pliocene 4.6–3.2 Ma (Marlow et al., 2000). Turnover is almost entirely absent in favor of nestedness-resultant processes in the BMP, but the AMP is defined entirely by a moderate rate of turnover with respect to both geographical and thermal distance. The dominance of the nestedness-resultant component of β-diversity within the BMP suggests that the complex, non-linear influence of the annual mean temperature selectively influences species richness from section to section in a pattern that is not spatially coherent, but this climatic variable (i.e., fraction [E], which has some spatial structure) explains only 39.1% of the βsne patterning in the region. Unaccounted-for non-climatic variables and perhaps non-environmental variables may explain some of the remaining variation. This points to habitat and/or environmental heterogeneity and neutral processes such as dispersal limitation, stochastic influences or habitat heterogeneity as the main structuring agents of communities within the region. Indeed, habitat heterogeneity generated through stochastic influences such as patchiness due to storms denuding portions of kelp beds (which is the dominant habitat type along the west coast) is well-known in kelp forests (Smale et al., 2011). Furthermore, the prevalence of two kinds of assembly processes to either side of the Cape Peninsula (17), coinciding with the point where a region of low α-diversity (BMP) transitions into a region of higher α-diversity (the Benguela-Agulhas Transition Zone, B-ATZ), suggests that the answer to our question is probably that some historical event was responsible for that disjunction, and that this barrier is still maintained today due to limited mixing between the Agulhas and Benguela Currents (although some of the sections along the west coast between 16/17 and 10 certainly have a reasonable probability of sharing species with sections at the western-most side of the AMP, e.g., 23, 24; see more on the B-ATZ, below).
Rapid rates of species turnover with respect to geographic and thermal distance along the east coast, and the steep decrease in species richness around the 16/17–22 transition region, suggest that relatively more endemic species should occur within the BMP and AMP. Seawater becomes warmer along a steep gradient along the east coast northward toward 54 where-after there is a further steepening of the gradient toward Mozambique (north of 58). This steepening of the already strong temperature gradient in the northern portion of the ECTZ supports the conclusion of Bolton and Anderson (2004) that it represents the transition from a tropical Indian Ocean seaweed flora in Mozambique to a temperate flora in the south. A similar pattern exists for the rocky intertidal biota in the region (Sink et al., 2005). The temperature of the coldest month of the year sufficiently accounts for the environmental gradient in all bioregions, except for in the BMP. This strong coupling between the thermal gradient and βsim suggests that a niche difference mechanism is the primary species compositional assembly process (Nekola et al., 1999). This implies species being sorted based on their physiological tolerances along a gradient, and the particularly high rate of species turnover along the east coast reflects the species' narrow thermal ranges. Species here have narrow local distributions, and there is only a 2% seaweed endemism in the region (Bolton and Anderson, 2004).
The overlap region between 22 and the southern tip of the Cape Peninsula (16/17), which we call the B-ATZ, is an area where aspects of both currents may be periodically seen (Largier et al., 1992), and it is not surprising that biogeographically this region includes biotic elements from the marine provinces on both sides of it. This is clearly visible in our analyses, which show the highest rates of species turnover (βsim) with respect to geographical distance, the mean temperature for August and the range in temperature in February. This is a form of regional biogeographical structuring called the temporal/spatial constraint model (Nekola et al., 1999). Similar conclusions with regards to the B-ATZ being an transitional area have been reached by Stephenson (1948), Bolton (1986), Bolton and Anderson (2004), and Mead et al. (2013). For marine macroalgae at least, this area cannot be recognized as a marine province as it does not meet the criteria for it to be classified as such (Spalding et al., 2007), i.e., it lacks cohesion, levels of endemism, and there is an absence of distinct abiotic features (with exception of such conditions periodically arising in False Bay due to prevailing mesoscale oceanographic conditions).
It was recently suggested that studies which associate marine macroalgal distributional patterns with broad-scale temperature gradients neglect to consider alternative or additional explanations, such as that offered by connectivity due to ocean currents (Wernberg et al., 2013). It is possible that seaweeds are similarly influenced by ocean currents around South Africa, and it is temping to suggest that the direction of turnover along the east coast is from north to south, and along the south coast from east to west, as this would coincide with the direction of the Agulhas Current. From a purely thermal point of view, the imprints of those ocean currents are certainly present at the coast in the zone that our seaweed data represent, but the physical processes that operate in that coastal zone are considerably more complex and decoupled from mesoscale influences (Schlegel et al., 2017). Personal observations by one of us (AS) recount dislodged individuals of the kelp, Ecklonia maxima, having washed ashore some 400 km east of the eastern-most distributional limit for the species in South Africa—this is against the direction of the current, which in this region is situated considerably far from the shore along the edge of the Agulhas Bank. Our data on the rate of species turnover as a function of the thermal difference between coastal sections (βsim vs. dE) within bioregions show an almost precise relationship between thermal difference and change in species composition (even though the relationship between dE and geographical distance differs between the east and south coasts), and provide very strong support for temperature, rather than connectivity due to currents, as an overriding driver of seaweed community composition over large spatial scales. Within the context of Beijerinck's “Law” that “everything is everywhere but the environment selects” (De Wit and Bouvier, 2006), the only role that ocean currents can have in setting up biogeographical patterns is to cause everything not to be everywhere. At these scales, with or without the influence of ocean currents, temperature still selects. In other words, at broad spatial scales within the region encompassed by the ECTZ and the AMP, the seaweed flora is structured predominantly by niche-assembly processes driven by a thermal gradient.
However, we should point out that ca. 13.7% of the seaweed flora's turnover is also explained by some other broad scale spatial influences (MEM2, MEM3 and MEM5), which represent exogenous (unmeasured) environmental variables and autocorrelation. The importance of these relative to each other is unknown. Future work may find that some of these environmental influences may indeed be that of ocean currents, which operate at those broad scales, but our current interpretation does not support this notion. Rather, we suggest that the unknown influence is more likely to be neutral influences (Hubbell, 2001) manifesting at a reasonably broad scale. These may include stochastic events (e.g., Smale et al., 2011) or biotic processes caused by differences in species demography, dispersal or the influence of grazers, etc. Such processes are usually finer scaled, but, having said that, the entire west coast region that makes up the BMP is dominated by nestedness-resultant β-diversity, which is also symptomatic of neutral assembly processes.
Our approach throughout this paper hinges upon the biogeographical provinces and overlap regions defined using both seaweed and thermal data. Together, thermal and species gradients/patterns provide the most parsimonious explanation for the processes assembling the seaweed flora within the study region. The explanation is that thermal gradients set up β-diversity based on species turnover as a result of environmental filtering within the southern and eastern coastal sections, while along the west coast a very different kind of assembly process, nestedness-resultant β-diversity, emerges. We suggest that this disjunction results from a historical event that is currently still limiting connectivity between the Indian and Atlantic Oceans around the Cape Peninsula.
We also show how thermal metrics other than the annual mean temperatures can aid in the delineation of macroecological patterns. Some species are limited by cold and some by warm temperatures, and the constraining factor may differ at a species' southern and northern limits (e.g., van den Hoek, 1982a; Breeman, 1988). If these are important determinants, as we indeed show there are, other measures of climate should be used in our repertoire of explanatory variables. For example, the mean temperature for the coldest month, standard deviation, variation in annual temperature extremes, or temperature ranges (e.g., as used by Qian and Ricklefs, 2007; Leaper et al., 2011) may readily be extracted from time series of daily temperatures. These approaches, as we apply them here, may greatly improve macroecological and biogeographic studies and aid our ability to analyse broad-scale communities patterns and the processes that assemble them.
AS conceptualized the scope of the research reported in this paper, undertook all the numerical and statistical analyses, made the first round of interpretation, and did the bulk of the writing. JB and RA collected all the samples over the last 30 years, compiled the database of species distribution records for the region included in this analysis, contributed in equal part to the conceptualization of the research and the interpretation of the findings, and provided significant editorial input into penultimate and final drafts of the document.
The research was partly funded by the South African National Research Foundation (http://www.nrf.ac.za) programme “Thermal characteristics of the South African nearshore: implications for biodiversity” (CPRR14072378735). Aside from providing funding, the funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Robert Schlegel is thanked for his assistance in the preparation of some of Figure 1.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2017.00404/full#supplementary-material
Anderson, M. J., Tolimieri, N., and Millar, R. B. (2013). Beta diversity of demersal fish assemblages in the North-Eastern Pacific: interactions of latitude and depth. PLoS ONE 8:e57918. doi: 10.1371/journal.pone.0057918
Barton, P. S., Cunningham, S. A., Manning, A. D., Gibb, H., Lindenmayer, D. B., and Didham, R. K. (2013). The spatial scaling of beta diversity. Glob. Ecol. Biogeogr. 22, 639–647. doi: 10.1111/geb.12031
Baselga, A., Gómez-Rodríguez, C., and Lobo, J. M. (2012). Historical legacies in world amphibian diversity revealed by the turnover and nestedness components of beta diversity. PLoS ONE 7:e32341. doi: 10.1371/journal.pone.0032341
Baselga, A., Orme, D., Villeger, S., De Bortoli, J., and Leprieur, F. (2017). betapart: Partitioning Beta Diversity into Turnover and Nestedness Components. R package version 1.4-1. Available online at: https://CRAN.R-project.org/package=betapart
Bolton, J. J. (1986). Marine phytogeography of the Benguela upwelling region on the west coast of southern Africa: a temperature dependent approach. Bot. Mar. 29, 251–256. doi: 10.1515/botm.1922.214.171.124
Bolton, J. J., and Anderson, R. J. (2004). “Marine vegetation,” in Vegetation of Southern Africa, eds R. M. Cowling, D. M. Richardson, and S. M. Pierce, (Cambridge, UK: Cambridge University Press), 348–370.
Bolton, J. J., Leliaert, F., De Clerck, O., Anderson, R. J., Stegenga, H., Engledow, H. E., et al. (2004). Where is the western limit of the tropical Indian Ocean seaweed flora? An analysis of intertidal seaweed biogeography on the east coast of South Africa. Mar. Biol. 144, 51–59. doi: 10.1007/s00227-003-1182-9
Breeman, A. M. (1988). Relative importance of temperature and other factors in determining geographic boundaries of seaweeds: experimental and phenological evidence. Helgoländer Meeresuntersuchungen 42, 199–241. doi: 10.1007/BF02366043
Broitman, B. R., Blanchette, C. A., Menge, B. A., Lubchenco, J., Krenz, C., Foley, M., et al. (2008). Spatial and temporal patterns of invertebrate recruitment along the west coast of the United States. Ecol. Monogr. 78, 403–421. doi: 10.1890/06-1805.1
Byrne, M., Ho, M., Selvakumaraswamy, P., Nguyen, H. D., Dworjanyn, S. A., and Davis, A. R. (2009). Temperature, but not pH, compromises sea urchin fertilization and early development under near-future climate change scenarios. Proc. R. Soc. B Biol. Sci. 276, 1883–1888. doi: 10.1098/rspb.2008.1935
Davidar, P., Rajagopal, B., Mohandass, D., Puyravaud, J.-P., Condit, R., Wright, S., et al. (2007). The effect of climatic gradients, topographic variation and species traits on the beta diversity of rain forest trees. Glob. Ecol. Biogeogr. 16, 510–518. doi: 10.1111/j.1466-8238.2007.00307.x
De Wit, R., and Bouvier, T. (2006). ‘Everything is everywhere, but, the environment selects’; what did Baas Becking and Beijerinck really say? Environ. Microbiol. 8, 755–758. doi: 10.1111/j.1462-2920.2006.01017.x
Dormann, C. F., Elith, J., Bacher, S., Buchmann, C., Carl, G., Carré, G., et al. (2013). Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography 36, 027–046. doi: 10.1111/j.1600-0587.2012.07348.x
Dray, S., Legendre, P., and Peres-Neto, P. R. (2006). Spatial modelling: a comprehensive framework for principal coordinate analysis of neighbour matrices (PCNM). Ecol. Model. 196, 483–493. doi: 10.1016/j.ecolmodel.2006.02.015
Dray, S., Pélissier, R., Couteron, P., Fortin, M. J., Legendre, P., Peres-Neto, P. R., et al. (2012). Community ecology in the age of multivariate multiscale spatial analysis. Ecol. Monogr. 82, 257–275. doi: 10.1890/11-1183.1
Fitzpatrick, M. C., Sanders, N. J., Normand, S., Svenning, J.-C., Ferrier, S., Gove, A. D., et al. (2013). Environmental and historical imprints on beta diversity: insights from variation in rates of species turnover along gradients. Proc. R. Soc. B Biol. Sci. 280:20131201. doi: 10.1098/rspb.2013.1201
Griffiths, C. L., Robinson, T. B., Lange, L., and Mead, A. (2010). Marine biodiversity in South Africa: an evaluation of current states of knowledge. PLoS ONE 5:e12008. doi: 10.1371/journal.pone.0012008
Hommersand, M. H., and Fredericq, S. (2003). “Biogeography of the marine red algae of the South African West Coast: a molecular approach,” in Proceedings of the 7th International Seaweed Symposium (Oxford: Oxford University Press), 325–340.
IPCC. (2013). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press UK; New York, NY, USA.
Largier, J. L., Chapman, P., Peterson, W. T., and Swart, V. P. (1992). The western Agulhas Bank: circulation, stratification and ecology. South Afr. J. Mar. Sci. 12, 319–339. doi: 10.2989/02577619209504709
Leaper, R., Hill, N. A., Edgar, G. J., Ellis, N., Lawrence, E., Pitcher, C. R., et al. (2011). Predictions of beta diversity for reef macroalgae across southeastern Australia. Ecosphere 2, 1–18. doi: 10.1890/ES11-00089.1
Lombard, A. T., Strauss, T., Harris, J., Sink, K., Attwood, C., and Hutchings, L. (2004). South African National Spatial Biodiversity Assessment 2004. Technical report, South African National Biodiversity Institute.
Marlow, J. R., Lange, C. B., Wefer, G., and Rosell-Mele, A. (2000). Upwelling intensification as part of the Pliocene-Pleistocene climate transition. Science 290, 2288–2291. doi: 10.1126/science.290.5500.2288
Mead, A., Griffiths, C. L., Branch, G. M., McQuaid, C. D., Blamey, L. K., Bolton, J. J., et al. (2013). Human-mediated drivers of change – impacts on coastal ecosystems and marine biota of South Africa. Afr. J. Mar. Sci. 35, 403–425. doi: 10.2989/1814232X.2013.830147
Melo, A. S., Rangel, T. F. L. V. B., and Diniz-Filho, J. A. F. (2009). Environmental drivers of beta-diversity patterns in New-World birds and mammals. Ecography 32, 226–236. doi: 10.1111/j.1600-0587.2008.05502.x
Nekola, J. C., White, P. S., Carolina, N., Nekola, C., and Curriculum, P. S. W. (1999). The distance decay of similarity in biogeography and ecology. J. Biogeogr. 26, 867–878. doi: 10.1046/j.1365-2699.1999.00305.x
Oksanen, J., Guillaume Blanchet, F., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2017). vegan: Community Ecology Package. R package version 2.4-5. Available online at: https://CRAN.R-project.org/package=vegan
Peres-Neto, P. R., and Legendre, P. (2010). Estimating and controlling for spatial structure in the study of ecological communities. Glob. Ecol. Biogeogr. 19, 174–184. doi: 10.1111/j.1466-8238.2009.00506.x
Peres-Neto, P. R., Legendre, P., Dray, S., and Borcard, D. (2006). Variation partitioning of species data matrices: estimation and comparison of fractions. Ecology 87, 2614–2625. doi: 10.1890/0012-9658(2006)87[2614:VPOSDM]2.0.CO;2
Potts, W. M., Henriques, R., Santos, C. V., Munnik, K., Ansorge, I., Dufois, F., et al. (2014). Ocean warming, a rapid distributional shift, and the hybridization of a coastal fish species. Glob. Change Biol. 20, 2765–2777. doi: 10.1111/gcb.12612
Schlegel, R. W., Oliver, E. C. J., Wernberg, T., and Smit, A. J. (2017). Nearshore and offshore co-occurrence of marine heatwaves and cold-spells. Prog. Oceanogr. 151, 189–205. doi: 10.1016/j.pocean.2017.01.004
Sink, K. J., Branch, G. M., and Harris, J. M. (2005). Biogeographic patterns in rocky intertidal communities in KwaZulu-Natal, South Africa. Afr. J. Mar. Sci. 27, 81–96. doi: 10.2989/18142320509504070
Smale, D., Kendrick, G., and Wernberg, T. (2011). Subtidal macroalgal richness, diversity and turnover, at multiple spatial scales, along the southwestern Australian coastline. Estuar. Coast. Shelf Sci. 91, 224–231. doi: 10.1016/j.ecss.2010.10.022
Smale, D. A., Kendrick, G. A., and Wernberg, T. (2010). Assemblage turnover and taxonomic sufficiency of subtidal macroalgae at multiple spatial scales. J. Exp. Mar. Biol. Ecol. 384, 76–86. doi: 10.1016/j.jembe.2009.11.013
Smit, A. J., Roberts, M., Anderson, R. J., Dufois, F., Dudley, S. F. J., Bornman, T. G., et al. (2013). A coastal seawater temperature dataset for biogeographical studies: large biases between in situ and remotely-sensed data sets around the coast of South Africa. PLoS ONE 8:e81944. doi: 10.1371/journal.pone.0081944
Spalding, M. D., Agostini, V. N., Rice, J., and Grant, S. M. (2012). Pelagic provinces of the world: a biogeographic classification of the world's surface pelagic waters. Ocean Coas. Manag. 60, 19–30. doi: 10.1016/j.ocecoaman.2011.12.016
Spalding, M. D., Fox, H. E., Allen, G. R., Davidson, N., Ferdaña, Z. A., Finlayson, M. A. X., et al. (2007). Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. BioScience 57, 573–583. doi: 10.1641/B570707
Straub, S. C., Thomsen, M. S., and Wernberg, T. (2016). “The dynamic biogeography of the anthropocene: the speed of recent range shifts in seaweeds,” in Seaweed Phylogeography, eds Z.-M. Hu and C. Fraser (Dordrecht: Springer), 63–93.
Stuart-Smith, R. D., Edgar, G. J., and Bates, A. E. (2017). Thermal limits to the geographic distributions of shallow-water marine species. Nat. Ecol. Evol. 1, 1846–1852. doi: 10.1038/s41559-017-0353-x
Tittensor, D. P., Mora, C., Jetz, W., Lotze, H. K., Ricard, D., Berghe, E. V., et al. (2010). Global patterns and predictors of marine biodiversity across taxa. Nature 466, 1098–1101. doi: 10.1038/nature09329
Tonial, M., Silva, H., Tonial, I., Costa, M., Silva Júnior, N. J., and Diniz-Filho, J. (2012). Geographical patterns and partition of turnover and richness components of beta-diversity in faunas from Tocantins river valley. Brazil. J. Biol. 72, 497–504. doi: 10.1590/S1519-69842012000300012
Turpie, J. K., Beckley, L. E., and Katua, S. M. (2000). Biogeography and the selection of priority areas for conservation of south african coastal fishes. Biol. Conserv. 92, 59–72. doi: 10.1016/S0006-3207(99)00063-4
Tuya, F., Cacabelos, E., Duarte, P., Jacinto, D., Castro, J., Silva, T., et al. (2012). Patterns of landscape and assemblage structure along a latitudinal gradient in ocean climate. Mar. Ecol. Prog. Ser. 466, 9–19. doi: 10.3354/meps09941
Tyberghein, L., Verbruggen, H., Pauly, K., Troupin, C., Mineur, F., and De Clerck, O. (2012). Bio-ORACLE: a global environmental dataset for marine species distribution modelling. Glob. Ecol. Biogeogr. 21, 272–281. doi: 10.1111/j.1466-8238.2011.00656.x
van den Hoek, C. (1982a). Phytogeographic distribution groups of benthic marine algae in the North Atlantic Ocean. A review of experimental evidence from life history studies. Helgoländer Meeresuntersuchungen 35, 153–214.
Verbruggen, H., Tyberghein, L., Pauly, K., Vlaeminck, C., Nieuwenhuyze, K. V., Kooistra, W. H. C. F., et al. (2009). Macroecology meets macroevolution: evolutionary niche dynamics in the seaweed Halimeda. Glob. Ecol. Biogeogr. 18, 393–405. doi: 10.1111/j.1466-8238.2009.00463.x
Waters, J. M., Wernberg, T., Connell, S. D., Thomsen, M. S., Zuccarello, G. C., Kraft, G. T., et al. (2010). Australia's marine biogeography revisited: back to the future? Aust. Ecol. 35, 988–992. doi: 10.1111/j.1442-9993.2010.02114.x
Wernberg, T., Thomsen, M. S., Connell, S. D., Russell, B. D., Waters, J. M., Zuccarello, G. C., et al. (2013). The footprint of continental-scale ocean currents on the biogeography of seaweeds. PLoS ONE 8:e80168. doi: 10.1371/journal.pone.0080168
Wieters, E. A., Broitman, B. R., and Branch, G. M. (2009). Benthic community structure and spatiotemporal thermal regimes in two upwelling ecosystems: comparisons between South Africa and Chile. Limnol. Oceanogr. 54, 1060–1072. doi: 10.4319/lo.2009.54.4.1060
Keywords: beta-diversity, species assembly, seaweed, turnover, nestedness-resultant, Benguela Current, Agulhas Current
Citation: Smit AJ, Bolton JJ and Anderson RJ (2017) Seaweeds in Two Oceans: Beta-Diversity. Front. Mar. Sci. 4:404. doi: 10.3389/fmars.2017.00404
Received: 13 July 2017; Accepted: 28 November 2017;
Published: 12 December 2017.
Edited by:Thomas Wilke, Justus Liebig Universität Gießen, Germany
Reviewed by:Jian-Wen Qiu, Hong Kong Baptist University, Hong Kong
Fernando Tuya, University of Las Palmas de Gran Canaria, Spain
Copyright © 2017 Smit, Bolton and Anderson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Albertus J. Smit, email@example.com