Burrowing benthic organisms promote water and solute fluxes across the sediment–water interface. Water and solutes penetrate the burrow walls and are transported into/out of the sediment when organisms flush their burrows with overlying water. Few studies have been done to investigate bioirrigation under shifting environmental conditions. We experimentally quantified bioirrigation by Chironomus plumosus larvae in the laboratory at 3 ranges of O2 saturation (low, medium, and high O2 concentrations), 2 temperatures (10 and 20°C), and over different seasons. We measured ventilation activities with O2 and flow-velocity microsensors, flow velocities during pumping periods with color tracers, pumping rates with conductivity exchange experiments, and rates of advective and diffusive water influx into the sediment by influx assays (NaCl was the tracer in both latter experiments). O2 saturations <12% extended pumping durations/h, whereas saturations <3% decreased pumping durations to ∼0. Flow velocities were 2× higher when O2 saturation was >50% than when it was <10%. Rising temperatures altered larval pumping (higher pumping frequency, lower pumping length) and increased flow velocity. Hence, pumping rate and rates of water influx were significantly higher at 20 than at 10°C. Seasonal variations in bioirrigation occurred despite constant laboratory conditions, i.e., the rate of water influx was significantly higher in spring/summer than in autumn. Our study shows that temporally varying environmental conditions should be considered when evaluating bioirrigation-mediated benthic fluxes across the sediment–water interface.
Burrowing benthic organisms enlarge the sediment surface area and promote active water transport across the sediment–water interface. Their ventilation activity, i.e., flushing their burrows with oxic water from the overlying water body (bioirrigation; Polerecky et al. 2006, Gallon et al. 2008) transports solutes and particulate constituents into the sediment pore water or back into the overlying water (Aller and Aller 1992, Aller 2001, Koretsky et al. 2002). Bioirrigation oxidizes the burrow walls, changes balances of redox reactions, and affects nutrient release or immobilization (Aller and Aller 1998, Lewandowski and Hupfer 2005, Meysman et al. 2006a). These changes influence microbial community structure in the sediment (Johnson et al. 1989, Stief et al. 2004, Bertics and Ziebis 2009) and chemical diagenesis (Andersson et al. 1988, Kristensen 2000, Stief and de Beer 2006). Thus, quantification of the benthic fluxes across the sediment–water interface induced by bioirrigating organisms is essential to understanding the microbiological and biogeochemical processes occurring in aquatic sediments (Boudreau and Marinelli 1994, Wang and van Cappellen 1996, Meysman et al. 2006b).
Quantification of bioirrigation-mediated benthic fluxes for a whole lake is difficult. First, the volumes of water pumped by diverse species differ greatly. Pumping rates between 1 and 16,000 mL h−1 individual−1 have been reported (Kristensen 1983a, Leuchs 1986, Osovitz and Julian 2002). Population-wide pumping rates and, thus, quantification of bioirrigation-mediated fluxes depend on population densities. Moreover, fluxes are quite variable because pumping rates are susceptible to environmental biotic and abiotic conditions, such as food availability, presence of predators, and weather (Seymour 1972, Smethie et al. 2003, Hölker and Stief 2005). The response of benthic organisms to changes in temperature and O2 concentration and the consequences for bioirrigation-mediated benthic fluxes, especially in freshwater sediments, have been investigated in few studies (Walshe 1948, Riisgård et al. 1992, Stief and Schramm 2010). Quantification of modified pumping rates under different environmental conditions is important because water temperatures and O2 concentrations in freshwater lakes are expected to change in the future (Gerten and Adrian 2002, Wilhelm and Adrian 2008, Kirillin 2010).
Season also has a strong effect on benthic organisms and their activity. Dipteran larvae (Neumann 1965), algae, fish, and mammals (Lee and McClintock 1986, Grey et al. 2004, van Dijk et al. 2005) show seasonal and annual cycles of activity, emergence, growth, or fecundity even when they are kept at constant conditions in the laboratory. Results of a few studies suggest seasonal variations of benthic fluxes across the sediment–water interface (Hammond et al. 1985, Martin and Sayles 1987, Schlüter et al. 2000). However, ventilation activity and bioirrigation commonly are investigated under laboratory conditions without taking season into account. Consequently, alteration of benthic fluxes in response to the season is not completely understood and, especially for freshwater lakes, poorly investigated.
The aim of our study was to use innovative measurement techniques to evaluate ventilation activity (Roskosch et al. 2011) and bioirrigation-mediated benthic fluxes under experimental conditions mimicking shifting environmental conditions (Fig. 1). We conducted experiments with Chironomus plumosus (Diptera: Chironomidae) larvae, which build U-shaped tubes in muddy sediment (burrow length = up to 40 cm, depth = up to 20 cm; McLachlan and Cantrell 1976, Roskosch et al. 2010, 2011) and periodically pump water through the tubes by undulating body movements (Roskosch et al. 2010). Chironomidae occur in high densities (100–1000 larvae/m2) in the profundal area of freshwater lakes (McLachlan 1977, Andersen and Jensen 1991, Kajak 1997). We hypothesized that changing O2 concentrations and temperatures would alter ventilation activity (pumping frequency, length, and duration, and flow velocity) of C. plumosus larvae and result in significantly altered bioirrigation-mediated benthic fluxes (pumping rate and influx rate). We also hypothesized that seasonal variation caused by annual cycles would occur even if bioirrigation were measured under constant conditions in the laboratory.
Fig.1.
Schematic diagram of the measurement techniques used for the 4 variables (experiment 1: pumping frequency, length, and duration; experiment 2: flow velocity; experiment 3: pumping rate; and experiment 4a, b: influx rate), and environmental variables that were manipulated in each experiment.

Methods
Organisms, sediment, and water
We collected C. plumosus larvae from sediments at 6-m water depth at a fixed location (lat 52°44′N, long 13°65′E) from the shallow (mean water depth = 4.9 m), eutrophic, polymictic Lake Müggelsee in Berlin, Germany, approximately every 2 mo over 3 y (November 2007 to November 2010). We used 4th-instar (final-instar) larvae that were ∼2 cm long and constructed burrows ∼1.7 mm in diameter (Roskosch et al. 2011). We used sieved Lake Müggelsee sediment (1-mm mesh size) topped with Lake Müggelsee water for all experiments and for housing the larvae. Larvae began constructing burrows as soon as they were placed on sediments. Before the experiments, we kept larvae up to 1 wk in aerated tanks in darkness at 10°C without food.
Experiments
Experiment 1. Pumping frequency, length, and duration.—
We measured the effects of O2 concentrations and temperature on pumping frequency (number of periods of constant undulations/h), pumping length (time of constant undulations/period [s]), and pumping duration (total time spent ventilating/h [min/h]) with microsensors (1st hypothesis). We filled U-shaped Perspex tubes (20-cm length, 7-mm inner diameter) with sediment and fixed them in water-filled plastic tanks (261 ×160 ×172 mm, 10 tubes each) as described in Roskosch et al. (2010, 2011). We placed larvae at one end of a tube to guide their digging of a burrow along the tube to the other end. In contrast to other studies in which water-filled tubes were used (Walshe 1950, 1951, Kristensen 1983b, c), this arrangement allowed us to know burrow length while providing larvae with relatively natural conditions that did not alter their ventilation activity (Kristensen 1983a, Roskosch et al. 2011). We used 60 larvae (n = 60, 6 different tanks) per range of O2 saturation and temperature. We acclimated larvae to the experimental conditions and allowed them to construct burrows for 2 d at the experimental temperature with constant aeration in darkness. Evaporation of water from the tanks was low, but if necessary, we replaced evaporated water with fresh lake water before the experiments were started.
We investigated the effects of low (0–3%), medium (>3–12%), and high (>12–100%) O2 saturation on pumping frequency, pumping length, and pumping duration at 10 and 20°C. We chose these ranges of O2 saturations to ensure that within-O2-group responses were similar (see Statistics below). We gassed the tanks constantly with a mixture of N2/CO2 based on the partial pressure of the ambient air so that the pH was not affected by gassing. We covered the tanks with cling film to limit the redissolution of O2 across the air–water interface. The different O2 saturations were realized by altering the degree of gassing and covering the tanks. We checked O2 concentration of the water continuously in the middle of the tank with O2 microsensors (described below). Measurements were started 2 h after O2 saturations were stable.
We investigated the effects of temperature at 10 and 20°C. We gassed the overlying water constantly with air (>>50% O2 saturation) because we did not want limited O2 to alter ventilation activity. The experimental temperatures corresponded to differences between spring/autumn (10°C) and summer (20°C) above-ground temperatures in Lake Müggelsee (R. Adrian, Leibniz-Institute of Freshwater Ecology and Inland Fisheries, personal communication).
We measured pumping frequencies and lengths with Clark-type O2 microsensors (OX-100, tip diameter = 50 or 100 µm; Unisense, Aarhus, Denmark) in burrow inlets and flow velocity microsensors (FS-25, tip diameter = 20 µm; Unisense Aarhus, Denmark; Revsbech et al. 1998, Brand et al. 2007) in burrow outlets as in Roskosch et al. (2011). Sensor tips were 2 mm deep in the middle of a burrow opening, and the start of a pumping period was indicated by simultaneous increases in O2 concentrations and flow velocities caused by pumping of the larva. When ventilation activity stopped, O2 concentrations decreased gradually, whereas flow velocities decreased instantly. We have no indication that the presence of microsensors affected ventilation activity behavior. Chironomus plumosus is a filter-feeder that usually dwells deep in its burrow and is not a surface-deposit feeder (Walshe 1951, Armitage et al. 1995). However, in few isolated instances, larvae detected the microsensors and built a new burrow tube beside the sensor tip. In such instances, only unaffected measurement periods were evaluated.
We calculated the arithmetic mean pumping frequency (f, /h) as the number of pumping periods divided by the length of the whole measurement time (average measurement time = 6 h). We estimated the arithmetic mean pumping length/period (l, /s) from all pumping periods during the whole measurement time. The mean pumping duration (D, min/h) was
Experiment 2. Flow velocity.—
We used color tracers to measure the effects of O2 concentrations and temperatures on flow velocities (mm/s) within pumping periods (1st hypothesis). We used the same sediment-filled U-shaped Perspex tubes, tanks, and larval collection/acclimation procedure described for experiment 1. We conducted the experiment analogously to experiment 1 except that temperature measurements were done at slightly different O2 saturations (low = 3–12%, medium = >12–30%, high = >30–100%) because pumping length was extremely short when O2 saturation was <3%, and reliable estimates of flow velocity were only possible for O2 saturations >3%. Furthermore, the effect of O2 concentration was investigated at 20°C only, whereas measurements in the range of high O2 concentration (>>50%) were done at 10 and 20°C. We used 40 larvae (n = 40, 4 tanks) per O2 group and temperature. We chose treatments based on statistical criteria as described for experiment 1.
We measured flow velocity with color tracers as described by Roskosch et al. (2011). We placed a dissolved color tracer (food color; Wusitta, Sitzendorf, Germany) with a pipette directly above the burrow inlet. When the larva started pumping, the tracer was sucked into the burrow. We immediately added a tracer of another color above the inlet. The time t (s) required for the 2nd tracer to become visible at the burrow outlet was timed with a stopwatch.
Assuming that the length l of the burrow (mm) built inside the Perspex tube was identical to the length of the Perspex tube, the flow velocity v (mm/s) of the color tracer during regular ventilation activity was
Experiment 3. Pumping rate.—
We evaluated the effects of temperature (10 and 20°C) on pumping rates (mL/s) with conductivity exchange experiments (1st hypothesis). We fixed a sediment-filled Perspex tube (22-cm length, 7-mm inner diameter) in a water-filled plastic tank separated into 2 equal chambers (500 mL/chamber, 1 tube/tank, each end of the tube in a different chamber; see Roskosch et al. 2011 for details). We used 20 larvae and setups (n = 20, 20 tanks) per temperature, and allowed larvae to acclimate to the experimental conditions as described for experiments 1 and 2. We replaced evaporated water with fresh lake water before starting the experiments.
The measurements of the pumping rates were done as described in Roskosch et al. (2011). We used color tracers to distinguish the inlet and outlet of a burrow. Repeated measurements confirmed that larvae were always oriented in the same direction during pumping periods. We added 50 mg/L NaCl to the chamber that housed the burrow inlet and ensured uniform mixing of the water by constant aeration of both chambers. The limit of salinity tolerance for C. plumosus (NaCl = 10,000 mg/L) is much higher than the concentration used in experiment 3 (50 mg/L NaCl corresponding to an increase in electrical conductivity of ∼100 µS/cm). The natural salinity of Lake Müggelsee is ∼750 µS/cm, so it is very unlikely that the increased salinity affected ventilation activity (Lauer 1969). Equilibration of the tracer concentration in the inlet chamber was recorded continuously by conductivity meters (GMH3430; Greisinger Electronic, Regenstauf, Germany; accuracy = 0.5% but at least ±2 µS/cm) for 2 d. We connected the inlet and outlet chambers with a water-filled tube (diameter = 1 mm) to ensure constant water levels in both chambers.
We calculated the pumping rate Q (mL/h) based on an equation for the exchange between 2 compartments (Mutz et al. 2007). We estimated Q by minimizing the root mean squared error between the conductivity measured in the outlet chamber and the theoretical conductivity Lf0 (µS/cm) at time t (h) based on the estimated Q. In the mass-balance equation 3, LfEq (µS/cm) is the conductivity in each chamber at equilibrium, Lfstart (µS/cm) is the conductivity in each chamber at the start of the experiment, and V (mL) is the equal water volume in both chambers:
Experiment 4a. Influx rate.—
We evaluated the effect of temperature on the influx rates ( = rate of advective and diffusive water influx of overlying water into the sediment [L/h]) (1st hypothesis). We filled round plastic columns (diameter = 59 mm) with sediment and water and added different densities of larvae (1 larva/column = 360 larvae/m2, height (h)sediment = 10 cm, hwater = 5 cm; 10 larvae/column = 3600 larvae/m2, hsediment = 15.5 cm, hwater = 15.5 cm). We chose this setup to assure almost natural conditions. The larvae were able to build their burrows freely in the sediment. The number of larvae corresponded to in situ abundances that are within the range normally found in eutrophic lakes (Svensson and Leonardson 1996, Kajak 1997).
We ran these experiments in batches (4 columns/date and temperature) in summer and autumn 2008 to 2010 (between July and November) (1 larva, 10°C, n = 24 columns; 1 larva, 20°C, n = 10 columns; 10 larvae, 10°C, n = 34 columns; 10 larvae, 20°C, n = 26 columns). We ran all influx experiments immediately after larvae, sediment, and water were sampled from Lake Müggelsee. We incubated columns for 2 d in darkness at 10°C so that larvae could build burrows and acclimate to the experimental conditions. On day 3, we raised temperatures to 20°C in ½ of the columns and allowed larvae to acclimate to the experimental temperatures for 1 d. We started measurements on day 4. We were focusing on the effect of temperature, but this effect could have been confounded by the effect of change in temperature. We aerated the overlying water in all columns constantly. We replaced evaporated water with fresh lake water before the experiments were started.
To start, we removed ∼10 mL of the overlying water and added NaCl dissolved in 10 mL of water (60 mg NaCl in experiments with hwater = 5 cm; 250 mg NaCl in experiments with hwater = 15.5 cm) to the laboratory column. The solution was mixed into the water column by aeration. We used the conductivity meters described above to record the decrease in NaCl in the overlying water caused by uptake of dissolved NaCl by the sediment for 7 d. We fitted the decrease in Lf0 (µS/cm) and the influx rate Qir (L/h) with the mass-balance Eq. 3 described above (Mutz et al. 2007). However, in this experiment the volume of the overlying water V (mL) and the volume of the pore water VPW (mL) had to be taken into account separately:
Experiment 4b. Influx rate under seasonal variation.—
We measured advective and diffusive influx rates (L/h) with influx experiments conducted seasonally over a period of 3 y (2nd hypothesis). We used methods identical to those described for experiment 4a (10 larvae; 3600 larvae/m2). We did the experiments at 10 and 20°C (n = 4 columns, 40 larvae/temperature and date, done on the same day) approximately every 2nd mo from November 2007 to November 2010 (18 dates for 10°C, 17 dates for 20°C). Experiment 4a for 10 larvae was a subset of experiment 4b (July to November).
Temperature coefficient.—
We calculated the temperature coefficient (Q10) from values measured at T1 = 10°C and T2 = 20°C as
The rate (R) represents any measurement at both temperatures, T1 and T2. The Q10-value illustrates the temperature dependency of pumping frequency, pumping length, pumping duration, flow velocity, pumping rate, and water influx. Q10 > 1 indicates that the variable increases, whereas Q10 < 1 indicates that the variable decreases with rising temperature.
Statistics
We measured pumping frequency, pumping length, flow velocity, and pumping rate at least twice per larva and obtained measurements for every larva that built tubes and ventilated them. By using a threshold approach we established boundaries for low, medium, and high O2 saturation. Therefore, observations between >0 and 100% were grouped so that the group variance was minimized (groups with low standard deviations were chosen).
We did all statistical analyses with the statistical software SPSS (version 14; IBM Corporation, Montauk, New York). We used Shapiro–Wilk tests to check for normality of all variables within O2 groups. We compared pumping frequency, pumping length, flow velocity, and pumping rate among O2 groups with 1-way analysis of variance (ANOVA) if data were normally distributed or the nonparametric Kruskal–Wallis H-test if they were not. We used post hoc (Least Significant Difference or Dunnett-T3) or repeated 2-sample tests to assess differences in the contributions between each pair of O2 groups. To compare 2 groups of values (e.g., between temperatures), we used Student's t-tests if distributions were normal and Mann–Whitney U-tests if they were not normal. The test and the p-value are reported when results were statistically significant (p ≤ 0.05). If >2 measurement techniques were compared, we adjusted the p-value for significance with a Bonferroni correction.
Results
Effect of O2
Experiment 1.—
Ventilation activities (pumping frequency, length, and duration) differed significantly among O2 groups (0–3% [low], >3–12% [medium], >12–100% [high]) at both temperatures (10 and 20°C). At 10°C, pumping frequency, length, and duration were significantly lower at low than at medium or high O2 and did not differ significantly between medium and high O2 (Table 1, Fig. 2A–C). However, pumping length and duration peaked at medium O2, whereas pumping frequency did not differ appreciably between medium and high O2. At 20°C, pumping frequency was significantly lower at medium than at high O2 and did not differ between medium and low or between low and high O2 (Fig. 2A). Pumping length and duration were significantly longer at medium than at low or high O2 and significantly longer at high than at low O2 (Fig. 2B, C). At low O2, larvae conducted minor pumping (Fig. 3), but when O2 saturation was <<2%, larvae almost stopped pumping.
Fig.2.
Box-and-whisker plots for pumping frequency (A), pumping length (B), and pumping duration (C) of Chironomus plumosus larvae at 0 to 3% (low), >3 to 12% (medium), >12 to 100% (high) O2 saturations at 10 and 20°C. Dotted lines indicate means, solid lines indicate medians, box ends are quartiles, whiskers show 2 standard deviations, and * indicates outliers. p-values indicate probabilities associated with 1-way analysis of variance done at each temperature (10 and 20°C) with O2 group as the main effect (see Table 1 for post hoc comparisons of means and results of tests for differences between temperatures). Bars with the same letters are not significantly different.

Fig.3.
Normal pumping at high O2 saturation (>20%) and minor pumping at low O2 saturation (<2%) measured by an O2 microsensor in a burrow inlet (2 mm deep) in experiment 1. Data were taken from separate experimental trials. Time is coded as days, hours, and minutes.

Table 1.
p-values for pairwise comparisons of pumping frequency, pumping length, and pumping duration among low, medium, and high O2 saturation groups at 10 and 20°C (experiment 1) and for flow velocity among low, medium, and high O2 saturation groups at 20°C (experiment 2). Superscripts indicate statistical test used to compare groups (a indicates analysis of variance, b indicates t-test, c indicates Mann–Whitney U-test). Bonferroni-corrected p-value for significance in all tests comparing a pair of groups was p ≤ 0.0167.

Experiment 2.—
Flow velocities measured within pumping periods at 20°C differed significantly among O2 groups (3–12% [low], >12–30% [medium], >30–100% [high]) and decreased with dropping O2 concentrations (Table 1, Fig. 4).
Fig.4.
Box-and-whisker plot for flow velocity of Chironomus plumosus larvae at 3 to 12% (low), >12 to 30% (medium), and >30 to 100% (high) O2 saturations measured at 20°C in experiment 2. Dotted lines indicate means, solid lines indicate medians, box ends are quartiles, whiskers show 2 standard deviations, and * indicates outliers. Data were analyzed with a Kruskal–Wallis H-test (see Table 1 for p-values). Bars with the same lower case letter are not significantly different.

Effect of temperature
Experiment 1.—
Ventilation activity differed between temperatures. Pumping frequency was significantly higher at 20 than at 10°C (Fig. 5A), whereas pumping length was significantly lower at 20 than at 10°C (Fig. 5B). Pumping duration did not differ significantly between temperatures (Fig. 5C).
Fig.5.
Box-and-whisker plots for pumping frequency (A), pumping length (B), pumping duration (C), flow velocity (D), and pumping rate (E) of Chironomus plumosus larvae at 10 and 20°C under high O2 concentrations (>>50%) in experiments 1 to 3. Dotted lines indicate means, solid lines indicate medians, box ends are quartiles, whiskers show 2 standard deviations, and * indicates outliers. Data were analyzed with t-tests. Bars with the same lower case letter are not significantly different.

Experiment 3.—
Constant pumping duration increased flow velocity at 20°C resulting in a significantly higher pumping rate at 20 than at 10°C (Fig. 5E). Q10-values for all variables except pumping length and duration were >1 (Table 2).
Table 2.
The temperature coefficient (Q10) for pumping frequency, pumping length, pumping duration, flow velocity, pumping rate, and rate of water influx at 10 and 20°C.

Rate of water influx affected by temperature and seasonal variation
Experiment 4a.—
The rate of advective and diffusive influx of overlying water into the sediment was affected by temperature at both larval densities. Q10-values were 1.7 or 1.6, respectively, for experiments with 1 and 10 larvae (Table 2). We calculated per capita influx rate to provide a basis of comparison among experiments. Per capita influx rates were 1.7× (at 10°C) to 1.9× (at 20°C) higher at the lower than at the higher larval density (Fig. 6A, B). Originally, no comparison was planned between experiments, so experimental designs with 1 and 10 larvae differed slightly in terms of starting concentration of NaCl, sediment volume, and water volume, but sediment surface area was identical.
Fig.6.
Box-and-whisker plots for rate of advective and diffusive water influx into the sediment per individual caused by pumping of Chironomus plumosus larvae at the individual (1 larva = 360 larvae/m2; experiment 4a) (A) and population level (10 larvae = 3600 larvae/m2; experiment 4b) (B) at 10 and 20°C. Experiment 4a was done in summer and autumn 2008, whereas experiment 4b was done in summer and autumn 2008, 2009, and 2010 (July to November each year). Dotted lines indicate means, solid lines indicate medians, box ends are quartiles, whiskers show 2 standard deviations, and * indicates outliers. Data were analyzed with t-tests. Bars with the same lower case letter are not significantly different.

Table 3.
Mean (±1 SD; median in parentheses) spring/summer (April–July) and autumn (September–October) advective and diffusive water influx (mL/h) into the sediment measured for 10 larvae at 10 and 20°C.

Experiment 4b.—
Rates of water influx measured over 3 y under constant laboratory conditions were significantly higher in spring/summer (April–July) than in autumn (September–October) (Fig. 7, Table 3). Across the entire period, the Q10-value was 1.7.
Fig.7.
Mean (±1 SD) rate of advective and diffusive water influx into the sediment caused by the pumping of a Chironomus plumosus population (10 larvae = 3600 larvae/m2) at 10 and 20°C during the 3-y study (experiment 4b). Experiments in grey boxes were conducted in summer and autumn 2008, 2009, and 2010 and were used to compare the rate of water influx into the sediment caused by pumping of C. plumosus larvae at 10 and 20°C.

Discussion
Effects of O2 concentration, temperature, and season
Our results supported hypothesis 1. Larvae modified their ventilation activity, and thus, bioirrigation, in response to the O2 concentration in the overlying water column (experiments 1 and 2). Our results suggest that larvae are able to balance a decrease in O2 by increasing their pumping duration. At very low O2 saturation, normal ventilation activity is minimized (Fig. 2A–C). We assume that in this case, O2 consumption, which is necessary to perform pumping, exceeds O2 gained by pumping. Some authors have described a switch to anaerobic metabolism for chironomid larvae in such almost-anoxic conditions (Augenfeld and Neess 1961, Nagell and Landahl 1978). At extremely low O2 saturation, larvae pump slowly, rarely, and briefly, but regularly (Fig. 3). Without occasional ventilation activity, larvae presumably would be unable to assess the O2 concentration of the overlying water, and thus, would not register a renewal of O2, should it occur. In our experiments, ventilation activity stopped when the O2 saturation in the overlying water decreased to <<2%. In addition, the O2 concentrations in the burrow adjacent to the larvae are probably even lower because O2 saturation decreases further when the water from the surrounding anoxic sediments flows through the burrow. Walshe (1950) also reported that C. plumosus pumping length increased as O2 concentrations decreased and that pumping length decreased drastically when O2 was <5% (cf. <3% in our study). Chironomus thummi, another bioirrigating larva occurring in shallow lakes, increased pumping length and pumping duration as O2 concentrations decreased (Leuchs 1986). In contrast to our study, the pumping frequency did not decrease with dropping O2 concentrations in either of these studies. Both studies were done with larvae in water-filled tubes with a diameter wider than natural burrows. Such tubes are known to alter bioirrigation (Kristensen 1983a, Roskosch et al. 2010). In our study, tubes were filled with sediment. Roskosch et al. (2010, 2011) showed that results from burrows built in sediment-filled tubes were not significantly different results from burrows built directly in the sediment.
Our results confirmed that water temperature affects bioirrigation of C. plumosus significantly and temperature changes result in modified ventilation activity and pumping rates (Experiments 1, 2, and 3). Presumably this change in activity with temperature is explained by the fact that the larvae are poikilotherms (Platzer 1967, Pinder 1986). Change in pumping frequency and length with temperature is consistent with results of studies on bioirrigating worms (Seymour 1972, Kristensen 1983c) and larvae (Walshe 1950, Leuchs 1986). However, Leuchs (1986) reported that the increased pumping rate of C. thummi with temperature was the result of a longer pumping duration, and not of increased flow velocity in the burrows as occurred in our study. The discrepancy may result from differences in the feeding behavior of the 2 species. Chironomus plumosus is a filter-feeder that catches particles from the water flowing through its burrow in a conical net (Walshe 1950, Armitage et al. 1995), whereas C. thummi is a deposit-feeder. Roskosch et al. (2010) showed that filter feeding increases the pumping rate substantially because a larva pumps for respiration and to meet its food requirements. Therefore, depending on the available food quantity, filter-feeders must pump a considerable amount of water through their burrows (Walshe 1951, Nielsen et al. 1995, Riisgård and Larsen 2005). When larvae ventilate their burrows with a faster flow velocity, the net is more quickly packed with particles. Hence, we assume that rising temperatures might result in shorter pumping periods of C. plumosus.
The rate of advective and diffusive water influx into the sediment was significantly altered by rising temperatures (Experiments 4a, b). The Q10-values (1.6–1.7; Table 2) confirmed that larval activity was modified by temperature changes (Naylor 1963). Our results concur with those from a previous study indicating that bioirrigation affects the pore-water composition more strongly in warm summer (up to 22°C) than in cold winter (down to 5°C) (Martin and Sayles 1987). Seasonal variation of bioirrigation-mediated benthic fluxes were assumed in a few in situ studies (Hammond et al. 1985, Schlüter et al. 2000, Smethie et al. 2003).
Our measurements of influx rates support the hypothesis that seasonal variation (circumannual cycling) occurs even if bioirrigation is measured when laboratory conditions, such as O2 supply or temperature, are constant (experiment 4b). We kept experimental conditions constant during the 3-y study, but rates of water influx into the sediment (spring/summer vs autumn) clearly differed between seasons. Such seasonal variability under constant conditions might occur if larvae had endogenous clocks responsible for the variations occurring in the laboratory. Circumannual cycles have been reported as a source of variability in other laboratory experiments. For instance, in a study of fish (juvenile roach, Rutilus rutilus) van Dijk et al. (2005) found significant differences in physical variables between summer and winter although laboratory temperatures were constant. Endogenous clocks or circumannual cycles also are known from studies of other insects, algae, cyanobacteria, and Daphnia populations kept at constant temperature, food availability, and photoperiod in the laboratory (Neumann 1965, Kao et al. 2010; N. Bauer, C. Dziallas, D. Baganz, and G. Staaks, Leibniz Institute of Freshwater Ecology and Inland Fisheries, personal communication). However, studies of the effects of endogenous clocks on ventilation activity and bioirrigation-mediated benthic fluxes, especially in freshwater sediments, are scarce (Saunders et al. 2004). Pumping rates are commonly determined under laboratory conditions without taking the season into account. Our results indicate that seasonal variations of bioirrigation should be considered carefully in future experiments. The major part of our experiments was conducted in summer to autumn 2008 because it was technically impossible to conduct all experiments at the same time. Thus, a part of the variability observed in our experiments might be attributable to seasonal variations.
Advanced environmental implications
A quantitative evaluation of bioirrigation-mediated benthic fluxes and their alterations is essential to understanding the biogeochemical processes in aquatic sediments and for applying dynamic diagenetic models (Boudreau and Marinelli 1994, Wang and van Cappellen 1996). Predictive modelling is becoming increasingly important because environmental conditions, such as water temperature and O2 concentrations, are shifting rapidly with global climate change. For example, long-term measurements of the water temperature of Lake Müggelsee show that its mean annual temperature has been increasing by ∼0.3°C per decade (Kirillin 2010). Moreover, Kirillin (2010) predicted that mixing patterns in the shallow lake eventually will change from polymictic (mixed water column) to dimictic (stratification periods in summer and winter) and, finally, to warm monomictic (no winter stratification, strong stratification in summer). Consequently, the hypolimnetic O2 concentration will decrease because warmer temperatures and earlier diatom spring blooms will lead to a higher mineralization rate, especially at the lake bottom (Wilhelm and Adrian 2008). In Lake Müggelsee, O2 saturation in the hypolimnion (measured directly above the sediment, 5-m water depth) usually does not decrease below 25% (R. Adrian, personal communication). However, during short-term stratification periods in summer, low O2 saturation (<2 mg/L) may occur (Wilhelm and Adrian 2008). In the future, stratification periods are expected to become longer (Kirillin 2010). Our results indicate that under such low O2 conditions, ventilation activity will be affected and pumping rate of C. plumosus larvae in Lake Müggelsee will decrease significantly.
Pumping rates and, thus, bioirrigation-mediated benthic fluxes are affected by population density. Mean larval density over 3 y in Lake Müggelsee was 514 ± 299 (n = 20) 4th-instar C. plumosus larvae/m2 sediment (Roskosch 2011). Assuming this abundance remains constant all year and that C. plumosus is the only bioirrigating species in the lake, a change in temperature from 10 to 20°C would yield a 40% increase (Q10 = 1.6–1.7; Table 2) in the population-wide pumping rate (from 430 to 600 L m−2 d−1, based on rates of 34.7 and 48.3 mL h−1 individual−1, respectively). Given the volume of Lake Müggelsee (3.65 ×109 L), this rate implies that the equivalent of the volume of the entire lake is pumped through the sediment in 10.4 d at 20°C instead of 14.4 d at 10°C (Roskosch et al. 2010). This extrapolation is only an approximation of the in situ pumping rate, but the values confirm the strong effects of bioirrigation on the aquatic environment.
Water influx into the sediment per individual decreased with increasing abundance of larvae (experiment 4a with 360 and 3600 larvae/m2 sediment). The efficiency per larva decreased with increased abundance by a factor of 0.58 in our experiments. McLachlan (1977) reported that C. plumosus kept at very high densities build shorter tubes. Consequently, water influx into the sediment depends not only on the pumping rate, but also on variables, such as community composition, burrow shape (I-, U-, Y- or branched tubes), and burrow dimension (effective exchange surface) (Graneli 1979, Koretsky et al. 2002). Hence, the influx rate per individual might be lower at densities >3600 larvae/m2. Bioirrigation-mediated fluxes at the sediment–water interface are highly sensitive to environmental conditions, population density, tube dimensions, and sediment physics (permeability and porosity, muddy or sandy sediment) (Meysman et al. 2006a, b). Bioirrigation-mediated benthic fluxes should be quantified under different environmental conditions and for diverse species and sediments.
O2 concentration, temperature, and season alter ventilation activity and advective and diffusive water influx into the sediment. As O2 concentration dropped, pumping duration increased and subsequently decreased, whereas flow velocity constantly decreased. As temperature rose, pumping duration remained the same, whereas flow velocity increased. Thus, changing conditions led to changes in pumping rates. Moreover, the rate of water influx was always higher in spring/summer than in autumn. This result shows that the season should be considered even when test organisms are brought into the laboratory and bioirrigation measurements are done under laboratory conditions. For meaningful quantification of bioirrigation-mediated benthic fluxes, measurements must be made carefully and under precise experimental conditions, especially when the data are being interpreted or extrapolated to population-scale rates.
Acknowledgements
The Leibniz-Institute of Freshwater Ecology and Inland Fisheries (IGB) and the German Research Foundation (DFG, LE 1356/3-1) are acknowledged for financial support. The authors thank S. Jordan, R. Hölzel, and B. Schütze (IGB) for data sampling and technical support, K. Pohlmann (IGB) for statistical advice, and A. Arthur-Goettig (bioXpress) for language editing. We thank T. Mehner and the participants in the course on scientific writing at IGB for helpful discussion on an earlier draft of this article.