The Impact of River Discharge and Water Temperature on Manganese Release from the Riverbed during Riverbank Filtration: A Case Study from Dresden, Germany

The climate-related variables, river discharge, and water temperature, are the main factors controlling the quality of the bank filtrate by affecting infiltration rates, travel times, and redox conditions. The impact of temperature and discharge on manganese release from a riverbed were assessed by water quality data from a monitoring transect at a riverbank filtration site in Dresden-Tolkewitz. Column experiments with riverbed material were used to assess the Mn release for four temperature and three discharge conditions, represented by varying infiltration rates. The observed Mn release was modeled as kinetic reactions via Monod-type rate formulations in PHREEQC. The temperature had a bigger impact than the infiltration rates on the Mn release. Infiltration rates of <0.3 m3/(m2·d) required temperatures >20 ◦C to trigger the Mn release. With increasing temperatures, the infiltration rates became less important. The modeled consumption rates of dissolved oxygen are in agreement with results from other bank filtration sites and are potentially suited for the further application of the given conditions. The determined Mn reduction rate constants were appropriate to simulate Mn release from the riverbed sediments but seemed not to be suited for simulations in which Mn reduction is likely to occur within the aquifer. Sequential extractions revealed a decrease of easily reducible Mn up to 25%, which was found to reflect the natural stratification within the riverbed, rather than a depletion of the Mn reservoir.


Introduction
Riverbank filtration (RBF) has been successfully used as a natural and cost-efficient water treatment method in many countries in Europe [1,2], the USA [3,4], Africa [5,6], and Asia [7,8].RBF can naturally occur or can be induced by pumping, whereby wells are placed adjacent to the river that creates a hydraulic potential gradient from the river towards the wells.RBF triggers a variety of natural attenuation processes that can largely improve the water quality of the bank filtrate (BF) and lower the post-treatment effort [9,10].For example, some organic micropollutants are effectively removed within the first meter of infiltration under oxic conditions [11].A series of redox processes along the flow path of the infiltrate can adversely affect the BF quality.Depletion of dissolved oxygen (DO), followed by denitrification and the reduction of manganese (Mn) minerals in the riverbed and the aquifer can cause elevated Mn concentrations, which require subsequent treatment [12][13][14].
However, the resulting Mn concentration depends on various factors, including residence time and water temperature [15].Assuming a steady well operation, the residence time can change due to river water level fluctuations, which are usually connected to river discharge variations (floods/droughts).For example, decreasing water levels usually lower the hydraulic gradient and prolong the distance between the river and the RBF wells.The lowered hydraulic gradient also affects the infiltration rate.Additionally, the infiltration rate is impacted by a usually smaller infiltration area at lower water levels.Water temperature variations affect the viscosity of the water and, therefore, also affect the residence times and infiltration rates.Additionally, according to the rule of Van 't Hoff [16], increasing temperatures lead to increased biological degradation rates, whereby Mn release is considered to be largely biologically mediated.
Current climate forecasts are thought-provoking at many RBF sites [17][18][19].Droughts potentially lower the river discharge, extend travel times, and promote anaerobic conditions along the flow path, while floods can shorten travel times or cause, for example, breakthroughs of pathogens and organic micropollutants.The current increase of climate extremes in Europe is expected to continue, with a higher frequency of heat waves, long-lasting droughts in some regions, heavy precipitation events and river floods [20].For the German state of Saxony, the air temperature is projected to increase by 3-3.5 • C and the mean summer rainfall is to decrease by 20-25% by the year 2100 [21].
The Waterworks (WW) Dresden-Tolkewitz (Saxony, Germany) was built in 1898 and is one of the oldest RBF schemes in Europe.This study investigated the impact of the climate-related variables, temperature and discharge, on the BF quality.Water quality data of a 10-year time span from a monitoring transect were examined to identify discharge-and temperature-related patterns.To assess the potential Mn release from the riverbed in Dresden-Tolkewitz, three columns filled with riverbed sediment from the Elbe river emulated three infiltration rates and four temperature conditions.In order to use the results from the column experiment for a planned modeling of the transect, the observed Mn release was reproduced by hydrogeochemical modeling with PHREEQC.Additionally, a sequential extraction procedure was applied to the riverbed sediment from the columns after the experiment.Based on that data, the implications of the results to the redox-related BF quality in a potentially changing climate are discussed with a focus on Mn.

Description of the RBF Waterworks Dresden-Tolkewitz
The WW Dresden-Tolkewitz is located at the upper Elbe river in Germany (Figure 1).Three siphon well galleries with 72 vertical wells abstract up to 1500 m 3 /h (36,000 m 3 /d).The portion of riverbank filtrate is around 83% during mean flow and 70-75% during low flow conditions [22].The focus of this study was a 95 m wide transect between the Elbe river and a production well (PW), which already was the focus of previous riverbed clogging studies [22].The PW fully penetrates the aquifer and the 4 m long filter screen is located directly above the aquitard.The transect has three observation wells (OW 1, 2 and 3).Each OW has sampling points at three depths (upper, middle, lower = OW i-1, i-2, i-3).During the mean flow, the nearest OW (OW 1) is around 21-30 m apart from the riverbank.OW 2 and OW 3 are around 40 and 80 m apart from the bank during mean flow.The average travel time along the transect is between 24 and 30 days [23].
The mean discharge of the Elbe river in Dresden is 332 m 3 /s (at 1.84 m river stage).The discharge varies during the mean low and high flow periods between 110 m 3 /s (0.75 m) and 1700 m 3 /s at a water level of 5.47 m [22].The climate in Dresden is humid continental with warm summers.The alluvial aquifer is unconfined, composed of gravel and coarse sand with a saturated thickness of 11-14 m and has a hydraulic conductivity of 1-2 × 10 −3 m/s [22].

Regular Monitoring in Dresden-Tolkewitz
In this study, the evaluation period of the water quality data from the WW Dresden-Tolkewitz was a 10-year time span from 1 January 2006 to 31 December 2016.Regular samples were taken twice a year from all sampling points of each OW by the waterworks staff.Additional event-based samples were taken during low flow periods and, if possible, during high flow periods.Sampling was carried out according to DVGW [24] and corresponding to earlier guidelines.Water quality data for the Elbe river were taken from the database of the Saxon state ministry [25].The relevant sampling point is located at river km 43.5, around 3 km upstream of the WW.Water levels in this study refer to the federal water level "Dresden Augustusbruecke" at river km 55.63 [26].

Set-Up of the Column Experiments
To understand the behavior of Mn at the RBF site in Dresden-Tolkewitz, three columns with riverbed sediments were set up in the laboratory of the University of Applied Sciences Dresden (Figure S1).Each column was 1 m long, had an inner diameter of 0.08 m, and was made up of galvanized steel.The filling material was riverbed sediment from the Elbe River.The riverbed material was collected in front of the investigated monitoring cross-section of the RBF Waterworks Dresden-Tolkewitz (Figure 1).The riverbed sediment was recovered around 20 m apart from the riverbank during a low discharge period.Due to the very low gradient of the riverbed towards the riverbank, the area around the sampling point was already flooded at slightly higher water levels, which still would occur during mean low flow conditions.Hence, the area around the sampling point can be considered as a potential infiltration area.Because of the relatively coarse riverbed, undisturbed sampling was not possible.Thus, the upper 5 cm of the riverbed where scratched first to represent the clogging layer.Subsequently, the deeper riverbed material was dug out layer-wise and sieved in place to a grain size <4 mm.Immediately after transporting it to the lab, the wet riverbed material was filled into columns in ≈0.1 m thick, separately compacted layers with the clogging layer on top.

Regular Monitoring in Dresden-Tolkewitz
In this study, the evaluation period of the water quality data from the WW Dresden-Tolkewitz was a 10-year time span from 1 January 2006 to 31 December 2016.Regular samples were taken twice a year from all sampling points of each OW by the waterworks staff.Additional event-based samples were taken during low flow periods and, if possible, during high flow periods.Sampling was carried out according to DVGW [24] and corresponding to earlier guidelines.Water quality data for the Elbe river were taken from the database of the Saxon state ministry [25].The relevant sampling point is located at river km 43.5, around 3 km upstream of the WW.Water levels in this study refer to the federal water level "Dresden Augustusbruecke" at river km 55.63 [26].

Set-Up of the Column Experiments
To understand the behavior of Mn at the RBF site in Dresden-Tolkewitz, three columns with riverbed sediments were set up in the laboratory of the University of Applied Sciences Dresden (Figure S1).Each column was 1 m long, had an inner diameter of 0.08 m, and was made up of galvanized steel.The filling material was riverbed sediment from the Elbe River.The riverbed material was collected in front of the investigated monitoring cross-section of the RBF Waterworks Dresden-Tolkewitz (Figure 1).The riverbed sediment was recovered around 20 m apart from the riverbank during a low discharge period.Due to the very low gradient of the riverbed towards the riverbank, the area around the sampling point was already flooded at slightly higher water levels, which still would occur during mean low flow conditions.Hence, the area around the sampling point can be considered as a potential infiltration area.Because of the relatively coarse riverbed, undisturbed sampling was not possible.Thus, the upper 5 cm of the riverbed where scratched first to represent the clogging layer.Subsequently, the deeper riverbed material was dug out layer-wise and sieved in place to a grain size <4 mm.Immediately after transporting it to the lab, the wet riverbed material was filled into columns in ≈0.1 m thick, separately compacted layers with the clogging layer on top.
During the filling of the columns, the sediment mass was measured using a balance to calculate the bulk density and assess the compaction of the material in the columns.The mean travel time (t a ) and effective porosity (n e ) for both columns were determined from electrical conductivity (NaCl) breakthrough curves from tracer experiments performed before start-up.
To adjust the temperature regiment, all three columns with riverbed material were stored in a thermostatic cabinet.A second thermostatic cabinet contained three storage containers with Elbe river water, collected in Dresden once per week.The outflow of the columns flowed into three additional containers within the second cabinet.The investigated temperatures were 10, 20, 30, and 35 • C (Table 1).Temperature (T), dissolved oxygen (DO), the pH-value, and electrical conductivity (EC) were determined using WTW Multi 3430 and appropriate electrodes (WTW, Weilheim, Germany) before the columns in the storage containers and after the columns in a flow-through cell.A series of valves allowed sending the outflow of each column separately through the cell.
At most RBF sites, the infiltration rates depend for example on abstraction rates of the wells, clogging of the riverbed, distance between the river and the wells and the infiltrating area, and are therefore very site specific.To represent the low, mean, and high infiltration rates of 0.3, 0.6 and 1.1 m 3 /(m 2 •d), the flow through the columns was adjusted to 1, 2, and 4 mL/min.The flow rate was adjusted individually for each column using ProMinent Beta diaphragm pumps (ProMinent, Heidelberg, Germany).Each of the 12 possible flow and temperature conditions run until 15 to 20 pore volumes (PV) of every column were exchanged.Sampling started after around 5 PV and continued until at least 10 and 15 PV.Up to three intermediate samples were taken if possible (e.g., weekends were skipped).At one sampling event, the water samples from all three columns were taken at separate sampling taps after the columns before the outflow container.Thus, the presented results for each temperature and flow rate represent the mean value of three similarly prepared, and independently operated columns.Samples from the storage containers (=inflow) were taken once per week.Alkalinity was determined at every sampling event by alkalimetric titrations with 0.1 M hydrochloric acid (HCl).

Sequential Extraction of the Riverbed Sediment
To estimate the mobilization behavior of Mn and to assess the mineralogical composition, a 4-step sequential extraction procedure was applied to the filling material of the columns after the experiment (Table S1).Rauret et al. [27] and Sutherland and Tack [28] described the applied procedure in detail.The total extractable Mn was determined by microwave acid digestion with HNO 3 for separate samples from the same sampling points.Samples were taken after 0.05 m (below the clogging layer), and at 0.3, 0.6, and 0.9 m along the columns, before being immediately filled into airtight sample containers and stored at 4 • C before analysis.

Estimation of Reduction Constants for the Elbe Riverbed with PHREEQC
Elevated Mn concentrations at many bank filtration sites are linked to the microbiological reduction of Mn minerals within the riverbed and the aquifer [1,29].The degradation (oxidation) of organic matter (OM, simplified CH 2 O) is the driving force for the associated redox reactions (Equations ( 1)-( 3)).
Aerobic respiration : Mn(IV) reduction : The results from the column experiments are considered to represent the potential Mn release from the riverbed.In order to use the results for a planned modeling of the transect, the observed Mn release was reproduced by chemical modeling with PHREEQC [30].By applying the approach of Henzler et al. [14], the relevant redox reactions were modeled as kinetic reactions using Monod-type rate formulations (Equations ( 4)-( 7)).Because neither increasing Fe concentrations nor decreasing sulfate concentrations were observed along the transect, additional redox reactions accounting for iron and sulfate reduction were excluded.
The parameters r ox , r nit and r mn denote the production and consumption rates (positive and negative) of dissolved O 2 , NO 3 − , and Mn 2+ .Rate constants for OM degradation under oxic, nitrate and manganese reducing condition are represented by k ox , k nit and k mn .C ox and C nit are the dissolved oxygen and nitrate concentrations.K ox and K nit denote Monod-half saturation constants.The inhibition of nitrate and manganese reduction under oxic conditions was implemented by the inhibition constants K ox inhb nit and K nit inhb mn .Accordingly, K nit inhb mn represents the inhibition constant for manganese reduction under nitrate-reducing conditions.The overall turnover rate of OM r OM (Equation ( 7)) is the sum of the reaction rates r ox , r nit , and r mn that are multiplied with the stoichiometric coefficients Y ox , Y nit , and Y mn corresponding to the redox reactions (Equations ( 1)-( 3)).Following a similar modeling approach of Greskowiak et al. [31], the parameter f reac was included to simulate a zone of increased reactivity at the first section of the infiltration path [14].The application of Equations ( 4)-( 7) implied two assumptions.First, MnO 2 was present in excess.Hence, MnO 2 was not rate limiting and the implementation of a Monod-half saturation constant for Mn was not necessary in Equation 6.Second, the OM content was assumed to be infinitely available (or redelivered by the river) and would not be exhausted during the simulation period.Similar to Diem et al. [32], Greskowiak et al. [31], and Sharma et al. [33], an additional temperature factor f T was implemented that accounted for the impact of temperature changes on the degradation rates (Equation ( 8)).
T opt denotes the optimal temperature for a maximal degradation rate and α as well as β are fitting parameters.Applying the results from Diem et al. [32], Henzler et al. [14], Greskowiak et al. [33] and Sharma et al. [33], none of the three parameters were to be changed from the initial data set.
In PHREEQC, a 1 m long column was represented by 50 cells with a length of each of them being 0.02 m (Table S2).Porosity and pore velocity were known from the tracer test.The dispersion and diffusion coefficients for the model were calibrated for the NaCl breakthrough curves from the tracer tests and non-reactive transport.Subsequently, the data from Henzler et al. [14] for the rate constants k ox , k nit , and k mn , as well as for the inhibition constants K ox inhb nit , K ox inhb mn , and K nit inhb mn , were used as the initial parameter set for reactive modeling.Calibration was initially carried out with PEST [34].Since the inhibition constants K ox inhb nit , K ox inhb mn and K nit inhb mn did not change during the initial calibration runs, and in order to speed up the calibration, the inhibition constants were held fixed at the initial values during further calibration.The following calibration of the rate constants k ox , k nit and k mn was first carried out by adjusting the parameters for best fit by hand (trial-and-error).Afterward, the trial-and-error results were checked with PEST.
The calibration targets were the determined median values of pH, DO, NO 3 − and Mn 2+ in the outflow water of the columns for each of the three flow and four temperature conditions.Hence, the calibration each resulted in 12 values for k ox , k nit , and k mn .

Seasonal Fluctuation of Redox-Sensitive Parameters Close to the Riverbank
Most relevant redox parameters of the Elbe river undergo strong seasonal fluctuations (Table S3).The median value for water temperature was 10.9 Along the transect from the Elbe river towards the PW seasonal fluctuations were also noticeable (Table S4).As expected, OW 1-1 and 2-1 in the upper aquifer showed the strongest variations.The temperatures at both observation points varied from winter to summer from 6.1 to 20.0 • C and 8.3 to 19.9 • C (Figures 2 and 3).Temperature fluctuations at the deeper observation points OW 1-2 and 2-2 were in the same order of magnitude.At the lowest observation points OW 1-3 and 2-3, as well as further along the flow path at OW 3, the temperature variations were ±3 • C compared to the median values.
DO showed similar patterns during the year.The higher DO concentration in the Elbe river in winter resulted in >6.5 mg/L at OW 1-1, >2.5 mg/L at OW 1-2 and >1.0 mg/L at OW 2-1.Further, along the flow path, DO was mostly found <0.5 mg/L.During summer, the DO was almost depleted at OW 1-1.Nitrate showed a similar behavior.During winter, NO 3 − concentrations were >20 mg/L at OW 1-1, >7 mg/L at OW 1-2 and >13 mg/L at OW 2-1, whereas in summer, the nitrate was with <1.0, <0.5 and <3.0 mg/L almost neglectable.
Mn concentration varied primary at OW 1-1 during the year.In winter, the Mn concentration was very low at OW 1-1, but increased up to 0.42 mg/L during summer time.At all other observation points, variations were found but without distinct patterns.Along the entire flow path after OW 1-1, the Mn concentration was always >0.1 mg/L and mostly below 0.35 mg/L (Table S4).S4).S4).

Mn Release Depending on the Temperature and Infiltration Rate During the Column Experiments
To investigate the effect of varying temperature to the Mn release from the riverbed, three columns filled with riverbed sediments run at 10, 20, 30, and 35 °C.Varying flow rates of 1, 2, and 4 mL/min represented low, mean and high infiltration rates for all four temperature regimes.

Mn Release Depending on the Temperature and Infiltration Rate During the Column Experiments
To investigate the effect of varying temperature to the Mn release from the riverbed, three columns filled with riverbed sediments run at 10, 20, 30, and 35 • C. Varying flow rates of 1, 2, and 4 mL/min represented low, mean and high infiltration rates for all four temperature regimes.
After changing temperature and/or flow, the Mn concentration changed within 5 pore volumes (PV) and was stable after 8-10 PV (e.g., Figure S4).At high infiltration rates (4 mL/min, pore velocity v a = 3.44•10 −5 m/s, residence time t R = 8.1 h) and water temperatures of 10 and 20 • C, the median Mn release was <<0.01 mg/L (n = 9 and 12) and almost neglectable (Figure 5, Table S6).At 30 • C, the Mn concentration increased in the outflow slightly to around 0.03 mg/L (median, n = 21).Only at 35 • C was Mn released (median of 0.51 mg/L).

Mn Release Depending on the Temperature and Infiltration Rate During the Column Experiments
To investigate the effect of varying temperature to the Mn release from the riverbed, three columns filled with riverbed sediments run at 10, 20, 30, and 35 °C.Varying flow rates of 1, 2, and 4 mL/min represented low, mean and high infiltration rates for all four temperature regimes.
The Mn concentration in the feed water maximally was 0.02 mg/L (median <LOD, limit of detection 0.005 mg/L Mn, n = 23).

Reduction Constants of O 2 , NO 3
− and Mn(IV) as Electron Acceptors The observed Mn release was reproduced by chemical modeling with PHREEQC using Monod-type rate formulations and the reduction rate constants k ox , k nit , and k mn as calibration parameters.The reduction rate constant k ox showed the highest value at high infiltration rates and at a temperature of 10 • C (1.8 × 10 −9 mol/(L•s)).With increasing temperature, k ox decreased down to 4.3 × 10 −10 mol/(L•s) (Table 2, Figure S5).For mean and low infiltration rates, a similar behavior was found.The lowest k ox of 9.7 × 10 −11 mol/(L•s) was determined at low infiltration rates and a temperature of 35 • C. For the reduction rate constant k nit , no distinct pattern was found.At high infiltration rates, k nit remained constant at 5.0 × 10 −12 mol/(L•s).For mean and low infiltration rates at 10 • C, k nit was in the same order of magnitude (7.0 and 2.0 × 10 −12 mol/(L•s)).With increasing temperatures of 20 • C and 30 • C, k nit increased up to 2.8 × 10 −10 mol/(L•s).At 35 • C, k nit decreased again at the mean and low infiltration rates.
The reduction rate constant k mn at 10 • C and 20 • C for high, mean and low infiltration rates were in the order of magnitude of 1.0 to 2.5 × 10 −10 mol/(L•s).At 30 • C and high infiltration rates, k mn remained in this range.For mean and low infiltration rates at 30 • C and 35 • C, k mn increased to 1.4-2.0× 10 −9 mol/(L•s).
With a percental error of −2.4% compared to the measured values (median, −6.7-0.2%,10-90%ile, n = 12), the simulated DO concentrations showed the errors of relevant redox parameters (Table S7).The error for NO 3 − with a median of 0.4% and a span of −0.6 to 1.1% (10 to 90%ile, n = 12) was lower.With an error of 0.2%, the simulated Mn 2+ concentrations showed the smallest median deviation but a comparable large span of −15.7 to 8.5% (n = 12).The median error for the simulated pH was −1.0% with a span of −8.7 to 1.7% (n = 12).

The Decrease of Easily Reducible Mn Along the Flow Path
To estimate the mobilization behavior of Mn, a 4-step sequential extraction procedure was applied to the filling material of the three columns after the experiment.All three columns showed qualitatively similar results and all the following values represent median values with n = 3 (Table S8).The total Mn mass (Mn tot , as the sum of all 4 extracted fractions) was around 270 mg/kg at the inlet after 0.05 m (Figure 6).Further along the flow path, the total Mn mass decreased to ≈150 mg/kg down to a minimum of ≈125 mg/kg the minimum at the outlet after 0.9 m.The total extractable Mn (microwave acid digestion) was around 112 mg/kg at the inlet and increased to 133, 250, and 473 mg/kg at the outlet (Table S8).With a percental error of −2.4% compared to the measured values (median, −6.7-0.2%,10-90%ile, n = 12), the simulated DO concentrations showed the largest errors of relevant redox parameters (Table S7).The error for NO3 − with a median of 0.4% and a span of −0.6 to 1.1% (10 to 90%ile, n = 12) was lower.With an error of 0.2%, the simulated Mn 2+ concentrations showed the smallest median deviation but a comparable large span of −15.7 to 8.5% (n = 12).The median error for the simulated pH was −1.0% with a span of −8.7 to 1.7% (n = 12).

The Decrease of Easily Reducible Mn Along the Flow Path
To estimate the mobilization behavior of Mn, a 4-step sequential extraction procedure was applied to the filling material of the three columns after the experiment.All three columns showed qualitatively similar results and all the following values represent median values with n = 3 (Table S8).The total Mn mass (Mntot, as the sum of all 4 extracted fractions) was around 270 mg/kg at the inlet after 0.05 m (Figure 6).Further along the flow path, the total Mn mass decreased to ≈150 mg/kg down to a minimum of ≈125 mg/kg the minimum at the outlet after 0.9 m.The total extractable Mn (microwave acid digestion) was around 112 mg/kg at the inlet and increased to 133, 250, and 473 mg/kg at the outlet (Table S8).The mass fraction of soluble and carbonate bound Mn ("Carb.")decreased from around 65 mg/kg to 30 mg/kg along the flow path.On a percentage base, the soluble and carbonate bound Mn always remained between 21-24% (Figure 7).Easily reducible Mn dropped from around 140 mg/kg (53% of the total mass) at the inlet to 36 mg/kg (<30%) after 0.9 m.Organically bound Mn ("Organ.")remained stable between 6-9 mg/kg (3-5%).The residual Mn fraction was found almost constant at 51-64 mg/kg, but the percentage increased from around 20 to >40%.The mass fraction of soluble and carbonate bound Mn ("Carb.")decreased from around 65 mg/kg to 30 mg/kg along the flow path.On a percentage base, the soluble and carbonate bound Mn always remained between 21-24% (Figure 7).Easily reducible Mn dropped from around 140 mg/kg (53% of the total mass) at the inlet to 36 mg/kg (<30%) after 0.9 m.Organically bound Mn ("Organ.")remained stable between 6-9 mg/kg (3-5%).The residual Mn fraction was found almost constant at 51-64 mg/kg, but the percentage increased from around 20 to >40%.

The Significance of the Calibrated Degradation Rate Constants
In order to use the results from the column experiment for a planned modeling of the transect, the observed Mn release was simulated using PHREEQC, focusing on relevant redox reactions.During an initial automated calibration with PEST, the inhibition constants K ox inhb_nit, K ox inhb_mn, and K nit inhb_mn did not change and were held fixed at the initial values during further calibration.Hence, only the reduction rate constants kox, knit and kmn were used for calibration.
The calibrated kox had an overall median value of 4.4 × 10 −10 mol/(L•s) (n = 12, Table 2).These results are in fair agreement with results of Henzler et al. [14], who determined a kox of 2.0 × 10 −10 mol/(L•s) for a bank filtration site at Lake Tegel in Berlin, Germany.Greskowiak et al. [31] used a kox of 1.52 × 10 −10 mol/(L s) to simulate the consumption of DO at an infiltration pond in Berlin.Sharma et al. [33] used a more than ten times higher degradation rate with kox = 1.30 × 10 −9 mol/(L•s) to simulate a transect at the RBF Waterworks Flehe in Düsseldorf, Germany.The average water temperature in their study was 13.5 °C and the BF traversed the ≈60 m wide stretch from the riverbank to the well within 60 days [33].This corresponds to a pore velocity (va) of ≈1.2 × 10 −5 m/s.Comparing this with the results of this study of kox = 8.65 × 10 −10 mol/(L•s) at 10 °C and va ≈ 1.72 × 10 −5 m/s, shows a good agreement, too.
For kmn, the overall median value was 2.3 × 10 −10 mol/(L•s) (n = 12, Table 2).This is ≈100-to 1000times higher than the reported data of kmn = 1.7 × 10 −12 mol/(L•s) [14], 8.9 × 10 −13 mol/(L•s) [31] or 3.0 × 10 −13 mol/(L•s) [35].Matsunga et al. [35] modeled a column experiment and included the precipitation of Rhodochrosite (MnCO3) as a sink for Mn 2+ .The precipitation of MnCO3 as a sink for Mn 2+ has also been shown at other bank filtration sites [38].Taking the average water quality of the outflow from the columns at 20 °C (pH 8, 50 mg/L Ca 2+ , 120 mg/L HCO3 − ) and mean infiltration rates for equilibrium in PHREEQC, the water would be supersaturated with respect to MnCO3 with around 0.1 mg/L Mn 2+ in equilibrium (saturation index 0.77).Thus, the precipitation ofMnCO3 would probably control the Mn concentration for longer residence times but did not along the columns due to slow reaction kinetics [39].As a consequence, the given kmn seem to be representative in simulating an Mn release from a (highly reactive) riverbed, but are not suited for simulations of longer transects, in which an Mn reduction is likely to occur within the aquifer and can be controlled by MnCO3 precipitation, for example.

The Significance of the Calibrated Degradation Rate Constants
In order to use the results from the column experiment for a planned modeling of the transect, the observed Mn release was simulated using PHREEQC, focusing on relevant redox reactions.During an initial automated calibration with PEST, the inhibition constants K ox inhb_nit , K ox inhb_mn , and K nit inhb_mn did not change and were held fixed at the initial values during further calibration.Hence, only the reduction rate constants k ox , k nit and k mn were used for calibration.The calibrated k ox had an overall median value of 4.4 × 10 −10 mol/(L•s) (n = 12, Table 2).These results are in fair agreement with results of Henzler et al. [14], who determined a k ox of 2.0 × 10 −10 mol/(L•s) for a bank filtration site at Lake Tegel in Berlin, Germany.Greskowiak et al. [31] used a k ox of 1.52 × 10 −10 mol/(L s) to simulate the consumption of DO at an infiltration pond in Berlin.Sharma et al. [33] used a more than ten times higher degradation rate with k ox = 1.30 × 10 −9 mol/(L•s) to simulate a transect at the RBF Waterworks Flehe in Düsseldorf, Germany.The average water temperature in their study was 13.5 • C and the BF traversed the ≈60 m wide stretch from the riverbank to the well within 60 days [33].This corresponds to a pore velocity (v a ) of ≈1.2 × 10 −5 m/s.Comparing this with the results of this study of k ox = 8.65 × 10 −10 mol/(L•s) at 10 • C and v a ≈ 1.72 × 10 −5 m/s, shows a good agreement, too.
For k mn , the overall median value was 2.3 × 10 −10 mol/(L•s) (n = 12, Table 2).This is ≈100to 1000-times higher than the reported data of k mn = 1.7 × 10 −12 mol/(L•s) [14], 8.9 × 10 −13 mol/(L•s) [31] or 3.0 × 10 −13 mol/(L•s) [35].Matsunga et al. [35] modeled a column experiment and included the precipitation of Rhodochrosite (MnCO 3 ) as a sink for Mn 2+ .The precipitation of MnCO 3 as a sink for Mn 2+ has also been shown at other bank filtration sites [38].Taking the average water quality of the outflow from the columns at 20 • C (pH 8, 50 mg/L Ca 2+ , 120 mg/L HCO 3 − ) and mean infiltration rates for equilibrium in PHREEQC, the water would be supersaturated with respect to MnCO 3 with around 0.1 mg/L Mn 2+ in equilibrium (saturation index 0.77).Thus, the precipitation ofMnCO 3 would probably control the Mn concentration for longer residence times but did not along the columns due to slow reaction kinetics [39].As a consequence, the given k mn seem to be representative in simulating an Mn release from a (highly reactive) riverbed, but are not suited for simulations of longer transects, in which an Mn reduction is likely to occur within the aquifer and can be controlled by MnCO 3 precipitation, for example.

Impact of the Discharge on the Observed Redox Patterns
River discharge indirectly affects the quality of the BF in multiple ways.Low river discharge and related low water levels are often associated with less dilution of wastewater effluent, higher loads of dissolved organic matter (DOM), prolonged travel times, and they are expected to promote DO consumption [19].Low discharge periods further result in reduced shear stress at the riverbed and more intense clogging at the bottom of the water body, which can also promote anoxic conditions [40,41].High discharge conditions/floods can cause a partial removal of the clogging layer, resulting in a better hydraulic conductivity of the riverbed and higher water levels lead to higher infiltration rates [18].Increased flow velocities combined with shorter flow paths cause shorter travel times.Furthermore, high flow conditions and shorter travel times were linked to an increased input of TOC (total organic carbon) [42,43].
To the authors' knowledge, very few researchers have addressed the impact of river discharge to Mn release from the riverbed or the Mn concentration in the BF.Previous work has focused on the impact of discharge on the degradation of natural organic matter (NOM) and DO consumption.In a column experiment, von Rohr et al. [44] evaluated the role of discharge on NOM degradation during RBF.The experiment was set up at 20 • C and with four flow rates, resulting in residence times (t R ) of 40, 20, 12, and 4 h.At residence times of 4 h, oxic conditions persisted through the entire column.With a residence time of 20 h, anoxic conditions were observed within 20 cm along the flow path.Results from the column experiment in this study confirm the flow/infiltration rate dependency of oxygen depletion.At 20 • C and at high infiltration rates (t R = 8.1 h), the columns remained oxic, but became anoxic at residence times of 16-32 h (Table S9).Diem et al. [32] found river discharge to be correlated with an enhanced POM input and higher DO consumption during flood events.Due to a generally higher DO consumption at temperatures >15 • C, the correlation was only found to be true for temperatures below 15 • C [17].Diem et al. [32] successfully modeled these observations using higher DO consumption for higher discharges.In this study, rates for DO consumption (k ox ) in PHREEQC were highest at high infiltration rates and decreased at mean and low infiltration rates, which is in agreement with the findings of Diem et al. [32].
Groffman and Crossey [45] found slightly increasing Mn concentrations within the upper aquifer section during periods with lower discharge at Rio Calaveras (New Mexico, USA).The observations from the transect in Dresden-Tolkewitz did not reveal a clear discharge/water level dependency.Figure 4 indicates an increasing Mn concentration at OW 1-1 during a long-lasting low discharge period of the Elbe River in 2015.In addition, Mn at OW 2-1 and 3-1 increased during this period (Figures S2 and S3).Contrary to this, at OW 1-1 Mn, already decreased when the water level was still low, whereas Mn further increased at OW 2-1 and 3-1.Hence, the field observations did not indicate a mere discharge-dependency of Mn along the transect.
Manganese release during the column experiment with riverbed material was low at 10 • C and 20 • C for the high, mean, and low infiltration rates (Figure 5).At 30 and 35 • C, Mn increased from high to mean and again to low infiltration rates.Correlation coefficients of −0.99 for 30 • C and −0.87 for 35 • C were very high and significant (p < 0.05, data not shown).
The relation between river discharge and the release of Mn from the riverbed is ambiguous.On the one hand, the literature data [44] and the results from the column experiments in this study suggest that periods of low discharge lead to elevated Mn concentrations due to lower infiltration rates and longer residence times.Contrary to that, e.g., Diem et al. [17] found higher river discharge to be correlated with higher DO consumption, which was also reinforced by modeling [32].Since higher DO consumption is often related to increasing Mn concentrations during RBF [1], periods of higher discharge may cause Mn release.Shorter travel times during periods of high river discharge would interfere with this.In addition, the results of the column experiment in this study showed that the temperature effect distorts this discharge dependency (see the following paragraph).
Manganese along the transect in Dresden-Tolkewitz varied the most at OW 1-1, 2-1, and 2-2 (Figures 4 and S2), which are closest to the Elbe river (≈20 m).Furthermore, those observation points showed the largest temperature fluctuations (Figures 2 and 3).Especially the Mn concentration at OW 1-1, which is closest to the riverbank, responded to higher water temperatures during low and mean flow conditions (Figures S6-S8).However, at the more distant observation points (OW 1-3, 2-3 and OW 3), Mn was comparably stable during the year.This can partially be attributed to extended travel times (compared to the average travel times), especially towards the deeper observation points.Similar observations were made at RBF sites at the Lot river in France [51] and at the Glatt river in Switzerland [50].Both found a seasonal trend for Mn, with elevated concentrations during the warm summer months and lower concentrations during the winter season.For the RBF site in Dresden-Tolkewitz, no seasonal trend for Mn was observed at any OW, which is probably due to the long sampling interval.Temperature seems to be a very important factor controlling initial Mn release from the riverbed in Dresden-Tolkewitz, but more research is needed in order to assess this parameter quantitatively.Since the assessment is only possible during temporary constant, low water levels [23], future research must focus on the rare long-lasting low discharge periods.
Results from the column experiments indicate a high-temperature dependency of Mn release.At constant infiltration rates (high, mean, or low), the Mn concentration increased with increasing temperature (Figure 5).At high infiltration rates (residence time t R = 8.1 h), the Mn concentration increased in the outflow sharply from 0.03 to 0.51 mg/L at higher temperatures of 30 and 35 • C. For mean and low infiltration rates (t R = 16.2 and 32.4 h), a similar Mn increase happened after water temperature increased from 20 to 30 • C. The calibrated k mn in this study reflects the sharp increases, showing a 20-fold increase from 30 to 35 • C at high infiltration rates as well as a 15-and 17-times higher k mn after temperature raised from 20 to 30 • C at mean and low infiltration rates (Table 2).
Bourg and Bertin [51] reported a threshold water temperature of 10 • C to trigger a microbiologically mediated Mn reduction.Hoehn et al. [52] observed extensive denitrification, stronger reducing conditions, and elevated Mn concentrations at water temperatures above 14 • C. The findings of the column experiment support these observations.At water temperatures of 10 • C, the Mn release was neglectable for all infiltration rates.In addition, the results from the column experiment and the determined reduction constants can expand these statements for Mn containing riverbeds.For high infiltration rates of 1.1 m 3 /(m 2 •d) and above, Mn release from the riverbed is unlikely for surface water temperatures that are typical in temperate climate zones (≤30 • C).If the BF infiltrates at infiltration rates in the order of magnitude of 0.6 m 3 /(m 2 •d) and below, water temperatures >20 • C are sufficient to trigger extensive Mn release from the riverbed.

Depletion of the Mn Reservoir Within the Riverbank
After the column experiments were finished, a sequential extraction procedure was applied to the riverbed sediment (Figures 6 and 7).
Manganese contents in the riverbed sediments are very site-specific and cannot be limited to regional or climatic differences [53].Moreover, the Mn content depends highly on the grain size.Jain and Ram [54] and Jain and Sharma [55] found 230-650 mg/kg Mn in the grain size fraction with the highest mass fraction (36%, d = 0.21-0.25 mm), whereas the grain size fraction with the lowest mass fraction of only 2% (d < 0.075 mm) contained up to 2800 mg/kg Mn.
However, the order of magnitude of ≈150-300 mg/kg Mn was low compared to many other riverbed sediments, e.g., 960 mg/kg Mn in the Rhine sediment [56] or 1700 mg/kg Mn in the Garonne [53].Nonetheless, the values are in agreement with the former analysis of the riverbed sediment from the Elbe river.Grischek et al. [57] determined ≈50 mg/kg Mn tot (n = 3) in riverbed sediment at Torgau around 100 km downstream of Dresden.In Dresden-Tolkewitz at the same transect, Paufler [58] determined around 250 mg/kg Mn tot (n = 3) for the upper 5 cm of the riverbed and only ≈80 mg/kg Mn tot (n = 6) for the depths 5-30 cm.
At the Glatt river in Switzerland, von Gunten et al. [13] hypothesized that the repetitive exhaustion of Mn deposits in the river sediments towards the fall season could be one reason for seasonal Mn variations.
Integrating the released Mn from one column over the entire experimental time in this study (all four infiltration rates and three temperature regimes), this results in an overall mass output of approx.140 mg Mn (data not shown).After the experiment, the total Mn mass was around 1600 mg within a column (Figure 7, ≈10 kg of sediment, calculated dry weight).Considering the large deviations from the total extractable Mn (Table S8), and the potential error of sequential extraction procedures [59], conclusions about a possible exhaustion of Mn during the column experiment would not be reasonable.Taking the observations of Paufler [58] into account, the observed steps along the columns in this study may reflect the natural stratification within the riverbed, rather than a depletion of the Mn reservoir.
With around 5 mg/L of filterable substances and a Mn content of 3500 mg/kg of the suspended matter (both median, n > 100 [25]), the Elbe river water contains approximately 0.04 mg/L Mn, which is bound to suspended solids.Considering the mean abstraction of 22,000 m 3 /d in Dresden-Tolkewitz and a bank filtrate portion of 70%, about 15,500 m 3 BF are abstracted per day.Assuming an even distribution of suspended solids in the infiltrating river water, the Elbe river delivers around 6 × 10 5 mg Mn per day into the riverbed, which is probably not entirely reducible.With the same BF portion and the observed Mn release of 0.1 to 0.8 mg/L, the Mn output from the riverbed into the BF would be 1.5 × 10 6 to 1.2 × 10 7 mg/d.Thus, a depletion of the Mn reservoir seems to be possible but more research is needed to evaluate this conclusively.

Implications for (River-)Bank Filtration Sites
Current climate forecasts show an increase in river discharge seasonality, with increasing high discharges and decreasing low discharges.Furthermore, global mean river water temperatures are projected to increase by 0.8-1.6 • C by 2100 with Europe, the United States, parts of southern Africa, Australia and eastern China facing the largest changes [60].
The results of this study showed stronger Mn release for low infiltration rates compared to high infiltration rates when the temperature was equal.RBF sites in regions with a colder climate probably will not notice different Mn concentrations in the future, since Mn release was neglectable at water temperatures of 10 • C for all infiltration rates.In regions with a temperate climate, Mn release from the riverbed is unlikely during high discharge events with infiltration rates ≥1.1 m 3 /(m 2 •d).At infiltration rates of around 0.6 m 3 /(m 2 •d) and below, water temperatures >20 • C are sufficient to trigger extensive Mn release from the riverbed.Thus, RBF sites in temperate climate zones at rivers with large seasonal discharge fluctuations may have to deal with increasing Mn concentrations in the future.Additionally, the effect of seasonal rivers showing the highest water temperatures during low discharge periods [60] may intensify this effect.Surface water temperatures of >30 • C are rare in the temperate climate zone.
Typically, rivers in dry/arid or tropical climates like the Nile river [61] or the Mekong river [62] show temperatures in this order of magnitude.Thus, the results of this study suggest that RBF sites in arid or tropical climate zones should be aware of an intensified Mn release from the riverbed in the future, which can distort the quality of the BF even at high infiltration velocities.In such cases, biological post-treatment could be a viable option to remove Mn from the BF [63].
Nonetheless, elevated temperatures are not necessarily causing elevated Mn concentrations during RBF.During the summer in 2003, maximum river water temperatures of 25 • C were observed at the RBF site at the Rhine river in Germany.No elevated Mn concentrations were detected, although anaerobic conditions developed within the aquifer [64].
The field observations in Dresden-Tolkewitz and at other RBF sites prove that an actual high Mn release from the riverbed must not necessarily lead to elevated Mn concentrations in the pumped raw water [1,47].Sorption, (re-)oxidation, and precipitation along the flow path are potential sinks for Mn released from a riverbed.Especially if the aquifer material contains Mn(hydr)oxides, Mn concentrations of the BF can change largely along the flow path due to the high affinity of Mn 2+ for Mn(hydr)oxides [63].Thus, a long distance between the riverbank and the pumping wells could be of advantage to buffer elevated Mn concentrations and temperature.However, such an advantage may have to be checked against a lower portion of abstracted bank filtrate and a higher portion of potentially Fe-and Mn-rich landside groundwater [65].
For the RBF site in Dresden-Tolkewitz, the temperature was found to be the driving force for Mn release from the riverbed and high discharge and infiltration rates limit the release at lower temperatures.In Dresden-Tolkewitz, the mean infiltration rate is around 0.2 m 3 /(m 2 •d), which corresponds to the infiltration rate that was found to be sustainable for RBF sites along the River Rhine, Elbe, and other European Rivers [66].After Soares [67], an infiltration rate of 0.32 m 3 /(m 2 •d) would be still sustainable for this RBF site.Ahrns [23] observed infiltration rates up to 0.95 m 3 /(m 2 •d) in Dresden-Tolkewitz.Thus, the investigated low infiltration rate of 0.3 m 3 /(m 2 •d) in this study represents the mean infiltration rate in Dresden-Tolkewitz.Such low infiltration rates must go in hand with temperatures above 20 • C to trigger a Mn release.With increasing temperatures, the infiltration rate becomes less important and at water temperatures around 30 • C, an extensive Mn release from the riverbed can be expected even at mean infiltration rates.

Conclusions
Current climate forecasts project increasing river discharge seasonality and water temperatures and thus, are thought-provoking at many (river)bank filtration sites.Water quality data of a 10-year time span from a monitoring transect at the RBF site Dresden-Tolkewitz and accompanying column experiments were used to assess the potential Mn release from the riverbed with respect to the climate-related variables, temperature and discharge.Temperature was found to be more important than discharge.Low infiltration rates 0.3 m 3 /(m 2 •d) required temperatures above 20 • C to trigger Mn release.With an increasing temperature, the discharge becomes less important and at 30 • C, the infiltration rates of ≈0.6 m 3 /(m 2 •d) can already cause an extensive Mn release from the riverbed.The subsequent modeling of the column experiment with PHREEQC resulted in degradation rates for DO that are applicable at other RBF sites for several water temperatures and flow velocities.The determined Mn reduction rate constants are appropriate to simulate an Mn release from riverbed sediments but are not suited for simulations in which Mn reduction is likely to occur within the aquifer and the Mn concentrations can be limited by precipitation, for example.S1: The chemical reagents and analytical conditions for the optimized BCR sequential extraction procedure [27,28]; Table S2: The input data of the column in PHREEQC;

•
Table S3: The statistical data of the Elbe river at Dresden for the entire observation period 2006-2016 [25]; • Table S4: The monthly and 10-year median values of selected parameters along all three OW's of the transect in the WW Dresden-Tolkewitz for the entire observation period 2006-2016; • Table S5: The statistical data for the low flow period in 2015; • Table S6: The statistical data for Mn release during the column experiment with riverbed sediment from the Elbe river in Dresden-Tolkewitz, Q 1 and Q 3 correspond to the 1st and 3rd quartile; • Table S7: The comparison of the measured parameter values during the column experiment and the modeled parameters from PHREEQC; • Table S8: The results and statistical data of the sequential extraction of the riverbed sediment after the column experiment (n = 3 for each statistical data); • Table S9: The median concentration of DO in the outflow during the column experiment.

Figure 1 .
Figure 1.(a) The location of the RBF Waterworks Dresden-Tolkewitz, (b) Observed transect between the Elbe river and the pumping well (PW), (c) Location of the observation wells within the transect and sampling point of the riverbed sediment for the column experiment, (d) cross-section of the transect and locations of the observation points for each observation well (OW).

Figure 1 .
Figure 1.(a) The location of the RBF Waterworks Dresden-Tolkewitz, (b) Observed transect between the Elbe river and the pumping well (PW), (c) Location of the observation wells within the transect and sampling point of the riverbed sediment for the column experiment, (d) cross-section of the transect and locations of the observation points for each observation well (OW).

19 Figure 2 .
Figure 2. The mean values (n = 2) of the relevant redox parameters along the transect in February (TableS4).

Figure 3 .
Figure 3.The mean values (n = 2) of the relevant redox parameters along the transect in September (TableS4).
water level of the Elbe river shows strong annual fluctuations.Long-lasting mean low discharge conditions are represented by a water level of ≤0.75 m and are rare.From 1998 to 2006, the Elbe water level was between 0.7 and 0.8 m only in 2003 (for almost 3 months).During the observation period from 2006 onwards, the Elbe river decreased to mean low discharge conditions in 2008, 2009, and 2016, which lasted for two weeks at maximum.In 2015, the latest low discharge period was observed that lasted for more than 3 months and the water level dropped to as low as 0.5 m [26].During this 153 day long low discharge period in 2015, the median water level was 0.74 m (0.62-1.05 m, 10-90%ile, n = 153) and the mean water temperature 20.1 °C (11.0-24.6 °C, n = 153, Table

Figure 2 .
Figure 2. The mean values (n = 2) of the relevant redox parameters along the transect in February (TableS4).
Water 2018, 10, x FOR PEER REVIEW 7 of 19

Figure 2 .
Figure 2. The mean values (n = 2) of the relevant redox parameters along the transect in February (TableS4).

Figure 3 .
Figure 3.The mean values (n = 2) of the relevant redox parameters along the transect in September (TableS4).
water level of the Elbe river shows strong annual fluctuations.Long-lasting mean low discharge conditions are represented by a water level of ≤0.75 m and are rare.From 1998 to 2006, the Elbe water level was between 0.7 and 0.8 m only in 2003 (for almost 3 months).During the observation period from 2006 onwards, the Elbe river decreased to mean low discharge conditions in 2008, 2009, and 2016, which lasted for two weeks at maximum.In 2015, the latest low discharge period was observed that lasted for more than 3 months and the water level dropped to as low as 0.5 m [26].During this 153 day long low discharge period in 2015, the median water level was 0.74 m (0.62-1.05 m, 10-90%ile, n = 153) and the mean water temperature 20.1 °C (11.0-24.6 °C, n = 153, Table

Figure 3 .
Figure 3.The mean values (n = 2) of the relevant redox parameters along the transect in September (TableS4).
The water level of the Elbe river shows strong annual fluctuations.Long-lasting mean low discharge conditions are represented by a water level of ≤0.75 m and are rare.From 1998 to 2006, the Elbe water level was between 0.7 and 0.8 m only in 2003 (for almost 3 months).During the observation period from 2006 onwards, the Elbe river decreased to mean low discharge conditions in 2008, 2009, and 2016, which lasted for two weeks at maximum.In 2015, the latest low discharge period was observed that lasted for more than 3 months and the water level dropped to as low as 0.5 m [26].During this 153 day long low discharge period in 2015, the median water level was 0.74 m (0.62-1.05 m, 10-90%ile, n = 153) and the mean water temperature 20.1 • C (11.0-24.6 • C, n = 153, Table Water 2018, 10, x FOR PEER REVIEW 8 of 19

Figure 4 .
Figure 4.The Mn concentration at OW 1 during a low discharge period in 2015.

Figure 5 .
Figure 5.The mean Mn concentration in the outflow of the columns with riverbed sediment for high, mean, and low infiltration rates (4, 2 and 1 mL/min) at different temperatures (see TableS6for no. of samples), Q1 and Q3 correspond to the 1st and 3rd quartile.

Figure 4 .
Figure 4.The Mn concentration at OW 1 during a low discharge period in 2015.

Water 2018 , 19 Figure 4 .
Figure 4.The Mn concentration at OW 1 during a low discharge period in 2015.

Figure 5 .
Figure 5.The mean Mn concentration in the outflow of the columns with riverbed sediment for high, mean, and low infiltration rates (4, 2 and 1 mL/min) at different temperatures (see TableS6for no. of samples), Q1 and Q3 correspond to the 1st and 3rd quartile.

Figure 5 .
Figure 5.The mean Mn concentration in the outflow of the columns with riverbed sediment for high, mean, and low infiltration rates (4, 2 and 1 mL/min) at different temperatures (see Table S6 for no. of samples), Q 1 and Q 3 correspond to the 1st and 3rd quartile.At mean infiltration rates (2 mL/min, v a = 1.72•10 −5 m/s, t R = 16.2 h), a minor Mn release was observed at 10 and 20 • C with median values of 0.01 and 0.02 mg/L, respectively.With increasing temperature, the Mn release increased to 0.49 mg/L at 30 • C (median, n = 9) and 0.63 mg/L at 35 • C (median, n = 15).

Figure 6 .
Figure 6.The median values (n = 3) of the absolute Mn mass fractions along the column.

Figure 6 .
Figure 6.The median values (n = 3) of the absolute Mn mass fractions along the column.

Figure 7 .
Figure 7.The median values (n = 3) of percentage Mn mass fractions along the columns.

Figure 7 .
Figure 7.The median values (n = 3) of percentage Mn mass fractions along the columns.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2073-4441/10/10/1476/s1,• Figure S1: The flow scheme of the column experiments inside the thermostatic cabinets; • Figure S2: The Mn concentration at OW 2 during a low discharge period in 2015; • Figure S3: The Mn concentration at OW 3 during a low discharge period in 2015; • Figure S4: The increase of the Mn concentration within the first 5 PV at 30 • C and 1 mL/min after lowering from 2 mL/min, Median, 10-and 90-%ile (each n = 3); • Figure S5: The calibrated reduction rate constants for the column experiment; • Figure S6: The Mn concentration at OW 1 during a low discharge period in 2007; • Figure S7: The Mn concentration at OW 2 during a low discharge period in 2007; • Figure S8: The Mn concentration at OW 3 during a low discharge period in 2007; • Table

Table 1 .
The experimental design of the column experiment.

Table 2 .
The calibrated reduction constants of this study compared to the literature data.