- 1Indiana University, O’Neill School of Public and Environmental Affairs, Bloomington, IN, United States
- 2St. Croix Watershed Research Station, Marine on St. Croix, MN, United States
- 3Washington State Department of Ecology, Lacey, WA, United States
River damming disturbs the transport and fate of nutrients, which alters biogeochemical cycling within reservoirs and alters the flux of nutrients downstream. In reservoirs, silicon (Si) is retained in larger proportions relative to nitrogen (N) and phosphorus (P) which can reduce dissolved Si (DSi) availability and push phytoplankton communities to be dominated by non-siliceous, potentially harmful taxa, such as cyanobacteria. Lake Monroe is the largest reservoir in Indiana and provides drinking water for more than 140,000 people, making it a critical water resource and exemplar for potential Si retention within the Mississippi River basin. To quantify the retention of DSi in Lake Monroe, we calculated an annual DSi budget using measured DSi inputs and outputs between April 2020 and March 2021. We also measured in-lake DSi and phytoplankton community composition between May and October 2020 as well as long-term biogenic Si fluxes to the sediment to better quantify the mechanisms controlling DSi retention. We found that Lake Monroe retained over half of its annual DSi inputs over the monitoring period and that retention was driven by diatom growth and sedimentation. As the construction of large dams continues globally, it is important to quantify how the biogeochemical cycling and transport of DSi is changing and the role of reservoirs in potentially shifting diatoms from N- or P- limited to Si-limited.
Introduction
The importance of silicon (Si) to inland and coastal ecosystems is well described, particularly in the context of eutrophication and harmful algal blooms (Officer and Ryther, 1980; Conley et al., 1993; Billen and Garnier, 2007). The Si depletion hypothesis posits that nutrient loading, particularly of phosphorus (P), can increase diatom production in lakes and the sequestration of biogenic Si in lake sediments (Schelske and Stoermer, 1971; Schelske et al., 1986). Over periods of years to decades, increased sequestration of biogenic Si in lake sediments disrupts the steady-state processes of dissolved Si (DSi) uptake and biogenic Si dissolution, thereby reducing DSi availability in the epilimnion of stratified lakes (Crowe, 2006; Carey and Fulweiler, 2013; Chen et al., 2014). Often, this results in nutrient conditions that limit diatom production and favor blooms of non-siliceous taxa, such as cyanobacteria. The process of long-term depletion of Si is distinct from the annual utilization of DSi by diatoms that temporarily reduces DSi concentrations (Conway et al., 1977; Schelske et al., 1986), though both phenomena result from the coupled cycling of nitrogen (N), P, and Si by phytoplankton (Hobbs et al., 2010; Scibona et al., 2022; Sethna and Royer, 2024).
Damming of rivers alters the quantity, timing, and stoichiometry of nutrient delivery to downstream water bodies which affects phytoplankton productivity and, subsequently, food web dynamics, carbon sequestration, and water quality (Smith, 2003; Paerl et al., 2006; Conley et al., 2009; Poff and Schmidt, 2016). Previous work has quantified global Si retention by river damming, which significantly alters the export of Si to downstream ecosystems and the stoichiometry of available nutrients (Beusen et al., 2009; Laruelle et al., 2009; Maavara et al., 2014). Additionally, reservoirs can alter the limiting nutrient in freshwater and marine systems, including potentially shifting diatoms from N- or P- limited to Si-limited (Paerl et al., 2006; Howarth and Marino, 2006), yet few studies have examined reservoir Si budgets and the role of Si stoichiometry in phytoplankton community succession.
The load of DSi transported by the Mississippi River to the northern Gulf of Mexico declined by approximately 50% during the middle of the 20th century (Turner et al., 1998; Turner and Rabalais, 2004). One proposed mechanism for the decline in DSi loads in the Mississippi River is sequestration of biogenic Si in nutrient-rich reservoirs on streams and rivers within the Mississippi River basin. Indeed, mass balance studies have shown reservoirs to be retention structures for Si (Maavara et al., 2014; Maavara et al., 2015a); however, there are very few published Si mass balance studies for lakes or reservoirs within the Mississippi River basin. Triplett et al. (2008) reported annual Si retention for two natural riverine lakes in Minnesota, Lakes Pepin and St. Croix, with higher retention occurring in eutrophic Lake Pepin. At present, the role of constructed reservoirs on DSi fate and transport in the Mississippi basin is largely unquantified.
Here we detail a 1-year DSi budget for Lake Monroe, a reservoir in southern Indiana and within the Mississippi River basin, using inputs from major tributaries and precipitation, outflow from the dam, drinking water withdrawal, in-lake changes in DSi concentrations and diatom abundance, and biogenic Si sequestration in reservoir sediments. Our study had three main objectives: (1) use a mass-balance approach to determine whether, and to what degree, Lake Monroe retained DSi, (2) explore the role of diatoms in driving DSi retention, and (3) calculate the relative retention of DSi and compare it with published values across a gradient of reservoir trophic status. We hypothesized Lake Monroe would show net retention of DSi during the sampling period; the main driver of this retention was expected to be uptake by diatoms, which can effectively sequester DSi through sedimentation (Turner et al., 1998; Frings et al., 2014; Wang et al., 2016; Liu et al., 2019). Diatoms are the most significant user of DSi in freshwater systems (Wetzel, 2001; Thamatrakoln and Kustka, 2009); therefore, the concentration of DSi in the epilimnion should correspond to the relative proportion of diatoms versus non-siliceous phytoplankton. If DSi concentrations decreased with increasing diatom abundance, this would indicate the growth of diatoms affected the availability of DSi for downstream transport. Finally, calculating the relative retention of DSi (RDSi, Maavara et al., 2014) in Lake Monroe quantified the proportion of DSi inputs retained in the reservoir and allowed for a comparison with other reservoirs and natural impoundments. As the rate of dam construction continues to increase globally (Zarfl et al., 2015), it is critical to quantify DSi retention in reservoirs, and this study provides estimates of retention within the Mississippi River basin where the load of DSi can strongly influence coastal food webs and ecological processes (Turner and Rabalais, 2004; Turner et al., 2007). For example, the timing of DSi delivery to the Gulf of Mexico can influence the growth of diatoms whose productivity provides the majority of the organic matter which drives the respiration that contributes to the formation of the hypoxic zone (Turner and Rabalais, 1994; Royer, 2020).
Methods
Site description
Lake Monroe is the largest lake in Indiana, with a surface area of approximately 44 km2 at normal pool and a maximum depth of 17 m (Jones et al., 1997). Globally, reservoirs ranging between 10 and 100 km2 are estimated to make up between 9%–12% of the total volume of impounded water (Downing et al., 2006; Lehner et al., 2011). The average volume and residence time calculated during the study period were 0.25 km3 and 186 days, respectively, based on the daily volume and outflow from the dam measured by the U.S. Army Corps of Engineers (ACOE). The lake is an impoundment on Salt Creek and was constructed in 1965 for the primary purpose of flood control and flow regulation during dry periods. The main tributaries to the lake are the North, Middle, and South Forks of the Salt Creek (Figure 1). The Lake Monroe watershed is approximately 1,095 km2 consisting of more than 86% forested land and about 12% agricultural land. Soils in the watershed are poorly to well drained and are all classified as “highly erodible lands” given the steep slopes within the watershed (Jones et al., 1997). The watershed is underlain primarily by limestone and siltstone of the Mississippian Borden Group (Jones et al., 1997). Today, Lake Monroe provides drinking water for more than 140,000 people, with withdrawals totaling 0.02 km3 between April 2020 and March 2021 (City of Bloomington Utilities). Lake Monroe is meso-eutrophic according to the Trophic State Index (Carlson, 1977), which was calculated using chlorophyll-a concentrations measured during our sampling period.

Figure 1. Map of the Lake Monroe Watershed. Lake sampling sites are numbered 1–5. Land cover data is a subset of the 2011 National Land Cover Database. Inset map shows the location of the Lake Monroe Watershed within the Mississippi River Basin.
Sampling regime
We monitored four tributaries and the dam outflow monthly for one calendar year between April 2020 and March 2021, collecting samples for DSi analysis and measuring instantaneous discharge. Each of the four tributaries drained sub-watersheds that, combined, account for about 55% of the Lake Monroe watershed area (Table 1; Figure 1). Stream water samples were filtered in the field using 0.45 µm cellulose filters (Fisherbrand), transported on ice, and frozen until analysis. Instantaneous discharge (Q) was measured using an electromagnetic water velocity probe (March-McBirney Model 2000 FloMate) as described in Hauer and Lamberti (2006). We also used instantaneous discharge measurements recorded by the U.S. Geological Survey at the North and South Forks of the Salt Creek (gage numbers 03371650 and 03371600, respectively). We sampled after one storm event in March 2021 in order to capture nutrient conditions and tributary discharge during high flows.

Table 1. Summary statistics for Lake Monroe tributary and outlet sampling sites, including mean DSi concentrations, instantaneous discharge (Q) descriptors, and instantaneous DSi loads during the sampling period (April 2020–March 2021). Each site was sampled monthly to construct a DSi budget for the watershed.
To characterize nutrient concentrations and phytoplankton community composition, we sampled five sites along a longitudinal transect through the lake monthly between late May and late October 2020 (Figure 1). From the three main sampling locations (Sites 1, 3, and 5), we collected epilimnetic samples using a 2-m integrated sampler and, during periods of stratification, hypolimnetic samples about 1 m above the lake bed using a Van Dorn sampler. Lake DSi samples were filtered and stored in the same way as tributary nutrient samples and additional, unfiltered samples were collected for total nitrogen (TN) and total phosphorus (TP) analysis. We used TN and TP to calculate stoichiometric ratios since bioavailable N and P in the water column are made up of particulate and organic forms which can rapidly remineralize (Dodds, 2003; Bergström, 2010; Carey et al., 2019). We also measured photosynthetically active radiation (Li-Cor Model LI-189 Quantum Radiometer-Photometer; Lincoln, NE) just below the water surface and collected temperature profiles (In-Situ Aquatroll 500 Multiparameter Sonde; Fort Collins, CO) at 1-m intervals to determine thermal stratification and metalimnion depth. At two midpoint sites between the main sites (Sites 2 and 4), we collected epilimnetic temperature and DSi samples. Lastly, we collected depth profiles of DSi concentration at 2-m intervals at Site 5, the deepest site located approximately 500 m upstream of the dam.
Phytoplankton samples were collected from the epilimnion of the three main sites, stored in dark plastic, and preserved with glutaraldehyde. Samples were imaged and identified to the genus level by PhycoTech, Inc. (St. Joseph, MI, United States) using the Imaging FlowCytobot (IFCB; McLane Research Laboratories, Inc., East Falmouth, MA, United States). The IFCB produced at least 2000 images per sample which were normalized to the volume of sample analyzed and identified to the genus level at 80% accuracy; however, classification at the functional group level is more accurate and provides an ecological assessment of the phytoplankton community. Biovolume of the phytoplankton was estimated for each individually imaged cell (Moberg and Sosik, 2012). In this analysis, we focus on the biovolume of diatoms (Bacillariophyta), which were dominated by the taxa Aulacoseira, Fragilaria, Nitzschia, and other less abundant centric and pennate diatoms.
All tributary and lake DSi samples were analyzed using the heteropoly blue method (Sultan, 2014) on a Lachat QuikChem flow injection analyzer (Model 8500; Hach Company; Loveland, CO). Samples were analyzed for TN and TP colorimetrically following an alkaline persulfate digestion (APHA, 2017) using an Alpkem FLOW Solution Autoanalyzer (Model 3570; OI Analytical; College Station, TX). For all nutrient analyses, we ensured the standard curves had coefficients of determination of at least 0.9, ran certified, external standards, and routinely calculated the method detection limits. Detection limits for DSi, TN, and TP were 0.06, 0.1, and 0.002 mg L−1, respectively. Samples below the detection limit included four TN samples from Crooked Creek and one TP sample from the epilimnion of Site 5. No DSi samples were below detection. Samples below detection were set to one half of the measured detection limit for purposes of statistical analysis.
We collected a sediment core from the deepest point in the reservoir, close to Site 5, in an area of the basin that was away from steep slopes that might be subject to slumping that was identified using sonar and bathymetric maps. The sediment core was collected using a piston corer with a 6.5-cm diameter polycarbonate core barrel and operated from the lake surface by metal drive rods (Wright, 1991). We recovered 23 cm of sediment which was transported vertically to the laboratory where the core was sectioned at 1-cm intervals and stored in polypropylene jars for the subsequent quantification of dry-density, water, organic, and mineral content, radiometric dating, and Si content.
The dry-density (dry mass per volume of fresh sediment), water content, and organic and mineral composition of each interval were determined using standard loss-on-ignition techniques (Dean, 1974). Lake diatoms were identified throughout the core, indicating the core represented sediment deposited after the reservoir was constructed. In order to accurately estimate the age of the sediment intervals, we used Gamma spectrometry to quantify the activity of lead (Pb) radioisotopes (214Pb and 210Pb) and construct a dating model for the core (Appleby, 2002). Briefly, we fit an exponential regression to the measured 210Pb activity to estimate the total inventory of unsupported 210Pb (Binford, 1990; Appleby, 2002). We used the activity of 214Pb as a measure of supported 210Pb; 214Pb is the daughter isotope of 226Ra which is in secular equilibrium with supported 210Pb (Appleby, 2002). We then used a constant rate of supply model to estimate the age and dry mass accumulation rates of each of the sediment core intervals (Krishnaswamy et al., 1971; Appleby, 2002; Supplementary Figure S1). Finally, we calculated a sediment focusing factor to correct for the spatial heterogeneity in sedimentation rates and to allow for a whole-basin estimate of Si burial (Hobbs et al., 2013; Engstrom and Rose, 2013).
DSi budget calculations
We quantified DSi inputs and outputs for Lake Monroe to construct an annual DSi mass balance from April 2020 through March 2021. Tributary DSi loads and atmospheric wet deposition were considered watershed inputs while the dam outflow and drinking water withdrawal were considered outputs. Total DSi retained in Lake Monroe during the sampling period was therefore calculated as the difference between total watershed inputs and outputs. The relative retention of DSi (RDSi) was calculated as:
where DSiin and DSiout are the total fluxes of DSi in and out of Lake Monroe, respectively (Maavara et al., 2014). We then compared the RDSi of Lake Monroe to the global dataset of reservoirs compiled by Maavara et al. (2014).
Daily loads in the tributaries and outflows were modeled using Loadflex, a package in R that allows for the simple, linear interpolation of solute loads between sampling events (Appling et al., 2015). Loads were interpolated between sampling events since DSi concentrations exhibited a chemostatic relationship with discharge at all sampling sites (Supplementary Figure S2). We used U.S. Geological Survey (USGS) continuous, daily discharge data, monthly DSi concentrations, and a simple linear interpolation model to estimate total daily DSi loads from the North and South Forks of the Salt Creek (Table 1). We used these daily loads to calculate daily DSi yields which were scaled to the outlet of the North and South Fork tributaries. We used the average of the North and South Fork yields to scale total loads from other sampled watersheds and the unmonitored portion of the Lake Monroe watershed. Similarities in watershed land use and DSi concentrations between sub-watersheds and sampled tributaries allowed us to reasonably assume proportional contributions of DSi from unmonitored portions of the watershed (David et al., 2006). Due to the inherent uncertainty associated with stream discharge measurements, we applied 10% error to each value measured at the U.S. Geologic Survey gaging stations (Harmel et al., 2006) and calculated the DSi loads using ±10% of the reported discharge values. Atmospheric deposition of DSi was calculated by multiplying the annual volume of precipitation over the lake by the average DSi concentration in precipitation (both rain and snow) collected during three distinct events during the monitoring period. Outflow DSi loads were modeled using the same Loadflex methods, but discharge data were provided by the U.S. ACOE. The City of Bloomington drinking water plant is located near Site 2 and monitors daily withdrawals. We estimated the total mass of DSi removed with drinking water using the average DSi concentration at Site 2 measured between April and October of 2020 multiplied by the total volume of water withdrawn.
We estimated the mass of DSi retained in the hypolimnion by multiplying the measured hypolimnetic DSi concentration by the volume of the hypolimnion before lake turnover, estimated based on the depth of the metalimnion and area of the lake. To quantify the mass of Si sequestration in reservoir sediments, we analyzed biogenic silica (BSi) concentrations using 30 mg subsamples of freeze-dried sediment representing the last 10 years. Samples were digested in 40 mL of 1% Na2CO3 solution heated at 85 °C in a reciprocating water bath for 5 hours (Conley and Schelske, 2002). After cooling and neutralization with 4.5 g of 0.021 N HCl solution, DSi was measured colorimetrically on a Unity Scientific SmartChem 170 discrete analyzer (SmartChem 2012a). The concentration of DSi in each sediment interval was converted to the flux of BSi over time using the dry mass accumulation rates.
Statistical analysis
All analyses and statistical tests were carried out in R (The R Foundation for Statistical Computer, Version 4.0.5, 2021) and data were visualized using the ggplot2 package (version 3.5.1; Wickham, 2016). We hypothesized that uptake by diatoms was the main driver of DSi retention in the lake; therefore, we expected a declining trend in DSi concentrations along the downstream gradient of the reservoir, higher DSi concentrations in the hypolimnion relative to the epilimnion, and a relationship between DSi concentrations and the relative abundance of diatoms that reflected uptake and sedimentation. We used simple linear regression to model the spatial trend in DSi concentrations as well as a one-way analysis of variance (ANOVA) to test for differences in DSi concentrations between lake sampling sites for each sampling date. We also used t-tests to detect significant differences in DSi concentrations between the epilimnion and hypolimnion. The relationship between DSi concentrations and the abundance of diatoms relative to the total phytoplankton community was modeled using simple linear regression. Prior to running each statistical analysis, we tested all data to ensure the assumptions for each test were met, including normality using the Shapiro-Wilk test and homoskedasticity by plotting model residuals against fitted results.
Results
Annual DSi budget
The primary source of DSi to Lake Monroe was tributary inputs, which varied through time and between sites. Between all sites and sampling dates, DSi concentrations in the tributaries ranged between 3.6 and 13.3 mg L−1 while discharge ranged over multiple orders of magnitude (Table 1). Due to the inherent uncertainty associated with stream discharge measurements, we applied 10% error to each value measured at the U.S. Geologic Survey gaging stations (Harmel et al., 2006). Annual input from the watershed was estimated to be 5,712 (±571) Mg SiO2. The total annual DSi input from precipitation onto the lake was approximately 7 Mg SiO2, which represents ∼0.1% of the total DSi inputs to the lake during the study period (Figure 2).

Figure 2. DSi budget for Lake Monroe for the period April 2020 through March 2021. Inputs include DSi loads from monitored and unmonitored tributaries as well as precipitation over the lake. The outputs from the watershed include the dam outflow and drinking water withdrawals. The ranges in tributary inputs and dam outflow represent the uncertainty applied to discharge measurements. The difference between the inputs and outputs, between 2,982 and 3,688 Mg SiO2, represents the net retention of DSi in the reservoir during the study period, of which 854 Mg are estimated to be buried annually in the sediments as biogenic Si.
Outflows from Lake Monroe include the dam outlet and drinking water withdrawals. DSi concentrations from the outlet fluctuated between 8 and 10 mg L−1 between April and October 2020 before declining to 3–5 mg L−1 through the winter and early spring of 2021. The total output of DSi from the dam was 2,181 (±218) Mg SiO2 over the period of record. The drinking water withdrawals totaled 0.02 km3 from April 2020 to May 2021 and we calculated the total mass of DSi removed with drinking water to be 203 Mg SiO2. There was intra-annual variation in DSi retention within the lake; loads in the outflow exceeded the total incoming load for most months in 2020 but were much less than the total input during the winter and early spring of 2021 (Figure 3). Based on total average input and output estimates, the annual retention of DSi was 3,335 (±335) Mg SiO2 and the annual RDSi was 0.58 (±0.08), meaning that during the study period, DSi outflow was equivalent to less than half of the DSi inputs. Even with the most conservative estimates of DSi retention (minimum inputs, maximum outputs), RDSi was 0.49, showing meaningful DSi retention in Lake Monroe.

Figure 3. Total monthly DSi loads (as megagrams [Mg] SiO2) from each sub-watershed (stacked bars) and the dam outflow (gray bars).
The dating model for the sediment core produced a continuous age-depth relationship for the full length of the sediment core. The base of the core was dated to approximately 2007 and the top 14 cm of core were found to represent the last decade of sedimentation in Lake Monroe. Focus-corrected, areal sediment BSi fluxes averaged 1.9 (±0.1) mg Si cm−2 year−1 and the consistent flux values within the last decade indicate Si sequestration has remained consistent over this time period (Figure 4). We calculated the annual flux of Si to the sediments within the entire Monroe basin by multiplying the average areal flux by the surface area of the lake (44 km2), yielding a whole basin burial of 854 Mg of Si annually and indicating that nearly 30% of the DSi retained in the reservoir is permanently buried as BSi in the lake sediments.

Figure 4. The annual, focus-corrected, areal flux of BSi measured over the last decade in Lake Monroe.
Variation in DSi concentrations in Lake Monroe
Contrary to our predictions, epilimnetic DSi concentrations did not significantly decline along the longitudinal gradient of the lake. In fact, epilimnetic DSi concentrations were not significantly different among any of the epilimnion sampling sites during the monitoring period. However, there were differences between epilimnion and hypolimnion concentrations that indicated changes in DSi storage through time within the lake (Figure 5). Specifically, in May, July, August, and September, DSi concentrations in the hypolimnion were significantly greater than epilimnion DSi (t-test, p < 0.01; Figure 6). Lake Monroe was stratified between May and September, limiting mixing between the epilimnion and hypolimnion; therefore, changes in DSi concentrations were likely a result of diatom uptake in the epilimnion and dissolution of biogenic Si in the hypolimnion. Based on the hypolimnetic DSi concentration measured at the Lower site (13.4 mg L−1) and the depth of the metalimnion (10.5 m) in late September, we estimated 2,653 Mg of DSi were recycled from the hypolimnion to the epilimnion during fall turnover, representing the other 70% of the DSi retained in the reservoir.

Figure 5. Isopleth diagram of DSi concentrations (mg SiO2 L−1) from Site 5 during the monitoring period. Samples were collected at 2-m intervals through the water column. Isopleths are spaced at 0.5 mg L−1 intervals.

Figure 6. Epilimnetic (light blue) and hypolimnetic (dark blue) DSi concentrations during the sampling period across all sampled sites. Epilimnetic concentrations were measured using a 2 m integrated sampling tube at the surface of the lake, and hypolimnetic samples were collected 1 m from the lake bed using a Van Dorn sampler. Data are grouped by depth and sampling month. Each box represents the 25th to 75th percentiles, the horizontal line indicates the median value, whiskers extend to ±1.5 Interquartile Range, and points indicate outliers in each group.
As we expected, DSi concentrations had a significant, negative relationship (simple linear regression, p < 0.001 for all sites) with diatom biovolume at sites 1 and 3 throughout the monitoring period (Supplementary Figure S3). We modeled the relationship between each sampling site individually to account for the differences between site behavior along the downstream gradient of the reservoir. DSi concentrations at sites 1 and 3 explained 84% and 66% of the variance in relative diatom biovolume, respectively, while there was no significant relationship between diatom biovolume and DSi concentrations at site 5. This suggests diatoms, as the most significant user of DSi, reduced epilimnetic DSi concentrations differently along the downstream gradient within the reservoir. The relationship between diatom biovolume and DSi concentrations is especially pronounced when comparing changes in diatom biovolume simultaneously with DSi concentrations through time (Figure 7). When diatoms peaked in July and October, there was a corresponding decline in DSi concentrations. While there was a significant relationship between diatom biovolume and DSi concentrations, there was no effect on diatom biovolume with changes in molar Si:TN, Si:TP, or TN:TP (Supplementary Figure S4).

Figure 7. The biovolume of diatoms, averaged across Sites 1, 3, and 5 (top) and the mean epilimnetic DSi concentration (bottom) during the sampling period in 2020. Error bars extend to the standard error of the mean and in some cases are too small to be visible.
Discussion
Understanding DSi cycling in inland waters is critical to resolving how modified landscapes influence nutrient delivery along the land-ocean continuum. Reservoirs represent a major gap in our knowledge of Si retention and transformation, despite their growing prevalence and known impacts on biogeochemical cycles (Maavara et al., 2020a; Maavara et al., 2020b). To our knowledge, our DSi budget in Lake Monroe is the first to quantify DSi retention within the Mississippi River basin, which contributes increasingly imbalanced loads of N, P, and Si to the Gulf of Mexico that can favor non-siliceous and potentially harmful algal taxa (Turner and Rabalais, 2004; Royer, 2020). We found that Lake Monroe retained approximately 58% of its annual DSi inputs, providing clear evidence that reservoirs with relatively short residence times can act as substantial Si sinks. Diatom uptake and sedimentation were key drivers of DSi retention, indicated by a strong, inverse relationship between DSi concentrations and diatom biovolume as well as the permanent burial of BSi in lake sediments. When compared to published estimates from reservoirs with similar water residence times, the DSi retention in Lake Monroe was higher than expected, underscoring the largely unquantified role of reservoirs in global Si retention. Taken together, our results highlight the important role of reservoirs in modulating Si transport to downstream ecosystems and coastal waters.
Lake Monroe was a net sink of DSi
Based on our measurements over the course of 1 year, Lake Monroe retained about 58% of the total DSi inputs. Compared with other reservoirs around the world, the RDSi of Lake Monroe is in the top 25th percentile but is similar to the median value of reservoirs with similar residence times (Harrison et al., 2012; Maavara et al., 2014). Globally, the average RDSi is 0.13, with values above 0.3 generally associated with water residence times greater than 5 years (Maavara et al., 2014; Maavara et al., 2015a). The average water residence time of Lake Monroe during the sampling period was about 0.5 years, indicating the DSi retention observed in Lake Monroe during our monitoring period was higher than expected based on the relationship between water residence times and DSi retention described by Maavara et al. (2014); however, this is not unexpected since the relationship in Maavara et al. (2014) was established using budgets for just 20 reservoirs from a global population of more than 75,000 reservoirs with a surface area >0.1 km2 (Lehner et al., 2011). The water retention time in Lake Monroe during the sample period was similar to the median value of historic water volume and dam outflow, indicating the DSi retention we measured is likely representative of annual Si retention during the last decade (Supplementary Table S1). Furthermore, sediment Si fluxes indicated relatively constant Si depositional rates, suggesting Si retention rates have not changed in the last 10 years and providing no evidence for Si depletion in Lake Monroe. The significant DSi retention observed in Lake Monroe provides further evidence of the impact of reservoirs on nutrient delivery to downstream waters and, ultimately, coastal ecosystems. In the midwestern U.S., occurrences of heavy precipitation are predicted to increase with climate change (Hamlet et al., 2020), which could increase inputs of DSi to the Mississippi River. However, high flow events will likely disproportionately increase N loads, reducing the availability of DSi relative to N in receiving waters (Leong et al., 2014; Sinha et al., 2017; Royer, 2020; Sethna et al., 2022).
Retention of DSi in reservoirs, as observed in Lake Monroe, has significant implications for the transport of DSi to the Mississippi River and, ultimately, the Gulf of Mexico. Two natural impoundments in the Mississippi River basin in Minnesota, Lakes Saint Croix and Pepin, had RDSi values of 0.04 and −0.11, respectively, with the lower of the two corresponding to Lake Pepin, which had an order of magnitude greater Si inputs and was also more eutrophic relative to Lake Saint Croix (Triplett et al., 2008; Maavara et al., 2014). Lake Monroe and the lakes studied by Triplett et al. (2008) all are in the upper portions of the Mississippi River basin and DSi released from these lakes is subject to retention and sedimentation before reaching the Gulf of Mexico. Interestingly, Downing et al. (2016) concluded that impoundments within several agricultural watersheds of the Mississippi River basin did not increase Si retention, rather, this was attributed to the low abundance of diatoms relative to cyanobacteria in those reservoirs. Together, Lakes Monroe, Saint Croix, and Pepin represent <0.001% of the Mississippi River basin area and <5% of the total Si load to the Gulf of Mexico, highlighting the significant uncertainty that exists in scaling up our measurements of Si retention in small reservoirs to Si retention at a continental scale.
Continued study of Si transport and availability in rivers and reservoirs throughout the basin are necessary to identify nutrient management practices that alleviate cultural eutrophication and hypoxia in the Gulf of Mexico. These types of studies are also critical in other large rivers such as the Yangtze River in China and the Amazon River in Brazil, where many large dams are currently under construction, or planned for in the near future (Yang et al., 2011; Zarfl et al., 2015; Anderson et al., 2018; Flecker et al., 2022). Increasing Si retention in these basins, coupled with increased anthropogenic loads of N and P, will likely contribute to more frequent and intense Si limitation in coastal systems, with implications for phytoplankton productivity, water quality, and carbon cycling (Garnier et al., 2010; Maavara et al., 2020a). These issues are compounded by dam-driven changes to the timing of nutrient loading, which could increase the potential for Si limitation during the growing season and contribute to increased blooms of harmful or nuisance algal taxa.
Contrary to our predictions, epilimnetic DSi concentrations did not decrease along the lake continuum suggesting watershed inputs maintained relatively constant Si concentrations. There were, however, significant differences in DSi concentrations between the epilimnion and hypolimnion, indicating that the probable mechanism facilitating retention is the growth and sedimentation of diatoms. Diatom growth reduces DSi concentrations during blooms and frustules sink after death or consumption, effectively sequestering the biogenic Si in lake sediments. As biogenic Si is more soluble than crystalline Si (Cornelis et al., 2011), some diatom frustules are recycled to DSi, thereby increasing the DSi concentration in the hypolimnion (Scibona et al., 2022). The storage and dissolution of diatom frustules in the hypolimnion explains why hypolimnetic DSi concentrations were higher than epilimnetic DSi.
The sequestration and dissolution of diatoms in lake sediments is the most likely fate of the retained DSi in Lake Monroe. Increased hypolimnetic DSi concentrations support the hypothesis that some of the diatom frustules dissolve in the hypolimnion and are recycled to the epilimnion at turnover (Figure 4); however, the flux of BSi to the sediments indicates that nearly a third of the DSi retained in the reservoir is buried in lake sediments (Triplett et al., 2008; Humborg et al., 2008; Maavara et al., 2015a). The other 70% of the Si was recycled back into the water column during fall turnover in late September (Figure 5). The recycling of DSi from the hypolimnion elevated epilimnetic DSi concentrations and likely facilitated the observed diatom bloom in late October (Figure 7). Further examination of the biogeochemical cycling of DSi between the water column and lake sediments could help us better describe the processes that facilitate BSi sequestration and recycling at lake turnover.
Silica utilization in Lake Monroe
The inverse relationship between the biovolume of diatoms and DSi concentrations indicated that DSi availability was influenced by diatom productivity. Our data suggest the elevated DSi concentrations created conditions suitable for diatom growth and their growth depleted the available DSi in the lake, thus creating conditions more favorable for the growth of non-siliceous phytoplankton groups (Beusen et al., 2009; Carey et al., 2019). The relationship between DSi concentrations and diatom biovolume likely differed between sites because Site 1, as the most upstream site, was most strongly influenced by riverine inputs, while conditions at Sites 3 and 5 were influenced by in-lake processes. Interestingly, the minimum DSi concentration in Lake Monroe was 5.09 mg SiO2 L−1, well above the limiting concentration of 0.1 mg SiO2 L−1 reported by Schelske and Stoermer (1971). In fact, Si:TN and Si:TP molar ratios indicated abundant DSi availability relative to both N and P even though non-siliceous phytoplankton taxa dominated the community throughout most of the monitoring period. Our monthly sampling frequency was likely not of fine enough resolution to capture the minimum DSi concentration at the peak of the diatom bloom; however, there was a clear decline in diatom biovolume after DSi concentrations decreased, and a large increase in diatom productivity once DSi concentrations rose above 8 mg SiO2 L−1 (Figure 7). This study is a clear example of annual Si utilization by diatoms, but quantitative studies of nutrient limitation are necessary to determine how N and P availability affect Si uptake and patterns in phytoplankton succession. As we only collected DSi data over the course of one calendar year, we do not have sufficient evidence to determine whether the occasional reduction in DSi concentrations due to diatom growth will cause long-term DSi depletion in Lake Monroe.
Concentrations of DSi explained much of the variation in diatom biovolume; however, nutrient stoichiometry (Supplementary Figure S4), temperature, or photosynthetically active radiation had little explanatory power (Supplementary Figure S5). This suggests DSi concentrations had the strongest control on diatom growth, despite previous work connecting decreasing Si:N and Si:P ratios with increasing occurrences of eutrophication and harmful algal blooms (Officer and Ryther, 1980; Justić et al., 1995; Turner et al., 1998; Teubner and Dokulil, 2002). Reservoirs have been shown to preferentially retain TN and TP over Si thereby increasing the relative Si availability to downstream waters (Maavara et al., 2020b; Maavara et al., 2020a). In other words, while reservoirs can induce DSi depletion and subsequent reductions in riverine DSi flux, in some cases they reduce TN and TP to an even greater extent, thus increasing Si:N and Si:P ratios to those favoring siliceous algal taxa in downstream waters (Humborg et al., 2000; Humborg et al., 2008). Increased DSi loads to the Gulf of Mexico could exacerbate hypoxia as diatom productivity provides much of the organic matter fueling hypoxic zone formation (Turner and Rabalais, 1994; Dortch et al., 2001). However, in human-dominated regions, reservoirs cannot buffer riverine stoichiometry against the ever increasing anthropogenic inputs of N and P relative to Si, which results in a declining trend of Si:N and Si:P in coastal receiving waters (Beusen et al., 2009; Laruelle et al., 2009; Maavara et al., 2014).
Implications and direction for future work
Reservoirs play a vital role in human civilization through flood protection, maintenance of water supplies, support for transportation and navigation, and as a source of renewable power generation (Lehner et al., 2011). At the same time, reservoirs have negative environmental impacts including declines in biodiversity, contributing significant quantities of greenhouse gas emissions to the atmosphere, and altering biogeochemical processing along the river continuum (Deemer et al., 2016; Grill et al., 2019; Maavara et al., 2020b). Empirical studies of N, P, and Si retention in reservoirs are lacking in published literature, limiting our ability to accurately model changes in nutrient delivery to coastal systems (Maavara et al., 2014; 2015b; Akbarzadeh et al., 2019). As dam construction continues to impede free-flowing rivers, an improved understanding of biogeochemical nutrient cycling is necessary for accurate environmental assessments and water quality projections (Friedl and Wüest, 2002; Maranger et al., 2018).
Our study is the first report of DSi dynamics in Lake Monroe and provides further evidence of DSi retention in reservoirs. Recent reviews of the global retention of DSi through reservoirs highlight the significant fraction of DSi sequestered in lentic systems which represents an important component of the global biogeochemical cycle of Si (Harrison et al., 2012; Frings et al., 2014; Maavara et al., 2014). Our results align with previous research that shows reservoirs can be significant sinks of DSi, with diatom-mediated uptake and sediment burial acting as dominant retention mechanisms, underscoring the importance of including reservoirs, across a wide range in water residence time, in regional and global Si budgets. Longer term analyses of DSi retention within Lake Monroe and other similar reservoirs would allow for better characterization of DSi dynamics under a range of hydrologic conditions and an examination of trends in retention through time, including explicit testing of the silica depletion hypothesis.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
LS: Conceptualization, Methodology, Writing – original draft, Investigation, Visualization, Supervision, Writing – review and editing, Formal Analysis. LR: Investigation, Writing – review and editing. SP: Investigation, Writing – review and editing. ZP: Formal Analysis, Writing – review and editing. AP: Writing – review and editing, Formal Analysis. WH: Supervision, Writing – review and editing, Conceptualization. TR: Resources, Writing – original draft, Supervision, Writing – review and editing, Conceptualization.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported in part by the Indiana Department of Environmental Management through a Section 319 Nonpoint Source Management Grant, the Monroe County Stormwater Management Board, and the City of Bloomington Utility Service Board.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvs.2025.1656708/full#supplementary-material
References
Akbarzadeh, Z., Maavara, T., Slowinski, S., and Van Cappellen, P. (2019). Effects of damming on river nitrogen fluxes: a global analysis. Glob. Biogeochem. Cycles 33, 1339–1357. doi:10.1029/2019gb006222
Anderson, E. P., Jenkins, C. N., Heilpern, S., Maldonado-Ocampo, J. A., Carvajal-Vallejos, F. M., Encalada, A. C., et al. (2018). Fragmentation of Andes-to-Amazon connectivity by hydropower dams. Sci. Adv. 4, eaao1642. doi:10.1126/sciadv.aao1642
APHA (2017). Standard method for the examination of water and wastewater. Washington, D.C., USA: American Public Health Association.
Appleby, P. G. (2002). “Chronostratigraphic techniques in recent sediments,” in Tracking environmental change using lake sediments. Editors W. M. Last, and J. P. Smol (Dordrecht: Kluwer Academic Publishers), 171–203.
Appling, A. P., Leon, M. P., and McDowell, W. H. (2015). Reducing bias and quantifying uncertainty in watershed flux estimates: the R package loadflex. Ecosphere 6, 1–25. doi:10.1890/es14-00517.1
Bergström, A.-K. (2010). The use of TN:TP and DIN:TP ratios as indicators for phytoplankton nutrient limitation in oligotrophic lakes affected by N deposition. Aquat. Sci. 72, 277–281. doi:10.1007/s00027-010-0132-0
Beusen, A. H. W., Bouwman, A. F., Dürr, H. H., Dekkers, A. L. M., and Hartmann, J. (2009). Global patterns of dissolved silica export to the coastal zone: results from a spatially explicit global model. Glob. Biogeochem. Cycles 23, 1–13. doi:10.1029/2008GB003281
Billen, G., and Garnier, J. (2007). River basin nutrient delivery to the coastal sea: assessing its potential to sustain new production of non-siliceous algae. Mar. Chem. 106, 148–160. doi:10.1016/j.marchem.2006.12.017
Binford, M. W. (1990). Calculation and uncertainty analysis of 210Pb dates for PIRLA project lake sediment cores. J. Paleolimnol. 3. doi:10.1007/bf00219461
Carey, J. C., and Fulweiler, R. W. (2013). Nitrogen enrichment increases net silica accumulation in a temperate salt marsh. Limnol. Oceanogr. 58, 99–111. doi:10.4319/lo.2013.58.1.0099
Carey, J. C., Jankowski, K., Julian, P., Sethna, L. R., Thomas, P. K., and Rohweder, J. (2019). Exploring silica stoichiometry on a large floodplain riverscape. Front. Ecol. Evol. 7, 346. doi:10.3389/fevo.2019.00346
Carlson, R. E. (1977). A trophic state index for lakes. Limnol. Oceanogr. 22, 361–369. doi:10.4319/lo.1977.22.2.0361
Chen, N., Wu, Y., Wu, J., Yan, X., and Hong, H. (2014). Natural and human influences on dissolved silica export from watershed to coast in Southeast China. J. Geophys. Res. Biogeosciences 119, 95–109. doi:10.1002/2013jg002429
Conley, D. J., and Schelske, C. L. (2002). “Biogenic silica,” in Tracking environmental change using lake sediments. Editors J. P. Smol, H. J. B. Birks, W. M. Last, R. S. Bradley, and K. Alverson (Netherlands, Dordrecht: Springer), 281–293.
Conley, D., Schelske, C., and Stoermer, E. (1993). Modification of the biogeochemical cycle of silica with eutrophication. Mar. Ecol. Prog. Ser. 101, 179–192. doi:10.3354/meps101179
Conley, D. J., Paerl, H. W., Howarth, R. W., Boesch, D. F., Seitzinger, S. P., Havens, K. E., et al. (2009). Controlling eutrophication: nitrogen and phosphorus. Science 323, 1014–1015. doi:10.1126/science.1167755
Conway, H. L., Parker, J. I., Yaguchi, E. M., and Mellinger, D. L. (1977). Biological utilization and regeneration of silicon in Lake Michigan. J. Fish. Res. Board Can. 34, 537–544. doi:10.1139/f77-085
Cornelis, J.-T., Delvaux, B., Georg, R. B., Lucas, Y., Ranger, J., and Opfergelt, S. (2011). Tracing the origin of dissolved silicon transferred from various soil-plant systems towards rivers: a review. Biogeosciences 8, 89–112. doi:10.5194/bg-8-89-2011
Crowe, A. M. (2006). The role of silica depletion in the eutrophication of Lac la Biche. Alberta: University of Alberta.
David, M. B., Wall, L. G., Royer, T. V., and Tank, J. L. (2006). Denitrification and the nitrogen budget of a reservoir in an agricultural landscape. Ecol. Appl. 16, 2177–2190. doi:10.1890/1051-0761(2006)016[2177:datnbo]2.0.co;2
Dean, W. E. (1974). Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: comparison with other methods. SEPM J. Sediment. Res. 44. doi:10.1306/74D729D2-2B21-11D7-8648000102C1865D
Deemer, B. R., Harrison, J. A., Li, S., Beaulieu, J. J., DelSontro, T., Barros, N., et al. (2016). Greenhouse gas emissions from reservoir water surfaces: a new global synthesis. BioScience 66, 949–964. doi:10.1093/biosci/biw117
Dodds, W. K. (2003). Misuse of inorganic N and soluble reactive P concentrations to indicate nutrient status of surface waters. J. North Am. Benthol. Soc. 22, 171–181. doi:10.2307/1467990
Dortch, Q., Rabalais, N. N., Turner, R. E., and Qureshi, N. A. (2001). Impacts of changing Si/N ratios and phytoplankton species composition. Coastal hypoxia: consequences for living resources and ecosystems. Am. Geophys. Union, 37–48.
Downing, J. A., Prairie, Y. T., Cole, J. J., Duarte, C. M., Tranvik, L. J., Striegl, R. G., et al. (2006). The global abundance and size distribution of lakes, ponds, and impoundments. Limnol. Oceanogr. 51, 2388–2397. doi:10.4319/lo.2006.51.5.2388
Downing, J. A., Cherrier, C. T., and Fulweiler, R. W. (2016). Low ratios of silica to dissolved nitrogen supplied to rivers arise from agriculture not reservoirs. Ecol. Lett. 19, 1414–1418. doi:10.1111/ele.12689
Engstrom, D. R., and Rose, N. L. (2013). A whole-basin, mass-balance approach to paleolimnology. J. Paleolimnol. 49, 333–347. doi:10.1007/s10933-012-9675-5
Flecker, A. S., Shi, Q., Almeida, R. M., Angarita, H., Gomes-Selman, J. M., García-Villacorta, R., et al. (2022). Reducing adverse impacts of Amazon hydropower expansion. Science 375, 753–760. doi:10.1126/science.abj4017
Friedl, G., and A. Wüest, (2002). Disrupting biogeochemical cycles – Consequences of damming. Aquatic Sciences 64:55–65. doi:10.1007/s00027-002-8054-0
Frings, P. J., Clymans, W., Jeppesen, E., Lauridsen, T. L., Struyf, E., and Conley, D. J. (2014). Lack of steady-state in the global biogeochemical Si cycle: emerging evidence from lake Si sequestration. Biogeochemistry 117, 255–277. doi:10.1007/s10533-013-9944-z
Garnier, J., Beusen, A., Thieu, V., Billen, G., and Bouwman, L. (2010). N:P:Si nutrient export ratios and ecological consequences in coastal seas evaluated by the ICEP approach. Glob. Biogeochem. Cycles 24, 2009GB003583–12. doi:10.1029/2009gb003583
Grill, G., Lehner, B., Thieme, M., Geenen, B., Tickner, D., Antonelli, F., et al. (2019). Mapping the world’s free-flowing rivers. Nature 569, 215–221. doi:10.1038/s41586-019-1111-9
Hamlet, A. F., Byun, K., Robeson, S. M., Widhalm, M., and Baldwin, M. (2020). Impacts of climate change on the state of Indiana: ensemble future projections based on statistical downscaling. Clim. Change 163, 1881–1895. doi:10.1007/s10584-018-2309-9
Harmel, R. D., Cooper, R. J., Slade, R. M., Haney, R. L., and Arnold, J. G. (2006). Cumulative uncertainty in measured streamflow and water quality data for small watersheds. Trans. ASABE 49, 689–701. doi:10.13031/2013.20488
Harrison, J. A., Frings, P. J., Beusen, A. H. W., Conley, D. J., and McCrackin, M. L. (2012). Global importance, patterns, and controls of dissolved silica retention in lakes and reservoirs. Glob. Biogeochem. Cycles 26. doi:10.1029/2011gb004228
F. R. Hauer, and G. A. Lamberti (2006). Methods in stream ecology. 2nd ed. (San Diego, California: Academic Press/Elsevier).
Hobbs, W. O., Lalonde, S. V., Vinebrooke, R. D., Konhauser, K. O., Paul Weidman, R., Graham, M. D., et al. (2010). Algal-silica cycling and pigment diagenesis in recent alpine lake sediments: mechanisms and paleoecological implications. J. Paleolimnol. 44, 613–628. doi:10.1007/s10933-010-9441-5
Hobbs, W. O., Engstrom, D. R., Scottler, S. P., Zimmer, K. D., and Cotner, J. B. (2013). Estimating modern carbon burial rates in lakes using a single sediment sample. Limnol. Oceanogr. Methods 11, 316–326. doi:10.4319/lom.2013.11.316
Howarth, R. W., and Marino, R. (2006). Nitrogen as the limiting nutrient for eutrophication in coastal marine ecosystems: evolving views over three decades. Limnol. Oceanogr. 51, 364–376. doi:10.4319/lo.2006.51.1_part_2.0364
Humborg, C., Conley, D. J., Rahm, L., Wulff, F., Cociasu, A., and Ittekkot, V. (2000). Silicon retention in river basins: far-reaching effects on biogeochemistry and aquatic food webs in coastal marine environments. AMBIO A J. Hum. Environ. 29, 45–50. doi:10.1639/0044-7447(2000)029[0045:srirbf]2.0.co;2
Humborg, C., Smedberg, E., Medina, M. R., and Mörth, C.-M. (2008). Changes in dissolved silicate loads to the Baltic Sea—the effects of lakes and reservoirs. J. Mar. Syst. 73, 223–235. doi:10.1016/j.jmarsys.2007.10.014
Jones, W. W., Jenson, M., Jourdain, E., Mitchell-Bruker, S., Strong, L., Bieberich, L., et al. (1997). Lake Monroe diagnostic and feasibility study. Bloomington: Indiana University, 352.
Justić, D., Rabalais, N. N., Turner, R. E., and Dortch, Q. (1995). Changes in nutrient structure of river-dominated coastal waters: stoichiometric nutrient balance and its consequences. Estuar. Coast. Shelf Sci. 40, 339–356. doi:10.1016/s0272-7714(05)80014-9
Krishnaswamy, S., Lal, D., Martin, J. M., and Meybeck, M. (1971). Geochronology of lake sediments. Earth Planet. Sci. Lett. 11, 407–414. doi:10.1016/0012-821x(71)90202-0
Laruelle, G. G., Roubeix, V., Sferratore, A., Brodherr, B., Ciuffa, D., Conley, D. J., et al. (2009). Anthropogenic perturbations of the silicon cycle at the global scale: key role of the land-ocean transition. Glob. Biogeochem. Cycles 23. doi:10.1029/2008gb003267
Lehner, B., Liermann, C. R., Revenga, C., Vörösmarty, C., Fekete, B., Crouzet, P., et al. (2011). High-resolution mapping of the world’s reservoirs and dams for sustainable river-flow management. Front. Ecol. Environ. 9, 494–502. doi:10.1890/100125
Leong, D. N. S., Donner, S. D., Hassan, M. A., Gabor, R., and Drummond, J. D. (2014). Sensitivity of stoichiometric ratios in the Mississippi River to hydrologic variability. J. Geophys. Res. Biogeosciences 119, 1049–1062. doi:10.1002/2013jg002585
Liu, D., Yuan, P., Tian, Q., Liu, H., Deng, L., Song, Y., et al. (2019). Lake sedimentary biogenic silica from diatoms constitutes a significant global sink for aluminium. Nat. Commun. 10, 4829. doi:10.1038/s41467-019-12828-9
Maavara, T., Dürr, H. H., and Van Cappellen, P. (2014). Worldwide retention of nutrient silicon by river damming: from sparse data set to global estimate. Glob. Biogeochem. Cycles 28, 842–855. doi:10.1002/2014gb004875
Maavara, T., Hood, J. L. A., North, R. L., Doig, L. E., Parsons, C. T., Johansson, J., et al. (2015a). Reactive silicon dynamics in a large prairie reservoir (Lake Diefenbaker, Saskatchewan). J. Gt. Lakes. Res. 41, 100–109. doi:10.1016/j.jglr.2015.04.003
Maavara, T., Parsons, C. T., Ridenour, C., Stojanovic, S., Dürr, H. H., Powley, H. R., et al. (2015b). Global phosphorus retention by river damming. Proc. Natl. Acad. Sci. U S A 112, 15603–15608. doi:10.1073/pnas.1511797112
Maavara, T., Akbarzadeh, Z., and Van Cappellen, P. (2020a). Global dam-driven changes to riverine N:P:Si ratios delivered to the coastal ocean. Geophys. Res. Lett. 47, e2020GL088288. doi:10.1029/2020gl088288
Maavara, T., Chen, Q., Van Meter, K., Brown, L. E., Zhang, J., Ni, J., et al. (2020b). River dam impacts on biogeochemical cycling. Nat. Rev. Earth Environ. 1, 103–116. doi:10.1038/s43017-019-0019-0
Maranger, R., Jones, S. E., and Cotner, J. B. (2018). Stoichiometry of carbon, nitrogen, and phosphorus through the freshwater pipe: stoichiometry of carbon, nitrogen, and phosphorus. Limnol. Oceanogr. Lett. 3, 89–101. doi:10.1002/lol2.10080
Moberg, E. A., and Sosik, H. M. (2012). Distance maps to estimate cell volume from two-dimensional plankton images: distance map cell volume algorithm. Limnol. Oceanogr. Methods 10, 278–288. doi:10.4319/lom.2012.10.278
Officer, C., and Ryther, J. (1980). The possible importance of silicon in marine eutrophication. Mar. Ecol. Prog. Ser. 3, 83–91. doi:10.3354/meps003083
Paerl, H. W., Valdes, L. M., Peierls, B. L., Adolf, J. E., and Harding, L. W. (2006). Anthropogenic and climatic influences on the eutrophication of large estuarine ecosystems. Limnol. Oceanogr. 51, 448–462. doi:10.4319/lo.2006.51.1_part_2.0448
Poff, N. L., and Schmidt, J. C. (2016). How dams can go with the flow. Science 353, 1099–1100. doi:10.1126/science.aah4926
Royer, T. V. (2020). Stoichiometry of nitrogen, phosphorus, and silica loads in the Mississippi-Atchafalaya River basin reveals spatial and temporal patterns in risk for cyanobacterial blooms. Limnol. Oceanogr. 65, 325–335. doi:10.1002/lno.11300
Schelske, C. L., and Stoermer, E. F. (1971). Eutrophication, silica depletion, and predicted changes in algal quality in Lake Michigan. Science 173, 423–424. doi:10.1126/science.173.3995.423
Schelske, C. L., Stoermer, E. F., Fahnenstiel, G. L., and Haibach, M. (1986). Phosphorus enrichment, silica utilization, and biogeochemical silica depletion in the Great Lakes. Can. J. Fish. Aquatic Sci. 43, 407–415. doi:10.1139/f86-051
Scibona, A., Nizzoli, D., Hupfer, M., Valerio, G., Pilotti, M., and Viaroli, P. (2022). Decoupling of silica, nitrogen and phosphorus cycling in a meromictic subalpine lake (Lake Iseo, Italy). Biogeochemistry 159, 371–392. doi:10.1007/s10533-022-00933-9
Sethna, L. R., and Royer, T. V. (2024). Response of dissolved organic matter to variation in nutrient availability and phytoplankton community composition in a mesotrophic, temperate-zone reservoir. Freshw. Sci. 43, 412–426. doi:10.1086/732260
Sethna, L. R., Royer, T. V., Speir, S. L., Trentman, M. T., Mahl, U. H., Hagemeier, L. P., et al. (2022). Silicon concentrations and stoichiometry in two agricultural watersheds: implications for management and downstream water quality. Biogeochemistry 159, 265–282. doi:10.1007/s10533-022-00927-7
Sinha, E., Michalak, A. M., and Balaji, V. (2017). Eutrophication will increase during the 21st century as a result of precipitation changes. Science 357, 405–408. doi:10.1126/science.aan2409
Smith, V. H. (2003). Eutrophication of freshwater and coastal marine ecosystems a global problem. Environ. Sci. Pollut. Res. 10, 126–139. doi:10.1065/espr2002.12.142
Sultan, S. M. (2014). Determination of silicon and silicates. in L. M. L. Nollet, and L. S. P. De Gelder editors. Handbook of water analysis. 3rd ed. Boca Raton: Taylor and Francis, 24.
Teubner, K., and Dokulil, M. T. (2002). Ecological stoichiometry of TN:TP:SRSi in freshwaters: nutrient ratios and seasonal shifts in phytoplankton assemblages. Arch. für Hydrobiol. 154, 625–646. doi:10.1127/archiv-hydrobiol/154/2002/625
Thamatrakoln, K., and Kustka, A. B. (2009). When to say when: can excessive drinking explain silicon uptake in diatoms? BioEssays 31, 322–327. doi:10.1002/bies.200800185
Triplett, L. D., Engstrom, D. R., Conley, D. J., and Schellhaass, S. M. (2008). Silica fluxes and trapping in two contrasting natural impoundments of the upper Mississippi River. Biogeochemistry 87, 217–230. doi:10.1007/s10533-008-9178-7
Turner, R. E., and Rabalais, N. N. (1994). Coastal eutrophication near the Mississippi river delta. Nature 368, 619–621. doi:10.1038/368619a0
Turner, R. E., and Rabalais, N. N. (2004). Suspended sediment, C, N, P, and Si yields from the Mississippi River Basin. Hydrobiologia 511, 79–89. doi:10.1023/b:hydr.0000014031.12067.1a
Turner, R. E., Qureshi, N., Rabalais, N. N., Dortch, Q., Justic, D., Shaw, R. F., et al. (1998). Fluctuating silicate:nitrate ratios and coastal plankton food webs. Proc. Natl. Acad. Sci. U S A 95, 13048–13051. doi:10.1073/pnas.95.22.13048
Turner, R. E., Rabalais, N. N., Alexander, R. B., McIsaac, G., and Howarth, R. W. (2007). Characterization of nutrient, organic carbon, and sediment loads and concentrations from the Mississippi River into the northern Gulf of Mexico. Estuaries Coasts 30, 773–790. doi:10.1007/bf02841333
Wang, B., Liu, C.-Q., Maberly, S. C., Wang, F., and Hartmann, J. (2016). Coupling of carbon and silicon geochemical cycles in rivers and lakes. Sci. Rep. 6, 35832. doi:10.1038/srep35832
Yang, S. L., Milliman, J. D., Li, P., and Xu, K. (2011). 50,000 dams later: erosion of the Yangtze River and its delta. Glob. Planet. Change 75, 14–20. doi:10.1016/j.gloplacha.2010.09.006
Keywords: silicon, river damming, reservoirs, silicon retention, silicon flux, silicon budget
Citation: Sethna LR, Rasnake LM, Powers SR, Plechaty ZI, Pouchak AH, Hobbs WO and Royer TV (2025) Dynamics of dissolved silicon in a large mesotrophic reservoir in south-central Indiana, U.S.A. Front. Environ. Sci. 13:1656708. doi: 10.3389/fenvs.2025.1656708
Received: 30 June 2025; Accepted: 04 August 2025;
Published: 29 August 2025.
Edited by:
Zimin Li, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Xiaodong Zhang, Tianjin University, ChinaNatalia Borrelli, Universidad Nacional de Mar del Plata, Argentina
Copyright © 2025 Sethna, Rasnake, Powers, Plechaty, Pouchak, Hobbs and Royer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lienne R. Sethna, bHNldGhuYUBzbW0ub3Jn
†ORCID: Lienne R. Sethna, orcid.org/0000-0003-1156-172X; Lindsey M. Rasnake, orcid.org/0009-0003-8498-4878; William O. Hobbs, orcid.org/0000-0001-7321-0779; Todd V. Royer, orcid.org/0000-0002-1945-1560