Introduction

Silicon is the second most abundant element in Earth’s crust and typically exists as silica (SiO2) or in silicate minerals. Dissolved silicate in a water environment is generally found as monomeric silicate (H4SiO4) and is characterized by a non-ionic form in neutral and weakly acidic solutions but is ionized in an alkaline solution (Philippini et al. 2006). The ionized ratio of dissolved silicate is approximately 0.016% at a pH of 7 and increases with increasing pH. This ratio is 50 and 94% at a pH of 9 and 11, respectively (Gottlieb and Meyers 2012). Overall, the average global concentration of dissolved silicate ranges from 0.07 mmol/L (2 mg-Si/L) in the ocean (Tréguer et al. 1995) to 0.14 mmol/L (4 mg-Si/L) in rivers (Conley 1997). High concentrations of dissolved silicate [e.g., > 1.42 mmol/L (40 mg-Si/L)] are also found in sediment pore waters (Hallmark et al. 1982).

Silicate is one of the essential nutrients used by aquatic organisms such as phytoplankton (e.g., diatoms) to form their skeletons (Perry and Keeling-Tucker 2000; Yee et al. 2003; Lacombe et al. 2007); this diatom accounts for 25% of the world’s entire net primary production (Willén 1991). Diatom growth competes with dinoflagellate growth, which causes red tide (Furumai 2012). Flagellates such as dinoflagellate, chrysophyte, and chlorophyte are potentially harmful phytoplankton in nature, and when diatoms disappear, they become dominant species (Schelske and Stoermer 1971; Conley et al. 1993). Thus, silicate is a critical nutrient that may influence primary production and dominant species change of phytoplankton (Officer and Ryther 1980; Conley et al. 1993), and a comprehensive understanding of the geochemical behavior of silicate is important in aquatic science and ecology.

Dissolved silicates are mainly supplied to lakes by rivers. A portion of this silicate is readily adsorbed onto ferric hydroxide [Fe(OH)3] which is distributed in an oxidative epilimnion to form relatively insoluble siliceous ferric hydroxide particles (Vempati and Loeppert 1989; Fischer and Knoll 2009; Gottlieb and Meyers 2012). These particles are then deposited in lake-bottom sediment. If the lake-bottom environment maintains oxidative properties, nearly all the particles precipitate and are buried semi-permanently in the bottom sediment. Therefore, this adsorption reaction influences the concentration of dissolved silicate in the water column (Edwards and Liss 1973). In addition, this causes a gradual depletion of dissolved silicate in some lake waters (Tallberg 2000).

At Lake Nakaumi, a brackish lake located in western Japan, concentrations of dissolved silicate increase seasonally with total hydrogen sulfide in the anoxic hypolimnion and pore water in bottom sediment, which is characterized by elevated values (i.e., a maximum of 6.96 mmol/L). These results suggest that dissolved silicate eluted from lake-bottom sediment in a highly anoxic-reducing environment.

In general, previous studies have proposed that dissolved silicate eluted from sediment originates from biogenic matter, such as diatom debris (van der Weijden 2007). However, Testa et al. (2013) noted that bacteria were not able to decompose diatom cell in an anoxic environment, and that elution of silicate adsorbed on ferric hydroxide would be main source for its supply to the pore water. In Lake Nakaumi, biogenic silica dissolution is also insufficient to explain the drastically increasing concentrations of silicate in an anoxic-reducing environment.

Silicate elution in reductive environments is considered to be one of the routes for the supply of silicate to a whole lake. Dissolved silicates, which are eluted from sediment, diffuse into the hypolimnion and circulate throughout the entire lake via vertical mixing processes. Aquatic organisms can reuse this silicate as a nutrient source (Lehtimäki et al. 2016) or the silicate may react with metal oxides and re-precipitates in bottom sediments, such that sediment plays an important role as both source and sink for silicate (Matisoff et al. 1985). Additionally, we hypothesized that hydrogen sulfide affects the elution of silicate adsorbed on ferric hydroxide in an anoxic-reducing environment. In this study, we investigated the mechanisms of silicate adsorbed on ferric hydroxide and silicate elution in the presence of hydrogen sulfide in the laboratory experiments.

Materials and methods

Sampling and chemical analysis

Sampling area and pretreatment for analysis

We investigated a dredged area in Lake Nakaumi, a brackish lake with an average depth of 5.4 m that lies between Japan’s Shimane and Tottori Prefectures (Fig. 1). In the 1960s, land reclamation projects here formed a dredged area with a maximum depth of approximately 17 m. As a result, this area is often stagnant and has anoxic-reducing hypolimnion throughout the year (Sugahara et al. 2015).

Fig. 1
figure 1

Sampling point in Lake Nakaumi. The dredged area is located at 35.43˚ N, 133.27˚ E

Water samples were collected with Kitahara-type water sampler (Rigo, type 5023-B) in July and August 2015. Sampling was performed between depths of 1 and 15 m at 2 m intervals. Water temperature, redox potential, dissolved oxygen (DO) content, turbidity, electric conductivity, and salinity were measured with a multi-water quality profiler (Hydrolab, Quanta and Datasonde5X). We measured pH with a pH meter (Horiba, type D-50).

Sediment cores were collected from the dredged area with a KK-type core sampler (diameter = 51 mm, length = 670 mm, Hashimoto Scientific Co., Ltd., Japan) in July 2015. All cores were sectioned into thicknesses of 1 or 2 cm. A pore water sample was extracted from each section based on the methods reported in Sugahara et al. (2010). 2.5 mL of sediment from each section was placed into a 10 mL polyethylene syringe. Each sample was mixed with 2.5 mL of anaerobic deionized water prepared via bubbling nitrogen gas. A 25 mm filtering unit (pore size: 0.45 μm, Whatman) was then attached to the tip of the syringe to filter the mixture. The filtrate was analyzed for the dissolved silicate and other constituent concentrations. The actual (Ca) and measured (in the filtrate; Cm) concentration were obtained using the following equations.

Water content (WC, %) of the wet sediment sample:

$$ WC = ~\frac{{TW - DW}}{{TW}} \times 100~\left( \% \right) $$
(1)

where TW is the total weight of a wet 10 mL sediment sample (g), and DW is the dry weight of a wet 10 mL sediment sample (g). The pore water dilution ratio (DR) in the syringe was calculated with the following equation:

$$ DR = ~\frac{{\frac{{WC}}{{100 - WC}} \times DW_{s} + 2.5}}{{\frac{{WC}}{{100 - WC}} \times DW_{s} }} $$
(2)

where DWs is the dry weight of the wet sediment sample collected in the syringe (g). The actual concentration in the pore water was calculated as follows:

$$ C_{a} = ~C_{m} ~ \times DR $$
(3)

The dry weight was measured after drying the wet sediment sample at 105 °C for 24 h. The concentrations of phosphate (PO43–) and total hydrogen sulfide (H2S + HS) were measured spectrophotometrically via the molybdenum blue (PO43–) and methylene blue methods (H2S), respectively (Sugahara et al. 2010). Dissolved silicate was analyzed using the modified molybdenum yellow method as described in the following sections. Dissolved iron concentration was determined with inductively coupled plasma atomic emission spectrometer (ICP-AES, Optima 5000Z, Perkin-Elmer). All water samples were analyzed after filtration with a GF/F filter (Whatman, UK).

Reagents

Highly purified Milli-Q water (Millipore, Tokyo, Japan) was used to prepare chemical standard and reagent solutions. Milli-Q water was also used to wash all bottles and experimental apparatus. All chemicals used were of guaranteed reagent grade.

A silicate standard solution (16.6 mmol/L = 1000 mg-SiO2/L) was prepared by fusing 0.5 g of silicon dioxide with 5 g of sodium carbonate anhydrate and dissolving the pellet in 500 mL of Milli-Q water. An ammonium molybdate solution (100 g/L) was obtained by dissolving ammonium heptamolybdate in Milli-Q water. A tartaric acid solution (200 g/L) was prepared by dissolving L + -tartaric acid in Milli-Q water.

Modified molybdenum yellow method for dissolved silicates

Total hydrogen sulfide reduces a portion of the Mo (VI) ions in a molybdosilicate complex to Mo (V), which changes its color from yellow to blue (Rich 2007). Total hydrogen sulfide is abundant in the anoxic bottom layer of a brackish lake. Lake Nakaumi’s hypolimnion typically contains approximately 0.2 mmol/L of dissolved total hydrogen sulfide, whose concentration is elevated in the sediment (Sugahara et al. 2015). Therefore, appropriate pretreatment requires the complete removal of total hydrogen sulfide. The procedure for an accurate spectrophotometric determination of dissolved silicate was as follows. Lake water samples were filtered with a GF/F filter (Whatman, UK). 1 mL of 2.4 mol/L HCl was added to 20 mL of water sample in a polypropylene bottle, and total hydrogen sulfide in the sample was removed by stirring for 10 min. 1 mL of 100 g/L ammonium molybdate solution was then added to the sample. After 15 min, 1 mL of 200 g/L tartaric acid solution was added to the mixture. After 5 min, analyte solution absorbance was measured at 420 nm with a Shimadzu model UV-1800 spectrophotometer (Japan).

For pore water sample analysis, we modified the volumes of sample and reagent solutions. 4 mL of 0.12 mol/L HCl was added to 0.5 mL of the filtered pore water sample. Total hydrogen sulfide was removed by stirring for 10 min. An ammonium molybdate solution (100 g/L, 0.25 mL) was then added to the sample. After 15 min, 0.25 mL of a 200 g/L tartaric acid solution was added to the mixture. Afterward, the procedures for sample treatment were identical to the lake water sample procedures.

Laboratory experiments

Adsorption and elution of silicate onto ferric hydroxide precipitates

Silicate and iron mixed solutions that contained 0.16 mmol/L SiO2 and 0 or 3.6 mmol/L FeCl3 were prepared (NaCl concentration = 3.5%). The solution pH was adjusted to 7.00 ± 0.05, which yielded a ferric hydroxide precipitate. After 1 day, 10 mL of the sample solution was filtered into a 50 mL polypropylene bottle with a GF/F filter (filtrate A), and the ferric hydroxide precipitate was collected on the filter. The Si concentration in filtrate A was measured using the molybdenum yellow method. The ferric hydroxide precipitate was placed into another 50 mL polypropylene bottle and 10 mL of a 3.1 mmol/L sodium sulfide solution with a pH of 7 and 3.5% NaCl was added to the precipitate. After standing for 1 day, the solution was filtered with a GF/F filter (filtrate B), and its Si concentration was also measured with the methods used for filtrate A.

Silicate elution from bottom sediment with sodium sulfide solution

In this experiment, we used sediment core samples at depths between 0 and 5 cm below the sediment surface. 35 mL of sediment sample was placed into a 50 mL polypropylene syringe. 15 mL of a sodium sulfide solution (0, 0.14, and 8.8 mmol/L) with a pH of 7 and 3.5% NaCl was injected into the syringe. The tip of the syringe was sealed to prevent the entry of atmospheric air and was placed in a water bath at 20 °C to inhibit airflow. After vigorous shaking for 10 min, 1 mL of the filtrate was obtained from the mixture for Si concentration analyses.

Results and discussion

Silicate concentrations in Lake Nakaumi

In July 2015, the halocline was observed at 2–4 m depth with dissolved oxygen depletion in layers deeper than 4 m (Fig. 2). The redox potential had a negative value below depths of 10 m. Total dissolved hydrogen sulfide was not detected in the epilimnion but increased with depth in the hypolimnion, where its concentration was 2.1 × 10−1 mmol/L at the bottom (15 m) (Fig. 3). The vertical distribution profile of dissolved silicate also had similar patterns as the total hydrogen sulfide. The dissolved silicate concentration was nearly constant in the epilimnion and increased with depth in the hypolimnion, which was characterized by a concentration of 0.078 mmol/L at a depth of 15 m.

Fig. 2
figure 2

Vertical profiles of (a) water temperature (WT, □) and salinity (open circle), and (b) dissolved oxygen (DO, ◇) and redox potential (ORP, open triangle) in the water column in the dredged area of Lake Nakaumi for July 2015

Fig. 3
figure 3

Vertical profiles of total hydrogen sulfide (filled triangle), phosphate (filled rectangle), and dissolved silicate (open circle) in a the water column and b sediments in the dredged area of Lake Nakaumi for July 2015

The total hydrogen sulfide concentration in sediment had a maximum concentration of 5.83 mmol/L at 2.5 cm depth, which began to decrease below 4.5 cm depth and was not detected below 20 cm depth (Fig. 3). In the upper sediment (0–2.5 cm), the dissolved silicate concentration increased from 0.83 to 1.51 mmol/L with depth. Highly concentrated dissolved silicate (3.29–3.86 mmol/L) was detected below 15 cm depth. Dissolved iron concentrations were around the quantification limit by ICP-AES and reliable data were not obtained. Its concentration showed very low value of 0.07 nmol/L at the depth of 1.5 cm with total hydrogen sulfide of 5.57 mmol/L. In August 2015, we observed even higher concentrations of dissolved silicate. Below 1.5 cm depth, the dissolved silicate concentration ranged from 3.25–6.96 mmol/L; the area weighted mean (4.22 mmol/L) was calculated from the vertically integrated silicate content based on the vertical distribution profile (Table 1). On the other hand, dissolved silicate concentration ranged 0.36–0.45 mmol/L in May 2015, and the area weighted mean was 0.41 mmol/L. Dissolved silicate concentration increased from May to August in the pore water.

Table 1 Concentration of silicates in the pore water in Lake Nakaumi sediments on 2015

The existence of high concentrations of dissolved silicate in the pore water of Lake Nakaumi has not yet been reported. Concentrations in pore water of similar marine sediments have been reported as 0.08–0.13 mmol/L for the Bermuda Rise (Fanning and Pilson 1971) and 0.04–0.12 mmol/L for the Atlantic Ocean (Schink et al. 1975), significantly lower than that of Lake Nakaumi.

Hypothesis regarding silicate elution in lake sediment

In this study, silicate concentrations increased with those of total hydrogen sulfide in the anoxic hypolimnion and elevated concentrations of dissolved silicate were detected in pore water sampled from the bottom sediment in the dredged area of Lake Nakaumi; summer conditions enhanced this phenomenon.

Dissolved silicate adsorbs onto ferric hydroxide precipitates at neutral pH values (Iler 1979; Swedlund and Webster 1999; Lehtimäki et al. 2016). In addition, ferric hydroxide is abundant in oxidative sediment found in brackish lakes along coastal sea areas (Kostka and Luther III 1994; Vigderovich et al. 2019). Based on this, we propose the following mechanisms as a hypothesis to explain the phenomenon observed in Lake Nakaumi. At oxidative conditions, iron exists as Fe3+(III), forms ferric hydroxide (Fe(OH)3), and precipitates in the bottom sediment. Then, dissolved silicate adsorbs onto ferric hydroxide. When bottom conditions change to anoxic, sulfate-reducing bacteria produce hydrogen sulfide. The newly-formed hydrogen sulfide reduces the ferric hydroxide, and the silicate within the hydroxide is eluted into the pore water. Ferrous iron, which forms via the reduction of ferric hydroxide, reacts with sulfide ions and ferrous sulfide precipitates in the hypolimnion and pore water. Okumura et al. (2009) reported that ferrous sulfide was saturated in the bottom layer of Lake Nakaumi. Since ferrous sulfide has an inferior affinity to silicate compared with ferric hydroxide (Philippini et al. 2006), eluted silicate does not adsorb onto ferrous sulfide and, subsequently, concentrates in the pore water.

Since seawater contains a large number of sulfate ions, sulfate reduction via microorganism activity occurs in the bottom layer of sea areas, even though there is not a complete depletion of dissolved oxygen. Rabalais et al. (2002) reported that sulfate reduction occurs when the dissolved oxygen concentration is less than 0.05 mg/L. When sulfate reduction produces sulfide ions, this results in the reduction of ferric hydroxide. Previous studies have reported that ferric hydroxides coated with ferrous sulfide were present in marine sediments (Roden and Urrutia 2002). These reports indicate that our hypothesis may satisfy the conditions described above.

Borch et al. (2007) reported that ferric hydroxide, whose surface is coated with strongly adsorbed oxyanions such as phosphate, is difficult to be reduced by biological processes. Silicate’s tetrahedral structure is the same as a phosphate molecule’s, and its properties (such as an adsorption reaction to ferric hydroxide) are very similar. Therefore, ferric hydroxide adsorbing silicates may stably remain in a reducing environment because of their resistance to microbial diagenesis. Also, Lehtoranta et al. (2009) reported that lake eutrophication accelerates the conversion of dominant species from iron-reducing bacteria to sulfate-reducing bacteria, allowing for the production of hydrogen sulfide from the unfinished reduction of iron oxide. This idea is supported by previous studies on the reaction between ferric hydroxide and hydrogen sulfide in anoxic reducing sediments of the brackish lakes and the oceans (Canfield 1989; Fukusawa et al. 1995; Weston et al. 2006). Furthermore, Anschutz et al. (1998) reported that iron oxide existed in sediment where biological sulfate reduction was progressing, and Thamdrup et al. (1994) reported that amorphous iron (III) coexisted with hydrogen sulfide below 5 cm sediment depth, indicating an anoxic environment. Therefore, iron reduction in the anoxic-reducing sediment of Lake Nakaumi may be caused by hydrogen sulfide because this is known to be a eutrophic lake with a hypoxic water layer (Kusunoki and Sakata 2018) and dissolved silicate concentrations in the pore water exceeded the solubility of diatoms in the same condition.

Consequently, we suggest that silicate adsorbed onto ferric hydroxide is eluted by a reaction between hydrogen sulfide and ferric hydroxide in an anoxic-reducing environment. As shown in Fig. 4. In the oxidative condition, dissolved silicate adsorbed onto ferric hydroxide, and this matter is accumulated in sediment. When the lake-bottom environment was changed to an anoxic-reducing condition, hydrosulfide (: HS, and also hydrogen sulfide) was produced by a sulfate-reducing bacteria (SRB) and react with ferric hydroxide adsorbing silicate. As a result, ferric iron was reduced to the black iron sulfide (ferrous sulfide), at the same time, silicate was released in the water.

Fig. 4
figure 4

Behavior of dissolved silicate with a change in a lake-bottom environment

Adsorption and elution of silicate onto ferric hydroxide in laboratory experiments

In the silicate and dissolved iron mixed solution, the dissolved silicate concentration in filtrate A collected after the adsorption experiment noticeably decreased with an increase in the iron concentration. The adsorbed silicate fraction was nearly 80% at the iron concentration (3.6 mmol/L) (Fig. 5). This result indicates that silicate is efficiently adsorbed onto ferric hydroxide at neutral pH values.

Fig. 5
figure 5

Adsorption of silicate onto ferric hydroxide precipitate for the concentration of adsorbed silicate (filled circle) and percentage of silicate adsorbed (open circle). The initial silicate concentration was 0.16 mmol/L; pH of each solution was adjusted to 7.00 ± 0.05

The precipitate (ferric hydroxide) collected on the filter was placed into 10 mL of a 3.1 mmol/L sodium sulfide solution (pH = 7.0), whose color changed from brown to dark gray. This result suggests that ferric hydroxide was reduced to ferrous ion by the addition of sodium sulfide and that ferrous sulfide was formed by the reaction between ferrous ions and bisulfide ions remaining in the solution.

Filtrate B was collected from a mixed solution of precipitate and sodium sulfide after standing for 10 min, after which the Si concentration in the filtrate was measured. Figure 6 shows the concentrations of silicate in filtrates A and B. The sum of the Si concentrations in the two filtrates was equal to 94% and 93% of the initial Si concentration (0.16 mmol/L) during the adsorption experiments at iron concentrations of 0.54 and 1.80 mmol/L, respectively. The silicate concentrations in filtrate B, i.e., the concentration of Si eluted from the ferric hydroxide precipitate by adding the sodium sulfide solution, were 82% and 92% of the adsorbed silicate (the difference between the initial concentration of 0.16 mmol/L and the concentration in filtrate A), respectively. These results agree with our hypothesis that silicate is eluted from the ferric hydroxide during its reaction with sulfide ions as we already explained in the former section using Fig. 4.

Fig. 6
figure 6

Elution of the silicate adsorbed onto ferric hydroxide upon the addition of sodium sulfide. Results for filtrate A (concentration of silicate not adsorbed onto ferric hydroxide) shown in black; results for filtrate B (concentration of silicate eluted from ferric hydroxide by the addition of sodium sulfide) shown in white

Silicate elution from Lake Nakaumi bottom sediment based on laboratory experiments

We examined the elution of silicate from sediment core samples collected from Lake Nakaumi by adding a sodium sulfide solution. The concentrations of dissolved silicate eluted from sediment increased with increasing sodium sulfide concentrations (Fig. 7). We also detected the silicate concentration by adding a solution without sodium sulfide because the sediment originally contained sulfide ions. However, the addition of sodium sulfide (8.8 mmol/L) increased the silicate concentration twofold from its level at 0 mmol/L. These results suggest that silicate, which was eluted to pore water from ferric hydroxide, remained in the bottom sediments due to hydroxide reduction via sulfide ions.

Fig. 7
figure 7

The amount of silicate eluted when sodium sulfide was added to the anoxic sediments collected from Lake Nakaumi, based on sample depth of 0–5 cm and pH of 6.94. The result was obtained at 10 min after adding the sodium sulfide solution, respectively

Other sources of dissolved silicate in pore water

In addition to ferric hydroxide adsorbing a great deal of silicate, biogenic silica is also a possible source of dissolved silicate in pore water (Van der Weijden 2007). However, Rickert et al. (2002) reported that the solubility of biogenic silica was about 1.4 mmol/L at a water temperature of 14 °C (this temperature is the average sediment temperature of Lake Nakaumi in the summer), which is much lower than the maximum concentration observed in Lake Nakaumi. In addition, the sediments are characterized by a decline in biogenic silica dissolution rate, and differences in dissolution rate based on the presence or absence of oxygen are relatively small (Lehtimäki et al. 2016). Rickert et al. (2002) also reported that inorganic and organic materials coat the surface of biogenic silica in sediments. These coatings significantly reduce its dissolution rate, which indicates that the complete dissolution of biogenic silica in sediment to pore or bottom water requires approximately 100 years (Hurd 1973). Furthermore, Testa et al. (2013) reported that biogenic silicate was not decomposed in an anoxic-reducing sediment because of bacteria’s inability to decompose diatom cell walls. Thus, biogenic silica may not contribute to the highly concentrated dissolved silicate and its rapid dissolution during the short periods that are applicable to this study of Lake Nakaumi.