- 1Microsensors Group, Max Planck Institute for Marine Microbiology, Bremen, Germany
- 2Biogeochemistry Department, Max Planck Institute for Marine Microbiology, Bremen, Germany
Nitrification rates are low in permeable intertidal sand flats such that the water column is the primary source of nitrate to the sediment. During tidal inundation, nitrate is supplied to the pore space by advection rather than diffusion, relieving the microorganisms that reside in the sand from nitrate limitation and supporting higher denitrification rates than those observed under diffusive transport. Sand flats are also home to an abundant community of benthic photosynthetic microorganisms, the microphytobenthos (MPB). Diatoms are an important component of the MPB that can take up and store high concentrations of nitrate within their cells, giving them the potential to alter nitrate availability in the surrounding porewater. We tested whether nitrate uptake by the MPB near the sediment surface decreases its availability to denitrifiers along deeper porewater flow paths. In laboratory experiments, we used NOx (nitrate + nitrite) microbiosensors to confirm that, in the spring, net NOx consumption in the zone of MPB photosynthetic activity was stimulated by light. The maximum potential denitrification rate, measured at high spatial resolution using microsensors with acetylene and nitrate added, occurred below 1.4 cm, much deeper than light-induced NOx uptake (0.13 cm). Therefore, the shallower MPB had the potential to decrease NOx supply to the deeper sediments and limit denitrification. However, when applying a realistic downward advective flow to sediment from our study site, NOx always reached the depths of maximum denitrification potential, regardless of light availability or season. We conclude that during tidal inundation porewater advection overwhelms any influence of shallow NOx uptake by the MPB and drives water column NOx to the depths of maximum denitrification potential.
Introduction
Coastal ecosystems receive massive inputs of nitrogen from land that can fuel algal primary production and cause harmful algal blooms and hypoxia (Rabalais, 2002; Gruber and Galloway, 2008). Denitrification removes bioavailable nitrogen from the biosphere, returning it to the atmosphere as inert nitrogen gas (N2). Denitrification rates are usually limited by the rate of nitrate supply by diffusion or advection (Aller, 1982; Joye and Anderson, 2008). In permeable coastal sediments, advection is driven by a combination of faunal bioirrigation and current-driven pressure gradients over ripples, transporting oxygen and nitrate centimeters deep into the sand (Huettel et al., 1996; Billerbeck et al., 2006; Werner et al., 2006; Ahmerkamp et al., 2017). As a result, advection can increase areal denitrification rates by an order of magnitude compared with diffusive transport in the same sediment (Gao et al., 2012). However, not only denitrifying microorganisms consume nitrate in these sediments. The microphytobenthos (MPB) also competes for available nitrate (Lomstein et al., 1990; Sundbäck et al., 1991; Lorenzen et al., 1998; Longphuirt et al., 2009; Stief et al., 2013). Diatoms are often an important component of the MPB (Longphuirt et al., 2009; Stief et al., 2013). They are known to take up and store nitrate at high concentrations within their cells (Lomas and Glibert, 2000; Kamp et al., 2011, 2015; Merz et al., 2020). The influence of the MPB on denitrification in coastal sediment has been investigated in the past under conditions where mass transport was driven by diffusion rather than advection (Nielsen et al., 1990; Risgaard-Petersen et al., 1994; Lorenzen et al., 1998; Sundbäck and Miles, 2000; Risgaard-Petersen, 2003). In general, algal ammonium uptake led to nitrogen limitation of nitrification, which in turn provided less nitrate substrate for denitrification, thereby slowing rates of coupled nitrification-denitrification (Risgaard-Petersen, 2003).
We studied how nitrate uptake by the MPB affects denitrification in an intertidal sand flat with permeable sediments and porewater advection during tidal inundation. Here, low nitrification rates limit the importance of coupled nitrification-denitrification (Marchant et al., 2014), hence denitrification is controlled by the nitrate supply from the water column. While the nitrate supply is ample in the winter when water column nitrate concentrations are typically in excess of 50 μmol L–1, it is possible that nitrate is limiting in the summer when water column concentrations are 1–2 μmol L–1 or less (Grunwald et al., 2010). Therefore, denitrification is most likely limited by nitrate availability in the summer. The decrease from winter maximum to summer minimum water column nitrate concentrations is driven by spring phytoplankton blooms, some of which are deposited on shallow sediments (Lomstein et al., 1990). This corresponds with the end of the time period, between fall and spring, when phytopigment and intracellular nitrate concentrations are highest in intertidal sand flats (Stief et al., 2013). Therefore, intense springtime nitrate uptake by the MPB may also decrease porewater nitrate concentrations. In summary, it is most likely that MPB nitrate uptake can limit the amount of nitrate available for denitrification during the spring and summer.
We tested the hypothesis that, at high tide when water column nitrate is driven into the sediment by advection, the MPB intercepts nitrate near the sediment surface and decreases its availability for denitrifying microorganisms located deeper into the sand. In laboratory experiments, we used NOx (nitrate + nitrite) and oxygen microsensor measurements to test our hypothesis. We measured the depths in the sediment at which net oxygen production in the light and potential denitrification rates were highest to show a vertical separation between the two processes. We then measured oxygen and NOx penetration depths in the light and the dark across a range of imposed nitrate fluxes into the sediment to directly demonstrate whether nitrate uptake activity at the sediment surface was sufficient to alter NOx availability deeper in the sediment.
Materials and Methods
Field Sampling
The Janssand intertidal sand flat is located in the German Wadden Sea, between the mainland and Spiekeroog Island. The site is characterized by fine quartz sand (mean grain size 176 μm, high permeabilities of 7.2–9.5 × 10–12 m2, and low sediment total organic carbon content between 0.04 and 0.11% (Billerbeck et al., 2006), and a porosity of 0.38 (Gao et al., 2012). All samples were collected from the upper portion of the flat near its eastern edge (53°44.236′N, 7°41.818′E). With a typical tidal range of around 3 m, during high tide the upper flat is inundated by 1.5–2 m of water while during low tide it is exposed for 6–8 h (Billerbeck et al., 2006). Sediment cores were collected from areas with dense brown diatom mats. On a preliminary sampling trip in March 2016, a 1 cm3 sediment subsample was collected from the top 2 cm of each core and mixed with 5 mL of fixative (2% formaldehyde, 1% potassium iodine in nitrate-free artificial seawater (ASW)) and stored at 4°C for diatom cell counting in 6 mL exetainers. The samples were mixed by inversion, allowed to settle, and the supernatant was decanted. The sediment was rinsed in this manner three more times with 2 mL of ASW, and the supernatant was combined. The number of diatoms in ten 1 μL volumes was counted using an inverted microscope at 320x magnification (Lund et al., 1958). The mean diatom abundance was 28 ± 7 × 104 cells cm–3 (n = 22), an order of magnitude higher than previously measured at this site (Marchant et al., 2014), indicating that it was possible to visually target areas with high diatom abundance for sampling.
The sediment cores used in the following measurements were collected during tidal flat exposure in April and August of 2016 and February of 2017, when water temperatures were 9, 21, and 3°C, respectively. Around 15 cm of sediment was collected in clear, 8.3 cm diameter and 30 cm tall core tubes. In each season, 3 sediment cores were collected haphazardly from within a ∼0.5 m2 area with visible microphytobenthos and no visible macrofauna or bioturbation activity. NOx concentrations (ΣNO2–, NO3) were determined using a chemiluminescence detector after reduction to NO with 90°C acidic Vanadium (III) chloride (Braman and Hendrix, 1989). In situ NOx concentrations were 15, 1.5, and 10 μmol L–1, in April 2016, August 2016, and February 2017, respectively. Sediment cores were immediately transported in a dark cooler to the laboratory. They were stored at in situ temperatures under an automated daily light/dark illumination cycle that was programmed to match the time of year at the field site, with a light intensity of 100 μmol photons m–2 s–1. This light intensity was close to in situ values during inundation (Billerbeck et al., 2007). Seawater collected from the study site was percolated down through each core twice per day (every 12 h, ∼8 a.m. and 8 p.m.) by opening a valve at the bottom of the core and gravity-draining the water. This replaced the stagnant porewater in the top 5 cm of the sediment, simulating natural, tidal patterns of porewater flushing that regularly expose the microbial community to oxygen and nitrate (Werner et al., 2006). Around 0.5 cm of seawater was maintained above the sediment surface at all times to prevent the core from drying out. Seawater NOx concentrations were monitored and kept at in situ values for all experiments. Experiments were completed between 1 day and 2 weeks after sample collection. This time period was chosen to balance the need for replication of time-consuming microsensor measurements with potential changes in the sediment microbial community. The careful control of laboratory conditions described above with respect to temperature, light/dark cycles, and porewater flushing were implemented to minimize these changes, but the biomass and composition of the microbial community were not tracked during these experiments. The use of clear core tubes may have allowed algal growth on the internal core tube walls, though no visible growth was observed in any sediment core. Moreover, all microsensor measurements were made in the center of these 8.3 cm diameter core tubes, at least 3 cm from the nearest tube wall, minimizing any influence on porewater chemistry at the measurement location.
All measurements made using these sediment cores, described in detail below, were completed during the daylight phase of the artificial day/night cycle, typically beginning 2–4 h after dawn. Because experimental manipulations involved percolating surface seawater through the core, they took the place of the regular, morning (8 am) percolation event meant to simulate tidal inundation. Therefore, measurements typically began 12–14 h following the previous (8 pm) percolation event. Because simulated tidal inundation in the laboratory was carried out by hand, it was not synchronized with the timing of tidal inundation at the field site as this would have required large shifts in the relative timing of light and inundation over a 2-week measurement period.
Microsensor Measurements of Volumetric NOx Consumption Rates
Oxygen and NOx consumption rates were measured using oxygen microsensors (Revsbech, 1989) and NOx microbiosensors as described previously (Larsen et al., 1996), using a microsensor set up with a motorized micromanipulator and computerized data acquisition. Initially, each sediment core was placed in the light (∼100 μmol photons m–2 s–1). The sediment surface with respect to the sensor tips was determined using a dissecting microscope. Two consecutive oxygen and NOx microprofiles were then measured ∼20 min apart to confirm that the profiles were at steady state. At 0.1–0.2 cm depth an oxygen maximum was observed in these steady state microprofiles, indicating the maximum photosynthetic activity.
The microsensors were then positioned at the depth of the oxygen maximum. Oxygenated overlying seawater containing ambient levels of NOx was rapidly gravity-drained through the core, until the oxygen and NOx concentrations at the target depth matched those in the overlying water. Water flow was then turned off, and volumetric NOx and oxygen production and consumption rates were obtained from linear regression of concentration versus time during the 10 min measurement window. All regressions were significant (p < 0.01) with r2-values greater than 0.95. These measurements were made on two independent cores each in the spring and summer.
Depth Profiles of Potential Denitrification Rates
Depth profiles of potential denitrification rates in the sediment were measured using a modified version of the acetylene block technique (Sørensen, 1978) in one core each in the summer and the winter. Acetylene blocks the last step in denitrification, the reduction of nitrous oxide (N2O) to N2, allowing N2O to accumulate (Yoshinari and Knowles, 1976). This technique was combined with the method for measuring oxygen consumption rate profiles in permeable sediment (Polerecky et al., 2005) by adding a N2O microsensor (Andersen et al., 2001). A volume of 160 mL of oxygenated site seawater, amended with nitrate to 20 μmol L–1, was added to the core above the sediment. The overlying water was gravity-drained down through the core in 3–4 min. Subsequently, repeated oxygen and N2O profiles from 0 to 4 cm depth at 0.1 cm resolution were measured approximately every 18 min. This process was repeated on the same core, this time with 20% (32 mL) of the overlying seawater volume (160 mL) replaced with acetylene-saturated seawater. Complete inhibition of N2O reduction has been observed at an acetylene concentration of 0.1 atmospheres (Yoshinari and Knowles, 1976), which can be achieved by replacing 10% of a given seawater volume with acetylene-saturated seawater (Seitzinger et al., 1993). Repeated oxygen and N2O microprofiles were measured as described above until consecutive microprofiles were nearly identical (around 2.5 h in the summer and 10 h in the winter).
Denitrification rates were calculated from the N2O profiles by solving Equation 1 for R
where dC/dt is the change in N2O concentration with respect to time, R is the denitrification rate, and Jout and Jin are the diffusive fluxes of N2O out of and into the volume of sediment centered on the point where the concentration measurement was made. dC/dt at each depth was calculated as the slope of the linear regression of N2O concentration versus time. A minimum of 3 N2O measurements was used to calculate dC/dt. The diffusive fluxes were calculated using Equation 2 as described by Deurer et al. (2008).
where φ is the dimensionless sediment porosity (0.38), Dw is the diffusion coefficient for N2O in water at 20°C, and dC/dz is the N2O concentration gradient with respect to depth. We assumed a diffusion coefficient of 6.58 × 10–6 m2 h–1 (Tamimi et al., 1994). dC/dz was estimated as the difference in concentration between the depth where the calculation was made and the depth above it (Jout) or below it (Jin) divided by the change in depth (0.001 m). Since multiple N2O profiles were measured during the time that dC/dt was constant, the net N2O flux out of each depth (Jout – Jin) was calculated for each profile and the mean flux was used to calculate denitrification rates. Denitrification rates are expressed on a per mole of nitrogen (not N2O) basis.
Potential Denitrification Rate Kinetics
This modified acetylene block technique was also used to estimate denitrification reaction kinetics. To do this, denitrification rates were calculated at each depth and time point where an N2O concentration measurement was made. This rate was estimated at each depth as the first derivative of N2O concentration versus time. The first derivative at each time was approximated as the slope of the curve’s tangent line at that point following Equation 3 below (Wood, 1982).
At any given depth in the profile, dC/dtn is the change in N2O concentration with time in profile n, N2O(n+1) is the N2O concentration measured in the next profile (profile n + 1), N2O(n–) is the N2O concentration measured in the preceding profile (profile n−1), t(n+1) is the time at which the N2O measurement was made in profile n + 1, and t(n–1) is the time at which the N2O measurement was made in profile n−1. Denitrification rates were calculated from dC/dtn using Equation 1. These denitrification rates are calculated for a single point in time. The rates measured in this section are referred to as instantaneous denitrification rates throughout the manuscript to differentiate them from the rates calculated in section “Depth Profiles of Potential Denitrification Rates.”
The apparent porewater NOx concentration was estimated at each depth and time point for comparison with the instantaneous denitrification rates. These concentrations were calculated based on the total amount of N2O produced and the original 20 μmL L–1 pool of NOx available in the porewater according to Equation 4.
Where NOx(app) is the apparent NOx concentration remaining in the porewater at a given time and depth, NOx(init) is the initial porewater NOx concentration (20 μmol L–1) added at the start of the incubation, N2O is the measured N2O concentration (multiplied by 2 to account for denitrification reaction stoichiometry), and f is a correction factor that accounts for the fraction of NOx that was consumed by processes other than denitrification (calculated according to Equation 5 below).
Where N2Omax is the maximum N2O concentration measured at a particular depth during any time point. This calculation assumes that all or nearly all the initial NOx available was consumed by some combination of processes over the course of the ∼2.5-h incubation.
Finally, denitrification reaction kinetics were calculated by fitting the instantaneous denitrification rate versus apparent porewater NOx concentration data using the Hill equation (Hill, 1913) (Equation 6 below).
Where R is the instantaneous denitrification rate derived from Equations 3 and 1 and NOx(app) is the apparent NOx concentration derived from Equation 4. Vmax is the maximum possible rate, KM is the half saturation constant, and n is the Hill coefficient. The Hill coefficient is a representation of the relative degree of interaction, or cooperativity, between multiple substrate molecules binding with an enzyme (Monod et al., 1965). The Hill equation has been used to describe the kinetics of allosteric enzymes (Weiss, 1997), and has recently been adapted to model microbial growth rates under substrate-limiting conditions (Arranz and Peinado, 2017). The Hill equation was used instead of the simpler Michaelis-Menten equation (Michaelis and Menten, 1913, 2013) because it yielded a better fit to our data (Akaike information criterion (AIC) values of 39.2 and 90.5, respectively).
Influence of the MPB on the Depth of NOx Penetration Into the Sediment
The MPB microbial community is heterogeneously yet broadly distributed across the intertidal sand flats at this study site (Chennu et al., 2013), making it impossible to collect control sediment cores that would be devoid of MPB activity. Instead, we opted to focus on sediments with the highest MPB density, identified as described in section “Field Sampling.” Our experimental approach centered on creating conditions that maximized the potential of the MPB to alter NOx availability at the depth of maximum denitrification activity. If NOx penetrated to that depth under these idealized conditions, then we consider it very unlikely that the MPB had any substantial influence on substrate availability for denitrification under either laboratory or in situ conditions.
NOx and oxygen penetration depths were measured using oxygen microsensors and NOx microbiosensors. An initial microprofile was measured with no porewater flow. Then overlying seawater containing oxygen and ambient NO3 was pumped down through the core, at either low or high flow rates. Low flow was defined as a porewater infiltration rate of 0.11 mL min–1 in all experiments, high flow was 0.36 mL min–1 in the spring and 1.3 mL min–1 in the summer. These porewater flow rates were equivalent to vertical seawater infiltration rates of 1.2, 4.0, and 14.4 L m–2 h–1, respectively. The imposed NOx flux into the sediment was 18–60 μmol m–2 h–1 in the spring and 2–22 μmol m–2 h–1 in the summer. All measurements were done in both the dark and at a light intensity of 100 μmol photons m–2 s–1, to test the influence of MPB nitrate uptake on nitrate penetration.
Continuous microprofiling began immediately after the onset of downward porewater flow and continued for 3–6 h with consecutive microprofiles measured at 30–60-min intervals. For each microprofile measured, the NOx and oxygen penetration depths were calculated as the first depth at which the concentration reached zero. The penetration depth appeared to gradually increase with ongoing percolation. Therefore, the penetration depth was plotted against the elapsed time since the onset of porewater flow and fit with an asymptotic curve of the form:
where d is the penetration depth, t is the time elapsed since the start of pumping, and b and c are dimensionless coefficients used to fit the curve, and a is the estimate of the asymptote. This asymptote is the maximum penetration depth, estimated using a least squares fit of Equation 7 (calculated using the non-linear least squares function in R) to the experimentally derived penetration depth versus time curves.
Results
NOx Consumption Rates by MPB
In the zone of highest photosynthetic activity, NOx consumption rates were higher in the light than in the dark, but only in the spring (Figure 1A). In spring, a switch from the dark to the light increased net NOx consumption rates from 13 mmol N (m3 sediment) –1 h–1 to 28 mmol N (m3 sediment) –1 h–1 in one core (a 1.8-fold increase) and from 23 mmol N (m3 sediment) –1 h–1 to 31 mmol N (m3 sediment) –1 h–1 in a duplicate core (a 1.3-fold increase). In contrast, in summer, net NOx consumption was not observed in either the light or the dark (Figure 1B), even when nitrate levels were increased to far above ambient summer concentrations (from 1.5 to 20 μmol L–1). Net oxygen production and consumption rates were an order of magnitude higher in the summer than the spring (Figures 1C,D).
Figure 1. Changes in net NOx consumption rates between dark and light conditions measured in the spring (A) and summer (B). Changes in net oxygen production rates between dark and light conditions measured in the spring (C) and summer (D). A line connects the dark and light measurements made on each core. Negative net consumption rates in (B) indicate net NOx production.
Depth Separation of Photosynthesis and Potential Denitrification
In the summer, maximum denitrification occurred far below the zone where the MPB were photosynthetically active (Figure 2). Whereas the photosynthetic oxygen maximum was found at 0.13 cm depth, denitrification peaked below 1.4 cm. Denitrification also occurred within the depth interval that was always exposed to oxygen (0–0.4 cm), but at a much lower level. The mean (± standard deviation) denitrification rate across 0–0.4 cm depth was 0.95 ± 0.19 mmol N (m3 sediment) –1 h–1. The potential denitrification rate increased steadily with depth down to about 1.4 cm after which it stabilized at values from 5.9 to 7.4 mmol N (m3 sediment) –1 h–1. The mean denitrification rate over this depth range (1.4–3.9 cm) was 6.9 ± 0.4 mmol N (m3 sediment)–1 h–1. In the winter the potential denitrification rates were lower, particularly deeper in the core. However, the pattern of low potential denitrification rates near the sediment surface which increased to a maximum by 1.5–2 cm depth was identical in both seasons (Figure 2B).
Figure 2. Representative summer steady state oxygen concentration depth profiles in the light (open circles) and in the dark (closed circles) (A). Potential denitrification rate profiles in the summer (black) and in the winter (gray) measured in the light (B). Error bars represent the linear regression standard error. The gap in the summer profile reflects two separate measurements made over different depth intervals in the same core. The shaded areas in (A,B) represent the depth interval over which oxygen was always present in steady state profiles. Repeated N2O concentration microprofiles following acetylene addition in the summer (C). The dashed line represents the maximum possible N2O concentration that could be generated from the available porewater NOx. The instantaneous potential denitrification rate plotted against the apparent amount of NOx remaining in the porewater (D) over the 1.4–2.0 cm depth range highlighted in (C). The black line is fit to the data using Equation 6.
Potential Denitrification Reaction Kinetics
Potential denitrification reaction kinetics were calculated from the same denitrification rate profiles (measured in summer in the light) described in section “Depth Separation of Photosynthesis and Potential Denitrification.” Focusing on a depth interval with stable, high denitrification rates (1.4–2.0 cm), it is apparent that when acetylene blocked the last step in denitrification, N2O accumulated rapidly during the first hour of incubation and slowly approached a maximum value thereafter (Figure 2C). The maximum amount of N2O that could be produced was 10 μmol L–1 based on the 20 μmol L–1 NOx added to the seawater medium. The mean N2O concentration across these depths at the last time point measured was 9.1 ± 0.1 μmol L–1, or around 91% of the maximum possible. Therefore, at the depth of maximum potential denitrification rates, more than 90% of the available NOx was consumed by denitrification.
Instantaneous potential denitrification rates increased non-linearly as a function of apparent porewater NOx concentrations (Figure 2C). These rates are called “instantaneous” to distinguish them from the potential denitrification rates from the preceding paragraph that were calculated as the change in N2O concentration over time at each depth. Instead, instantaneous rates were estimated at a given point in time from the derivatives of individual N2O concentration profiles. The non-linear increase in instantaneous potential denitrification rates as a function of substrate concentration was modeled using the Hill equation (Equation 6). Fitting this curve to the data yielded a half-saturation constant (KM) of 3.1 μmol L–1, a maximum instantaneous denitrification rate (Vmax) of 6.8 mmol N (m3 sediment) –1 h–1, and a Hill coefficient of 2.4. The upper inflection in the curve in Figure 2C occurs at a rate of around 6.1 mmol N (m3 sediment) –1 h–1, which is 90% of Vmax, and corresponds with an apparent NOx concentration of ∼7.8 μmol L–1. Above this NOx concentration, the instantaneous denitrification rate increased very slowly as a function of increasing NOx concentration.
Influence of the MPB on the Depth of NOx Penetration Into the Sediment
When overlying seawater containing oxygen and NOx was pumped down through a sediment core at a constant rate, these electron acceptors penetrated gradually deeper but their concentration gradients consistently remained steep (see representative profiles in Figure 3). There was ample NOx available in the porewater just a few millimeters above the NOx penetration depth regardless of the amount of time elapsed since the onset of pumping. Therefore, the NOx penetration depth serves as a proxy for the depth range over which denitrification was possible.
Figure 3. Time series of oxygen (A) and NOx (B) depth profiles after onset of downward water flow through a sediment core. These representative profiles are from the spring, in the dark, with a low flow rate of 0.11 mL min–1.
NOx penetration depths in the light and dark were directly compared in three separate cores, once at low flow in the spring, once at high flow in the spring, and once at high flow in the summer. In all cases there was no difference in NOx penetration between the light and dark treatments (Figures 4A,B). However, flow rate did have an effect; NOx penetration depths were much higher at high flow rates than at low flow rates in both the spring (Figure 4A; different flow rates measured in independent cores) and the summer (Figure 4B; measured in the same core). Similarly, flow rate had a much stronger influence on the NOx saturation depth than light (Figure 4C). The NOx saturation depth is the depth at which the porewater NOx concentration fell below the saturating concentration of 7.8 μmol L–1 required to drive potential denitrification rates at 90% of Vmax (see section “Potential Denitrification Reaction Kinetics” for details). Under all conditions tested in the spring, saturating NOx concentrations reached the depth of maximum denitrification potential (Figure 4C). This deep NOx penetration occurred within 4 h after the onset of porewater flow.
Figure 4. The change in NOx penetration depths (A,B), NOx saturation depths (the depth at which the NOx concentration falls below 7.8 μmol L–1) (C), and oxygen penetration depths (D,E) as a function of time following the onset of porewater flow. Data from the spring are displayed in the left column (A,C,D) and data from the summer are displayed in the right column (B,E). All summer profiles were measured in the same core. In the spring, low flow rate profiles were measured in one core and high flow rate profiles in a second, independent core. The lines represent asymptotic fits to the data according to Equation 7. The gray areas in (A–C) represent the deep depth interval over which potential denitrification rates reached their maximum (below 1.4 cm; Figure 2B).
As with NOx penetration depths, oxygen penetration depths were similar at the same flow rate regardless of light availability in both spring (Figure 4D) and summer (Figure 4E). The oxygen penetration depth was consistently shallower than the NOx penetration depth during porewater advection (mean ± standard deviation = 0.8 ± 0.4 cm) across all individual profile measurements in all experimental conditions (n = 20 profiles after time 0 in Figure 4).
Discussion
Light-Induced NOx Uptake by the MPB
In the spring, NOx was consumed more rapidly in the light than in the dark where the MPB were photosynthetically active (Figure 1A). Stimulation of net NOx consumption in the light is consistent with nitrate uptake and storage by diatoms grown in the light (Kamp et al., 2011), and with assimilation into biomass (Middelburg and Nieuwenhuize, 2000). The MPB may further stimulate bacterial nitrogen uptake in the light by exuding organic carbon that can serve as a substrate for heterotrophic activity (Evrard et al., 2008). NOx assimilation was more important than dissimilatory respiration in shallow, oxic sediment during the spring. Near the sediment surface, light-induced NOx uptake (the difference between net NOx uptake measured in the light and the dark) accounted for 29–33% of total net NOx consumption. In contrast, denitrification in oxic sediment in any season only accounted for at most 3–4% of total net NOx consumption. As DNRA rates at this study site were less than 20% of denitrification rates (Marchant et al., 2014) these were not considered further. The remaining balance of the observed net NOx uptake can most likely be attributed to microbial assimilation, by either autotrophs or heterotrophs, in the dark (e.g., Middelburg and Nieuwenhuize, 2000). In summary, dissimilatory NOx respiration was less important than MPB-influenced NOx uptake near the sediment surface in the spring. No net NOx consumption was observed in photosynthetic sediment during the summer in either the light or the dark, even when exposed to an unseasonably high concentration of 20 μmol L–1 (Figure 1A). This was despite very high net oxygen production rates (Figure 1B). The most likely explanation for this is that high rates of organic matter remineralization in the summer (Billerbeck et al., 2006) released sufficient ammonium into the porewater to meet the nitrogen demand of the MPB.
Denitrification Kinetics
Our half-saturation constant (KM) of 3.1 μmol L–1 and maximum instantaneous denitrification rate (Vmax) of 6.8 mmol N (m3 sediment) –1 h–1 fell within the ranges of 1.5–19.8 μmol L–1 and 0.9–7.5 mmol N (m3 sediment) –1 h–1 reported for permeable sediments across six sites in coastal Australia (Evrard et al., 2013). However, we observed a low sensitivity of denitrification rates to increasing NOx at low concentrations (Figure 2B) that is not apparent in the data of Evrard et al. (2013). It is possible that at very low porewater NOx concentrations, denitrification rates are limited by diffusion-driven mass-transport of NOx from the porewater to the surface of sand grains inhabited by microorganisms (Probandt et al., 2018; Ahmerkamp et al., 2020) that are likely the dominant sites of denitrification within the bulk sediment.
Controls on NOx Penetration Depth During Tidal Inundation
Despite high nitrate consumption activity in the top few millimeters of sediment (Figure 1A), when nitrate-bearing seawater was pumped into sediment cores to simulate in situ porewater advection (n = 7), some NOx always penetrated to depths with high denitrification potential regardless of season, light availability, or porewater flow rate (Figures 4A,B and Table 1). NOx saturation depths were at most 0.5 cm shallower than NOx penetration depths. Moreover, saturating NOx concentrations (7.8 μmol L–1) penetrated more than 1 cm deep into the sand under all conditions tested. This was the case even though these experiments were designed to maximize the potential influence of the MPB. The range in experimental NOx fluxes into the sediment was 2–60 μmol m–2 h–1, at least an order of magnitude lower than the 700–43,000 μmol N m–2 h–1 range of nitrate fluxes into the sediment that have been modeled for this site (Gao et al., 2012). Even so, the advective transport of nitrate into the sediment was always higher than nitrate consumption by the MPB such that the MPB could not limit nitrate supply to the denitrifying microbial community.
Table 1. Theoretical maximum NOx and oxygen penetration depths and NOx fluxes, under different light and flow conditions.
The observation that light-induced MPB activity did not influence the depth of nitrate penetration can be validated by a simple mass balance of nitrate in the zone of benthic photosynthetic activity. This zone extended over a depth range of at most 0.8 cm (the maximum steady state oxygen penetration depth; Figure 2A). Assuming that light induced an increase in the maximum net volumetric NOx consumption rate of 8 mmol N (m3 sediment) –1 h–1 (the difference between light and dark net NOx consumption rates in one core in the spring; Figure 2A) in this zone, the maximum increase in the areal net NOx consumption rate that can be attributed to the MPB was 64 μmol m–2 h–1. Even if we apply the maximum net NOx consumption rate measured in the spring (31 mmol N (m3 sediment) –1 h–1) to this entire depth zone, the areal net NOx consumption rate comes to 248 μmol m–2 h–1. This maximum possible rate is only 35% of the low end of the 700–43,000 μmol N m–2 h–1 range of nitrate fluxes into the sediment that have been modeled for this site (Gao et al., 2012). In summary, we find no evidence that the MPB limited the amount of nitrate available in the sediment for denitrification within our very conservative experimental framework, and we also have no reason to expect that the MPB would limit nitrate availability for denitrification in situ.
Coupling Between NOx and Oxygen Penetration Depths
For the porewater flow rates of 0.11, 0.36, and 1.3 mL min–1 used in this experiment, water would, in theory, penetrate to 2.6, 8.4, and 30 cm depth, respectively, over the 8 h represented in Figure 4. Measured NOx and oxygen penetration depths were much shallower than these values, demonstrating the importance of consumption throughout the sediment in controlling the penetration depth of terminal electron acceptors. During advection, oxygen and NOx penetrated easily through shallow, oxidized sand where consumption rates were low, but higher consumption rates in deeper, more reduced sands limited the depth to which they could penetrate. For example, potential denitrification rates were 6–7 times higher in anoxic than oxic conditions (Figure 2A). Therefore, NOx penetrated easily through the oxic zone of the sand, but was consumed completely within a few millimeters of the oxygen penetration depth.
Despite differences in penetration depth, oxygen and NOx coexisted over a depth range that was typically several centimeters deep during porewater advection (Figure 4). The presence of oxygen does not prevent this NOx from being used for denitrification. We measured denitrification rates in the oxic zone similar to those measured in other oxic permeable intertidal and subtidal sediments (Rao et al., 2008; Gao et al., 2010), and denitrifying microorganisms from our study site are capable of simultaneously respiring oxygen and nitrate (Marchant et al., 2017). Moreover, porewater oxygen concentrations ranged between 20 μmol L–1 to more than 200 μmol L–1 in the oxic zone of our sediment cores in which we detected potential denitrification, similar to the > 100 μmol L–1 oxygen background in which denitrification has been observed by others in whole sediment (Marchant et al., 2017) and bacterial cultures (Qu et al., 2016). Therefore, if oxygen and NOx profiles overlap substantially during sand flat inundation as they did in our experiments, then denitrification in the oxic zone may play a very important role in the overall nitrogen metabolism of the flat.
Conclusion
Based on our results, we can reject the hypothesis that MPB covering permeable intertidal sandflats intercepts nitrate near the sediment surface and decreases its availability for denitrifying microorganisms. Though MPB-inhabited sediments do indeed take up nitrate, this uptake is overwhelmed by advective nitrate transport. Therefore, it appears that MPB limitation of denitrification only occurs in sediments where diffusion is the dominant transport process (Lorenzen et al., 1998; Risgaard-Petersen, 2003). Previous studies have found no systematic differences in MPB abundance between diffusive and advective sediments that could account for this difference in the influence of MPB on denitrification (Sundbäck and Miles, 2000; Cartaxana et al., 2006). Our findings emphasize that intertidal permeable sediments are important sites of nitrogen removal via denitrification (Gao et al., 2012; Huettel et al., 2014; Marchant et al., 2014).
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The data and associated computer code are publicly available as a computing capsule at the Code Ocean repository (doi: 10.24433/CO.1921276.v1).
Author Contributions
CS led the data analysis and writing efforts. CS and PH-V carried out the experiments. CS, PH-V, and HM analyzed the samples. All authors contributed to the experimental design and to writing and editing the manuscript.
Funding
Funding for this work was provided by the Max Planck Society.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Gabi Eickert, Karin Hohmann, Vera Hübner, Anja Niclas, Ines Schröder, and Cäcilia Wigand for building the microsensors used in this study. We would also like to thank Manfred Schlösser, Tomas Wilkop, Anja Niclas, and Cäcilia Wigand for assisting with sample collection in the Wadden Sea. Finally, we thank Ugo Marzocchi for his advice in making the denitrification rate measurements, and Gabi Klockgether for her assistance with nitrogen isotope analyses.
References
Ahmerkamp, S., Marchant, H. K., Peng, C., Probandt, D., Littmann, S., Kuypers, M. M., et al. (2020). The effect of sediment grain properties and porewater flow on microbial abundance and respiration in permeable sediments. Sci. Rep. 10, 1–12. doi: 10.1016/j.ecss.2011.10.024
Ahmerkamp, S., Winter, C., Krämer, K., De Beer, D., Janssen, F., Friedrich, J., et al. (2017). Regulation of benthic oxygen fluxes in permeable sediments of the coastal ocean. Limnol. Oceanogr. 62, 1935–1954. doi: 10.1002/lno.10544
Aller, R. C. (1982). “The effects of macrobenthos on chemical properties of marine sediment and overlying water,” in Animal-Sediment Relations, eds P. L. McCall and M. J. S. Tevesz (New York, NY: Springer), 53–102. doi: 10.1007/978-1-4757-1317-6_2
Andersen, K., Kjaer, T., and Revsbech, N. P. (2001). An oxygen insensitive microsensor for nitrous oxide. Sens. Actuat. B 81, 42–48. doi: 10.1016/s0925-4005(01)00924-8
Arranz, F., and Peinado, J. (2017). A mesoscopic stochastic model for the specific consumption rate in substrate-limited microbial growth. PLoS One 12:e0171717. doi: 10.1371/journal.pone.0171717
Billerbeck, M., Roy, H., Bosselmann, K., and Huettel, M. (2007). Benthic photosynthesis in submerged Wadden Sea intertidal flats. Estuar. Coast. Shelf Sci. 71, 704–716. doi: 10.1016/j.ecss.2006.09.019
Billerbeck, M., Werner, U., Polerecky, L., Walpersdorf, E., De Beer, D., and Huettel, M. (2006). Surficial and deep pore water circulation governs spatial and temporal scales of nutrient recycling in intertidal sand flat sediment. Mar. Ecol. Prog. Ser. 326, 61–76. doi: 10.3354/meps326061
Braman, R. S., and Hendrix, S. A. (1989). Nanogram nitrite and nitrate determination in environmental and biological materials by vanadium(III) reduction with chemiluminescence detection. Analytic. Chem. 61, 2715–2718. doi: 10.1021/ac00199a007
Cartaxana, P., Mendes, C., Van Leeuwe, M., and Brotas, V. (2006). Comparative study on microphytobenthic pigments of muddy and sandy intertidal sediments of the Tagus estuary. Estuar. Coast. Shelf Sci. 66, 225–230. doi: 10.1016/j.ecss.2005.08.011
Chennu, A., Färber, P., Volkenborn, N., Al-Najjar, M. A., Janssen, F., De Beer, D., et al. (2013). Hyperspectral imaging of the microscale distribution and dynamics of microphytobenthos in intertidal sediments. Limnol. Oceanogr. Methods 11, 511–528. doi: 10.4319/lom.2013.11.511
Deurer, M., Von Der Heide, C., Böttcher, J., Duijnisveld, W. H. M., Weymann, D., and Well, R. (2008). The dynamics of N2O near the groundwater table and the transfer of N2O into the unsaturated zone: a case study from a sandy aquifer in Germany. Catena 72, 362–373. doi: 10.1016/j.catena.2007.07.013
Evrard, V., Cook, P. L., Veuger, B., Huettel, M., and Middelburg, J. J. (2008). Tracing carbon and nitrogen incorporation and pathways in the microbial community of a photic subtidal sand. Aquat. Microb. Ecol. 53, 257–269. doi: 10.3354/ame01248
Evrard, V., Glud, R. N., and Cook, P. L. M. (2013). The kinetics of denitrification in permeable sediments. Biogeochemistry 113, 563–572. doi: 10.1007/s10533-012-9789-x
Gao, H., Matyka, M., Liu, B., Khalili, A., Kostka, J. E., Collins, G., et al. (2012). Intensive and extensive nitrogen loss from intertidal permeable sediments of the Wadden Sea. Limnol. Oceanogr. 57, 185–198. doi: 10.4319/lo.2012.57.1.0185
Gao, H., Schreiber, F., Collins, G., Jensen, M. M., Kostka, J. E., Lavik, G., et al. (2010). Aerobic denitrification in permeable Wadden Sea sediments. ISME J. 4, 417–426. doi: 10.1038/ismej.2009.127
Gruber, N., and Galloway, J. N. (2008). An Earth-system perspective of the global nitrogen cycle. Nature 451, 293–296. doi: 10.1038/nature06592
Grunwald, M., Dellwig, O., Kohlmeier, C., Kowalski, N., Beck, M., Badewien, T. H., et al. (2010). Nutrient dynamics in a back barrier tidal basin of the Southern North Sea: time-series, model simulations, and budget estimates. J. Sea Res. 64, 199–212. doi: 10.1016/j.seares.2010.02.008
Hill, A. V. (1913). The combinations of haemoglobin with oxygen and with carbon monoxide. I Biochem. J. 7:471. doi: 10.1042/bj0070471
Huettel, M., Berg, P., and Kostka, J. E. (2014). Benthic exchange and biogeochemical cycling in permeable sediments. Annu. Rev. Mar. Sci. 6, 23–51. doi: 10.1146/annurev-marine-051413-012706
Huettel, M., Ziebis, W., and Forster, S. (1996). Flow-induced uptake of particulate matter in permeable sediments. Limnol. Oceanogr. 41, 309–322. doi: 10.4319/lo.1996.41.2.0309
Joye, S. B., and Anderson, I. C. (2008). “Nitrogen cycling in coastal sediments,” in Nitrogen in the Marine Environment, 2nd Edn, eds D. G. Capone, D. A. Bronk, M. R. Mulholland, and E. J. Carpenter (New York: Elsevier), 686–902.
Kamp, A., De Beer, D., Nitsch, J. L., and Stief, P. (2011). Diatoms respire nitrate to survive dark and anoxic conditions. Proc. Natl. Acad. Sci. U.S.A. 108, 5649–5654. doi: 10.1073/pnas.1015744108
Kamp, A., Høgslund, S., Risgaard-Petersen, N., and Stief, P. (2015). Nitrate storage and dissimilatory nitrate reduction by eukaryotic microbes. Front. Microbiol. 6:1492. doi: 10.3389/fmicb.2015.01492
Larsen, L. H., Revsbech, N. P., and Binnerup, S. J. (1996). A microsensor for nitrate based on immobilized denitrifying bacteria. Appl. Environ. Microbiol. 62, 1248–1251. doi: 10.1128/aem.62.4.1248-1251.1996
Lomas, M. W., and Glibert, P. M. (2000). Comparisons of nitrate uptake, storage, and reduction in marine diatoms and flagellates. J. Phycol. 36, 903–913. doi: 10.1046/j.1529-8817.2000.99029.x
Lomstein, E., Jensen, M. H., and Sørensen, J. (1990). Intracellular NH4+ and NO3– pools associated with deposited phytoplankton in a marine sediment (Aarhus Bight, Denmark). Mar. Ecol. Pro. Ser. 61, 97–105. doi: 10.3354/meps061097
Longphuirt, S. N., Lim, J.-H., Leynaert, A., Claquin, P., Choy, E.-J., Kang, C.-K., et al. (2009). Dissolved inorganic nitrogen uptake by intertidal microphytobenthos: nutrient concentrations, light availability and migration. Mar. Ecol. Pro. Ser. 379, 33–44. doi: 10.3354/meps07852
Lorenzen, J., Larsen, L. H., Kjaer, T., and Revsbech, N.-P. (1998). Biosensor determination of the microscale distribution of nitrate, nitrate assimilation, nitrification, and denitrification in a diatom-inhabited freshwater sediment. Appl. Environ. Microbiol. 64, 3264–3269. doi: 10.1128/aem.64.9.3264-3269.1998
Lund, J. W. G., Kipling, C., and Le Cren, E. D. (1958). The inverted microscope method of estimating algal numbers and the statistical basis of estimations by counting. Hydrobiologia 11, 143–170. doi: 10.1007/bf00007865
Marchant, H. K., Ahmerkamp, S., Lavik, G., Tegetmeyer, H. E., Graf, J., Klatt, J. M., et al. (2017). Denitrifying community in coastal sediments performs aerobic and anaerobic respiration simultaneously. ISME J. 11, 1799–1812. doi: 10.1038/ismej.2017.51
Marchant, H. K., Lavik, G., Holtappels, M., and Kuypers, M. M. (2014). The fate of nitrate in intertidal permeable sediments. PLoS One 9:e104517. doi: 10.1371/journal.pone.00104517
Merz, E., Dick, G. J., De Beer, D., Grim, S., Hübener, T., Littmann, S., et al. (2020). Nitrate respiration and diel migration patterns of diatoms are linked in sediments underneath a microbial mat. Environ. Microbiol. 23, 1422–1435. doi: 10.1111/1462-2920.15345
Michaelis, L., and Menten, M. L. (1913). Die Kinetic der Invertin-wirkung. Biochemische Zeitschrift 49, 333–369.
Michaelis, L., and Menten, M. L. (2013). The kinetics of invertin action: translated by TRC Boyde submitted 4 February 1913. FEBS Lett. 587, 2712–2720. doi: 10.1016/j.febslet.2013.07.015
Middelburg, J. J., and Nieuwenhuize, J. (2000). Nitrogen uptake by heterotrophic bacteria and phytoplankton in the nitrate-rich Thames estuary. Mar. Ecol. Prog. Ser. 203, 13–21. doi: 10.3354/meps203013
Monod, J., Wyman, J., and Changeux, J.-P. (1965). On the nature of allosteric transitions: a plausible model. J. Mol. Biol. 12, 88–118. doi: 10.1016/s0022-2836(65)80285-6
Nielsen, L. P., Christensen, P. B., Revsbech, N. P., and Sørensen, J. (1990). Denitrification and photosynthesis in stream sediment studied with microsensor and whole-core techniques. Limnol. Oceanogr. 35, 1135–1144. doi: 10.4319/lo.1990.35.5.1135
Polerecky, L., Franke, U., Werner, U., Grunwald, M., and De Beer, D. (2005). High spatial resolution measurement of oxygen consumption rates in permeable sediments. Limnol. Oceanogr. 3, 75–85. doi: 10.4319/lom.2005.3.75
Probandt, D., Eickhorst, T., Ellrott, A., Amann, R., and Knittel, K. (2018). Microbial life on a sand grain: from bulk sediment to single grains. ISME J. 12, 623–633. doi: 10.1038/ismej.2017.197
Qu, Z., Bakken, L. R., Molstad, L., Frostegård, Å, and Bergaust, L. L. (2016). Transcriptional and metabolic regulation of denitrification in P aracoccus denitrificans allows low but significant activity of nitrous oxide reductase under oxic conditions. Environ. Microbiol. 18, 2951–2963. doi: 10.1111/1462-2920.13128
Rao, A. M. F., Mccarthy, M. J., Gardner, W. S., and Jahnke, R. A. (2008). Respiration and denitrification in permeable continental shelf deposits on the South Atlantic Bight: N2:Ar and isotope pairing measurements in sediment column experiments. Continent. Shelf Res. 28, 602–613. doi: 10.1016/j.csr.2007.11.007
Revsbech, N. P. (1989). An oxygen microsensor with a guard cathode. Limnol. Oceanogr. 34, 474–478. doi: 10.4319/lo.1989.34.2.0474
Risgaard-Petersen, N. (2003). Coupled nitrification-denitrification in autotrophic and heterotrophic estuarine sediments: on the influence of benthic microalgae. Limnol. Oceanogr. 48, 93–105. doi: 10.4319/lo.2003.48.1.0093
Risgaard-Petersen, N., Rysgaard, S., Nielsen, L. P., and Revsbech, N. P. (1994). Diurnal variation of denitrification and nitrification in sediments colonized by benthic microphytes. Limnol. Oceanogr. 39, 573–579. doi: 10.4319/lo.1994.39.3.0573
Seitzinger, S. P., Nielsen, L. P., Caffrey, J., and Christensen, P. B. (1993). Denitrification measurements in aquatic sediments: a comparison of three methods. Biogeochemistry 23, 147–167. doi: 10.1007/bf00023750
Sørensen, J. (1978). Denitrification rates in a marine sediment as measured by the acetylene inhibition technique. Appl. Environ. Microbiol. 36, 139–143. doi: 10.1128/aem.36.1.139-143.1978
Stief, P., Kamp, A., and De Beer, D. (2013). Role of diatoms in the spatial-temporal distribution of intracellular nitrate in intertidal sediment. PLoS One 8:e73257. doi: 10.1371/journal.pone.0073257
Sundbäck, K., Enoksson, V., Granéli, W., and Pettersson, K. (1991). Influence of sublittoral microphytobenthos on the oxygen and nutrient flux between sediment and water: a laboratory continuous-flow study. Mar. Ecol. Prog. Ser. 74, 263–279. doi: 10.3354/meps075263
Sundbäck, K., and Miles, A. (2000). Balance between denitrification and microalgal incorporation of nitrogen in microtidal sediments, NE Kattegat. Aquat. Microb. Ecol. 22, 291–300. doi: 10.3354/ame022291
Tamimi, A., Rinker, E. B., and Sandall, O. C. (1994). Diffusion coefficients for hydrogen sulfide, carbon dioxide, and nitrous oxide in water over the temperature range 293-368 K. J. Chem. Eng. Data 39, 330–332. doi: 10.1021/je00014a031
Weiss, J. N. (1997). The Hill equation revisited: uses and misuses. FASEB J. 11, 835–841. doi: 10.1096/fasebj.11.11.9285481
Werner, U., Billerbeck, M., Polerecky, L., Franke, U., Huettel, M., Van Beusekom, J. E. E., et al. (2006). Spatial and temporal patterns of mineralization rates and oxygen distribution in a permeable intertidal sand flat (Sylt, Germany). Limnol. Oceanogr. 51, 2549–2563. doi: 10.4319/lo.2006.51.6.2549
Wood, G. A. (1982). Data smoothing and differentiation procedures in biomechanics. Exerc. Sport Sci. Rev. 10, 308–362.
Keywords: denitrification, microphytobenthos, sand flat, nitrate uptake, permeable sediment
Citation: Schutte CA, Huanca-Valenzuela P, Lavik G, Marchant HK and de Beer D (2021) Advection Drives Nitrate Past the Microphytobenthos in Intertidal Sands, Fueling Deeper Denitrification. Front. Microbiol. 12:556268. doi: 10.3389/fmicb.2021.556268
Received: 27 April 2020; Accepted: 10 May 2021;
Published: 17 June 2021.
Edited by:
Peter Stief, University of Southern Denmark, DenmarkReviewed by:
Sokratis Papaspyrou, University of Cádiz, SpainMartina Herrmann, Friedrich Schiller University Jena, Germany
Copyright © 2021 Schutte, Huanca-Valenzuela, Lavik, Marchant and de Beer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Charles A. Schutte, schutte@rowan.edu
†Present address: Charles A. Schutte, Department of Environmental Science, Rowan University, Glassboro, NJ, United States Paulina Huanca-Valenzuela, Horn Point Laboratory, University of Maryland Center for Environmental Science, Cambridge, MD, United States