Authorized (original) (raw)

Long-term Recovery - Review of sediment condition at Marine Farm lease No. 76 (Gunpowder Jetty), North West Bay

Macleod, C. K., Eriksen, R.S., Davey, A., Kelly, B. and Ross, D. J.
img-0.jpeg

IMAS - in collaboration with the Tasmanian State Government
img-1.jpeg
translatingnatureintoknowledge
www.imas.utas.edu.au

Title

Long-term Recovery -Review of sediment condition at Marine Farm lease No. 76 (Gunpowder Jetty), North West Bay

Authors

Catriona K. Macleod, Ruth S. Eriksen, Adam Davey, Bronagh M. Kelly and D. Jeff Ross

Disclaimer

The authors do not warrant that the information in this document is free from errors or omissions. The authors do not accept any form of liability, be it contractual, tortious, or otherwise, for the contents of this document or for any consequences arising from its use or any reliance placed upon it. The information, opinions and advice contained in this document may not relate, or be relevant, to a reader’s particular circumstance. Opinions expressed by the authors are the individual opinions expressed by those persons and are not necessarily those of the Institute for Marine and Antarctic Studies or the University of Tasmania.

This work is copyright. Except as permitted under the Copyright Act 1968 (Cth), no part of this publication may be reproduced by any process, electronic or otherwise, without the specific written permission of the copyright owners. Information may not be stored electronically in any form whatsoever without such permission.
©The Institute for Marine & Antarctic Studies (IMAS)2013.

Enquires should be directed to:

Dr Catriona Macleod
Institute for Marine and Antarctic Studies
University of Tasmania
Private Bag 49, Hobart, Tasmania 7001, Australia
Catriona.Macleod@utas.edu.au
Ph: (03) 62277277 Fax: (03) 62278035

Citing of this report

Macleod, C. K., Eriksen, R.S., Davey, A., Kelly, B.and Ross, D. J. (2014).Long-term Recovery - Review of sediment condition at Marine Farm lease No. 76
(Gunpowder Jetty), North West Bay. IMAS Technical Report 32pp.

Contents

Contents … 3
Non-technical summary … 4
Acknowledgements … 5

  1. Introduction … 6
  2. Methods … 7

2.1 Sample Site Location … 8
2.2 Sediment Characterisation(Physical-Chemistry Properties) … 9
2.3 Faunal Characterisation … 11
2.4 Video Footage … 11
2.5 Mesocosm Incubations … 11
2.6 Statistical Analysis … 12
3. Results and Discussion … 13
3.1 Particle Size Distribution … 13
3.2 Organic Matter Measurement … 14
3.3 Sulphide … 14
3.4 Macrofaunal Diversity … 15
3.5 Nutrient Flux Measures … 16
3.6 Macrofauna … 19
3.7 Video Assessment … 22
3.8 Copper and Zinc loadings … 23
4. Conclusions … 25
5. References … 26

Non-technical summary

This study was undertaken to assess the Gunpowder Jetty site(Marine Farm lease No.76) 10 years after it was finally vacated to determine whether the benthic, visual and physicalchemical conditions were consistent with background conditions and to what extent the system has recovered. The results suggest that, based on the selected environmental condition variables assessed, the sediments around the Gunpowder lease have largely recovered. Whilst some parameters still showed small differences between locations these do not appear to be having a major effect on the ecology or nutrient mineralisation processes and consequently there is no evidence that the farming activity has had any permanent impact on the ecology in this area.

KEYWORDS: aquaculture impacts, environmental assessment, benthic infaunal assessment, sediment biogeochemistry, nutrients dynamics

Acknowledgements

We would like to thank the staff at DPIPWE for their help in compiling the data, in analysing the video data and in providing feedback and suggestions along the way (particularly Graham Woods and Eric Brain) and Marc Santos and John Preston for assistance with field sampling.

1. Introduction

Management of the environmental impact of finfish farms requires a good understanding of the local physical, chemical and biological processes, and a reliable understanding of enrichment and recovery status. Localised aquaculture impacts have been well studied (e.g. Brown et al., 1987; Ritz et al., 1989; Weston, 1990; Woodward et al., 1992; Holmer and Kristensen, 1992; Findlay et al.,1995; Cheshire et al., 1996; Hargrave, et al., 1997 . Buschmann et al., 2006; Hargrave et al., 1997; Janowicz and Ross, 2001; Karakassis et al., 2002). Organic enrichment is the main effect of cage fish farming and areas of localised organic enrichment can arise when fish feed and faeces settle on the seabed (Brown et al., 1987; Lumb, 1989; Holmer and Kristensen, 1992; Hargrave et al., 1993; Karakassis et al. 2000). The ability of sediments to assimilate such increases in organic load is variable and depends on both abiotic and biotic factors (Lumb, 1989; Chang and Thonney, 1992; Wu et al., 1994). Typically, the flora and fauna of impacted sediments adapt to utilise this new nutrient source, resulting in changes in benthic community structure. However, if the sediment’s capacity to assimilate organic inputs is exceeded and the sediment becomes anoxic, the sediment biogeochemistry will be altered towards a system dominated by anaerobic forms of metabolism and toxic degradation products (hydrogen sulphide and ammonia) can be released into the environment affecting farm production and aquatic ecosystem health (Gowen and Bradbury, 1987, Gowen et al., 1988; Rosenthal and Rangeley, 1988; Holmer and Kristensen, 1992, Hargrave et., 2008). As this organic matter degrades the sediment oxygen demand increases and the underlying sediments become deoxygenated (Hargrave et al., 2008). This results in a lowered redox (oxidation-reduction) potential and increased sediment sulphide concentrations (Hargrave et al., 2008). Sediments under cages can also be enriched with zinc and copper sourced from feed or antifouling agents. These conditions will affect the prevailing benthic communities, changing the species mix and in the worst case scenarios where the sediments are anoxic, azoic conditions can result (Gowen and Bradbury, 1987).

Whilst it is important to have a good understanding of the extent to which the sediments are impacted, it is also important to understand the potential for recovery, as this determines how long an impact will persist and the potential for reversibility of any effects. To date estimates of recovery vary greatly, from less than six months (Brooks et al., 2003; Ritz et al., 1989) to more than 5 years (Brooks et al., 2004). This variability in the estimates of recovery probably reflects differences in both underlying ecosystem conditions and farm management but may also be related to the specific recovery end-point employed, as the benchmark against which recovery is being judged can have a significant impact on recovery time estimations (Macleod et al., 2004b, Keeley et al., 2013). That said, most studies tend to target short-term recovery, with only a very few addressing longer-term remediation. There are still relatively few studies which have been conducted over long enough timeframes to

identify “complete” recovery i.e. to pre-impact condition or to conditions consistent with those of selected control sites (e.g. Karakassis et al., 1999; Macleod et al., 2004b).

Against this background, a study was initiated in 1999 jointly by the aquaculture industry (Aquatas Tasmania Australia Pty Ltd) and the Department of Primary Industries, Water and Environment (DPIWE) with the aim of assessing the rate of sediment recovery after removal of all salmon cages from a lease located in North West Bay, South-East Tasmania. It is clear that there is a very broad range of factors which may influence the sediment recovery rate, but one of the principal objectives of this research was to identify the changes/stages in both the benthic faunal community structure and the physical/ chemical status of the sediments over time associated with long term fallowing of an intensively farmed marine finfish cage site in the cool temperate waters of Tasmania. The results of this research were published in Macleod et al (2004) and the key findings suggesting that sediments under the cage were highly impacted at the time that the cages were removed but that the extent of impact diminished with both time and distance from the cages. However, the results suggested that although most of the biogeochemical parameters assessed had returned to conditions equivalent to those at the reference after 2 years, the benthic community structure still differed significantly.

This project follows on from that original study to evaluate the benthic, visual and physicalchemical condition of the sediments at the Gunpowder lease almost 13 years after it was first assessed, and 10 years after it was finally vacated. These results provide a valuable insight into the rate and extent to which the sediment has recovered and inform our understanding of ecosystem resilience, providing valuable information for long-term management of aquaculture lease areas to ensure the sustainability of ongoing operations.

Study Aims:

2. Methods

Sampling was undertaken at the same sites and in a manner consistent with the protocols described in the original study (Macleod et al., 2002). In the initial study all samples were collected by diver however, more recent safety precautions now prohibit diving under these conditions and so the sediment samples in the current study were all collected by core or grab, with the video footage collected by ROV.

Included in this temporal comparison are the data collected one week after removal of cages, the data from 24 months after removal and the most recent data collected almost 13 years after initial cage removal (which is approximately 10 years after complete vacation of the farm in July 2003).

2.1 Sample Site Location

The Gunpowder lease (lease 76), in North-West Bay, Tasmania, was first granted to Aquatas Pty Ltd in 1985 and operated as a commercial salmon farming operation until August 1999. When in operation, the farm was a smolt ongrowing site, occupying an area of 3.12 ha , with up to 16 production pens including both 60 and 80 m circumference polar circles. The lease area is relatively shallow ( 14−20 m14-20 \mathrm{~m} ) and can be subject to elevated summer temperatures (DPIF, 1997). In the year prior to closure the farm stocked approximately 200 - 300 tonnes of fish; however stocking levels were markedly reduced in the 3-4 months prior to the site’s closure as stock were transferred from the site (S.Percival (Aquatas), pers. comm.). Over the preceding four years this site had been stocked more or less continually with little or no fallowing, adjustment of stocked biomass according to the time of year (water temperatures) was the only response to environmental conditions (S.Percival (Aquatas), pers. comm.). Relatively good fish performance was achieved at the site and farm management considered the growing conditions of this lease to be “average”. The site was restocked briefly after the initial study, from mid-October 2002 to early-April 2003, and was permanently vacated shortly thereafter.

Three cage positions were randomly selected prior to the removal of the cages, in late August 1999 for inclusion in the initial study, and fixed transects were positioned on the seabed running from directly beneath the cage ( -10 m ) to a distance 35 m from the cage boundary (Fig. 1, Table 1). Reference stations were established 150 m from the cage position. The positions of all stations were fixed using a differential global positioning system (DGPS).In this follow-up study only the -10 m and 150 m stations were re-assessed.

In the final recovery assessment samples were also collected from reference and compliance positions originally sampled as part of the DPIPWE original baseline survey at this farm site. Consequently, the full suite of samples surveyed includes; 6 sites consistent with the original University of Tasmania study (i.e. the -10 m and 150 m for each of the original transect positions) and 6 sites collected from DPIPWE’s original baseline survey (Fig. 1). In combination these sites reflect 4 previously farmed sites (1-4), 4 boundary/ compliance sites (5-8) and 4 unfarmed reference sites (9-12); the location of the sites is shown in Fig. 1 and their position information and identification in each of the previous studies is shown in Table 1.

Table 1. Site identification details - showing current and previous labels, eastings and northings and relevant impact category.

Site ID Original ID Easting/ Northing Category
1 MFB 1 524242/5234269524242 / 5234269 Farm
2 TAFI T1 (-10m) 524250/5234204524250 / 5234204 Farm
3 TAFI T2 (-10m) 524299/5234159524299 / 5234159 Farm
4 TAFI T3 (-10m) 524348/5234111524348 / 5234111 Farm
5 MFB 2 524208/5234349524208 / 5234349 Compliance
6 MFB 6 524214/5234129524214 / 5234129 Compliance
7 MFB 3 524396/5234312524396 / 5234312 Compliance
8 MFB 5 524549/5234032524549 / 5234032 Compliance
9 MFB 7 523918/5234687523918 / 5234687 Reference
10 TAFI T1 (150m) 524139/5234105524139 / 5234105 Reference
11 TAFI T2 (150m) 524199/5234051524199 / 5234051 Reference
12 TAFI T3 (150m) 524260/5233990524260 / 5233990 Reference

img-2.jpeg

Figure 1. Location of sampling sites in relation to the original Gunpowder Jetty lease. (Reference - O, Compliance - O, Farm - O).

2.2 Sediment Characterisation(Physical-Chemistry Properties)

At each sample station three replicate core samples were collected using perspex tubes 250 mm long and 45 mm internal diameter. Samples were collected for analysis of particle size composition, organic content, loss on ignition, (LOI), redox potential, sulphide concentration. Redox potential and sulphide concentration were measured at two points ( 1 cm and 4 cm ) in a single core from each triplicate.

Half of the top 4 cm from two cores was collected for sediment particle size analysis. A subsample of each was passed wet through a graded series of sieves ( 4 mm,2 mm,1 mm4 \mathrm{~mm}, 2 \mathrm{~mm}, 1 \mathrm{~mm},

500μ m,250μ m,125μ m500 \mu \mathrm{~m}, 250 \mu \mathrm{~m}, 125 \mu \mathrm{~m} and 63μ m63 \mu \mathrm{~m} ). The sediment retained on each sieve was dried and weighed and the percentage of the total sample weight calculated. The fraction smaller than 63μ m63 \mu \mathrm{~m} was determined by calculation of the difference between the initial sample weight and the combined weight of the retained fractions.

Total organic matter was determined by a modification of the loss on ignition technique (Greiser and Faubel, 1988). Samples collected from the top 4 cm of each core were homogenised and a sub-sample of approximately 2-5 grams taken. In order to remove excess carbonate from the samples, samples were sieved to remove large shell fragments and any remaining carbonate was neutralised by acidification with 1 M HCl . The samples were oven dried for 24 hours at 60∘C60^{\circ} \mathrm{C} before being transferred to a muffle furnace for 4 hours at 500∘C500^{\circ} \mathrm{C}. The weight of organic material was calculated as the difference between the oven dried and final furnace ashed weights.

In the laboratory sub-samples ( 4 ml ) were taken from each replicate core for sulphide analysis at depths of 1 cm and 4 cm using a cut-off 5 ml syringe. Sulphide was measured using a Cole-Parmer 27502-40 silver/sulfide electrode following the technique described by Wildish et al. (1999). Sulphide standards were prepared before each sampling event to establish electrode calibration curves.

Redox potential data was not analysed in this study as the original data from 1999 and 2002 was considered doubtful due to contamination of the redox probe during the study, and it would appear that a similar issue may have occurred with the 2012 data.

Sediment samples for metal (copper and zinc) analysis were collected using a multicore sampling device. Upon retrieval the cores were sealed using latex bungs. Overlying water was carefully removed, and the top 30 mm of sediment extruded using a plunger; the sediment samples from all 3 mini cores collected at each position were combined and homogenised in a 120 ml jar. Samples were refrigerated prior to analysis for total recoverable copper by CSIRO.

Samples of porewaters for evaluation of copper and zinc loading were collected from the bulk sediment samples. Duplicate tubes were centrifuged at 3500 rpm for 15 minutes, at 4∘4^{\circ} C. The supernatant from each of the duplicate sub-samples was combined to give sufficient volume for analysis, and filtered through 0.45μ m0.45 \mu \mathrm{~m} filters into 10 mL acid washed vials. This process was undertaken in a nitrogen filled glove box to avoid introducing additional oxygen to the samples. Samples were acidified with Merck Suprapure Nitric Acid to 0.2%0.2 \% v/v. A procedural blank (deionised water) was prepared with every 10 samples. Samples were sent to CSIRO for analysis. Where results are reported as below detection, the limit of reporting has been halved for the purpose of plotting and statistical analysis.

2.3 Faunal Characterisation

In the initial study five replicate samples were collected by diver for assessment of the benthic macrofaunal community structure at each sample sites using hand held 150 mm diameter PVC pipe corers to a depth of 100 mm (sampling area of 0.0177 m20.0177 \mathrm{~m}^{2} ). Samples were placed immediately into mesh bags (0.875 mm)(0.875 \mathrm{~mm}), and then rinsed on the boat prior to being transferred to containers of 4%4 \% formalin for fixation of the fauna. In the laboratory each sample was sieved to 1 mm , sorted and the animals retained were identified to the lowest possible taxonomic level and enumerated. In 2012 the same processing and identification protocols were followed but only 3 replicate samples were collected at each sample position using a Van Veen grab (sampling area of 0.0675 m20.0675 \mathrm{~m}^{2} ).

2.4 Video Footage

Camera drops/ ROV footage was collected at each location in 1999 and 2012 and analysed by DPIPWE according to the criteria outlined in Marine Farm Environmental Monitoring protocols and described by Crawford et al. (2001).

2.5 Mesocosm Incubations

Three replicate large mesococm core samples were collected for incubations from each of the original TAFI sites at the -10 m and 150 m positions (sites 2,3,4,10,112,3,4,10,11 and 12), a total of 18 cores. Incubations were profiled for oxygen and nutrient flux rates assessed (start/ end).

Sediments for mesocosms were collected by box corer. A mesocosm core ( 300 mm×150300 \mathrm{~mm} \times 150 mm diameter) was pushed into the sediment to a depth of 100 mm , and capped at the bottom. The mesocosm core was gently filled with bottom water collected from the site, and then the top cap put in place. Cores were gently transferred to a Nally bin filled with site water, ice-blocks were added to keep the samples cool during transport and cores were kept upright and as stable as possible to allow sediments to settle. Cores were transported to the laboratory as soon as practicalable. At the laboratory, cores were carefully transferred to the incubation system, large round tubs filled with Taroona seawater, and allowed to equilibrate; chiller units kept the incubation temperatures equivalent to that of the site bottom water.

In addition sufficient site water to incubate all cores (approx. 400 L ) was collected from just above the sediments. Site water was collected from about 1 m off the bottom (to avoid resuspending sediments) at one sample location at each farm using a length of tubing and a bilge pump. The tubing was flushed with at least 60 L to ensure all water was from the bottom. Carboys were triple rinsed then filled with site bottom water. Site water was kept at a temperature equivalent to that of the relevant site bottom water (determined from CTD profile data) until used in incubation of mesocosms.

Sediments were allowed to equilibrate ( 24 hours) in the laboratory before incubations commenced. Mesocosm cores were flushed with site water for 20-30 minutes before incubations commenced. Incubations were undertaken in the dark, as the depths from which the samples were collected would preclude phytoplankton or microphytobenthic activity. Water samples were collected for dissolved nutrients (ammonia, nitrate + nitrite, phosphate) and pH analysis from each individual core at the start and end ( ∼4\sim 4 hours) of the incubation to estimate sediment nutrient fluxes. Dissolved oxygen was monitored at the start and end of each incubation using a Hach HQ40d meter with LDO101-15 probe. Data loggers continuously monitored O2\mathrm{O}_{2} levels ( 3 minute intervals) in one or two cores from each treatment. Temperature was maintained at a level approximately equivalent to field conditions using programmable heater/chiller units. Nutrient samples were analysed by Analytical Services Tasmania, with a minimum reporting limit of 2μ g/L2 \mu \mathrm{~g} / \mathrm{L} as N or P . Nutrient flux rates in μmol/m2/hr\mu \mathrm{mol} / \mathrm{m} 2 / \mathrm{hr} were calculated from the incubation time, the measured height of water in each core, the surface area of the core and the concentration changes during incubation (Eriksen et al., 2012).

2.6 Statistical Analysis

Univariate data were analysed by Analysis of Variance (ANOVA) with homogeneity of variances being checked using box-plots. Data were untransformed. A two-way fixed effects model ANOVA with factors location and time was used to assess variations in particle size, organic content, sulphide concentration, macroinvertebrate abundance, number of species and Shannon diversity. One-way ANOVA to test for differences in sediment fluxes and metal loadings between locations.

Multivariate analyses were conducted on the community assessment results to identify patterns within the data, with the results displayed as ordination plots using multidimensional scaling (MDS). Other than where indicated, replicates were combined and the data was subjected to a square root transformation in order to adjust the importance of species dominants. Two-way crossed analysis of similarities (ANOSIM) was conducted on the a priori groupings of farm, compliance and reference at each time to test within and between these a priori groups. SIMPER analysis was used to determine if any particular species were indicative of these patterns (Clarke and Warwick, 2001). All the multivariate analysis techniques were included in the Plymouth Routines in Multivariate Ecological Research (PRIMER) software package.

Sample analyses were specifically designed to address the specified study aims:

  1. Determine whether sediments within and around the Gunpowder Jetty lease have recovered such that they are consistent with reference conditions.
  1. Determine whether there is any significant difference in community composition between any of the sampled sites
  1. Determine whether there is any significant difference in functional performance between the reference and previously farmed cage sites

3. Results and Discussion

3.1 Particle Size Distribution

The sediment particle size was generally fine at all sites, with the silt-clay fraction ( <0.063 mm<0.063 \mathrm{~mm} ) comprising more than 50%50 \% of the sediment at all sites and times. At the end of farming (1999) there was no significant difference in the percentage silt-clay in the sediments from any sites, even though levels at the compliance sites were lower than in the reference sediments ( p=0.07,df2,7, F=3.95p=0.07, \mathrm{df}_{2,7}, \mathrm{~F}=3.95 ) (Fig. 2). The percentage silt clay in sediments has not changed markedly at the farm sites since over time, but would appear to have declined slightly at the reference sites over this period, with the average dropping from 98%98 \% (±0.5)( \pm 0.5) to 62%(±2.7)(p<0.01,df2,6, F=8.0962 \%( \pm 2.7)\left(p<0.01, \mathrm{df}_{2,6}, \mathrm{~F}=8.09\right.
img-3.jpeg

Figure 2. Average percentage silt/clay ( <0.063 mm<0.063 \mathrm{~mm} ) (s.d.) for farm, compliance and reference sites in 1999, 2002 and 2012.

3.2 Organic Matter Measurement

There was no significant difference between organic matter levels at farm, compliance or reference within each year, nor was there any difference in farm levels between 1999 and 2002 (Fig. 3). However, over time there has been a marked decline in organic matter content at the farm sites since 1999 ( p<0.01,df3,6, F=18.02p<0.01, \mathrm{df}_{3,6}, \mathrm{~F}=18.02 ).
img-4.jpeg

Figure 3. Percentage organic matter content (s.d) for farm, compliance and reference sites in 1999, 2002 and 2012.

3.3 Sulphide

Sulphide levels were very high at the farm stations (both at 1 and 4 cm depth) at the point at which farming ceased (1999) but dropped at the end of the initial recovery period (2002) and have remained low ever since (Fig.4).Sulphide levels at the reference stations have consistently been much lower, levels being negligible in 2002 and 2012. Consequently the interaction between location and time was highly significant ( p<0.01,df2,6 F=147.55p<0.01, \mathrm{df}_{2,6} \mathrm{~F}=147.55 ).

Whilst the sulphide concentrations at the farm stations remained higher than that recorded at the reference stations at all times, it is important to note that levels have decreased markedly since farming ceased and are consistent with levels reported from un-farmed areas elsewhere (Macleod et al., 2004a,b).

img-5.jpeg

Figure 4. Sulphide levels (s.d) at a) 1 cm and b) 4 cm for farm and reference sites in 1999, 2002 and 2012 (Note that sulphide was not sampled at compliance stations).

3.4 Macrofaunal Diversity

There were differences in the average number of species recorded from the farm, compliance and reference sites over time ( p<0.01,df2,66 F=11.68p<0.01, \mathrm{df}_{2,66} \mathrm{~F}=11.68 ) (Fig. 5). This appears to be associated with the compliance site; in 1999 more species were recovered from the compliance locations than in 2012(p<0.05,df1,22 F=6.57)2012\left(p<0.05, \mathrm{df}_{1,22} \mathrm{~F}=6.57\right) and at the end of the study in 2012 there were significantly more species at the compliance site than either the farm or reference sites ( p<0.01,df2,66 F=11.68p<0.01, \mathrm{df}_{2,66} \mathrm{~F}=11.68 ). However, it is important to note that these differences were still relatively minor with the lowest average number of species reported being 10 and the highest 25 .
img-6.jpeg

Figure 5. Average number of species (s.d.) at farm, compliance and reference sites in 1999, 2002 and 2012.

Total abundance has decreased markedly at farm sites since the cessation of farming ( p<0.01,df2,33 F=6.93p<0.01, \mathrm{df}_{2,33} \mathrm{~F}=6.93 ) (Fig.6). Similarly abundances have also declined significantly at the compliance sites ( p<0.01,df1,22 F=11.68p<0.01, \mathrm{df}_{1,22} \mathrm{~F}=11.68 ). Although abundances were very much lower at the

reference sites than at either compliance or farm sites in 1999 and 2002 average abundances at these sites were still slightly elevated relative to the results for 2012 ( p<0.05p<0.05, df2,27 F=3.48\mathrm{df}_{2,27} \mathrm{~F}=3.48 ). However, in 2012 there was no significant difference between any of the sampling positions ( p=0.13,df2,33 F=2.15p=0.13, \mathrm{df}_{2,33} \mathrm{~F}=2.15 ).
img-7.jpeg

Figure 6. Average abundance per m3\mathrm{m}^{3} (s.d) at farm, compliance and reference sites in 1999, 2002 and 2012.

Total diversity improved markedly at the farm sites between 1999 and 2002 ( p<0.01,df2,27p<0.01, \mathrm{df}_{2,27} F=15.41\mathrm{F}=15.41 ) and has remained at a similar level ever since (Fig. 7), although the variability observed at the farm site in 2002 appears to have stabilised. Once again, in 2012 there was no significant difference between any of the sampling positions ( p=0.21,df2,33 F=1.60p=0.21, \mathrm{df}_{2,33} \mathrm{~F}=1.60 ).
img-8.jpeg

Figure 7. Shannon diversity (s.d) at farm, compliance and reference sites in 1999, 2002 and 2012.

3.5 Nutrient Flux Measures

Where waste food and faeces are present in the benthos, they will be mineralised in the sediments and released directly back into the water column as soluble forms that can fuel algal growth (eutrophication) in the water column. The transformation and cycling of nutrients (oxygen, nitrogen, phosphorous) between sediments and the water column is therefore intimately linked to the rate of organic enrichment and degradation.
Consequently the rate at which these processes occur can provide some insight into the condition and performance of benthic habitats.

Generally, high ammonia and phosphorus fluxes from the sediment reflect organically enriched and oxygen depleted sediments. Similarly, fluxes of nitrate in the opposite direction, i.e. from the water column into the sediments, may reflect organically enriched and oxygen depleted sediments. This is because the process of nitrification that creates and supplies nitrate for denitrifying microbes in the sediments is limited in oxygen depleted sediments and the “denitrifiers” must therefore rely on nitrate sourced from the water column.

Although the average NH4\mathrm{NH}_{4} flux at the site 150 m from the cage was less than half that under the cage the replicate results were highly variable and so not significantly different (Fig. 8). Furthermore, these fluxes are very low (by several orders of magnitude) in comparison to fluxes previously recorded under or near operational sea-cages (e.g. Holmer and Kristensen, 1994; Christensen et al., 2000; Holmer et al., 2002, 2003;Bissett et al., 2009; Erickson et al., 2012), and within the range of fluxes reported from reference conditions (e.g. Christensen et al., 2000; Morrisey et al., 2000;Bissett et al., 2009; Erickson et al., 2012). Similarly, the fluxes of NOX,NO2,NO3\mathrm{NO}_{X}, \mathrm{NO}_{2}, \mathrm{NO}_{3} and PO 4 measured at the farm and reference sites are consistent with fluxes reported from reference conditions elsewhere (e.g. Christensen et al., 2000; Morrisey et al., 2000;Bissett et al., 2009; Erickson et al., 2012), including previously reported fluxes from other areas remote from fish farming within NW Bay (Macleod et al., 2008). There was no significant difference between the flux rates for any of the nitrogen species at the previously sampled cage and reference positions. Nitrate, ammonia and NOx\mathrm{NO}_{x} all indicate positive effluxes of nutrients from the sediments whilst nitrite was generally absorbed.

There were significant differences in phosphate (PO4) concentrations (df 1,5, F=23.80, p=0.008\mathrm{p}=0.008 ), but the overall rate was once again very low; the observed difference well within the flux range that might be expected from normal, oxygenated fine sediments. Fluxes from actively farmed sediments easily being an order of magnitude higher (Ericksen et al., 2012.
img-9.jpeg

Figure 8. Nutrient flux rates (s.d) in 2012 associated with sediments collected from the original cage transect sampling positions ( -10 m and 150 m ).

Oxygen flux rates provide a result that is consistent with the interpretation of the other nutrients (Fig. 9). Oxygen is being drawn into sediments but not at an excessive rate, which suggests a healthy respir benthic ecology, and there is no significant difference between the levels at the farm and at the reference locations.
img-10.jpeg

Figure 9. Oxygen flux rates (s.d) in 2012 associated with sediments collected from the original cage transect sampling positions ( -10 m and 150 m ).

The water overlying the sediments was well oxygenated and although there was a fairly rapid drop off in oxygen levels within the first few millimetres of the sediment, this is as might be expected in such fine (silt-clay) sediments (Fig. 10). There is no evidence that the farm sediments differed significantly from those of the reference locations.
img-11.jpeg

Figure 10. Oxygen penetration depths (s.d) for farm ( -10 m ) and reference ( 150 m ) sites in 2012. Sediment interface is shown at 0 mm .

3.6 Macrofauna

The ordination of the sites shows that there were significant differences in the community structure of the farm sites in 1999 compared to all other sites and times (Fig. 11). The ordination plot and associated ANOSIM analysis shows that in 2002 (after 3 years recovery) there had been considerable recovery at the farm positions such that for the most part the farm sites were indistinguishable from the reference, and that 10 years on the community structure and variability at the farm compliance and reference sites was indistinguishable but that these communities were as a whole different to those sampled on previous years, such that there appears to have been a general temporal shift.
img-12.jpeg

Figure 11.MDS ordination plot (Stress=0.17) showing the relative positions of farm, compliance and reference sites based on the benthic community composition of all replicate samples collected in 1999, 2002 and 2012.

Examining the nature of these changes overall shows that there were significant differences between years (Global R=0.478,p<0.05R=0.478, p<0.05 ) as well as between the different categories of site (Global R=0.234,P<0.05R=0.234, P<0.05 ), but the differences between years were greater (Global RR values). The communities at the farm sites at the end of production (1999) were dominated by opportunistic species indicative of organic enrichment (i.e. Capitella capitata, Malacoceros tripartitus and Nebalia longicornis), which together accounted for more than 86%86 \% of the overall similarity of communities at these sites and were quite different to those at reference locations, where the community was characterised by the brittle star Amphiura elandiformis, previously established as a useful indicator of clean sediment conditions and a number of other filter and surface deposit feeding polychaetes, bivalves and crustaceans more commonly associated with uncontaminated environmental conditions such as (Table 2). In contrast there was a lot of overlap in the communities at the farm, compliance and

reference locations in 2012, with differences for the most part being associated with changes in abundances of these key species rather than in species replacements. The species mix also suggests relatively good ecological conditions, with no obvious dominants and a broad range of functional types represented.

Table 2. SIMPER output for the full community assessment indicating the average within group similarity and the average abundance per group, variability (st.dev. ratio), individual % contribution to the group similarity and cumulative % similarity of the key species in each of the main groups over time.

Group 1999 Farm
Average similarity: 31.73

Species Av.Abund Sim/SD Contrib% Cum.%
Capitella sp.(MoV 2558) 38.14 0.92 51.21 51.21
Malacoceros tripartitus 21.73 0.70 25.82 77.03
Nebalia longicornis 7.90 0.51 9.64 86.67
Corbula gibba 4.23 0.27 3.28 89.95
Nassarius spp 3.43 0.39 2.70 92.65

Group 1999 Comp
Average similarity: 37.37

Species Av.Abund Sim/SD Contrib% Cum.%
Corbula gibba 12.00 1.34 19.06 19.06
Lysilla jennacubinae 9.22 1.03 13.84 32.90
Euphilomedes spp. 8.59 1.19 13.42 46.33
Ostracod sp. 8.69 1.03 10.62 56.95
Theora lubrica 7.71 0.99 10.00 66.95
Capitella sp.(MoV 2558) 9.17 0.22 3.81 70.76
Nephtys spp. 2.61 0.66 2.83 73.59
Nassarius spp 2.71 0.66 2.82 76.41

Group 1999 Ref
Average similarity: 28.11

Species Av.Abund Av.Sim Sim/SD Contrib% Cum.%
Lysilla jennacubinae 7.50 5.28 0.80 18.79 18.79
Amphiura elandiformis 7.14 4.57 0.96 16.26 35.05
Theora lubrica 7.96 4.25 0.56 15.13 50.18
Corbula gibba 7.78 3.81 0.82 13.54 63.73
Nucula pusilla 7.43 3.64 0.60 12.95 76.67
Lumbrinereidae spp. 4.49 1.72 0.64 6.11 82.78
Euphilomedes spp. 6.06 1.27 0.41 4.53 87.31
Nassarius spp 3.39 0.97 0.42 3.47 90.77

Group 2012 Farm
Average similarity: 55.81

Species Av.Abund Av.Sim Sim/SD Contrib% Cum.%
Callianassa limosa 7.52 11.80 7.47 21.14 21.14
Theora lubrica 7.52 11.80 7.47 21.14 42.27
Phoronid sp. 5.64 6.50 1.06 11.65 53.93
Euchone limnicola 5.64 6.35 1.07 11.38 65.31
Lysilla jennacubinae 5.01 4.77 0.84 8.55 73.86
Corbula gibba 5.01 4.64 0.84 8.32 82.17
Nassarius spp 4.38 3.37 0.68 6.04 88.21
Melitid sp. 3.39 1.76 0.42 3.15 91.36

Group 2012 Comp
Average similarity: 44.09

Species Av.Abund Av.Sim Sim/SD Contrib% Cum.%
Theora lubrica 6.89 7.90 1.98 17.91 17.91
Euchone limnicola 6.26 6.12 1.39 13.87 31.78
Callianassa limosa 5.64 5.70 1.06 12.92 44.69
Nassarius spp 5.01 4.02 0.81 9.11 53.81
Corbula gibba 4.38 3.71 0.66 8.41 62.21
Phoronid sp. 4.38 3.43 0.64 7.77 69.98
Nemertean sp. 4.38 3.09 0.64 7.01 77.00
Lysilla jennacubinae 4.38 2.84 0.66 6.44 83.43
Kalliapseudes sp. 3.13 1.15 0.41 2.61 86.04
Mediomastus australiensis 3.13 1.13 0.42 2.57 88.61
Oedicerotid spp. 2.51 0.77 0.31 1.76 90.37

Group 2012 Ref
Average similarity: 57.52

Species Av.Abund Av.Sim Sim/SD Contrib% Cum.%
Callianassa limosa 7.52 11.50 7.56 19.99 19.99
Theora lubrica 7.52 11.50 7.56 19.99 39.99
Corbula gibba 6.26 7.98 1.41 13.87 53.86
Nassarius spp 6.26 7.67 1.42 13.33 67.19
Amphiura elandiformis 5.01 4.91 0.84 8.53 75.71
Euchone limnicola 4.38 3.47 0.67 6.03 81.74
Lysilla jennacubinae 3.76 2.66 0.53 4.63 86.37
Phoronid sp. 3.76 2.49 0.54 4.34 90.71

Comparison of the community composition data for each of the locations sampled in 2012 shows a high degree of overlap between the sites, such that there was no significant difference between the farm, compliance and reference communities (Global R=0.113R=0.113, p=0.17\mathrm{p}=0.17 ) (Fig.12).
img-13.jpeg

Figure 12.MDS ordination plot (Stress=0.13) showing the relative positions of farm, compliance and reference sites based on the benthic community composition at each of the locations sampled in 2012.

3.7 Video Assessment

The principal components analysis shown in Figure 13 shows the separation of the sites based on the video scores. The first two axes account for approximately 74%74 \% of the variation in the dataset with 41%41 \% of the variation associated with PC1. The key features separating the stations along PC1 are an increased presence of both Brittlestars and Small Fish. Differentiation of sites along PC2 is not as well defined and generally results from a reduction in Brittlestars along with declines in abundance of Molluscs and Burrow Density. Whilst there does seem to be some evidence that the video footage in 1999 might differ from that in 2012, this is at a very low level and may simply be a data artefact (i.e. a result of improved video capture techniques and clarity of images in 2012), consequently it may be wise not to read too much into this. In addition it is important to note that whilst the presence of brittlestars in a video survey will tend to infer good environmental conditions, the absence of brittlestars does not mean that conditions are in any way degraded.

When the 2012 data is considered in isolation there is little to suggest any major differentiation between the site categories (reference, compliance and farm) (Fig. 13). None

of the usual farm impact assessment criteria (e.g. gas bubbles, sediment colour, Beggiatoa, pellets/ faeces) suggested any significant impacts at any of the positions sampled in 2012.
img-14.jpeg

Figure 13. Principal components analysis (PCA) showing the relative positions of farm, compliance and reference sites based on the video assessment at each of the locations sampled in 1999 and 2012.

3.8 Copper and Zinc loadings

As copper based antifoulants were used in the latter years of production at the Gunpowder site in the follow-up study we assessed total sediment load and porewater concentrations of both copper and zinc.

Total loadings for surface sediments were well below the ANZECC guidelines for sediment copper, with the highest level recorded being 43mg/kg43 \mathrm{mg} / \mathrm{kg} (i.e. low level guideline of 65mg/kg65 \mathrm{mg} / \mathrm{kg} ), and there was no significant difference between levels at the farm and reference locations ( p=0.11,df1,5 F=3.52p=0.11, \mathrm{df}_{1,5} \mathrm{~F}=3.52 ) (Fig.14). Sediment levels were within the range of expectation, being consistent with concentrations reported for reference sites through ongoing salmon farm monitoring (Macleod and Eriksen, 2009) and in a similar range to levels previously reported from SE Tasmania (HEST 2000, Jones et al., 2003, Butler, 2005, Macleod and Helidoniotis, 2005). Porewater copper levels at the reference location were effectively below the detection level of the analysis ( <1ug/L<1 \mathrm{ug} / \mathrm{L} ) and although the farm levels may appear a little higher the variability was such that these levels were not significantly different from the

reference sites ( p=0.10,df1,19 F=2.91p=0.10, \mathrm{df}_{1,19} \mathrm{~F}=2.91 ) and were also very low and consistent with previously reported levels.
img-15.jpeg

Figure 14. Copper levels (s.d) in a) sediment and b) porewater at the original cage transect sampling positions ( -10 m and 150 m ) in 2012.

Sediment zinc loadings were consistently above the ANZECC low level guideline (i.e. 200 mg/kg\mathrm{mg} / \mathrm{kg} ), with no significant difference between the farm and reference locations ( p=0.52p=0.52, df1,5 F=0.45\mathrm{df}_{1,5} \mathrm{~F}=0.45 ) (Fig.15). The overall average zinc concentration in the sediments sampled being 263mg/kg263 \mathrm{mg} / \mathrm{kg} (s.d. 16.68). However, this is consistent with a previous report from Jordan et al. (2002) which noted that both lead and zinc were elevated in NW Bay such that 22%22 \% of sites sampled in NW Bay exceeded the low level guideline for zinc. Jordan et al. suggested that the sites with higher zinc levels were generally those with finer sediments from the deeper regions of the Bay and concluded that the source was most likely associated with the local catchment.

Porewater zinc levels were low and at levels unlikely to elicit a biological effect (Stuart Simpson CSIRO-CERC, pers. comm.), and again did not differ significantly between the reference and farm sites ( p=0.40,df1,19 F=0.72p=0.40, \mathrm{df}_{1,19} \mathrm{~F}=0.72 ).
img-16.jpeg

Figure 15. Zinc levels (s.d) in a) sediment and b) porewater at the original cage transect sampling positions ( -10 m and 150 m ) in 2012.

4. Conclusions

In summary it would appear from analysis of the selected environmental condition variables that in 2012 the sediments around the Gunpowder lease have largely recovered.

Sediment geochemical parameters show some small differences but these do not appear to translate to ecological differences or any major changes in sediment biogeochemistry or function. Levels of organic matter are broadly equivalent, with no evidence of any residual organic loadings. Sediments in this region are very fine (predominantly silt-clay), and the farm and reference locations were indistinguishable. The compliance locations had a slightly lower proportion of fines than observed at either the reference or farm locations, but this was quite minimal and consistent within the range of particle sizes that might be expected in this region. Previous sampling in NW Bay would suggest that sediment particle sizes are generally fine with a high proportion of silt-clay but that the composition can vary markedly with depth, with a greater proportion of coarser material in shallower regions (Macleod et al, 2008). Sulphide was still detectable in the sediments collected at the farm locations, but the levels were very low (close to the detection limit for the analysis) and do not appear to have any major impact on the biological community. In terms of abundance, number of species, diversity and overall community structure, the farm sites were largely indistinguishable from the reference or compliance sites.

There is no evidence that any major or substantive changes have occurred at the farm sites following cessation of farming, nor is there evidence that the farming activity has had any permanent impact in this area. Although there does appear to have been a very slight change in community composition over time, this change was evident at all of the sampling stations, was at a very low level and would appear more likely to be a result of natural ecological dynamics. There does not appear to be any significant functional difference in the sediments at the farm sites relative to nearby compliance or reference sites, and the nutrient flux rates are consistent with expectations, suggesting a relatively healthy biological response for this type of sediment. It is also worth noting that over the period the intervening 10 year period since the farm site has been fully vacated there have been a number of other major changes to the local catchment and potential inputs to NW Bay (i.e. significant urban development, and upgrading and relocation of sewage treatment facilities).

5. References

211-214.