Next Article in Journal
ISSR-Based Genetic Diversity Assessment of Genus Jasminum L. (Oleaceae) from Pakistan
Next Article in Special Issue
Aquatic Macrophytes Occurrence in Mediterranean Farm Ponds: Preliminary Investigations in North-Western Sicily (Italy)
Previous Article in Journal
Identification and Functional Characterization of Genes Encoding Phenylacetaldehyde Reductases That Catalyze the Last Step in the Biosynthesis of Hydroxytyrosol in Olive
Previous Article in Special Issue
Wetland Restoration with Hydrophytes: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wind Exposure Regulates Water Oxygenation in Densely Vegetated Shallow Lakes

1
EA 4592 Géoressources & Environnement, F-33600 Pessac, France
2
INRAE, UR EABX, F-33612 Cestas, France
3
Pôle R&D Écosystèmes Lacustres (ECLA), F-13100 Aix-en-Provence, France
*
Author to whom correspondence should be addressed.
Plants 2021, 10(7), 1269; https://doi.org/10.3390/plants10071269
Submission received: 14 May 2021 / Revised: 14 June 2021 / Accepted: 18 June 2021 / Published: 22 June 2021
(This article belongs to the Special Issue Macrophytes in Inland Waters: From Knowledge to Management)

Abstract

:
The presence of dense macrophyte canopies in shallow lakes locally generates thermal stratification and the buildup of labile organic matter, which in turn stimulate the biological oxygen demand. The occurrence of hypoxic conditions may, however, be buffered by strong wind episodes, which favor water mixing and reoxygenation. The present study aims at explicitly linking the wind action and water oxygenation within dense hydrophytes stands in shallow lakes. For this purpose, seasonal 24 h-cycle campaigns were carried out for dissolved gases and inorganic compounds measurements in vegetated stands of an oligo-mesotrophic shallow lake. Further, seasonal campaigns were carried out in a eutrophic shallow lake, at wind-sheltered and -exposed sites. Overall results showed that dissolved oxygen (DO) daily and seasonal patterns were greatly affected by the degree of wind exposure. The occurrence of frequent wind episodes favored the near-bottom water mixing, and likely facilitated mechanical oxygen supply from the atmosphere or from the pelagic zone, even during the maximum standing crop of plants (i.e., summer and autumn). A simple model linking wind exposure (Keddy Index) and water oxygenation allowed us to produce an output management map, which geographically identified wind-sheltered sites as the most subjected to critical periods of hypoxia.

1. Introduction

In lentic shallow water bodies, the diel and seasonal oxygen balance is given by the interplay between the photosynthetic activity of primary producers (net production of dissolved oxygen, DO), their respiration (net DO consumption) and heterotrophic respiration of bacteria and animals (net DO consumption). When present, submerged aquatic vegetation (SAV) induces significant diel fluctuations in oxygen levels [1]. During the day, a supersaturation of oxygen (>100%) is observed due to photosynthesis, while at night oxygen is no longer produced, and consumption processes are predominant due to respiration. This type of nycthemeral variation is particularly emphasized in summer, when plant photosynthetic rate and heterotrophic respiration are at their maximum, and whose net effect largely exceeds the contribution of temperature-dependent oxygen solubility [2,3]. Primary production releases high amounts of oxygen in the water column, allowing for the oxidation of methane (CH4) by methanotrophic epiphytic bacteria [4,5]. In addition, radial oxygen loss in the rhizosphere [6,7,8] contributes to reduce benthic CH4 flux through benthic methanotrophy or oxidation of nitrate [9]. The synthesis of large quantities of biomass occurs through the assimilation of nutrients (including N-compounds, phosphate and carbon dioxide, CO2); SAV is thus able to uptake gases and nutrients coming from the sediment and the atmosphere, and synthesize them in biomass [10]. However, oxygen dynamics can be altered in densely vegetated stands, such as those dominated by invasive macrophytes: here, due to a fast cycle of growth and decay, fresh organic matter continuously replenishes the organic bulk in the sediment [11]. In those conditions, elevated biological oxygen demand (BOD) and hypoxia are coupled to accumulation and stratification of nutrients, such as carbon dioxide, methane, ammonium and reactive phosphorous, also during the day [12,13,14,15].
Dense plant canopies are known to locally generate bottom shading and modify water circulation, that impeding convective cooling, even within small depths [16,17,18,19,20]. This can affect or, inversely, exacerbate the thermally driven exchange flow of nutrients and DO between pelagic and littoral zones of the lake [21]. Nevertheless, it has been assessed that wind episodes may induce surface and internal flows even within stratified waters [22]. Shallow lakes are continuously subjected to mixing and wave-breaking in function of wind speed [21]; in function of the lake size, slope and bathymetry, wind-exposed lakes are also concerned by periodic seiche events which contribute to water mixing [23,24]. Still, only few studies infer that wind action may induce local turbulent mixing and reaeration even within dense submerged canopies [14,18]. Indeed, the explicit interaction between wind exposure and ecosystem functioning in shallow plant-dominated lakes remains almost undescribed [16,25], despite the growing demand for precision by stakeholders (users, managers, politicians) in the domain of biological invasions.
Invasive aquatic plants management is one of the main issue concerning global changes in freshwaters [26]. In the past, manager’s decisions on invaded environments were mainly driven by socio-economical questions relating to tourism, boating and swimming [27]. Recently, invasive plants management started to progressively embrace ecological sciences, with the primary goal of understanding if invaded sites are de facto degraded or imperiled regarding to their functioning. As a result, aquatic weeds are more and more in question about their overall role on ecosystem metabolism. Coherently, managers are demanding to ecologists to produce effective and readable tools for ameliorating their interventions [28,29]. We herein report a study specifically addressing this issue, which employs DO saturation as a reliable indicator of net ecosystem metabolism, related to autotrophic and heterotrophic processes [30]. Firstly, we hypothesize that the impact on DO levels within invasive macrophyte stands significantly differs in function of the degree of wind exposure. Secondly, we hypothesize that DO levels are inversely correlated to plant densities and sedimentary organic matter content. Thirdly, we propose a quantitative tool to spatially identify sites that are more risky for hypoxic events, and thus need intervention by managers.

2. Results

2.1. Diel Variations in Vegetated Stands

Sampling sites in Lacanau Lake (hereafter, LAC Lake) were homogeneously distributed within the largest invasive macrophyte stands of the lake, which developed in the most sheltered zones of the lake (Figure 1).
Results from seasonal 24 h-cycle campaigns showed that most of the sites were hypoxic (DO saturation <100%), and that DO depletion also occurred during daylight (Figure 2). Concomitantly, CO2 was mostly supersaturated and pH acid (pH < 7), with some exceptions during daylight in summer and autumn (Figures S1 and S2). CH4, NH4+ and NO3 buildup in the water column appeared both during the night and day (Figures S3–S5). Water temperature ranged from 11.1 ± 0.2 to 26.7 ± 0.3 °C (in spring and summer, respectively—Figure S6) and DOC averaged 13.1 ± 0.2 mg L−1 on an annual basis.
ANOVA test revealed that, with temperature and ammonium as solely exceptions (day > night), dissolved gases and inorganic compounds measured within plant stands did not differ between day and night; all parameters varied seasonally (Table 1).

2.2. Wind-Sheltered vs. Wind-Exposed Sites

The choice of sampling sites in Parentis-Biscarrosse Lake (hereafter, PAR Lake) was based on two co-occurring conditions: the presence of densely vegetated areas and the difference in wind exposure (Figure 1). Results from seasonal campaigns in PAR Lake showed that dissolved gases and inorganic compounds significantly changed in function of wind exposure (ANOVA, Table 2).
Differences between sheltered and exposed sites were significant for every physicochemical parameter, yet only at vegetated sites and in function of the season (in summer and in autumn). Significant differences between vegetated and plant-free sites occurred only at sheltered sites. Tukey’s HSD test indicated that water temperature was lower at sheltered sites than at exposed ones (Figure S7). pH values were lower at vegetated and sheltered sites than at exposed ones, only during summer (Figure S8). DO was lower at vegetated and sheltered sites than at exposed ones (Figure 3); CO2 and CH4 were higher at vegetated and sheltered sites than at exposed ones (Figures S9 and S10); NH4+ and NO3 values differed seasonally between sheltered and exposed sites, with no differences between vegetated and plant-free areas (Figures S11 and S12). DOC averaged 6.2 ± 0.2 mg L−1 on an annual basis.

2.3. Dependence of DO Saturation on Plant Biomass and Sedimentary OM

Total biomass varied seasonally at both lakes, with values comprised between 319 ± 245 and 668 ± 414 gDW m−2 at LAC Lake (in spring and summer, respectively), and between 1626 ± 132 and 4528 ± 2413 gDW m−2 at PAR Lake (in spring at exposed and in autumn at sheltered sites, respectively). OM content in vegetated sediments ranged from 0.7 ± 0.2 to 71 ± 3% and from 0.7 ± 0.1 to 1.2 ± 0.1% as LOI, for LAC Lake and PAR Lake, respectively. Linear mixed-effects model, calculated on the two lakes dataset, showed that DO saturation was not dependent on OM sedimentary content on the total plant biomass; only DO values measured in LAC Lake during summer resulted in being negatively correlated to biomass (p-value < 0.01).

2.4. Dependence of DO on Wind Exposure and Hypoxia Risk Map Production

The regression of DO saturation against wind exposure, identified with the segmented function in R, showed a structural breakpoint at Keddy Index = 2.9 (Figure 4). We considered this breakpoint as a threshold of hypoxia risk, i.e., low risk above this value and high risk below. This threshold is assumed to be the minimum wind exposure which would be able to decrease the risk of hypoxia in dense submerged plant stands.
Further, in order to produce a hypoxia risk map, Keddy Index was calculated for each 4 h-long period (n = 2190) on each pixel cell (n = 4031 for LAC Lake and n = 14,438 for PAR Lake) matching with densely vegetated areas, presenting biomass >50 gDW m−2, mapped at the lake scale (1.19 km2 and 4.17 km2 in LAC and PAR Lakes, respectively, from 31) (Figure 5). Hypoxia risk was above 50% in 70 ha of plants stands (corresponding to 60% of the total vegetated surface) in LAC Lake and in 50 ha in PAR Lake (12% of the total vegetated surface). This risk was above 75% in 11 ha of plants stands (9% of the total vegetated surface) in LAC Lake and in 11ha in PAR Lake (3% of the total vegetated surface).

3. Discussion

In vegetated stands, diel variations of inorganic compounds typically reflect plants photosynthetic activity, with the lowest dissolved carbon and nitrogen concentrations measured in the water at late afternoon, corresponding to nutrients depletion by plants uptake, then an accumulation during the night, with a peak just before dawn. At the same time, DO and pH follow an exactly inverse pattern. In our study, the described nycthemeral shape was detectable only at some sites and mostly during summer. At other sites, heterotrophic activity, stimulated by temperature increase during summer and autumn, exceeded net oxygen release during the day, that resulting in hypoxia/anoxia events and buildup of CO2, CH4 and NH4+ in the water column. This observation is recurrent in dense stands formed by invasive macrophytes, where the sedimentation of organic matter generates an elevated benthic BOD during the period of senescence of plants; this implies a permanent DO deficit [15,31,32]. In dense hydrophyte stands, DO input from the atmosphere can be limited to the surficial layer of the water column, as long stems constitute a physical barrier, as floating-leaved macrophytes do [13]. In our case, the occurrence of a vertical “plant wall” at the external boundaries of vegetated stands may also lead to the annihilation of the horizontal flow of nutrients and DO from the pelagic to the littoral zones [22].
Hypoxic events and inorganic compounds buildup can be however contrasted by the wind action, which may induce local turbulent mixing and reaeration even within dense submerged canopies [14,18]. Coherently, some of the diel variations measured in our study showed a flattened shape, with constant values along the 24 h-cycle. On one hand, elevated DO values during the night could be attributable to convective mixing due to air temperature nightly decrease [17,33]. On the other hand, the maintaining of constant DO values along a diel cycle may be an indicator of stationary wind conditions and turbulent mixing; this supposition is supported by the second part of our study. Seasonal campaigns at wind-sheltered and -exposed sites showed that, ecosystem functioning was not ascribable to the plant presence/absence or to the seasonal biomass variation only. Indeed, DO and CO2 saturation at wind-exposed sites hovered at about 100% all year round, indicating that wind-driven diffusion continuously outreached net production and consumption within the water column, even in invaded areas of the lake. Overall results show thus that the presence of invasive hydrophytes does not systematically promote water hypoxia, if local wind conditions allow an efficient water mixing by wind.
When considering the whole dataset, only DO values measured in LAC Lake during summer resulted in being dependent on plant density; moreover, vegetated stands in this lake mainly developed at sheltered sites [34]. Prevailing winds oriented from the northwest created low hydrodynamic conditions, because of the natural barrage formed by sand dunes [35]. Elevated plant biomass matching with shallow depths in wind-sheltered areas seemed to generate favorable conditions for water hypoxia, a phenomenon exacerbated by an elevated turnover of biomass during summer. In contrast, extremely elevated biomass measured in PAR Lake, largely exceeding values reported until now for Egeria spp. invaded sites [3,31], did not generate an extreme DO deficit even at wind-sheltered sites. The difference between the two lakes is evident also from a thermic point of view: at LAC Lake, a previous study had showed that water temperature measured in vegetated stands was significantly lower than that measured in plant-free areas, irrespective of the season [19]. The present study on PAR Lake shows instead that no significant difference exists between vegetated and plant-free areas, irrespective of the season (Figure S7). As for the DO and CO2, the divergence in temperature results among the two lakes could be due to the different size of the lake, the second being larger and permitting fetch length—and thus, water mixing—to be more important.
The hypoxia risk map shows that elevated hypoxia probability is associated with wind-sheltered areas of the lakes, and that oxygenation shortage can affect a large total surface of several tens of hectares. Hypoxia risk is at its maximum in both lakes at enclosed and wind-sheltered areas, like small marinas and public boat launches, which are known to be important drivers of aquatic plant spread [36,37]. On the other hand, large surfaces of the lake invaded by elevated plant densities would not be affected by hypoxia and would thus not necessitate intervention. The hypoxia risk map we produced represents a preliminary and concrete tool, coupling field measurements and modelling, which can reduce plant management costs, as it indicates precisely where invasive plants constitute a problem for ecosystem functioning. A similar approach providing reproducible management tools has been recently published [38], that coupling lake depth or bathymetry to anoxia probability in the hypolimnion of deep lakes. Our model should be, however, calibrated site-specifically, because the intrinsic sedimentary features and the trophic status of the lake could affect the magnitude of hypoxia level and nutrients flux. The two lakes we studied presented different DOC values, sedimentary OM content and resulted in very different concentrations of CO2 and CH4. Also, due to a different fetch length, the reaeration strongly varied even within comparable wind velocity. A possible improvement of our method could have been to introduce the local bathymetry in the model. Indeed, waves induce vertical upward forces acting on the water column movements and sediment resuspension [34]; furthermore, wind-induced circulation in nearshore zones appears to be crucial in littoral plant-free areas [24]. We can expect an increase of wind effect on water mixing in shallow zones due to orbital movements translated to the lake bottom. Nevertheless, SAV also reduce waves action and current velocities within beds [39,40]. Future modelling works should thus focus on integrating vegetation in the photic region to better define how the cross-shore water circulation works. Another possible improvement in the future would be the use of automatic oxygen probes, in order to obtain a finer resolution scale of diel and seasonal variations, and perfect the calculation of hypoxia risk probability on a long temporal scale.
Our results highlight the need to consider local hydrodynamics in lake management decisions. Wind exposure should be used for spatially organizing management plans and prioritizing zones where invasive biomass control actions are needed. Mapping hypoxia risk in densely vegetated stands is a promising tool for the management of invasive hydrophytes in shallow lakes.

4. Materials and Methods

4.1. Study Area

Lacanau Lake and Parentis-Biscarrosse Lake are shallow lakes located in the southern Atlantic coast of France. Those lakes are characterized by sandy acidic substrate and classed as oligo-mesotrophic (Lacanau, 16 km2) and eutrophic (Parentis-Biscarrosse, 32 km2). Within the two lakes, large submerged stands of Egeria densa Planch. and Lagarosiphon major (Ridl.) Moss develop between 1 and 7 m deep, with dense stands being preferentially located at shallow and wind-sheltered sites, or at deep and wind-exposed sites [34].

4.2. Field Campaigns

Between June 2013 and November 2015 at Lacanau Lake, seasonal 24 h-cycle campaigns were carried out at 15 sites. Sampling sites were homogeneously distributed within the largest invasive macrophyte stands of the lake [34]. Water was collected within plant canopy at depths ranging from 100 to 330 cm, with a frequency of four times a day (two samplings during the day, between 11 a.m. and 3 p.m.; two samplings during the night, between 9 p.m. and 6 a.m.). Water temperature (T, °C), pH, dissolved oxygen (DO, expressed as saturation %), dissolved carbon dioxide (CO2, %), dissolved methane (CH4, µM), nitrate (NO3, µM), ammonium (NH4+, µM) and dissolved organic carbon (DOC, mg L−1) were measured according methods reported in [19]. Finally, we tested the influence of the sampling time on the biogeochemistry of the water column by a two-way ANOVA with interactions among factors. The diel variation (two levels: day vs. night) and the season (three levels: spring vs. summer vs. autumn) were considered as fixed factors, while the sampling site (fifteen levels) was considered as a random factor.
Between January 2016 and January 2017 at Parentis-Biscarrosse Lake, seasonal sampling campaigns were carried out, during the day only, at vegetated (3 wind-sheltered and 3 wind-exposed sites) and at plant-free sites (3 wind-sheltered and 3 wind-exposed). The degree of wind exposure was estimated by previous modeling of wind exposure Keddy Index [41]. Water was collected within plant canopy at depths ranging from 150 to 300 cm, between 11 a.m. and 3 p.m. Water samples collection, treatment and analyses are the same as those adopted in Lacanau Lake and reported in [19]. Finally, we tested the influence of spatial exposure to wind on the biogeochemistry of the water column by a three-way ANOVA with interactions among factors. The degree of wind exposure (two levels: exposed vs. sheltered), the plant presence (two levels: vegetated vs. plant-free) and the season (four levels: winter vs. spring vs. summer vs. autumn) were considered as fixed factors, while the sampling site (twelve levels) was considered as a random factor.
Normal distribution (Shapiro–Wilk Test) and homoscedasticity (Levene’s Test) were tested before running ANOVAs. Post hoc analyses were performed by Tukey’s Honestly Significant Difference (HSD) test. Statistical analyses were performed with R Program [42]. Mean values are reported with their standard deviation.
Macrophytes sampling was carried out by rake for total biomass (gDW m−2) measurements, immediately after water samplings, as reported in [19]. Concomitantly, within the plant stands, sediment samples were collected by grabber, as described in [34], for sedimentary organic matter (OM, as loss of ignition, % LOI) measurements. In order to test the dependence of DO saturation on plant biomass and OM content, a linear mixed-effects model fit by maximum likelihood was performed on the whole dataset (DO measurements from Lacanau and Parentis-Biscarrosse Lakes), with the sampling site as a random factor.

4.3. Wind Exposure Calculations

Wind exposure was calculated according [41] for both lakes by using a fetch matrix (i.e., the distance over which waves can build up) obtained from lake open-water raster grid cells (resolution of 17 m) for each wind compass direction (10–360°, in 10° increments). Wind data (hourly and daily mean speed and direction) were provided by Météo France in Cap-Ferret (44°37′54″ N, 1°14′53″ O) and Biscarrosse (44°25′54″ N, 1°14′51″ O) weather stations for Lacanau and Parentis-Biscarrosse Lakes, respectively. It is possible to generate values which should be related to the effect of wind at a given point (here, a grid cell) by using fetch and wind velocity. For a given compass direction, one measure of exposure is the product of mean wind speed and direction and the percent frequency of the wind blowing in that direction.
In order to position wind-sheltered and -exposed sampling sites in Parentis-Biscarrosse Lake, daily mean wind speed and direction were used to build a wind exposure map during 1-year period (2014). One measure of exposure was calculated for each grid cell over 36 compass wind directions according to a fetch matrix. A cell’s total exposure is given by the sum of values calculated for all of the compass directions during the 1-year period. Sampling sites were chosen within lake areas identified as low- or highly-exposed to wind action.

4.4. Coupling DO and Wind Exposure

Keddy Index was calculated for the 4 h before the exact timing of water sampling. Then, each DO value was coupled to the sum of Keddy Index values for this period. This is the duration estimated being necessary for the water mixing at low depth [17,43]. In order to test the dependence of DO saturation on wind exposure, a Chow test was performed to determine the presence of a structural break at some point of the data series [44]. We used the sctest function from the strucchange package in R software to perform a Chow test, which resulted in F = 10.7, p-value = 2.7 × 10−5. The significance of the test indicates that a structural breakpoint is present in the regression. Or else, that two regression lines can fit the pattern in the data more effectively than a single regression line. Finally, we applied the segmented function in R to analyze segmented relationships in the regression, in order to obtain a breakpoint value.

4.5. Hypoxia Risk Map Production

We calculated 4 h-long Keddy Index values each day during one year (2014 and 2016 for Lacanau and Parentis-Biscarrosse Lakes, respectively) for each raster cell corresponding to densely vegetated areas of the lake, and presenting a biomass >50 gDW m−2 [34]. Each 4 h-long period and each grid cell in which wind exposure was under the breakpoint value, indicating a high risk of hypoxia, was classified as “1”, whereas 4 h-long period with low risk of hypoxia (>2.9) were classified as “0”. The probability of hypoxia was expressed as the percentage (0–100%) of 4 h-long periods where wind exposure was below the hypoxia threshold during one year. Finally, this probability was reported on raster grid cells to map the hypoxia risk at densely vegetated areas scale.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10071269/s1, Figure S1: CO2 results from seasonal 24 h-cycle campaigns in Lacanau Lake; Figure S2: pH results from seasonal 24 h-cycle campaigns in Lacanau Lake; Figure S3: CH4 results from seasonal 24 h-cycle campaigns in Lacanau Lake; Figure S4: NH4+ results from seasonal 24 h-cycle campaigns in Lacanau Lake; Figure S5: NO3 results from seasonal 24 h-cycle campaigns in Lacanau Lake; Figure S6: Water temperature results from seasonal 24 h-cycle campaigns in Lacanau Lake; Figure S7: Water temperature results from seasonal campaigns in Parentis-Biscarrosse Lake; Figure S8: pH results from seasonal campaigns in Parentis-Biscarrosse Lake; Figure S9: CO2 results from seasonal campaigns in Parentis-Biscarrosse Lake; Figure S10: CH4 results from seasonal campaigns in Parentis-Biscarrosse Lake; Figure S11: NH4+ results from seasonal campaigns in Parentis-Biscarrosse Lake; Figure S12: NO3 results from seasonal campaigns in Parentis-Biscarrosse Lake. In S7, S8, S9, S10, S11 and S12, Tukey’s HSD results are not reported for the seasonality factor. *** indicates p-value < 0.001, ** p-value < 0.01, * p-value < 0.05.

Author Contributions

C.R., V.B. and J.T.-R. conceptualized the project and methodology, M.E. and G.J. contributed to the lab and the field work. C.R. prepared the original draft; V.B., J.T.-R., M.E. and G.J. reviewed and edited the final version. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by AEAG (Agence de l’Eau Adour-Garonne), within the conventions #310330122, #310330085, and #310330109, by Irstea (Institut National de Recherche en Sciences et Technologies pour l’Environnement et l’Agriculture), by the Regional Council Nouvelle-Aquitaine (CCRRDT) and by the Regional Research and Development Program (PSDR) within the project AquaVIT (#2015-1R20603—00005190-00005191).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article or supplementary material.

Acknowledgments

We would like to thank SIAEBVELG and Police du Lac de Lacanau, for facilitating night activities in the field. We are also grateful to B. Delest, A. Moreira, and K. Madarassou for nutrients analyses and to D. Poirier (EPOC, University of Bordeaux) for assisting in CH4 analyses. J.-L. Keller, S. Moreira, T. Huguet, G. Ducasse, A. Poli, H. Guerreiro, J. Chabanne, J. Vedrenne, I. Halgand, R. Le Barh, J. Gueguen, N. Dagens and Y. Bourguignon considerably helped during field and laboratory activities and we sincerely acknowledge them. We are grateful to T. Leboucher and S. Boutry for contributing to the implementation of the hypoxia risk map. We are grateful to students A. Poli, V. Mauro and N. Naulin for preliminary data analysis. We acknowledge Météo France for wind data supplying. We are grateful to H. Ulloa and another anonymous reviewer for the useful comments on the original draft of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Caraco, N.F.; Cole, J.J. Contrasting impacts of a native and alien macrophyte on dissolved oxygen in a large river. Ecol. Appl. 2002, 12, 1496–1509. [Google Scholar] [CrossRef]
  2. Pokorný, J.; Rejmánková, E. Oxygen regime in a fishpond with duckweeds (Lemnaceae) and Ceratophyllum. Aquat. Bot. 1983, 17, 125–137. [Google Scholar] [CrossRef]
  3. Pelicice, F.M.; Agostinho, A.A.; Thomaz, S.M. Fish assemblages associated with Egeria in a tropical reservoir: Investigating the effects of plant biomass and diel period. Acta Oecologica 2005, 27, 9–16. [Google Scholar] [CrossRef]
  4. Sorrell, B.K.; Downes, M.T.; Stanger, C.L. Methanotrophic bacteria and their activity on submerged aquatic macrophytes. Aquat. Bot. 2002, 72, 107–119. [Google Scholar] [CrossRef]
  5. Yoshida, N.; Iguchi, H.; Yurimoto, H.; Murakami, A.; Sakai, Y. Aquatic plant surface as a niche for methanotrophs. Front. Microbiol. 2014, 5, 30. [Google Scholar] [CrossRef]
  6. Sand-Jensen, K.; Prahl, C.; Stokholm, H. Oxygen release from roots of submerged aquatic macrophytes. Oikos 1982, 38, 349–354. [Google Scholar] [CrossRef]
  7. Sorrell, B.K.; Dromgoole, F.I. Oxygen transport in the submerged freshwater macrophyte Egeria densa Planch. I. Oxygen production, storage and release. Aquat. Bot. 1987, 28, 63–80. [Google Scholar] [CrossRef]
  8. Ribaudo, C.; Bartoli, M.; Racchetti, E.; Longhi, D.; Viaroli, P. Seasonal fluxes of O2, DIC and CH4 in sediments with Vallisneria spiralis: Indications for radial oxygen loss. Aquat. Bot. 2011, 94, 134–142. [Google Scholar] [CrossRef]
  9. Haroon, M.F.; Hu, S.; Shi, Y.; Imelfort, M.; Keller, J.; Hugenholtz, P.; Yuan, Z.; Tyson, G.W. Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature 2013, 500, 567–570. [Google Scholar] [CrossRef] [PubMed]
  10. Reddy, K.R.; Tucker, J.C.; Debusk, W.F. The role of Egeria in removing nitrogen and phosphorus from nutrient enriched waters. J. Aquat. Plant Manag. 1987, 25, 14–19. [Google Scholar]
  11. Drexler, J.Z.; Khanna, S.; Lacy, J.R. Carbon storage and sediment trapping by Egeria densa Planch, a globally invasive, freshwater macrophyte. Sci. Total Environ. 2020, 755, 142602. [Google Scholar] [CrossRef]
  12. Barko, J.W.; Godshalk, G.L.; Carter, V.; Rybicki, N.B. Effects of Submersed Aquatic Macrophytes on Physical and Chemical Properties of Surrounding Water; Aquatic Plant Control Research Program Technical Report A-88-1 I; Depart Army, U.S. Army Corps of Engineers: Washington, DC, USA, 1988; p. 28.
  13. Frodge, J.D.; Thomas, G.L.; Pauley, G.B. Effects of canopy formation by floating and submergent aquatic macrophytes on the water quality of two shallow Pacific Northwest lakes. Aquat. Bot. 1990, 38, 231–248. [Google Scholar] [CrossRef]
  14. Sand-Jensen, K.; Andersen, M.R.; Martinsen, K.T.; Borum, J.; Kristensen, E.; Kragh, T. Shallow plant-dominated lakes–extreme environmental variability, carbon cycling and ecological species challenges. Ann. Bot. 2019, 124, 355–366. [Google Scholar] [CrossRef] [Green Version]
  15. Moore, T.P.; Clearwater, S.J.; Duggan, I.C.; Collier, K.J. Invasive macrophytes induce context-specific effects on oxygen, pH, and temperature in a hydropeaking reservoir. River Res. Appl. 2020, 36, 1717–1729. [Google Scholar] [CrossRef]
  16. Herb, W.R.; Stefan, H.G. Dynamics of vertical mixing in a shallow lake with submersed macrophytes. Water Resour. Res. 2005, 41. [Google Scholar] [CrossRef]
  17. Coates, M.J.; Folkard, A.M. The effects of littoral zone vegetation on turbulent mixing in lakes. Ecol. Model. 2009, 220, 2714–2726. [Google Scholar] [CrossRef]
  18. Andersen, M.R.; Sand-Jensen, K.; Woolway, R.I.; Jones, I.D. Profound daily vertical stratification and mixing in a small, shallow, wind-exposed lake with submerged macrophytes. Aquat. Sci. 2017, 79, 395–406. [Google Scholar] [CrossRef]
  19. Ribaudo, C.; Tison-Rosebery, J.; Buquet, D.; Jan, G.; Jamoneau, A.; Abril, G.; Anschutz, P.; Bertrin, V. Invasive aquatic plants as ecosystem engineers in an Oligo-Mesotrophic shallow lake. Front. Plant Sci. 2018, 9, 1781. [Google Scholar] [CrossRef] [PubMed]
  20. Vilas, M.P.; Marti, C.L.; Oldham, C.E.; Hipsey, M.R. Macrophyte-induced thermal stratification in a shallow urban lake promotes conditions suitable for nitrogen-fixing cyanobacteria. Hydrobiologia 2018, 806, 411–426. [Google Scholar] [CrossRef]
  21. Zhang, X.; Nepf, H.M. Thermally driven exchange flow between open water and an aquatic canopy. J. Fluid Mech. 2009, 632, 227. [Google Scholar] [CrossRef]
  22. Ulloa, H.N.; Constantinescu, G.; Chang, K.; Horna-Munoz, D.; Hames, O.; Wüest, A. Horizontal transport under wind-induced resonance in stratified waterbodies. Phys. Rev. Fluids 2020, 5, 054503. [Google Scholar] [CrossRef]
  23. Lin, Y.T. Wind effect on diurnal thermally driven flow in vegetated nearshore of a lake. Environ. Fluid Mech. 2015, 15, 161–178. [Google Scholar] [CrossRef]
  24. Ulloa, H.N.; Constantinescu, G.; Chang, K.; Horna-Munoz, D.; Steiner, O.S.; Bouffard, D.; Wüest, A. Hydrodynamics of a periodically wind-forced small and narrow stratified basin: A large-eddy simulation experiment. Environ. Fluid Mech. 2019, 19, 667–698. [Google Scholar] [CrossRef] [Green Version]
  25. Torma, P.; Wu, C.H. Temperature and circulation dynamics in a small and shallow lake: Effects of weak stratification and littoral submerged macrophytes. Water 2019, 11, 128. [Google Scholar] [CrossRef] [Green Version]
  26. Havel, J.E.; Kovalenko, K.E.; Thomaz, S.M.; Amalfitano, S.; Kats, L.B. Aquatic invasive species: Challenges for the future. Hydrobiologia 2015, 750, 147–170. [Google Scholar] [CrossRef]
  27. Matzek, V.; Covino, J.; Funk, J.L.; Saunders, M. Closing the knowing–doing gap in invasive plant management: Accessibility and interdisciplinarity of scientific research. Conserv. Lett. 2014, 7, 208–215. [Google Scholar] [CrossRef] [Green Version]
  28. Young, S.L.; Kettenring, K.M. The social-ecological system driving effective invasive plant management: Two case studies of non-native Phragmites. J. Environ. Manag. 2020, 267, 110612. [Google Scholar] [CrossRef] [PubMed]
  29. Le Floch, S.; Ginelli, L. The victorious battles of the lost war against aquatic invasive plants. ‘Fluid’ categorisation and multiple forms of ordinary commitment. Trans. Inst. Br. Geogr. 2021. [Google Scholar] [CrossRef]
  30. Cole, J.J.; Pace, M.L.; Carpenter, S.R.; Kitchell, J.F. Persistence of net heterotrophy in lakes during nutrient addition and food web manipulations. Limnol. Oceanogr. 2000, 45, 1718–1730. [Google Scholar] [CrossRef] [Green Version]
  31. Bini, L.M.; Thomaz, S.M.; Carvalho, P. Limnological effects of Egeria najas planchon (Hydrocharitaceae) in the arms of Itaipu Reservoir (Brazil, Paraguay). Limnology 2010, 11, 39–47. [Google Scholar] [CrossRef]
  32. Wang, H.; Lu, J.; Wang, W.; Yang, L.; Yin, C. Methane fluxes from the littoral zone of hypereutrophic Taihu Lake, China. J. Geophys. Res. Atmos. 2006, 111. [Google Scholar] [CrossRef]
  33. Brothers, S.; Kazanjian, G.; Köhler, J.; Scharfenberger, U.; Hilt, S. Convective mixing and high littoral production established systematic errors in the diel oxygen curves of a shallow, eutrophic lake. Limnol. Oceanogr. Methods 2017, 15, 429–435. [Google Scholar] [CrossRef]
  34. Bertrin, V.; Boutry, S.; Jan, G.; Ducasse, G.; Grigoletto, F.; Ribaudo, C. Effects of wind-induced sediment resuspension on distribution and morphological traits of aquatic weeds in shallow lakes. J. Limnol. 2017, 76, 84–96. [Google Scholar] [CrossRef] [Green Version]
  35. Tastet, J.P.; Lalanne, R.; Maurin, B.; Dubos, B. Geological and archaeological chronology of a late Holocene coastal enclosure: The Sanguinet lake (SW France). Geoarchaeology 2008, 23, 131–149. [Google Scholar] [CrossRef]
  36. Tamayo, M.; Olden, J. Forecasting the vulnerability of lakes to aquatic plant invasions. Invasive Plant Sci. Manag. 2014, 7, 32–45. [Google Scholar] [CrossRef]
  37. Bertrin, V.; Boutry, S.; Alard, D.; Haury, J.; Jan, G.; Moreira, S.; Ribaudo, C. Prediction of macrophyte distribution: The role of natural versus anthropogenic physical disturbances. Appl. Veg. Sci. 2018, 21, 395–410. [Google Scholar] [CrossRef]
  38. Deeds, J.; Amirbahman, A.; Norton, S.A.; Suitor, D.G.; Bacon, L.C. Predicting anoxia in low-nutrient temperate lakes. Ecol. Appl. 2021. [Google Scholar] [CrossRef]
  39. Carper, G.L.; Bachmann, R.W. Wind resuspension of sediments in a prairie lake. Can. J. Fish. Aquat. Sci. 1984, 41, 1763–1767. [Google Scholar] [CrossRef]
  40. Schutten, J.; Dainty, J.; Davy, A. Wave-induced hydraulic forces on submerged aquatic plants in shallow lakes. Ann. Bot. 2004, 93, 333–341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Keddy, P.A. Quantifying within-lake gradients of wave energy: Interrelationships of wave energy, substrate particle size and shoreline plants in axe lake, Ontario. Aquat. Bot. 1982, 14, 41–58. [Google Scholar] [CrossRef]
  42. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 23 May 2021).
  43. Jalil, A.; Li, Y.; Du, W.; Wang, W.; Wang, J.; Gao, X.; Khan, H.O.S.; Pan, B.; Acharya, K. The role of wind field induced flow velocities in destratification and hypoxia reduction at Meiling Bay of large shallow Lake Taihu, China. Environ. Pollut. 2018, 232, 591–602. [Google Scholar] [CrossRef] [PubMed]
  44. Chow, G.C. Tests of equality between sets of coefficients in two linear regressions. Econom. J. Econom. Soc. 1960, 591–605. [Google Scholar] [CrossRef]
Figure 1. Keddy Index calculated on annual basis for Lacanau Lake (on the left) and Parentis-Biscarrosse Lake (on the right). The windrose is calculated on wind speed and direction hourly data on an annual basis. Lake bathymetry and sampling sites for seasonal 24-cycle campaigns (LAC Lake), as well as sampling sites for seasonal campaigns at wind-sheltered and -exposed sites (PAR Lake) are reported.
Figure 1. Keddy Index calculated on annual basis for Lacanau Lake (on the left) and Parentis-Biscarrosse Lake (on the right). The windrose is calculated on wind speed and direction hourly data on an annual basis. Lake bathymetry and sampling sites for seasonal 24-cycle campaigns (LAC Lake), as well as sampling sites for seasonal campaigns at wind-sheltered and -exposed sites (PAR Lake) are reported.
Plants 10 01269 g001
Figure 2. DO results from seasonal 24 h-cycle campaigns in LAC Lake. Measurements were carried out within densely vegetated areas, presenting biomass >100 gDW m−2 during samplings. Hatched color indicates night periods.
Figure 2. DO results from seasonal 24 h-cycle campaigns in LAC Lake. Measurements were carried out within densely vegetated areas, presenting biomass >100 gDW m−2 during samplings. Hatched color indicates night periods.
Plants 10 01269 g002
Figure 3. DO results from seasonal campaigns in PAR Lake. Measurements were carried out at wind-sheltered and -exposed sites, in vegetated and plant-free areas. For a better readability, Tukey’s HSD results are not reported for the seasonality factor. *** indicates p-value < 0.001.
Figure 3. DO results from seasonal campaigns in PAR Lake. Measurements were carried out at wind-sheltered and -exposed sites, in vegetated and plant-free areas. For a better readability, Tukey’s HSD results are not reported for the seasonality factor. *** indicates p-value < 0.001.
Plants 10 01269 g003
Figure 4. Dependence of DO saturation on wind exposure (Keddy Index). Dashed line is the breakpoint indicating the threshold between low-risk and high-risk of hypoxia/anoxia events in densely vegetated areas.
Figure 4. Dependence of DO saturation on wind exposure (Keddy Index). Dashed line is the breakpoint indicating the threshold between low-risk and high-risk of hypoxia/anoxia events in densely vegetated areas.
Plants 10 01269 g004
Figure 5. Hypoxia risk maps calculated for one year in densely vegetated areas of LAC Lake (on the left) and PAR Lake (on the right). Lakes bathymetry is also reported.
Figure 5. Hypoxia risk maps calculated for one year in densely vegetated areas of LAC Lake (on the left) and PAR Lake (on the right). Lakes bathymetry is also reported.
Plants 10 01269 g005
Table 1. Summarized results of the two-way ANOVA on physicochemical parameters (diel variation and season as fixed factors; sampling site as random factor). Results refer to seasonal 24 h-cycle campaigns carried out at 15 vegetated sites of LAC Lake.
Table 1. Summarized results of the two-way ANOVA on physicochemical parameters (diel variation and season as fixed factors; sampling site as random factor). Results refer to seasonal 24 h-cycle campaigns carried out at 15 vegetated sites of LAC Lake.
T (°C)pH UnitsDO (%)
Sourcedf, Residualsp-Valuedf, Residualsp-Valuedf, Residualsp-Value
Diel variation1, 204<0.0011, 204n.s.1, 383n.s.
Season2, 204<0.0012, 204<0.0012, 383<0.001
Diel × Season2, 204n.s.2, 204n.s.2, 383n.s.
CO2 (%)CH4 (µM)NH4+ (µM)NO3 (µM)
df, Residualsp-Valuedf, Residualsp-Valuedf, Residualsp-Valuedf, Residualsp-Value
Diel variation1, 347n.s.1, 415n.s.1, 410<0.051, 394n.s.
Season2, 347<0.0012, 415<0.0012, 410<0.052, 394<0.001
Diel × Season2, 347n.s.2, 415n.s.2, 410n.s.2, 394n.s.
Table 2. Summarized results of the three-way ANOVA on physicochemical parameters (wind exposure, plant presence and season as fixed factors; sampling site as random factor). Results refer to seasonal campaigns carried out at 12 vegetated and plant-free sites of PAR Lake. Please refer to figures for Tukey’HSD test differences between treatments.
Table 2. Summarized results of the three-way ANOVA on physicochemical parameters (wind exposure, plant presence and season as fixed factors; sampling site as random factor). Results refer to seasonal campaigns carried out at 12 vegetated and plant-free sites of PAR Lake. Please refer to figures for Tukey’HSD test differences between treatments.
T (°C)pH UnitsDO (%)
Sourcedf, Residualsp-Valuedf, Residualsp-Valuedf, Residualsp-Value
Wind exposure1, 8<0.0011, 8<0.0011, 8<0.001
Plant presence1, 8n.s.1, 8n.s.1, 8n.s.
Season3, 24<0.0013, 24<0.0013, 72<0.05
Wind × Plant1, 8n.s.1, 8<0.051, 8<0.05
Wind × Seas3, 24<0.0013, 24<0.0013, 72<0.001
Plant × Seas3, 24n.s.3, 24n.s.3, 72<0.05
Wind × Plant × Seas3, 24n.s.3, 24n.s.3, 72<0.001
CO2 (%)CH4 (µM)NH4+ (µM)NO3 (µM)
df, Residualsp-Valuedf, Residualsp-Valuedf, Residualsp-Valuedf, Residualsp-Value
Wind exposure1, 8<0.0011, 8<0.0011, 8n.s.1, 8<0.05
Plant presence1, 8n.s.1, 8n.s.1, 8n.s.1, 8n.s.
Season3, 72<0.0013, 72<0.0013, 72<0.0013, 72<0.001
Wind × Plant1, 8<0.051, 8n.s.1, 8n.s.1, 8<0.05
Wind × Seas3, 72<0.0013, 72<0.0013, 72<0.0013, 72<0.001
Plant × Seas3, 72<0.0013, 72<0.053, 72n.s.3, 72n.s.
Wind × Plant × Seas3, 72<0.0013, 72<0.053, 72n.s.3, 72n.s.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ribaudo, C.; Tison-Rosebery, J.; Eon, M.; Jan, G.; Bertrin, V. Wind Exposure Regulates Water Oxygenation in Densely Vegetated Shallow Lakes. Plants 2021, 10, 1269. https://doi.org/10.3390/plants10071269

AMA Style

Ribaudo C, Tison-Rosebery J, Eon M, Jan G, Bertrin V. Wind Exposure Regulates Water Oxygenation in Densely Vegetated Shallow Lakes. Plants. 2021; 10(7):1269. https://doi.org/10.3390/plants10071269

Chicago/Turabian Style

Ribaudo, Cristina, Juliette Tison-Rosebery, Mélissa Eon, Gwilherm Jan, and Vincent Bertrin. 2021. "Wind Exposure Regulates Water Oxygenation in Densely Vegetated Shallow Lakes" Plants 10, no. 7: 1269. https://doi.org/10.3390/plants10071269

APA Style

Ribaudo, C., Tison-Rosebery, J., Eon, M., Jan, G., & Bertrin, V. (2021). Wind Exposure Regulates Water Oxygenation in Densely Vegetated Shallow Lakes. Plants, 10(7), 1269. https://doi.org/10.3390/plants10071269

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop