Experimental setup
The experiment was performed in a stream-side mobile mesocosm facility (MOBICOS) (Fink et al. 2020) located at the forested, upstream part of the Holtemme River (51° 49′ 00.7″ N, 10° 43′ 29.26″ E, Germany). This river section is characterized by a closed canopy and a pristine hydromorphology, moderate NO3 concentrations, and low soluble reactive phosphorus (SRP) concentrations compared to reaches further downstream, which are influenced by urban centers and agriculture (Weitere et al. 2021). Surface water was continuously supplied to the flumes by a submersible pump (type XJ 40, ABS Sulzer, Germany). The inflowing water was filtered with a 50 µm self-cleaning filter (GEFA Process Technik GmbH, Model AP3 DA) to prevent the inflow of suspended particles and macroinvertebrates. The mesocosm setup comprised 24 flumes (72 cm length, 14 cm height, 8 cm width) (Online resource 1, Fig. 1) filled with commercially available sediment previously washed and dried at 120°C. The sediment height was 8 cm allowing for a benthic (0–2 cm depth) and a hyporheic compartment (2–8 cm depth). The height of the water column was 4 cm. Four ceramic tiles (6 x 1 x 0.5 cm) were equidistantly placed in each flume perpendicularly with the flow to facilitate water infiltration into the sediment. Discharge was set to a constant flow rate of 0.17 ± 0.006 L \({s}^{-1}\) in each flume and continuously monitored with a magnetic inductive flow meter (IFM, model: SM8000). Water passed through a twin-wall polycarbonate honeycomb with an opening diameter of 5 mm to ensure laminar flow. The light was provided using one elongated LED lamp (Solar Stinger daylight 16W) placed above each flume, and light: dark cycles were set to 10:14 h.
Treatments
We assigned the 24 flumes to four treatments: control, fine sediment, eutrophication, and fine sediment + eutrophication (Table 1). In the control treatment, we mimicked the environmental conditions found in the Holtemme River by using similar-sized gravel (2–4 mm) and comparable light intensities (i.e., 7.5 ± 1.5 µmol photons \({m}^{-2}{s}^{-1}\)). The fine sediment treatment was created by using fine sand (grain size 0.1–0.4 mm). The eutrophication treatment was created by increasing the light intensity from 7.5 ± 1.5 µmol photons \({m}^{-2}{s}^{-1}\) to 51.2 ± 6.5 µmol photons \({m}^{-2}{s}^{-1}\) and increasing soluble reactive phosporous (SRP) concentrations from below the quantification limit (< 0.003 \(mg {l}^{-1}\)) to 0.04 mg \({l}^{-1}\) by constantly supplying a KH2PO4 solution with a peristaltic pump (Watson-Marlow 205U, Falmouth, UK). To test the effect of stressors interaction, we simultaneously manipulated substrate, SRP (from now on “P”), and light in the eutrophication + fine sediment treatment (Table 1). Each treatment was replicated with six flumes. Flumes were operated under these conditions in a constant flowing mode for two weeks.
In the hyporheic compartment, not only P but also DOC limitation can strongly influence heterotrophic microbial growth because local, autotrophic organic carbon (OC) production can be negligible (Stutter et al. 2020). We expected this to be of little relevance to our experiment due to the high bulk DOC concentrations in the stream. To make sure that DOC limitation was not the case, we tested in an auxiliary experiment with laboratory microcosms to which extent N-NO3 uptake was potentially limited by DOC availability. We used a protocol similar to Graeber et al. 2021. Briefly, we sampled water from each flume outlet once after two weeks of flume colonization and measured the uptake of DOC, dissolved N, and dissolved P fractions in a bioavailability experiment in the dark for an incubation time of eight days. Based on this, we calculated the reactive DOC, reactive dissolved N, and reactive dissolved P concentrations and calculated the C : N : P ratio of reactive, hence microbially available macronutrients. To show the potential reactive DOC or P limitation of heterotrophic NO3 uptake, we used a ternary plot approach using the ggtern package in R (Graeber et al. 2021), in which we normalized the C : N : P to the median of a large analysis of bacterial C : N : P ratios (Godwin and Cotner 2018). We also assessed potential changes in dissolved organic matter (DOM) molecular composition via fluorescence spectroscopy and parallel factor analysis using the staRdom toolbox in R (Pucher et al. 2019). A detailed protocol and results can be found in Online resource 2.
Table 1
Summary of conditions found in the flumes assigned to the different treatments. Values are mean (n = 5) and standard deviation except for sediment grain size (min-max)
| |
Control
|
Eutrophication
|
Fine sediment
|
Eutrophication+
fine sediment
|
|
Sediment grain size (mm)
|
2–4
(Gravel)
|
2–4
(Gravel)
|
0.1–0.4
(Fine sand)
|
0.1–0.4
(Fine sand)
|
|
Light intensity
(µmol photons \({m}^{-2}{s}^{-1}\))
|
7.5 ± 1.5
|
51.2 ± 6.5
|
7.5 ± 1.5
|
51.2 ± 6.5
|
|
DOC (\({mgl}^{-1}\))
|
9.4 ± 0.23
|
9.38 ± 0.21
|
9.41 ± 0.15
|
9.39 ± 0.23
|
|
SRP (\({mgl}^{-1}\))
|
≤ 0.003 ± 0.00
|
0.04 ± 0.00
|
≤ 0.003 ± 0.00
|
0.04 ± 0.00
|
|
N-N-NO3 (\({mgl}^{-1}\))
|
1.05 ± 0.01
|
1.04 ± 0.00
|
1.05 ± 0.01
|
1.05 ± 0.01
|
Tracer addition, sampling, and analyses
After two weeks of colonization, one flume per treatment was sampled to determine ambient isotope values (15N atom fraction, x(15N)) of benthic and hyporheic organic matter (OM). Subsequently, we added an enriched K15N-NO3 solution constantly for 24 h into the turbulent part of the flumes’ inlet section with syringe pumps (NE-1200 Twelve Channel Programmable Syringe Pump, New Era Pump Systems, Inc.). The addition was adjusted to achieve a target x(15N) of 30. After the addition, we collected sediment samples from the benthic and hyporheic compartments of each flume. To prevent sediment mixing between the benthic and hyporheic compartments, we first shoveled the upper 2 cm of sediment and stored it in a sealed container (benthic sample). We then collected a sediment sample from the hyporheic compartment and stored it in another sealed container (hyporheic sample). To separate OM from the mineral substrate, we rinsed the samples multiple times with water and passed them through a 100 µm sieve. The benthic samples exhibited visibly higher OM content compared to the hyporheic samples. Therefore, a larger volume of sediment from the hyporheic samples had to be processed to obtain sufficient OM for further analysis. The processed sediment was transported to the laboratory, where it was dried at 60°C for a minimum of 24 h and then weighed. To determine the volume of the sediment, the weight was divided by the specific density of the sediment, which was measured in the lab using a pycnometer (Blaubrand™).
In the laboratory, OM was freeze-dried (Delta 2–24 LSCplus) for at least 24 h, homogenized, and three replicates of each sample was weighed into tin caps (HEKAtech, Germany) for isotope analysis. Stable nitrogen isotopic composition was analyzed on an elemental analyzer (Flash 2000 Organic Elemental Analyzer, Thermo Fisher Scientific, Germany), directly connected via an open split system (ConFlo IV, Thermo Fisher Scientific, Germany) to an isotope ratio mass spectrometer (Delta V Advantage IRMS, Thermo Fisher Scientific, Germany).
During the analysis, only a minor enrichment was found (i.e., maximum x(15N) = 6.38). Thus, the normalization was done by analyzing the reference materials IAEA-311 (x(15N) = 2.05) (Parr and Clements 1991) and AS-3 (x(15N) = 0.36) (Kornexl et al. 1999) and applying a two-point calibration approach. The analytical precision was below ± 0.2.
Particulate organic carbon (POC) and particulate nitrogen (PN) concentrations of the OM were measured with an elemental analyzer (Elementar Analysen Systeme GmbH, Hanau).
Calculation of N uptake
N-NO3 uptake by benthic and hyporheic OM was calculated from sediment samples collected before and after the isotope tracer addition following Mulholland (2004b). Briefly, we corrected all 15N values by the background values and calculated the excess atom fraction (AFE). Then, we calculated the areal uptake rate \({U}_{Areal}\) and biomass uptake \({U}_{Biomass}\). For that, we first calculated the amount of tracer in the OM (\(15{N}_{OM}\), g 15N):
$$15{N}_{OM}=OM\times \left(\frac{\%N}{100}\right)\times {AFE}_{OM}$$
where OM (g dry weight) is the organic matter stock within the samples, %N is the nitrogen elemental concentrations (%) of the OM, and \({AFE}_{OM}\) the excess 15N fraction of OM.
Then, we scaled the result by sediment volume (\({15N}_{Sed}\), g 15N \({cm}^{-3}\)):
$${15N}_{Sed} = \frac{{15N}_{OM}}{V}$$
Where V is the volume of sediment from which the OM was collected (\({ cm}^{-3}\)). The volume of sediment was calculated by dividing the dry weight (W \(g\)) by the sediment density (d \(g{ cm}^{-3}\)) as:
Where d, i.e. sand density = 2.76 \(g {cm}^{-3}\), gravel density = 2.54 \(g {cm}^{-3}\).
Subsequently, we calculated the N-NO3 uptake rate of the sediment (\({U}_{sed}\) g 15 N \({cm}^{-3}\) \({d}^{-1}\)):
$${{U}_{sed}}_{}=\frac{{15N}_{Sed}}{\left(\frac{{15N}_{Flux}}{{N}_{flux}}\right)}$$
where
$$15NFlux=\left({AFE}_{w,a}\text{*}Q\text{*}C)-({AFE}_{w,b}\text{*}Q\text{*}C\right)$$
$${N}_{flux}=Q\text{*}C$$
\(15NFlux\) is the tracer mass flux, where we subtracted from the total mass flux \(\left({AFE}_{w,a}\right)\) the background mass flux \(\left({AFE}_{w,b}\right)\), C is the N-NO3 concentration (mg \({l}^{-1}\)) and Q is flume discharge (\({l s}^{-1}\)).
We expressed uptake on an areal basis (\({U}_{Areal}\), g N-NO3 \({m}^{-2}{d}^{-1}\)), i.e. N-NO3 uptake of a 1 \({m}^{2}\) surface with a height of 1 cm by multiplying the \({U}_{sed}\) by 10,000.
We calculated biomass uptake (\({U}_{Biomass}\) g N-NO3 \({{g}^{-1}}_{OM}\) \({d}^{-1}\)), i.e. N-NO3 uptake of 1 gram of OM, by dividing the areal uptake by the amount of OM found in that volume as:
$${U}_{Biomass}= \frac{{U}_{Areal}}{OM}$$
We finally summed the areal uptake of the benthic and hyporheic compartments to obtain the total areal uptake (\({U}_{tot}\) g N-NO3 \({m}^{-2}{d}^{-1}\)).
Chlorophyll-a and microbial density
Chlorophyll-a concentrations in the sediment were measured from benthic samples after the two weeks of colonization. Sediment samples were immediately frozen at -20°C and chlorophyll-a was analyzed via high-performance liquid chromatography with a Thermo Scientific UltiMate 3000 HPLC System (Dionex, Thermo Fisher Scientific Corporation, Waltham, MA, USA). Bacterial abundance was measured by collecting a sediment sample from the benthic and hyporheic compartments of each flume in sterile Falcon tubes. To detach cells from the sediment, we used established protocols as described by Chen et al. (2021). Then cells were stained with DAPI and counted using an epifluorescence microscope at 100x magnification. A microscopic investigation was performed using a Zeiss AxioImager.Z2 epifluorescence microscope equipped with an HXP R 120W/45C UV Hg-vapor lamp, Colibri.2 LED illuminations and the following fluorescence filters: DAPI (365/10 nm excitation, 420 LP emission, FT 395 Beam Splitter), Alexa488 (472/30 excitation, 520/35 emission, 495 Beam Splitter), and Alexa594 (562/40 excitation, 624/40 emission, 593 Beam Splitter). The filter sets discriminated clearly between DAPI (365 nm) and natural autofluorescence of the diatom cells. Imaging was done with 100X oil objective numerical aperture N:A 1.4. The software for image acquisition allows for overlapping images acquired successively with different filter sets (Zen software from Carl Zeiss). Images are reported in Online resource 3.
Physical-chemical parameters monitoring
Light intensity was measured with a frequency of five minutes using a light intensity data logger MX2202 (Onset, Bourne, Massachusetts, USA) placed on the sediment of three flumes per treatment. Planar oxygen sensors (SP-PSt6-YAU, PreSens Precision Sensing GmbH, Germany) were glued to the inner wall of three flumes for each treatment to monitor oxygen concentration. Oxygen was measured every 25 seconds at two different depths: 2 cm (benthic) and 6 cm (hyporheic). Some of the cables connecting the planar sensors to the oxygen readers detached from the wall during the measurements, resulting in incomplete temporal data for certain treatments (specifically, 3 out of 9 incomplete datasets for the benthic and 2 out of 9 incomplete datasets for the hyporheic compartment). To ensure comparability of data among treatments, we selected the oxygen data from two flumes per treatment for each compartment where the oxygen data was complete and uninterrupted for a minimum of seven consecutive days (Online resource 1, Fig. 2 and Fig. 3). Water samples for N-NO3, nitrite (N-NO2), ammonium (N-NH4), reactive phosphorous (SRP) were collected and analyzed as described in Online resource 1.
Statistical analysis
To assess treatment effects on N-NO3 areal uptake, biomass uptake, OM, POC, PON, bacterial abundance, and chlorophyll-a, we applied a type I analysis of variance (ANOVA package car, function aov) and a Tukey HSD post-hoc test separately for the benthic and hyporheic compartment. Data were ln(x) transformed where necessary to approximate assumptions of normality and homogeneity of variance of residuals. Values exceeding three times the standard deviation were classified as outliers and removed from the analysis. This happened for four variables: benthic bacterial abundance (2 out of 20), hyporheic POC content (2 out of 20), hyporheic PON (2 out of 20), and benthic OM (2 out of 20). All tests were performed in R (R Core Team 2021).