Organic-carbon-rich sediments: benthic foraminifera as bio-indicators of depositional environments

Fjords have been described as hotspots for carbon burial, potentially playing a key role within the carbon cycle as climate regulators over multiple timescales. Nevertheless, little is known about the long-term fate of the carbon that may become stored in fjordic sediments. One of the main reasons for this knowledge gap is that carbon arriving on the seafloor is prone to post-depositional degradation, posing a great challenge when trying to discriminate between an actual change in the carbon deposition rate and postdepositional carbon loss. In this study, we evaluate the use of modern benthic foraminifera as bio-indicators of organic carbon content in six voes (fjords) on the west coast of Shetland. Benthic foraminifera are known to be sensitive to changes in organic carbon content in the sediments, and changes in their assemblage composition therefore reflect synchronous variations in the quantity and quality of carbon reaching the seafloor. We identified four environments based on the relationship between benthic foraminiferal assemblages and organic carbon content in the sediments: (1) land-locked regions influenced by riverine and/or freshwater inputs of organic matter, namely the head of fjords with a restricted geomorphology; (2) stressed environments with a heavily stratified water column and sediments rich in organic matter of low nutritional value; (3) depositional environments with moderate organic content and mild or episodic current activity; and (4) marginal to coastal settings with low organic content, such as fjords with an unrestricted geomorphology. We conclude that foraminifera potentially provide a tool to disentangle primary organic carbon signals from post-depositional degradation and loss of organic carbon because of their environmental sensitivity and high preservation potential in the sedimentary record.


Introduction
Globally, fjords sequester about 18 Mt of organic carbon (OC) annually in their sediments, partly acting to reduce the 25 remineralisation of OC into the environment and thus buffering the release of carbon dioxide (CO2) to the atmosphere (e.g. Smith et al. 2015;Cui et al., 2016). A recent study by  estimated that the sediments of Scotland's 111 major fjords hold 295 ± 52 Mt OC. In Shetland, these fjordic systems are different from their mainland counterparts, notably in that they are shallower and less glaciated (Edwards and Sharples, 1986) and peat or peaty soils dominate their catchments (Soil Survey of Scotland, 1981). Erosion of terrestrial organic matter (OM) is known to be a major source of OC in fjordic sediments (Cui et al., 2017) and Scottish peatlands hold 1620 Mt of OC (Chapman et al., 2009). Shetland peatlands cover over one third of the archipelago's total area and their erosion accounts for a loss of 1 to 4 cm yr -1 from peat surfaces (Birnie, 1993). 5 Eroding peatlands may therefore hold the potential to contribute significantly to OC burial and accumulation in Shetland's voes and to the overall coastal OC store of this region. Nevertheless, little is known about the long-term fate of the OC that may accumulate and become stored in these sediments. One of the main reasons for this knowledge gap is that OC is prone to post-depositional modification and degradation. It is therefore difficult to discern whether a down-core decrease in OC content is due to an actual change in carbon deposition over time or simply to this post depositional loss (Hülse et al., 2017;Stolpovsky 10 et al., 2017).
Benthic foraminiferal assemblages have long been used as indicators of environmental change and to reconstruct the intrinsic characteristics of marine ecosystems. A recent study (Duffield et al., 2017) explored the relationship between foraminiferal assemblages and variations in OC content in the surficial and recent sediments of Norwegian fjords, highlighting temporal trends in OC fluxes and ecological quality status. In this study, we evaluate the use of modern benthic foraminifera as bio-15 indicators of OC content in six voes (fjords) on the west coast of Shetland whose sedimentation dynamics are likely influenced by the erosion of OC rich peat from the surrounding catchments (Fig. 1). As benthic foraminifera are sensitive to changes in OC content in the sediments (Alve et al., 2016 and references therein), changes in their assemblage composition should reflect synchronous variations in the amount of OC reaching the seafloor. Given the high preservation potential of their tests in marine sediments (e.g. Murray, 2003;Thomas et al., 1995), foraminifera potentially offer a means to discriminate primary OC signals 20 from post-depositional OC degradation and may provide a tool to validate studies and models reconstructing past changes in the carbon cycle (Hülse et al., 2017;Stolpovsky et al., 2017). Nevertheless, the relationship between OC and foraminiferal assemblages is not straightforward and warrants further investigation.
Early studies developed conceptual models explaining benthic foraminiferal microhabitat preferences in deep sea settings based on changes in the organic flux reaching the seafloor and bottom water oxygenation (Jorissen et al., 1995;Van Der Zwaan 25 et al., 1999). These authors found that in well-oxygenated bottom waters, food is usually the limiting factor controlling the distribution of benthic foraminifera, whereas under eutrophic conditions, bottom water oxygenation may become the dominant control on the composition of foraminiferal assemblages (Jorissen et al., 1995;Van Der Zwaan et al., 1999). Typically, food availability is strongly connected to OC fluxes and the transport of organic matter to the seafloor. In this context, it becomes crucial to understand not only food availability in terms of the quantity of OM reaching the seafloor, but also to understand 30 the quality of this food source. For example, refractory organic matter (ROM), often derived from reworked terrestrial sources, is known to be of low nutritional value for benthic foraminifera compared to the more bio-available form of labile organic matter (LOM) (e.g. Fontanier et al., 2008;Jorissen et al., 1998). It should also be noted that OC accumulation rate and bottom water oxygen concentrations are often related, with low oxygen concentrations reducing the reworking/remineralisation of OM and thus favouring its preservation in marine sediments. In certain environments, this implies that similar OC fluxes could result in different OC accumulation rates depending on bottom water oxygenation, the latter defining the composition of benthic foraminiferal assemblages (Van Der Zwaan et al., 1999).
To investigate the relationships between sedimentary OC in six west Shetland voes and the associated changes in benthic foraminiferal assemblages, this study combines geochemical and micropalaeontological techniques to: 1) Fingerprint the 5 source (terrestrial vs. marine) and quality (refractory vs. labile) of organic matter and the form (organic vs. inorganic) of sedimentary carbon. 2) Establish benthic foraminiferal distribution in Shetland's voes from recent surficial sediments. 3) Investigate the use of benthic foraminifera as bio-indicators of OC content in coastal sediments and their potential for palaeo-OC reconstruction purposes.

Area of study and sampling strategy
Shetland is a subarctic archipelago covering a total area of 1466 km 2 , lying about 170 km north of mainland Scotland, UK ( Fig. 1). Several voes dissect Shetland's coastline, providing an extended route for the transport of OC rich material from land to sea (Cui et al., 2016). The current study focuses on six voes on Shetland's west coast: Clift Sound, Sand Sound, Olna Firth, Aith Voe, Busta Voe and Vaila Sound (Fig. 1, Table 1). These six voes are characterised by different geomorphologies, 15 freshwater inputs, proximity to the open sea and intensity of local currents; together these characteristics determine the amount of carbon deposited in the sediments. Clift Sound lies south of Shetland, the islands of Trondra and East Burra to the west, has steep-sides and a relatively unrestricted geomorphology (Fig. 1). Two main types of soil are present in the catchment of Clift Sound: humus-iron podzols cover the islands of East Burra and Trondra, while peaty soils dominate the eastern coast of Clift Sound ( Supplementary Fig. 1). Winds are generally from the southwest, can be channelled and result in intensified currents 20 due to the geometry of the sound and the surrounding highlands. Tidal energy can also be channelled into more vigorous currents (Cefas, 2007a).
Sand Sound is somehow T-shaped, located in the southwest of Shetland and comprises three areas: the head of the voe, the inner basin and the outer basin (Fig. 1). The head of the voe, almost perpendicular to the inner and outer basins, is a sheltered area that receives freshwater from a significant number of rivers draining the surrounding land. Three main types of soil 25 characterise this area: peaty gleys, peaty podzols and peaty rankers ( Supplementary Fig. 1). Here, the voe is very shallow and drying can occur in its eastern and western periphery (Cefas, 2007b(Cefas, , 2008b. In contrast, the rest of the voe is much deeper. The inner basin gently slopes away reaching more than 20 m at its maximum depth; drying can still occur where this basin meets the head of the voe (Cefas, 2007b). The outer basin has steep sides, is 42 m deep at its deepest point and completely open to sea. A shallow sill divides the inner from the outer basins (Cefas, 2007b) 30 Busta Voe, Olna Firth and Aith Voe are part of a major inlet on the southern coastline of St Magnus Bay on the west coast of Shetland (Fig. 1). Olna Firth represents the eastern branch of the inlet and is roughly oriented eastwest, it is the furthest away from the open sea (Fig. 1). A large area of this voe exceed 30 m water depth; the northern coastline gradually slopes into the voe, whereas the south side shows a steeper gradient. Olna Firth is classified as micro-tidal and due to the generally low level of energy in this system, stratification may occur, especially during summer time (Cefas, 2013). 5 Aith Voe is the southern component of the inlet, with a northsouth orientation (Fig. 1). Due to the hilly landscape (up to 100 m) surrounding the voe, Aith Voe is exposed to winds from the north, which can significantly alter local circulation and prevent stratification of the water column. Two types of soil surround Aith voe: peaty and organic soils ( Supplementary Fig. 1); few streams drain these soils, with the largest of them discharging on the east coast of Aith Voe (Cefas, 2010).
Busta Voe lies in the northern part of the inlet and is oriented northsouth, sheltered and with a maximum water depth of 39 10 m (Fig. 1). Three main types of soil surround this area: peaty gleys, peaty podzols and peaty rankers ( Supplementary Fig. 1); two rivers drain a relatively small catchment area into Busta Voe. Tidal flow is weak in the voe and wind generated currents are predominant; stratification of the water column especially during warm periods can occur (Cefas, 2008a).
Vaila Sound lies on the western coast of Shetland, with the island of Vaila to the southwest and the isle of Linga in the middle ( Fig. 1). Three types of soil dominate the catchment of this voe: peaty gleys, peaty podzols and peaty rankers (Supplementary 15 Fig. 1). Most streams discharge in the northern part of the sound; hence, locations in this area might become more affected by terrigenous inputs. Additionally, Linga offers protection from wind and currents, which could facilitate localised accumulation of terrigenous material. Shelter from strong winds is also provided by the island of Vaila to the south. Here the tidal range is small and the associated energy weak; however, Vaila Sound is connected to the Atlantic Ocean west and east of the island of Vaila (Cefas, 2009). 20 The MV Moder Dy 2015 cruise in west Shetland surveyed the six voes in August 2015 (Table 1). Marine surface sediments were sampled using a Duncan and Associates Van Veen grab with a sampling area of 0.1 m 2 . Twenty-three surface sediment samples were obtained by scraping the top layer (~ 1 cm thick) of each grab with a domestic spoon, samples were then stored in a cold-box in sealed plastic bottles.
An earlier field survey of Shetland voes carried out in August 2009 measured bottom water temperature (BWT), salinity (BWS) 25 and oxygen (O2) at the same locations as this study (Fig. 2); however, no data were collected in Vaila Sound. Oceanographic measurements were made using a SeaBird SBE 19 plus conductivity-temperature-depth (CTD) profiler with a dissolved oxygen probe. Bottom water temperature varies between 14.14 and 12.20 °C being on average 12.74 °C (Fig. 2a). Bottom water salinity remains almost constant across these voes ranging between 35 and 35.2 p.s.u. with the exception of MD15-06 where BWS reached a minimum of 34.7 p.s.u (Fig. 2b). Similarly, bottom waters are mostly well-oxygenated in all the voes 30 (O2 > 5mg l -1 ) with the exception of Olna Firth where lower O2 were observed (Fig. 2c).

Particle size measurement
Particle size analyses were performed on a sub-set of surface sediments using a Coulter LS230 particle size analyser to quantify the volume (%) of the various grain sizes in each sample. We identified three categories based on the International Organisation for Standardisation (ISO) scale: 'Clay' includes particles < 2 µm, 'Silt' particles between 2-63 µm, and 'Sand' particles > 63 µm (upper limit of 2000µm). Before analysis, sediment samples were digested using 30 % hydrogen peroxide to remove 5 organic matter, and 10 % hydrochloric acid to remove carbonate, following a similar methodology to that outlined by Austin and Evans (2010). Digestions were carried out to obtain an insoluble residue useful for comparing sediments having very different carbonate and organic matter content, which otherwise might skew the measurements because of in situ processes.

Loss on ignition (LOI) analysis
Loss on ignition measurements were carried out to quantify the percentage of total organic matter (TOM), and the respective 10 amounts of labile (LOM) and refractory organic matter (ROM). The initial combustion temperature was set to 250˚ C, as significant mass loss of LOM has been recorded at this temperature (Mook and Hoskin, 1982). Sediments were successively heated to 550˚ C to quantify ROM, as refractory terrestrial and aquatic OM will likely burn off at this temperature (Kristensen, 1990). About 1 g of dried sediment was precisely weighted into crucibles of known weight for each surface sample (M0). Where Ms is the sample weight minus the crucible weight. Replicate samples were measure at stations MD15-01A and 01B 20 and MD15-05A and 05B which resulted in a mean relative error of ± 0.07 % for LOM, ± 0.06 % for ROM and ± 0.03 % for TOM, pointing to very good data reproducibility and representation of local conditions.

Bulk carbon elemental & stable isotope analyses
Surface sediment samples were analysed to determine bulk elemental C and stable isotope ( 13 C) values. Each sediment sample was freeze dried and homogenized and approximately 12 mg of milled sediment were placed into silver capsules and a further 25 10 mg into tin capsules.
The samples encapsulated in silver underwent an acid fumigation step (Harris et al., 2001) to remove inorganic carbon (IC).
After drying for 24 hours at 40 °C, both OC and  13 C were measured using an elemental analyser coupled to an isotope ratio mass spectrometer (IRMS) at the NERC Life Science Mass Spectrometer Facility (Lancaster, UK). The standard deviation of  13 C triplicate measurements was 0.07 ‰ and  13 C values are reported in standard delta notation relative to Vienna Pee Dee Belemnite (VPDB).
The samples in tin capsules were analysed for total carbon (TC) using a Elementar Elemental Analyser (EA) at the School of Geography and Sustainable Development, University of St Andrews (Verardo et al., 1990). Precision of the analysis is calculated based on repeat measurements of standard reference material B2178 (medium Organic content standard from 5 Elemental Microanalysis, UK) with C = 0.08%.
These results were combined together to calculate the quantity of IC in each sample; the percentage of OC was subtracted from the percentage of TC to obtain the percentage of IC.

Binary mixing model
To discriminate between marine-sourced and terrestrially-sourced OC, we used a two end-members (binary) mixing model 10 assuming a two-point source of OC, marine and terrestrial, following Thornton and McManus (1994). We used  13 C measurements in the sediment samples as tracer of the source of OC in west Shetland voes and calculated the fraction of terrestrially-sourced OC (OCterr) as follows: OCterr = ( 13 Cmar- 13 Cterr) / ( 13 Cmar- 13 Csample) OCterr + OCmar = 1 15 where the end-member value for marine-sourced OC ( 13 Cmar) was taken from Smeaton and Austin (2017) based on phytoplankton, zooplankton, macroalgae, and benthic microalgae carbon stable isotope values, while the end-member for terrestrially sourced OC ( 13 Cterr) was based upon the work of Thornton et al. (2015) who studied the distribution of carbon stable isotope in Scotland's topsoil; OCmar indicates the fraction of marine-derived OC. We assumed no isotope discrimination between the source and the surface sediment samples. 20

Benthic foraminiferal counts and statistics
Aliquots of surface sediments were stained with Rose Bengal to allow the identification of living foraminifera, as detailed in Schönfeld et al. (2012). However, it should be noted that foraminiferal counts are 'total' (live + dead) because the main objective of this study is to provide a tool for the interpretation of fossil foraminiferal assemblages and their relationship with changes in OM and OC content in sediments over time (e.g.: Conradsen, 1993). Our priority is to understand how benthic 25 foraminifera respond to prevailing (long-term) environmental conditions, rather than to seasonal variability. Additionally, scraping the top layer (~ 1 cm thick) of each grab with a domestic spoon may lead to a potential underestimation of living fauna if mixing of the top layer occurred during sampling. In this scenario, the number of living foraminifera at the sediment surface will be "diluted" due to the presence of dead/fossil specimens from deeper sediments. Having said this, a recent study by (Rillo et al., 2019) reported that historical sediment samples collected in a way that could have caused disturbance of the 30 sediment surface (sounding and dredge) are still representative of surface conditions and their foraminiferal assemblages can be used to reconstruct environmental changes reliably. At two sites, MD15-01 and MD15-05, replicate samples were taken (-01A and -01B, -05A and -05B, respectively) to check measurement reproducibility and seafloor heterogeneity.
As part of an unpublished work, the rose Bengal stained surface sediments were wet sieved over a 63 m mesh; the residues oven-dried at < 60 °C, then weighed and dry-sieved into 63-150 m and >150 m size fractions. As part of an associated work aiming to understand the influence of different size fractions on the composition of benthic foraminiferal assemblages (Lo 5 Giudice Cappelli et al., in review in Frontiers), both size fractions were analysed independently, and the results demonstrated that a more holistic picture of environmental change is obtained when benthic foraminiferal assemblages of both size fractions are compiled together. This improves the representation of the entire assemblage and results in more statistically robust environmental reconstructions. Similar findings were also reached by Weinkauf and Milker (2018 and references therein).
Depending on sample volume, each sample was divided into a number of splits using a standard dry-splitter and, when possible, 10 at least 300 specimens were dry-picked to ensure statistical significance when discussing benthic foraminiferal assemblages' diversity and composition (Schönfeld et al., 2012). Species were identified following Austin (1991) and counted using a tally sheet; notes on tests preservations were made when signs of etching were observed.
Relative benthic foraminiferal abundances were calculated for each sample (Supplementary Table 1). It should be noted that we grouped under the name E. excavatum both forma selseyense and forma clavata despite them being recognised as 15 genetically different species (Darling et al., 2016). It was not possible to consistently identify and discriminate between these two species, hence the use of a common identifier.
The Palaeontological Statistics software package PAST (version 3.16, Hammer et al., 2001) was used to analyse relationships between benthic foraminiferal assemblages through cluster analysis, canonical component analysis (CCA) and non-metric multidimensional scaling (MDS) based on the Bray-Curtis similarity index and relative abundance data. Cluster analysis 20 showed good reproducibility of benthic foraminiferal assemblages ( Fig. 4b) with the two pairs of replicate stations exhibiting assemblage reproducibility of 90% (MD15-01A and 01B) and 94% (MD15-05A and -05B).

Carbon content and source in west Shetland bulk sediments
The percentage of total carbon (TC) in west Shetland voes varies between 4.31 and 14.32% (Fig. 3d), with the highest percentage found at station MD15-05 (Clift Sound) and the lowest at station MD15-09 (Sand Sound). Organic carbon (OC) content ranges between 2.0-8.51%, peaking at station MD15-11 (Olna Firth) and reaching a minimum at site MD15-17 (Busta 10 Voe) ( Fig. 3d; Table 3). The percentage of inorganic carbon (IC) varies between and 0.35-7.51% and is the highest in station MD15-03 (Clift Sound) and the lowest in station MD15-14 (Aith Voe) (Fig. 3d). We excluded from the data set site MD15-15, as this sample contains coal, which is believed to have been shed accidentally from a commercial vessel. In general, there is a tendency for high OC and low IC at the head of the voes and the reverse occurs where the voes connect to the sea (Fig.   3d). 15 In terms of carbon source, OCterr is generally higher at sites close to land than in proximity to the open sea (Fig. 3b). It reaches a peak of 72% at station MD1-06 in Sand Sound, and a minimum of 21% at station MD15-03 (Fig. 3b).

Benthic foraminiferal assemblages in west Shetland voes
We observed a total of 64 species of benthic foraminifera in the 21 stations analysed in this study, including the two replicates samples (Supplementary Table 1 Other common agglutinated species are Haplophragmoides sp. and S. wrightii, with relative abundance ranging from 0.1 % (Vaila Sound) to 2.9 % (Busta Voe), and 0.1% (Olna Firth) and 3.9 % (Clift Sound), respectively (Supplementary Table 1).
The distribution of the remaining agglutinated species is patchy, and their relative abundance varies around the average value of 1 % (Supplementary Table 1 Table 1).
The remaining 49 species are hyaline foraminifera of which 9 species have relative abundances of > 10 % in at least one of the 21 stations studied: Ammonia spp., Bulimina marginata, Buliminella elegantissima, Cibicides spp., Elphidium gerthi, 10 Elphidium margaritaceum, Elphidium excavatum, Rosalina spp. and Stainforthia fusiformis (Table 2). Additionally, of all the hyaline species, 20 are rare having relative abundances < 1 %, while 18 species are frequent having a relative abundances between 1-5 % and only 2 species are common having relative abundances ranging between 5-10 % (Supplementary Table 1 Overall, ten taxa of benthic foraminifera were observed in each one of the six voes and had relative abundances of more than 15 10 % in at least one of the 21 stations studied (Table 2): Ammonia spp., Bulimina marginata, Buliminella elegantissima, Cibicides spp., Eggerelloides scaber, Elphidium gerthi, Elphidium margaritaceum, Elphidium excavatum, Rosalina spp. and Stainforthia fusiformis. These ten taxa were categorised as dominant and imaged with a scanning electron microscope (SEM) at the Scottish Oceans Institute, University of St Andrews (Fig. 5). Altogether, these ten dominant taxa account for more than 70% of the total assemblage at each station (Table 2). In general, the most abundant taxa observed in west Shetland's six voes 20 is the hyaline Cibicides spp. relative abundance, followed by the agglutinated E. scaber relative abundance (Table 2 and Fig.   5).

Statistical analyses of benthic foraminifera data and environmental parameters
We used cluster analysis of environmental parameters (BWT, BWS, O2, % Clay, OC, IC, ROM, LOM,  13 C and WD) based on Euclidean distance and of benthic foraminifera relative abundance based on the Bray-Curtis index to illustrate similarities 25 between locations (Figs. 4a and 4b), and performed CCA on the combined dataset (benthic foraminifera + environmental parameters) to examine the response of benthic foraminiferal assemblages to environmental gradients (Fig. 4c). Additionally, we run non-metric MDS based on the Bray-Curtis similarity index and foraminifera relative abundance data to analyse for similarities between the populations of the 21 studied stations, independently from environmental forcing (Fig. 5). To note that temperature, salinity and oxygen measurements are from 2009 and may not fully represent the conditions in 2015 when the 30 sediment samples were collected.
Based on both CCA and non-metric MDS, we defined the same four groups with 75% confidence; in the CCA plot axis 1 explains 57 % of total variability and axis 2 the 26 % based on the 10 environmental parameters (Figs. 4c and 5).  5). Station MD15-17 in Busta Voe falls just outside the 75% confidence interval defining Group 4; however, cluster analysis 10 ( Fig. 4b) shows ~82% similarity between the assemblage of this site and station MD15-20 and more than the 70% similarity with the rest of the locations belonging to Group 4. Therefore, for the purposes of our discussion, we include station MD15-17 in Group 4.

Particle size distribution in west Shetland voes 15
The particle size data from the voes of Shetland exhibit a large range with the silt and sand fractions, but little variation in the proportion of clay (Fig. 3a). These samples suggest a significant size sorting which is consistent with the tidal and estuarine dynamics of these voes. In Vaila Sound, intensified bottom currents and winnowing of finer sediments owing the influx of Atlantic Ocean waters (Cefas, 2009) appear to control local grain size distribution (Fig. 3a). Here, stations exposed to the open sea and intensified bottom currents are characterised by coarser sediments (MD15-20 and 21), whereas more sheltered and 20 calm locations  are characterised by finer sediments (Fig. 3a). In contrast, sediments in Clift Sound are siltdominated ( Fig. 3a) despite the unrestricted geomorphology ( Fig. 1) and vigorous bottom currents (Cefas, 2007a). It appears that a different mechanism than winnowing of finer sediments by intensified bottom currents drives particle size distribution here. In Clift Sound, the overall low percentage of OM together with the relatively high IC content and  13 C values point to a reduction in terrigenous input as the main driver of particle size distribution (Fig. 3). Conversely, in Sand Sound, an increased 25 terrigenous input appears to be responsible for the coarser sediments found in stations MD15-06, -07 and -09 (Fig. 3). The head of the voe is a sheltered area that receives freshwater from a significant number of rivers draining the surrounding land (Cefas, 2007b(Cefas, , 2008b, resulting in an overall high percentage of OM and OC, low  13 C values (Figs. 3c and 3d) and coarser sediments (Fig. 3a). In contrast, station MD15-08 at the head of Sand Sound is characterised by finer sediments (Fig. 3a) despite the high terrigenous supply at this location (Figs. 3c and 3d). It is possible that, due to the restricted geomorphology of 30 Sand Sound (Fig. 1), weak currents are unable to transport the coarse material to station MD15-08, the deepest station at the head of the voe (11 m vs. 6 m). In Sand Sound inner basin (MD15-10), silty sediments dominate (Fig. 3a). Here, a combination of reduced terrigenous input and/or weaker currents could be driving particle size distribution. Nørgaard-Pedersen et al. (2006) found similar drivers of particle size distribution in Loch Etive, a sea loch in mainland western Scotland. Vigorous bottom currents and/or enhanced terrigenous supply would result in locally coarse sediments, whereas weaker currents and/or limited terrigenous supply would result in finer sediments. Due to the duality of drivers determining particle size distribution in west 5 Shetland voes (bottom currents intensity vs. terrigenous supply), it is challenging to consistently group locations based on granulometry data (Fig. 3a) or identify a relationship between grain size distribution, OM and OC content in the sediments ( Supplementary Fig. 3).

Spatial distribution, quantity, quality and source of organic matter, and carbon in west Shetland voes.
A clear pattern is evident with regard to the quantity, quality and source of organic matter and carbon measured in the six west 10 Shetland voes (Fig. 3). Typically, carbon is present in the sediments as IC and OC and the distribution of these two forms of carbon tends to mirror one another, with low IC and high OC percentages towards land while the reverse occurs close to the open sea (Fig. 3d). Organic matter follows the same pattern of OC distribution regardless of its quality (refractory vs. labile); however, ROM content in the sediments is always higher than LOM (Fig. 3c) owing its resilience to degradation.
About the source of OM and OC, low  13 C and high OCterr values characterise the head of the voes where high ROM and OC 15 are accumulated, indicating that the source of this material is terrestrial (Fig. 3). Terrestrial OM is more stable and resistant to decay than marine OM (Batten, 1996) and it is enriched in 12 C from soils resulting in lower  13 C values than marine OM derived from algae (e.g. Balasse et al., 2005;Ficken et al., 1998;Marconi et al., 2011;Schiener et al., 2014;Schmidt and Gleixner, 2005). Sand Sound is dominated by organic rich and terrestrially sourced material as it receives freshwater and terrestrial OM from the numerous streams draining the surrounding land (Fig. 3). In contrast, Clift and Vaila Sounds are 20 characterised by marine sourced material and high IC content in their sediments (Fig. 3) as the inflow of Atlantic Ocean waters improves the voes' potential for IC storage (Burrows et al., 2017;Smeaton et al., 2016. Olna Firth has the highest LOM percentage as this relatively deep and micro-tidal basin (Cefas, 2013) is prone to water mass stratification which can lower bottom water oxygenation and slow down the reworking/remineralisation of OM favouring the preservation of the more labile material (Fig. 3). 25 A slight deviation from this overall trend was recorded at station MD15-05. We speculate that the unusually high TC content at station MD15-05 at the mouth of Clift Sound (Fig. 3d)

Group 2 -E. scaber, B. marginata and S. fusiformis 15
Eggerelloides scaber, B. marginata and S. fusiformis are species typical of organic rich sediments, stratified waters and possibly low oxygen concentrations (Fig. 4c), often suggesting stressed environments (Alve, 1994Evans et al., 2002;Fontanier et al., 2002;Murray, 1992;Nørgaard-Pedersen et al., 2006). Bulimina marginata (~9%) seems to thrive in regions with high ROM (9.6%) due to its ability to use OM of low nutritional value as food (Klitgaard-Kristensen & Buhl-Mortensen 1999) (Figs. 4c and 5). Stainforthia fusiformis (~21%) is a known opportunistic species which follows high OM gradients at 20 stratification fronts Alve and Murray, 1997;Scott et al., 2003). In general, high organics in the sediments can result in low oxygen concentrations at the seafloor (Fig. 4c) due to eutrophication and increased respiration . Additionally, the production of humic acids as a by-product of organic matter biodegradation may lower the pH (Bauer and Bianchi, 2012), overall making the environment more hostile for foraminifera to calcify, which could result in the dominance of agglutinated species like E. scaber (~31%) (Fig. 5). Etching of Ammonia spp. shells ( Supplementary Fig. 2

) in 25
Olna Firth supports local acidification of bottom waters, which is substantiated by the low IC content (2.6%) of stations MD15-11 and -12 (Fig. 3). Therefore, assemblages characterised by E. scaber, B. marginata and S. fusiformis reflect a stressed environment with a heavily stratified water column -possibly low oxygen levels, and organic rich sediment especially in ROM of low nutritional value. Despite seawater properties (T, S and O2) were measured in 2009, they seem to be representative of seawater characteristic in 2015 when the sediment samples were collected, as the two sets of measurements (seawater and sediments) appear to be in agreement, reflecting long-term of environmental dynamics.

Group 3 -no obvious dominant taxa
The assemblage of Group 3 is characterised by high relative abundance of taxa with affinities for OC (4.9%) and OM (6.3%) like Elphidium (8%), E. scaber (13%) and B. elegantissima (5%) (e.g.: Mendes et al. 2012;Alve & Murray 1999) as well as 5 high energy taxa like Cibicides (22%) and Ammonia (5%) (Scott et al., 2003) (Figs. 4c and 5). The low relative abundance of B. marginata (2%) and S. fusiformis (5%) additionally suggests that OM at these locations has good nutritional value and current activity is present at these sites (reduced stratification, possibly vertical mixing). Therefore, assemblages showing this composition are likely to reflect depositional environments with mild or episodic current activity and moderate organic content.

Group 4 -Cibicides spp. and Rosalina spp. 10
In Group 4, the dominant taxa is Cibicides (~57%) followed by Rosalina (8%) (Fig. 5 A clear relationship is evident between low organics (OC 4.2%, ROM 3.2% LOM 1.5%) (Fig. 4c) and the dominance of Cibicides (Figs. 4c and 5). All locations falling in Group 4 are on the path of Atlantic Ocean waters entering the fjords meaning that this assemblage composition is most likely found in areas dominated by sea processes. 15

Other drivers influencing modern benthic foraminifera distribution
Seawater properties can compete with OC and OM in influencing the distribution of benthic foraminifera across the six voes, especially BWS and O2 concentration may exert a second-order control on the composition of foraminiferal assemblages (Jorissen et al., 1995;Van Der Zwaan et al., 1999). A highly variable BWS generally influences the distribution of stenohaline (strictly marine species) and euryhaline species (able to adapt to a wide range of salinity) along BWS gradients; however, in 20 Shetland, BWS is almost constant across the voes (Fig. 2) having very minimal influence on benthic foraminiferal assemblages (Fig. 4c). Similarly, bottom waters are mostly well-oxygenated ([O2] > 5 mg l -1 ) in all the voes (Fig. 2), meaning that food availability is the main driver of changes in the composition of benthic assemblages (Jorissen et al., 1995;Van Der Zwaan et al., 1999). It should be noted that seawater properties reported in Fig. 2 were measured in 2009 and may differ from the actual conditions in 2015 when sediment samples were collected; nevertheless, our results show that similar dynamics were most 25 likely in place in 2015, as the two sets of measurements compare very well in depicting local environmental dynamics.
Food availability is deeply connected to organic fluxes and transport of organic matter to the seafloor, hence the relationships found between benthic foraminiferal assemblages and depositional environments (Figs. 4c and 5). Only one station in Olna Firth (MD15-11) has poorly oxygenated bottom waters (O2 < 5 mg l -1 ), which could represent the onset of hypoxic conditions (Vaquer-Sunyer and Duarte, 2008) potentially affecting the distribution of benthic foraminifera (Fig. 2). The benthic 30 foraminiferal assemblage of this location reflects not only organic rich sediments, but also stratified waters and low oxygen concentrations, pointing to a stressed environment (Fig. 4c). However, this is far from the dysoxic/anoxic conditions required for a switch from food-limited to oxygen-limited conditions (Fig. 2). Therefore, the quantity and quality (labile vs. refractory) of organic matter reaching the seafloor is still the major control on the benthic foraminiferal assemblage found at station  4.4 Comparison with other fjords.

Spatial distribution and sources of OM and C.
The processes governing OM and C inputs in west Shetland voes are comparable with those reported for other Scottish fjords by  who found a relationship between the OC content in Scottish fjords and their physical characteristics such as tidal range, precipitation, catchment area and runoff. However, as catchments in west Shetland transport sediment from 10 the adjacent C rich peatlands into the marine system, west Shetland voes overall store more TC (Supplementary Fig. 1) compared to other Scottish fjords (Smeaton et al., 2016). On average, the TC content of the six west Shetland voes is 8.8%, whereas TC in other Scottish fjords ranges between 0-8% (Smeaton et al., 2017 in Fig. 5). Nevertheless, as stated in , there are similarities between Scotland's mainland fjords and these six voes in west Shetland. For example, fjords like Loch Broom and Little Loch Broom in Scotland have relatively unrestricted geomorphologies which makes them highly 15 efficient in capturing IC, similar to Clift Sound and Vaila Sound in Shetland (Fig. 3). Transient low oxygen concentrations in Scotland's mainland fjords seem to be an important factor in the burial and preservation of OC (Gillibrand et al., 2006).

Benthic foraminiferal assemblages
Our knowledge of benthic foraminiferal assemblages in Scottish fjords is limited to studies based on fauna from Loch Etive Nørgaard-Pedersen et al., 2006) and their response to changes in sea level and terrigenous input entering 25 the catchment. Species in common with west Shetland voes include Cibicides lobatulus, Ammonia batavus and E. scaber.
Assemblages dominated by C. lobatulus and A. batavus in Loch Etive indicate marine conditions and intensified bottom currents possibly associated with deep water renewal events (Nørgaard-Pedersen et al., 2006). This agrees well with the assemblages dominated by Cibicides spp. and Rosalina spp. observed in Clift Sound and Vaila Sound suggestive of highenergy environments (Fig. 5). The dominance of agglutinated species such as E. scaber occurs in Loch Etive when OM 30 deposition increases eventually resulting in low oxygen levels in bottom waters and carbonate dissolution Nørgaard-Pedersen et al., 2006). We observed similar environmental conditions and assemblages in Olna Firth, the middle of Busta Voe and inner basin of Sand Sound (Figs. 4c and 5).
Similarities can also be observed between assemblages from the Scottish coast and shelf and west Shetland voes. Klitgaard-Kristensen and Sejrup (2002) and Murray (2003) found high abundance of C. lobatulus and Rosalina sp. in areas with strong 5 currents, low organics and coarse-grain sediments. We do not see a correlation between grain size and Cibicides or Rosalina in our samples due to the multiple forcings driving particle size distribution in Shetland (bottom current intensity vs. terrigenous supply); however we do observe a clear correlation between high relative abundance of Cibicides and Rosalina and low organics and high energy environments (Figs.4c and 5). Klitgaard-Kristensen and Sejrup (2002) and Murray (2003) also found a relationship between the dominance of B. marginata and low oxygen levels, high organics and fine sediments. Additionally, 10 Murray (2003) also observed S. fusiformis dominated assemblages in association with high availability of LOM and low oxygen levels. Overall, this compares really well with the assemblage we described in Group 2 (Figs. 4c and 5).
A recent study by Murray and Alve (2016) synthesise benthic foraminiferal biogeography in NW European fjords, mostly focusing on Norwegian fauna. Similarities with west Shetland voes include C. lobatulus as an "Atlantic species" pointing to the influence of Atlantic Ocean waters on local estuarine/fjordic circulation. Assemblages dominated by E. scaber have been 15 observed in Oslofjord and Drammensfjord in Norway (Alve, 1995;Alve and Nagy, 1986) and Gullmar Fjord in Sweden (Qvale et al., 1984), generally associated with organic rich sediments and shallow waters. Additionally, Klitgaard-Kristensen and Buhl-Mortensen (1999) found high relative abundance of the opportunistic species S. fusiformis in Norwegian fjords with a high per cent of sand and organic carbon, and assemblages characterised by B. marginata in environments rich in ROM.
In our sample set, we had no means to distinguish in situ tests from potentially advected tests when picking and counting 20 unstained specimens, and the close match between our benthic foraminiferal assemblages and published data suggests that advection of material from other locations, especially in the more energetic environments, is negligible. Overall, benthic foraminiferal assemblages do respond to variations in the amount of OM and OC becoming deposited in fjord sediments (Fig.   4c), and comparisons can be drawn between fjordic systems. We recognise that this is a qualitative approach and we are not able to quantify changes in OM or OC content reaching the seafloor based on foraminiferal assemblages. However, this study 25 provides an initial insight into the use of foraminifera to discriminate between post-depositional OC degradation and actual OC burial and accumulation in fjords, as foraminifera would only preserve the latter information in their assemblage composition. We reckon however that extreme conditions at the seafloor may result in post depositional dissolution of benthic foraminiferal tests, which would obviously compromise the assemblage composition. For this reason, we recommend to carefully consider the level of preservation of each assemblage and, when necessary, make appropriate remarks on etching of 30 tests and other possible signs of tests dissolution (e.g. Khanna et al., 2013). Additionally, assemblages dominated by agglutinated species like in Group 2 might have a lower diagnostic potential in their application down-core compared to their hyaline and porcelaneous counterpart due to a tendency of being poorly preserved in the fossil record. Nevertheless, foraminifera are generally abundant in the fossil record and this approach shows down-core potential to support reconstructions of changes in OC burial and accumulation over time and could be paired with statistical and spatial model studies to tease out post depositional OC degradation.

Conclusions
Our records from west Shetland voes provide new insights into sediments provenance and type, modern benthic foraminiferal 5 spatial distribution across environmental gradients, and their use as indicators for OC enrichment in fjordic systems. Our results support a terrestrial source of both OM and OC in land-locked areas and calm depositional settings in west Shetland voes.
Additionally, in locations characterised by low bottom water oxygenation, LOM is present in higher percentages than in well ventilated sites. Conversely, marine processes govern IC storage, with higher concentrations in marginal to open marine settings. Benthic foraminiferal assemblages reliably reflect changes in OC content reaching the seafloor and show great 10 potential as bio-indicators to trace these changes over time. We were able to identify four groups (75% confidence) based on foraminiferal relative abundance data, which reflect four types of depositional settings: 1) Land-locked with high organics and very weak currents; 2) Stressed environment with high organics, heavily stratified water column and low oxygen concentrations; 3) Transitional environment characterised by moderate concentrations of organics and mild or intermittent currents; 4) Marginal to coastal settings with low organics and a high-energy environment. 15 Therefore, studies aiming to reconstruct past variations in the amount of OC reaching the seafloor and becoming stored therein could benefit from the use of benthic foraminifera assemblages as a qualitative approach to understand OC preservation in the sediment column and, for example, ground truth model predictions.
Author contribution. WENA and ELGC designed the study. ELGC, JC and CS analysed sediment physical properties, while 20 ELGC and WENA worked on benthic foraminiferal assemblages. KD provided hydrographic data. ELGC prepared the manuscript with contributions from all co-authors.
Competing interests. The authors declare that they have no conflict of interest.

25
Acknowledgements. The crew of the MV Moder Dy are gratefully acknowledge for support of the work during the cruise in Shetland. We are thankful to undergraduate students Laura Smith, Leonie Strobl and Bethan Hudd for their help in the microscopy lab assisting with benthic foraminifera picking. This project was supported by BBSRC/NERC (ref. BB/M026620/01). Isotope analysis was supported by NERC Life Science Mass Spectrometry Facility (CEH_L_098_11_2015). KD was funded by the EU INTERREG IV Northern Periphery project: Warning of Algal Toxin 30 Events to support Aquaculture in the NPP coastal zone Region "WATER". We thank two anonymous referees and the Handling Associate Editor Dr Hiroshi Kitazato for valuable and constructive comments that improved the manuscript. fjords, Geophys. Res. Lett., 44(3), 1457-1465, doi:10.1002/2016GL072260, 2017 The genetic diversity, phylogeography and morphology of Elphidiidae (Foraminifera) in the Northeast Atlantic, Mar. Micropaleontol., 129, 1-23, doi:10.1016/j.marmicro.2016.09.001, 2016 -4-4005-2007, 2008. 50 doi:10.1111/ejss.12289, 2015: Application of organic carbon and nitrogen stable isotope and C/N ratios as source indicators of organic matter provenance in estuarine systems: Evidence from the tay estuary, Scotland, Estuar. Coast. Shelf Sci., doi:10.1006/ecss.1994.1015, 1994