Low Benthic Oxygen and High Internal Phosphorus-Loading are Strongly Associated With the Invasive Macrophyte Nitellopsis obtusa (starry stonewort) in a Large, Polymictic Lake

Lake Scugog is an important headwater to the Trent Severn-Waterway in Ontario, Canada. In recent years, notable ecosystem-level changes have occurred coinciding with the emergence of the non-native invasive charophyte Nitellopsis obtusa. Despite N. obtusa arriving in North America in the early 1970s, studies documenting the impact of N. obtusa on invaded ecosystems are scarce. Given the increasing dominance of N. obtusa in inland waters of the Great Lakes basin, we investigated the ecosystem-level impacts of N. obtusa in Lake Scugog over a 3-year period. We show for the first time a strong association between N. obtusa occurrence and biomass with benthic anoxia in this shallow, polymictic lake. Benthic dissolved oxygen concentrations were significantly lower (p-value < 0.001) at sites with N. obtusa compared to sites without N. obtusa. Additionally, N. obtusa biomass was a negative predictor of near-bed oxygen concentration (R2 = 0.59, p-value < 0.001). Knowing that anoxia can promote the internal loading of phosphorus, we measured soluble reactive phosphorus (SRP) in the pore-water of sediments at each site, and found N. obtusa biomass explained 90% of sediment pore-water SRP (R2 = 0.90, p-value < 0.001). These notable associations between N. obtusa and key lake elements indicates that N. obtusa may be acting as an ecosystem engineer in invaded lakes by altering the biogeochemical fate of oxygen and phosphorus.

CHANGE OF HEAT FLUX AT THE WATER-BOTTOM SEDIMENTS BOUNDARY IN A SMALL LAKE DURING THE YEAR

The features of the variability of the heat flux at the boundary of the water column with bottom sediments in the central part of the small polymictic Lake Vendyurskoe (southern Karelia) were analyzed according to the data of year-round temperature measurements from November 2018 to October 2019. It was shown that the heat flux directed from the bottom sediments into the water reached 6–8 W/m2 in the pre-ice period and at the moments of ice formation and breaking, and did not exceed 1 W/m2 during most of the winter. A week after the ice-off (May 3–4, 2019) the water column was completely mixed, the heat flux changed its sign, and heat accumulation in the upper layer of bottom sediments began. The heat flux directed to bottom sediments was 4–6 W/m2, periodically increasing to 8-12 W/m2 from mid-May to late July and from the second decade of August to the first decade of September. Against the background of a prolonged cooling during the first ten days of August, the heat flux changed sign with daily frequency, the amplitude of daily fluctuations reached ± 2.5 W/m2. Autumn cooling of the Lake Vendyurskoe began in early September 2019, the heat flux sign changed on September 14–16, 2019.

Autocalibration of a one-dimensional hydrodynamic-ecological model (DYRESM 4.0-CAEDYM 3.1) using a Monte Carlo approach: simulations of hypoxic events in a polymictic lake

Automated calibration of complex deterministic water quality models with a large number of biogeochemical parameters can reduce time-consuming iterative simulations involving empirical judgements of model fit. We undertook autocalibration of the one-dimensional hydrodynamic-ecological lake model DYRESM-CAEDYM, using a Monte Carlo sampling (MCS) method, in order to test the applicability of this procedure for shallow, polymictic Lake Rotorua (New Zealand). The calibration procedure involved independently minimizing the root-mean-square error (RMSE), maximizing the Pearson correlation coefficient (r) and Nash–Sutcliffe efficient coefficient (Nr) for comparisons of model state variables against measured data. An assigned number of parameter permutations was used for 10 000 simulation iterations. The “optimal” temperature calibration produced a RMSE of 0.54 ∘C, Nr value of 0.99, and r value of 0.98 through the whole water column based on comparisons with 540 observed water temperatures collected between 13 July 2007 and 13 January 2009. The modeled bottom dissolved oxygen concentration (20.5 m below surface) was compared with 467 available observations. The calculated RMSE of the simulations compared with the measurements was 1.78 mg L−1, the Nr value was 0.75, and the r value was 0.87. The autocalibrated model was further tested for an independent data set by simulating bottom-water hypoxia events from 15 January 2009 to 8 June 2011 (875 days). This verification produced an accurate simulation of five hypoxic events corresponding to DO < 2 mg L−1 during summer of 2009–2011. The RMSE was 2.07 mg L−1, Nr value 0.62, and r value of 0.81, based on the available data set of 738 days. The autocalibration software of DYRESM-CAEDYM developed here is substantially less time-consuming and more efficient in parameter optimization than traditional manual calibration which has been the standard tool practiced for similar complex water quality models.

Auto-calibration of a one-dimensional hydrodynamicecological model (DYRESM 4.0-CAEDYM 3.1) using a Monte Carlo approach: simulations of hypoxic events in a polymictic lake

Automated calibration of complex deterministic water quality models with a large number of biogeochemical parameters can reduce time-consuming iterative simulations involving empirical judgements of model fit. We undertook auto-calibration of the one-dimensional hydrodynamic-ecological lake model DYRESM-CAEDYM, using a Monte Carlo Sampling (MCS) method, in order to test the applicability of this procedure for shallow, polymictic Lake Rotorua (New Zealand). The calibration procedure involved independently minimising the root-mean-square-error (RMSE), maximizing the Pearson correlation coefficient (r) and Nash-Sutcliffe efficient coefficient (Nr) for comparisons of model state variables against measured data. An assigned number of parameter permutations was used for 10,000 simulation iterations. The optimal temperature calibration produced a RMSE of 0.54 °C, Nr-value of 0.99 and r-value of 0.98 through the whole water column based on comparisons with 540 observed water temperatures collected between 13 July 2007–13 January 2009. The modeled bottom dissolved oxygen concentration (20.5 m below surface) was compared with 467 available observations. The calculated RMSE of the simulations compared with the measurements was 1.78 mg L−1, the Nr-value was 0.75 and the r-value was 0.87. The autocalibrated model was further tested for an independent data set by simulating bottom-water hypoxia events for the period 15 January 2009 to 8 June 2011 (875 days). This verification produced an accurate simulation of five hypoxic events corresponding to DO < 2 mg L−1 during summer of 2009–2011. The RMSE was 2.07 mg L−1, Nr-value 0.62 and r-value of 0.81, based on the available data set of 738 days. The auto-calibration software of DYRESM-CAEDYM developed here is substantially less time-consuming and more efficient in parameter optimisation than traditional manual calibration which has been the standard tool practiced for similar complex water quality models.