the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Lake hydrochemistry and aquatic plant diversity across permafrost landscapes of Yakutia, Eastern Siberia
Abstract. Freshwater ecosystems are a major feature of the northern landscapes that are expected to experience significant future changes due to climate change and land-use alterations. In Central Yakutia, abundant lakes in topographic permafrost-thaw depressions, named ‘alaas’, define the traditional cultural landscape that is home to the indigenous Sakha people, with critical ecosystem services like freshwater supply, meadows for cattle breeding, as well as fishing and hunting grounds. In contrast, lakes in the Verkhoyansk mountain region east of Central Yakutia are of glacial origin or developed on glacial moraines and represent deeper and more oligotrophic lake systems much less used as human resources.
Here, we analyse the hydrochemistry, sedimentary DNA (sedDNA)-derived aquatic plant diversity, geomorphology, and adjacent land cover of sixty-six lakes across the Central Yakutian lowland permafrost landscape and the Verkhoyansk Oymyakon high mountain plateau to understand their characteristics and environmental drivers. Our hydrochemical analysis reveals a clear distinction between the low-mineralised mountain lakes and the highly variable hydrochemistry of the lowland thermokarst lakes. The lake developmental stage within the thermokarst lake sequence seems to be a key driver of lake hydrochemistry in the lakes of the Central Yakutian lowland. Specifically, the lake’s developmental stage is reflected by dissolved organic carbon (DOC), pH, its stable isotopic composition, and the hydrochemical facies of alkali and earth alkali elements. New thermokarst lakes have a depleted stable isotopic composition, possibly due to contributions from meltwater of adjacent permafrost ground-ice. This thermokarst lake stage is typically located within forest and has the highest DOC. In contrast, the hydrologically open thermokarst lake systems, typically located in large connected alaas systems with settlements and managed land use, have lower DOC and fewer mineralisation than recently formed thermokarst lakes or old alaas lakes. The dilution in the hydrologically connected alaas lakes occurs due to flushing, mainly during high discharge events such as the regular snowmelt. Old alaas lakes show an enriched oxygen isotope composition and have high salinity and mineral content, suggesting processes of evaporation and highlighting their vulnerability to future warming. However, low chloride together with an enriched isotopic composition and elevated fluoride characterise several of the sampled high-salinity lakes. This points to an additional process beyond the current evaporation, such as fluoride leakage from lacustrine sediments or salt deposits.
SedDNA-derived macrophyte diversity reflects lake types and reveals the dominance of brackish water-tolerant cosmopolitan submerged macrophytes, particularly Stuckenia and Potamogeton, across all lake types. The macrophytes Myriophyllum and M. verticillatum are exclusively found in freshwater lakes in the lowlands and the mountain regions, supporting their indicator value for freshwater conditions. Our results provide a detailed examination of lake systems in modern conditions within highly climate-sensitive lowland and mountain permafrost landscapes.
- Preprint
(4399 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
- RC1: 'Comment on egusphere-2025-4114', Anonymous Referee #1, 07 Oct 2025
-
RC2: 'Comment on egusphere-2025-4114', Anonymous Referee #2, 04 May 2026
This is an interesting paper on lake hydrochemistry of two contrasting sites in Republic of Sakha. Authors examine drivers of lake inorganic and organic chemistry by considering landscape and climate characteristics. The idea is clear; the research is scientifically sound and data rich. Lakes are discussed from different perspectives, but at the same time, in a quite complex way. Here are several suggestions on how to improve the quality of the manuscript.
Authors are advised to improve M&M, Results, and Discussion in a way that the narrative better connects RQ stated, Results obtained, and the Discussion of these results.
This can be done by 1) better explaining reasons for assigning thermokarst development stage to lowland lakes and the choice of lakes for sampling; 2) better explaining the rationale of choosing the parameters involved (buffer zones, LC dataset, and climatic variables; 3) better explaining the reason of MPA.
Further, although the laboratory analytics and data post-processing and Methods seem very convincing, the Discussion of the obtained Results is a little speculative.
Abstract requires some numbers and clearer statements about drivers: higher/lower seams non-informative; hydrochemical analysis clearly distinguished mountain lakes from lowland lakes: mountain lake hydrochemistry is determined by the climate and topography; lowland thermokarst lakes are diverse: this is explained by the lake development stage, these groups show different ranges of hydrochemical characteristics: DOC, pH etc.
Please explain how sampling sites were chosen: sample sizes for each type are inequal (15-15-8-28). Did Authors consider environmental drivers (LULC, climate) while choosing sampling sites? How did Authors account for geomorphology? Through lake shore elevation?
L30-32 – this sentence is purely formulated: Do you intend to understand environmental drivers of lake hydrochemistry by land cover, plant diversity, land use, and climate? Or what? It must be rephrased.
L46-50 – this part must be introduced into the main narrative (L35-45) about environmental drivers. Otherwise, it stays aside of the main results.
Since the impact of lake type is among the main findings of a paper, the characteristic features of each type (lowland thermokarst lakes) must clearly be specified. How did Authors differentiate visually between tympy lake and mature alaas lake, that are of different stages of development? Do we connect these types to ages?
Alaas lakes, alaas landscapes, lake systems (L96) – all are very misleading terminologies.
L99 - CAL – these might be connected during one season and “disconnected” during another. More robust explanation of CAL is required.
L128 – thermokarst/alaas development stage?
L178-192 – do Authors present maximal lake depth, or depth of sampling? How did Authors use transparency (Secchi disk), and chromophoric DOM in the analysis? (these parameters are introduced)
L276-279 - how did Authors use lake morphology in the analysis? (round, oval-elongated etc)
In general, sub-Section 2.4 describes the UAV data processing and GIS analysis. A little is described regarding the lake type classification (lowland thermokarst lakes)
F.e. Fig. 2D to my subjective understanding may be considered as NTL, CTL, and OTL (depending on the position in the alas)
Sub-Section 2.5. I do not understand why Authors have chosen to additionally utilize 10-m ESA WC land cover, if drone images are available?
At the same time, important topographic catchment characteristics have not been included (i.e. from ArcticDEM).
Please explain the rationale of choosing 50- and 360-m buffers for extracting LC fractions. Why not 100 or 200 m?
L312-313, L315-316 What is the rationale of merging grasslands with croplands from WC for Central Yakutia, if Authors intend to analyze the impact of the land use? Same with Built-Up class
L330-337 Why not to choose data from nearest meteorological stations? What about permafrost relevant parameters: ground temperature, ALT?
L369-370 This is a repetition
L384 – ML TBS is 49 mg per L, but Table 1 shows the value of 61. Similar to OAL.
L414-419 Authors write that in 60 lakes, alkaline earths exceed alkalis, but at the same time, an opposite effect is observed in 9 lakes. How could that be?
What is the idea of the lake clustering? How is answers the RQs stated at the beginning? (Intro)
L512 Forest land cover fraction?
L555 – what is the Elev variable? Ellipsoid elevation of the shoreline? Or mean catchment elevation?
L575-580 This is the repetition. I suggest removing this part.
L579 what is the “hydrochemical lake status”?
L581-584 It is obvious that permafrost lakes receive an input from surface runoff, but Authors do not directly measure the hydrochemical parameters of the active layer or catchment streams.
L594-595 One could add stable water isotope values for ground ice and add these numbers to the Figure 5 with isotopic composition with relevant references.
Section 4.3 In the title, Authors declare “thermokarst lakes”, but it is obvious that these are of different origin
L650 – it is not a glacial type, it is a glacial origin
L656-659 – this is a repetition as in 4.1
L660-667 how this paragraph is connected to findings of this paper?
Fig. 8 – orange and yellow colors can hardly be distinguished
4.5. Yakutian lakes as a resource – L731-751 (this is an irrelevant information). I suggest removing this part.
L1324 – forming
Regarding the speculations about sources of fluoride in samples considered: Authors are advised referring to a research paper by Pavlova & Fedorova 2022 “Fluoride distribution in sub-permafrost groundwater in Central Yakutia” published in Earth’s Cryosphere.
The impact of thermokarst development stage on lake hydrochemistry has already been discussed in Manasypov et al. 2014 (‘Thermokarst lake waters across the permafrost zones of western Siberia’)
Overall, the Discussion is the weakest part of the manuscript, that must largely be improved by replacing the irrelevant parts with truly important findings of other researchers.
Citation: https://doi.org/10.5194/egusphere-2025-4114-RC2
Data sets
Summer anorganic hydrochemistry (cations and anions) of thermokarst lakes in the Central Yakutian Lowland and mountain lakes in the Verkhoyansk Mountain Range in Eastern Yakutia [dataset]. I. Baisheva et al. https://doi.org/10.1594/PANGAEA.959459
Stable water isotope data of 66 lakes from a summer field campaign in Central and Eastern Yakutia, Siberia in 2021 (RU-Land_2021_Yakutia) [dataset] A. Stieg et al. https://doi.org/10.1594/PANGAEA.950688
53 orthomosaics processed from Unoccupied Aerial Vehicle (UAV) image data of lake shorelines sampled during a field campaign in Central and Eastern Yakutia, Siberia in 2021 (RU-Land_2021_Yakutia) A. Stieg et al. https://doi.org/10.1594/PANGAEA.956223
Lake information on orthomosaics created from lake image data sampled during a field campaign in Central and Eastern Yakutia, Siberia in 2021 (RU-Land_2021_Yakutia) [dataset]. A. Stieg et al. https://doi.org/10.1594/PANGAEA.955723
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 11,958 | 1,831 | 230 | 14,019 | 234 | 220 |
- HTML: 11,958
- PDF: 1,831
- XML: 230
- Total: 14,019
- BibTeX: 234
- EndNote: 220
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
In this study, the authors present a comprehensive comparative analysis of 66 lakes across Central Yakutia and the Verkhoyansk–Oymyakon mountain plateau to investigate how geomorphological setting, hydrochemistry, and sedimentary DNA (sedDNA)-derived macrophyte diversity reflect lake development stages and environmental drivers. The study distinguishes low-mineralised glacial mountain lakes from highly variable thermokarst lakes in the Central Yakutian lowlands and demonstrates that lake developmental stage is a key factor controlling hydrochemical characteristics such as DOC concentration, pH, isotopic composition, and ionic facies. The authors further show that newly formed thermokarst lakes influenced by permafrost meltwater exhibit depleted isotopic signals and high DOC, while older alaas lakes display increased salinity and enriched isotopic compositions due to evaporation processes. SedDNA results reveal clear macrophyte community differentiation among lake types, with freshwater taxa serving as indicators of specific hydrological and chemical conditions. Overall, the study provides valuable insights into the functioning and vulnerability of permafrost-affected freshwater systems under ongoing climatic and land-use change. This study is very well structured, and its content is highly relevant to the scientific advancement of this field. I thoroughly enjoyed reading it. I only have a few minor comments/suggestions.
Line 108 to 109. I suggest modifying this sentence as follows: “Inorganic and organic hydrochemistry (e.g., concentrations of major and minor ions and dissolved organic carbon (DOC)) can reflect natural or anthropogenic processes in lakes such as land cover and human practices (Hu et al., 2020; Du et al., 2023).”
Line 108 to 109. Focusing on DOC is appropriate for these lakes, but the authors should clarify why POC was not included.
Line 173 I suggest including additional lake characteristics in Table A1, such as area, length, width, mean depth, maximum depth, lake shape.
Line 186 “The water samples were collected from the upper water layer at a sampling depth of 0.5 m into 2 L sample containers”. Could the authors clarify why, for the deeper lakes, samples were not also taken from the middle or bottom layers to capture potential vertical variability in hydrochemistry?
Line 215-216. I believe the isotope notation in this sentence is incorrect. Please replace it with: The compositions of hydrogen and oxygen isotopes are represented using the δ notation of δD and δ18O in per mille (‰) versus”.
Line 286 Legend of Figure 2: I believe the lake numbers should be corrected. “In panels (A, D), two examples of new thermokarst lakes are EN21448 and EN21418 (purple outlines).”
Line 1187, Table A3: Could the authors clarify why P and PO4 values are not reported for Mountain lakes and Deep Mountain lakes?