1. Introduction
As climate change continues to exacerbate the effects of anthropogenic impacts across global land and seascapes, restoration practitioners seek new and novel approaches to restore ecosystem functioning in degraded habitats [
1,
2]. Among the numerous aquatic systems that require restoration interventions, temperate estuaries, largely through excessive nutrient input, are the most degraded of all marine ecosystems [
3]. This is particularly true for those estuaries that have been the foci for human expansion, resulting in extensive modifications to catchments for agriculture, industry and other human developments including the expansion of ports, leading to the introductions of non-indigenous species (NIS) and the loss of important biogenic habitats, such as shellfish reefs [
4,
5,
6,
7].
Shellfish reefs, comprising mussels and oysters, were once commonplace in temperate estuaries throughout the world but were largely lost or degraded through a range of anthropogenic influences, including declining water quality, disease, and overharvesting for food and construction material [
4,
5]. As shellfish reefs provide a multitude of ecosystem services, shellfish reef restoration projects have accelerated in the last decade, and this includes the use of artificial hard structures to provide suitable habitat for larval settlement [
4,
5,
6]. For example, in lieu of natural oyster shell, restoration in Chesapeake Bay after 5 years concrete modular reefs deployed subtidally (7 m depth) were found to contain >8000 mussels m
−2 of surface area [
6]. The addition of hard substrates has also been employed in temperate regions of Australia, where large-scale restoration efforts have occurred at over 40 locations since 2015 [
8].
In addition to shoreline erosion protection during storms and their enormous water filtration capacity, shellfish reefs also support fishery production through the provision of habitat and food resources [
9,
10]. This is particularly important in estuaries, which largely comprise unstructured soft sediments, with the protrusion of the three-dimensional hard structure interspersed with interstitial spaces, of variable sizes, providing a predator refuge. Thus, shellfish reefs often support greater fish densities than unstructured habitats [
11]. The three-dimensional structures also promote the settlement, growth and survival of shellfish and other filter feeders by providing better access to feeding opportunities. The elevation provides a buffer zone from the sediment–water interface where oxygen concentrations are typically lower and also reduces exposure to benthic detritus that can become resuspended by waves [
12,
13,
14,
15]. As seston derived from resuspension can contain high concentrations of silt and is thus of low quality for suspension feeders, it can lead to increased processing costs (pseudofeces production) and low absorption efficiencies and thus reduced survival [
7,
16,
17]. For example, the condition of the Horse Mussel
Atrina zelandica was found to decrease after exposure to elevated sediment concentrations after three days and their natural distribution was governed by suspended sediment loads [
17]. The habitat heterogeneity provided by the reefs also promotes biodiversity primarily through the presence of structurally complex microhabitats, comprising matrices of living shellfish, shells and sediment, which support a diverse range of invertebrate species [
18,
19]. For example, beds of the Mediterranean Mussel
Mytilus galloprovincialis in South Africa were found to support a rich assemblage of 80 macroinvertebrate species [
20].
Restoration initiatives are particularly appropriate for estuaries in south-western Australia as these systems have been heavily modified since European settlement [
21,
22]. However, this region is also recognised as a global warming and drying hotspot [
23], where sea surface temperatures have increased at a faster rate (0.1–0.3 °C per decade) than 90% of the ocean globally [
24], and rainfall declines have been among the most pronounced of any Australian region [
25,
26,
27]. The effects of these modifications have been particularly apparent in the microtidal Swan-Canning Estuary, which flows through Western Australia’s capital city, Perth [
28].
The Swan-Canning Estuary is among the most extensively modified estuary in Western Australia and was recognised as one of the most hypereutrophic of 131 estuaries and coastal ecosystems examined globally by Cloern et al. [
29]. The estuary comprises an upper riverine region (the Swan and Canning rivers), a large central basin and an entrance channel, which has been widened and deepened to facilitate Western Australia’s largest and busiest container port. The presence of the port accounts for the estuary containing 46 of the States’ 60 marine NIS (non-indigenous species) [
30]. The basin of the Swan-Canning Estuary also once supported reefs of the native Flat Oyster
Ostrea angasi, as was evident from the banks of shell removed in the first half of the twentieth century [
22]. While the only remaining shellfish in high abundance is
M. galloprovincialis, the larvae of this species has a preference to settle on hard structures, which are now non-existent in deeper water.
The declines in streamflow and increased seawater intrusion have led to salinities in bottom waters throughout the basin of the Swan-Canning Estuary remaining close to seawater (e.g., 32–37) throughout the year [
31]. This is particularly the case during ‘dry years’ when annual rainfall is below average (759 mm), which has occurred in 15 out of the 20 years prior to 2022 [
32]. While the increased marinization of the estuary has resulted in habitat compression of estuarine species upstream, such as the Black Pygmy Mussel
Xenostrobus securis, it has also provided a more stable saline environment for species with marine affinities [
31,
33]. The estuary has also had a long history of hypoxia (<2 mg L
−1 of dissolved oxygen), particularly in the bottom waters of the upper riverine region, causing major impacts for fish and invertebrate fauna [
28,
34]. During ‘wet years’, however, when rainfall is above (or approaches) average, the water column in the estuary basin can become stratified and a hypoxic zone can persist in bottom waters for several weeks at the time [
35].
Large-scale shellfish reef construction projects in Australia have typically involved the deployment of a limestone reef base, which elevates the shellfish off the typically soft and low oxygen bottom waters [
8]. While shellfish reef restoration has largely focused on hatchery-sourced oysters (
Saccostrea glomerata and/or
Ostrea angasi),
M. galloprovincialis, which in Australia comprises both northern and southern hemisphere halotypes [
36], has been receiving increased attention [
7,
8]. Furthermore, as mytilid mussels continue to produce byssal threads throughout their life, if dislodged, they can reattach to new structures.
Mytilids are an important commercial aquaculture species and a range of approaches have been employed to maximise harvests [
37]. These typically involve suspending the mussels in the water column to provide better feeding opportunities and less burial from sediments as well as to reduce mortality from benthic predators, such as crabs and rays [
38]. While the use of longlines or raft cultures, where ropes are suspended from the surface, are more typical to facilitate these requirements, a style of farming developed on the French Atlantic coast, referred to as bouchot (shellfish bed) farming, employs mussels grown on wooden poles staked into the substrate. Although the protrusion of wooden poles from the benthos may not appear natural, these structures, which will decompose over time, have the potential to form the foundations of shellfish reefs; thus, this type of approach is potentially transferable as a restoration option.
The aim of this study was to assess the efficacy of bouchot-style reefs as a restoration option for M. galloprovincialis in the Swan-Canning Estuary. This was undertaken by quantifying the performance of translocated M. galloprovincialis (growth, body condition and survival) on three constructed bouchot reefs, the presence of newly recruited M. galloprovincialis and an assessment of the fish and invertebrate species that colonised these habitats. The results of this study are of interest to environmental managers and will assist shellfish reef practitioners in assessing appropriate approaches and locations for future shellfish reef restoration projects.
2. Materials and Methods
2.1. Study Location
The Swan-Canning Estuary (32.0° S and 115.9° E) in south-western Australia is ~50 km long, has a surface area of ~55 km
2, and remains permanently open to the Indian Ocean (
Figure 1) [
20]. Although most of the estuary is shallow (<2 m in depth), it reaches a maximum depth of ~20 m in the entrance channel, and up to ~17 m in the main basin. The composition of the sediment differs markedly between the shallow marginal shoals (<2 m deep) and the deeper waters (>2 m deep). Nearshore substrates contain a greater proportion of medium-grain-size sand than the offshore sediments, which comprise higher percentages of particulate organic matter and finer inorganic particles (i.e., silt and fine sand) [
39].
South-western Australia experiences a Mediterranean climate [
23], with warm, dry summers and cool, wet winters, and ~80% of the annual rainfall occurs between June and September. This highly seasonal rainfall pattern, combined with the microtidal tidal regime of the Swan-Canning Estuary (tidal range < 2 m), results in marked seasonal variations in water physical–chemical conditions in this salt-wedge estuary [
33]. Salinities are typically stable and relatively high throughout much of the estuary basin during the austral summer (December to February), but during winter, may vary markedly along the estuary following substantial freshwater discharge [
33], leading to marked stratification of the water column and hypoxia [
28,
33].
The Swan-Canning Estuary and its catchment have been heavily modified since European settlement and contains the City of Perth, the capital of Western Australia (~2 million people). Such modifications include substantial land clearing for urbanization and agriculture, hardening of shorelines, blasting and dredging of the mouth of the estuary to permanently connect it to the ocean and allow the construction of Fremantle Port [
21]. As such, the system was categorised as “extensively modified” in an audit of the condition of ~1000 estuaries across Australia [
40].
2.2. Environmental Variables
Long-term average monthly rainfall (1944–2023) and total monthly rainfall (2019–2023) recorded at Perth Airport (Station 009021) were extracted from the Australian Government Bureau of Meteorology (2023) database [
31]. Long-term average monthly streamflow (1970–2023) and total monthly streamflow (2019 and 2023) recorded for the Swan River (Station 616011) were extracted from the Water Information Reporting (WIR) database [
41]. Weekly values (2019–2023) for water temperature, salinity and dissolved oxygen concentration at the surface (0.5 m depth) and bottom of the water column (6.0 m depth) recorded at Heathcote (site 6161869;
Figure 1) were also extracted from the WIR database [
41].
2.3. Reef Construction and Sampling Regime
Individuals of M. galloprovincialis (shell length 40–70 mm, 0+ age class) were collected from mooring ropes in the entrance channel of the Swan-Canning Estuary, and their external epibionts were removed by a combination of aerial exposure and manual agitation. Three liters of M. galloprovincialis was then used to seed the top 1 m of each of 300 wooden (jarrah; Eucalyptus marginata) stakes (1500 (H) × 25 (W) × 25 (D) mm). For this, a PCV pipe was used to apply a commercial cotton sock around the mussels on each stake. Thus, the use of established mussels negated the use of coconut fibre ropes to collet spat, as is used in traditional bouchot culture. Instead, this study employed a similar process to that used in commercial aquaculture where established mussels are reattached to grow-out ropes at a desired density that promotes growth, and this has been highly successful for food production.
The collection and socking of mussels onto the stakes were undertaken over three consecutive weekends and involved community volunteers. On each weekend between late March and early April 2019, 100 stakes were deployed at 3 locations in the main basin of the Swan-Canning Estuary (
Figure 1). The reefs were constructed in ~4 m depth by SCUBA divers hammering the stakes ~0.5 m into the substrate approximately 1 m apart in a grid pattern covering an area of ~100 m
2. These sites were selected based on their depth, with 3 m being the minimum vertical clearance allowed for structures in navigable waters, and their requirement to be located outside of recreational zones and marine parks.
Three stakes were retained at the time of deployment and a further three stakes were removed each month from each of the three reefs by snorkelling until February 2020, with a further three stakes collected from each reef in March 2021 and January 2023. Stakes were individually placed in large plastic bags, which were sealed and transported to the laboratory, where the above substrate portion of each stake was divided into three ~300 mm sections, with the top of the stake between 700 and 1000 mm from the sediment, the middle section 350 to 650 mm from the sediment and the bottom section of the stake 0 to 300 mm from the sediment. For each stake and section, each M. galloprovincialis was counted, its shell length measured to the nearest 0.5 mm using vernier callipers and its tissue removed and dried at 60 °C for 72 h or until a constant dry weight was reached. These samples were employed to explore patterns in growth, recruitment, body condition and survival.
Fish assemblages at the reef sites and adjacent unstructured control sites, located 0.5–1 km away (
Figure 1) were monitored using Remote Underwater Videos (RUVs) in accordance with Animal Ethics Permit O2996/17. RUVs were chosen over baited systems as the spread of the bait plume has the potential to attract fish from outside of the reef area and, thus, those fish that are not using the reef directly. Each RUV comprised a camera (GoPro Hero 2018; GoPro, Inc. San Mateo, CA, USA) attached to a weighted rectangle steel frame (340 mm in height and 310 × 310 mm in length and width) and lowered to the substrate. All footage was recorded at a resolution of 1920 × 1080 pixels at 60 frames s
−1 and for at least one hour.
In each month between October 2018 and June 2020 (except for January, March and October 2019), one RUV was deployed at each reef and adjacent control site. On each sampling occasion, one of the three reefs was chosen at random and a second deployment was conducted. RUV videos were collected prior to the deployment of the reefs (October 2018 to February 2019) to determine whether there were any underlying differences in the characteristics of the fish fauna among site and then post-deployment (April 2019 to June 2020) to investigate differences between the habitats.
The relative abundance of each fish and large crustacean species, e.g., Blue Swimmer Crabs (
Portunus armatus), which is a key recreationally targeted species [
42], was determined by counting the maximum number of individuals (MaxN) of each species in each frame commencing from the moment the RUV frame hit the benthos and continuing for 60 min [
43]. Each of the taxa recorded was allocated to a number of functional guilds, i.e., estuarine usage, feeding mode and habitat (
Table A1), following the rationale given by [
44,
45].
The invertebrates directly attached to the stakes and their associated epibiont community were removed from the top section of each stake (where they were most abundant), from samples in May, June, September and December 2019, and March 2021 and January 2023. Each sample was wet sieved using a 2000 μm mesh, and the invertebrates were stored in a 70% ethanol solution. Invertebrates were then separated from any material and sediment retained from the sieving and subsequently enumerated and identified to their lowest taxonomic level.
2.4. Statistical Analyses
The growth patterns of M. galloprovincialis were initially explored using traditional and seasonally adjusted von Bertalanffy growth models. While the seasonally adjusted provided a better fit (p < 0.001), the model still yielded low R2 values when applied to each reef, i.e., R2 = 0.18–0.73. Instead, as all translocated M. galloprovincialis were assumed to be the same age on each sampling occasion, shell length (mm) differences among reefs and stake sections in each month were analysed using Analysis of Variance (ANOVA). The assumption that all mussels belonged to the same cohort was consistent with the shell lengths in the adjoining embayment of Cockburn Sound where this species is grown commercially and spat collection ropes are replaced each year, the location of the mussels towards the top of the mooring ropes and, compared to the larger mussels that were present, largely lacked biofouling.
The body condition of
M. galloprovincialis was explored using the relationship between soft tissue dry mass and shell length, an approach used in previous assessment of bivalve condition (e.g., [
46,
47]). To facilitate comparisons among month, reef site and stake section Analysis of Covariance (ANCOVA) was employed, which standardised for shell length (to the mean value 61 mm) to help account for spatial and temporal differences in shell length. Prior to analyses, data for shell length and dry mass were log transformed. Shell length and body condition data were analysed using IBM SPSS (Version 28.0).
The MaxN of each species in each replicate RUV sample was used to construct a data matrix and calculate the number of species and total MaxN in each sample. The data for these two univariate variables were examined using a Draftsman plot, to visually assess if the distribution of values were skewed and, if so, the transformation required to ameliorate this effect. This showed that total MaxN required a fourth-root transformation. The data for both variables before and after the deployment of the bouchot-style reefs were used to construct four Euclidean distance matrices. In turn, each was subjected to a two-way univariate Permutational Multivariate Analysis of Variance (PERMANOVA) [
48] test to determine whether they varied significantly (
p < 0.05) with Habitat (two levels; reef and control) and Time (two levels; pre-deployment, i.e., October-December 2018 and February 2019 or post-deployment, i.e., the 14 months between April 2019 and June 2020).
The data matrix of the MaxN abundance of each species in each replicate sample was divided to separate those samples collected at each site before and after the deployment of the bouchot-style reefs. The data in each matrix were dispersion-weighted [
49] to down-weight the abundances of heavily-schooling species whose numbers were erratic, relative to those that return more consistent values, and then square-root transformed to balance the contributions of relatively abundant species, compared to those with lower MaxN values [
50]. The transformed data in each matrix were then used to construct a Bray–Curtis resemblance matrix. These were separately subjected to the same two-way PERMANOVA design above. If required, post hoc testing was conducted using pairwise PERMANOVA.
Trends in the fish faunal community data were explored using a range of multi-dimensional scaling (MDS) approaches and Canonical Analysis of Principal coordinates (CAP) [
50]. The averages of repeated bootstrap samples (bootstrapped averages) were used to construct a mMDS ordination plots for the main effects of Habitat and, if significantly different, also Time. Superimposed on to each plot were points representing the group average (i.e., the average of the bootstrapped averages) and the associated, smoothed and marginally bias-corrected 95% bootstrap region, in which 95% of the bootstrapped averages fall [
51]. Note that the Time factor for the post-deployment analyses contained 14 levels, which, for clarity, was too complex to show on a bootstrapped mMDS plot and so was visualised using a centroid nMDS ordination. The same Bray–Curtis resemblance matrix was used to create a distance among centroids matrix, which was, in turn, used to generate a nMDS ordination plot [
52]. A CAP ordination plot was created to show how the fish faunal communities changed among habitats before and after deployment. Superimposed onto the CAP are vectors for species whose abundance changes in a linear direction (Pearson correlation > 0.5) relative to the CAP axes.
Shade plots [
53] were used to illustrate trends in the transformed abundance of each species across Habitats and Time separately, when significant, for the pre- and post-deployment data. A shade plot is a visualization of the averaged data matrix, where a white space for a species demonstrates that it was not recorded, while the depth and colour of shading, ranging from light grey shades through to black, represents increasing values for the abundance of that species. Species were seriated based on their Bray–Curtis similarities and placed in optimum order.
Data on the invertebrate community were analysed using the same suite of analyses as for the fish communities, although in the PERMANOVA tests the Habitat factor was replaced by Site (three levels; Reef 1, 2 and 3) and, for visual clarity, the CAP analyses only considered Time rather than the significant but far less influential Time × Site interaction. All statistical analyses for fish and invertebrate communities were performed using the PRIMER v7 multivariate statistics software package [
51] with the PERMANOVA+ add-on [
50].
4. Discussion
This study was the first to assess the efficacy of employing a bouchot-style approach as an option for shellfish reef restoration. This was assessed at three locations in the Swan-Canning Estuary through quantitative assessment of the biological characteristics (growth, body condition and survival) of translocated M. galloprovincialis, the presence of naturally recruited M. galloprovincialis and the suitability of the reef structure for fish and other macroinvertebrates. The study demonstrated that distance from sediment was an important factor when considering shellfish reef restoration approaches, with mussel survival, growth, body condition and recruitment greater on the top section of the stake. Furthermore, more fish species inhabited the bouchot-style reefs than unstructured habitats and the reefs were also colonised by a range of invertebrate species, including 11 native and 6 non-indigenous species; thus, this type of an approach has potential as a shellfish reef restoration option. However, due to bottom water hypoxia that can occur in wet years or following unseasonal rainfall, the study also found that the deeper waters of the estuary may not be conducive to sustaining populations of M. galloprovincialis over multiple years, the results of which can provide guidance for future shellfish restoration initiatives.
4.1. Overall Performance of Mytilus galloprovincialis
The performance (recruitment, growth, body condition and survival) of M. galloprovincialis differed among reef locations and particularly on stake section, i.e., distance from the substrate. Stake top sections, which were 700–1000 mm above the sediment, consistently outperformed those on the middle and bottom sections, which were 350–650 and 0–300 mm above the sediment, respectively. While the factors driving these marked differences among stake sections may be widely varying, the consistency in results implies that distance from substrate plays an important role, and this may be related to biotic and abiotic factors.
While the microtidal Swan-Canning Estuary contains weak tidal currents, its large width-to-depth ratio allows wind-generated waves to resuspend the fine sediments that are characteristic of estuary basins [
7,
13]. Although sediments can be displaced and distributed by turbulence throughout the water column, turbidity typically reaches its maximum near the seabed [
54]. Thus, in close proximity to the seabed, particulate inorganic matter (PIM) concentrations are likely to exceed a threshold that mussels stop feeding, due to the metabolic cost of pseudofeces production and/or their sorting capacity becoming saturated [
55,
56,
57,
58]. It is thus relevant that in a recent study ~5 km upstream of Reef 1, where the performance of
M. galloprovincialis was consistently poor, PIM in waters close to the seabed reached 305 mg L
−1 (A. Cottingham unpublished data) and thus far greater than normal mussel feeding conditions (PIM < 10 mg L
−1) [
7,
59,
60]. Furthermore, these conditions are also more likely to be accompanied by lower dissolved oxygen concentrations [
14].
4.2. Body Condition of Translocated Mytilus galloprovincialis
In most cases, the body condition of
M. galloprovincialis decreased following translocation, particularly at Reef 1. This indicates that mussels are using energy reserves to survive, which is unsustainable over long durations [
61,
62]. This conclusion is consistent with recent
M. galloprovincialis feeding behaviour experiments using water collected 0.5 m above the substrate in the basin of the Swan-Canning Estuary during summer, which yielded absorption efficiencies that were likely to equate to a net energy loss and thus progressive declining mass if conditions persist [
7]. The low values for body condition during winter may have also be exacerbated by spawning [
63,
64] and these have been observed in other marine bivalves, such as the New Zealand green lipped mussel
Perna canaliculus [
65], and this may have further reduced glycogen reserves leading to increased vulnerability to stressful conditions [
66]. The conclusion that spawning occurs in winter in the Swan-Canning Estuary was consistent with the presence of new recruits in samples in August.
4.3. Factors Influencing Reef Sustainability
Despite the lifespan of
M. galloprovincialis typically ranging between two and four years, only 2% of translocated mussels survived until the end of the study, i.e., 10 months [
67,
68]. This rate was far lower than that recorded in a study in the Aegean Sea in which 77.5% of
M. galloprovincialis in nets survived for a similar time [
69]. While the low survival rate of translocated mussels is likely to be at least partly explained by sequential declines in body condition following translocation, only four mussels were recorded in all samples in March 2021, and those include those naturally recruited onto the reefs. The scarcity of
M. galloprovincialis was accompanied by marked declines in other invertebrate species, which declined from 4495 individuals belonging to 14 species in December 2019 to 35 individuals belonging to 4 species in March 2021. Such consistency in declines in the number of species and individuals is important when considering potential drivers.
As the majority of species (both NIS and native) recorded on the bouchot-style reefs were temperate species, summer temperature may have exceeded critical thermal maximums for some of these. For example, the abundant native species
A. variegatus has a similar geographic distribution to
M. galloprovincialis and thus also likely has a similar temperature tolerance. It is thus relevant that the water temperature in the Swan-Canning Estuary was >24 °C for several months of the year, a temperature that has been associated with mass mortality of
M. galloprovincialis [
70]. While some of the abundant NIS species were introduced from cooler waters in Europe, e.g.,
A. aspersa, and thus also have a limited tolerance to elevated temperatures (26 °C) [
71],
A. senhousia has an upper thermal limit of 31 °C and is capable of surviving in salinities between 17 and 37 [
72,
73]. Thus, the temperatures and salinities in the current study remained well within the tolerance range for
A. senhousia.
While mytilids generally have a better survival when exposed to low oxygen than most other invertebrate species, this tolerance has evolved from oversized lamellibranchiate gills that aid in filter feeding and the ability of mussels to survive during emersion over a tidal cycle [
74,
75]. For example, the blue mussel
Mytilus edulis only showed a slight reduction in respiration rates when subjected to dissolved oxygen concentrations decreasing from 9 to 2 mg L
−1 and demonstrated excellent survival (98%) when exposed to 2 mg L
−1 for 16 days [
75]. Mytilids also contain a robust mechanism to cope with these stressors, which involves isolating their living tissue through valve closure [
7,
76]. For example, the valves of
M. galloprovincialis have been shown to close in response to declining salinities (~10) and hypoxic conditions [
75,
76]. During valve closure, however, oxygen concentrations in the mantle cavity can decline rapidly, resulting in a switch to anaerobic respiration [
77]. Although this metabolic pathway still yields adequate amounts of ATP and glycogen stores [
78,
79], prolonged closure can result in substantial mortality [
76]. For example, exposing
M. galloprovincialis to a salinity of 5, thus inducing valve closure, resulted in 77% mortality after ~4 days [
76]. It was thus relevant that dissolved oxygen concentrations during 2021 dropped below 2 mg L
−1 for several weeks, conditions that are also considered lethal for most benthic invertebrates [
80,
81].
4.4. Faunal Composition
While the average number of species and total MaxN remained relatively constant throughout the four months prior to the reef deployment, a discernible temporal shift in faunal composition occurred following the first sampling occasion in October 2018, i.e., mid-spring, in which the fish community was depauperate. This is consistent with the well-established seasonal changes in the fish fauna of permanently-open temperate estuaries [
82,
83,
84,
85]. During summer and autumn months, when salinities and temperatures are elevated, the number of fish species with marine affinities are typically far greater than in winter, when salinities and temperatures decline [
86,
87]. Furthermore, it is relevant that the annual rainfall during that year (2018) of 744 mm would have produced water column stratification and reduced oxygen condition in bottom waters.
Following reef deployment, fish faunal composition differed significantly between the reef and control sites, with these differences largely attributed to the increased numbers of species utilising the reef habitat. Comparisons between structurally complex shellfish reef sites vs. unvegetated controls undertaken in other parts of the world have yielded similar results, which are explained by the provision of shelter and additional prey sources [
10,
88,
89,
90,
91]. Almost all species found exclusively at reef sites, namely the boxfish
A. amygdaloides, the seahorse
H. angustus and several species of leatherjacket, e.g.,
Monacanthus chinensis,
Acanthaluteres vittiger and
Eubalichthys mosaicus, were marine species with a preference for complex structures [
43,
92,
93,
94]. While the number of benthopelagic species remained similar between reef and control habitat, there was a marked increase in abundance at the reef sites. This difference was largely driven by the apogonid,
O. ruppelli, which formed large schools at the reef sites, capitalising on the increased structure and availability of its small crustacean prey [
95].
The number of species, total MaxN and faunal composition of invertebrates all differed significantly between months, with successive increases in mean species richness reflecting colonisation of the reefs over time, and these were assisted by the provision of new surfaces as the numbers of
M. galloprovincialis declined. Of the invertebrate species, those of barnacles and ascidians were the most abundant. The ascidian
Ciona intestinalis, which is known to feed on bivalve larvae [
96], was one of the most abundant species during winter (May–September 2019), which coincided with the spawning season of
M. galloprovincialis. This ascidian then declined markedly in December, which may be attributed to the inability of this species to function at temperatures >21 °C [
97]. Although this may not directly contribute to its mortality, it would have reduced its ability to compete, and this was evident as the presence of the colonial ascidian
S. plicata increased during the warmer months. Although
S. plicata was not abundant, its individuals were typically large. This species is also more tolerant to warmer waters and prefers higher salinities, which could explain its appearance in December and not in the cooler months when
C. intestinalis was more abundant [
98].
There was a marked shift in invertebrate composition in March 2021 following the prolonged period of hypoxia with species such as
M. galloprovincialis, the barnacle
A. variegatus and the ascidians
C. intestinalis, A. aspersa and
S. pilcata all disappearing and being replaced by sponges, which were only recorded in one of the three reefs and even then, in low quantities prior to hypoxia. Ascidians are considered fairly sensitive to low oxygen conditions, with several field studies demonstrating that
C. intestinalis suffers mortality while oysters survived [
96,
99]. While
M. galloprovincialis and a congener of
A. variegatus have been shown to have considerable tolerance to hypoxia, they have a lower threshold than those of sponges [
100], thus explaining why the latter taxa dominated at all reef sites in March 2021. However, their low abundance before March 2021 and in January 2023 suggests they are poor competitors, and their occurrence was likely linked to their ability to withstand prolonged hypoxia.
Although not contributing substantially to changes in faunal composition, it is important to note the impact of the non-native solitary ascidian
A. aspersa and barnacle
A. variegatus, which colonised the stakes rapidly upon their deployment. Both these species were recorded in high abundances throughout months and on all reefs. The rapid colonisation of the stakes particularly by
A. aspersa may be attributed to its biology. This species can live for 18 months and adapt their reproductive biology from gonochoristic to a hermaphrodite to allow spawning to occur throughout the year [
101]. Furthermore,
A. aspersa is well suited to urbanised and highly eutrophic estuaries and deeper slow-moving waters [
102], which are characteristic of the basin region of the Swan-Canning Estuary. These biological traits likely explain the lack of temporal change in the abundance of this NIS ascidian.
4.5. Contribution to Future Knowledge
Although M. galloprovincialis is found throughout the Swan-Canning Estuary, it is almost exclusively restricted to a ~1 m depth and is particularly abundant on hard substrates in the intertidal zone, the majority of which are man-made (e.g., jetties and pylons). The success of mytilids in the intertidal zone, compared to many other types of fauna (and flora), is largely attributed to the ability of mytilids to prevent desiccation by closing their valves during emersion. Thus, shellfish reef restoration using these types of species would be far better suited to nearshore shallow/intertidal zones.
In contrast, as supported by this study, subtidal hard structures can attract competition from a range of faunal groups, some of which were able to attain abundances far greater than
M. galloprovincialis. While the majority of the first colonisers were native species, by December (after about six months), this trend reversed and NIS species dominated, which included
A. aspersa, the tunicate
Molgulidae sp. and
A. senhousia. Thus, potential locations for shellfish reef restoration efforts should consider the types of species that may colonise these habitats and whether they are likely to outcompete the target species. Moreover, the timing of the deployment should, where possible, not coincide with the availability of larvae from NIS, which could settle on the newly available habitat. In the case of the Swan-Canning Estuary, and despite its environmental, cultural, and societal importance and its historic loss of shellfish reefs, the estuary contains 46 of the States’ 60 marine NIS [
30], and deployment of hard structures in the vicinity of busy ports thus needs to be given consideration. While seeding densities, size at seeding and timing of reef constructions could potentially result in improved
M. galloprovincialis performance, the ability of mussels (and other fauna) to be sustained over longer durations appears limited, particularly during wet years in bottom water, which can become hypoxic for long durations.