Diverse marine benthic communities and reduced anthropogenic contaminants near Scott Base (Hut Point Peninsula, Ross Island, Antarctica)

Antarctic marine ecosystems are largely thought to be among the planet’s least impacted, yet habitats adjacent to research stations can be heavily polluted. Despite long-term monitoring and remediation being high priorities for international environmental protection, the ecological responses to contaminants and stressors remains poorly characterised, limiting our ability to manage and reduce human impacts. This study compares epifaunal community composition at two sites close to Scott Base with a reference site further away. We couple these data with environmental characteristics, including current data, sediment properties, and contaminant concentrations within the sediment and in the tissues of two epifaunal species, both from this survey and those previously reported. Previously high concentrations of polychlorinated biphenyls and polyaromatic hydrocarbons are now undetectable and concentrations of heavy metals were mostly similar or reduced compared to previous data from 2002. High within-site variance suggests localised contamination footprints from being situated within a deposition zone and/or from the erosion of anthropogenic debris. Despite the persistence of some metals (arsenic, copper and lead) at one site, our study revealed high biodiversity at all three sites (22–28 taxa per 0.25 m2). Benthic community structure was influenced by a combination of factors, including sea ice characteristics, sediment type and habitat complexity. Overall, our study clearly highlights the influence of human activities on the benthos in adjacent marine habitats. The established monitoring protocols coupling diver and remote sampling will enable regular monitoring, filling a critical need for time-series data in order to detect long-term trends and interactions with climate drivers.


Introduction
Humans have affected all of Earth's ecosystems to varying degrees (Western 2001;Halpern et al. 2008;Ellis et al. 2021). Arguably among the planet's least human-impacted are high latitude marine ecosystems in Antarctica that are covered by ice for most of the year (Halpern et al. 2008;Aronson et al. 2011). The region's relative remoteness and inaccessibility has resulted in limited pollution, resource extraction (minerals, fishing) and habitat fragmentation making coastal marine biological assemblages in Antarctica important bellwether systems for studies of climate-related change (Aronson et al. 2011;Dayton et al. 2019).
Yet, habitats adjacent to current and former research stations are far from pristine. Sites adjacent to some stations have been polluted by shipping, spills, accidents, intentional dumping, fuel soot, seal hunting, sled dogs, and treated and untreated sewage discharge (Lenihan et al. 1990;Kennicutt II et al. 1995, 2010Stark et al. 2005Stark et al. , 2016Palmer et al. 2021). In addition, non-point source pollutants can be introduced through increases in sediment-laden freshwater runoff and wind-blown dust, which can further deliver bound anthropogenic contaminants to the marine benthos (Griscom et al. 2004;Birch et al. 2008). Remediation of environmental damage and long-term monitoring have been identified as high priorities by The Antarctic Treaty's Protocol on Environmental Protection (Kennicutt II et al. 2010;CEP 2021). However, ecological responses to many of these human contaminants and stressors remain poorly characterised, limiting our ability to manage and reduce human impacts in coastal areas.
Seventy-six research stations are operational in Antarctica, with the majority of these within coastal areas (COMNAP 2017). One of the largest is the United States' McMurdo Station (MCM) located on Hut Point Peninsula, which houses > 1000 people in the summer. Owing to its size and diverse range of activities including scientific research and related logistical support (shipping, aviation, ground vehicles, fuel storage, power generation, construction, solid waste processing, sewage treatment, etc.), coastal marine soft-sediment benthic habitats near MCM in Winter Quarters Bay have been heavily polluted (e.g., Kennicutt II et al. 1995;Negri et al. 2004Negri et al. , 2006. There, concentrations of hydrocarbons are comparable to those in some of the world's busiest and most contaminated shipping ports (Lenihan 1992).
The nearest research station to MCM in terms of geographical distance is New Zealand's Scott Base (SB). SB has a maximum summer population of ~ 80 people, a size that is more typical of the majority of research stations in Antarctica. Whilst habitats near larger stations such as MCM have been extensively monitored, much less is known about the levels of anthropogenic impact in the marine areas adjacent to smaller stations, such as SB. The only previous assessment of SB marine benthos focused on sediment contamination and showed elevated levels of hydrocarbons and heavy metals (Negri et al. 2004(Negri et al. , 2006. Contamination was lower than at nearby MCM sites (Anderson and Chague-Goff et al. 1996), and similar to levels reported elsewhere in Antarctica (Alam et al. 1993;Giordano et al. 1999;Santos et al. 2005). However, there has been no ongoing marine benthic contaminant monitoring near SB since the early work of Negri and colleagues.
Seafloor macroinvertebrate assemblages are regularly used for environmental monitoring because many of the species are sedentary, long-lived and sensitive to changes in water and sediment quality. Changes in Antarctic macrobenthic communities have been observed in response to natural disturbances such as variations in sea ice cover (Norkko et al. 2007;Lohrer et al. 2013) and ice scour (Gutt et al. 1996), as well as local anthropogenic impacts (Dayton et al. 1970;Conlan et al. 2004;Tin et al. 2009) and global climate change (Aronson et al. 2011). In the marine environment adjacent to SB, wastewater is discharged and seafloor observations have provided evidence of past anthropogenic pollution (glass, metal, wood, vehicle parts, etc.) (Negri et al. 2006;Hale et al. 2008). Understanding the impacts of SB activities on the local marine environment is imperative, particularly as sewage and metal and chemical contaminants have been shown to modify benthic communities elsewhere in multiple locations Palmer et al. 2021), reducing biodiversity, changing species abundances and altering the prevalence of dominant species (Conlan et al. 2004;Stark et al. 2014;Palmer et al. 2021).
The degree to which marine assemblages in the receiving environment are affected by contaminants can depend on many factors, including their proximity to the source, the relative sensitivities of the biota to differing concentrations and combinations, dispersion characteristics and the initial concentration of the pollutant/ contaminant. The latter two factors are considerably influenced by local currents and hydrodynamics and are therefore important elements for determining potential contamination footprints, but are currently not well characterised at SB.
The objectives of this research were to (i) compare epifaunal community composition and biodiversity at two sites close to SB and at a reference site further away, and (ii) contextualise the findings with current speed and direction data from the SB sites, elucidating contaminant transport potential and their relative vulnerabilities to anthropogenic stressors. We also (iii) compare the 2019 data to contamination information recorded previously by Negri et al. (2006) as an indication of potential improvement (or lack thereof).

Site locations and descriptions
Sampling took place at three coastal marine sites around Hut Point Peninsula, southern Ross Island during October-November 2019 ( Fig. 1; Table 1). Two of the sites were adjacent to SB and were chosen to coincide with potential contamination entry points from freshwater runoff and with sites previously sampled for sediment contamination (Negri et al. 2004(Negri et al. , 2006. These sites were Scott Base 1 (SB1) on the southern coast of Pram Point approximately 300 m westsouthwest of the SB sewage outfall, and Scott Base 3 (SB3) on the northern coast of Pram Point, northeast of the outfall. An initial aim was to establish up to five monitoring sites at increasing distances away from Scott Base. However, establishing additional sites around SB was impeded by practicalities including (1) a limited field work window, (2) very steep underwater terrain and thus limited seafloor area in the targeted 15-25 m depth zone, (3) sea ice conditions in front of Scott Base in 2019, including large pressure ridges directly over suitable working depths to the east of Pram Point, and (4) the Ross Ice Shelf to the north of SB which prevented the establishment of sites in that direction. Ultimately, Arrival Heights Site 1 (AH1) was chosen as a reference site. The establishment of a reference site between SB3 and AH1 was decided against due to past and recent , with information on summertime freshwater flows (yellow arrows) and positions of the RO intake/ discharge and sewage discharge (purple arrows, with * denoting the sewage outfall). Predominant current flow directions for the two sites with ADCP current meters are also shown for reference (see methods and results for more details), and the prevailing current direction on the western side of Cape Armitage in dark blue (Littlepage 1965;Barry et al. 1988). Note that the pressure ridges were much closer to shore in 2019 than they are in this image from a prior year Table 1 Locations and descriptions of sites Sea ice conditions includes thickness of land-fast sea ice, the snow layer atop the sea ice, and the platelet ice layer beneath the sea ice. The 2019 notes are based on the personal observations of the authors and 2002 notes were copied verbatim from Negri et al. (2006)  reports of high contamination within the sediment and benthic fauna in front of MCM in Winter Quarters Bay and towards Cape Armitage ( Fig. 1; Lenihan et al. 1990;Kennicutt II et al. 1995;Negri et al. 2004Negri et al. , 2006Palmer et al. 2022), and thick spicule mats covering the sediments at sites near Cape Armitage that made these areas unsuitable for contaminant sampling (Klein et al. 2014).

Site characterisation and sampling
Two acoustic doppler current profilers (ADCP; Nortek Signature 500) were deployed at both SB sites from October to December 2019 to capture information on current speed and direction. Both ADCPs were deployed in a downward orientation and set with a 2 min sampling period in 1 m vertical depth bins from the underside of the ice to the seabed. Custom frames and tripods were used to suspend the instruments so that the instrument transducer heads were ~ 0.5 m beneath the under-surface of the ice and clear of any brash/ platelet ice. Extra non-magnetic weights were added to the frames to maintain each instrument's vertical orientation throughout the deployment. Information on benthic habitats and epifaunal assemblage composition were collected using high resolution video. At each site, two 25 m transects were established at ~ 22 m depth. Along each transect, diver collected video was recorded from a height of 0.5 m above the seabed. Additional video footage of the seafloor surrounding each transect was taken remotely from the surface using a cabled submersible remotely operated vehicle (ROV). Multiple overlapping passes were made across the seabed transects at ~ 0.5 m depth contours between ~ 20 and 26 m in order to create a 2D orthomosaic image of each site.
To determine characteristics of surface sediments at each site, divers used sterile plastic jars to collect scrapes of surface sediment (upper 2 cm). Each replicate was a ~ 500 g composite of five individual sediment scrapes (jars). After collection, all replicates were kept in the dark and frozen until analysis. Four replicate samples were analysed for sediment contaminant concentrations (metals, PCBs, and hydrocarbons) and one composite replicate per site was analysed for particle size distribution, organic matter content, and pigment concentration (chlorophyll a and phaeophytin).
Samples of two sessile suspension-feeding taxa, the sponge Sphaerotylus antarcticus and the bivalve Laternula elliptica, were collected for tissue contaminant concentration (body burden) analyses. Data from sessile benthic organisms are likely to be more reflective of a site's history of exposure to contaminants than analyses of mobile species (e.g., fish such as Trematomus spp., Notothenia spp.) where site fidelity is less certain. Long-lived suspension-feeding taxa living on (e.g., sponges) and in (e.g., infaunal bivalves) seafloor sediments, and that filter large volumes of water over their lifespans, make excellent bio-indicator species due to the greater probability of concentrating diffuse contaminants (de Mestre et al. 2012;Padovan et al. 2012;Batista et al. 2014;Orani et al. 2022). In addition, S. antarcticus and L. elliptica were common to abundant at all three of the sites sampled and were previously collected for contaminant analyses by Negri et al. (2006), enabling temporal comparisons. Divers attempted to collect two other sponges (Homaxinella balfourensis and Mycale sp. in addition to S. antarcticus), but both additional species were absent at one or more of the three sites and were thus not analysed. At each site, a total of four replicate samples of 20-30 g blotted tissue wet weight were collected for S. antarcticus and L. elliptica. Sponge tissue was removed in situ using a plastic knife and L. elliptica specimens were removed whole from the sediment. After collection, each sample was placed into a double ziplock bag and frozen until analysis.

Current meter processing
ADCP data were downloaded, extracted from their raw formats, and averaged into 10 min intervals. A magnetic declination of 141.09° E was applied to the measured current direction to correct the readings to reflect true north and a pressure offset was applied to standardise depths relative to ambient air pressure at the seawater surface. Data were then assessed for quality based on internal measurements of instrument tilt, signal correlation and signal strength and removed if quality thresholds were not met. Although the ADCPs were able to record for 4.5 weeks, the measuring transducers on each instrument became progressively obstructed by platelet ice growth, leading to reduced sampling range followed by complete loss of reliable data. Site SB1 produced 21 days of reliable data, with the first 17 days spanning the entire water column and the remaining 4 days recording at increasingly shallower depths. Site SB3 produced 12 days of data for the entire water column before becoming fouled.

Video analysis
Analysis of the diver-collected video was done using individual frames rather than continuous footage (Cummings et al. 2006;Thrush et al. 2011;Cummings et al. 2018). The video along each transect was 'divided' into 10 equal time segments and still frames were grabbed at random from the first, third, fifth, seventh and ninth segments. This ensured the independence of replicates and that the entire length of each transect was represented in the sampling. If any of the frames were of insufficient quality for analysis (e.g., blurry), a new one was selected. Eight video frames were analysed per transect (i.e., n = 8 per transect and n = 16 per site) by one individual to minimise observer bias. The still frames were imported into Adobe Photoshop software and scaled to encompass a standard field of view of 56 cm wide by 45 cm high. A computerised grid consisting of 80 equal sized squares was then superimposed onto each image. Each image was analysed in several stages. First, each of the 80 grid cells was classified as being dominated by boulder, rock, cobble, pebble, gravel, sand, shell hash, or obscured, with sizes classes following Wentworth (1922). The proportional representation of each substrate class in each image was calculated by dividing the number of cells where it was dominant by the total available cells (i.e., 80 minus the "obscured" cells). For large sessile animals (that could cover multiple grid cells in an image) and abundant small sessile animals (difficult to count individually), a technique similar to that used for the base substrate was used to estimate their presence and relative abundance. That is, the number of cells where they were present was divided by the total available cells. Most of the large taxa (mobile and sessile) in the images were able to be identified to genus or species levels, whereas smaller taxa were more difficult to identify (the maximum number of taxa identified per frame was 11). Mobile animals that were visible in the images (e.g., seastars, brittle stars, urchins, sea spiders) were counted individually rather than using the proportional representation method.

Sediment and animal tissue analysis
To determine sediment characteristics, sediment samples were analysed for particle size distribution, organic matter content, and pigment concentrations using previously published standard methods (e.g., Lohrer et al. 2004) within one month of collection. Samples for particle size distribution were immersed in a 9% hydrogen peroxide solution to digest any organic material and grain size fractions were determined by a wet sieving/pipetting method. Organic matter content was assessed via combustion of dry sediment (48 H at 60 °C) in a muffle furnace (5.5 h at 400 °C) and expressed as percent dry weight loss on ignition. Frozen chlorophyll a and phaeophytin content samples were freeze-dried, with pigments then extracted in ethanol and analysed spectrophotometrically (Turner Designs 10-AU) (Sartory 1982), yielding pigment concentrations.
Analysis of contaminants was conducted on sediment and animal tissue, as per Negri et al. (2006), using standard methods in an accredited laboratory in Hamilton, New Zealand. Sediments were sieved and the fraction < 2 mm retained, before being dried to a moisture content of 2-5% and mechanically ground to produce homogenous samples for analysis. Bivalve tissue was removed from the shell, and all animal tissue samples were chopped, minced and blended to produce homogenous samples. Individual tests were conducted on each sample for six metals (Arsenic (As), Cadmium   were examined using sonication extraction, solid-phase extraction clean up followed by GC-MS analysis. Sediments were additionally analysed for seven petroleum hydrocarbons (C7-C9, C10-C11, C12-C14, C15-C20, C21-C25, C26-C29 and C30-C44) using sonication extraction, silica clean up and GC-FID analysis, and animal tissues were additionally analysed for total lipid content using gravimetric analysis.

Statistical analyses
The seafloor biodiversity data from SB1, SB3 and AH1 were compiled into a sample-by-taxon matrix and imported into PRIMER (version 7) software for analysis. Differences in seafloor community composition were assessed from Bray-Curtis similarities following a fourth-root transformation of the raw abundance data (to down-weight the influence of abundant taxa on the results). Ordination plots were created (non-metric multidimensional scaling, nMDS) and tests for significant differences in community composition among sites were conducted (Permutational Multivariate Analysis of Variance, PERMANOVA, with "Site" as a fixed factor in post-hoc tests; Similarity Percent, SIMPER). Environmental data collected at each site were also compiled and used to describe differences among sites (e.g., Principal Components Analyses, PCA) (Clarke 1993;Clarke et al. 2006).

Site characteristics
There were differences in sea ice characteristics between sites (Table 1). The sea ice was thickest at SB3 and thinnest 1 3 at AH1 (Table 1). AH1 had almost no snow above the sea ice and almost no platelet/brash ice underneath (Table 1), resulting in higher underwater light availability.
Orthomosaic maps of the seafloor at each site produced from composite images from the ROV (Fig. 2), also revealed differences between sites. The largest differences were observed between the reference site AH1 and the two SB sites (SB1, SB3). The AH1 seabed was dominated by cobbly volcanic rock interspersed with scoria rubble and gravel and few large boulders. Anthropogenic debris was found within the site, consisting of old rope and floats and a gangplank 300 m south of the transect which fell off a ship in 1960 (Dayton et al. 2016). The two SB sites were dominated by a mixture of moderately unconsolidated volcanic scoria rubble and gravels with interspersed rocky outcrops and were steeply sloped (~ 40° compared to ~ 30° at AH1). Both sites contained bivalve shell hash, predominantly from the infaunal bivalve, Limatula hodgsoni, and high cover of an unidentified filamentous 'fluffy' turf, which was likely comprised in part from the silica spicules of sponges. SB3 had larger patches of sediment and slightly more ambient light than SB1 (diver obs.). The biggest difference between the two SB sites was the large amount of anthropogenic debris occurring at SB3. This included glass bottles, bamboo flag poles, Fig. 2 Orthomosaic of one 25 m transect at a Arrival Heights (AH1), c Scott Base 1 (SB1) and e Scott Base 3 (SB3), along the 22 m depth contour. Images b, d, and f shows an enlarged image from each transect at AH1, SB1 and SB3, respectively. The top of each image is towards the shore, and because of steep sloping of the seabed, is slightly shallower than the bottom of the image rusting metal drums, wooden boards, and food waste (corn cobs). Rusting and disintegrating metal turned the sediment an orange/red colour in patches at SB3. Sediment characteristics also differed by site (Table 2). For example, sediment mud content was 5-7 times higher at SB3 (62%) compared to the other two sites (9-12%). Sediment organic matter content, which is often positively correlated with mud content, was highest at SB3 (Table 2). Sediment chlorophyll a and phaeophytin content were up to 10 times higher in sediments at AH1 relative to SB1 and SB3 (Table 2).

Current strength and direction
Hut Point Peninsula, where two ADCPs were deployed, has diurnal tides with a range of ~ 1 m on a spring tide. Currents at SB1 were relatively weak (median and mean < 6 cm/s; Online Resource 1) and uniform from the surface to the seabed (Fig. 3C). However, brief pulses of strong flow were recorded (18-20 cm/s). An oscillating eastward-westward flow regime was observed (Fig. 3A), though the easterly currents were generally stronger than the westerly currents, suggesting the predominant residual flow pattern is from SB1 towards the SB wastewater discharge point (Fig. 1). Currents at SB3 did not follow a bi-directionally pattern but were instead predominately in a south-westerly direction. Median, mean and near maximum currents were weaker, on average, than those at SB1, with bottom water currents almost half of those at SB1 (Online Resource 2). Additionally, a more vertical current structure was observed at SB3, with currents faster underneath the ice, and slower near the seafloor (Fig. 3D).

Contaminant concentrations within the sediment
With very few exceptions, concentrations of all PAH and PCB congeners in sediments at all sites were below the detection thresholds of the analytical procedures used. Total PAH and PCB concentrations (i.e., all congeners combined) were also below detection thresholds at all sites, indicating very little existing organic contamination at the study sites in 2019. The only exception was petroleum hydrocarbons, which, while below detection limits at SB1 and AH1 (< 70 mg/kg dry weight), were present in moderate concentrations at SB3 (average 157.5 ± 52.2, range 90-300 mg/kg dry weight).
Heavy metals were detected in the sediments at all three sampling sites (Fig. 4, Online resource 3). Average concentrations of As, Cu, Pb, Zn, and Cd were generally highest at SB3. In contrast, the concentration of Hg was highest at AH1, with all replicates at this site exceeding the indicative sediment toxicity Default Guideline Value for Hg of 0.15 mg/kg dry weight (Fig. 4).
Considerable among-replicate variation was detected at SB3. One of the four replicates had substantially higher concentrations of all metal species tested, with concentrations of some metals at (Cu) or above (As, Pb) indicative sediment toxicity guideline values (ANZECC & ARMCANZ 2000). The Pb concentration in this sample was 100 times higher than that of the other samples at the site, while As was ~ 10 times higher, and Cu, Cd and Hg were ~ 3 times higher.

Contaminant concentrations in animal tissues
Similar to the sediment analysis, concentrations of PAHs and PCBs within animal tissues were generally below detection limits. While some of the individual PCB congeners  were just above the detection threshold in L. elliptica tissues sampled at AH1, the concentrations were still very low, with total PCBs < 0.02 mg/kg tissue in all replicates of both species at all three sites.
Heavy metal contaminants were detected in the tissues of both suspension feeding species (Online resource 3, Fig. 4). Cd, which was in very low concentrations in sediment, was relatively concentrated in the tissues of both species (Online resource 3). Across all metal species, concentrations in L. elliptica tended to be highest at SB1, intermediate at SB3, and lowest at AH1 (Fig. 4). No site-related pattern was apparent for S. antarcticus.

Epifaunal community structure
The ecological community data gathered from frame grabs of diver-collected video revealed all sites had high diversity, with a total of 22-28 taxa per site (see Online Resource 4 for all species and abundances). Although 28 individual taxa were recorded at SB3, the site with the highest average richness, evenness and diversity of taxa per frame was SB1. AH1 had the lowest average richness, abundance and diversity per frame (Table 3).
There were distinct differences in community assemblages between all sites (P perm < 0.001 for all pairwise contrasts, PERMANOVA; Fig. 5), however, the two SB sites were more similar to each other (31% similarity) than either one was to AH1 (13 and 17% similarity, SIMPER; Table 4).

Discussion
Marine benthic habitats adjacent to Antarctic research stations have the potential to be highly impacted. Heavy metals, PAHs and PCBs were detected in sediments and in animal tissues in the vicinity of SB in 2002 (Negri et al. 2006), presumably the consequence of base activities in years prior. Base practices, including sewage treatment and pollution management, have improved substantially since 2002, and contaminant levels within animal tissues were found to be much lower in 2019. For example, PAHs and PCBs were below detection limits in 2019, and sediment heavy metal concentrations were mostly present at similar or lower levels compared to 2002, except for an increase in some metal species (As and Pb) at SB3. However, to our knowledge, there have been no assessments of seafloor community characteristics and diversity until now. Our study revealed diverse benthic communities at 22 m depth at all three sites (a total of 22-28 taxa per sample per site), including SB3, which was once a dump site and which current flows and grain size indicate may be a contaminant depositional zone. The benthic communities at both SB sites were considerably different to AH1, being characterised by relatively high abundances of the brittle star O. antarctica, cone sponges P. invaginata, and sea spiders (Pycnogonida). Comparatively, these species were rare to absent at AH1, and instead the soft coral C. frankliniana was the most abundant epifaunal taxon. Differences in benthic community composition between sites can be partially attributed to variation in benthic substates  (2000)) which "indicate the concentrations below which there is a low risk of unacceptable effects occurring, and should be used, with other lines of evidence, to protect aquatic ecosystems". DGVs for Zn and Cd are off scale and therefore not shown. N.B. the difference in scale on the y-axis between sediment and tissue data, and Cd and Pb are in a different order for sediment and tissue concentrations to allow for the lower concentrations to be visible and habitats between the three sites, as highlighted by the orthomosaic maps (Fig. 2). For example, the two SB sites had much steeper slopes and were dominated by a mixture of unconsolidated volcanic scoria rubble compared to cobbly volcanic rock at AH1. Differences and changes to seafloor complexity (both abiotic and biotic) can influence macrobenthic communities in multiple ways. The mechanisms may involve reducing the erosion of fine sediment from the bed, creating biogenic habitat for encrusting organisms, and providing protection to other organisms against predation Smith et al. 2015). Analysis of the surface sediment further differentiated the three sites, especially between the two SB sites. Sediments at SB3 had high proportions of silt and clay, with 63% mud content compared to < 12% at the other two sites (Table 2). High rates of suspended fine sediment deposition leading to elevated muddiness are likely to influence benthic communities both directly (e.g., by smothering small organisms and interfering with their food gathering) and indirectly (e.g., by altering near-bed oxygen, ammonium, and sulphide concentrations and sediment biogeochemistry) (Ellis et al. 2002;Marinelli et al. 2002;Lohrer et al. 2004;Cummings et al. 2009). When coupled with current data, the high sediment mud content and organic enrichment observed at SB3 is unlikely to have come from the sewage outfall. The large-scale circulation pattern in McMurdo Sound moves water southward past Cape Evans towards Cape Armitage, whereupon it turns east north-east after Cape Armitage to flow under the Ross Ice Shelf (Robinson et al. 2012), and would result in SB3 largely being in the lee of the prevailing north-easterly flow and potentially in a recirculation eddy. This and the weaker currents recorded at SB3, especially near the sea floor, suggest that it may be a fine sediment deposition zone, which is further supported by the high proportion of silt and clay.
The deposition of silt and clays can additionally alter benthic macrofaunal communities via the transportation of bound contaminants (Affleck et al. 2014). Previous determination of contamination by Negri et al. (2004; found evidence of sediment and animal tissue contamination at both SB sites including PAHs, PCBs and metals (Cu, Cd, Zn, Pb and As). Our assessment at three sites in 2019 revealed PAHs and PCBs are now below detection limits or present in very low concentrations, however, the presence of metals was evident at all sites (Fig. 4). Within the sediment, metal concentrations have decreased at SB1 compared to Negri et al. (2006), but are the same (Cu, Cd, Hg), decreased (Zn) or increased by at least an order of magnitude (Pb and As) at SB3 (Online resource 3). These sediment contaminant concentrations at SB are comparable to those reported adjacent to similar sized research stations and much lower than those reported at nearby MCM sediments (e.g., Lenihan 1992;Kennicutt II et al. 1995;Stark et al. 2003;Santos et al. 2005;Kennicutt II et al. 2010;Stark et al. 2014;Palmer et al. 2021). However, the increased concentration of Pb and As in SB3 sediments and the high variance between samples at this site (Online resource 3) suggests, similar to other studies, the differential deposition of contaminant bound sediments and/or the degradation of anthropogenic debris in situ can result in persistent and patchy contamination hot spots (Kennicutt II 2003;Negri et al. 2006).
Contaminants can be harmful to marine organisms and communities, causing pathological and molecular abnormalities (Evans et al. 2000;Corbett et al. 2014), affecting behaviour and survival (Lenihan 1992), and ultimately have been linked to reduced species diversity (Lenihan et al. 1995). However, a reduction in the accumulation of metals within the tissues of long-lived benthic species was observed at both SB1 and SB3 compared to Negri et al. (2006). These reductions in bioaccumulation along with the high diversity of benthic communities observed at both sites suggest that the current contamination impact is likely to be low and localised around Scott Base, and that Antarctic benthic communities may have a level of resilience to moderate levels of pollution. Whilst more extensive monitoring (including benthic infauna) is needed to explore this, high benthic abundance and diversity (especially of sponges and bryozoans) has recently been recorded close to MCM despite its level of pollution .
Benthic community structure is influenced by a complex interaction of multiple biotic and abiotic factors, such as sea ice characteristics, current flows, sediment type and food availability, as well as species interactions (Dayton et al. 1969;Gutt et al. 1996;Raguá-Gil et al. 2004;Cummings et al. 2018). Differences in land-fast sea ice characteristics were particularly evident within this study and may explain the distinct differences in benthic communities between both SB sites and AH1 (Fig. 5, Table 4; Clark et al. 2013). Thicker sea ice, as observed at both SB sites, attenuates more light, reducing microalgal abundance, as shown by the low sediment chlorophyll a concentrations (< 2 µg/g), and can result in currents providing the dominant food source to benthic communities (Dayton et al. 1986;Norkko et al. 2007;Kim et al. 2019). Comparatively, AH1 was characterised by thinner ice cover, greater light availability and thus sediment chlorophyll a concentrations that were 20-40 times higher (28 µg/g). Food supply is known to affect biodiversity in nearshore Antarctic communities (Norkko et al. 2007;Thrush et al. 2011;Lohrer et al. 2013) and therefore is likely to be a significant factor in the dissimilarity observed in benthic communities between AH1 and the two SB sites. Altered food and climate patterns are capable of driving relatively rapid change at the scale of the entire McMurdo Sound (Dayton et al. 2016;, with AH1 being one of the sites exhibiting remarkable change since 1967. However, as a reference site to understand the anthropogenic impacts around SB, AH1 was not ideal. AH1 was on the opposite flank of Hut Point Peninsula from the SB sites, increasing the possibility of different water masses flowing past these sites, increasing their physical separation, and resulting in different levels of food supply and lower ecological connectivity (for example, dissimilar supply and limited exchange of larvae). Differences in food availability have previously been reported along the eastern and western coasts of McMurdo Sound (Barry et al. 1988;Cummings et al. 2003), attributed to proximity and connectedness to productive open water and sea ice thickness. Whilst the establishment of a reference site closer to SB than AH1 was initially confounded by accessibility and known contamination in front of MCM, future exploration of more suitable control sites can now be supported by additional knowledge of current patterns (Fig. 3) and proven utilisation of the ROV which could be used for preliminary site investigations. Our study clearly highlights the influence of SB activities over past years on the benthos and sediments in the adjacent marine area. The sites we have established and sampled (SB1 and SB3), in combination with the historical information available from the 2002 investigations by Negri et al. (2004;, are a valuable basis on which to monitor and detect future changes, and particularly those associated with the current Scott Base rebuild (ANZ 2023). Additionally, the orthomosaic images generated from the ROV footage will enable detailed comparisons of the wider area surveyed over time (Piazza et al. 2019). The infrastructure and protocols we have established will enable regular routine monitoring at future dates using combinations of diver-collected samples and remotely operated cameras, filling a critical need for time-series data gathering for detecting long-term trends and interactions with climate drivers.
Acknowledgements This work was funded by Antarctica New Zealand. Manuscript preparation was supported by the New Zealand Government's Strategic Science Investment Fund (SSIF) to the National Institute for Water & Atmospheric Research (NIWA). Thank you to Vonda Cummings for reviewing a draft of the manuscript, and two anonymous reviewers.
Author contribution AML was principal investigator, designed the study and compiled and analysed the data. AML, PM, RB, DB and LWT conducted the field work. PM and LWT collected and processed the video data. DB collected and BG processed the oceanographic data. LWT processed photogrammetry data. SM and AML wrote the manuscript. All authors read and approved the manuscript.

Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author upon request.

Competing interests
The authors have no competing interests to declare that are relevant to the content of this article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.