Suzdalevo Lake (Central Siberia, Russia)—A Tunguska Event-Related Impact Crater?

In 1908, a massive explosion known as the Tunguska Event (TE) occurred in Central Siberia. However, its origin remains widely discussed and environmental impacts are not known in detail. We investigated evidence of the TE in sediments of Suzdalevo Lake, which is located near the explosion epicenter. According to local nomads (Evenkis), Suzdalevo Lake did not exist before the TE and was considered as a possible impact-origin water body. However, apart from oral testimony, there is no evidence of the lake formation process. Two short sediment cores (SUZ1 and SUZ3) were retrieved from the lake and dated using 210Pb and 137Cs. The sedimentary record was characterized using magnetic susceptibility, X-ray fluorescence, and the screening for melted magnetic microspherules. To study possible effects of the TE on the lake ecosystem, we performed diatom and freshwater fauna remains analyses. Results indicate that the lake contains sediments that originated before the TE and thus its formation was not related to the impact. Also, the depth to diameter ratio of the lake basin is too low (<1/100) for a young impact crater. In one of the two cores (SUZ1), we documented distinct changes in the lake-catchment ecosystem that occurred within a 5-cm-thick depth interval calculated for the best fit depths for the year 1908 using three alternative age-depth models (CRS, CIC, CFCS), namely, increases in terrestrial matter input (abundant fine plant macroremains, peaks in magnetic susceptibility and the Sr to Rb ratio) and taxonomic diversity and relative abundance of benthic taxa. The shifts in aquatic biota assemblages were likely caused by nutrient supply and improved water column mixing following a catchment disturbance. Nevertheless, precise timing of the observed abrupt changes in relation to the TE is not clear due to uncertainty of the 210Pb dating method and absence of melted magnetic microspherules or an event layer. The disturbance signals in the proxy data may postdate the TE. Our results demonstrate potential usefulness of the paleolimnological approach to understand the possible environmental consequences of the TE and similar events elsewhere.


INTRODUCTION
Impact and impact-like events are considered potential triggers of abrupt climatic and environmental changes, as their direct destructive power influences atmospheric and geochemical processes on local as well as global levels (Covey et al., 1994). Their potential consequences include not only earthquakes, tsunami, and fires, but also processes such as dust loading, ozone depletion, sulfate aerosol formation, disturbances of Earth's magnetic field, and nitric acid rain (Turco et al., 1981;Kolesnikov et al., 1997;Toon et al., 1997;Shishkin 2007;Wünnemann and Weiss 2015;Artemieva and Shuvalov 2016). On the morning of 30th June 1908, a massive explosion occurred in Central Siberia, near the Podkamennaya Tunguska River, in the Evenkiysky District, Russia (e.g., Gladysheva, 2020a;Gladysheva, 2020b and references therein). The explosion affected an area of about 2,000 km 2 of continental taiga and uprooted more than 80 million trees with an estimated power equivalent to 20-30 million tons (Mt) of TNT (Rosanna et al., 2015;Artemieva and Shuvalov 2016;Robertson and Mathias 2019). It was likely the most destructive impact-like event witnessed in the modern history of humankind. Fortunately, the Tunguska Event (TE) epicenter was located in a very remote and sparsely populated region. At least three people (local Evenki nomads) died as a direct consequence of the TE, and an unknown number of people were injured from the effects of the shockwave, thermal radiation, or glass damage (Jenniskens et al., 2019). A much smaller explosion (~0.5 Mt of TNT) from a meteorite about 19 m in diameter caused injuries of~1,600 people on 15th February 2013 in a more populated area in the Chelyabinsk Oblast, Russia (Brown et al., 2013;Kletetschka et al., 2015;Kartashova et al., 2018), documenting the destructive potential of such extreme events. TE-related glass damage was reported over a wide area, especially around the Angara River, even 400 km away from the explosion epicenter (Jenniskens et al., 2019). The night sky shone brightly for a few days over an even larger area, including Europe. The origin of this phenomenon was explained by the same mechanism as the formation of noctilucent clouds, with water and dust creating ice crystals at heights from 70 to 300 km in the atmosphere. The reflective surface of these TE noctilucent clouds was 10 4 × greater than usual, probably caused by a two to three orders of magnitude larger crystal size (Gladysheva 2012).
The TE is considered to have been caused by the impact of a cosmic body, but neither an impact crater nor any stony fragments have been found (Artemieva and Shuvalov 2016). An exotic quarzitic boulder, known as John's Stone, was hypothesized by some authors as a potential meteorite, possibly of Martian origin, but research by Bonatti et al. (2015) suggested its relation to the local Permian-Triassic Trap magmatism-hydrothermalism. In recent years, several lines of evidence have been published in support of the hypothesis that Cheko Lake (~10 km NW of the inferred TE epicenter) fills an impact crater created by a fragment of the parent TE cosmic body. These include geophysical data (e.g., Gasperini et al., 2007;Gasperini 2015), the presence of tree trunks and branches below lake sediments (Gasperini et al., 2014), and physical modelling (Foschini et al., 2019). However, this hypothesis was questioned by Collins et al. (2008) and Rogozin et al. (2017) based on theoretical models and sediment dating, respectively. According to local Evenki reindeer-herding nomads (Evgeniya Karnouchova, pers. comm.), a TE-related origin is also assumed for Suzdalevo Lake, a small water body located in the southern part of the tree fall area, as it was first mapped by expeditions after the TE (Vasilyev 1998;Jenniskens et al., 2019, and references therein). Today, more than 110 years after the explosion, many doubts and uncertainties persist about the TE. As the original asteroid hypothesis was called into question due to the absence of a meteorite fragment, a cometary origin of TE has been suggested by some authors (e.g., Kresák 1978;Gladysheva, 2020a;Gladysheva, 2020b). However, other causes (e.g., the volcanic ejection of natural gas, a dark matter collision, solar activity) have been discussed as well (Kundt 2001;Froggatt and Nielsen 2015;German 2019). In a recent paper, Khrennikov et al. (2020) argue that the TE was caused by an iron asteroid body, which passed through the atmosphere at a minimum altitude of 10-15 km with trajectory length about 3,000 km and continued to the near-solar orbit.
Overall knowledge of the TE is limited due to the lack of scientific observations at the time of the explosion and almost no research on the epicenter performed in the following two decades. Moreover, most of the later investigations were conducted to determine what phenomenon caused the TE, with much less attention paid to the topic of environmental damage. Natural archives, such as tree rings, peat, and especially lake sediments, have the great potential to provide important new evidence of this kind.
Only a few studies on these natural archives in the Tunguska region have been performed so far. Tree-ring research by Kletetschka et al. (2017) revealed a directional chemical response in the xylem of larch (Larix) trees that survived the TE inside the tree fall area. Badyukov et al. (2011) described the presence of metal microspherules consisting of Ni(Cr) bearing wüstite and magnetite with Ni-rich metal inclusions and glassy silicate microspherules in soils above the floodplain of the Chunya River, which may indicate the possible presence of cosmic body material in various sedimentary natural archives. Similar findings of microspherules were reported much earlier by Zaslavskaya et al. (1964). Lake sediments have been rarely studied, as the location is very remote, and lakes are not common in the area. An important contribution was provided by investigations of the Zapovednoe Lake sedimentary record. Despite its location outside of the tree fall area and the radiant burn area (60 km from the TE epicenter; Johnston and Stern 2019 and references therein), the lake sediment contained a characteristic yellowish (clayey) TE layer possibly caused by an earthquake and likely related catchment erosion (Kletetschka et al., 2019a;Darin et al., 2020). A considerable increase in magnetic susceptibility followed immediately after the deposition of this TE layer, pointing to intense fires near the TE epicenter. These fires (e.g., Svetsov 2002;Jenniskens et al., 2019) could converted non-magnetic iron in organic material into magnetite particles with high magnetic susceptibility that were then incorporated into ash and deposited in soils and various sediments (Kletetschka and Banerjee 1995;Kletetschka et al., 2019a). Therefore, such a characteristic signal should also be present in other lake sediment records and peat sequences in the region. A similar clayey layer of TE age has already been found in Cheko Lake sediments . If combined with a paleoecological approach, the presence of characteristic TE layers should enable studies on the effects of the event on related ecosystems. The potential of such paleoecological methods has already been demonstrated by Tositti et al. (2006) who discovered TE-related anomalies in tree pollen productivity and the onset of a phase of more humid conditions at the Raketka peat bog site (8 km from the TE epicenter). Other peat profiles (sites located 6-65 km from the TE epicenter; plant material of Sphagnum fuscum) have revealed probable traces of acid rains as a reaction to the TE, as evidenced by a positive anomaly in δ 15 N (Kolesnikov et al., 1998).
In this study, we present a paleoenvironmental record of Suzdalevo Lake sediments. As Suzdalevo Lake has been reported by local Evenki people to have appeared just after the TE, we assumed that the lake could be a TE impact crater and/or an important natural archive to investigate the disturbance of the taiga forest and related changes in the aquatic environment. Reconstructing local environmental changes caused by the TE is crucial to understanding the processes triggered by this unique phenomenon.

REGIONAL GEOLOGICAL SETTING AND THE STUDY SITE
The TE tree fall area of butterfly shape (axis along the 115°a zimuth; Figure 1A) is located in the Tunguska Basin, which is in the western part of the Siberian craton that is composed of highgrade Archeanic metamorphic complexes and granites (Gladkochub et al., 2006;Artemieva and Shuvalov 2016). At the Permian-Triassic boundary (~251 Ma), the area was affected by the extensive eruption of Siberian flood-volcanic rocks, the largest subaerial volcanic event known (~7 × 10 6 km 2 buried by basalts) (Kamo et al., 2003). Recent local river basins are filled with Quaternary sediments (fluvial and aeolian deposits) and soils, which also applies to the surroundings of the study site. Permafrost is discontinuous in this area with thickness up to 25 m (Gasperini et al., 2007;Ponomarev et al., 2019). Distinct climatic changes have been observed in overall central Siberia over the last century. The climate has become much warmer, and in some regions dry summer conditions have prolonged the "no-rain" period. This increases the likelihood of extremely dry periods and the risk of extreme fire events (Tchebakova et al., 2011). At the nearest meteorological station, that is placed in Vanavara (60°20′45″N 102°16′54″E; Figure 1A), the current (period 2000-2020) mean annual temperature is −4.5°C, the mean January temperature −27.7°C, and the mean July temperature 17.8°C. The "normal" mean annual precipitation (period 1933-2020) is 421 mm, however, the minimum mean annual precipitation of 271 mm was observed in 2016 (Weather and Climate 2021).
Suzdalevo Lake is a shallow water body located very close to the Chamba River (60°39′29.05″N, 102°3′3.36″E) a tributary of the Podkamennaya Tunguska River, 20 km SW from the TE epicenter (i.e., inside the tree fall area; Figure 1A). The lake is separated from the Chamba River by a natural dike ≥45 m wide and ≥4 m high. No bathymetric measurements or limnological characterizations of this site have been previously published. The lake was allegedly named after K. I. Suzdalev, a merchant from Vanavara, who visited the site shortly after the TE (Evgeniya Karnouchova, pers. comm.). Vegetation in the lake catchment is characterized by mixed taiga consisting mostly of pine (Pinus), birch (Betula), spruce (Picea), willow (Salix), and Dahurian larch (Larix gmelini) trees, with undergrowth of blueberries (Vaccinium myrtillus), cranberries (Vaccinium vitis-idaea), mosses and lichens. No historical logging at this site is known.

Morphobathymetric Characterization of the Study Site
In May 2019, we surveyed the lake with a close-spaced grid of echosounding lines, with the aim to compile a morphobathymetric map of its floor and identify a suitable position for sediment coring. We used techniques and instruments like those employed by Gasperini et al. (2007) at Cheko Lake. The bathymetric measurements were performed with a 200 kHz echosounder tied to a wooden stick aside of a boat while a transducer was located at about 30 cm below the water surface. Instead of acquiring only the depth value provided by the echosounder through the serial port, we collected the entire echogram by means of a 16 bit A/D converter driven by SwanPRO acquisition software. The acquisition parameters of the echosounder signal were as follows: sample rate 0.856 µs, trace length 25 ms, and 2,140 trace samples (seven traces per second). The collected data were stored in XTF format, converted to SEGY format, consistency was checked using the Segy-change preprocessing software program (Stanghellini and Carrara 2017), and processed using SeisPrho software program (Gasperini and Stanghellini 2009) to display echograms and semi-automatically digitize the sediment-water interface. The lake-floor profiles were subsequently corrected for the vertical offset between the water surface and the transducer vertical position, and then exported in ASCII (lon, lat, z) format. The file was then processed with GMT to produce both a regular grid file and the final map. The original data were affected by spatial noise due to GPS precision being limited to about 1 m. To reduce this spatial noise, the final morphobathymetric map was compiled using a 2D median filter with a searching radius of 10 m.

Core Retrieval, Sediment Lithology
During the same field survey, two short lake sediment cores, SUZ1 and SUZ3, were collected using a Kajak gravity corer (sampling tube diameter of 5.8 cm and length of 50 cm) from the deepest middle part of Suzdalevo Lake. The lengths of the obtained cores SUZ1 and SUZ3 were 42 and 46 cm, respectively. Both cores were immediately sliced in steps of 1 cm, with the uppermost 1 cm (0-1 cm) divided into two layers-0-0.5 cm and 0.5-1 cm. The color, grain size, and amount of plant macroremains (qualitative data) were recorded during slicing and during analyses of aquatic biota remains.

Short-Lived Isotope Dating
Both sediment cores were measured using gamma spectroscopy for the specific activity of 210 Pb, 137 Cs, and 226 Ra isotopes in the Radiometry Laboratory (Institute of Geochemistry, Mineralogy and Mineral Resources Charles University, Prague, Czechia). Individual samples (1-cm-thick layers; the two uppermost layer were merged) were measured by well-in-well geometry in a SILAR ® low-background anti-Compton-anti-coincidence gamma spectrometer with a specially designed 40 × 40 mm Na(Tl) well-type LEADMETER ® detector with a total efficiency of 46.2% (for 210 Pb line of 47 keV) placed in a well-type guard NaI(Tl) detector 160 × 125 mm in 10-cm-low-background lead shielding (Hamrová et al., 2010). A Canberra DSA 2000 multichannel analyzer controlled by GENIE 2000 software was used to determine the specific activities of the abovementioned isotopes. The measuring time for individual samples was 2 days, and 6 days for background. Special inhouse standards with a light matrix were used for radionuclide quantification in the same geometry of the 8 ml vials. The IAEA-447 standard (moss-soil) was used as 210 Pb reference material with a result of 340 ± 6 Bq kg −1 for the recommended value of 338 Bq kg −1 and corrected for decay from the reference date (15/11/2009). Due to low weights of the SUZ1 samples affecting the results, we decided to measure isotope activities also for the SUZ3 core by merging four adjacent 1-cm-thick layers in a row. The output from the heavier 4-cm-merged samples (~4 g of dry sediment) was expected to be more accurate. The age of individual layers was obtained using the Constant Rate of Supply (CRS), the Constant Initial Concentration (CIC), and the Constant Flux Constant Sedimentation rate (CFCS) agedepth models that are implemented in an R package "serac" (Bruel and Sabatier 2020 and references therein). These alternative models were used because none of them can be considered as a universal method of age-depth model construction.

Sediment Magnetism and Geochemistry
All sediment samples from both cores were placed into plastic cups and magnetic susceptibility was measured using a SM30 Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 777631 magnetic susceptibility meter (ZHinstruments Inc., Czechia) at an oscillation frequency of 8 kHz and generated magnetic field amplitude of 40 A/m. The magnetic susceptibility values [X (SI), average of three measurements] were mass normalized. The total content of selected elements was determined by means of X-ray fluorescence (XRF) spectrometry using a handheld ED-XRF analyzer VANTA VMR with a Silicon Drift Detector (Olympus, United States). The ED-XRF analyzer was coupled with a programmable moving core holder and autonomously run by a computer (for details Kletetschka et al., 2018;Kletetschka et al., 2019b). National Institute of Standards and Technology (NIST) standard reference materials 2711a Montana II Soil and 2710a Montana I Soil were used for quality control. Finally, we focused on an interpretation of concentrations and ratios of individual elements that are likely relevant to tracking the history of the lake's productivity as well as the erosional activity in the lake catchment. The strontium to rubidium ratio (Sr/Rb) and titanium (Ti) concentration were selected as proxies characterizing the input of material from the catchment (Cuven et al., 2010;Xu et al., 2010), phosphorus (P) and the silicon to zirconium ratio (Si/Zr) as proxies of lake productivity (Cuven et al., 2011), and the iron to manganese ratio (Fe/Mn) with iron (Fe) concentration control as proxy of near bottom redox conditions (Mackereth 1966;Davison 1993).

Separation and Analyses of Magnetic Microspherules
For separation of magnetic microspherules from the sediment samples (all SUZ1 layers), we performed a standard separation technique described by Israde-Alcántara et al. (2012). The separated material was analyzed under a dissecting optical microscope at 25-70× magnification. Suspected spheric objects (diameter >3 µm) were manually placed on aluminum stubs with nonmagnetic tools (sharpened wooden sticks). Finally, the spheric objects were identified with a scanning electron microscope (SEM) TESCAN Vega in back-scattered electron (BSE) and secondary electron (SE) modes, and subsequent elemental microanalysis was conducted using an energy dispersive X-ray spectroscope (EDS, detector X-Max 50; Oxford Instruments, United Kingdom) at the Laboratory of Scanning Electron Microscopy (Institute of Petrology and Structural Geology, Faculty of Science, Charles University, Prague, Czechia).

Diatoms and Freshwater Fauna Remains
Diatom analysis and the analysis of freshwater fauna remains were conducted at 1 cm resolution throughout the SUZ1 sediment core. Diatom samples were prepared following the method described by van der Werff (1955). The material was cleaned by adding 37% H 2 O 2 and heating to 80°C for about 1 h, then by the addition of KMnO 4 . Following digestion and centrifugation, the resulting clean material was diluted with distilled water to avoid excessive concentrations of diatom valves, which may hinder reliable observations. Known quantities of Lycopodium spores were added to estimate diatom concentrations.
Cleaned diatom valves were mounted in Naphrax ® , a highrefractive index medium. In each sample, 400 diatom valves were identified and enumerated on random transects at 1,000× magnification (oil-immersion), using an Olympus BX53 microscope equipped with Differential Interference Contrast (Nomarski) optics and an Olympus UC30 Imaging System. Diatoms were categorized into five ecological groups according to their life form (modified according to Buczkó et al., 2013): 1) aerophyton (living in subaerial and terrestrial habitats), 2) benthos (living at the bottom and near shore, mainly in soft sediments), 3) periphyton (attached to submerged substrata or organisms; more sessile than benthic diatoms), 4) plankton (living in the water column), and 5) tychoplankton (taxa frequently encountered both in the water column and in surface sediments). In addition, the category of benthos was further subdivided into three subcategories: epipsamon (living attached to sand grains), benthos sensu stricto (lake benthos besides epipsamon), and river benthos. Furthermore, the D to C ratio of diatom valves to chrysophyte cysts was calculated. The chrysophyte cyst sums are expressed relative to the number of diatom frustules (two valves = one frustule) plus chrysophyte cysts, with the following formula: D:C = [number of diatom frustules/(number of chrysophyte cysts + number of diatom frustules)] × 100 (Douglas and Smol 1995).
Samples for zoological indicators were washed with distilled water through a 100 µm sieve and transferred into a modified plastic Sedgewick-Rafter counting cell. All identified animal remains, mainly chironomid (non-biting midge) head capsules (HCs), were picked with either fine forceps or a steel needle using a dissecting optical microscope (15× magnification). Samples were then dehydrated in ethyl alcohol (~96%) and mounted on glass microscope slides in Euparal ® mounting medium. For taxa identification, we used an optical microscope (100-400× magnification) and followed several identification keys-mainly Wiederholm (1983), Rieradevall and Brooks (2001), Brooks et al. (2007), and Szeroczyńska and Sarmaja-Korjonen (2007). Due to low abundances of all zoological indicators, including chironomid HCs, in most of the samples, we show only concentrations of individual taxa (i.e., fossils per 1 g of dry sediment).
Statistically significant changes in the diatom record were determined using the optimum sum of squares partitioning on percentage data implementing the ZONE program (Birks and Gordon 1985;Lotter and Juggins 1991). The results from optimum sum of squares partitioning on percentage data were tested for significance using a broken-stick model to derive a stratigraphical zonation (Bennett 1996). Due to the presence of zero values in the chironomid record in some of the analyzed layers (which do not allow independent zonation), we use the diatom zonation when the zoological indicator data and other proxies are discussed. Figure 1C displays the morphobathymetric map of Suzdalevo Lake obtained after data processing and compilation, while Supplementary Figure S2 shows echographic (depth) profiles with the depth magnified by a factor of 10. The lake is very shallow, with a maximum depth of about 2.3 m and a surface area of~1.8 ha. The depth-to-diameter ratio is about 0.015, and its irregular inverted-cone morphology is slightly elongated in the N-S direction. The depocenter is located close to the N shore, while the E and W shores show morphological irregularities. Both sediment cores (SUZ1 and SUZ3) were recovered from the depocenter ( Figure 1C) and consisted of homogeneous dark brown gyttja with no distinct lamination or a clayey layer.

Sediment Age Determination
In the SUZ1 core, 210 Pb, 226 Ra, and 137 Cs activities were measured only for the upper 26 cm. Due to their low values and large error bars (Supplementary Figure S3, Supplementary Table S1), the same measurements were not performed for the deeper layers. Excess ("unsupported") 210 Pb ( 210 Pb ex ) was detected even in the deepest analyzed sample (depth of 25-26 cm), therefore the 210 Pb ex inventory could not be calculated. The 137 Cs specific activities are relatively low and irregular at depths below 16 cm. The maximum 137 Cs activity of 93 Bq kg −1 was detected at the depth of 9.5 cm. Weights of some SUZ1 samples were <1 g of dry sediment (average 0.9 g).
The problem of low sample weights in the SUZ1 core was resolved using the parallel core SUZ3, where four adjacent layers were merged to produce larger samples (3.1-5.6 g of dry sediment). In nine merged (4-cm-thick) layers, the specific activities of 210 Pb decrease exponentially from the topmost layers to background ("supported") activities in the lower part of the core, showing a relatively smooth decreasing curve (Supplementary Figure S3). 210 Pb ex was successfully detected in upper eight merged samples (Figure 2, Supplementary Table  S1). The trend in the SUZ3 relation between 210 Pb ex and depth in core is similar to the SUZ1 record, but the curve is smoother and shows lower error bars (Supplementary Figure S4). The 137 Cs specific activities are relatively low (maximum activity of 61 Bq kg −1 ), demonstrating low contamination by radioactive fallout. An increase in 137 Cs activity, that should be equivalent with the global fallout due to nuclear weapon tests before its maximum in 1963 CE, was found between depths of 14 and 18 cm (  Table S2). However, this 137 Cs peak is flat and low 137 Cs activities are dispersed in all samples of the upper part of the core, probably due to continuous flushing out from the surrounding contaminated area, occasional sediment-water interface mixing (redeposition), and younger atmospheric fallout (later nuclear tests and power plant accidents).
Calculated sediment accumulation rates for the SUZ3 core differ among the age-depth models (Supplementary Table S2).
Based on the CRS model, the accumulation rates decrease FIGURE 2 | Short-lived radionuclide activities and age-depth models for the SUZ3 core. From left to right: 210 Pb ex activity (semilogarithmic plot), 137

The Magnetic and Geochemical Record
Mass normalized magnetic susceptibility values [X (SI)] in the SUZ1 core show a peak at the depth of 27.5 cm and an increasing trend between the depths 15 and 0 cm ( Figure 3). For the same core, the concentrations of selected elements (ppm) and their ratios are plotted with depth and the diatom zonation (see details in Section 3.5 and Figure 3) to better interpret the lake sedimentation history. Two proxies characterizing catchment material input to the lake are used. The Sr to Rb ratio has a similar trend as magnetic susceptibility (a peak at the depth of 26.5 cm and an increasing trend from 16 cm), except for a drop in the two uppermost layers that correspond to an increased concentration of Ti. Ti concentrations are otherwise relatively stable in the profile. Both proxies of primary productivity in the lake environment, P concentration and the Si to Zr ratio, show peaks at 37.5 cm, low values around 26.5 cm, and increased values between 24 and 12 cm (max. P concentration of~2,250 ppm and max. Si/Zr of~5,800) followed by low values (min. P concentration of~1,070 ppm and max. Si/Zr of~1850). Trends of the Fe to Mn ratio curve correspond well with Fe concentrations, showing increased values around 27.5 cm and a sharp increase with maxima in the two uppermost layers (shifts from 85 to 133 for Fe/Mn and from~50,500 to~137,200 ppm for the Fe concentration). The Fe to Mn ratio, a potential proxy for changes in the near-bottom redox conditions, was therefore determined by fluctuations in the Fe concentration. A similar magnetic susceptibility and geochemical record was found for the SUZ3 core ( Figure 3). Nevertheless, it shows two noticeable differences: 1) absence of the magnetic susceptibility peak at the depth of 27.5 cm followed by the Sr/Rb peak at the depth of 26.5 cm, and 2) absence of the magnetic susceptibility, Ti, Fe/Mn, and Fe peaks in the two uppermost layers. Based on this comparison, the SUZ1 core provided more dynamic records and was selected for extraction of iron-rich microspherules and analyses of biological indicator remains.

Iron-Rich Microspherule Presence
During the SEM observation of the magnetic extract from the Suzdalevo Lake sediment, we found reduced (oxygen content <50%) titanomagnetite (TM) grains in the depth interval of 26-30 cm (Supplementary Figure S5). However, their concentration is very low, not exceeding three grains per 1 g of dry sediment. Only one magnetic microspherule larger than 3 µm in diameter was found in the same samples with an original wet volume of 28 cm 3 (14 cm 3 for the two youngest layers). This iron oxide microspherule (10 µm in diameter; Fe 67.6%, O 30.4%, Cr <1%, Cu <1%) was separated from the sample 3-4 cm (~2011 CE). Its outer surface exhibits a distinctive skeletal texture characteristic of rapid melting-quenching processes (Supplementary Figure S6).

Biological Indicators
Altogether, 156 diatom taxa representing 55 genera were identified in the sediment studied. The most abundant taxa include Staurosira costruens var. venter (Ehrenberg) P. B.  Figure  S7). Most of the taxa belong to benthos sensu stricto, tychoplankton, and periphyton, whereas the epipsamon, river benthos, and aerophyton ecological groups are rare (Figure 4). Although Tabellaria flocculosa is usually considered to be a planktonic diatom, periphytic forms can be also common. Therefore, we decided to keep this taxon as a separate category. The diatom record was divided into four significant diatom assemblage zones (DAZ 1-4) based on the relative abundances. DAZ 1 was further divided into two subzones (DAZ 1a, DAZ 1b), because of striking changes in the diatom valve concentration per 1 g and diversity (number of species per 400 valves) at the depth of 28 cm (Figure 4).

DAZ 1a (42-28 cm)
This (sub)zone is characterized by high abundances of benthic and tychoplanktonic diatom taxa, together with lower diversity and diatom valve concentration values. The zone is further characterized by fluctuation in the relative abundances of tychoplanktonic A. laevissima and Aulacoseira italica (Ehrenberg) Simonsen. Typical for those lowermost layers are high abundances of benthic colonial Fragilariaceae taxa including S. construens var. venter and P. brevistriata. Notable is the presence of acidophilous Eunotia species (Eunotia glacialis Meister, Eunotia tenella (Grunow) Hustedt or Eunotia minor (Kützing) Grunow) and a gradual increase in the D to C ratio. Chironomids, dominated by taxa associated with water macrophytes (e.g., Dicrotendipes nervosus-type and Psecrocladius sordidellus-type) and lake bottom sediment (Chironomus anthracinus-type and Tanytarsus lugens-type), reach their highest abundance and highest taxonomic diversity (number of taxa per samples), especially at the depth interval of 39-32 cm, peaking at 35.5 cm (131 HCs per 1 g of dry sediment). On the other hand, outside this depth interval, findings of chironomid HCs are rare (0-4 HCs per 1 g of dry sediment). The same trends were observed in abundances of bryozoan statoblasts and microturbelarian resting-egg capsules. Findings of ceratopogonid HCs and remains of planktonic cladocerans The onset of this (sub)zone is characterized by a sharp decrease of tychoplankton, whereas the abundances of benthos, periphyton and aerophyton increase markedly. Furthermore, remarkable changes in the diatom valve concentration and diversity occur in this section, with a significant peak of diversity and valve concentration at 27.5 cm, followed by a peak in chrysophyte cyst concentration (low D:C ratio) at 26.5 cm. The significant increase in periphyton is most apparent in the higher abundances of taxa indicative of greater habitat availability, such as Cocconeis placentula Ehrenberg, Gomphonema parvulum Kützing, and Epithemia adnata (Kützing) Brébisson. This zone is further characterized by the first occurrence of species of the genus Psammothidium, typically bound to sand grains (epipsamon), and river benthic taxa. In addition, several aerophilous diatoms, indicative of inwash from the lake catchment, such as Hantzschia amphioxys (Ehrenberg) Grunow, Pinnularia borealis Ehrenberg, and Luticola acidoclinata Lange-Bertalot are present. Chironomid abundance and taxonomic diversity as well as abundances of bryozoan statoblasts and microturbelarian resting-egg capsules reach the second highest peaks in DAZ 1b (up to 30 chironomid HCs per 1 g of dry sediment). The most dominant taxon, Chironomus anthracinus-type, is accompanied mainly by other eurytopic taxa living in bottom sediments (e.g., Cladotanytarsus mancus-type and Microtendipes pedellus-type). Tanytarsus lugens-type, an inhabitant of well-oxygenated cold lakes, is not present in the assemblages. Within the planktonic cladocerans, Bosmina is replaced by Daphnia and Ceriodaphnia.
In DAZ 1b, we also observed fine plant macroremains (length 1-5 mm) that were more abundant than in the other zones (Supplementary Figure S9).

DAZ 2 (25-9 cm)
The relative abundances of tychoplanktonic taxa markedly increase (up to 48%), with a gradual increase of A. laevissima from 15 to 33% at the onset of this zone, accompanied by a decrease in benthic taxa, diatom valve concentration, and diversity. Valve concentration and diversity are low at first, then rise gradually till the top of the zone. DAZ 2 is also characteristic by a rise and a subsequent steep decline in abundances of A. italica, and by extremely low abundances of chironomids HCs (0-3 HCs per 1 g of dry sediment). Chironomids associated with water macrophytes and remains of planktonic cladocerans disappear. In addition to bottom sediment associated chironomids, only bryozoan statoblasts of the genus Plumatella document the continuous presence of benthic invertebrates in DAZ 2.

DAZ 3 (9-1 cm)
A marked decrease in the abundance of diatom tychoplankton occur at the base of this section and benthic S. construens var. venter become the most dominant taxon. A gradual decrease of tychoplanktonic A. alpigena and A. laevissima is accompanied by a complete disappearance of A. italica. On the other hand, benthos, periphyton, aerophyton, and epipsamon increase, leading to relatively high taxonomic diversity of the assemblages. Sellaphora saugerressi (Desmazières) C.E.Wetzel and D.G.Mann, Pseudostaurosira trainorii E.A.Morales, P. brevistriata, and Achnanthidium minutissimum (Kützing) Czarnecki form minor but characteristic elements of this zone, and the relative abundance of Tabellaria flocculosa increase markedly. Water macrophyte chironomids appear again, but total abundances of chironomid HCs increase only slightly (1-10 HCs per 1 g of dry sediment). This increase is more pronounced in bryozoan statoblasts and microturbelarian resting-egg capsules. The dominant chironomid taxa are similar to DAZ 1b-Chironomus anthracinus-type, Dicrotendipes nervosus-type, and Microtendipes pedellus-type; however, several Tanytarsus lugens-type specimens were also detected. Findings of Daphnia and Ceriodaphnia ephippia are rare and situated in the younger half of the zone. Bosmina remains are absent.

DAZ 4 (1-0 cm)
A distinct shift in diatom assemblages occur in the two top layers of the core, showing a decrease in valve concentration and diversity with a simultaneous increase in chrysophyte cyst concentrations (high D:C ratio). The most prominent change in this zone is a gradual increase in abundances of T. flocculosa and a sudden disappearance of diatom tychoplankton, together with high abundances of river benthic (Planothidium frequentissimum (Lange-Bertalot) Lange-Bertalot) and aerophytic taxa (Microcostatus krasskei (Hustedt) J. R. Johansen and J. C. Sray, H. amphioxys). High abundances of Aulacoseira species are replaced by an increase in the abundance of species that inhabit shallow water environments. The benthic diatom Staurosirella pinnata (Ehrenberg) D. M. Williams and Round, first occurring in the lowermost layers of the core, suddenly becomes the dominant taxon. Apparent changes were also observed in the zoological indicators. The proportion of water macrophyte chironomids increases due to the first occurrence of Endochironomus impar-type in the assemblages, and a single finding of Eukiefferiella claripennistype, a rheophilic taxa, was observed. However, chironomid abundances remain low, similarly to in DAZ 3 (5-12 HCs per 1 g of dry sediment). The taxonomic composition of planktonic cladocerans also changes, with ephippia of the genus Ceriodaphnia replaced by Bosmina remains and abundant ephippia of g. Daphnia (63-92 ephippia per 1 g of dry sediment).

The Origin of Suzdalevo Lake in Relation to the Tunguska Event
The cores SUZ1 and SUZ3 were collected at the same location within the Suzdalevo Lake basin and show relatively similar geochemical, magnetic, and short-lived isotope activity records (Figure 3, Supplementary Figure S3). Due to this similarity, although partial differences are apparent, we used the 210 Pb dating of the SUZ3 core for interpretation of SUZ1 proxy data. In both cores, the observed peaks in 137 Cs (onset of the peaks at the depth of~16 cm; Figure 2, Supplementary Figure  S4) were interpreted as signals of the global fallout from atmospheric nuclear weapons with the 137 Cs maximum in 1963 CE (e.g., Appleby 2008; Bruel and Sabatier 2020), including the tests carried out by the Soviet Union at Novaya Zemlya from 1957 to 1962 (Khalturin et al., 2004), which is in good agreement with the results of the CRS age-depth model (1958)(1959)(1960)(1961)(1962)(1963) CE at the depth of 16 cm) for the SUZ3 core. This model is likely the most relevant one for dating of the younger half of the sedimentary record as it reflects well sediment compaction in middle parts of the cores. The signal of the Chernobyl accident in 1986, that was reported from many regions of the Northern Hemisphere, is not known from lake sediments investigated in the area of the TE (Rogozin et al., 2017;Darin et al., 2020). However, the maximum 137 Cs peaks in Suzdalevo Lake cores are not very prominent. In SUZ1, low 137 Cs activities were detectable also below the depth of 16 cm, likely due to partial mixing and/or remobilization and downward diffusion within the sediments (Matisoff 2017). Despite the fact that the "no sediment mixing" assumption for construction of the 210 Pb-based age-depth models is not fully met, all three models (CRS, CIC, CFCS) show that the base samples of the studied sediment cores are older than 1908 CE (Figure 2, Supplementary  Table S2). Since even the oldest layers contain valves of lake diatoms and remains of lake invertebrates, sedimentation in a lake environment is apparent throughout the whole record (Figures 4, 5). Moreover, the diatom and zoological indicator assemblages below the first layer with no 210 Pb ex activity (>36 cm, calculated ages >1908 CE) are not uniform (Figures 4, 5,  Supplementary Figures S7, S8) which contradicts potential very high accumulation rate in a lake formed by the TE. Therefore, Suzdalevo Lake was evidently formed before the TE. Due to the nature of the studied material (short cores sampled by a gravity corer) it is not possible to determine the time or the process of lake formation. It can only be speculated that it may have been associated with the selective melting of permafrost that is still present in the region (Gasperini et al., 2007;Ponomarev et al., 2019). Despite the oral testimony of the Evenki people, Suzdalevo Lake is not a TE impact crater filled with water, as has been proposed for the basin of Cheko Lake (Collins et al., 2008;Rogozin et al., 2017). Note that this conclusion is supported by the depth to diameter ratio being lower than 1/100 ( Figure 1C and Supplementary Figure S2), while typical simple impact craters have values between 1/5 and 1/7 (Melosh 1989;Włodarski et al., 2017).
On the other hand, the existence of the lake before the TE is not supported by identification of a specific "TE layer". No distinct clayey TE layer was observed, in contrast to those documented in the sediments of Cheko Lake and Zapovednoe Lake (thickness from 2 to 8 mm) by Kletetschka et al. (2019a) and Darin et al. (2020), respectively (Supplementary Figure S10). This is probably due to the fact that Zapovednoe and Cheko lakes are deep riverine lakes with large catchments, while Suzdalevo is a small shallow water body separated from the Chamba River ( Figure 1C, Supplementary Figures S1, S2). Sediments of riverine lakes and lakes with large tributaries are often very sensitive to erosion in a wider area compared to lakes with small catchments (e.g., Xu et al., 2010;Wilhelm et al., 2017). The most conspicuous, but not strong, clastic material input event in Suzdalevo Lake occurred in DAZ 4, based on the simultaneously increased values of mass normalized magnetic Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 777631 susceptibility, the Sr to Rb ratio, and Ti concentration in the SUZ1 core (Rapuc et al., 2020). However, DAZ 4 represents the youngest part of the sedimentary record and no similar erosion signal was found in the SUZ3 core. Another evidence of disturbance to the lake catchment was observed between the depths 25 and 28 cm (DAZ 1b), i.e., 14 cm above the base of the SUZ1 core (Figure 3). The DAZ 1b (SUZ1 core) contains abundant fine plant remains (Supplementary Figure S9), and its onset is characterized by increases in mass normalized magnetic susceptibility and the Sr to Rb ratio. The mass normalized magnetic susceptibility values reflect the concentration of potential magnetic carriers in the sediment (magnetite, pyrrhotite, maghemite, titanomagnetite, titanomaghemite, greigite); therefore, their enhancement in DAZ 1b reflects the presence of higher concentrations of at least one of these carriers. Since the area was exposed to fires due to the explosion and extreme radiation from the sky (Svetsov 2002;Jenniskens et al., 2019;Johnston and Stern 2019), the burning of organic material (e.g., grass, herbs, tree leaves and branches) has been associated with the conversion of ironbearing organic matter to magnetite (Kletetschka and Banerjee 1995). The increase in the Sr to Rb ratio also reflects the higher content of plant remains in the sediment samples rather than a decline in terrigenous clastic input (Xu et al., 2010). This is supported by the relatively stable Ti concentration values, indicating no distinct change in the proportion or grain size of clastic material (Cuven et al., 2010). The above mentioned geochemical and magnetic anomalies in DAZ 1b were, however, observed only in the SUZ1 core (Figure 3). The absence of similar shifts in the core SUZ3 implies variability of the disturbance signal strength within the lake basin. The DAZ 1b onset at the depth of 28 cm is in agreement with the CRSmodel-calculated depth for 1908 CE (27.75 cm), however, this "TE-depth" is placed ca. 5 and 4.3 cm deeper according to the CIC and CFCS models, respectively. It is important to notice that the time of the TE (i.e., 111 years before the sediment coring) is close to the limit of the 210 Pb dating method and the related models' applicability. Moreover, when the CRS model is applied, a "tooold" age error is always present for the deeper core layers due to underestimation of 210 Pb ex (Binford 1990).
We found no increase in concentration of melted magnetic microspherules in layers of potential TE age. Possible existence of a layer containing these objects would have been used for validating the age-depth models. The lack of melted magnetic microspherules in the SUZ1 sediments means that the TE explosion produced less than 354 magnetic microspherules per 1 m 2 , assuming the diameter of our gravity corer. This agrees with the known maximum concentration of these objects near the TE epicenter (90 magnetic microspherules per 1 m 2 , when objects with diameter from 10 to several hundred µm were counted) as reported by Badyukov et al. (2011). They studied material collected in the early 1960s by an expedition led by K. P. Florenskii, and distinguished magnetic iron-oxide microspherules and magnetic silicate microspherules with metallic droplets, both subdivided based on the presence/ absence of a Ni admixture. The single microspherule we found in the SUZ1 core (Supplementary Figure S6) corresponds to the category of iron-rich spherules without a Ni admixture that was, according to Badyukov et al. (2011), likely created by a hightemperature process from terrestrial rocks (i.e., melted and evaporated due to sudden heat pulse, and later nucleated and/ or quenched due to subsequent cooling). The same type of microspherules can be, however, of anthropogenic origin and transported for long distances if sufficiently small (<150 µm in diameter) (Zhang et al., 2011). The maximum calculated age of the extracted microspherule of 2003 CE (CFCS model, depth of 4 cm; Supplementary Table S2) indicates that it likely originated from the background micrometeorite flux, redeposition of an older TE sediment or soil, or anthropogenic activities. It seems that, due to their low concentration, these objects (diameter >3 µm) cannot be effectively used as markers of the TE in lake sediment cores from the Tunguska region. Also, the somewhat reduced TM grains that were found at depths of 30-26 cm (Supplementary Figure S5) cannot be interpreted as clear evidence of the TE explosion. The concentration of these grains is low (<3 grains per 1 g of dry sediment) and could have originated from the erosion of local volcanic rocks (Kamo et al., 2003). Such erosion in the lake catchment and subsequent inwash of the TM grains to the lake depocenter could have rather been associated with the TE-related forest disturbance. In future studies focused on lake sediments in the Tunguska region, the detection of increased concentrations of Ir and other Pt group elements may be more promising as follows from the "Northern" peat bog record presented by Hou et al. (2004).

Lake Ecosystem Changes at the Time of the Tunguska Event
In this study, we present the first diatom and freshwater zoological indicator records from the area exposed to the TE explosion in 1908 CE (Figures 4, 5). Throughout its recent history, Suzdalevo Lake lacked a significant planktonic diatom community and was dominated by tychoplanktonic Aulacoseira species together with small benthic alkaliphilic species of the Fragilariaceae family, commonly found in nutrient limited lakes with long periods of ice cover in alpine (Lotter et al., 1997;Karst-Riddoch et al., 2005), arctic (Douglas and Smol 1999;Laing and Smol 2000), and tundra regions (Westover et al., 2006). Similarly, the remains of planktonic cladocerans are rare in the sediments, demonstrating a greater suitability of the shallow lake habitats for benthic and littoral taxa (Szeroczyńska and Sarmaja-Korjonen 2007). This is supported by frequent findings of various benthic zoological indicators, such as chironomids, ceratopogonids, mayflies, bryozoans, or microturbelarians. All invertebrate taxa that were identified are known from the lakes in southern Siberia (42-65°N) or have a Holarctic distribution (e.g., Massard and Geimer 2008;Nazarova et al., 2008;Kotov 2016;Biskaborn et al., 2019).
According to the alternative age-depth models (Figure 2, Supplementary Table S2), the "TE-layer" and related signals of ecosystem changes should be located in DAZ 1a around the depth interval of 32-33 cm (CIC and CFCS models) or at the DAZ 1a/DAZ 1b transition, i.e., near the depth of 28 cm (CRS model). The first scenario is more likely in the case of the underestimation sediment accumulation rate in deeper layers mentioned by Binford (1990), however, the second scenario is also realistic when considering the error bars and thickness of the sediment layers used for the dating (4 cm). Moreover, only the second scenario corresponds to a distinct abrupt shift in lake biota (Figures 4, 5). The reality of the first scenario would imply an absence of a unique signal in the diatom and zoological indicator records. Such a conclusion would, for example, be consistent to the observations from sediments of Lake Strzeszyńskie (W Poland) that experienced the Morasko iron meteorite shower 4.9-5.3 kyr.cal. BP (Pleskot et al., 2018). The estimated entry mass of the Morasko meteoroid was between 600 and 1,100 tons and fragments of this cosmic body formed seven impact craters with diameters ranging from 20 to 90 m (Bronikowska et al., 2017).
If the age of the DAZ 1b onset is~1908 CE as suggested by the CRS model, the pre-TE status of Suzdalevo Lake is represented by the bottom zone (DAZ 1a) that is dominated by benthic colonial Fragilariaceae taxa (Figure 4, Supplementary Figure S7). These ecological generalists are common in lakes at high elevations and high latitudes, where ice-free intervals are short and nutrient supplies are low, although they also occur across a wide ecological gradient (Wolfe 2003). Furthermore, they inhabit the littoral zone, which is the first area to become ice-free (Westover et al., 2006). Thus, these taxa may quickly form blooms and outcompete larger diatoms during the short growing season (Smol 1988;Lotter et al., 1999;Grönlund and Kaupilla 2002). Fluctuations in abundances of all benthic zoological indicators ( Figure 5 and Supplementary Figure S8) imply the alternating presence of at least one adverse factor. Such severe suppression of benthic fauna in lakes is often caused by oxygen depletion (hypoxia) or even anoxia (e.g., Moravcová et al., 2021). Also, the presence of Ceriodaphnia ephippia and Plumatella and Cristatella mucedo statoblasts indicates limiting near-bottom oxygen concentrations caused by the intensive decomposition of organic matter accompanied by pronounced thermal stratification, limiting oxygen replenishment (Ursenbacher et al., 2020). However, the stratification and oxygen depletion were not permanent, as documented by the high abundance and high taxonomic diversity of chironomid HCs at depths from 39 to 32 cm, where both taxa that are tolerant (e.g., Chironomus anthracinus-type) and sensitive (Tanytarsus lugens-type) to low oxygen concentrations (Brooks et al., 2007;Tichá et al., 2019) were observed. Short periods of more intense water column mixing (affecting the sediment-water interface) could have occurred during the recent lake history as suggested by detectable low 137 Cs activities bellow the depth of 16 cm in both cores (Figure 2, Supplementary Figure S3), but such short-term changes cannot be reconstructed due to the time resolution of our records. The hypothesized changes in the mixing regime during DAZ 1a were likely triggered by an external factor, such as a climate change or a forest disturbance. A potential impact of floods is not supported by any erosion proxies (Figure 3) or presence of taxa indicative of running water environment (Figures 4, 5). Therefore, if a flood occurred, it did not have an erosive effect on the lake shore. Another possible explanation comes from the alternative age-depth models. If the "TE-depth" calculation using CIC and CFCS models is more realistic then in the case of the CRS model, an effect of the TE explosion cannot be ruled out as well.
As mentioned above, the assemblages in DAZ 1b could represent a direct response of the lake-catchment ecosystem to the TE (Figures 4, 5). In diatoms, the occurrence of aerophyton indicates the influence of the substantial transport of terrestrial material from the lake catchment, whereas the high abundances of periphyton and benthos suggest greater substrate variability. These changes therefore correspond to the observed increase in the amount of fine plant remains and the peak in mass normalized magnetic susceptibility (SUZ1 core), documenting a disturbance of the surrounding taiga. Diatom assemblages are known to also be structured by organic matter content and/or availability (Grimes et al., 1980), which significantly increased during this period. The diversified conditions, together with the input of organic matter, most probably led to increased diatom diversity and valve concentration (Adrian et al., 1999). Our data do not imply any conditions promoting the growth of acidophilous taxa, possibly connected with acidification or acid rains, as previously suggested by Kolesnikov et al. (1998). This also applies to the diatom assemblages in DAZ 1a. However, the peak in diatom valve concentration in DAZ 1b could have been supported by mild eutrophication after the increased deposition of nitrogen (N), estimated by Kolesnikov et al. (2003) to 200,000 tons across the 2,000 km 2 of tree fall area. This extra N supply was interpreted as a result of the high-temperature oxidation of nitrogen in the atmosphere with formation of nitrogen oxides during the impact (Kolesnikov et al., 1998(Kolesnikov et al., , 2003; however, it could also have arisen from another unconventional process. Boslough and Crawford (2008) modelled the TE explosion in the atmosphere and calculated an overpressure reaching 1 GPa at 6-8 km above the Earth's surface. Such a pressure wave allows the compression of N in air to reach near the liquid phase (Alkhaldi and Kroll 2019), and thus may have allowed the oxidation and creation of N species over the impacted area. The deposition value published by Kolesnikov et al. (2003) (1 ton per 1 ha) may be overestimated or apply only to sites located close the TE epicenter (distance <10 km); however, even a smaller deposition could have affected Suzdalevo Lake, as by shown by Bergström and Jansson (2006). Moreover, an input of P-rich dust after the TE explosionrelated thermal radiation and shockwave could be expected, transporting this key limiting nutrient to the studied lake and its catchment. No increase in P concentration was observed in the studied cores near the DAZ 1b onset or around the depth interval 32-32.84 cm (Figure 3), but trends in total P concentration sometimes differ from trends in biologically available P as shown by Norton et al. (2011). The taxonomic diversity and abundance of zoological indicators increase in DAZ 1b, with the dominance of Chironomus anthracinus-type and the absence of Tanytarsus lugens-type (Supplementary Figure S8) suggesting an improvement in near-bottom oxygen conditions but to a lesser extent than in the case of the period around the chironomid abundance and diversity peak at 35.5 cm (Johnson and Wiederholm 1989;Gąsiorowski and Sienkiewicz 2010). The potential destruction of the shoreline forest by the TE explosion could have caused increased wind speeds and shortened the periods of thermal stratification via wind mixing of the water column, resulting in improved dissolved oxygen concentrations (Klaus et al., 2021). On the other hand, the concurrent inwash of organic material could have accelerated respiration and oxygen consumption rates, which had the opposite effect and prevented the reappearance of the sensitive Tanytarsus lugens-type (Luoto 2013). In any case, the observed changes in DAZ 1b do not indicate any major disturbance or lake ecosystem collapse.
After the DAZ 1b, the ecosystem did not return to the pre-TE state. A new equilibrium, represented by DAZ 2, was established, and lasted until the 1990s (i.e., ca. 50-70 years) (Figure 4). We interpret this zone as a period of low water transparency and pronounced oxygen depletion, indicated by the substantial increase in the tychoplanktonic diatom taxa and extremely low abundances of all benthic zoological indicators, respectively. We assume that both processes were driven by the presence of additional dead plant biomass in the catchment after the TE. This led to a subsequent increase in dissolved organic carbon (DOC) production, and later to DOC leaching, DOC transport to the lake, and lake water brownification (Kopáček et al., 2018). The DOC leaching may have been amplified by reduced evapotranspiration in the disturbed forest, leading to increased soil moisture (Brothers et al., 2014;Kopáček et al., 2018). Brownification is a fundamental factor influencing lake ecosystem structure and function through more stable thermal stratification caused by increased absorption of solar radiation and limited light availability to benthic primary producers due to decreased water column transparency (Brothers et al., 2014;Solomon et al., 2015;Vasconcelos et al., 2016). These conditions likely caused the decline in benthic diatoms and especially in all benthic invertebrates, such as chironomids, ceratopogonids, bryozoans, and microturbelarians (Figures 4,  5). The chironomid abundances are so low in DAZ 2 (0-3 HCs per 1 g of dry sediment) that even episodic near-bottom anoxic conditions at the deepest part in the lake cannot be ruled out (Johnson and Wiederholm 1989;Houfková et al., 2017;Ursenbacher et al., 2020). These conditions would have supported Aulacoseira species (Supplementary Figure S7), which are adapted to low light availability (Kilham 1990;Bradbury et al., 1994), but most of them also indicate turbulent, unstable environments, likely caused by mixing of the water column. Such ecological characteristics are related to the production of heavy resting cells during the Aulacoseira life cycle that requires resuspension from the sediments into the lake water to establish (tycho)planktonic populations (Round et al., 1990). In addition, Aulacoseira species also create heavy silicified vegetative cells and therefore require an increased degree of turbulence to prevent their sedimentation (Saunders et al., 2009;Buczkó et al., 2013). In view of these facts, the stratification could not be long-term, but occurred frequently enough to disadvantage the benthic fauna. A similar record showing a long-term Aulacoseira presence in a shallow lake affected by anoxia was published by Houfková et al. (2017). We can speculate that a severe winter stratification due to long ice cover period was the key factor that suppressed the benthic invertebrates but enabled a summer bloom of Aulacoseira.
While the highest diatom valve concentration per 1 g of dry sediment was observed in DAZ 1b, the highest P concentrations and the highest values of the Si to Zr ratio, the other two potential proxies for algal primary production (e.g., Cuven et al., 2011), were observed in DAZ 2 ( Figure 3). We explain this discrepancy by the above-mentioned presence of large and heavy silicified Aulacoseira valves in DAZ 2, which could have increased the total diatom biomass and the concentration of biogenic silica (i.e., the Si to Zr ratio) in the sediment samples. Alternatively, the increased P content in DAZ 2 could be explained by an increase in P immobilization in the sediments caused by a higher aluminum hydroxide to iron hydroxide ratio in highly organic sediments during the brownwater phase, as shown by Kopáček et al. (2007). In the latter case, P concentration values would not reflect the availability of this limiting nutrient for aquatic organisms.
Since about the 1990s (DAZ 3), the lake biota in Suzdalevo Lake began to resemble those of the DAZ 1a (Figures 4, 5). The overall dominance of benthic diatoms (around 70% of the total diatom assemblage) may be attributed to processes linked to increased water transparency and with more stable conditions after recovery of the taiga forest, when DOC concentration in the lake water decreased and mature trees on the lake shore started to attenuate water mixing by wind (Kilham, 1990;Westover et al., 2006). The higher diatom diversity could be explained by the decreased dominance of Aulacoseira taxa, leading to more even diatom communities. These new conditions were also more favorable for benthic invertebrates, including taxa associated with aquatic macrophytes, that were absent in DAZ 2.
We also attempted to document changes in oxygen availability and related redox conditions in the sedimentary record using the Fe to Mn ratio ( Figure 3). This proxy is based on the fact that Mn is more easily mobilized and passes into solution more readily than Fe under reducing conditions (Davison, 1993). However, in Suzdalevo Lake, trends in the Fe to Mn ratio coincide with the Fe concentrations, indicating that the ratio was mainly affected by changes in Fe supply from the lake catchment (Mackereth 1966), and thus could not be used to reconstruct changes in redox conditions.

Recent Environmental Changes in the Suzdalevo Lake Sedimentary Record
The SUZ1 core provided an anomalous multi-proxy record for the uppermost two layers (DAZ 4), representing approximately the years 2016-2019 CE. In spite of the absence of these anomalies in the SUZ3 core that indicates variability of the signal strength within the lake basin or partial loss of the surface layer during the coring, the shifts observed in most of the proxies in the SUZ1 core are conspicuous and deserve a separate comment (Figure 3). The low Sr to Rb ratio and the sharp increase in Ti concentrations in DAZ 4 document an increased input of clastic material (Cuven et al., 2010;Xu et al., 2010), which was not observed in the other zones, including the depth interval where "TE-depth" is expected according to the age-depth models. Such a change in the proportion of inorganic particles could have been caused by Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 777631 increased erosion in the catchment or by flood activity of the nearby Chamba River (Supplementary Figure S1). Also, river benthic diatom taxa, such as P. frequentissimum (Figure 4, Supplementary Figure S7), reach their maximum abundance in DAZ 4, indicating a possible riverine influx and, together with higher abundances of aerophyton, may reflect the increased input of allochthonous material. In addition, the low diatom diversity and valve concentration and the dominance of S. pinnata and T. flocculosa suggest the persistence of relatively unproductive conditions throughout this period (Dam et al., 1994;Hofmann et al., 2011). The same zone is characterized by the presence of abundant remains of two planktonic cladocerans, Daphnia and Bosmina, while the abundances of the benthic zoological indicators remain rather low ( Figure 5). Among these benthic invertebrates, the only exception is the slight increase in chironomids associated with aquatic macrophytes and the solitary finding of the rheophilic chironomid Eukiefferiella claripennis-type (Brooks et al., 2007) Figure S8). Thus, the zoological record also suggests the influence of an external factor that changed conditions in the lake. DOC concentrations in boreal streams are often the lowest after high discharge events arising from snowmelt (e.g., Jonsson et al., 2007). If the factor responsible for the observed changes in lake biota was a snowmelt flood (or several floods during the period 2016-2019 CE) on the Chamba River, a drop in DOC accompanied with an increase in water transparency, along with other changes in water chemistry, could have occurred in Suzdalevo Lake. Such a shift would, for example, support the growth of aquatic plants and increase the abundances of the associated chironomid taxa ( Figure 5). Due to the presence of Bosmina, a cladoceran well adapted to fish predation, shorttime colonization of the lake by fish after a flood event cannot be ruled out either (Johnsen and Raddum 1987); however, we did not observe any fish during our field survey.

(Supplementary
We do not have data on floods on the Chamba River and we do not know of any observations of changes in the connectivity between Suzdalevo Lake and the river, as the site is very remote. Nevertheless, based on available data from the nearest meteorological station in Vanavara (Weather and Climate 2021), severe floods caused by rapid snowmelt or an ice jam near the study site could have occurred during the given period. At this station, the average annual temperature and precipitation between the years 1933 and 2020 were -5.5°C and 421 mm, respectively, while between the years 2016 and 2019 the average annual temperature increased to −4.1°C and average annual precipitation increased to 625 mm. This change in the local climate would have been favorable for stronger snowmelt floods and corresponds to trends observed in the wider area of Central Siberia (e.g., Tchebakova et al., 2011).

CONCLUSION
The TE in 1908 CE was an unusual extreme event that became part of the Evenki oral heritage. We performed this study to assess the potential for a TE-related origin of Suzdalevo Lake, a shallow water body located in the TE tree fall area, and to understand the effects of the catastrophic explosion on the lake-catchment ecosystem. Our data show that the lake basin did not originate from the impact, but lake sediments from the site provide a multiproxy record of a disturbance that occurred within a 5-cm-thick depth interval calculated using three alternative age-depth models for the position of the year 1908 CE. However, a clear link between this disturbance and the TE was not proven as the Suzdalevo Lake sediments do not contain an erosional clayey TElayer, that was observed in the Zapovednoe and Cheko riverine lakes, or any melted magnetic microspherules produced by the TE explosion. We interpret the observed abrupt changes in the proxy records as the consequences of the inwash of dead plant biomass to the lake, increased wind mixing of the water column possibly caused by deforestation of the shoreline, and mild eutrophication by elevated nutrient input. After the phase of higher productivity (i.e., increased concentrations and taxonomic diversity of diatom valves and zoological indicator remains), the lake did not return to its original state for about 50-70 years, likely due to the increased DOC supply from its catchment and associated truncation of the photic zone and near-bottom hypoxia. Such changes in a lake environment that may be related to the TE are reported here for the first time, however, more lakes in the TE tree fall area should be studied to reveal the real impact of the TE on ecosystems.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in the  and ensuring the smooth running of the entire expedition, Nicholas Hasson for help during the coring, Evžen Stuchlík for helpful comments on interpretation of our results, Petra Pokorná for help with R software, and David Hardekopf for language corrections. We also thank both reviewers for improving the manuscript via their comments.