Next Article in Journal
Settling of Road-Deposited Sediment: Influence of Particle Density, Shape, Low Temperatures, and Deicing Salt
Next Article in Special Issue
Benthic Metabolism in Fluvial Sediments with Larvae of Lampetra sp.
Previous Article in Journal
Remote Estimation of Water Quality Parameters of Medium- and Small-Sized Inland Rivers Using Sentinel-2 Imagery
Previous Article in Special Issue
Monthly Abundance Patterns and the Potential Role of Waterbirds as Phosphorus Sources to a Hypertrophic Baltic Lagoon
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Water
Volume 12
Issue 11
10.3390/w12113125
water-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Contrasting Effects of Bioturbation Studied in Intact and Reconstructed Estuarine Sediments

by
Marco Bartoli
1,2,*,
Sara Benelli
1,
Monia Magri
1,2,
Cristina Ribaudo
3,
Paula Carpintero Moraes
4 and
Giuseppe Castaldelli
5
1
Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, 4314 Parma, Italy
2
Marine Research Institute, University of Klaipeda, 92294 Klaipeda, Lithuania
3
Bordeaux INP, University Bordeaux Montaigne, 33600 Bordeaux, France
4
Laboratory of Biogeochemistry of Nutrients, Micronutrients and Traces in Oceans—Oceanographic Institute, São Paulo University, São Paulo 05508-120, Brazil
5
Department of Life Sciences and Biotechnology, University of Ferrara, 44121 Ferrara, Italy
*
Author to whom correspondence should be addressed.
Water 2020, 12(11), 3125; https://doi.org/10.3390/w12113125
Submission received: 2 October 2020 / Revised: 1 November 2020 / Accepted: 3 November 2020 / Published: 7 November 2020
(This article belongs to the Special Issue The Role of Macrobiota in Aquatic Nutrient Cycling)
Browse Figures
Versions Notes

Abstract

:
Macrofauna can produce contrasting biogeochemical effects in intact and reconstructed sediments. We measured benthic fluxes of oxygen, inorganic carbon, and nitrogen and denitrification rates in intact sediments dominated by a filter and a deposit feeder and in reconstructed sediments added with increasing densities of the same organisms. Measurements in reconstructed sediments were carried out 5 days after macrofauna addition. The degree of stimulation of the measured fluxes in the intact and reconstructed sediments was then compared. Results confirmed that high densities of bioturbating macrofauna produce profound effects on sediment biogeochemistry, enhancing benthic respiration and ammonium recycling by up to a factor of ~3 and ~9, respectively, as compared to control sediments. The deposit feeder also increased total denitrification by a factor of ~2, whereas the filter feeder activity did not stimulate nitrogen removal. Moreover, the effects of deposit feeders on benthic fluxes were significantly higher (e.g., on respiration and ammonium recycling) or different (e.g., on denitrification) when measured in intact and reconstructed sediments. In intact sediments, deposit feeders enhanced the denitrification coupled to nitrification and had no effects on the denitrification of water column nitrate, whereas in reconstructed sediments, the opposite was true. This may reflect active burrowing in reconstructed sediments and the long time needed for slow growing nitrifiers to develop within burrows. Results suggest that, in bioturbation studies, oversimplified experimental approaches and insufficient preincubation time might lead to wrong interpretation of the role of macrofauna in sediment biogeochemistry, far from that occurring in nature.

1. Introduction

Bioturbation by macrofauna is a fascinating topic that has fueled a diversified body of scientific research [1,2,3,4,5,6]. Macrofauna couple the pelagic and benthic environments in various ways. For instance, feces and pseudofeces deposition by filter feeders, containing viable and active phytoplankton, displaces pelagic primary production at the sediment level [7,8]. Filter feeders activity, such as that of clams, may process the entire water column of shallow water bodies many times per day, assuring water transparency, increased light penetration, and benthic primary production [9,10,11]. Filter feeders also excrete large amounts of the mineralized phytoplankton biomass back into the water column [12,13,14]. Biodeposits can increase the organic content of sediments and stimulate both aerobic and anaerobic metabolic pathways [15,16]. Filtration and excretion thus have the potential of maintaining pelagic and benthic primary production elevated. Deposit feeders, such as polychaete worms, explore surficial or deep sediment horizons to extract food and avoid predation. As oxygen (O2) penetration in sediments is extremely limited, burrowing macrofauna have evolved different ways to cope with anoxia and potentially toxic concentrations of compounds such as ammonia or sulfide. Burrow architectures may vary among species, but most deposit feeders ventilate their burrows by injecting and circulating oxic bottom water or by extracting pore water from sediments [17,18,19]. Such a need to get rid of metabolic end-products and to import O2 for their own respiration produces different biogeochemical side-effects on sediments [20]. For example, the volume of oxidized sediments in bioturbated areas is much higher than that in their not bioturbated counterparts, which means that the oxidation status of different elements is different depending on the presence of macrofauna. Ecosystem-level implications are multiple as oxidized sediments may retain elements avoiding their release to the water column or favor their permanent loss via coupled oxic–anoxic processes [21,22].
Macrofauna, through bioturbation but also feeding strategy and ecological interactions, affect both the physical and biological environments in sediments and in the water column. During their life span, macrofauna contribute to the creation and become part of a complex network of micro- and macro-organisms and of chemical and biological processes [23,24,25]. With this respect, little is known about the complex and intimate mutual interactions within holobionts, which are associations between microbes and macrofauna, and that affect benthic functioning [6,26].
Research activities in the field of bioturbation, with the aim of understanding density-dependent effects of macrofauna on selected processes, have largely adopted experimental approaches based on the addition of a single species to reconstructed sediments (e.g., sieved or freeze-thaw, in order to remove other macrofauna species, then homogenized and packed back into microcosms) [27,28,29,30]. Results from this kind of studies have allowed us to understand how important macrofauna is for a number of biogeochemical processes, including aerobic and anaerobic microbial respiration rates, nutrient recycling, reworking, and mixing of new and old organic matter and oxidation of chemically reduced metal pools [31,32,33]. The level of enhancement or inhibition of specific processes operated by macrofauna and measured with such approaches has been likely used to calibrate models, to upscale processes under varying macrofauna density or to propagate the effects to the water column [34]. However, such an approach may lead to a partial understanding of the true effect of macrofauna due to multiple factors. One of them is the absence of other macrofauna groups that may compete with, facilitate, or inhibit the targeted species. Another one is the time taken for an organism living in sediment, which is necessary to engineer the environment and, therefore, to change the chemical conditions and the biological communities (e.g., meiofauna and microbes). Due to these factors, we hypothesized that processes measured under in situ conditions may deeply differ from those measured shortly after sediment manipulation and macrofauna addition. We also hypothesized that such differences may vary along different macrofauna functional groups.
To test these hypotheses, we measured aerobic, anaerobic (denitrification), total respiration, and nutrient fluxes in intact and reconstructed sediments collected from two estuarine areas dominated by a filter feeder and a deposit feeder. Specific objectives were i) to analyze the effects of bioturbation on different benthic processes in intact and reconstructed sediments and ii) to evaluate whether differences between the two experimental conditions (i.e., intact versus reconstructed) depend upon macrofauna functional groups.

2. Materials and Methods

2.1. Study Sites

Samples were collected from the Sacca di Goro Lagoon (Figure 1), a shallow microtidal lagoon within the Po River Delta, Northern Italy (44.78–44.83° N, 2.25–12.33° E) [35]. The lagoon is approximately triangular in shape, has a surface area of 27 km2, an average depth of 1.5 m, and it is connected to the sea by a large mouth. The lagoon presents a western portion influenced by the nutrient-rich freshwater inlet of the Po di Volano, a central portion with marine influence and an eastern, confined area cultivated with clams (Ruditapes philippinarum). The average annual temperature and precipitation in this area are 17 ± 1 °C and 631 ± 117 mm, respectively. For the aims of this study, two stations were selected as representative of a pristine and of a heavily impacted lagoon area and located in the western and in the eastern portion of the lagoon, respectively. Station 1 has turbid water due to the Po di Volano freshwater inputs and has muddy, highly bioturbated sediments. During summer, the macrofauna community is dominated by the polychaete Alitta succinea. Station 2 has transparent water and sandy sediments and is located within a licensed area for clam farming; the dominant macrofauna is represented by the clam R. philippinarum.

2.2. Sampling Activities

2.2.1. Reconstructed Sediment Cores

Sediments and water were collected in July 2016. On 2 July 2016, nearly 50 L of sediment and 100 L of water were collected from the two stations. As both stations were shallow (<1 m depth), the water was collected by hand in 25 L tanks, whereas the sediment was collected with a shovel and transferred into 20 L buckets. Water and sediments were brought to the laboratories at Ferrara University within 2 h. Here, sediments from the two stations were (i) sieved (0.5 mm) in two tanks containing in situ water to remove large macrofauna, (ii) homogenized, and (iii) left to settle for 3 days. The tanks with sieved sediments and in situ water were maintained in the dark at in situ temperature, and the overlying water was kept oxic by air bubbling. Individuals of A. succinea and of R. philippinarum retrieved during the sieving were maintained alive in separate small aquaria with a few cm of sediments and 10 cm of well aerated water phase. On 5 July 2016, 9 cores (inner diameter = 8 cm, height = 30 cm) were subsampled from each tank. Water and sediment phases were subsampled by pushing vertically transparent Plexiglas cores in the tanks. Once retrieved, all cores were leveled in order to have a sediment height of 10 cm and a water column of 17 cm. Cores with muddy sediments from Station 1 were treated as follows: 3 cores were used as control, 3 cores were added with 3 individuals of A. succinea, and 3 cores were added with 6 individuals of A. succinea, to reproduce densities (600 and 1200 ind m−2, respectively) comparable to those measured in situ. In a few minutes, all added polychaetes burrowed into the sediments. Cores with sandy sediments from Station 2 were treated as follows: 3 cores were used as control, 3 cores were added with 2 clams (R. philippinarum), and 3 cores were added with 4 clams, to simulate densities at farmed areas (400 and 800 ind m−2, respectively). Within 30 min from clams addition all organisms burrowed within the sand, and only siphons were visible. All 9 + 9 cores were provided with stirring units (a rotating magnet suspended 5 cm above the sediment–water interface driven by an external motor at 40 rpm) and submersed in separate incubation tanks containing aerated and well mixed water from the two stations inside a temperature-controlled room maintained in the dark at the in situ temperature (25 °C).

2.2.2. Intact Sediments

On 9 July 2016, a second sediment and water sampling was carried out during which at each station 6 Plexiglas cores (inner diameter = 8 cm, height = 30 cm) were randomly collected by hand for flux and denitrification measurements and 5 Plexiglas cores (inner diameter = 4 cm, height = 20 cm) were collected for sediment characterization (for the methods see [14,36,37]). Cores with disturbed vertical profiles or with visually broken clam shells along walls were immediately discarded. Besides intact sediments, nearly 100 L of water from each site was collected for preincubation and incubation periods. Water column temperature, pH, and salinity were measured at the two sites by means of a multiple probe (mod 556, YSI Instrument, Yellow Spring, OH, USA). All the material was brought to the laboratories at Ferrara University within 2 h. Here, intact sediment cores from the two sites were submersed into two 50 L tanks containing in situ aerated and well mixed water inside a temperature-controlled room maintained in the dark at the in situ temperature. Intact cores for flux measurements were provided with a rotating magnet, suspended 5 cm above the sediment–water interface driven by an external motor at 40 rpm.

2.3. Dark Incubations for Flux and Denitrification Measurements

On 10 July 2016, 5 days after macrofauna addition in reconstructed cores and after an overnight preincubation of the intact sediment cores, two sequential dark incubations were performed for the measurements of net solute fluxes across the sediment–water interface and of denitrification rates. During the first incubation, measurements of benthic fluxes of dissolved O2, dissolved inorganic carbon (DIC), molecular nitrogen (N2), and dissolved inorganic nitrogen (DIN = NH4+ + NO2− + NO3−) were performed with a start–end approach [38]. The incubation lasted between 2 and 3 h in order to keep O2 drops from the initial values between 25 and 40%. Incubation time was set according to regular check (0.5 h) of dissolved O2 in the water phase via a microelectrode (OX-50, Unisense A/S, Aarhus, Denmark) through the cores lid. At time, zero water samples (n = 4) were collected from each incubation tank. At time, final water samples were collected from each core water phase through a valve in the cores lid. Incubation started when gas-tight lids were positioned on the top of each core. Water samples were collected via 60 mL glass syringes; 20 mL was disposed in 12 mL exetainers (Exetainer®, Labco Limited, High Wycombe, UK) and added with Winkler reagents for O2 measurements, 20 mL was disposed in 12 mL exetainers for DIC determination, which was measured with 0.1 N HCl titration with an automatic titration unit (ABU91, Radiometer, Copenhagen, Denmark) and a combined pH sensor [39]. An aliquot of 20 mL was filtered (Whatman GF/C, Cytiva, Marlborough, MA, USA) and transferred to scintillation vials for nutrient analysis. Dissolved inorganic nitrogen was determined spectrophotometrically. Ammonium (NH4+) was determined using salicylate and hypochlorite in the presence of sodium nitroprussiate [40]. Nitrate (NO3−) was determined after reduction to nitrite (NO2−) in the presence of cadmium, and NO2− was determined using sulphanilamide and N-(1-naphthyl)ethylenediamine [41,42]. Fluxes were defined as the difference between final and initial concentrations of solutes in the water column multiplied by the volume of the core water phase and divided by the area of the sediment and the incubation time (as reported in [38]). The respiratory quotient RQ, which is the ratio of dissolved inorganic carbon to oxygen fluxes (absolute values), was calculated.
After the incubation for flux measurements, lids were removed, the water in each tank was renewed, and the cores were submersed. Cores were maintained open and stirred in in situ oxygenated and well mixed water in the dark and at the in situ temperature. Approximately 5 h later, denitrification measurements started. The water level in each tank was lowered just below the cores opening, in order to separate the cores from the tank water. Labelled NO3− from a stock solution (20 mM, sodium nitrate-15N, ≥98 atom%, Sigma-Aldrich, Merk KGaA, Darmstadt, Germany) was added to the water column of each core in order to have a final 15NO3− concentration of 25 µM (Isotope Pairing Technique -IPT- [43]). Before and after 5 min from the addition of the labelled NO3− a subsample (5 mL) of the water phase was collected, filtered, and dispensed into 10 mL plastic vials for NO3− analysis. Labelled NO3− was allowed to diffuse into the anoxic horizon and then a gas-tight lid was positioned on the top of each core and the incubation started. The incubation length of the denitrification experiment overlapped that of the flux measurements. At the end of the incubation, the cores were opened, slurred, and samples were collected for labelled N2 analyses via membrane inlet mass spectrometry (MIMS, Bay instrument, Easton, MD, USA) [44]. Water samples for labelled N2 analyses were poisoned with ZnCl2 to stop microbial activity. Rates of total denitrification, denitrification coupled to nitrification (Dn), and water column nitrate-driven denitrification (Dw) were calculated, as described in [43]. Denitrification efficiency (DE) was calculated as the ratio between denitrification and cumulative flux of N from sediment (N2 + DIN) and expressed as a percentage [45]. The slurry in each core was then sieved (0.5 mm) in order to collect all macrofauna; macrofauna were recognized, and dry weight was quantified after drying the soft tissue at 70 °C to a constant weight.

2.4. Statistical Analysis

Differences among net solute fluxes or denitrification rates measured in the intact cores of the two stations and in the reconstructed cores of C, L, and H treatments were tested via one-way analysis of variance (ANOVA) followed by the post hoc Holm Sidak test. The dependency of solute fluxes or denitrification rates on macrofauna biomass were tested using linear regression analysis. Normality of residuals was confirmed numerically with a Shapiro–Wilks test.
The comparison of regressions between the biomass and measured process rates or calculated respiratory quotients and denitrification efficiency in intact and reconstructed cores were performed via analysis of covariance (ANCOVA). The level for statistical significance was set at 0.05. Statistical analyses were performed with R software v. 3.5.1 (The R Foundation for Statistical Computing, Vienna, Austria) [46].

3. Results and Discussion

3.1. General Features of the Study Sites

Station 1 had soft, organic sediments, with lower density and higher porosity, organic matter, C and N contents, and chlorophyll a as compared to Station 2 (Table 1). The biomass of macrofauna was one order of magnitude higher at Station 2 mostly due to high size of the cultivated R. philippinarum (Table 1). At Station 1, the dominant organism retrieved from sediments was the deep burrower A. succinea, with a density of 870 ± 530 ind m−2. Besides A. succinea, the macrofauna community at Station 1 included other small polychaetes and occasional individuals of Monocorophium insidiosum, Chironomus salinarius, Hydrobia spp., and Capitella capitata, whose pooled biomass was <10% of the total. At Station 2, R. philippinarum dominated the macrofauna community, with a density of 630 ± 490 ind m−2. Besides clams, sediments also hosted a few individuals of M. insidiosum and Echinogammarus spp., whose pooled biomass was <1% of the total. Station 1 was influenced by the freshwater inputs of the Po di Volano and had lower salinity and higher concentrations of DIN as compared to Station 2, with NO3− as dominant chemical species (Table 1).

3.2. Benthic Respiration, DIN Fluxes, and Rates of Denitrification in Intact Sediment Cores

Respiration rates tended to be higher at Station 2, dominated by R. philippinarum, however differences were not significant (Figure 2). Similar net O2 uptake and DIC release at the two stations, despite large difference in organic matter content, were probably the result of much higher contribution of macrofauna metabolic activity in the sandy Station 2, sustained by clams. Both stations recycled large and not statistically different amounts of NH4+. At the muddy, organic-rich Station 1, such amounts are likely sustained by high microbial ammonification, whereas at Station 2, such amounts are likely excreted by clams [47]. Net inorganic carbon to NH4+ efflux ratios averaged 6.9 ± 0.7 and 13.9 ± 4.7 at Station 1 and 2, respectively, suggesting in general the labile nature of the available organic matter. The low ratio measured at Station 1, very close to that of phytoplankton, may be due to much higher rates of NH4+ oxidation within sediments via nitrification [22]. This is confirmed by higher net NO2− and NO3− fluxes, averaging 119.35 ± 5.35 and 122.27 ± 15.11 µmol N m−2h−1, respectively, and by denitrification rates of NO3− produced within sediments, all significantly higher at Station 1. Elevated nitrification rates at Station 1 are probably supported by the activity of burrowers, increasing O2 availability within NH4+-rich pore waters [27,48,49].
Rates of total denitrification were higher at Station 1 due to higher denitrification coupled to nitrification and to higher denitrification of water column NO3− (Table 1). Denitrification efficiency was higher by a factor of 3 at Station 1 (15%) than at Station 2, due to similar NH4+ regeneration but much higher denitrification. Despite such difference, at both sites the major dissolved inorganic N flux was towards the water column as recycled NH4+, whereas the fraction lost to the atmosphere was minor.

3.3. Benthic Respiration, DIN Fluxes and Rates of Denitrification in Reconstructed Cores

Increasing densities of A. succinea and R. philippinarum in reconstructed sediments significantly stimulated benthic respiration, including the consumption of O2 and the production of DIC (Figure 3). The degree of stimulation of O2 and DIC fluxes was different along with the density gradient, resulting in an increase in the respiratory quotient, approaching the unit in the high density treatments. Both macrofauna functional groups largely enhanced NH4+ regeneration, whereas their net effects on NO2− and NO3− fluxes were less marked likely due to simultaneous stimulation of multiple microbial processes (e.g., nitrification and denitrification, Figure 3) [50]. The isotope pairing technique revealed a significant effect of A. succinea on the removal of water column NO3− via denitrification (Dw), whereas the same was not apparent with R. philippinarum [28]. Both macrofauna groups did not stimulate denitrification coupled to nitrification (Dn) (Figure 3). In reconstructed sediments, denitrification efficiency tended to decrease along with increasing macrofauna densities (p = 0.088 and p = 0.152 for A. succinea and R. phlippinarum, respectively) confirming a major effect of both deposit-feeding and filter-feeding macrofauna biomass on N recycling than on N losses [5]. Denitrification efficiencies overlapped values determined in intact sediments and averaged 15% (range 11–18%) and 5% (range 2–8%) in the A. succinea and R. philippinarum treatments, respectively.
These results allow to perform simple linear regressions between measured process rates and macrofauna biomass to calculate the degree of enhancement of a certain process (Table 2). The biomass of added A. succinea (10.3 ± 3.7 and 21.4 ± 4.8 gdw m−2 for L and H, respectively) was an order of magnitude lower than that of added R. philippinarum (222.2 ± 37.5 and 410.7 ± 68.8 gdw m−2 for L and H, respectively), however the effects produced on benthic biogeochemistry was comparable, except for denitrification rates. Burrowing A. succinea contributed less to respiration and excretion than R. philippinarum, but produced a much larger effect on microbial activity, compensating the different biomass. This interpretation is supported by the higher degree of stimulation of Dw by A. succinea. Denitrification in fact is not performed by deposit-feeding macrofauna but by the microbes growing along burrows that take advantage from the macrofauna-mediated bioirrigation, transporting NO3−-rich water within sediments [8,19,31]. Burrowing and ventilation activities moreover mobilize pore water solutes as NH4+, which is released in large amounts in A. succinea L and H treatments. We speculate that fluxes measured in sediments with R. philippinarum are instead mostly sustained by clams excretion [51,52].

3.4. Macrofauna Stimulate Differentially Biogeochemical Processes in Intact and in Reconstructed Sediments

Table 2 reports the slopes of linear regressions between measured fluxes or calculated rates and biomass of A. succinea and R. philippinarum in intact and reconstructed sediment cores, whereas Table 3 reports the results of the ANCOVA where such regressions were compared. As a general outcome, the enhancement of specific processes by the two functional groups of macrofauna (e.g., oxygen uptake, NH4+ efflux and denitrification of water column NO3− –Dw–) suggested much larger effects of the deposit feeder as compared to the filter feeder (Table 2). For example, the increase in benthic oxygen consumption per gram of dry weight of the burrower varied between 81 and 156 µmol gdw−1h−1. This is likely due to the construction and irrigation of subsurface burrow structures, to the increased sediment–water exchange surface and to the stimulation of microbial activity within sediments [28].
Interestingly, the stimulatory effects of A. succinea on most measured benthic fluxes, with NO3− as the only exception, were generally different in intact and reconstructed sediments, whereas this was not true for R. philippinarum (Table 2 and Table 3). In reconstructed sediment cores A. succinea produced much higher enhancement of O2 respiration and of water column NO3− consumption via denitrification (by a factor of 2) and of NH4+ recycling (by a factor 5) as compared to intact sediments. In unmanipulated sediments, A. succinea stimulated denitrification of NO3− produced via nitrification (Dn), which is expected, whereas this was not apparent in reconstructed sediments. In the latter, A. succinea stimulated the denitrification of water column NO3−, likely due to recently constructed burrows with sharp redox gradients and limited numbers of nitrifying bacteria on their walls.
Besides differences in the degree of stimulation of certain processes, measurements in intact and reconstructed bioturbated sediments led sometimes to contrasting effects of macrofauna. For example, A. succinea in intact cores enhanced DIC retention in sediments, whereas the same organism in reconstructed cores stimulated DIC release to the water (Table 2 and Table 3). Opposite effects by burrowing fauna were apparent also for NO2− fluxes.
Interestingly, the effects of R. philippinarum on solute fluxes in intact and reconstructed cores were not statistically different (Table 2 and Table 3). We speculate that clams drive fluxes mostly through their metabolic activities but produce a limited effect, at least in the short-term, on sedimentary variables and on microbial communities within sediments.

4. Conclusions

Results of this study confirm the key role of bioturbating macrofauna in benthic biogeochemistry, with deposit feeders simulating the removal of nitrogen via denitrification and increasing the denitrification efficiency of the benthic system. Deposit feeders also increase the rates of oxygen consumption, partly due to chemical or biological oxidation of anaerobic metabolism end products. As such, they may buffer the toxicity of sulphides and act upon redox-sensitive processes such as P retention, ultimately producing positive feedbacks on benthic biodiversity and negative feedbacks on eutrophication. In the Sacca di Goro lagoon, clams are cultivated, and their density is orders of magnitude higher than under natural conditions. Clams are demonstrated to favor benthic-pelagic coupling, but with such densities, they may determine a too rapid (re)cycling of excreted nutrients associated to ingested phytoplankton. As clams produce a scarce stimulation of denitrification, they maintain low the denitrification efficiency, and they may promote (or sustain) eutrophic conditions.
The present study suggests that results of experiments with sediments and macrofauna manipulations need to be carefully considered. We demonstrated here that the degree of stimulation of commonly measured biogeochemical processes as aerobic, anaerobic (e.g., denitrification), and total respiration or nutrient fluxes by macrofauna can be substantially different when calculated in intact, not manipulated natural sediments compared to reconstructed microcosms. Reasons are probably multiple and include, first of all, the consequences of sediment sieving and packing, which deeply alter vertical gradients, organic matter distribution, the sediment redox, the pore water chemistry and likely affect important fractions of reactive components of the organic matter pools. Sieving also removes all the macrofauna in order to highlight the biogeochemical effects of a single species. This helps to have very clear, often linear, and highly significant results; on the other hand, the simplification of the community removes an unaccountable number of ecological interactions that shape the true functioning of the benthic ecosystems. Results obtained from manipulated sediments are real, as they are measured, but are likely never found in the field, due to oversimplification of the experimental approach, including the removal of chemical gradients and multiple ecological interactions. Aquatic ecologists often use manipulative approaches, where they remove species to analyze the effect of a few targeted organisms. Results from this work, and in particular, fluxes measured in manipulated sediments, suggest that clear and easy-to-interpret results are obtained from these manipulations. The intensity of process stimulation calculated in manipulative studies can then be coupled to data on macrofauna abundances in order to upscale the effects of bioturbation at the whole system level. In doing so, one can get results that are quite far from the reality. Something similar likely occurs when incubating a deep burrower or a clam in a small vial to calculate their respiration and excretion rates; they will try to dig through the glass walls during the incubation and their metabolic activity will be much higher as compared to when they are laying in their burrows, surrounded and protected by sediments. A major fraction of the effects of macrofauna on sediments are indirect, and associated to the way they affect microbial communities and metabolic activities and, with them, chemical microgradients within sediment, which are flux drivers. Measurements in manipulative experiments should also consider these aspects and consider or quantify how long it takes for macrofauna to produce a steady state in terms of microbial communities, activities, and chemical gradients in manipulated sediments. The latter require a long conditioning time before measurements, as the appropriate preincubation period is probably closer to weeks than to days.

Author Contributions

Conceptualization, M.B.; validation, M.B., S.B., M.M., and C.R.; investigation, M.B. and G.C.; resources, M.B. and G.C.; writing—original draft preparation, M.B., S.B., and M.M.; writing—review and editing, C.R., P.C.M., and G.C.; funding acquisition, G.C. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Fabio Vincenzi is kindly acknowledged for assistance during lab activities and for the analysis of nutrients.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aller, R.C. Bioturbation and remineralization of sedimentary organic matter: Effects of redox oscillation. Chem. Geol. 1994, 114, 331–345. [Google Scholar] [CrossRef]
  2. Welsh, D.T. It’s a dirty job but someone has to do it: The role of marine benthic macrofauna in organic matter turnover and nutrient recycling to the water column. Chem. Ecol. 2003, 19, 321–342. [Google Scholar] [CrossRef]
  3. Mermillod-Blondin, F.; Rosenberg, R.; Fluviaux, H.; Claude, U.; Lyon, B. Ecosystem engineering: The impact of bioturbation on biogeochemical processes in marine and freshwater benthic habitats. Aquat. Sci. 2006, 68, 434–442. [Google Scholar] [CrossRef]
  4. Laverock, B.; Gilbert, J.A.; Tait, K.; Osborn, A.M.; Widdicombe, S. Bioturbation: Impact on the marine nitrogen cycle. Biochem. Soc. Trans. 2011, 39, 315–320. [Google Scholar] [CrossRef] [Green Version]
  5. Stief, P. Stimulation of microbial nitrogen cycling in aquatic ecosystems by benthic macrofauna: Mechanisms and environmental implications. Biogeosciences 2013, 10, 7829–7846. [Google Scholar] [CrossRef] [Green Version]
  6. Zilius, M.; Bonaglia, S.; Broman, E.; Chiozzini, V.G.; Samuiloviene, A.; Nascimento, F.J.A.; Cardini, U.; Bartoli, M. N2 fixation dominates nitrogen cycling in a mangrove fiddler crab holobiont. Sci. Rep. 2020, 10, 13966. [Google Scholar] [CrossRef]
  7. Gergs, R.; Rinke, K.; Rothhaupt, K.O. Zebra mussels mediate benthic-pelagic coupling by biodeposition and changing detrital stoichiometry. Freshw. Biol. 2009, 54, 1379–1391. [Google Scholar] [CrossRef] [Green Version]
  8. Benelli, S.; Bartoli, M.; Ribaudo, C.; Fano, E.A. Contrasting effects of an alienworm on benthic N cycling in muddy and sandy sediments. Water 2019, 11, 465. [Google Scholar] [CrossRef] [Green Version]
  9. Swanberg, I.L. The influence of the filter-feeding bivalve Cerastoderma edule L. on microphytobenthos: A laboratory study. J. Exp. Mar. Biol. Ecol. 1991, 151, 93–111. [Google Scholar] [CrossRef]
  10. Vaughn, C.C.; Hakenkamp, C.C. The functional role of burrowing bivalves in freshwater ecosystems. Freshw. Biol. 2001, 1431–1446. [Google Scholar] [CrossRef] [Green Version]
  11. Sandwell, D.R.; Pilditch, C.A.; Lohrer, A.M. Density dependent effects of an infaunal suspension-feeding bivalve (Austrovenus stutchburyi) on sand flat nutrient fluxes and microphytobenthic productivity. J. Exp. Mar. Biol. Ecol. 2009, 373, 16–25. [Google Scholar] [CrossRef]
  12. Asmus, R.M.; Asmus, H. Mussel beds: Limiting or promoting phytoplankton? J. Exp. Mar. Biol. Ecol. 1991, 148, 215–232. [Google Scholar] [CrossRef]
  13. Nizzoli, D.; Welsh, D.T.; Viaroli, P. Seasonal nitrogen and phosphorus dynamics during benthic clam and suspended mussel cultivation. Mar. Pollut. Bull. 2011, 62, 1276–1287. [Google Scholar] [CrossRef] [PubMed]
  14. Benelli, S.; Bartoli, M.; Racchetti, E.; Carpintero, P.; Mindaugas, M.; Lubiene, I.; Anna, E. Rare but large bivalves alter benthic respiration and nutrient recycling in riverine sediments. Aquat. Ecol. 2017, 51, 1–16. [Google Scholar] [CrossRef]
  15. Atkinson, C.L.; First, M.R.; Covich, A.P.; Opsahl, S.P.; Golladay, S.W. Suspended material availability and filtration—Biodeposition processes performed by a native and invasive bivalve species in streams. Hydrobiologia 2011, 667, 191–204. [Google Scholar] [CrossRef]
  16. Murphy, A.E.; Anderson, I.C.; Smyth, A.R.; Song, B.; Luckenbach, M.W. Microbial nitrogen processing in hard clam (Mercenaria mercenaria) aquaculture sediments: The relative importance of denitrification and dissimilatory nitrate reduction to ammonium (DNRA). Limnol. Oceanogr. 2016, 1589–1604. [Google Scholar] [CrossRef] [Green Version]
  17. Kristensen, E.; Kostka, J.E. Macrofaunal burrows and irrigation in marine sediment: Microbiological and biogeochemical interactions. In Interac-Tions between Macro- and Microorganisms in Marine Sediments; Kristensen, E., Haese, R.R.K.J., Eds.; American Geophysical Union: Washington, DC, USA, 2005; pp. 125–158. [Google Scholar]
  18. Murphy, E.A.K.; Reidenbach, M.A. Oxygen transport in periodically ventilated polychaete burrows. Mar. Biol. 2016, 163, 1–14. [Google Scholar] [CrossRef]
  19. Benelli, S.; Bartoli, M.; Zilius, M.; Vybernaite-Lubiene, I.; Ruginis, T.; Petkuviene, J.; Fano, E.A. Microphytobenthos and chironomid larvae attenuate nutrient recycling in shallow-water sediments. Freshw. Biol. 2018, 63, 187–201. [Google Scholar] [CrossRef]
  20. Volkenborn, N.; Woodin, S.A.; Wethey, D.S.; Polerecky, L. Bioirrigation in Marine Sediments. In Reference Module in Earth Systems and Environmental Sciences; Elias, S.A., Ed.; Elsevier Inc.: Amsertam, The Netherlands, 2016; pp. 1–9. ISBN 9780124095489. [Google Scholar]
  21. Hölker, F.; Vanni, M.J.; Kuiper, J.J.; Christof, M.; Grossart, H.-P.; Stief, P.; Adrian, R.; Lorke, A.; Dellwig, O.; Brand, A.; et al. Tube-dwelling invertebrates: Tiny ecosystem engineers have large effects in lake ecosystems. Ecol. Monogr. 2015, 85, 333–351. [Google Scholar] [CrossRef] [Green Version]
  22. Moraes, P.C.; Zilius, M.; Benelli, S.; Bartoli, M. Nitrification and denitrification in estuarine sediments with tube-dwelling benthic animals. Hydrobiologia 2018, 819, 217–230. [Google Scholar] [CrossRef]
  23. Mermillod-blondin, F.; Lemoine, D.G. Ecosystem engineering by tubificid worms stimulates macrophyte growth in poorly oxygenated wetland sediments. Funct. Ecol. 2010, 24, 444–453. [Google Scholar] [CrossRef]
  24. Herren, C.M.; Webert, K.; Drake, M.D.; Vander Zanden, M.J.; Einarsson, A.; Ives, A.R.; Gratton, C. Positive feedback between chironomids and algae creates net mutualism between benthic primary consumers and producers. Ecology 2017, 98, 447–455. [Google Scholar] [CrossRef]
  25. Magri, M.; Benelli, S.; Bondavalli, C.; Bartoli, M. Benthic N pathways in illuminated and bioturbated sediments studied with network analysis. Limnol. Oceanogr. 2018, 63, S68–S84. [Google Scholar] [CrossRef]
  26. Van der Heide, T.; Govers, L.L.; De Fouw, J.; Olff, H.; Van Der Geest, M.; Van Katwijk, M.M.; Piersma, T.; Van De Koppel, J.; Silliman, B.R.; Smolders, A.J.P.; et al. A three-stage symbiosis forms the foundation of seagrass ecosystems. Science 2012, 336, 1432–1434. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Svensson, J.; Enrich-Prast, A.; Leonardson, L. Nitrification and dentrification in a eutrophic lake sediment bioturbated by oligochaetes. Aquat. Microb. Ecol. 2001, 23, 177–186. [Google Scholar] [CrossRef]
  28. Nizzoli, D.; Bartoli, M.; Cooper, M.; Welsh, D.T.; Underwood, G.J.C.; Viaroli, P. Implications for oxygen, nutrient fluxes and denitrification rates during the early stage of sediment colonisation by the polychaete Nereis spp. in four estuaries. Estuar. Coast. Shelf Sci. 2007, 75, 125–134. [Google Scholar] [CrossRef]
  29. De Backer, A.; Van Coillie, F.; Provoost, P.; Van Colen, C.; Vincx, M.; Degraer, S. Bioturbation effects of Corophium volutator: Importance of density and behavioural activity. Estuar. Coast. Shelf Sci. 2011, 91, 306–313. [Google Scholar] [CrossRef]
  30. Bonaglia, S.; Bartoli, M.; Gunnarsson, J.S.; Rahm, L.; Raymond, C.; Svensson, O.; Yekta, S.S.; Brüchert, V. Effect of reoxygenation and Marenzelleria spp. bioturbation on Baltic Sea sediment metabolism. Mar. Ecol. Prog. Ser. 2013, 482, 43–55. [Google Scholar] [CrossRef] [Green Version]
  31. Pelegrì, S.P.; Blackburn, T.H. Bioturbation effects of the amphipod Corophium volutator on microbial nitrogen transformations in marine sediments. Mar. Biol. 1994, 121, 253–258. [Google Scholar] [CrossRef]
  32. Stocum, E.T.; Plante, C.J. The effect of artificial defaunation on bacterial assemblages of intertidal sediments. J. Exp. Mar. Biol. Ecol. 2006, 337, 147–158. [Google Scholar] [CrossRef]
  33. Kauppi, L.; Bernard, G.; Bastrop, R.; Norkko, A.; Norkko, J. Increasing densities of an invasive polychaete enhance bioturbation with variable effects on solute fluxes. Sci. Rep. 2018, 65, 1–12. [Google Scholar] [CrossRef]
  34. Braeckman, U.; Provoost, P.; Gribsholt, B.; Van Gansbeke, D.; Middelburg, J.J.; Soetaert, K.; Vincx, M.; Vanaverbeke, J. Role of macrofauna functional traits and density in biogeochemical fluxes and bioturbation. Mar. Ecol. Prog. Ser. 2010, 399, 173–186. [Google Scholar] [CrossRef]
  35. Viaroli, P.; Giordani, G.; Bartoli, M.; Naldi, M.; Azzoni, R.; Nizzoli, D.; Ferrari, I.; Comenges, J.M.Z.; Bencivelli, S.; Castaldelli, G.; et al. The Sacca di Goro lagoon and an arm of the Po River. In Estuaries; Springer: Berlin/Heidelberg, Germany, 2006; Volume 5, pp. 197–232. [Google Scholar]
  36. Bartoli, M.; Benelli, S.; Lauro, M.; Magri, M.; Vybernaite-Lubiene, I.; Petkuviene, J. Variable oxygen levels lead to variable stoichiometry of benthic nutrient fluxes in a hypertrophic estuary. Estuaries Coasts 2020. [Google Scholar] [CrossRef]
  37. Magri, M.; Benelli, S.; Bonaglia, S.; Zilius, M.; Castaldelli, G.; Bartoli, M. The effects of hydrological extremes on denitrification, dissimilatory nitrate reduction to ammonium (DNRA) and mineralization in a coastal lagoon. Sci. Total Environ. 2020, 740, 140169. [Google Scholar] [CrossRef]
  38. Dalsgaard, T.; Nielsen, L.P.; Brotas, V.; Viaroli, P.; Underwood, G.J.C.; Nedwell, D.B.; Sundbäck, K.; Rysgaard, S.; Miles, A.; Bartoli, M.; et al. Protocol Handbook for NICE-Nitrogen Cycling in Estuaries: A Project under the EU Research Programme: Marine Science and Technology (MAST III); Ministry of Environment and Energy National Environmental Research Institute, Denmark© Department of Lake and Estuarine Ecology: Silkeborg, Denmark, 2000; ISBN 8777725352. [Google Scholar]
  39. Anderson, L.G.; Hall, P.O.J.; Iverfeldt, A.; Rutgers van der Loejf, M.M.; Sundby, B.; Westerlund, S.F.G. Benthic respiration measured by total carbonate production. Limnol. Oceanogr. 1986, 31, 319–329. [Google Scholar] [CrossRef] [Green Version]
  40. Bower, C.E.; Holm-Hansen, T. A salicylate-hypochlorite method for determining ammonia in seawater. Can. J. Fish. Aquat. Sci. 1980, 37, 794–798. [Google Scholar] [CrossRef]
  41. Golterman, H.L.; Clymo, L.S.; Ohnstand, M.A.M. Methods for physical and chemical analysis of freshwaters. In Handbook, Blackwell Science, 2nd ed.; I.B.P: Oxford, UK, 1978; Volume 8. [Google Scholar]
  42. APHA (American Public Health Association). Standard Methods for the Examination of Water and Wastewaters, 18th ed.; American Public Health Association: Washington, DC, USA, 1992. [Google Scholar]
  43. Nielsen, L.P. Denitrification in sediment determined from nitrogen isotope pairing. FEMS Microbiol. Lett. 1992, 86, 357–362. [Google Scholar] [CrossRef]
  44. Kana, T.M.; Darkangelo, C.; Hunt, M.D.; Oldham, J.B.; Bennett, G.E.; Cornwell, J.C. Membrane Inlet Mass Spectometer for rapid high-precision determination of N2, O2, and Ar in environment water samples. Anal. Chem. 1994, 66, 4166–4170. [Google Scholar] [CrossRef]
  45. Eyre, B.D.; Ferguson, A.J.P. Comparison of carbon production and decomposition, benthic nutrient fluxes and denitrification in seagrass, phytoplankton, benthic microalgae- and macroalgae- dominated warm-temperate Australian lagoons. Mar. Ecol. Prog. Ser. 2002, 229, 43–59. [Google Scholar] [CrossRef] [Green Version]
  46. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  47. Nizzoli, D.; Bartoli, M.; Viaroli, P. Oxygen and ammonium dynamics during a farming cycle of the bivalve Tapes philippinarum. Hydrobiologia 2007, 587, 25–36. [Google Scholar] [CrossRef]
  48. Pelegri, S.P.; Blackburn, T.H. Effect of bioturbation by Nereis sp., Mya Arenaria and Cerastoderma sp. on nitrification and denitrification in estuarine sediments. Ophelia 1995, 42, 289–299. [Google Scholar] [CrossRef]
  49. Bosch, J.A.; Cornwell, J.C.; Kemp, W.M. Short-term effects of nereid polychaete size and density on sediment inorganic nitrogen cycling under varying oxygen conditions. Mar. Ecol. Prog. Ser. 2015, 524, 155–169. [Google Scholar] [CrossRef] [Green Version]
  50. Bartoli, M.; Nizzoli, D.; Viaroli, P.; Turolla, E.; Castaldelli, G.; Fano, E.A.; Rossi, R. Impact of Tapes philippinarum farming on nutrient dynamics and benthic respiration in the Sacca di Goro. Hydrobiologia 2001, 455, 203–212. [Google Scholar] [CrossRef]
  51. Welsh, D.T.; Nizzoli, D.; Fano, E.A.; Viaroli, P. Direct contribution of clams (Ruditapes philippinarum) to benthic fluxes, nitrification, denitrification and nitrous oxide emission in a farmed sediment. Estuar. Coast. Shelf Sci. 2015, 154, 84–93. [Google Scholar] [CrossRef] [Green Version]
  52. Murphy, A.E.; Nizzoli, D.; Bartoli, M.; Smyth, A.R.; Castaldelli, G.; Anderson, I.C. Variation in benthic metabolism and nitrogen cycling across clam aquaculture sites. Mar. Pollut. Bull. 2018, 127, 524–535. [Google Scholar] [CrossRef]
Water 12 03125 g001 550
Figure 1. (A) shows the location of the two sampling stations in Sacca di Goro Lagoon (Po River Delta, Northern Italy). (B) shows reconstructed sediments from the two stations, added with no (C, control) or increasing densities of Alitta succinea (L, low, 3 ind core−1, and H, high, 6 ind core−1) or Ruditapes philippinarum (L, low, 2 ind core−1, and H, high, 4 ind core−1). (C) shows randomly collected intact sediments from the two stations, with natural macrofauna community, dominated by the species indicated. (D) shows a reconstructed sediment core with burrowing individuals of Ruditapes philippinarum.
Figure 1. (A) shows the location of the two sampling stations in Sacca di Goro Lagoon (Po River Delta, Northern Italy). (B) shows reconstructed sediments from the two stations, added with no (C, control) or increasing densities of Alitta succinea (L, low, 3 ind core−1, and H, high, 6 ind core−1) or Ruditapes philippinarum (L, low, 2 ind core−1, and H, high, 4 ind core−1). (C) shows randomly collected intact sediments from the two stations, with natural macrofauna community, dominated by the species indicated. (D) shows a reconstructed sediment core with burrowing individuals of Ruditapes philippinarum.
Water 12 03125 g001
Water 12 03125 g002 550
Figure 2. Dark fluxes of dissolved O2 and dissolved inorganic carbon (DIC), dissolved inorganic nitrogen (DIN), and rates of denitrification measured in intact sediment cores at the two stations. In the upper panel, the y-axes on the right reports the respiratory quotient (RQ), calculated as the ratio between DIC and O2 fluxes (absolute value). Where present, letters above bars indicate statistically significant differences. Averages ± standard errors (n = 6) are reported.
Figure 2. Dark fluxes of dissolved O2 and dissolved inorganic carbon (DIC), dissolved inorganic nitrogen (DIN), and rates of denitrification measured in intact sediment cores at the two stations. In the upper panel, the y-axes on the right reports the respiratory quotient (RQ), calculated as the ratio between DIC and O2 fluxes (absolute value). Where present, letters above bars indicate statistically significant differences. Averages ± standard errors (n = 6) are reported.
Water 12 03125 g002
Water 12 03125 g003 550
Figure 3. Dark fluxes of dissolved O2 and DIC, DIN, and rates of denitrification measured in reconstructed sediments added with Alitta succinea (Station 1) and Ruditapes philippinarum (Station 2). C indicates the control treatment, L the low density treatment (600 and 400 ind m−2 for polychaetes and clams, respectively), and H the high density treatment (1200 and 800 ind m−2 for polychaetes and clams, respectively). In the upper panels, the y-axes on the right reports the respiratory quotient (RQ), calculated as the ratio between DIC and O2 fluxes (absolute value). Where present, letters above bars indicate statistically significant differences. Averages ± standard errors (n = 3) are reported.
Figure 3. Dark fluxes of dissolved O2 and DIC, DIN, and rates of denitrification measured in reconstructed sediments added with Alitta succinea (Station 1) and Ruditapes philippinarum (Station 2). C indicates the control treatment, L the low density treatment (600 and 400 ind m−2 for polychaetes and clams, respectively), and H the high density treatment (1200 and 800 ind m−2 for polychaetes and clams, respectively). In the upper panels, the y-axes on the right reports the respiratory quotient (RQ), calculated as the ratio between DIC and O2 fluxes (absolute value). Where present, letters above bars indicate statistically significant differences. Averages ± standard errors (n = 3) are reported.
Water 12 03125 g003
Table
Table 1. Sediment descriptors, including the biomass of the dominant macrofauna (Alitta succinea and Ruditapes philippinarum at Station 1 and 2, respectively) and water physico-chemical characteristics.
Table 1. Sediment descriptors, including the biomass of the dominant macrofauna (Alitta succinea and Ruditapes philippinarum at Station 1 and 2, respectively) and water physico-chemical characteristics.
Sediment Features
DensityPorosityOrganic MatterCNChl aMacrofauna
(g mL−1) (%)(%)(%)(mg m−2)(gdw m−2)
Station 11.14 ± 0.110.74 ± 0.066.03 ± 0.071.79 ± 0.040.20 ± 0.0137.3 ± 9.816.0 ± 9.1
Station 21.74 ± 0.020.47 ± 0.011.61 ± 0.080.35 ± 0.020.04 ± 0.004.7 ± 0.8220.8 ± 292.1
Granulometry
<0.063 mm0.063 < x < 0.125 mm0.125 < x < 0.5 mm0.5 < x < 2 mm>2 mm
(%)(%)(%)(%)(%)
Station 124.422.050.43.10
Station 210.735.251.32.70
Water features
TemperatureSalinitypHNH4+NO2−NO3−
(°C) (µM)(µM)(µM)
Station 12557.771.86.631.6
Station 225177.9614.61.33.1
Table
Table 2. Slope, intercept, level of significance, and r2 values of linear regressions between solutes fluxes (µmol m−2h−1), respiratory quotient, denitrification rates (µmol m−2h−1), and efficiency and macrofauna biomass of Alitta succinea and Ruditapes philippinarum (gdw m−2) retrieved in intact and reconstructed cores from the two stations. For slopes and intercepts, standard errors are reported in brackets.
Table 2. Slope, intercept, level of significance, and r2 values of linear regressions between solutes fluxes (µmol m−2h−1), respiratory quotient, denitrification rates (µmol m−2h−1), and efficiency and macrofauna biomass of Alitta succinea and Ruditapes philippinarum (gdw m−2) retrieved in intact and reconstructed cores from the two stations. For slopes and intercepts, standard errors are reported in brackets.
Alitta succineaRuditapes philippinarum
Intact SedimentsReconstructed SedimentsIntact SedimentsReconstructed Sediments
SlopeInterceptpr2SlopeInterceptpr2SlopeInterceptpr2SlopeInterceptpr2
O2−81.75 (25.67)−4031.78 (461.74)0.0330.72−158.53 (13.25)−2283.70 (291.41)<0.0010.95−12.80 (2.94)−5417.04 (1016.76)0.0120.83−6.98 (2.27)−4149.12 (969.48)0.0180.57
DIC−137.11 (75.55)7417.94 (1359.12)0.1440.45264.35 (45.38)−144.03 (997.82)<0.0010.8312.57 (2.26)4150.56 (781.78)0.0050.8912.09 (1.75)112.07 (748.16)<0.0010.87
RQ−0.04 (0.02)1.71 (0.29)0.0550.640.04 (0.01)−0.09 (0.28)0.0150.59−0.0001 (0.0006)0.918 (0.207)0.8500.010.001 (0.0002)0.067 (0.082)<0.0010.88
NH4+5.90 (3.86)736.22 (69.37)0.2010.3731.43 (2.71)169.36 (59.60)<0.0010.950.87 (0.69)514.09 (239.52)0.2760.281.47 (0.40)145.15 (169.82)0.0080.66
NO2−1.13 (0.45)−101.37 (8.12)0.0670.61−2.66 (0.71)145.65 (15.53)0.0070.67−0.03 (0.02)30.42 (7.14)0.2860.27−0.15 (0.05)149.04 (19.52)0.0140.60
NO3−−2.46 (1.63)161.53 (29.26)0.2050.36−2.19 (4.68)279.11 (102.81)0.6540.03−0.05 (0.02)5.76(8.53)0.1030.53−0.09 (0.11)148.70 (46.32)0.4580.08
DIN4.57 (4.49)999.12 (80.79)0.3670.2126.59 (4.77)594.12 (104.93)<0.0010.820.79 (0.68)550.26 (234.96)0.3070.261.24 (0.40)442.88 (168.65)0.0170.58
Dw1.28 (0.93)93.19 (16.65)0.2390.322.48 (0.29)71.43 (6.30)<0.0010.910.004 (0.005)9.33(1.75)0.4320.160.001 (0.001)5.66 (0.59)0.4340.09
Dn2.61 (0.75)35.31 (13.51)0.0250.750.06 (0.30)34.99 (6.53)0.8520.01−0.007 (0.017)17.90 (5.91)0.7050.04−0.009 (0.007)19.98 (2.86)0.2130.21
Dtot3.89 (0.61)128.50 (10.88)0.0030.912.53 (0.33)106.42 (7.29)<0.0010.89−0.003 (0.016)27.22 (5.67)0.8850.01−0.008 (0.008)25.64 (3.27)0.3290.14
DE0.21 (0.08)11.73 (1.52)0.0680.61−0.11 (0.05)15.20 (1.09)0.0540.43−0.005 (0.006)6.21(2.08)0.4360.16−0.007 (0.003)6.29 (1.20)0.0420.47
Table
Table 3. Results of the ANCOVA between the slopes of fluxes or calculated respiratory quotients and denitrification efficiency and macrofauna biomass in intact and reconstructed sediments from the two stations. While for Ruditapes philippinarum, most comparisons are not significant, for Alitta succinea the two conditions lead to rather different slopes.
Table 3. Results of the ANCOVA between the slopes of fluxes or calculated respiratory quotients and denitrification efficiency and macrofauna biomass in intact and reconstructed sediments from the two stations. While for Ruditapes philippinarum, most comparisons are not significant, for Alitta succinea the two conditions lead to rather different slopes.
Alitta succineaRuditapes philippinarum
FpFp
O26.5590.0272.4860.143
DIC16.7230.0020.0290.869
RQ10.4390.0087.8090.017
NH4+20.4240.0010.6620.433
NO2−7.9090.0174.1150.067
NO3−0.0010.9760.0570.816
DIN5.5500.0380.3710.555
Dw2.1490.1710.5860.459
Dn11.7400.0060.1210.735
Dtot3.4090.0920.1180.737
DE9.0410.0120.0890.771
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bartoli, M.; Benelli, S.; Magri, M.; Ribaudo, C.; Moraes, P.C.; Castaldelli, G. Contrasting Effects of Bioturbation Studied in Intact and Reconstructed Estuarine Sediments. Water 2020, 12, 3125. https://doi.org/10.3390/w12113125

AMA Style

Bartoli M, Benelli S, Magri M, Ribaudo C, Moraes PC, Castaldelli G. Contrasting Effects of Bioturbation Studied in Intact and Reconstructed Estuarine Sediments. Water. 2020; 12(11):3125. https://doi.org/10.3390/w12113125

Chicago/Turabian Style

Bartoli, Marco, Sara Benelli, Monia Magri, Cristina Ribaudo, Paula Carpintero Moraes, and Giuseppe Castaldelli. 2020. "Contrasting Effects of Bioturbation Studied in Intact and Reconstructed Estuarine Sediments" Water 12, no. 11: 3125. https://doi.org/10.3390/w12113125

APA Style

Bartoli, M., Benelli, S., Magri, M., Ribaudo, C., Moraes, P. C., & Castaldelli, G. (2020). Contrasting Effects of Bioturbation Studied in Intact and Reconstructed Estuarine Sediments. Water, 12(11), 3125. https://doi.org/10.3390/w12113125

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop