Silica in Lake Superior: mass balance considerations and a model for dynamic response to eutrophication

The dissolved silica concentration in waters of Lake Superior probably is in a steady state because it is not influenced significantly by man, and the climate, topography and vegetation in the drainage area of the lake have been stable for the past 4000 years. Therefore the rate at which dissolved silica is introduced to the lake should equal the output rate.The primary inputs are: tributaries (4.1–4.6 × 10⁸kgSiO₂/yr), diffusion from sediment pore waters (0.21?0.78 × 10⁸kgSiO₂/yr) and atmospheric loading (0.26 × 10⁸kgSiO₂/yr). Silica is lost from the lake waters by: outflow through the St. Marys River, diatom deposition, adsorption onto particulates in the sediments, and authigenic formation of new silicate minerals. Tributary outflow accounts for less than one half the annual input of silica, and diatom deposition and silica adsorption withdraw less than 10% of the annual input. Therefore the formation of new silicate phases must be the dominant sink for dissolved silica in Lake Superior. The specific phases formed are not identified in the bottom sediments. X-ray diffraction studies suggest that smectite is one product, and amorphous ferroaluminum silicates may be another product.Mathematical modeling of the dissolved silica response to lake eutrophication suggests that the phosphate loading to Lake Superior would have to increase by about 250-fold to cause a silica depletion rate equal to that reported for Lake Michigan, assuming no change in the rate of upwelling of deep waters.