The legacy of industrial pollution in estuarine sediments: spatial and temporal variability implications for ecosystem stress

The direct impacts of anthropogenic pollution are widely known public and environmental health concerns, and details on the indirect impact of these are starting to emerge, for example affecting the environmental microbiome. Anthropogenic activities throughout history with associated pollution burdens are notable contributors. Focusing on the historically heavily industrialised River Clyde, Scotland, we investigate spatial and temporal contributions to stressful/hostile environments using a geochemical framework, e.g. pH, EC, total organic carbon and potentially toxic elements: As, Co, Cr, Cu, Ni, Pb and Zn and enrichment indicators. With regular breaches of the sediment quality standards in the estuarine system we focused on PTE correlations instead. Multivariate statistical analysis (principle component analysis) identifies two dominant components, PC1: As, Cr, Cu, Pb and Zn, as well as PC2: Ni, Co and total organic carbon. Our assessment confirms hot spots in the Clyde Estuary indicative of localised inputs. In addition, there are sites with high variability indicative of excessive mixing. We demonstrate that industrialised areas are dynamic environmental sites dependant on historical anthropogenic activity with short-scale variation. This work supports the development of ‘contamination’ mapping to enable an assessment of the impact of historical anthropogenic pollution, identifying specific ‘stressors’ that can impact the microbiome, neglecting in estuarine recovery dynamics and potentially supporting the emergence of antimicrobial resistance in the environment.


Introduction
Industrial activities, whether contemporary or historical, have often occurred globally along major watercourses, particularly in estuaries where transport hubs and habitation sites have developed over civilisation. These features are also naturally chemically and physically important at the interface between terrestrial and marine systems where high flux of natural geochemical transport occurs, stimulating biodiversity and supporting critical ecosystems. The consequence has often impaired water quality, habitat loss and diminished ecosystem services, which result in deleterious changes in ecosystem structure with risks to human and aquaculture health. These changes could result in toxicologically stressful environments that, in turn, could be related to changes in microbiome and/or resistance-e.g. development of antimicrobial resistance (AMR) (Rodgers et al. 2018). Recovery of estuaries from long-term anthropogenic disruption is understood in some detail for higher trophic levels and within that context shows great variation from months to decades depending on biological species or system components (Borja et al. 2010). Very little understanding exists on the microbiome response to pollution and its impact on other characteristics.
The Clyde Estuary, Scotland, is an optimal example site to determine environmental 'stress', historically being one of the most contaminated estuarine environments in the UK (Balls et al. 1997). Through its historical development in Industrial Revolution, trade and wider urban development has seen massive anthropogenic disruption over relatively short spatial scales. Despite the recent improvement in water column quality and positive response of higher trophic levels (Whitfield and Elliot 2002), sediment contamination is retained and slower to recover (Borja et al. 2010). Consequently, an extensive list of historical pollutants has been thoroughly studied Edgar et al. 2003;Hursthouse 2001;Fernandes et al. 2009;Lass-Evans et al. 2012;Rowan 2002). The River Clyde and its tributaries have been important vectors and sink for contaminated waste from succession of industries located in the conurbation of Glasgow and along the banks of the estuary during the industrial development phase in Western Europe in the nineteenth and twentieth centuries. Consequently, the river has received pollution from the onset of the Industrial Revolution (eighteenth century) to the present day, changing with urban and industrial development in its catchment (Vane et al. 2007(Vane et al. , 2010Edgar et al. 2006). Metals such as As, Co, Cr, Cu, Ni, Pb and Zn (and other potentially toxic elements-PTEs) can be linked to anthropogenic and natural inputs (Caccia et al. 2003;Rodriguez-Iruretagoiena et al. 2016;Birch et al. 2015;Lenart-Boroń and Boroń 2014) which are now managed through a range of regulatory instruments, e.g. Water Framework Directive (2000/60/EC) (Cuculić et al. 2009;Khan et al. 2017;Larrose et al. 2010). The bioavailability and environmental influence of these pollutants depend on many processes such as (i) mobilisation in interstitial water, (ii) chemical species, (iii) transformation, e.g. methylation, (iv) sediment composition, e.g. oxides of Fe and organic matter, (v) competition with other metals for complexation, (vi) spatial distribution and (vii) influence of bioturbation, salinity, redox or pH on these processes (Bryan and Langston 1992;Weil and Brady 2016;Peng et al. 2008;Konhausera et al. 2002;Akcil et al. 2015;Zaaboub et al. 2015;Caccia et al. 2003;Berner and Berner 2012;Petit et al. 2015).
Being able to identify key correlations between physicochemical and PTE parameters is crucial in order to identify 'stressors' within the local environment. With past discharges to the River Clyde being reworked by estuarine processes and exposure conditions having a dynamic temporal context, there is a unique ecosystem to test the prevalence of pollutionrelated impacts. For example, the link between AMR and metal content which has previously been identified as a promoter of genetic dissemination by cross and/or co-resistance (Ashbolt et al. 2013;Baker-Austin et al. 2006;Berg et al. 2010;Seiler and Berendonk 2012;Perry and Wright 2013;Martinez et al. 2009;Knapp et al. 2017).
In this study, we explore metal levels and other geochemical characteristics along the River Clyde, Glasgow (Scotland), with an aim to support our baseline hypothesis that specific 'stressors' within the environment contribute risks to public health. Through identification of pollution conditions targeting locations with demonstrable exceedances of toxicity screening levels, we report the statistical significance of physicochemical parameters (pH, Ec, OM % and total organic carbon (TOC)) and key PTEs: As, Co, Cr, Cu, Ni, Pb and Zn within an industrialised estuary to identify key relationships (components). This research acts as a precursor to a larger study which will incorporate and correlate these results with microbial data.

Sample collection
Sediment samples were collected from 19 sites along the inner Clyde Estuary (Fig. 1) during sampling campaigns in Autumn 2017 to Spring 2018, to a maximum depth of 30 cm (handheld coring device; AMS, Inc.). Three cores were collected at each site: one from the high tide line, one from the low tidal front and one-half way between the two; each core was approximately 5 m apart. Four exceptions: Rothesay Dock, King George V Dock, White Cart Estuary and the Kelvin Confluence, were grab samples collected from the centre of the river. Samples were returned to the laboratory, and each core was subdivided into 10 cm sections, before homogenisation, with subsamples taken for microbiological analysis. The remaining sample was air-dried until constant weight and then sieved \ 500 lm for normalisation. This provided nine subsamples per location.

Physical and chemical properties
Organic matter (OM) was determined using loss on ignition with total organic carbon (TOC) using van Bemmelen factor (Dean 1974;Salehi et al. 2011); pH and electrical conductivity were measured using calibrated Mettler Toledo FiveGo meters (Corwin and Rhoades 1982;Van Reeuwijk 2002). Moisture content was determined using a Mettler Toledo HB43-5 m. All analysis was carried out in replicate (n = 5).

Elemental analysis
Replicate aliquots (n = 5) of 0.25 g ± 0.05 g of sample were treated with 5 ml 30% H 2 O 2 overnight to remove organic matter in digestion cups (Environmental Express, USA) and then digested using aqua regia at 95°C in a hot block until 5 ml in volume. They were then made up to 50 ml using UHP water and filtered in situ using 0.45-lm PTFE filter (FilterMate, Environmental Express, USA). A Certified Reference Material (CRM) was included in the digestion method (BCR667 Estuarine Sediment BCR Ò ; Sigma Aldrich).
Elements were determined by inductively coupled plasma mass spectrometry (ICP-MS) (Thermo X-Series II). A calibration series of multi-element standard (BCR320R-40G-channel sediment) was determined regularly (every 60 samples ? linearity r 2 [ 0.99) with indium used as the internal standard. Each determination was made in quadruplicate, and dilutions were made if necessary. Isotopes used for measurement were 75 As, 59 Co, 52 Cr, 63/65 Cu, 60 Ni, 206/207/208 Pb and 66/67 Zn. Results were subjected to the geoaccumulation index in order to gauge the degree of anthropogenic influence.

Results
Limits of detection were (mg/kg, dry weight): As: 0.10, Co: 0.012, Cr: 0.024, Cu: 0.082, Ni: 0.056, Pb: 0.01 and Zn 0.45 mg/kg with recoveries between 85 and 110%. Physical properties and PTE concentrations are given in Table 1, with comparable data for sediment quality criteria shown in Table 2. Statistical analysis was carried out using SPSS, with data tested for normal distribution using Kolmogorov-Smirnov and Shapiro-Wilk tests; PTE concentrations were logtransformed for analysis of variance (ANOVA) and principal component analysis (PCA).

Geoaccumulation index (I geo )
To provide further indication of combined contamination status, the I geo was calculated using the methods of Muller (1979) and Abrahim and Parker (2008): where C n is the concentration of the metal of interest in sediment sample, B n is the background level of the same metal, and 1.5 is the factor to minimise the effect of possible variations in the background (Al-Haidarey et al. 2010;Okedeyi et al. 2014;Hasan et al. 2013).
And evaluated against seven degrees of metal pollution in terms of seven contamination classes (see Fig. 2), with \ 0 demonstrating no contamination and [ 5 extremely contaminated.
IGO data show that all sediment samples (including at varying depths) are contaminated with low to extreme levels of contamination for all PTEs. However, this varied widely between PTE and distance/ depth.
As and Co show low levels of contamination along the hole of the Clyde historically, i.e. with depth.
Chromium predominantly shows an increase in contamination throughout history, whereas PTEs Cu, Ni, Pb and Zn show no discernible trend. Industrial trends throughout history, however, cannot be discounted without radiometric dating. Zn is seen to be the most extremely contaminated with only moderate levels seen in Port Glasgow. This site shows the lowest level contamination trend for all PTEs. There is varying evidence, which suggests both an accumulation in PTEs over time and an insight where contamination may have occurred. We can see that the high levels of contamination are seen between 25 and 33 km from the weir and again 16.44 km to the weir. These include GNK, CAR, KBP and DBW at the estuaries widest and again BSP, ESK, NSI, WCB and CUN, demonstrating zonal relationships (Fig. 1).
Analysis of variance (ANOVA) with Kruskal-Wallis testing was used to determine whether there were any significant differences in PTE concentration with sample location, depth and distance from the high tide mark (core number) (Table 3). There are significant differences (p \ 0.001) between the locations for all PTEs, but there is no difference in the distance from the high tide line, i.e. tidal influence does not significantly affect PTE concentration. With regard to depth, only Pb showed significant difference with sample depth.
Multiple correlation analysis (Table 4) shows a significant positive correlation for most variables with the exception of core number, suggesting that the distance to the high tide level does not influence the PTE concentrations. Note the negative correlation between location and pH (highlighting the changing salinity of estuarine sediments) and depth with Pb and Zn.
PCA was carried out using the correlation matrix and varimax rotation to elucidate the principal components (Alexander et al. 2008). Four principal components (PCs) with an eigenvalue [ 1 were returned (Table 5) and moisture content was not included due the nature of the matrix, i.e. at or near saturation. The dominant factors (41% of total variance) were anthropogenic PTEs: As, Cr, Cu, Pb and Zn which are typical 'urban' elements ) which were all positively related to each other. The remaining three PCs account for 41.9% of the total variance: PC2 shows relationships between Co, Ni and total organic content (19.6%), PC3 is pH (11.8%), and PC4 is conductivity (10.5%).  TEL threshold effect level, PEL probable effect level a G-Baseline refers to the median concentrations (mg/kg) measured for European geochemical stream sediment baseline maps (Ccme 2001) Interestingly, Ni exceeds the G-base, threshold effect limit and probable effect limits (Table 2); this potentially highlights the importance of their roles as essential nutrients.

Discussion
PTE concentrations (Table 1) were compared to a number of previous studies, including the European Geochemical Baseline, Clyde waters and estuary sediments (Jones et al. 2018;Krom et al. 2009), the IUGS/IAGC Global Geochemical Baseline and the FOREGS (Forum of European Geological Survey) (Salminen et al. 2005). The Global Geochemical  Baselines are the median levels of PTEs we expect to see in stream sediments. FOREGS produced a standardised method of sampling, chemical analysis and data management across 26 countries to produce national baseline datasets for comparable measures. These do not regulate or highlight parameter impact; rather, it allows a comparison with broad expected values for terrestrial surface materials (Table 2) (Salminen et al. 2005). The comparison showed that values from this study exceeded the geochemical baseline values and were more comparable to the Clyde Estuary baseline study results, suggesting that levels exceeding the baseline are likely affected by anthropogenic inputs. As previously reported, estuarine sediments in industrialised regions create environmentally stressful conditions for the sediment microbiome (Rodgers et al. 2018); we therefore focussed on 'stressors' and concentrations above the threshold effect level (TEL) and probable effect level (PEL). The TEL represents the concentration below which adverse effects are unlikely to occur, whereas the PEL defines the level in which adverse effects are expected to occur (Macdonald et al. 2000).
It can be seen in Table 1 that As, Co and Cr median values are within the expected geochemical baseline, while Cu, Ni, Pb and Zn are above typical levels. The physicochemical parameters (pH and TOC) show that pH sits within typical baseline values, whereas organic carbon content varies greatly. Water content, conductivity and organic matter show a vast variety in concentrations and vary depending on the sampling site.
Newshot Island, Rothesay Dock, King George V Dock and the Kelvin Confluence exceed the TEL for all, if not most, of the PTEs. This could be an artefact of the different sampling techniques employed at the latter three; these were grab samples collected from the centre of the river rather than within the tidal area; or it could be due to the heavy industrialised presence at these locations. Newshot Island is on the opposite bank to one of Scotland's largest wastewater treatment works which is a potential source of pollution in this area (Rodriguez-Iruretagoiena et al. 2016;Edgar et al. 2003). The PEL is exceeded at a number of locations: Greenock (Cr, Cu and Pb), Dumbarton Castle (Cr), Newshot Island (As, Pb), Rothesay Dock (As, Cr, Ni, Pb), King George V Dock (Cr, Ni, Zn) and Kelvin Confluence (Ni).
There is no spatial trend in As distribution; however, there are a number of hot spots, e.g. Newshot Island and Rothesay Dock where levels exceed the PEL. Historically, As was used to treat wood at timber and shipyards, both of which are known to have been near these locations (Clydewaterfront 2014). Elevated Cr levels at Dumbarton Castle (as well as a I geo level of 'extremely polluted'), Rothesay Dock and King George V Dock may also be linked to shipyards and timber yards as well as leather tanning (Muirhead 2014), particularly associated with known historic pollution levels in the Vale of Leven. Cu is predominantly below the TEL for most sample sites with little variation; the greatest variability is seen at Greenock and Newshot Island, with the former showing up to 13% variability, half of which exceed the PEL.
Elevated levels of Ni at the same site (average of 33.9 mg/kg) can be linked to Cu-Ni alloys typically used in ship building. With Ni there are also hot spots in areas of similar industrial heritage: Rothesay Dock, White Cart Bridge, King George V Dock and Kelvin Confluence. It is also worth noting that the Ni threshold effect level is exceeded at all but five locations, in conjunction with being moderatelyextremely polluted (I Geo , Fig. 2) which will impact the microbiome. There is no clear trend with Co, but there is high variability in concentrations (2.3-15 mg/ kg), with only the higher [ 11.2 mg/kg to be considered above the natural baseline. The highest concentration of Zn was found at King George V Dock, and this was the only location to exceed the PEL. Zn, however, has the most polluted Igeo factor suggesting a strong human influence. There appears to be no pattern in PTE concentration as the locations progress from saline conditions (Helensburgh) to freshwater conditions (Kelvin Confluence). Our control site at Cuningar Loop exceeded the TEL for Ni and three sample sites did not exceed the TEL for any PTE: Helensburgh, Cardross and the White Cart Estuary. These locations are all areas of low-current industrial activity compared to other sites. Changes in concentration are also reflected in the minimum and maximum values, with the maximum values exceeding regulatory and typical values. A comparison with a recently published study (Jones et al. 2018) which reported on summary data for samples collected in campaigns nearly two decades earlier from the Firth of Clyde (saline conditions) found that their concentrations are higher than our studies; however, our most comparable site (Greenock) exceeds the TEL for Ni and Zn and the PEL for Cr, Cu and Pb. Greenock also has the highest concentrations of Cr and Cu in our study (Table 1). A major difference between our study and Jones et al. (2018) is in relation to chromium; Jones et al. reported the majority of locations exceeded the PEL, whereas although twelve of our sites exceeded, the TEL only four exceeded the PEL.
As the Clyde Estuary is an 'open system', it is possible that there has been natural attenuation of historical pollution in the upper sediment column, particularly as Jones et al. (2018) collected samples from a greater depth. A greater column depth would not only reflect an earlier period in time, but also reflect more clay condition and possibly more 'stable' conditions, i.e. sediments that are not as affected by tidal changes. In our study the uppermost sediments were 'jelly-like' in their composition and are constantly being washed which could reduce higher PTE concentrations in areas where there are reduced industrial inputs. This therefore has implications for a potential increase in antimicrobial resistance due to exposure of historical pollution at depth (Rodgers et al. 2018). The time difference between our campaign and previous reports, particularly for Cr, shows the dynamic nature of the system with regard to complexity and elemental movement. It is widely reported that there are localised upstream sources of Cr contamination (Whalley et al. 1999); in the last decade, extensive engineering associated with transport infrastructure and international athletics events has transformed some of the most significant sources of Cr pollution with subsequent knock on effect downstream.
PCA analysis highlighted two prominent groupings relating to anthropogenic and natural sources; Component 2 (Co, Ni and TOC) we suggest is the natural and biological requirements within ecosystems because they are less commonly a consequence of anthropogenic causes. However, Co and Ni are typically distributed uniformly within sediments, their concentration increases with depth, and certain locations are likely to be as industrial depositions have accumulated over time (Leyssens et al. 2017;Cempel and Nikel 2006).
These results suggest that changes in salinity and proximity to high tide level have little impact on concentration; however, there are a number of pollution hot spots in the Clyde Estuary (Greenock, Newshot Island, Rothesay Dock, King George V Dock and the Kelvin Confluence) suggestive of sitespecific inputs. In addition, there are sites with high variability which can be attributed to many factors including sediment mixing with historical inputs. Observations of the potential impact of upstream interventions, in particular on the input of chromium, highlight the dynamic nature of the system and still the system acts as a net pollution sink.

Conclusion
Understanding the distribution of PTEs and geochemical factors in an industrialised estuary is important to elucidate the contributions to stressful/hostile environments. We see levels of PTEs often exceed TEL/ PEL levels at various locations and depths along the Clyde, although they are starting to reduce in comparison with historical levels. Our PCA analysis identified four components with PC1 and PC2 accounting for the highest variance and associations: As, Cr, Cu, Pb and Zn as well as Ni and Co with TOC. Compared to extensive assessment of sediment quality in the region over more than a decade ago, our data show there has been a gradual reduction in sediment contamination levels, and for example impact of upstream source control (extensive remediation of Cr contamination in Glasgow). However, significant exceedances of sediment quality standards still occur in the estuary and are spatially variable.
It is evident that the estuarine sediments are still a reservoir of contaminants resulting in persistent exposure and disruption to the microbial consortia. Further research is required in order to determine whether these components can be identified as 'stressors' within an environmental microbiome. If identified correctly, does this subsequently exacerbate the prevalence of genetic mutations and more specifically antimicrobial resistance? The 'contamination map' created within this study provides the basis to enable a detailed study of the impact on the microbiome and causes the emergence of antimicrobial resistance in the environment.