the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Thermal stratification and meromixis in four dilute temperate zone lakes
Chris Harding
Sajjad Akam
Ioan Lascu
Gabrielle Ledesma
Pratik Poudel
Heeyeon Sun
Samuel Duncanson
Karly Bandy
Alex Branham
Liza Bryant-Tapper
Tanner Conwell
Omri Jamison
Lauren Netz
Abstract. Four adjacent lakes (Arco, Budd, Deming, and Josephine) within Itasca State Park in Minnesota, USA are reported to be meromictic in the scientific literature. However, seasonally persistent chemoclines have never been documented. We collected seasonal profiles of temperature and specific conductance and placed temperature sensor chains in two lakes for ~ 1 year to explore whether these lakes remain stratified through seasonal mixing events, and what factors contribute to their stability. The results indicate that all lakes are predominantly thermally stratified and are prone to mixing in isothermal periods during spring and fall. Despite brief, semi-annual erosion of thermal stratification, Deming Lake showed no signs of complete mixing from 2006 to 2009 and 2019–2022 and is likely meromictic. Geochemical data indicate that water in Budd Lake, the most dilute lake, is predominantly sourced from precipitation. The water in the other three lakes is calcium-magnesium bicarbonate type, reflecting a source of water that has interacted with the landscape. δ18OH2O and δ2HH2O measurements indicate the lakes are supplied by precipitation modified by evaporation. The water residence time in meromictic Deming Lake is short (100 days), yet it maintains a large reservoir of dissolved iron. Josephine, Arco, and Deming lakes sit in a valley with likely permeable sediments and may be hydrologically connected through wetlands, and recharged with shallow groundwater, as no streams are present. All four lakes develop subsurface chlorophyll maxima layers during the summer. All lakes also develop subsurface oxygen maxima that may result from oxygen trapping in the spring by rapidly developed thermoclines. Documenting the mixing status and general chemistry of these lakes enhances their utility and accessibility for future biogeochemical studies.
- Preprint
(3715 KB) - Metadata XML
-
Supplement
(12404 KB) - BibTeX
- EndNote
Elizabeth Swanner et al.
Status: open (until 09 Nov 2023)
-
RC1: 'Comment on egusphere-2023-1764', Anonymous Referee #1, 21 Sep 2023
reply
This manuscript investigates the thermal stratification and meromixis of four adjacent lakes (Arco, Budd, Deming, and Josephine) within the Itasca State Park in Minnesota, USA. Its relevance and main motivation are to evaluate if the reported lakes can be classified as meromictic. The authors want to assess if these lakes remain stratified through seasonal mixing events, and what factors contribute to their stability. Results show that all Arco, Budd, and Josephine Lakes may be holomictic or even dimictic and that Deming Lake is likely meromictic. I think this study will make a relevant contribution to the field of limnology. I recommend the publication of this manuscript after the following comments are addressed.
Major comments:
It was not easy to follow the manuscript results because the Materials and Methods section is very incomplete. Additionally, after reading the manuscript. I think that it would be relevant to understand what type of monitoring procedure would ensure a more solid conclusion regarding the lake’s classification. I think that this aspect of the study needs to be discussed.
Specific comments:
Line 39-50: “The stability of a lake against mixing is conferred by density differences between… meromixis in lakes of the temperate zone (Boehrer et al., 2017).” Why do you have this text in two different locations?
Line 120- “ The goals of this study are to 1) determine whether these four lakes are meromictic, 2) investigate the water type, sources, …to global biogeochemical cycles that may result from climate change increasing stratification in lakes.” The same with this text. Why do you have this text in two different locations? This text belongs to the introduction. In fact, the majority of the text that is included in the Material and Methods section belongs to the introduction (I´m not saying that you must include all of text in the introduction section). In this section (Material and Methods) you must describe for example the lake’s location, the location of the sampling points (profiles of temperature and specific conductance; temperature sensor chains) and the mathematical concepts considered in the analysis (e.g. Brunt-Vaisala or buoyancy frequency (N) equation; The dimensionless lake number equation). This section should also identify all the dataset’s sources. For example, you only mention the water colour datasets in the results section. In my opinion this section needs to be completely reformulated.
Line 150: Can you include the lakes sampling points location in Figure 3)?
Line 212: (Supplementary Figures 7-10). I suggest considering the same scale range in all figures.
Line 251 Please replace MAMSL with: Meters above mean sea level (MAMSL). This is the first time the acronym appears in the text.
Line 253 – I think you mean (Figures 1 and Supplementary Figure 2).
Line 212 - Figure 6. Caption. Crosses are spring or bog water data from Itasca State Park (Supplementary Table 1) This caption is correct? Table 1 shows water color in mgPt L-1
Line 275: “The Nicollet Creek spring sample lies closest to the intersection of the LEL and LMWL, whereas the Deming bog sample lies closest to the lakes but is more enriched than the lakes (Figure 6).” I suggest including these samples in Figure 6.
Line 317. “Deming Lake rapidly develops a thermocline after ice-off (Supplementary Figure 6)”. I think that the figure number is not correct.
Supplementary Figure 19. Can you please describe the meaning of the gray area?
Citation: https://doi.org/10.5194/egusphere-2023-1764-RC1 -
AC1: 'Reply on RC1', Elizabeth Swanner, 26 Sep 2023
reply
Response to major comments:
I can add text to the conclusion about best practices for monitoring procedures to determine mixing classifications.
Response to comments on lines 39-50 and lines 120-:
Regarding the two above points, it appears that when I transferred the manuscript to the template provided by the journal, I copied the Introduction twice – once into section 1 and the second time into section 2. The Materials and Methods are missing. I have pasted them below. I apologize for this oversight.
“Formation of the four study lakes occurred during the late-Wisconsin glaciation ~12,000 years ago (marine isotope stage 2; Jennings and Johnson 2011). They occupy a tunnel valley that was formed beneath the Wadena lobe of the Laurentide Ice Sheet. Following glacial retreat, the melting of stagnant ice blocks within the tunnel valley left depressions in the landscape now occupied by lakes and wetlands (Wright Jr. 1993).
Today, the four lakes investigated in this study sit in the HUC-12 watershed that sources the headwaters of the Mississippi River (U.S. Geological Survey 2017). Budd (478.6 meters above mean sea level; MAMSL) is the highest elevation, while Arco (465.8 MAMSL) and Josephine (465.4 MAMSL) lie at similar elevations (Figure 1). Deming is the lowest elevation of the lakes (464.8 MAMSL).
Lake depth measurements were collected using a Garmin Striker 4 dual-beam transducer (sonar) attached to a rowboat or canoe. Depth and GPS measurements were taken every six seconds while the boat was in motion. A Garmin GLO 2 GPS receiver and ArcGIS Collector app was used to navigate, track the boat’s course, and ensure even coverage. The shoreline of the lakes was obtained by walking along accessible areas of the shore with the Garmin GLO 2 GPS receiver, or from Lidar-derived digital elevation models. Bathymetry rasters (1 m resolution) were generated from the depth measurements in ArcGIS Pro 3.0 using a 3rd-degree Local Polynomial Interpolation. These rasters were used to calculate lake volumes and contour maps. Rasters and volume data have been deposited with the Environmental Data Initiative (Swanner et al. 2022).
Chemical, physical, and biological parameters measured on the four lakes included depth, temperature, specific conductance, salinity, turbidity, pH, oxidation-reduction potential, dissolved oxygen, photosynthetically active radiation, chlorophyll-a, and phycocyanin. Major cations, anions, and isotopes of water (δ2H-H2O and δ18O-H2O) were determined on lake water retrieved from different depths within the four lakes. Taxon-specific chlorophyll-a fluorescence was collected with a Fluoroprobe (BBE Moldaenke). The data and description of methods are available in the Environmental Data Initiative (Swanner et al. 2022). Measurements were made and samples were collected from a boat anchored within the deepest basin of each lake.
A string of temperature loggers (HOBO Water Temp Pro v2) placed at different depths were deployed into the deep areas of Deming, Arco, and Budd Lakes for one year. A conductance logger (HOBO Conductivity Logger) was added near the bottom of the strings in Arco and Budd after six months. These sensors measured temperature every thirty minutes and specific conductivity every 2 hours. The sensor string was not retrieved from Budd, as it could not be located in May 2022. Conductance measurements with a Yellow Spring Instruments ProDSS temperature/conductivity sensor on deployment and removal were used to check for drift in the HOBO conductance logger. Hourly wind speed data for the duration of sensor deployment utilized the ITCM5 (47.2400, -95.1900, 1480 feet elevation) weather station in Itasca State Park. Data was downloaded from MesoWest (https://mesowest.utah.edu/). Plots and analyses were produced in Python or R Studio (2022.07.2) using the RLakeAnalyzer package v.1.11.4.1 (Winslow et al. 2019).
Major anions (CO32-, HCO3-, Cl-, and SO42-) and cations (Na+, K+, Ca2+, Mg2+) were used to produce a Piper diagram in Geochemist’s Workbench 15.0. The concentration of cations and anions was calculated as the percentage of total cations and anions in meq L-1.
The isotopes of water (δ2HH2O and δ18OH2O) were measured on spring or seep water that had been filtered with 0.45 micron nylon syringe filters and stored at 4 °C with minimal headspace until analysis. Samples were analyzed with a Picarro L1102-i Isotopic Liquid Water Analyzer at the Stable Isotope Laboratory at Iowa State University. The analytical uncertainty and average correction factor for δ18OH2O are ± 0.05 ‰ and ± 0.30 ‰ for δ2HH2O relative to V-SMOW.
Samples for microscopy and water color were collected into amber bottles with a Van Dorn sampler from three different depths in each lake, including the SCML, if present, as determined with the YSI ProDSS. Water color was determined on water filtered through a GF-75 (Advantec) glass fiber filter (Cuthbert and del Giorgio 1992). Absorbance was measured at 440 nm and 750 nm. The absorption coefficient (g; m-1) was calculated by subtracting the absorbance at 750 nm from the absorbance at 440 nm and dividing by the path length (m):
g440 = (2.303 * A440 – A750)/(path length) (2)
A conversion was necessary to determine the color (mg Pt L-1) of the lake water:
Color = 18.216*(g440) - 0.209 (3)
Water sampled from the SCML was preserved with 1% Lugol’s solution upon returning to the laboratory. Fixed samples were settled in the dark for three to seven days.
Student reports from courses taking place over several decades at the Itasca Biological Station and Laboratories (IBSL) (Knoll and Cotner 2018), formerly the Itasca Biological Station, were acquired from the library at the University of Minnesota, Twin Cities.”
Response to comments on line 150:
We did not sample from a fixed mooring. The methods now include this statement, “Measurements were made and samples were collected from a boat anchored within the deepest basin of each lake.”
Response to comments on line 212:
Adjusting the scale to each dataset allows the trends to be visibly resolvable, and makes for easy visual comparison of one time point to another. I could keep the scales the same but would need to change the number of plots in the paper width, so it would be harder to compare. It is just a trade-off.
Response to comments on line 251:
This acronym is defined in the second paragraph of the methods.
Response to comments on line 253:
Supplementary Figure 2 is the drought record. Supplementary Figure 1 contains a cross-sectional lake profile showing the lake level referred to in this figure, and corresponds to the cross-sections identified in Figure 1.
Response to comments on Figure 6:
Thanks for catching this. It should be Supplementary Table 2.
Response to comments on line 275:
The last sentence of the caption has been modified to clarify, “The cross closest to the lake data points is the Deming bog sample, and the cross closest to the intersection of the LEL and LMWL is the Nicollet Creek sample.”
Response to comments on line 317:
This should be Supplementary Figure 5.
Response to comments on line Supplementary Figure 19:
This is the 95% confidence interval. I have added this information to the caption.
Citation: https://doi.org/10.5194/egusphere-2023-1764-AC1
-
AC1: 'Reply on RC1', Elizabeth Swanner, 26 Sep 2023
reply
Elizabeth Swanner et al.
Elizabeth Swanner et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
152 | 37 | 9 | 198 | 22 | 3 | 5 |
- HTML: 152
- PDF: 37
- XML: 9
- Total: 198
- Supplement: 22
- BibTeX: 3
- EndNote: 5
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1