Carbon bioavailability in a high arctic fjord influenced by glacial meltwater, NE Greenland

The land-to-ocean ﬂux of organic carbon is increasing in glacierized regions in response to increasing temperatures in the Arctic (Hood et al., 2015). In order to understand the response of the coastal ecosystem metabolism to the organic carbon input it is essential to determine the bioavailability of the different carbon sources in the system. We quantiﬁed the bacterial turnover of organic carbon in a high Arctic fjord system (Young Sound, NE Greenland) during the ice-free period (July-October 2014) and assessed the quality and quantity of the 3 major organic carbon sources; (1) local phytoplankton production (2) runoff from land-terminating glaciers and a lowland river and (3) inﬂow from the ocean shelf. We found that despite relatively low concentrations of DOC in the rivers, the bioavailability of the river–DOC was signiﬁcantly higher than in the fjord, and characterized by high cell-speciﬁc bacterial production and low C:N ratios. In contrast, the DOC source entering via inﬂow of coastal shelf waters had high DOC concentrations with high C:N and low speciﬁc bacterial production. The phytoplankton production in the fjord could not sustain the bacterial carbon demand, but was still the major source of organic carbon for bacterial growth. We assessed the bacterial community composition and found that communities were speciﬁc for the different water types i.e., the bacterial community of the coastal inﬂow water could be traced mainly in the subsurface water, while the glacial river community strongly dominated the surface water in the fjord.


INTRODUCTION
Carbon consumption and mineralization by pelagic heterotrophic bacteria play a key role in marine ecosystems. Increasing temperatures, with pronounced effects at high latitudes, have raised questions about how reduced ice-cover and increased runoff can affect the ecosystem metabolism i.e., the balance between respiration and primary production. Bacterial carbon turnover was traditionally suggested to be limited by low temperature in high latitude systems, and consequently play only a minor role in turning over primary production (Pomeroy and Deibel, 1986;Pomeroy et al., 1991). Several high latitude studies during the past decade have however reported bacterial production rates similar to those reported in lowlatitude non-oligotrophic systems (Børsheim, 2000;Rysgaard and Nielsen, 2006;Sejr et al., 2007). Two seasonal studies found annual ranges in bacterial production of 5-42 mg C m −2 d −1 in Kobbefjord (64 • N) (Middelboe et al., 2012) and 90-165 mg C m −2 d −1 in Kongsfjorden, Svalbard (78 • N) (Iversen and Seuthe, 2011), with maximum values during spring, coinciding with the spring phytoplankton bloom.
Estimates of bacterial carbon cycling is often based on measurements of net bacterial production, however variability in the factors used to convert radioisotope (e.g., 14 C-leucine or 3 H-thymidine) incorporation to carbon production affects the growth estimates and potentially complicates comparison of studies. Measurements of the bacterial respiration (BR) are required in order to estimate the bacterial carbon demand (BCD) and growth efficiency (BGE). Such measurements in Arctic systems are few and the BGE reported are highly variable, but all in the low end of those reported in other aquatic systems (del Giorgio and Cole, 1998). It has been hypothesized that low temperatures limit substrate uptake and consequently argued that Arctic bacteria need relatively higher concentrations of carbon to grow at low temperatures (Pomeroy and Wiebe, 2001). Other studies have demonstrated an inverse relationship between BGE and temperature (Rivkin and Legendre, 2001;Apple et al., 2006) leading to speculations that the low efficiency found in the Arctic may instead be a result of poor quality carbon sources (Middelboe et al., 2012). While the concentration of organic matter alone does not reflect the carbon quality (Kirchman et al., 2005), the elemental ratios of the dissolved organic matter (DOM) i.e., the C:N ratio has provided insight on the DOM bioavailability (del Giorgio and Cole, 1998;Pradeep Ram et al., 2003;Kragh and Søndergaard, 2004).
The activity of bacteria is tightly coupled to DOM bioavailability (Amon and Benner, 1996;Kragh and Søndergaard, 2004). Bioavailable dissolved organic carbon (BDOC) has been estimated to constitute on average <1% of the oceanic DOC pool, however elevated in the surface waters (Hansell, 2013). Studies in the Greenland Sea, Fram Strait, and Kobbefjord have shown that BDOC constitutes 13-36% of total DOC in surface water (Middelboe and Lundsgaard, 2003;Middelboe et al., 2012;Jørgensen et al., 2014). Glacial meltwater from both Alaskan (Hood et al., 2009) and Alpine glaciers (Singer et al., 2012) contained highly bioavailable (>60%) DOC. As the Greenland Ice Sheet is melting at record speed (Nghiem et al., 2012) and the melt is projected to continue increasing (Keegan et al., 2014) it poses the question whether coastal bacterial carbon turnover will increase and drive the fjord systems toward more heterotrophic conditions in the future.
The bioavailability of different DOM types is influenced by a number of factors including composition of substrate, availability of mineral nutrients and the bacterial communities and their enzymatic capabilities (Middelboe and Lundsgaard, 2003;Kritzberg et al., 2010;Traving et al., 2016). Only few studies have tried to directly connect specific bacterial groups to different types of DOM (Kirchman et al., 2007;Baña et al., 2013;Osterholz et al., 2016). Meltwater from glaciers has been found to significantly modify the structure of microbial communities in the connected fjord (Gutiérrez et al., 2015). Consequently, increased runoff associated with warming climate will not only affect the transport of organic matter, it may also change the dynamics of coastal bacterial communities toward a higher influence of riverine bacteria and thus potentially changes in BGE and DOM degradation of the coastal bacterial community (Fortunato et al., 2013). Exploring links between the bacterial community composition and the various DOM sources, are therefore highly relevant in Arctic environments. Young Sound receives most of its runoff from the Greenland Ice Sheet via land terminating glaciers (Citterio et al., 2017), resulting in a clear gradient of allochthonous sources of both organic matter and silt throughout the fjord . The organic carbon sources in the fjord comprise two allochthonous carbon sources 1) meltwater from land-terminating glaciers and a lowland river and 2) coastal water that contains traces of DOM from the Arctic Ocean (Amon and Budéus, 2003).
Based on previous findings we hypothesize that the glacial runoff in Young Sound contains highly bioavailable DOM compared to the local production and inflowing coastal water. In order to understand the response in ecosystem metabolism we evaluate the importance of each of the 3 carbon sources as substrate for bacterial carbon degradation and examine associations between the bacterial DOM degradation and the genetic diversity.

Study Site and Sampling
The study was conducted in the high arctic fjord Young Sound, NE  • W). A sill of 45 m depth separates the deeper parts of the fjord from the coastal shelf waters, which are influenced by the East Greenland Current (for more info see Rysgaard et al., 2003). Sampling was conducted at four stations located along a length section from the inner fjord (St. 1) to the shelf waters on the outer side of the sill (St. 4) (Figure 1). The stations 1, 2, 3, and 4 are located according to stations monitored yearly by the Greenland Ecosystem Monitoring (GEM) MarinBasis Zackenberg programme in which they are named Tyro 05, YS 3.18, Standard St. and GH 05, respectively. The fjord stations were each sampled approximately every 10th during the early ice-free period (July 15-August 7) and the late summer period before new ice formation (September 4-October 4) ( Figure S1).
The first sampling was conducted prior to the sea-ice breakup through a hole in the ice at St. 3, when only the central part of the fjord was still ice covered (see satellite photos Figure 1). The ice broke up in the central part on July 15 and the fjord rendered ice-free within 24 h. The remaining sampling was carried out from the research vessel Aage V. Jensen using mini rosette with 12 × 1.7 L Niskin bottles from 6 standard depths (1,10,20,30,40, and 100 m) and 1-2 additional depths at the deep chlorophyll maximum (DCM) when this did not overlap with one of the standard depths. The DCM was determined FIGURE 1 | Study area in Northeast Greenland (74 • N, 21 • W) showing the 4 fjord stations (black dots) and the 3 rivers (blue squares). Above are Satellite TERR images from June-August 2014 illustrating the ice and snow conditions as well as the gradient of riverine input (available at http://ocean.dmi.dk/arctic/daneborg.uk.php).
prior to every sampling using a Satlantic Free-falling Optical Profiler . A Seabird SBE 19+ CTD profiler was deployed at every sampling occasion and recorded vertical profiles of temperature ( • C), salinity (ratio; no units), chlorophyll fluorescence (flu chl , relative; no units), turbidity (FTU), and photosynthetically active radiation (PAR, µmol m −2 s −1 ) at every sampling occasion. The light attenuation was estimated from the CTD-profiles using a two-phase Weibull function as described in Murray (2015).
Glaciers cover ca. 33% of the drainage area of the fjord and land-terminating glaciers contribute 50-80% of the annual terrestrial runoff with highest contribution in the inner fjord (Bendtsen et al., 2014;Citterio et al., 2017). Three of the major rivers discharging into the fjord were sampled as long as sufficient water was flowing (last sampling was on September 10). The meltwater in the Tyroler river (R1) and Clay Bay river (R2) flows from glaciers trough rocky sediment basins with close to zero vegetation for a distance of ca. 0.5 and 2 km, respectively, before they reach the fjord. A model study estimate the residence time of river water in the fjord to be about 2 weeks in July and up to a month in August (Bendtsen et al., 2014). R1 receives water directly from the Greenland Ice Sheet and has the largest catchment area (Bendtsen et al., 2014), while R2 receives meltwater from smaller local glaciers (Figure 1). The largest river, R3 (Zackenberg river) has the second largest catchment area and is connected to 2 lakes. It flows through lowland permafrost soils covered with vegetation types like dwarf shrub heath (Salix arctica) and grasses (e.g., Arctagrostis latifolia), however the riverbed is rocky and without vegetation (Elberling et al., 2008). River water was collected just below surface in 5 L plastic bottles. A 10-year time series of temperature and the organic and inorganic particulate biomass and the dissolved organic carbon (DOC) recorded from the Zackenberg river by the GeoBasis programme by Greenland Ecosystem Monitoring are included in Figure S3.
A total of 25 profiles were sampled at the fjord stations and the rivers were sampled each 3 times. Based on the salinity and temperature (Figure 2) five water types were defined ( Table 1). Bacterial abundance, production and chemical parameters (nutrients, DOC, DON, and chl a) were measured in total 174 times each. A total of 42 samples were collected in the rivers and at 1 m and DCM for "extra" analysis of carbon bioavailability, particle associated bacterial production and community composition analysis. These "focus samples" are marked with large symbols in Figure 2. Environmental data associated with these 42 focus samples are given in Table S1. Note no focus samples were collected from the water mass defined as Shelf water.

Chlorophyll a and Primary Production
Concentrations of chlorophyll a (chl a) were determined according to Jespersen and Christoffersen (1987). Triplicates of 250 mL water was filtered onto GF/F, 2 and 10 µm polycarbonate filters and chl a was extracted in 5 mL 96% ethanol for 12-24 h and analyzed on a Turner Design Fluorometer calibrated against a chl a standard. The measurements were done in FIGURE 2 | Potential temperature (θ) and salinity plot for the open-water season in Young Sound 2014 with water types divided into five categories following the criteria presented in Table 1 and selected based on salinity and temperature profiles ( Figure S1): River water, River plume, Surface water, Subsurface fjord water, and Shelf water. The season is split into two periods: July-August (orange) and September-October (blue). All data points from CTD casts throughout the fjord is shown as small triangles, while larger symbols (circles and triangles) represent the "focus depths" which include measurements of bioavailability, particle associated production and bacterial community composition. The river samples (triangles) are labeled according to origin; Tyroler river (R1), Clay Bay river (R2), and Zackenberg river (R3) and their salinity is set to 0. triplicates. Primary production (PP) was measured as 14 C-uptake (Nielsen, 1952) according to Markager et al. (1999) and measured for 3 size fractions; dissolved (<0.7 µm), small phytoplankton (0.7-10 µm) and lager phytoplankton (>10 µm). Samples were collected at 1 m depth and at one or two additional depth with a notable DCM (26 samples in total). The areal primary production was calculated according to Lyngsgaard et al. (2014). The daily area production was estimated by integrating over 24 h and with depth. The light intensity at each depth was calculated from the light attenuation and the surface light measured at the nearby Zackenberg research station as part of the GEM Programme.

Flow Cytometry
The abundance of bacteria was determined on an Attune R Acoustic Focusing Flow Cytometer (Applied Biosystems by Life technologies) with a syringe-based fluidic system and a 20 mW 488 nm (blue) laser. Samples were fixed with glutaraldehyde (0.5% final conc.) and kept dark at 4 • C until analysis within 12 h. Samples were stained with SYBR Green I (Molecular Probes, Eugene, Oregon, USA) for min 0.5 h at low flow rate of 25 µL min −1 following the protocol of Marie et al. (1999). It should be noted that the bacterial abundance in the rivers was not easily counted by flow cytometry as the inorganic particle signal was high and obscured the counts of free-living bacteria.
To reduce the problem the samples were diluted x10 with TEbuffer.

Nutrients
Unfiltered seawater was filled directly from the Niskin bottles into 30 mL acid washed HDPE bottles and stored at −20 • C. Nitrite and nitrate (NO − 2 +NO − 3 ), phosphate (PO 3− 4 ) and silicic acid (H 4 SiO 4 ) were measured on a Smartchem200 (by AMS Alliance) autoanalyser following procedures as outlined in Wood et al. (1967) for NO − 3 +NO − 2 , Murphy and Riley (1962) for PO 3− 4 and Koroleff (1983) for the determination of H 4 SiO 4 . Concentration of NH + 4 was determined directly in fresh samples using orthophthaladehyde according to Holmes et al. (1999).

Organic Matter Concentration
DOC and DON samples for determining the initial concentration of DOC were collected in 60 mL acid washed HDPE (highdensity polyethylene) bottles and stored frozen (−20 • C) until analysis. DOC is here considered to equal total organic carbon (TOC) according to (Anderson, 2002). DOC concentrations were determined by high temperature combustion (720 • C) using a Shimadzu TOC-V CPH-TN carbon and nitrogen analyser calibrated using a standard series of acetoanilide and the accuracy of the instrument was evaluated using seawater reference material provided by the Hansell CRM (consensus reference material) program. DON was calculated by subtraction of inorganic nitrogen. For particulate organic carbon and nitrogen (POC and PON) a total of 85 samples were collected. 15 L samples from the rivers and 30 L samples from the fjord (collected at 1 m, DCM and 100 m) were filtered onto 47 mm pre-combusted GFF filters using a peristaltic pump. The filters were placed in a desiccator containing concentrated 37% HCl for 12-14 h to remove inorganic carbon. POC and PON was measured using a Carlo Elba NC1500 (Milan, Italy) CHN elemental analyser following the method of Hedges and Stern (1984).

Bacterial Production
Bacterial production was estimated from incorporation of 3 Hthymidine (Riemann et al., 1982). From each water sample, four replicates of 10 mL unfiltered seawater samples were transferred to 20 mL plastic vials. One replicate was immediately amended with 500 µL of 100% trichloroacetic acid (TCA) and served as control. Samples were incubated with 10 nM 3 H-thymidine (final concentration) for 3-5 h at in situ temperature and stopped by addition of 500 µL 100% TCA. Samples were filtered onto 0.2 µm cellulose-nitrate filters, which were subsequently washed 10 times with ice-cold 5% TCA. Filters were transferred to 6-mL plastic vials and stored at −20 • C until analysis. In the laboratory 5 mL of scintillation liquid was added and the radioactivity was counted on a Perkin Elmer Liquid Scintillation Analyzer Tri-Carb 2800TR. The measured thymidine incorporation was converted to cell production assuming 2.0 x10 18 cells produced per mole 3 H thymidine incorporated (Fuhrman and Azam, 1980). Bacterial population growth rate (GR) was calculated as cell production (cells mL −1 d −1 ) divided by the cell abundance (cells mL −1 ). Cell production was converted to bacterial carbon production (BP) assuming 2.0 × 10 −14 g C cell −1 (Lee and Fuhrman, 1987). For the calculation of area-integrated bacterial carbon consumption, the bacterial production was depth-integrated across the upper 100 m, and divided by the BGE measured in July (see below). Additionally, at the "focus depths" and in the rivers the particulate bacterial production was measured as the fraction associated with particles larger than 3 µm and calculated as the total BP subtracted the <3 µm fraction.

Bacterial Respiration and Growth Efficiency (BGE)
Bacterial respiration (BR) was measured as oxygen consumption for ∼48 h at constant temperature in water from 1 m and DCM in triplicate 12 mL gas tight Exetainers equipped with an optical sensor. Water was pre-filtered through a 3 µm-polycarbonate filter to reduce bias from eukaryotic cell respiration and grazing. A sample with 20 µL HgCl served as control. Oxygen was measured every 5 min for >24 h using a 4-channel Fiber-Optic Oxygen Meter (FireSting, Pyroscience) using the program Pyro Oxygen Logger Software version 2.37 (PyroScience). Respiration rates were calculated as the decrease in oxygen concentration (µM) over the incubation time after subtracting control values. Conversion from oxygen consumption to carbon respiration was done assuming a respiratory quotient (RQ) of 0.82 (Søndergaard and Middelboe, 1995). BGE was calculated from measurements of net bacterial production (NBP) and bacterial respiration (BR) in the 3 µm-filtered samples as: BGE (%) = BP/ (BP+BR) × 100. BGE was measured at all four stations at 1 m depth and DCM, and from the three rivers, however only BR measurements that fulfilled the three following criteria were used: (1) constant incubation temperature (± 0.1 • C) for > 24 h, (2) low abundance (<5,000 cells mL −1 ) of small phytoplankton (<3 µm) (as we saw indications of possible respiration contamination from these), and finally (3) the measured respiration rate should not exceed the total respiration rate of the system (measured by T. Dalsgaard unpublished).

Long Term BDOC Experiments
The quantity of DOC available for bacterial degradation (BDOC) was measured in long term (126-148 d) oxygen consumption experiments. A total of 36 incubations were established from the focus depths and rivers (Table S1). 2 L were 0.22 µm filtered (Millipore R Sterivex) into acid washed HDPE bottles (the filters were later used for extraction of nucleic acids). Note we thus did not measure the bioavailability of the particulate organic matter. The bacterial inoculum was prepared by GF/F filtering 100 mL into 2 × 50 mL falcon tubes. The 0.22 µmfiltered water and bacterial inoculum was stored cold (2 • C) until experimental set-up. The oxygen consumption experiment was set up with five replicate 65 mL Winkler glass bottles equipped with an optical oxygen sensor for each sample and incubated in dark at 8 • C for 148 days (samples collected in July and August) or 126 days (samples collected September and October). Prior to incubation NO − 3 and PO 3− 4 (final conc. 5 µM and 1 µM, respectively) were added along with the bacterial inoculum (10% vol.) to ensure that N or P was not limiting C-degradation during incubation. Consequently, these measurements do not reflect in situ conditions neither the bacterial communities at the time of sampling, but are rather quantitative measure of the bioavailable DOC pool. All Winkler bottles contained a magnet to ensure mixing, and were incubated in a water bath in order to minimize oxygen contamination. In addition, parallel bottles containing 100 and 0% air saturated seawater were measured to correct oxygen measurements for deviations in the 100% control. As control incubations, triplicates of sample water were incubated without bacterial inoculum as well as a sample with 20 µL HgCl. The change in % air saturation over time was monitored every 7-14 days using a Fibox 3 fiber optic patch oxygen sensor (Presense) calibrated with 100% and 0% air saturation using the program OxyView-PST3-V6.02. As for the BGE measurements the change in oxygen concentration over time was converted to carbon consumption assuming an RQ of 0.82 (Søndergaard and Middelboe, 1995). DOC concentration was measured initially and by the end of the incubation.

Nucleic Acids Extraction, Amplification, and Amplicon Sequencing
All environmental samples for molecular analysis were collected by filtering water onto 0.22 µm pore size Millipore R Sterivex filters (the filtrate was used for bioavailability measurements described above). Note that samples were not prefiltered i.e., also the particle associated bacteria are included in this analysis. The filters were immediately frozen and stored at -80 • C until nucleic acid extraction. DNA and RNA were extracted simultaneously using the AllPrep DNA/RNA Mini Kit (Qiagen, Hilden, Germany) according to manufacturer's instructions with modifications for extraction from Sterivex filters as in Paulsen et al. (2016). RNA was subsequently treated with the DNAfree DNA Removal kit (Invitrogen, CA, USA) and reverse transcribed using the SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen). Amplification of cDNA and DNA was performed using a two-step nested PCR approach with primers 519F (CAGCMGCCGCGGTAA; Øvreås et al. (1997) and 806R (GGACTACHVGGGTWTCTAAT; Caporaso et al. (2011) targeting the bacterial (and Archaeal) 16S rRNA gene V4 hypervariable region. For the first PCR step triplicate samples were amplified in reaction volumes of 20 µL including 10 ng DNA or cDNA, 10 µL HotStarTaq Master Mix (Qiagen), 500 nM of each primer and nuclease free water. PCR cycles consisted of an initial denaturation of 15 min at 95 • C, followed by 25 cycles of 95 • C for 20 s, 55 • C for 30 s and 72 • C for 30 s and a final extension step of 72 • C for 7 min. Triplicate PCR products were pooled, purified using the DNA Clean & Concentrator-5 kit (Zymo Research Corporation, CA, USA) and quantified using the Qubit 3.0 Fluorometer. For the second PCR step, 10 ng of pooled PCR product was used in a reaction mixture containing 25 µL HotStarTaq Master Mix, 500 nM of each nested primer with an unique eightnucleotide barcode (total of 96 combinations) and nucleasefree water to bring the mixture to the total volume of 50 µL. Thermal cycles had an initial denaturation for 15 min at 95 • C, followed by 15 cycles at 95 • C for 20 s, 62 • C for 30 s, 72 • C for 30 s, and a final extension step of 72 • C for 7 min. PCR products were purified using Agencourt AMPure XP Beads (Beckman Coulter Inc., CA, USA) and prepared for sequencing by pooling the amplicons in equimolar amounts.
The quality and concentration of the amplicon pool were assessed by agarose gel electrophoresis and by using a Qubit 3.0 Fluorometer, before sending to the Norwegian Sequencing Centre (Oslo, Norway) for High-Throughput Sequencing on a MiSeq platform (Illumina, CA, USA) using the MiSeq Reagent Kit v2 (Illumina). Sequencing data is available at "The European Bioinformatics Institute" under study accession number PRJEB16067 (http://www.ebi.ac.uk).

16S rRNA Gene Sequence Analysis
Paired-end sequences were processed using different bioinformatic tools incorporated on a qiime-processing platform (Caporaso et al., 2010). FASTQ files were quality end-trimmed at a phred quality score ≥24 using Trimmomatic (Bolger et al., 2014) and merged using PANDAseq (Masella et al., 2012), while all reads <200 bp were removed. Prokaryotic OTUs were selected at a sequence similarity threshold of 97% using a de novo uclust (Edgar, 2010) OTU clustering method with default parameters and taxonomy assigned using the Silva 111 reference database (Quast et al., 2013). OTUs with a taxonomic identification were assembled to an OTU table providing abundances for each sample excluding singletons and rarefied to the number of sequences of the smallest samples (5,000 sequences). A total of 4,096,371 sequences were retrieved from the Illumina sequencing of the 16S rRNA gene V4 hypervariable region from total RNA across 52 samples. After removal of singletons, unassigned OTUs and chloroplast reads, sequences were rarefied to 5,000 reads per sample, with a total of 15,922 unique OTUs at 97% sequence identity. Multivariate statistical analysis was performed on basis of the rarefied OTU matrix to explain variations in the data and test for multivariate environmental correlation with the prokaryotic community structure. Bray-Curtis resemblance, ANOSIM, principal component analysis and redundancy analysis were calculated using primer-e version 6 (Plymouth, UK) and Canoco 5 (Ter Braak and Šmilauer, 2012).

Source Tracker Analysis
To illustrate the spread of bacterial communities in the fjord the SourceTracker 0.9.5 software (Knights et al., 2011) was applied in QIIME. It is designed to track the relative contribution of predefined microbial sources in sink samples using a Bayesian approach, as done in Storesund et al., in review. The relative abundance of 16S rRNA genes of all focus samples was used as input data. The OTU table, comprising all OTUs with a taxonomic identification (excluding singletons and chloroplast reads) was rarefied to 1,000 sequences and filtered to only include OTUs that were abundant in more than 3 samples. The rivers (R1, R2 and R3) and as a proxy for inflowing coastal water St. 4 DCM samples were used as "source populations." All remaining stations during the sampling period from July until October were defined as the "sink samples." The result is given as % likely origin from the 4 defined sources. The remaining are categorized as "unknown source". The result were visualized as fjord transects with weightedaverage extrapolation between points using Ocean Data View (Schlitzer, 2016).

Hydrography
During the first period (July-August) the fjord was stratified by a strong halocline at 5-6 m depth with low salinity (<20) in the surface and more saline bottom water (>30). The stratification was strongest in the inner fjord and the salinity of the surface water increased eastwards toward the shelf ( Figure S1). The later period (September-October) was in general colder (surface water <2 • C) and frequent storms mixed the upper layer, which deepened to 30 m at the two inner stations and down to 50 and 80 m at St. 3 and 4, respectively ( Figure S1). Five water types were defined for this study based on the thermohaline properties (Table 1 and Figure 2). This was done to facilitate the interpretation of the results as the dominant carbon sources can be expected to differ across these water types. The "Shelf water" represented the mixing of waters from the East Greenland Current (T∼1.5 • C, S < 34) with warmer and more saline Atlantic water (S > 34.4). "Surface water" is influenced by runoff and the seasonal surface heating, while "Subsurface fjord" water represents the shelf water that has entered the fjord passing over the sill and gradually being mixed with the "Surface water." We defined the "River plume waters" to be the fjord surface waters that were under direct influence of the river discharge. The change between the two periods is apparent in the T-S plot (Figure 2). The runoff was strongest in July and during the last river sampling on September 10 the flow from the glacial rivers (R1 and R2) had almost terminated, while Zackenberg river (R3) was still flowing until end-September. The inner stations were strongly affected by the river silt, with a high turbidity of surface water and the photic zone was therefore initially shallow at St. 1 and 2 (<5 and 10 m, respectively), however later in September it deepened to 25-35 m as turbidity decreased. The opposite trend was observed at the outermost station where the photic depth decreased from 35 to 10 m from July to October, due to the decreasing irradiation.

Bacterial Abundance and Chlorophyll a
While the central fjord was still ice-covered (Figure 1), the inner part was ice-free and a reduction in nutrients measured below the ice at St. 3; NO − 2 +NO − 3 , PO 3− 4 , SiO 4 , in the surface (0.01, 0.33, 1 µM) compared to 100 m values (3.2, 0.7, 6.4 µM) indicated that a phytoplankton spring bloom had already initiated in the inner fjord prior to current sampling program. Chl a was highest (up to 3.1 µg L −1 ) in the outer region of the fjord in July and August with a deep maximum at 20-40 m (Figure 3D), while in the late period chl a was highest in the inner fjord ( Figure 3A). The average bacterial abundance (BA) in the surface water was 3.2 ± 1.0 × 10 5 mL −1 and generally peaked at 0-20 m in the first period, with a maximum abundance in the outer fjord. BA was significantly higher in the upper 20 m for 19 out of 25 profiles, thus the BA maxima were decoupled from the chl a, especially in the first period where the chl a max was deep (20-40 m) (Figure 3). Overall only a weak significant linear correlation was found between BA and chl a (r 2 = 0.03, p = 0.03, n = 164). At the innermost station BA did not correlate with chl a at any time, but there was a positive correlation between BA and turbidity (r 2 = 0.757, p < 0.01, n = 31). At the outermost station chl a correlated with BA during the study (r 2 = 0.44, p < 0.01). At the two mid-fjord stations (2 and 3), correlation between BA and chl a was found only after the runoff had ceased (St. 2: r 2 = 0.83, p < 0.01 and St. 3: r 2 = 0.64, p < 0.01). BA in the rivers was significantly lower than in the fjord surface waters (Table 3). Abundances ranged from a minimum of 1.3 × 10 5 in the glacial rivers (R1 and R2) to 6.4 × 10 5 cells mL −1 in R3 ( Table 2). Note that only the free-living bacteria were enumerated.

Bacterial Production and Particle Association
Bacterial production (BP) was highest in the beginning of the sampling period (17-21 July) with a maximum of 2 µg C L −1 d −1 in the river plume water (St.1 + 2, 1 m) (Table 2, Figures 4A,B). The third highest measure of BP (1.7 µg C L −1 d −1 ) was found in the fresh water layer just below the ice sampled on July 11. BP was lowest within the Shelf water at St. 4 where it remained below 0.09 µg C L −1 d −1 throughout the entire study period, despite BA being relatively high. BP in the rivers was on average 0.39 ± 0.07 µg C L −1 d −1 and 0.32 ± 0.06 µg C L −1 d −1 in the first and the second period, respectively, and thus generally higher than those measured in the fjord despite low BA (Table 3). We found a larger fraction of BP in the R3 to be free-living (Table 2), which is possibly due to the presence of lakes acting as sedimentation basins. The contribution of particle-associated (>3 µm) bacterial production was considerable in both the fjord and rivers. In 17 out of 25 cases particle association was higher in the DCM sample than in the associated 1 m sample, however there was no significant difference between the %particle associations at the two depths (Figure 4).

Bacterial Carbon Demand and Primary Production
Only two of the BR measurements fulfilled all criteria for solid measurements of oxygen consumption throughout the incubation. These gave growth efficiency values of 7.3 ± 1.0 and 6.4 ± 2.0%, during incubation at the in situ temperature of −1.1 • C (St. 3, July) and at 3 • C (St. 4), respectively. The average BGE of 6.9% was applied for determination of bacterial carbon demand (BCD) to allow comparison with the total amount of carbon fixed by planktonic primary production (PP). When integrated over the photic zone the estimated BCD:PP was on average 1.7 ± 1.2 across the sampling period, suggesting that bacterial carbon demand could not be sustained by the local phytoplankton production. There were no clear spatial or temporal trends in the BCD:PP ratio (Figure 5), and estimates were similar when integrated to 100 m. PP was highest initially and could support a high BCD at St. 1 and 3, while bacteria were not sustained by fresh PP initially at St. 2 and 4 in this period. Toward the end of the open water period different patterns developed as bacteria in the inner fjord could be sustained by fresh PP, while BCD was decoupled (up to 5 times higher) from PP at the outermost station (Figure 5). The dissolved fraction (<0.7 µm) of PP was in general high in the fjord and contributed 39-52% in the first period and less, 27-36%, in the late period. The largest fraction contributed most at the outer station (max 32%) and in the last period it was never higher than 13%, whereas the fraction 0.7-10 µm became dominant in the late period ( Figure 5).

Organic Matter Concentrations and C:N Ratios
DOC concentrations showed little systematic vertical variability except for the surface samples at 1 m, which had significantly lower concentrations in 19 out of 25 profiles ( Figure S2). On average, the DOC concentration decreased from the outer part to the inner fjord from 130 ± 16 µM at St. 4 to 106 ± 30 µM at St. 1 (Figure 6). From the first to the second sampling period where the mixed layer deepened, the DOC concentration at 1 m increased, while at the DCM there was a slight decrease. The DOC concentration was highest in the deep water types characterized as Subsurface fjord water (101 ± 20 µM, n = 49) and slightly lower in the Surface water (97 ± 27 µM, n = 97). The River water had significantly lower DOC (40 ± 13 µM, n = 24), and thus the River plume water was diluted to an averaged concentration of Where there is no SD, only one sample was available. 67 ± 13 µM, n = 7. DOC:DON was also lowest in the River water (8 ± 2) and the River plume water (11 ± 2), while the Shelf water had a significantly higher C:N ratio (21 ± 8) ( Table 3).
Maximum POC concentration was found in the rivers with up to 14 µM, while in the fjord a maximum of 4 µM POC was found at the DCM. In the fjord the particulate fraction (%) of the total organic matter was in general minor (avg. 2.3 ± 1, max = 5%), FIGURE 4 | Vertical distribution (upper 100 m) at Station 1 (A), Station 2 (B), Station 3 (C) and Station 4 (D) of bacterial production ± SD, n = 3 (µg C L −1 d −1 ). The horizontal bars indicate the relative distribution between particle-associated bacterial production and free-living production (%) at 1 m and at the deep chlorophyll max (DCM) from July to October at the 4 stations.
while in the rivers the contribution of POC was significant (avg. 13 ± 8, max = 35%) (Tables 2, 3). This gave a significant negative relationship between salinity and POC (r 2 = 0.3, p < 0.001, n = 86). POC did not correlate with chl a, turbidity or to particulate bacterial production. The C:N ratio of the particulate matter was significantly lower in the rivers than in the fjord (p = 0.047, Ftest). In the fjord, the C:N ratio increased from the first period to the end of the open water season (p = 0.0023, F-test), whereas no trend was observed in the rivers.
Inorganic nitrogen (NO − 2 +NO − 3 ) was the limiting inorganic nutrient for primary production and was reduced to 0.4 µM in the surface water during the entire period, while the background level in the Shelf water was ca. 6 µM (Table 3, Figure S2). Ammonium (NH + 4 ) was only measured at St. 3 and ranged between 0.05 and 0.4 µM, with a maximum at 40 m (below DCM).

Bioavailability of DOM
Despite the increase in average DOC concentration from St. 1 to 4, there was a slight decrease in the concentration of BDOC from 19 ± 10 to 11 ± 8 µM from St. 1 and 4, respectively (when values were averaged for 1 m and DCM over the entire period). The fraction of BDOC relative to total DOC (%BDOC) thus decreased from the inner to the outer part of the fjord (Figure 6). BDOC concentrations in the rivers were relatively high (avg. 18 ± 8 µM), resulting in significantly higher %BDOC in the rivers than in the fjord (p < 0.005, t-test) (Figure 6, Table 3). Averaging %BDOC within each water type revealed a decrease as river water was mixed with the fjord water, i.e., River water : River plume : Surface water : Subsurface water ( Table 3).
There was a significant negative correlation between the relative abundance of C to N (DOC:DON) and %BDOC (r 2 = 0.17, p = 0.08, n = 42). C:N ratios were generally higher in the last period and the highest C:N of 29 ± 3 was found in the Shelf water.

Bacterial Community Composition
In general, Proteobacteria were the most abundant bacteria phylum in both fjord and river samples (∼86% of the bacterial phylum). Differences in community composition were observed at class level, with fjord DCM samples containing more Alphaproteobacteria than the rivers and the Subsurface water being strongly dominated by Gammaproteopbacteria (80%), opposed to surface samples and rivers samples with a higher share of Betaproteobacteria (Table 3). Cyanobacteria were present in all rivers and constituted up to 2.7 ± 1.7% of the microbial community in the innermost stations. Alpha diversity (described by the Shannon index) ranged from 4.9 to 9 (Table S2) indicating a greater diversity in the rivers and lowest in fjord samples from July and August (Figure 7). A redundancy analysis (including spatial and temporal parameters) was performed to identify factors that significantly affect bacterial community composition. This explained 47.2% of the total variation in the OTU data, with "water type" (19.5%; river (13%), p = 0.002) and time of sampling expressed as "month" (10.8%; September (8.5%), p = 0.002) being the most significant variables. As "station" generally explained <5% this variable was not included. Samples from all three rivers with overall higher species richness clustered together. The majority of the September fjord samples clustered tightly together, independent of their respective water types and differences between surface and DCM were minimal. Samples from July and August however clustered according to sample depth (i.e., surface or DCM) (Figure 7). Explanatory variables (water type, month and station) explain in total less than 50% of the variation in diversity for the entire data set.
To elucidate changes in community composition throughout the fjord the complexity of the dataset was reduced by constraining a phylogenetic analysis to only include the most abundant bacterial taxa (relative abundance >1%). Further the temporal factor was reduced by restricting the analysis to cover a 10-day period in early August (Figure 8). The heat map shows that all three rivers were very different from the fjord community, and that each river had unique taxa as the most abundant (blue = high, red = low relative abundance). While bacteria found in R1 (with the closest connection to the Greenland Ice Sheet) all could be found in the two other rivers, the R3 (that runs through vegetation-covered catchment area) included some unique families e.g., Granulosicoccaceae, Alcaligenaceae, and an unknown family from the order vadinha64 (the latter presented as "uncultured_bacterium" from the class Opitutae in Figure 8). In early August the majority of the most abundant river taxa were absent in the fjord samples, with the exception of R1 that shared a great number of its most abundant taxa with the nearby Station 1 surface sample which had more unique taxa than any other station (Figures 1, 8). Stations 2, 3, and 4 showed a higher number of shared taxa indicating a gradient from the inner to the outer fjord. The opposite was observed for the deeper DCM samples, which showed a gradual decrease of certain taxa from the outer St. 4 toward the innermost St. 1.
The SourceTracker analysis (including the full dataset at OTU level) showed clear contribution of especially the R1 community to the surface fjord stations with maxima of 75, 63, and 66% similarity at St. 1, 3, and 4, respectively in July-August and ca. 40% in September (Figure 9). The coastal community (source set as St. 4 DCM) entering the outer part of the fjord dominated strongly in the deeper fjord samples throughout the entire study (Figure 9 and Table S1). However, the coastal communities also dominated at the surface waters of St. 2 in the first period, where R1 only contributed 15%. As the R1 community was also present in R2 and R3, overall river contribution is incorporated in the R1 plot (Figures 9A,B). The community unique to R2 and R3 only contributed to minor degree (7-15%) to the fjord surface community. In order to investigate the influence of environmental parameters exerted on the bacterial community structure, we performed canonical correspondence analysis (CCA) between numerical factors, such as bacterial production, chl a fluorescence, temperature, salinity, DOC, and %BDOC, (Figure S2). Due to the complexity of the dataset no strong correlations were found with this analysis. However, when single taxa were correlated to specific environmental parameters we found that especially genera from the class Gammaproteobacteria showed correlations with specific DOM characteristics e.g., a strong positive correlation was found between the genus Glaciecola (order: Alteromonadales) and bacterial production (r 2 = 0.5283; p < 0.0001) with a maximum relative abundance of 25% when BP was highest. Another unknown genus from the order Alteromonadales showed a positive correlation with %BDOC (r 2 = 0.2720; p < 0.0075).

Bioavailability of Allochthonous DOM Sources
River input in marine systems are often a source of high DOM concentrations and therefore DOM often correlate negatively to salinity (Cauwet, 2002). The Young Sound system however deviates from this trend due to two factors. Firstly the freshwater input in fjord has low DOC concentrations due to the dominance of glacial meltwater and limited catchment vegetation. Secondly, we hypothesize that the coastal shelf waters entering the fjord are characterized by high levels of terrestrial organic matter that originates predominantly from Siberian rivers (we did however not sample sufficiently deep at St. 4 to capture the pure Polar water). These rivers discharge into the Arctic Ocean (Amon and Budéus, 2003) and the terrestrial DOC is retained in the surface waters exiting the Arctic Ocean via the Fram Strait as part of the East Greenland Current (Granskog et al., 2012). A large fraction of the terrestrially derived DOM transported by the major rivers Ob and Yenisei is found to be refractory (Meon and Amon, 2004), which explain the low bacterial activity and low DOC bioavailability we found in the Shelf water despite high concentrations. These conditions therefore result in a positive correlation between DOC and salinity (r = 0.51, p < 0.001, df = 175) in Young Sound.
In addition to the quantitative difference, the qualitative C:N ratio of the two allochthonous DOM sources differ. The ratio between bioavailable DOC and inorganic nitrogen exceeded the Redfield ratio (C:N:P = 106:16:1) by thousand fold, whereas the BDOC:PO 3− 4 ratio was 45 ± 30. This emphasizes the importance of DON as the main source of nitrogen for bacteria. The low C:N ratios in the river water resemble those reported from Alaskan glacial rivers, where the relatively high source of DON is explained by microbial production of protein-rich DOM in the subglacial environment (Hood and Scott, 2008). The concentration of bioavailable DOM in the rivers entering Young Sound was slightly lower than values obtained in Alaskan glacier outflow (Hood et al., 2009), but very similar to other studies from the Greenland ice sheet meltwater (Lawson et al., 2014). In contrast to the river DOC, the allochthonous DOC entering the fjord from the open ocean had high C:N ratios, similar to Arctic surface water (Benner et al., 2005). Together, our results demonstrate that open ocean DOC is less labile than the DOC produced in the fjord and that supplied from the rivers.
As expected the load of particulate organic carbon was relatively high in the rivers (Hood et al., 2015). However a high POC-signal was not traceable in the surface water of the inner fjord stations as was the case e.g., for the silt particles (measured as turbidity) and silicate. This indicates that the POC has been lost from water column either by sedimentation, dissolution or bacterial degradation. Bioavailability of POC was not determined, and BDOC values thus potentially represent an underestimation of labile organic carbon concentration. However, since POC on average accounted for 2.3 ± 1.0% and 13 ± 8% of total DOC in the fjord and rivers, respectively, the contribution of bioavailable POC to total bioavailable carbon is probably relatively small. The spatial gradient in %BDOC along the fjord transect suggests a gradual consumption of the labile DOC entering the fjord via the rivers. Consumption of riverine DOC in the fjord is supported by the high BP in the river and river plume water (Table 3) and the negative correlation of BP to salinity. The differences in DOM concentration and composition has been suggested as a driver for diversification of bacterial communities (Crump and Hobbie, 2005;Blanchet et al., 2016;Roiha et al., 2016). Our study FIGURE 8 | Heat map displaying the highly represented (relative abundance>1%) bacterial taxa at family level (the class level is indicated from the left hand color-code) across each of the 3 rivers and the 4 fjord stations (1 m and DCM) all sampled from August 1 to 10. The highest abundance of each of the 45 bacterial families (including cyanobacteria) was used as respective maximum and the remaining samples are given as a percentage relative to the maximum. The darkest blue illustrates maximum, yellow medium, and red lowest abundances.
suggests that the labile character of the river-DOM may have a role stimulating and shaping the activity and structure of bacterial communities in the fjord, by favoring fast growing bacteria. In samples where the community structure was analyzed, high BP was negatively correlated to salinity (r 2 = 0.2743; p < 0.0103). Further, high bacterial production was associated with growth of specific taxa, such as a positive correlation between BP and the relative abundance of the genus Glaciecola.

Bacterial Community Composition
In general, our results suggest that the bacterial community is largely structured by the different water sources to the fjord and changes along the salinity gradient, as also found in other coastal environments (del Giorgio and Bouvier, 2002;Crump and Hobbie, 2005;Gutiérrez et al., 2015). Especially in the early period large differences were observed between 1 m and DCM bacterial communities, which were explained by the  ) is plotted for the two periods using SourceTracker analysis. Transects are plotted using weighted-average extrapolation between points using Ocean Data View (Schlitzer, 2016). strong stratification at 5-7 m in this period ( Figure S1). This resulted in lower species richness during stratification in July and August (Figure 8), possibly because bacterial communities exhibit environmental niche partitioning when the water masses remain separated in this period (Morris et al., 2005;Delong et al., 2006;Chow et al., 2013;Salter et al., 2014). As river runoff decreased and mixing events increased (September) the resemblance between communities at 1 m and DCM naturally increased. One specific strong storm that lasted for several days (21-26 September) may explain the high similarity of 1 m surface samples from St. 1 and 2 to the river samples (Figure 8), as samples were collected immediately after the storm where the disturbance apparently caused transport of terrestrial bacteria to the fjord (Crump and Hobbie, 2005).
In temperate and Polar Regions Flavobacteriia often dominate during phytoplankton bloom (Wilson et al., 2017). We also found this group to dominate when chl a concentrations were high (i.e., DCM samples at St. 4), while they contributed only marginally (<0.1%) to the total bacterial community in the remaining fjord and particularly little in the rivers (Figure 8). Cyanobacteria were relatively abundant in the glacial runoff and can be considered as freshwater tracer, since marine cyanobacteria are usually not found in these Arctic regions (Paulsen et al., 2016). The bacterial communities in the rivers were highly specific to each of them (Figure 8), which we suggest is due to their difference in origin and catchment area; close connection (0.5 km) to the Greenland Ice Sheet (R1), longer distance (2 km) from smaller local glacier (R2), and the lowland vegetation rich river with lake connection (R3). This is in agreement with a recent study from West Greenland that similarly show distinct communities in rivers and a proglacial lake over a 2 km distance (Hauptmann et al., 2016), and find less terrestrial species in river samples with a more direct connection to glaciers.
In the present study R1 has the closest connection to the glacier and all the bacterial families found in R1 were also found in the two other rivers. The SourceTracker analysis (Figure 9) revealed that the R1 community had the far most influence on the surface water, despite not being the largest river. The species unique to R3 hardly contributed to the fjord community, despite being the largest river in terms of volume. This suggests that there is a higher potential for the glacial bacterial community to persist in the fjord than terrestrial communities. Gutiérrez et al. (2015) also found persistence of specific bacterial communities in the surface water of a glacier influenced fjord in Patagonia, and suggested that the glacial meltwater community in the surface was maintained by the competitive advantage of tolerating the cold fresh conditions. Given the averaged doubling time of bacteria in the surface water of Young Sound of 11 ± 7 days (fastest doubling time of 0.9 days) and a transport time of 14-30 days from innermost to the mouth of the fjord, the persistence of the specific communities may be a result of fast transport within the surface water and limited mixing. The differentiation between surface and subsurface communities was lost when the thermal stratification broke down in present study and in Gutiérrez et al. (2015).

Annual Estimations of Carbon Production and Turnover in Young Sound
Due to the strong stratification, high turbidity and input of allochthonous carbon sources, relatively heterotrophic conditions were expected to prevail in the fjord system as concluded in Nielsen et al. (2007). The seasonal and spatial resolution of bacterial carbon demand and primary production allowed a rough, but more robust than previous estimates of the annual carbon budget in Young Sound (Rysgaard and Nielsen, 2006;Nielsen et al., 2007), as previous studies cover a shorter period of the productive season and use a literature BGE of 33%. However, the present study did not cover the early productive season during ice cover in June. As an attempt to account for that we extrapolate the carbon uptake and primary production during the ice-covered productive period using the values obtained on July 11 at St. 3 during ice-cover and by assuming a productive period of 125 days. Based on these assumptions, the annual bacterial production amounted to 1.3 g C m −2 year −1 , corresponding to a carbon demand of 18 g C m −2 year −1 in the upper 100 m. This BP value is ∼3 times lower than a previous estimate of net annual bacterial production in Young Sound of 4.2 g C m −2 year −1 , based on measurements conducted during ice-cover at St. 3 (Rysgaard and Nielsen, 2006), emphasizing that the conditions at St. 3 in June are not representative for the late open water period (August-October) or for the entire Fjord system.
The assessment of bacterial carbon demand rely greatly on estimates of BGE, which is known to depend on a number of factors including DOC composition and lability, bacterial community and temperature, thus it is important to emphasize that BGE is likely to vary across time and space in Young Sound. The BGE we found (6.9%) is low, but still in line with other Arctic studies ranging from 2.2 ± 2.1% in the Chucki Sea (Cota et al., 1996), 6.3 ± 3% in the Fram Strait (Kritzberg et al., 2010), 19.1 ± 9.5% in the rivers Ob, Yenisei and the adjacent Kara Sea (Meon and Amon, 2004) and 9.5 ± 8.7% in Kobbefjord, Greenland (Middelboe et al., 2012).
The bacterial carbon demand on average exceeded primary production during the study period ( Figure 5) and indicated that only about 30% of the PP is dissolved and thus potentially available for bacterial uptake. Further, protist and copepod grazing of phytoplankton consumed a substantial fraction of the particulate PP (Middelbo et al. submitted;Arendt et al., 2016). Even though some of the particulate primary production are eventually degraded by pelagic bacteria as detrital matter, our results strongly indicated that primary production could not sustain bacterial carbon demand in the fjord within the productive 125 days and much less on an annual basis. Consequently, the allochthonous DOM sources contributes significantly to the bacterial carbon turnover in Young Sound.
The current study presents an overview of the organic carbon sources and their turnover in the high Arctic Young Sound. Further it highlights that the meltwater associated bacterial community from the glacial rivers persists and is actively transforming the river-DOM within the freshwater lens. The calculations of annual carbon production and turnover are obviously associated with large uncertainties regarding extrapolations across time and space and use of factors for converting thymidine incorporation to carbon production and net bacterial production to carbon demand. However, by applying a high temporal and spatial resolution of sampling and on-site estimation of BGE, this study provides a relatively solid estimation of the annual carbon budget in a high Arctic Fjord compared to previous studies. The fjord is net heterotrophic and in future scenarios with increasing temperatures the relative contribution of riverine DOC is expected to increase, driving the system toward more heterotrophic conditions. However, more measurements of bacterial growth efficiency and factors controlling this are required to provide more solid budget estimates of bacterial carbon consumption on an ecosystem scale.

AUTHOR CONTRIBUTIONS
MP and SN led the collection and analysis of data and the writing of the paper with equal contribution. MM led the planning of the study, and contributed with sample collection, data analysis and contributing greatly to the writing of the paper. OM led the collection, laboratory work and analysis of the molecular data. SM led the planning and calculations of primary production. AD led the collection and measurements of POC and PON. All other authors contributed to writing the paper. In addition EM, AL, TJ, and MS helped collecting data and CS helped analysing the organic matter samples.

ACKNOWLEDGMENTS
Big thanks to Egon Frandsen, Kunuk Lennert, and Ivali Lennert for organization, expertise in the area and technical assistance in the field and to Johan the polar fox for entertainment. We are grateful to be able to include data from the Greenland Ecosystem Monitoring Programme were provided by the Department of Bioscience, Aarhus University, Denmark in collaboration with Department of Geosciences and Natural Resource Management, Copenhagen University, Denmark.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmars. 2017.00176/full#supplementary-material Figure S1 | Profiles of salinity, temperature ( • C), DOC (µM), DOC:DON and the concentration of nutrients NO − 2 +NO − 3 and Si (µM) at the 4 stations in the periods from July-August (red colors) and September-October (blue colors). Figure S2 | Canonical correspondence analysis (CCA) for sequencing data of 16S rRNA gene fragments at RNA level (most abundant OTUs >1%). Arrows indicate selected environmental variables: bacterial production, dissolved organic carbon (DOC), percent bioavailable dissolved organic carbon (%BDOC), fluorescence, temperature and salinity. The OTUs selected for further analysis are labeled with their genus classification. Figure S3 | Time-series measurements from Zackenberg River of temperature ( • C), Sediment (mg L −1 ), LOI (mg L −1 ) (LOI is organic matter determined by Loss On Ignition when heated at 105 • C ≈ particulate organic matter), DOC (µM). Data is collected by Greenland Ecosystem Monitoring, the GEOBASIS Programme.