Next Article in Journal / Special Issue
Plant Part Age and Size Affect Sessile Macrobenthic Assemblages Associated with a Foliose Red Algae Phycodrys rubens in the White Sea
Previous Article in Journal
Do Local Environmental Factors and Lunar Cycle Influence Timing and Synchrony of Oviposition of a Turtle with Strict Nocturnal Nesting?
Previous Article in Special Issue
Biogeographical Patterns of Endolithic Infestation in an Invasive and an Indigenous Intertidal Marine Ecosystem Engineer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Canopy Cover on the Ecological Function of A Key Autogenic Ecosystem Engineer

1
Department of Marine Invertebrates, Museum Victoria, Carlton, VIC 3053, Australia
2
School of BioSciences, University of Melbourne, Parkville, VIC 3010, Australia
3
School of Life and Environmental Sciences, Centre for Integrative Ecology, Deakin University, Warrnambool Campus, PO Box 423, Warrnambool, VIC 3280, Australia
*
Author to whom correspondence should be addressed.
Diversity 2019, 11(5), 79; https://doi.org/10.3390/d11050079
Submission received: 8 April 2019 / Revised: 12 May 2019 / Accepted: 15 May 2019 / Published: 17 May 2019
(This article belongs to the Special Issue Diversity of Ecosystem Engineers in the World Coasts and Oceans)

Abstract

:
Intertidal fucoid algae can function as ecosystem engineers across temperate marine regions. In this investigation, we assessed the function of the alga dominating rocky reefs in temperate Australia and New Zealand, Hormosira banksii. Invertebrate and algal species assemblages were examined within areas of full H. banksii canopy, areas where it was naturally patchy or absent (within its potential range on the shore) and areas where the intact canopy was experimentally disturbed. Differences in species assemblages were detected between areas with natural variation in H. banksii cover (full, patchy, negligible), with defined species associated with areas of full cover. Differences were also detected between experimentally manipulated and naturally patchy areas of canopy cover. Species assemblages altered in response to canopy manipulations and did not recover even twelve months after initial sampling. Both light intensity and temperature were buffered by full canopies compared to patchy canopies and exposed rock. This study allows us to predict the consequences to the intertidal community due to the loss of canopy cover, which may result from a range of disturbances such as trampling, storm damage, sand burial and prolonged exposure to extreme temperature, and further allow for improved management of this key autogenic ecosystem engineer.

1. Introduction

Canopy-forming seaweeds can function as autogenic ecosystem engineers on rocky seashores, providing a suitable habitat for species that would otherwise be excluded or only occur in low abundance [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]. Many studies referring to species as autogenic ecosystem engineers [16,17] (or similar theoretical terms such as dominant species [18], facilitators, habitat-forming species [9] and biogenic habitat) have investigated the associated loss of species when this one species is lost through a major disturbance. What hasn’t been examined in full is whether the fucoid intertidal alga Hormosira banksii functions as an ecosystem engineer under varying densities in the absence of disturbance. Is it dependent on spatial dominance, or can it function as part of a mosaic of habitat density and type? Here we demonstrate through both mensurative and manipulative experiments how H. banksii functions as an autogenic ecosystem engineer, and how its ability to modify conditions, and thus habitat provision, differs under varying canopy cover. We also show that manipulative experiments alone cannot demonstrate that H. banksii is an autogenic ecosystem engineer (defined as changing the environment via its physical structure, [3]) and habitat provider on Victorian rocky shores.
Throughout its range, the fucoid Hormosira banksii is the dominant canopy-forming species on intertidal rocky shores, and unlike studies of fucoids in the northern hemisphere, it has no functional equivalent or ecological successor [2,5,6,9,12,19,20,21,22,23,24,25,26,27,28,29], likely due in part to its specialised ability to cope with harsh environmental conditions associated with Austral summers during tidal emersion [30]. Few studies examining the response of H. banksii understorey communities to disturbance have investigated the response of more than a few species or components (e.g., functional groups) of the system [9,12,17,25,31]. These studies highlight the potential for H. banksii to be categorised as an autogenic ecosystem engineer, though none have yet examined whether H. banksii also shows this function under naturally reduced densities.
Disturbance studies have shown that damage to H. banksii through trampling or experimental removal causes a decline in the abundance of species usually associated with the canopy [2,5,6,9,12]. H. banksii canopy-removal experiments in New Zealand [9,27] found that temperature was lower within the canopy in comparison to areas where it was removed and this was most apparent during the summer months when surface temperatures in removal plots were at least 5 °C higher than in control plots [9]. Lilley and Schiel [9] found that understorey algae declined in abundance, and an influx of ephemeral algae occurred (e.g., Ulva spp., Colpomenia sinuosa). The abundance of invertebrates also changed where some species declined in abundance (e.g., Micrelenchus tesselatus, as Cantharidella tesselata) and some increased (e.g., Notoacmea spp.). Trampling experiments in south-eastern Australia have shown that bare rock habitat and its associated fauna weren’t detrimentally impacted by trampling, yet H. banksii was damaged even at low trampling intensities with a lack of recovery even after 400 days in heavy trampling treatments [2]. Invertebrates common to bare areas were found to colonise the heavily damaged areas of H. banksii indicating the possibility of a shift in the community from trampling [2]. In contrast, King [21] found that variable responses to trampling and suggested that seasonal changes potentially obscure the impacts of trampling disturbance. Although both of these trampling studies found low power was a problem in testing the responses of individual species to trampling disturbance [2,21], they did, however, show that different species assemblages occupy different parts of the shore and that high variability is indicative of complexity beyond simple facilitative relationships.
The ecological function of H. banksii may vary in space and time and given the results of disturbance studies [2,5,6,9,12], the magnitude of its habitat provision could depend on its physical structure or canopy ‘cover’. Previous studies, however, have not examined the species relationships that occur with naturally reduced ‘patchy’ canopies. Here we examine whether the ability of H. banksii to modulate the understorey environment (e.g., light, temperature, desiccation stress) for associated species [3] varies in relation to canopy cover which is naturally variable on Victorian rocky shores. This study investigated the relationship between H. banksii and co-occurring species using mensurative surveys examining natural patterns and experimental manipulations at multiple sites, to determine whether the amount of canopy cover influences the autogenic ecosystem engineering function of H. banksii on Victorian rocky shores.

2. Materials and Methods

2.1. Study Sites

All sites were moderately wave-exposed sites facing Bass Strait [32]. On Victorian shores (south-eastern Australia), areas within the vertical range of H. banksii can be separated into: full beds (continuous canopy cover with little to no understorey visible), patchy areas (canopy cover is interspersed with ‘bare’ rock (devoid of observable biota), sand or alternate biota such as mussels or turfing algae), and areas where it is rare or absent (see Figure A1 for examples). All the study reefs were flat calcarenite reef platforms with a mixture of habitat types including extensive Hormosira banksii beds and sections of platform without canopy, comprising turf-forming algae, mussels and ‘bare’ rock (hereafter referred to as ‘negligible’).
Although substratum temperatures often exceed 30 °C during low tide in the summer months, the study species Hormosira banksii has a high tolerance to desiccation due to its unique fluid-filled beaded morphology and the presence of UV filters (physodes) in the cell walls [30]; thus H. banksii is the only fucoid capable of forming a canopy in the intertidal of mainland temperate Australia. The Victorian coastline has a micro-tidal range (below 2 m, [33]); all sites were sampled when tides were predicted to be below 0.4 m above mean low water (MLW).
Four sites were surveyed including two sites selected for manipulative experiments on canopy cover hereafter referred to as “experiment” sites and two sites selected to compare to the controls at experiment sites hereafter referred to as “background” sites. Experiment Sites: Bridgewater Bay, Mornington Peninsula, Victoria, Australia and 13th Beach 32W (Barwon Heads), Bellarine Peninsula, Victoria, Australia are both backed by steep headlands (Figure 1). Bridgewater Bay is a narrower platform (30 m) than 13th Beach (60 m) though both extend for hundreds of metres along the shoreline. Background Survey Sites: Sorrento Back Beach, Mornington Peninsula, Victoria, Australia (see control area in Reference [29]) and Point Roadknight, Anglesea, Victoria, Australia are both very large platforms hundreds of meters long, and approximately 100-m wide (Figure 1). Sorrento Back Beach runs parallel along the coastline, whereas Point Roadknight forms a southerly point.

2.2. Investigation of Species Associations with Natural Variations in Canopy Cover of Hormosira banksii

To determine the species assemblages associated with the H. banksii canopy, sampling was undertaken at the two background sites and two experimental sites prior to experimental manipulation (see experimental investigation) using haphazardly placed quadrats (0.09 m2) within the area of H. banksii distribution on the platform. Sampling occurred within three distinguished habitat types: negligible (<5%), patchy (~40%–80% configured as multiple large clumps separated by ‘bare’ rock without H. banksii), and full (>90%) areas of H. banksii canopy (see Figure A1). These discrete habitats were highly variable in their connectivity to one another. In each habitat, replicate quadrats (0.09 m2) were quantitatively sampled for percentage cover (using 49-point intercepts) of all components (e.g., sand, rock), sessile aggregating species (e.g., Capreolia implexa, Xenostrobus pulex both ≤ 1 cm). The abundance of all mobile and sedentary macroinvertebrate species was individually counted both on fronds of Hormosira and on the substrate beneath the canopy (background sites n = 10, experiment sites n = 5; see full species list Table A1). Percentage cover records included quantifying the primary cover layer which includes all dominant components and biota (referred to throughout as ‘primary’) and the layer below the canopy (referred to throughout as ‘understorey’). In negligible and removed plots where there is no canopy the ‘understorey’ is equivalent to the primary cover and where described as understorey, it is presented to serve comparison. Organisms were identified to the lowest taxonomic level possible. The depth of sand over rock was measured within quadrats at three random points to the nearest 2.5 mm. The lengths of the longest H. banksii fronds were recorded within quadrats from 3 randomly selected fronds to the nearest 2.5 mm. Sampling was undertaken during two visits per site to account for temporal variability in species abundances (Point Roadknight: January and December 2007 (summer); Sorrento: May and December 2007 (autumn and summer respectively)).

2.3. Experimental Investigation of the Effect of Variations in Canopy Cover of Hormosira banksii on the Physical Environment and Associated Species Assemblages

To examine whether algal and macroinvertebrate communities are influenced by the presence of full, patchy or absent canopies; an experiment was carried out at two intertidal rock platforms (13th Beach and Bridgewater Bay). Twenty-five fixed 0.25 m2 plots were established haphazardly across a reef area approximately 10,000 m2 within three categories of canopy cover and marked with labelled cattle tags fixed with stainless steel screws. Five plots were in areas where H. banskii was naturally absent (0%–5%), 5 in areas where H. banksii had naturally a patchy cover (40%–80%) and 15 were established in areas with full canopies (>95%, full control plots). Ten of these full plots (five replicates each) were randomly assigned to the experimental manipulations: thinned (40%–80%) and removed (0%–5%); the five remaining plots remained as full controls. Hormosira banksii cover was reduced by trimming the algal fronds in the experimental plots just above the holdfast using garden secateurs (to prevent long term disturbance by allowing regrowth from holdfasts, Underwood [24]). Any fronds that could overlay the plots from the outside edge were also trimmed. Furthermore, any regrowth above 5 mm was retrimmed as necessary creating a press disturbance (recovery was not a question of interest). Sampling of plots was undertaken as for the background sites using 0.09 m2 gridded quadrats placed in the centre of each plot. Taxa were identified to species in situ where possible (see Table A1 for full species list). Depth of sand to rock was measured within quadrats at three random points to the nearest 2.5 mm. The longest frond of each of three randomly chosen H. banksii plants was measured in each quadrat. To examine whether H. banksii modifies the physical environment to the understorey; temperature was recorded on the substratum of each plot using a waterproof thermometer to one decimal place, and relative light intensity on the substratum was recorded using a custom-built meter (JR Pocklington) including a light sensitive diode connected to a standard electrical meter, which measured the degree of shading. Temperature and light readings were taken beneath the canopy within plots (as they were being sampled for species assemblage). All measurements and surveys for all plots were conducted during on low tide, at each sampling time. Sampling intervals roughly followed a logarithmic time series; at 13th Beach: 1, 5, 15, and 379 days after manipulation (23/2/2007, 27/2/2007, 9/3/2007, 13/3/2008; time n = 4); at Bridgewater Bay: 1, 15, 19, and 371 days after manipulation (23/3/2007, 27/3/2007, 10/4/2007, 28/3/2008; time n = 4). The experiment start dates and sampling intervals differed slightly between sites on two occasions due to logistical difficulties. Long intervals across the winter period were due to low tide predictions being too high to expose the H. banksii zone, or only occurring during the night when species assemblages may differ from day records and accurate sampling couldn’t be guaranteed.

2.4. Data Analysis

Comparisons of species assemblages associated with both different natural levels of Hormosira cover and experimental manipulations of Hormosira canopy cover were assessed with multivariate analyses using PRIMER 6 version 6.1.1 and PERMANOVA + version 1.0.5 (PRIMER-E Ltd, 2013, Plymouth, UK). Differences in primary cover, understorey cover and mobile invertebrate assemblages between Hormosira-cover treatments were visualised separately for natural and manipulated covers using Non-metric Multidimensional Scaling (NMDS) [34] based on Bray-Curtis dissimilarities. Ordinations were accepted if stress values were < 0.20 and ideally < 0.10 (as recommended by Reference [35]). Shepherds Plots were examined to confirm the goodness of fit. Two-factor permutational multivariate analyses of variance (PERMANOVAs) were performed on fourth-root transformed data with 999 permutations for both natural (baseline) and manipulated covers, examining primary cover, understorey cover and mobile invertebrate assemblages separately. For comparisons of assemblages associated with different natural levels of Hormosira cover prior to manipulation, treatments (three levels: full, patchy, negligible; fixed) and sites (four levels: 13th Beach, Bridgewater Bay, Point Roadknight, Sorrento; random) were compared by PERMANOVA, but because logistical constraints prevented starting experiments and sampling all sites at the same times, any temporal differences in assemblages is incorporated into the unexplained residual variation in analyses. To examine the influence of experimental manipulation of Hormosira canopy cover on associated assemblages, treatments (five levels: full, patchy, negligible, thinned, removed; fixed) and sites (two levels: 13th Beach, Bridgewater Bay; random) were compared by PERMANOVA at each of days 1, 5, 15/19, and 371/379 since experimental manipulation. Tests for homogeneity of dispersions within treatment groups were performed using PERMDISP [36] with distance to centroids. When the treatment × site interaction was significant another PERMANOVA was run on the interaction term comparing treatments individually by site. In these cases, the P-value based on Monte Carlo random draws, P(MC) was selected as it is more robust with the reduced number of possible permutations in pairwise tests [36]. Similarity Percentages (SIMPER) was used to assess the contributions (> 10%) of individual taxa/cover-groups to the separation of treatments. Repeated measures analyses of variance (rmANOVA) were then used to determine the differences in abundances of individual taxa/cover-groups identified by SIMPER amongst treatments and between sites over time (at days 1, 5, 15/19, and 371/379 since experimental manipulation). The physical factors of temperature and light intensity; and the covariates sand depth and Hormosira frond length were also analysed using rmANOVA amongst treatments and between sites over time. All univariate analyses were conducted using SYSTAT version 10 and α = 0.05 for all statistical analyses. Assumptions of rmANOVA were assessed by examination of probability and residual plots; the abundance of individual mobile invertebrates was square root transformed, otherwise no data transformations were necessary.

3. Results

3.1. Investigation of Species Associations with Natural Variations in Canopy Cover of Hormosira banksii

The differences in assemblages associated with natural variation in canopy cover of Hormosira were not consistent amongst sites for the primary cover, understorey cover or mobile invertebrates (PERMANOVA treatment × site interactions; P(perm) = 0.001 for all analyses). However, for all sites, there were significant differences in the primary cover assemblages between full, patchy and negligible natural Hormosira-cover categories as expected (PERMANOVA on category × site interaction terms; P(perm) < 0.05 for all pairwise comparisons for each site; Figure 2a). The understorey and mobile invertebrate assemblages associated with these natural cover categories (full, patchy and negligible Hormosira) showed some spatial variability (and potentially temporal variability due to differences in sampling times of sites), but both understorey and mobile invertebrate assemblages differed significantly amongst all natural cover categories for each site (PERMANOVA on category × site interaction terms; P(perm) < 0.05 for all pairwise comparisons for each site), except full vs. patchy at 13th Beach (P(perm) = 0.138 and 0.227 for understorey cover and mobile invertebrate assemblages, respectively; Figure 2b,c) and patchy vs. negligible cover for mobile invertebrate assemblages at both Bridgewater Bay and Point Roadknight (P(perm) = 0.165 and 0.075, respectively; Figure 2c). The differences in all assemblages (primary cover, understorey cover and mobile invertebrates) between natural Hormosira-cover categories were, at least in part, due to differences in dispersion amongst Hormosira-cover categories (PERMDISP by categories: P(perm) = 0.001 for all analyses; Figure 2).
As expected, the differences in the primary cover assemblages amongst Hormosira-cover categories identified by SIMPER were primarily driven by differences in the percentage cover of Hormosira, with substrata devoid of visible biota (hereafter ‘bare rock’) and the presence of the mussel Xenostrobus pulex, further contributing to separation of patchy and negligible Hormosira-cover categories, respectively. For the understorey cover, the presence of sand in both full and patchy Hormosira-cover categories contributed 35% and 59%, respectively, to the separation amongst categories (mean ± SE% cover of sand across all sites: 65% ± 5.2 and 24% ± 6.1, respectively), whereas the presence of bare rock, X. pulex, the serpulid worm Galeolaria caespitosa and the rhodophyte Capreolia implexa contributed to 88% of the separation of the negligible Hormosira-cover categories (mean ± SE% cover across all sites: 41% ± 4.9, 31% ± 5.2, 3% ±1.0, 12% ± 3.5, respectively). Mobile invertebrate assemblages were defined by the presence of the gastropods Chlorodiloma adelaidae, Lunella undulata, Cominella lineolata and Dicathais orbita contributing 88% of the separation of the full cover category; the herbivorous gastropods Bembicium nanum, Austrocochlea constricta and a few pulmonates and prosobranch limpet species contributing 86% of the separation of the patchy cover category; and the presence of B. nanum and several pulmonate and prosobranch limpet species contributed to 91% of the separation of the negligible Hormosira-cover category (Figure 3).

3.2. Experimental Investigation of the Effect of Variations in Canopy Cover of Hormosira banksii on the Physical Environment and Associated Species Assemblages

Canopy cover of Hormosira banksii remained within the parameters of the experimental treatments and controls over the course of the experiment (Table 1, Figure 4). However, after experimental manipulation of the Hormosira canopy cover, the differences in both understorey cover and mobile invertebrates assemblages associated with canopy-cover treatments were not consistent amongst sites for any time point (PERMANOVA treatment × site interactions; P(perm) < 0.05 for all analyses). In several cases differences between treatments were, at least in part, due to differences in dispersion amongst treatment groups at various time points for both understorey cover (PERMDISP by treatments P(perm) < 0.005 for days 5, 15/19 and 371/379 post manipulation) and mobile invertebrate assemblages (PERMDISP by treatments P(perm) < 0.03 for days 5 and 371/379 post manipulation; Figure A2, Figure A3).
One day after canopy manipulations, the understorey cover assemblages were significantly different from each other for all pairwise treatment comparisons except for naturally full, patchy and experimentally thinned treatments at 13th Beach, and the plots with the experimentally removed Hormosira canopy at Bridgewater Bay did not significantly differ from either the patchy or thinned plots (Table 2a, Figure A2). By the end of the experiment at 13th Beach (day 379), significant differences in understorey assemblages amongst all treatments except naturally full, patchy and experimentally thinned treatments remained (Table 2a, Figure A2). At the second experimental site (Bridgewater Bay) however, by the end of the experiment (day 371), the understorey cover assemblages in plots with naturally full canopy differed significantly from all treatments except the experimentally thinned canopy, which had converged on the full treatment by day 5 of the experiment (Table 2, Figure A2). At this site (Bridgewater Bay), the understorey cover assemblages in the patchy treatment differed significantly from only the full treatment by day 371, but the treatments largely devoid of Hormosira cover (negligible natural cover and experimentally removed) had significantly different understorey cover components to all other treatments (Table 1, Figure A2). Overall the understorey components that contributed most to the separation of treatment groups (as identified by SIMPER, Table A2) at the end of the experiment were the cover of the turfing C. implexa and small mussel X. pulex (individually and in association with each other as an intertwined complex), the amount of sand covering the substratum and bare rock devoid of visible biota (Figure 5).
The abundance of turfing alga Capreolia implexa and small mussel Xenostrobus pulex (combined in the analysis due to common occurrence as an intertwined complex) showed a consistent pattern throughout the experiment (Table 1, Figure 5). Both species were consistently most abundant in the negligible plots (75% at 13th Beach and ~55% at Bridgewater Bay, Table 1, Figure 5 at ~20% cover in patchy plots and either absent or < 3% cover in all other treatments at both sites (Table 1, Figure 5). Sand cover varied between treatments and sites throughout time, with higher abundance found in full plots when differences occurred (Table 1, Figure 5). Sand depth did not correspond with cover and on many occasions, deeper sand was recorded in plots with a lower overall cover of sand (Table 1 and Figure A4). Deeper sand (0.8–0.4 cm 13th Beach, 1–0.3 cm Bridgewater Bay) was most frequently recorded in full and patchy treatments (Table 1 and Figure A4). Bare rock cover varied through time, between treatments and sites (Table 1, Figure 5). Removed treatments had the most consistent cover of bare rock throughout the experiment (~30% 13th Beach, ~60% Bridgewater Bay).
For the mobile invertebrate assemblages associated with the different levels of Hormosira-canopy cover, there were significant differences in assemblages amongst all treatments one day after manipulation, except patchy vs. both full and thinned at 13th Beach and thinned vs. both patchy and full at Bridgewater Bay (Table 2, Figure 6). However by the end of the experiment, whilst there were no significant differences in the mobile invertebrate assemblages associated with full, patchy and thinned treatments at 13th Beach, full and patchy treatments were significantly different with respect to mobile invertebrates at Bridgewater Bay (Table 2, Figure 6). Mobile invertebrate assemblages associated with the naturally negligible Hormosira canopy and experimentally removed canopy were significantly different from most other canopy treatments over the duration of the experiment for both sites Table 2, Figure 6 and was principally driven (SIMPER, Table A3) by the abundance of the turbinid Lunella undulata, the trochid Chlorodiloma adelaidae and predatory whelk Dicathais orbita in treatments with a full or intermediate Hormosira canopy present; and juvenile littorinids, B. nanum and abundant limpets, in treatments lacking a Hormosira canopy (Figure 6, juvenile littorinids and limpets not shown).
The turbinid Lunella undulata was consistently most abundant in the full canopy plots at both sites through time (Table 1, Figure 6), with approximately double the average abundance at Bridgewater Bay compared with 13th Beach. Lunella undulata was found in patchy and thinned canopy plots at Bridgewater Bay but not at 13th Beach (Table 1, Figure 6). None were recorded in negligible or removed canopy treatments at either site (Table 1, Figure 6). The trochid Chlorodiloma adelaidae, was most abundant in the full canopy plots at both sites through time, though more were found in patchy and thinned treatments at 13th Beach during day 15 sampling (Table 1, Figure 6). The predatory whelk Dicathais orbita varied in abundance between treatments and time, though was most often more abundant in patchy, full and thinned canopy treatments and in higher overall abundance at 13th Beach (Table 1, Figure 6). The grazing gastropod Bembicium nanum varied in abundance between treatments through time at both sites, mostly in highest abundance in negligible canopy treatment plots (Table 1, Figure 6).
Species richness ranged from an average of 4.5–9 at 13th Beach and 4.5–10.5 at Bridgewater Bay, differing between treatments and sites through time with no consistent pattern (Table 1, Figure 7). The temperature remained lower under full canopies (average 20 °C) compared with all other treatments on all but one occasion at one site (day 5, Bridgewater Bay) during the experiment (Table 1, Figure 8). Negligible and/or removed canopy treatments had the highest temperatures (average 24 °C) with patchy and thinned canopy plots having intermediate temperatures (average 22 °C and 21 °C respectively, Table 1, Figure 8). Light intensity was consistently lower in full canopies compared with all other treatments throughout the experiment (Table 1, Figure 9). Negligible canopy plots usually had the highest relative levels of light, though this wasn’t always statistically significant and is likely influenced by intermittent cloud cover during sampling (Table 1, Figure 9). Light levels in patchy and thinned treatments were generally intermediate between the two extremes of canopy cover (Table 1, Figure 9). Hormosira banskii frond length was found to be longer in full canopy plots at both sites (average 14 cm 13th Beach, 11 cm at Bridgewater Bay, Table 1, Figure A5). Frond length averaged 9 cm in patchy and thinned plots at 13th beach; 10 and 7 cm at Bridgewater Bay, respectively (Table 1, Figure A5). The fronds found in negligible and removed canopy plots never exceeded 2 cm throughout the course of the experiment (Figure A5).

4. Discussion

4.1. Hormosira banksii as An Ecosystem Engineer

The results of this study support the suggestion of Keough and Quinn [5] that ‘established’ H. banksii canopies function as ecosystem engineers on Victorian rocky shores by providing environmentally buffered refuges (cool and shaded) for particular species (see also References [2,9,12,21,27]). We found differences in assemblages associated with different natural states of the canopy (full, patchy, negligible cover) across four sites; and like other studies [2,5,9,12,21] a loss of species occurred with the removal and part-removal of full canopies at two sites. Canopies that were naturally patchy did not function in the same way as full canopies, the environmental parameters measured under these canopies and the associated species assemblages were generally intermediate between those found under full canopies and in the absence of a canopy. This study builds on the findings of studies undertaken in Victoria and New Zealand [2,5,9,12,17,21] by examining the role of H. banksii canopies not only in response to disturbance but also on areas of the shore that have naturally patchy canopies. This demonstrates that the assemblages associated with patchy canopies result from the reduced canopy cover and not disturbance per se.
The understorey and macroinvertebrate assemblages in areas of naturally negligible canopy cover (which were between or adjacent to canopy habitats with no differences in rock type, exposure or height on shore should) differed to those where the canopy was removed. Naturally, negligible cover areas on the shore were dominated by grazing snails including limpets and periwinkles, and often aggregations of the mussel Xenostrobus pulex and the turfing alga Capreolia implexa. Since the disturbance treatments did not have assemblages that matched the natural controls (negligible areas) it suggests that it takes longer than 1 year for an alternative community such as the Xenostrobus pulex and Capreolia implexa complex to establish. This is consistent with studies showing variable recruitment (e.g., [37]).
The physical conditions and associated assemblages within thinned canopy treatments and patchy canopies varied throughout time and between sites with negligible and full canopy areas varying from discrete to similar (consistent with studies in South Australia [26] and in New Zealand [12]). As found in New Zealand [9] no replacement of H. banksii occurred in this study. In contrast to similar studies in the northern hemisphere (e.g., [7]), no other dominant canopy forming perennial algae extend into the intertidal region of Victorian rocky shores. In contrast to Lilley and Schiel [9], other fucoids were rarely recorded in the understorey of H. banksii and were generally only observed in very low areas of the shore, rarely exposed by tides (e.g., Cystophora torulosa). This may be due to differences in H. banskii frond length (up to 18 cm in our study, up to 25 in Lilley and Schiel [9], see also Reference [38]) or climatic differences in extreme temperatures during tidal emersion between Victorian and southern New Zealand shores. The results of this study demonstrate that dense and patchy canopies support different species assemblages and have differing effects on the physical environment. Thus although H. banksii dominates reef platforms, the magnitude of its ecosystem engineer function for associated communities is influenced by the amount of canopy cover.

4.2. Species Associations with Cover

This study found that particular species were found to be associated with full canopies, patchy canopies and areas where canopies were absent from the shore. Some of these species were the same as those found in other studies such as Lunella undulata (as Turbo undulatus) in full canopies [21], and Bembicium nanum in patchy canopies [24]. Lunella undulata was generally only found within canopies or in submerged areas of the rock platform during low tide throughout the study, which suggests that they are positively associated with the buffered conditions supplied by a full canopy. Interestingly species richness did not vary consistently between treatments during the experiment despite the differences observed in species assemblages between treatments. This demonstrates the power of recording species-specific data and using multivariate analysis when investigating autogenic ecosystem engineers [10]. Our results differed from a Victorian study examining trampling disturbance and its effect on H. banksii assemblages [21]. King’s study found background variation in the abundance of L. undulata and Chlorodiloma odontis (as Austrocochlea odontis) was stronger than the trampling intensity [21]. In our study, L. undulata and C. adelaidae consistently differed in abundance between treatments indicating disturbance and canopy cover investigations are not equal. Although H. banksii canopies studied by Underwood [24] correspond to the ‘patchy’ areas in this study (similar frond length and cover), he found no difference in the abundance of the grazing gastropod B. nanum between canopy and non-canopy treatments. In our study, they were most abundant in negligible areas on all except one sampling day suggesting either site-specific/latitudinal effects or that this species has different habitat preferences in the absence of disturbance.
Despite the lack of coralline turf, the fleshy turf-forming alga Capreolia implexa was common in both patchy (all sites) and negligible-cover plots (at both experimental sites) and was observed to form large mats discrete from H. banksii often in conjunction with the mussel Xenostrobus pulex. Hormosira banksii was rarely observed growing within the C. implexa turf suggesting this species may be competitively excluding H. banksii from growing on some areas of the shore (as found for coralline turf elsewhere, [39]), though this would need to be tested experimentally.

4.3. Alternate Habitat Types and Stable States

What was evident during this study was that discrete habitats are occurring on Victorian rocky shores, made up of stable, full-canopy H. banksii beds, areas of patchy (possibly disturbed) areas of canopy mixed with a range of other habitats, and ‘bare’ (negligible canopy) areas that consist of rock with abundant grazers, and substrate dominated by aggregating mussels, or turfing algae. Despite such habitats commonly referred to in the literature as vertical ‘zones’, the results of this study support the earlier refutation of vertical zonation on temperate Australian shores by Underwood [40], as all these discrete habitats occurred without variation in exposure, height on shore or other observable abiotic factors. The multiple habitats occurring on the shores investigated in this study and their different species assemblages support the concept of multiple stable states (though planned manipulative experiments are necessary to test this [41]). The differences in assemblages associated with patchy or dense H. banksii cover demonstrates that it functions as an autogenic ecosystem engineer at least at the covers investigated. Since removed and negligible-cover plots supported alternative assemblages, it is evident that not all species on the shore are dependent on canopy cover as habitat, and its possible species such as turfing algae and aggregating mussels could also be functioning as ecosystem engineers (singularly or providing cumulative effects e.g., Reference [42]), although not at the same scale as the canopy alga e.g., mussels sampled during this study were ≤ 1 cm in height. It was clear that the presence of these multiple habitats across a shore allows a higher overall shore biodiversity than would exist with only one or few habitats, even extensive and dense H. banksii. Since predictable assemblages were associated with easily identifiable habitat areas, rapid assessment of reef health could be possible if applied to areas of high conservation value.

4.4. The Influence of Sand

One aspect of the study that added a level of complexity to the ecosystem engineering function of H. banksii, in comparison to the literature on intertidal canopy algal assemblages, was the abundance of sand accumulated under H. banksii canopy. King [21] found that sand accumulated into H. banksii mats during winter and spring and suggested it could be a cause of natural disturbance decreasing H. banksii cover. Airoldi and Hawkins [43] found that sediment accumulation on rocky shores in the United Kingdom reduced the grazing and even survival of the limpet Patella vulgata (commonly associated with algal canopies). Given the findings of Airoldi and Hawkins [43], the accumulation of sediment under H. banksii canopies in this study may explain the lack of limpets found in the understorey in contrast to studies in the northern hemisphere. Sand accumulation may have also led to the underestimation of some understorey species such as encrusting coralline algae. Encrusting coralline algae wasn’t found to differ in abundance between treatments but was recorded in the removed treatments at both sites in the short term following disturbance before being, lost presumably due to high light and/or temperature exposure. In many of the full canopies during this study sand was a major component of the understorey with depths often 3–4 mm at both experimental sites, precluding identification of encrusting and very small species buried by the sand. In terms of negative effects on H. banksii, sand may cause loss of H. banksii cover due to the smothering of lower vesicles (observed at several sites) making it susceptible to breakage [44], sand has been shown to also prevent the attachment of H. banksii recruits [45]. It is possible high sand accumulation in H. banksii canopies may be contributing to the comparatively lower abundances of understorey algae and invertebrates recorded in observations of other canopies such as Ascophyllum nodosum in the northern hemisphere where sand cover was not found in significant amounts (see [46]). Further well-planned investigations into the role of sand on intertidal rocky reefs [47] may well provide valuable insight into this finding.

5. Conclusions

Full canopies (within the vertical range of H. banksii) provided habitat for species that are not found elsewhere on the reef at four sites. Low canopy cover (both natural and disturbed) resulted in different species than found within an intact canopy, correlating with a reduction in the ability of the canopy to buffer temperature and light stress. Patchy canopies (between 40%–60% cover) have species assemblages and physical conditions intermediate between full canopies and areas where H. banksii was absent on the shore.
Previous studies into canopy-forming H. banksii on intertidal rocky shores have highlighted the problem of background variation preventing clarity in their conclusions [9,10,12,21]. By examining the background variation in support of the response of the associated assemblage to disturbed areas of H. banksii, this study clearly demonstrates that H. banksii canopies function as autogenic ecosystem engineers on Victorian rocky intertidal shores, increasing diversity at the shore scale. Despite overall similar outcomes to that found in other studies of H. banksia [2,9,12,21,25], the species assemblages found here did differ from those reported elsewhere supporting the suggestion by Lilley and Schiel [9] that “superficially similar communities in different geographic regions may have different ecological mechanisms affecting structure”. These results suggest that site-specific evaluation and monitoring of intertidal communities are required to adequately detect the stability and health of rocky shores in relation to canopy algae. Our results demonstrate that to maintain high biodiversity at the seascape level on these shores, a combination of dense and patchy H. banksii canopy in addition to areas devoid of H. banksii canopy is ideal.

Author Contributions

Conceptualization, J.B.P., A.B.; Methodology, J.B.P.; Software, J.B.P., T.D.O. and A.B.; Validation, J.B.P., A.B., M.J.K., T.D.O.; Formal Analysis, J.B.P. and A.B.; Investigation, J.B.P.; Resources, T.D.O., M.J.K. and A.B.; Data Curation, J.B.P. and A.B.; Writing-Original Draft Preparation, J.B.P.; Writing-Review and Editing, J.B.P., A.B., T.D.O., M.J.K.; Visualization, J.B.P.; Supervision, A.B., T.D.O.; M.J.K.; Project Administration, J.B.P. and T.D.O.; Funding Acquisition, T.D.O.

Funding

This research was funded by an Australian National Heritage Trust grant awarded to T.D.O, project no: 202244. JBP was funded by a studentship under this grant and a Melbourne University Graduate Research scholarship.

Acknowledgments

Acknowledgement is given to Rebecca Citroen, Sera Cutler, Trudy Costa, Robert Pocklington and John Pocklington for field assistance. Gerry Kraft assisted with the identification of algae. The Marine Research Group (The Field Naturalists Club of Victoria) assisted with invertebrate identification. Dane Balodis created the map of study sites. Trudy Costa, Prue Addison, and Liz Greaves provided discussion and feedback throughout the study. Drafts of this manuscript were improved by feedback from Jorge Gutierrez and two anonymous reviewers. Research in National Parks was carried out under permits (10003820 & 10003879) pursuant to the National Parks Act 1975, research at other sites was carried out under permit (RP699) pursuant to the Fisheries Act 1995.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Photo pane showing dense (top), patchy (middle) and absent/negligible (bottom) areas occurring on rocky seashores in Victoria and as sampled. The scale of each image is approximately 20 cm across.
Figure A1. Photo pane showing dense (top), patchy (middle) and absent/negligible (bottom) areas occurring on rocky seashores in Victoria and as sampled. The scale of each image is approximately 20 cm across.
Diversity 11 00079 g0a1
Table A1. Species and component list by site and type (> 78 species).
Table A1. Species and component list by site and type (> 78 species).
Bridgewater Bay13th BeachSorrento Back BeachPoint Roadknight
Primary AlgaeCapreolia implexa
Cladophora sp.
Coralline Turf
Hormosira banksii
Scytosiphon lomentaria
Ulva sp.
Capreolia implexa
Ceramium sp.
Cladophora sp.
Colpomenia sinuosa
Coralline Turf
Encrusting Coralline
Hormosira banksii
Sargassum sp.
Laurencia filiformis
Laurencia spp.
Ulva sp.
Unidentified Reds
Capreolia implexa
Cladophora sp.
Coralline Turf
Hormosira banksii
Laurencia sp.
Rivularia sp.
Ulva sp.
Capreolia implexa
Coralline turf
Diplocladia patersonis
Encrusting Coralline
Hormosira banksii
Notheia anomala
Scytosiphon lomentaria
Ulva sp.
Understorey Algae
‘Coralline Turf’ is likely Corallina officinalis but due to compact form was not able to be positively identified in the field.
Caulocystis cephalornithos
Ceramium sp.
Chaetomorpha sp.
Cladophora sp.
Colpomenia sinuosa
Coralline Turf
Diplocladia patersonis
Encrusting Coralline
Gelidium australe
Laurencia botryoides
Laurencia filiformis
Laurencia spp.
Notheia anomala
Scytosiphon lomentaria
Sargassum sp.
Ulva sp.
Unidentified Encrusting Black
Unidentified Epiphytes
Unidentified Reds
Encrusting green
Capreolia implexa
Caulerpa brownii
Ceramium sp.
Chaetomorpha sp.
Cladophora sp.
Colpomenia sinuosa
Coralline Encrusting
Coralline Turf
Cystophora spp.
Laurencia botryoides
Laurencia filiformis
Laurencia spp.
Notheia anomala
Pachydictyon sp.
Sargassum sp.
Ulva sp.
Unidentified Reds
Capreolia implexa
Ceramium sp.
Cladophora sp.
Cystophora sp.
Coralline Turf
Diplocladia patersonis
Encrusting Coralline
Laurencia sp.
Rivularia sp.
Ulva sp.
Capreolia implexa
Ceramium sp.
Cladophora sp.
Coralline Turf
Diplocladia patersonis
Encrusting Coralline
Laurencia sp.
Rivularia sp.
Splachnidium rugosum
Unidentified Reds
Invertebrates
BivalvesBrachidontes rostratus
Xenostrobus pulex
Unidentified Bivalve
Xenostrobus pulex
Unidentified Bivalve
Xenostrobus pulexXenostrobus pulex
PolycheatesGaleolaria caespitosaGaleolaria caespitosaGaleolaria caespitosaGaleolaria caespitosa
Limpets
(and false limpets)
Cellana tramoserica
Lottia mixta
Montfortula rugosa
Notoacmea alta
Notoacmea flammea
Onchidella nigricans
Patelloida alticostata
Patelloida latistrigata
Siphonaria diemenensis
Siphonaira zelandica
Juvenile limpets
(unidentified)
Cellana tramoserica
Lottia mixta
Notoacmea flammea
Notoacmea sp.
Onchidella nigricans
Patelloida alticostata
Siphonaria diemenensis
Siphonaria zelandica
Juvenile limpets
(unidentified)
Cellana tramoserica
Lottia mixta
Notoacmea alta
Notoacmea flammea
Notoacmea petterdi
Notoacmea sp.
Patelloida alticostata
Patelloida latistrigata
Siphonaria diemenensis
Siphonaria funiculata
Siphonaria zelandica
Juvenile limpets
(unidentified)
Cellana tramoserica
Montfortula rugosa
Notoacmea sp.
Notoacmea flammea
Notoacmea petterdi
Patella peroni
Patelloida alticostata
Patelloida latistrigata
Siphonaria diemenensis
Siphonaria zelandica
Siphonaria funiculata
Juvenile limpets
(unidentified)
GastropodsAustrocochlea constricta
Austrocochlea porcata
Austrolittornia unifasciata
Bembicium nanum
Cantharidius pulcherimus
Chlorodiloma adelaidea
Chlorodiloma odontis
Cominella lineolata
Dicathais orbita
Littorinid juvenile
(unidentified)
Mitre sp.
Phasianotrochus eximus
Thalotia conica
Lunella undulata
Austrocochlea constricta
Bembicium nanum
Cantharidius pulcherimus
Chlorodiloma adelaidea
Chlorodiloma odontis
Cominella lineolata
Cystophora sp.
Dentrimitrella pulla
Dicathais orbita
Lepsellia vinosa
Littorinid juvenile
(unidentified)
Mitre cabonaria
Phasianotrochus eximus
Thalotia conica
Lunella undulata
Waimatea obscura
Austrocochlea constricta
Austrolittorina unifasciata
Bembicium nanum
Chlorodiloma adelaidae
Cominella lineolata
Dicathais orbita
Littorinid juvenile
(unidentified)
Nerita atramentosa
Lunella undulata
Austrocochlea constricta
Austrocochlea juvenile
Bembicium nanum
Chlorodiloma adelaidae
Cominella lineolata
Dicathis orbita
Lepsiella vinosa
Littorinid juvenile
(unidentified)
Mitre sp.
Nerita atramentosa
Lunella undulata
EchinodermsMeridiastra calcar
Parvulastra exigua
Meridiastra calcarMeridiastra calcar
Parvulastra exigua
CrustaceansChthalaumus antennatus
Paragrapsus quadradentatus
Guinusia chabrus
Naxia aurita
Unidentified Amphipods
Unidentified Copepods
Unidentified Harpactocoids
Unidentified Isopods
Chthalaumus antennatus
Naxia aurita
Nectocarcinus turbculosa
Paragrapsus quadradentatus
Guinusia chabrus
Unidentified Amphipods
Unidentified Copepods
Unidentified Harpactocoids
Unidentified Isopods
Cthalaumus antennatus
Paragrapsus quadradentatus
Paragrapsus sp.
Unidentified Isopods
Unidentified Amphipods
Chthalaumus antennatus
Paragrapsus quadradentatus
Unidentified Isopods
AnemonesActinia tenebrosa
Anthothoe albocinta
Aulactinia verartra
Oulactis mucosa
Oulactis sp.
Unidentified anemones
Anthothoe albocinta
Oulactis mucosa
Oulactis sp.
Oulactis mucosa
Oulactis sp.
Unidentified Anemone
Actinia tenebrosa
Aulactinia verartra
ChitonsUnidentified ChitonUnidentified ChitonUnidentified Chiton
Figure A2. Non-metric multidimensional scaling (nMDS) plots of understorey cover assemblages associated with different levels (full cover = dark blue, patchy cover = aqua, negligible cover = green, thinned = pink, removed = grey) of Hormosira cover (a) 1 day, (b) 5 days, (c) 15/19 days and (d) 371/379 days following manipulation at each of two sites (13th Beach = triangles, Bridgewater Bay = circles). All data were fourth-root transformed. 3D nMDS plots based on Bray-Curtis dissimilarities are shown where stress exceeded 0.20 for 2D plots.
Figure A2. Non-metric multidimensional scaling (nMDS) plots of understorey cover assemblages associated with different levels (full cover = dark blue, patchy cover = aqua, negligible cover = green, thinned = pink, removed = grey) of Hormosira cover (a) 1 day, (b) 5 days, (c) 15/19 days and (d) 371/379 days following manipulation at each of two sites (13th Beach = triangles, Bridgewater Bay = circles). All data were fourth-root transformed. 3D nMDS plots based on Bray-Curtis dissimilarities are shown where stress exceeded 0.20 for 2D plots.
Diversity 11 00079 g0a2
Figure A3. Non-metric multidimensional scaling (nMDS) plots of mobile invertebrate assemblages associated with different levels (full cover = dark blue, patchy cover = aqua, negligible cover = green, thinned = pink, removed = grey) of Hormosira cover (a) 1 day, (b) 5 days, (c) 15/19 days and (d) 371/379 days following manipulation at each of two sites (13th Beach = triangles, Bridgewater Bay = circles). All data were fourth-root transformed. 3D nMDS plots based on Bray-Curtis dissimilarities are shown where stress exceeded 0.20 for 2D plots.
Figure A3. Non-metric multidimensional scaling (nMDS) plots of mobile invertebrate assemblages associated with different levels (full cover = dark blue, patchy cover = aqua, negligible cover = green, thinned = pink, removed = grey) of Hormosira cover (a) 1 day, (b) 5 days, (c) 15/19 days and (d) 371/379 days following manipulation at each of two sites (13th Beach = triangles, Bridgewater Bay = circles). All data were fourth-root transformed. 3D nMDS plots based on Bray-Curtis dissimilarities are shown where stress exceeded 0.20 for 2D plots.
Diversity 11 00079 g0a3
Table A2. Summary SIMPER table: Average dissimilarity of treatment pairs for understorey cover combining Bridgewater Bay and 13th Beach. Contributions > 10% given by understorey component/species.
Table A2. Summary SIMPER table: Average dissimilarity of treatment pairs for understorey cover combining Bridgewater Bay and 13th Beach. Contributions > 10% given by understorey component/species.
Negligible, FullNegligible, PatchyNegligible, RemovedNegligible, ThinnedFull, PatchyFull, RemovedFull, ThinnedPatchy, RemovedPatchy, ThinnedRemoved, Thinned
Day 1
% contribution
Average dissimilarity85.6986.7992.1995.8071.5093.7986.4689.4486.4679.29
Sand31.0624.75 11.3941.6434.3235.0323.1135.0334.13
Bare rock24.28 25.1014.4916.6223.9213.0523.92
Xenostrobus pulex22.0721.2220.5633.77
Capreolia implexa 10.43 12.93 12.93
Capreolia implexa and Xenostrobus pulex complex15.0613.7114.1113.99
Hormosira holdfasts 18.10 17.87 18.77 34.13
Day 5
% contribution
Average dissimilarity83.0668.0776.8983.5067.4865.4256.8263.2058.4458.43
Sand 10.01 12.00
Bare rock24.2626.68 22.7423.8026.5414.6923.8629.5427.27
Xenostrobus pulex 11.70 10.59
Capreolia implexa 10.21
Capreolia implexa and Xenostrobus pulex complex11.00 10.24
Hormosira holdfasts 14.32 14.17 14.59 15.47
Encrusting Coralline 10.68 11.11
Day 15/19
% contribution
Average dissimilarity89.4289.2989.1187.7577.8776.0663.5082.1079.8070.73
Sand 10.92 16.44 12.3513.83
Bare rock44.7628.5133.7343.1937.6152.3555.8635.5949.5555.25
Xenostrobus pulex20.1424.3120.6123.00
Capreolia implexa
Capreolia implexa and Xenostrobus pulex complex15.3114.1515.42 14.66
Hormosira holdfasts 15.20 17.96 17.70 15.68
Encrusting Coralline 10.9910.7013.59
Day 371/379
% contribution
Average dissimilarity74.7765.2174.8470.4466.2382.6743.7273.2061.1768.96
Sand15.7913.54 13.5623.9116.27 11.5021.0816.26
Bare rock15.2010.61 11.5016.2015.4910.7013.5613.3613.76
Xenostrobus pulex10.5812.59
Capreolia implexa 12.29 10.50 13.04
Capreolia implexa and Xenostrobus pulex complex18.1420.8912.8117.40
Hormosira holdfasts 10.21 11.71 17.40 12.54
Encrusting Coralline 11.18 15.70 11.15
Cladophora sp. 10.37 13.85
Galeolaria caespitosa 13.07
Coralline Turf 11.72
Table A3. Summary SIMPER table: Average dissimilarity of treatment pairs for mobile macroinvertebrates combining Bridgewater Bay and 13th Beach. Contributions > 10% given by species.
Table A3. Summary SIMPER table: Average dissimilarity of treatment pairs for mobile macroinvertebrates combining Bridgewater Bay and 13th Beach. Contributions > 10% given by species.
Negligible, FullNegligible, PatchyNegligible, RemovedNegligible, ThinnedFull, PatchyFull, RemovedFull, ThinnedPatchy, RemovedPatchy, ThinnedRemoved, Thinned
Day 1
% contribution
Average dissimilarity10087.0895.1510087.8297.6084.5194.2278.4192.31
Bembicium nanum13.3714.7324.1317.73
Juvenile littorinid 13.5719.8414.48
Chlorodiloma adelaidae 11.9813.3521.6915.7713.3210.33
Dicathais orbita 10.37 14.6419.2215.7135.69
Lunella undulata 15.0614.35 12.40
Siphonaria diemenensis 13.24
Cominella lineolata 12.95
Day 5
% contribution
Average dissimilarity97.2287.2096.0410086.4099.4690.9397.9896.60100
Bembicium nanum14.7616.9530.9325.31 12.1610.15
Juvenile littorinid15.5218.2635.5527.24
Chlorodiloma adelaidae 12.3310.02
Dicathais orbita 12.99 12.0112.3310.5023.1919.25
Lunella undulata 10.9812.4011.00 13.02
Oulactis mucosa 15.9010.3414.59
Amphipoda 23.80
Anthothoe albocinta 11.90
Day 15/19
% contribution
Average dissimilarity99.5189.8397.0394.9081.6596.3985.7295.7386.6995.48
Bembicium nanum 10.0919.6513.43 16.17
Juvenile littorinid14.0614.7424.8719.41
Chlorodiloma adelaidae 11.8218.1816.94 10.86
Dicathais orbita 11.30 15.3515.3515.90
Lunella undulata 10.41 14.4415.9914.9620.6116.87
Juvenile limpet 18.0913.77
Cominella lineolata 12.4513.20 11.4718.27
Aulactinia veratra 10.2010.45
Day 371/379
% contribution
Average dissimilarity10084.8483.7897.9987.8197.1579.5188.7785.5289.96
Bembicium nanum14.2415.9613.4715.11
Juvenile littorinid21.2525.3621.0721.59
Chlorodiloma adelaidae11.43 17.0112.9723.57 15.6210.04
Dicathais orbita 11.19 11.5612.2916.7512.7314.11
Lunella undulata12.95 18.6614.5520.70
Siphonaria diemenensis 11.00 10.40 11.55 11.67
Cominella lineolata 11.70 13.89 15.47
Juvenile limpet 14.67 11.15 14.71 12.26
Thalotia conica 10.95
Figure A4. Average sand depth (cm ± se) between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Figure A4. Average sand depth (cm ± se) between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Diversity 11 00079 g0a4
Figure A5. Average Hormosira frond length (cm ± se) between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Figure A5. Average Hormosira frond length (cm ± se) between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Diversity 11 00079 g0a5

References

  1. Menge, B.A. Predation Intensity in a Rocky Intertidal Community: Effect of an algal canopy, wave action and desiccation on predator feeding rates. Oecologia 1978, 34, 17–35. [Google Scholar] [CrossRef]
  2. Povey, A.; Keough, M.J. Effects of trampling on plant and animal populations on rocky shores. Oikos 1991, 61, 355–368. [Google Scholar] [CrossRef]
  3. Jones, C.G.; Lawton, J.H.; Shachak, M. Organisms as ecosystem engineers. Oikos 1994, 69, 373–386. [Google Scholar] [CrossRef]
  4. Alper, J. Ecosystem ‘Engineers’ Shape Habitats for Other Species. Science 1998, 280, 1195–1196. [Google Scholar] [CrossRef]
  5. Keough, M.J.; Quinn, G.P. Effects of periodic disturbances from trampling on rocky intertidal algal beds. Ecol. Appl. 1998, 8, 141–161. [Google Scholar]
  6. Schiel, D.R.; Taylor, D.I. Effects of trampling on a rocky intertidal algal assemblage in southern New Zealand. J. Exp. Mar. Biol. Ecol. 1999, 235, 213–235. [Google Scholar]
  7. Jenkins, S.R.; Norton, T.A.; Hawkins, S.J. Long term effects of Ascophyllum nodosum canopy removal on mid shore community structure. J. Mar. Biol. Assoc. UK 2004, 84, 327–329. [Google Scholar] [CrossRef]
  8. Dudgeon, S.; Petraitis, P.S. First year demography of the foundation species, Ascophyllum nodosum, and its community implications. Oikos 2005, 109, 405–415. [Google Scholar] [CrossRef]
  9. Lilley, S.A.; Schiel, D.R. Community effects following the deletion of a habitat-forming alga from rocky marine shores. Oecologia 2006, 148, 672–681. [Google Scholar] [CrossRef]
  10. Schiel, D.R. Rivets or bolts? When single species count in the function of temperate rocky reef communities. J. Exp. Mar. Biol. Ecol. 2006, 338, 233–252. [Google Scholar] [CrossRef]
  11. Hastings, A.; Byers, J.E.; Crooks, J.A.; Cuddington, K.; Jones, C.G.; Lambrinos, J.G.; Talley, T.S.; Wilson, W.G. Ecosystem engineering in space and time. Ecol. Lett. 2007, 10, 153–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Schiel, D.; Lilley, S. Gradients of disturbance to an algal canopy and the modification of an intertidal community. Mar. Ecol. Prog. Ser. 2007, 339, 1–11. [Google Scholar] [CrossRef] [Green Version]
  13. Benedetti Cecchi, L.; Cinelli, F. Canopy removal experiments in Cystoseira-dominated rockpools from the Western coast of the Mediterranean (Ligurian Sea). J. Exp. Mar. Biol. Ecol. 1992, 155, 69–83. [Google Scholar] [CrossRef]
  14. Cervin, G.; Åberg, P. Do littorinids affect the survival of Ascophyllum nodosum germlings? J. Exp. Mar. Biol. Ecol. 1997, 218, 35–47. [Google Scholar] [CrossRef]
  15. Pavia, H.; Carr, H.; Åberg, P. Habitat and feeding preferences of crustacean mesoherbivores inhabiting the brown seaweed Ascophyllum nodosum (L.) Le Jol. and its epiphytic macroalgae. J. Exp. Mar. Biol. Ecol. 1999, 236, 15–32. [Google Scholar] [CrossRef]
  16. Coleman, M.A.; Clark, J.S.; Doblin, M.A.; Bishop, M.J.; Kelaher, B.P. Genetic differentiation between estuarine and open coast ecotypes of a dominant ecosystem engineer. Mar. Freshw. Res. 2018. [Google Scholar] [CrossRef]
  17. Wootton, H.F.; Keough, M.J. Disturbance type and intensity combine to affect resilience of an intertidal community. Mar. Ecol. Prog. Ser. 2016, 560, 121–133. [Google Scholar] [CrossRef]
  18. Tait, L.W.; Schiel, D.R. Legacy Effects of Canopy Disturbance on Ecosystem Functioning in Macroalgal Assemblages. PLoS ONE 2011, 6, e26986. [Google Scholar] [CrossRef] [PubMed]
  19. Raffaelli, D. The grazer-algae interaction in the intertidal zone of New Zealand rocky shores. J. Exp. Mar. Biol. Ecol. 1979, 38, 81–100. [Google Scholar] [CrossRef]
  20. Underwood, A.J. The vertical distribution and seasonal abundance of intertidal macroalgae on a rocky shore in New South Wales. J. Exp. Mar. Biol. Ecol. 1984, 78, 199–220. [Google Scholar] [CrossRef]
  21. King, A. Monitoring and Management of Human Activity on Rocky Shores. Ph.D. Thesis, Melbourne University, Melbourne, Australia, 1992. [Google Scholar]
  22. Bellgrove, A. Recruitment of Intertidal Macroalgae on a Wave-Exposed Rocky Coast; Monash Clayton: Melbourne, Australia, 1998. [Google Scholar]
  23. Ralph, P.J.; Morrison, D.A.; Addison, A. A quantitative study of the patterns of morphological variation within Hormosira banksii (Turner) Decaisne (Fucales: Phaeophyta) in south-eastern Australia. J. Exp. Mar. Biol. Ecol. 1998, 225, 285–300. [Google Scholar] [CrossRef]
  24. Underwood, A.J. Grazing and disturbance: An experimental analysis of patchiness in recovery from a severe storm by the intertidal alga Hormosira banksii on rocky shores in New South Wales. J. Exp. Mar. Biol. Ecol. 1998, 231, 291–306. [Google Scholar] [CrossRef]
  25. Underwood, A.J. Physical disturbances and their direct effect on an indirect effect: Responses of an intertidal assemblage to a severe storm. J. Exp. Mar. Biol. Ecol. 1999, 232, 125–140. [Google Scholar] [CrossRef]
  26. Pocklington, J.B. Positive Interactions on Rocky Seashores. Are They Important, or Can They Be Ignored? Honours Thesis, Flinders University of South Australia, Adelaide, Australia, 2003. [Google Scholar]
  27. Lilley, S.A. Removal of Habitat-Forming Species and the Consequences on Community Biodiversity in New Zealand Rocky Shore Ecosystems. Master’s Thesis, Canterbury University, Christchurch, New Zealand, 2004. [Google Scholar]
  28. Macinnis-Ng, C.M.O.; Morrison, D.A.; Ralph, P.J. Temporal and spatial variation in the morphology of the brown macroalga Hormosira banksii (Fucales, Phaeophyta). Bot. Mar. 2005, 48, 198–207. [Google Scholar] [CrossRef]
  29. Addison, P.F.E.; Koss, R.S.; O’Hara, T.D. Recreational use of a rocky intertidal reef in Victoria: Implications for ecological research and management. Aust. J. Environ. Manag. 2008, 15, 169–179. [Google Scholar] [CrossRef]
  30. Schoenwaelder, M.E.A. Physode distribution and the effect of ‘Thallus sunburn’ in Hormosira banksii (Fucales, Phaeophyceae). Bot. Mar. 2002, 45, 262–266. [Google Scholar] [CrossRef]
  31. Bishop, M.J.; Morgan, T.; Coleman, M.A.; Kelaher, B.P.; Hardstaff, L.K.; Evenden, R.W. Facilitation of molluscan assemblages in mangroves by the fucalean alga Hormosira banksii. Mar. Ecol. Prog. Ser. 2009, 392, 111–122. [Google Scholar] [CrossRef]
  32. Short, A.D. Beaches of the Victorian Coast & Port Phillip Bay: A Guide to Their Nature, Characteristics, Surf and Safety; University of Sydney & Surf Life Saving Victoria: Sydney, Australia, 2006. [Google Scholar]
  33. Bird, E. The Coast of Victoria: The Shaping of Scenery; Melbourne University Press: Carlton, Australia, 1993. [Google Scholar]
  34. Kruskal, J.B. Non-parametric multidimensional scaling: A numerical method. Psychometrika 1964, 29, 115–129. [Google Scholar] [CrossRef]
  35. Clarke, K.R. Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 1993, 18, 117–143. [Google Scholar] [CrossRef]
  36. Anderson, M.J.; Gorley, R.N.; Clarke, K.R. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods; Primer-E Limited: Plymouth, UK, 2008. [Google Scholar]
  37. Bellgrove, A.; Clayton, M.N.; Quinn, G.P. An integrated study of the temporal and spatial variation in the supply of propagules, recruitment and assemblages of intertidal macroalgae on a wave-exposed rocky coast, Victoria, Australia. J. Exp. Mar. Biol. Ecol. 2004, 310, 207–225. [Google Scholar] [CrossRef]
  38. Gemelli, F.; Johnson, C.R.; Wright, J.T. Gastropod communities associated with different morphologies of the intertidal seaweed Hormosira banksii. Mar. Freshw. Res. 2019, 70, 280–291. [Google Scholar] [CrossRef]
  39. Bellgrove, A.; McKenzie, P.; McKenzie, J.; Sfiligoj, B. Restoration of the habitat-forming fucoid alga Hormosira banksii at effluent-affected sites: Competitive exclusion by coralline turfs. Mar. Ecol. Prog. Ser. 2010, 419, 47–56. [Google Scholar] [CrossRef]
  40. Underwood, A.J. A refutation of critical tidal levels as determinates of the structure of intertidal communities on British shores. J. Exp. Mar. Biol. Ecol. 1978, 33, 185–202. [Google Scholar] [CrossRef]
  41. Petraitis, P.S.; Latham, R.E. The Importance of Scale in Testing the Origins of Alternative Community States. Ecology 1999, 80, 429–442. [Google Scholar] [CrossRef]
  42. Gutiérrez, J.; Bagur, M.; Palomo, M. Algal Epibionts as Co-Engineers in Mussel Beds: Effects on Abiotic Conditions and Mobile Interstitial Invertebrates. Diversity 2019, 11, 17. [Google Scholar] [CrossRef]
  43. Airoldi, L.; Hawkins, S.J. Negative effects of sediment deposition on grazing activity and survival of the limpet Patella vulgata. Mar. Ecol. Prog. Ser. 2007, 332, 235–240. [Google Scholar] [CrossRef] [Green Version]
  44. McKenzie, P.F.; Bellgrove, A. Dislodgment and attachment strength of the intertidal macroalga Hormosira banksii (Fucales, Phaeophyceae). Phycologia 2009, 48, 335–343. [Google Scholar] [CrossRef]
  45. Schiel, D.R.; Wood, S.A.; Dunmore, R.A.; Taylor, D.I. Sediment on rocky intertidal reefs: Effects on early post-settlement stages of habitat-forming seaweeds. J. Exp. Mar. Biol. Ecol. 2006, 331, 158–172. [Google Scholar] [CrossRef]
  46. Pocklington, J.B.; Jenkins, S.R.; Bellgrove, A.; Keough, M.J.; O’Hara, T.D.; Masterson-Algar, P.E.; Hawkins, S.J. Disturbance alters ecosystem engineering by a canopy-forming alga. J. Mar. Biol. Assoc. UK 2018, 98, 687–698. [Google Scholar] [CrossRef]
  47. Airoldi, L. The effects of sedimentation on rocky coast assemblages. In Oceanography and Marine Biology, An Annual Review; Gibson, R.N., Atkinson, R.J.A., Eds.; Aberdeen University Press/Allen & Unwin: London, UK, 2003; Volume 41, pp. 161–236. [Google Scholar]
Figure 1. Map showing study sites (black filled circles) including latitudes and longitudes. Bridgewater Bay and 13th Beach are experiment sites, Sorrento and Point Roadknight are background sites. All sites on the coast of the Southern Ocean, Bass Strait is directly south (not shown).
Figure 1. Map showing study sites (black filled circles) including latitudes and longitudes. Bridgewater Bay and 13th Beach are experiment sites, Sorrento and Point Roadknight are background sites. All sites on the coast of the Southern Ocean, Bass Strait is directly south (not shown).
Diversity 11 00079 g001
Figure 2. Non-metric multidimensional scaling (nMDS) plots of (a) primary cover, (b) understorey cover and (c) mobile invertebrates assemblages associated with different natural levels of Hormosira cover prior to manipulation (full cover = dark blue, patchy cover = aqua, negligible cover = green) at each of four sites (13th Beach = triangles, Bridgewater Bay = circles, Point Roadknight = squares, Sorrento = diamonds). All data were fourth-root transformed. 3D nMDS plots are shown where stress exceeded 0.20 for 2D plots.
Figure 2. Non-metric multidimensional scaling (nMDS) plots of (a) primary cover, (b) understorey cover and (c) mobile invertebrates assemblages associated with different natural levels of Hormosira cover prior to manipulation (full cover = dark blue, patchy cover = aqua, negligible cover = green) at each of four sites (13th Beach = triangles, Bridgewater Bay = circles, Point Roadknight = squares, Sorrento = diamonds). All data were fourth-root transformed. 3D nMDS plots are shown where stress exceeded 0.20 for 2D plots.
Diversity 11 00079 g002
Figure 3. (a) Average ± SE densities of mobile invertebrates (per 0.09 m2 quadrat) amongst natural Hormosira cover treatments for each of four sites: (a) Chlorodiloma adelaidae, (b) Lunella undulata, (c) Cominella lineolata, (d) Dicathais orbita, (e) Bembicium nanum, (f) Limpets (Siphonaria spp., Cellana tramoserica, Patelloida spp. combined). N = 5 for all treatments at 13th Beach and Bridgewater Bay and N = 10 for Point Roadknight and Sorrento.
Figure 3. (a) Average ± SE densities of mobile invertebrates (per 0.09 m2 quadrat) amongst natural Hormosira cover treatments for each of four sites: (a) Chlorodiloma adelaidae, (b) Lunella undulata, (c) Cominella lineolata, (d) Dicathais orbita, (e) Bembicium nanum, (f) Limpets (Siphonaria spp., Cellana tramoserica, Patelloida spp. combined). N = 5 for all treatments at 13th Beach and Bridgewater Bay and N = 10 for Point Roadknight and Sorrento.
Diversity 11 00079 g003
Figure 4. Average ± SE percentage cover of Hormosira banksii canopy over time and between treatments at 13th Beach (left) and Bridgewater Bay (right). N = 5 for all treatments.
Figure 4. Average ± SE percentage cover of Hormosira banksii canopy over time and between treatments at 13th Beach (left) and Bridgewater Bay (right). N = 5 for all treatments.
Diversity 11 00079 g004
Figure 5. Average ± SE percentage cover of (a) C. implexa and X. pulex, (b) sand (c) bare rock for each treatment over time since manipulation of the canopy. 13th Beach on left hand side of panel, Bridgewater Bay on right hand side. N = 5 for all treatments.
Figure 5. Average ± SE percentage cover of (a) C. implexa and X. pulex, (b) sand (c) bare rock for each treatment over time since manipulation of the canopy. 13th Beach on left hand side of panel, Bridgewater Bay on right hand side. N = 5 for all treatments.
Diversity 11 00079 g005
Figure 6. Mean ± SE percentage cover of (a) Lunella undulata (b) Chlorodiloma adelaidae, (c) Dicathais orbita (d) Bembicium nanum, for each treatment over time since manipulation of the canopy. 13th Beach on left hand site of the panel, Bridgewater Bay on right hand side. N = 5 for all treatments.
Figure 6. Mean ± SE percentage cover of (a) Lunella undulata (b) Chlorodiloma adelaidae, (c) Dicathais orbita (d) Bembicium nanum, for each treatment over time since manipulation of the canopy. 13th Beach on left hand site of the panel, Bridgewater Bay on right hand side. N = 5 for all treatments.
Diversity 11 00079 g006
Figure 7. Average species richness ± se between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Figure 7. Average species richness ± se between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Diversity 11 00079 g007
Figure 8. Average temperature (°C) ± se between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Figure 8. Average temperature (°C) ± se between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Diversity 11 00079 g008
Figure 9. Average light intensity ± se as relative units between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Figure 9. Average light intensity ± se as relative units between treatments and over time at 13th Beach (left hand side) and Bridgewater Bay (right hand side).
Diversity 11 00079 g009
Table 1. Summary repeated measures analysis of variance (rmANOVA) results comparing factors of interest and abundances of individual taxa and understorey cover components contributing significantly (identified by Similarity Percentage (SIMPER) to differences amongst Hormosira-canopy treatments, sites and time. Only times when sampling of the two sites were close to each other were included for analysis: days 1, 5, 15/19 and 371/379. MSResidual (as MSR) and degrees of freedom (df) are provided to allow reconstruction of the full ANOVA table. * denotes P-values that were adjusted to Greenhouse-Geisser Epsilon, † denotes data Square-Root transformed to meet assumptions of normality, n.s = not significant (p > 0.05).
Table 1. Summary repeated measures analysis of variance (rmANOVA) results comparing factors of interest and abundances of individual taxa and understorey cover components contributing significantly (identified by Similarity Percentage (SIMPER) to differences amongst Hormosira-canopy treatments, sites and time. Only times when sampling of the two sites were close to each other were included for analysis: days 1, 5, 15/19 and 371/379. MSResidual (as MSR) and degrees of freedom (df) are provided to allow reconstruction of the full ANOVA table. * denotes P-values that were adjusted to Greenhouse-Geisser Epsilon, † denotes data Square-Root transformed to meet assumptions of normality, n.s = not significant (p > 0.05).
FactorTime × Site × Treatment
df = 12
Time × Treatment
df = 12
Time × Site
df = 12
Time
df = 3
Site × Treatment
df = 4
Treatment
df = 4
Site
df = 1
Hormosira banksii canopy coverMSR = 179.473
p = 0.011
MSR = 588.142
p = 0.0001 *
n.sn.sn.sMSR = 76,844.5
p = 0.0001
n.s
Species Richnessn.sMSR = 12.151
p = 0.0001 *
n.sMSR = 63.880
p = 0.0001 *
n.sn.sn.s
UNDERSTOREY COVER
Xenostrobus pulex & Capriola implexan.sn.sn.sn.sn.sMSR = 30,174.2
p = 0.0001
n.s
Bare RockMSR = 1563.12
p = 0.002 *
MSR = 2127.67
p = 0.0001 *
n.sMSR = 7219.06
p = 0.0001 *
n.sn.sMSR = 12,059.0
p = 0.009
SandMSR = 1275.09
p = 0.037 *
MSR = 2983.68
p = 0.0001 *
n.sMSR = 10,801.3
p = 0.0001 *
n.sMSR = 10,821.3
p = 0.0001
n.s
MOBILE INVERTEBRATES
Lunella undulata†MSR = 0.611
p = 0.0001 *
MSR = 0.445
p = 0.003 *
MSR = 0.657
p = 0.009 *
n.sMSR = 4.386
p = 0.001
MSR = 8.106
p = 0.001
MSR = 2.948
p = 0.001
Chlorodiloma adelaidae†n.sn.sn.sn.sMSR = 5.445
p = 0.002
MSR = 15.755
p = 0.0001
n.s
Dicathais orbita†n.sMSR = 0.596
p = 0.004 *
n.sn.sn.sMSR = 3.880
p = 0.0001
MSR = 1.348
p = 0.035
Bembicium nanum†n.sMSR = 0.866
p = 0.004 *
n.sn.sn.sMSR = 8.752
p = 0.0001
n.s
PHYSICAL
TemperatureMSR = 9.834
p = 0.0001 *
MSR = 11.296
p = 0.0001 *
MSR = 589.470
p = 0.0001 *
MSR = 87.398
p = 0.0001 *
MSR = 27.602
p = 0.021
MSR = 156.467
p = 0.0001
MSR = 112.350
p = 0.001
Light intensityn.sn.sMSR = 5111.26
p = 0.0001 *
MSR = 9461.34 p = 0.0001 *n.sMSR = 40,416.1 p = 0.0001MSR = 10,353.6
p = 0.0001
COVARIATES
Frond lengthn.sMSR = 6.528
p = 0.015 *
n.sMSR = 22.376
p = 0.0001 *
n.sMSR = 1371.70 p = 0.0001MSR = 31.113
p = 0.037
Sand depthMSR = 0.512
p = 0.02 *
n.sn.sMSR = 1.377
p = 0.004 *
MSR = 1.149
p = 0.016
MSR = 1.653
p = 0.002
n.s
Table 2. Pairwise comparisons of (a) understorey cover and (b) mobile invertebrate assemblages between Hormosira-canopy treatment groups by sites following Permutational multivariate analysis of variance (PERMANOVA) on the significant treatment × site interaction terms. P(MC) is the P-value based on Monte Carlo random draws provides a more robust significance value with the reduced number of possible permutations in pairwise tests; statistical significance (α = 0.05) is indicated in bold.
Table 2. Pairwise comparisons of (a) understorey cover and (b) mobile invertebrate assemblages between Hormosira-canopy treatment groups by sites following Permutational multivariate analysis of variance (PERMANOVA) on the significant treatment × site interaction terms. P(MC) is the P-value based on Monte Carlo random draws provides a more robust significance value with the reduced number of possible permutations in pairwise tests; statistical significance (α = 0.05) is indicated in bold.
a. Understorey Cover
13th BeachDay 1Day 5Day 15Day 379
GroupstP(MC)tP(MC)tP(MC)tP(MC)
Negligible, Full2.73550.0053.10760.0013.47610.0033.07160.002
Negligible, Patchy2.2610.0112.06890.0162.5560.0071.63550.048
Negligible, Removed3.72940.0014.15160.0023.82150.0024.16950.001
Negligible, Thinned2.9990.0022.32860.0042.79930.0052.70440.002
Full, Patchy1.320.1751.71430.0611.52540.0991.45540.112
Full, Removed1.79880.0372.41510.0072.6840.0055.37620.001
Full, Thinned1.43720.1241.21420.231.42120.1441.32160.178
Patchy, Removed2.27440.011.9540.0261.9740.0172.39820.007
Patchy, Thinned1.58980.0940.288470.9561.48530.1081.00350.4
Removed, Thinned2.13270.021.61140.0822.02280.0353.49190.001
Bridgewater BayDay 1Day 5Day 19Day 371
GroupstP(MC)tP(MC)tP(MC)tP(MC)
Negligible, Full3.99280.0013.4720.0034.30750.0013.31290.003
Negligible, Patchy2.26830.011.7910.0552.97450.0070.996080.416
Negligible, Removed2.62920.0022.84060.0055.05780.0021.81230.029
Negligible, Thinned3.63410.0033.43660.0033.53780.0012.01450.017
Full, Patchy3.00280.0041.93090.0431.99640.0312.48890.009
Full, Removed3.15640.0022.55950.0072.07380.042.85310.005
Full, Thinned4.34950.0011.31480.1991.8570.081.00170.399
Patchy, Removed1.40450.141.50030.1111.77360.0531.35590.127
Patchy, Thinned2.36060.0131.29670.1871.72550.0721.37920.15
Removed, Thinned1.06780.3471.51060.0990.471410.8481.75020.033
b. Mobile Invertebrates
13th BeachDay 1Day 5Day 15Day 379
GroupstP(MC)tP(MC)tP(MC)tP(MC)
Negligible, Full2.10680.0042.47020.0022.350.0032.76480.004
Negligible, Patchy2.32730.0032.34180.0052.70970.0012.55550.001
Negligible, Removed3.03790.0026.0460.0013.74860.0012.88490.002
Negligible, Thinned3.62490.0012.71320.0032.46780.0052.79190.005
Full, Patchy1.15990.2631.32680.1280.806820.6590.877740.495
Full, Removed1.9930.0153.09280.0031.97190.0082.03660.015
Full, Thinned1.9820.0161.64990.0580.600270.8170.950890.452
Patchy, Removed2.05580.023.26910.0022.0050.0411.75230.04
Patchy, Thinned1.38350.1411.80720.031.02760.4071.13340.315
Removed, Thinned3.06430.0061.5120.1541.67270.0741.65860.076
Bridgewater BayDay 1Day 5Day 19Day 371
GroupstP(MC)tP(MC)tP(MC)tP(MC)
Negligible, Full3.83730.0023.32960.0022.85130.0022.87350.001
Negligible, Patchy1.74350.0282.02750.0161.58830.0520.843610.554
Negligible, Removed2.21460.022.27480.011.83330.0281.05690.366
Negligible, Thinned2.20840.0032.16930.0051.37210.1111.79790.02
Full, Patchy2.52340.0031.94010.0212.04770.0131.95660.016
Full, Removed2.44720.0062.33790.0032.2910.0191.96590.017
Full, Thinned1.4360.1111.61960.0691.85220.0151.34020.174
Patchy, Removed1.86080.0261.52590.1071.71040.030.926690.524
Patchy, Thinned1.36340.1321.41080.1231.19690.2310.989850.419
Removed, Thinned1.54590.0481.22010.2261.19470.241.20610.209

Share and Cite

MDPI and ACS Style

Pocklington, J.B.; Keough, M.J.; O’Hara, T.D.; Bellgrove, A. The Influence of Canopy Cover on the Ecological Function of A Key Autogenic Ecosystem Engineer. Diversity 2019, 11, 79. https://doi.org/10.3390/d11050079

AMA Style

Pocklington JB, Keough MJ, O’Hara TD, Bellgrove A. The Influence of Canopy Cover on the Ecological Function of A Key Autogenic Ecosystem Engineer. Diversity. 2019; 11(5):79. https://doi.org/10.3390/d11050079

Chicago/Turabian Style

Pocklington, Jacqueline B., Michael J. Keough, Tim D. O’Hara, and Alecia Bellgrove. 2019. "The Influence of Canopy Cover on the Ecological Function of A Key Autogenic Ecosystem Engineer" Diversity 11, no. 5: 79. https://doi.org/10.3390/d11050079

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop