Stages of soil development in the coastal zone of a disappearing lake—a case study from central Poland

The aim of this study was to estimate the impact of lake disappearance, groundwater table fluctuations, and some aspects of human pressure on stages of soil cover development in the catchment of Lake Rakutowskie. Research was conducted in soil chronosequence composed of eight soil profiles located in 785-m-long transect. Soils were sampled every 10 cm. Physical and chemical properties were determined with the use of disturbed and undisturbed samples. Radiocarbon dates were obtained for the bottom part of peat in each soil profile. The mineral composition of selected soil samples was determined using the X-ray diffraction (XRD) method. Additionally, the hydrological conditions were continuously monitored from July 27, 2013 to September 23, 2014. An increased rate of shallow lakes overgrowth and final disappearance is a commonly observed process in the European lowlands. This paper presents and discusses the influence of changes in the range of the Lake Rakutowskie wetlands complex (central Poland) and water level fluctuations on soil cover transformation in the immediate vicinity of the lake. Five major phases of soil cover development (from Subaquatic Histosols (Limnic) to Hypereutric Regosol (Drainic, Humic)) were hypothetically reconstructed based on detailed studies of soil physical and chemical properties, mineral composition, water level fluctuations, and radiocarbon dating. The degree of peat decomposition in the studied profiles increased with the distance from the contemporary lake surface. Along with the progressive dehydration and mineralization of organic soils, their physical and chemical properties deteriorated, which is visible, e.g., in soil structure, and mineral and chemical composition (especially TOC and TN). Obtained radiocarbon ages of selected samples indicate that the oldest peat layers developed about 5856 cal BCE (6970 ± 70 14C BP), while the youngest were deposited about 1074 cal CE (980 ± 40 14C BP)–1573 cal CE (300 ± 35 14C BP). Hydrological studies indicate a progressive process of lake disappearance. The process of soil cover transforming of the studied area takes place relatively quick, mainly recently due to the influence of human activity. Based on these results, it can be predicted that the lake will most likely disappear in a relatively short time, and the soil cover will be entirely in the fifth phase (P5) of evolution.


Introduction
The origin of most lakes occurring in northern and central Poland is associated with the ice sheet activity which occupied this area during the Pleistocene period. The results of paleolake reconstruction of Polish lakelands indicate that about 67.4% of postglacial lake surface has already disappeared (Kalinowska 1961), with 11.2% of this occurring in the last 70 years (Choiński et al. 2016). There are an increasing number of sites that either are or will be in the near future in the transition phase.
The process of lake disappearance worldwide is a stage of landscape evolution and its intensity depends on many environmental and anthropogenic factors. However, in the case of shallow lakes, the impact of human activity, mainly agriculture, can be the major factor determining their existence. The stages of lake basin evolution also depend on morphometric parameters and water characteristics. Shallow lakes are highly unstable ecosystems, sensitive to climate changes, in particular the climatic water balance (Scheffer 2004). Lake disappearance can be treated as a sum of two factors: surface area declines and shallowing followed by the lake succession process (Birks and Birks 1980;Cohen 2003;Rydin and Jeglum 2006). In addition to natural factors, the impact of anthropogenic processes, especially over the past decades, plays a significant role in lake area reduction (Scheffer 2004;Rydin and Jeglum 2006;Choiński et al. 2011;Markiewicz et al. 2017). The first phase of lake disappearance were climate changes which occurred in the Late Glacial and Holocene and caused lake condition modifications which had an impact on water level and thereby on sediment deposition. Global warming, which in recent decades has greatly accelerated, has a very large impact on the functioning of ecosystems. Climate projections for Poland foresee warming in all seasons and changes in the distribution and quantity of precipitation (Szwed 2015). Analyses of climatic water balance (the difference between precipitation P and evapotranspiration E) conducted by Bartczak et al. (2019) for the investigated area showed a negative trend for the years 1967-2015. The projected further decrease in the summer climatic water balance (Szwed et al. 2010) may adversely affect shallow lake ecosystems, such as Lake Rakutowskie, causing increased vegetation overgrowth and transformation into wetland. The succession of wetland vegetation in this area, in combination with an oxygen-poor water environment, leads to accumulation of poorly decomposed plant remains and peat formation, which is the last stage of lake functioning (Tobolski 2000;Succow and Joosten 2001;Rydin and Jeglum 2006).
Peatland characteristics are mainly determined by plant cover and the geomorphological and hydrological situation and are modified by climatic conditions (Charman 2002). Hydroclimatic instability most often leads to a lowering of the groundwater table (Słowińska et al. 2010;Lamentowicz et al. 2016;Samson et al. 2018), which influences changes in acrotelm processes (Morris et al. 2011;Jassey et al. 2013;Marcisz et al. 2014;Reczuga et al. 2018). Rivers, lakes, and peatlands are also susceptible to agricultural land use and catchment deforestation (Marcisz et al. 2015;Hernandez-Almeida et al. 2016;Kruczkowska 2016;Ott et al. 2017;Poraj-Górska et al. 2017). Most water bodies and adjacent areas are subject to drainage, which results in water level lowering, the uncovering of lacustrine sediments, and soil formation (Gonet et al. 2010;Mendyk et al. 2016) often with the participation of peat-forming vegetation. However, peat undergoes a mineralization process as a result of dehydration, which in turn leads to its degradation. All processes (natural and anthropogenic) leave a trace in Histosols, which can be treated as archives of natural peatland development, human activity, and the history of the catchment environment (Shotyk et al. 1998). The issue of organic soil evolution as a result of the impact of water lowering or anthropogenic processes has been widely studied (e.g., Renger et al. 2002;Pawluczuk and Gotkiewicz 2003;Becher et al. 2013;Długosz et al. 2018;Bogacz et al. 2019). However, reconstructions of soil development stages on partially or totally drained areas, including changes in soil typology, are relatively rarely represented in literature. Moreover, the majority of such studies are focused on physical and chemical property transformation as an effect of drainage, while there is no information about the mineral composition of Histosols and/or murshic soils that develop in coastal zones of disappearing lakes.
The aim of this study was to estimate the impact of lake disappearance, groundwater table fluctuations, and some aspects of human pressure on soil formation, morphology, and various properties in the catchment of Lake Rakutowskie. This lake is located in central Poland where precipitation deficit often occurs. The lake ecosystem is in a transitional phase as a result of vegetation overgrowth (Kramkowski et al. 2014) and is a good area for studies of pedogenesis in terrains adjacent to disappearing lakes.

Study area
Lake Rakutowskie is located in the wetlands complex in the immediate neighborhood of the Great Poland-Kuyavian Lake District (Fig. 1), where a rapid decrease in the lake area in Poland has been observed from the 1920s (Choiński et al. 2011). It is a specific reservoir with a large and shallow lake basin with an average depth of 1 m, where the surface area has decreased by two-thirds over the last century (documented range of the lake was 329.3 ha in 1930 and 101.8 ha in 2009) (Kramkowski e t al . 2 0 14 ; Bartczak et al. 2 01 9) . N o w , L a k e Rakutowskie is a eutrophic lake with a maximum depth of about 2.5 m and long, gentle slopes located within a valley in the southern part of the Płock basin (central Poland) within a young glacial landscape. The study area is a remnant of a past extensive paleolake developed during the Vistulian Glaciation (Glazik 1978). The lake is located within the largest complex of wetland ecosystems in the Gostynin Lake District and surrounded by flooded magnocaricetum grasslands and extensive biogenic plains with 2-m deep peat cover (Kramkowski et al. 2014). Quaternary sediments in the Płock basin are represented by a 35-40 m thick sandy-gravelly fluvioglacial series locally separated by marginal clayey mursh (Baraniecka 1993;Roman 2011;Kramkowski et al. 2014), which was remodeled aeolian during Late Glacial (Kruczkowska et al. 2020;Manikowska, 1982;Manikowska 1991).
The current shoreline of the lake is covered with rushes, predominantly Phragmites australis, Typha latifolia, and Schoenoplectus lacustris. Almost the whole lake bottom is covered with Chara sp. and in its bays Nuphar lutea and Nymphaea alba L. occur. In the rushes zone, the protected species Cladium mariscus has also been encountered. To the south Lake Rakutowskie is bordered with wet alder forest, while Fraxino-Alnetum and Ficario-Ulmetum minoris forests also occur. In the immediate vicinity of the lake, there are also small areas of crop fields that are extensively used for agricultural purposes due to the poor quality of the soils (Kramkowski et al. 2014).
The average sum of annual precipitation  within the Płock basin and surrounding areas measured at three meteorological stations (Baruchowo, Brześć Kujawski, Duninów) was 523 mm, which is lower than average for Poland (about 600 mm) (Bartczak 2007;Bartczak et al. 2019). The study area is located within the rain shadow (e.g., latitudinal zone of low precipitation (< 550 mm)) which extends from the eastern part of Poland through the center to the western borders of the country.
Water level fluctuations are controlled by natural factors, mainly precipitation and evapotranspiration, and human activity. Intensified deforestation, primarily adaptation of forest areas for agricultural purposes, is the major factor that affected the contemporary range of Lake Rakutowskie.
The percentage of arable lands in the vicinity of Lake Rakutowskie is small. The immediate vicinity of the lake is covered with hay meadows and pastures. The northern part of Lake Rakutowskie is covered with forest.

Hydroclimatic monitoring
The hydrological conditions were continuously monitored by the authors from July 27, 2013 to September 23, 2014. In Fig. 1 Location of the study area. a Poland in comparison with Europe, b Central and northern part of Poland, c Lake Rakutowskie in neighborhood of the Great Poland-Kuyavian Lake District order to record groundwater level fluctuations, we installed three PVC wells along the hydrological gradient within the peatland (from dry to wet sites) and one in the lake (Fig. 2). Each of these was equipped with a Leveloggers Junior Edge (Solinst Canada Ltd.). Water levels were measured every 3h, and they were compensated using atmospheric pressure and averaged to daily values. Meteorological records for the years 1967-2014 were derived from regional meteorological stations of the Institute of Meteorology and Water Management (Baruchowo station-precipitation, Płock-air temperature, and humidity), both of which were representative for the research area. We used standardized anomalies of monthly average air temperature, sum of precipitation, and climatic water balance from long-term (1967-2014) monthly averages. Thermal conditions were classified using the methodology proposed by Lorenc (1994) and precipitation and climatic water balance using the Guttman (1999) classification.

Soil sampling and analysis
Eight soil profiles were analyzed along a 785-m-long transect constituting the chronosequence of soils crossing varied plant communities comprising ecotone cover edge of willow bushes through a young alder forest to an old alder forest (Fig. 2).
The distance of the first profile from the contemporary lake surface was 1200 m, while distances between each soil profile were 65 m (1 to 2), 160 m (2 to 3), 120 m (3 to 4), 100 m (4 to 5), 110 m (5 to 6), 130 m (6 to 7), and 100 m (7 to 8). The profiles were located on a relatively flat terrain, and the height difference between the highest and lowest point was 0.8 m (altitude in m a.s.l.: profile 1-72.25, 2-72.20, 3-72.35, 4-72.50, 5-72.70, 6-72.85, 7-73.05, 8-73.00). The studied soil profiles were characterized by increasing degrees of soil organic material transformation resulting from lake drainage. The soils were described and classified according to the WRB soil system (FAO 2006;IUSS Working Group WRB 2015) and sampled every 10 cm. Two samples of undisturbed structure were collected from each layer/horizon using 100 cm 3 steel rings to determine physical characteristics of the soils, including bulk density and porosity. Moreover, one disturbed soil sample was taken from each horizon/layer to analyze chemical properties. The peat type and decomposition degree were determined in the field using a 3-point scale (R1-R3) according to Okruszko (1976a) and Myślińska (2001). Bulk density, total porosity, and actual moisture content were analyzed using undisturbed samples with the dryingweight method. Disturbed samples were dried at 40°C and milled into powder. The following properties were determined from these samples: loss on ignition at 550°C, pH potentiometrically in a suspension with water in a soil:water ratio of 1:2.5 (pH-meter Elmetron CPC-401), the content of carbonates according to Scheibler's volumetric method, the content of total organic carbon (TOC) with the Alten method (Dziadowiec and Gonet 1999), and total nitrogen (TN) with the Kjeldahl method (Van Reeuwijk 2002). Samples were digested in a mixture of 40% HF and 60% HClO 4 , and the total content of phosphorus in solution was analyzed using the molybdenum-blue method, whereas contents of potassium, Fig. 2 The study area. a Soil profiles location. b Climatic conditions and Lake Rakutowskie water level changes calcium, magnesium, iron, chromium, copper, manganese, nickel, lead, and zinc were determined using flame atomic absorption spectrometry (Perkin Elmer AAS 2100).
The mineral composition of soil samples was determined using the X-ray diffraction (XRD) method. Samples were dried at 40°C and ground into powder. The powders were analyzed by means of a Bruker AXS D5005 diffractometer equipped with a KRISTALLOFLEX® 760 X-ray generator, a vertical goniometer, 1-mm divergence slit, 2-mm anti scatter slit, 0.6-mm detector slit, and a graphite diffracted-beam monochromator. CoKα radiation was used with the applied voltage of 40 kV and 30 mA current. Random amounts of the ground materials were scanned at a counting time of 2 s per 0.01°step from 3 to 70°2θ.

Radiocarbon dating
Radiocarbon dates were obtained for the bottom part of peat in each soil profile. Dating was performed in the Laboratory of Absolute Dating, Cianowice. Samples for radiocarbon dates were taken from transition between lake sediment and peat deposits. Conventional dates were calibrated using OxCal 4.2.3 software (Bronk Ramsey and Lee 2013) and the IntCal13 calibration curve (Reimer et al. 2013). In this paper, dates are also presented in calibrated years CE/BCE.

Hydroclimatic characteristics of the Rakutowskie wetlands complex in 2013-2014
Most of the months in the investigated period July 2013-September 2014 were characterized by normal air temperature fluctuations relative to long-term averages   (Fig.  2). However, March and July 2014 were very warm, and November-December 2013 and April 2014 were warm, which caused the hydrological year 2014 to be classified as very warm. Based on classification with the standardized precipitation index (SPI) by Guttman (1999), the hydrological year 2013 was neutral and 2014 was moderately dry. The same situation was noted for climatic water balance.
Groundwater level fluctuations at all locations were controlled by seasonal changes in hydrometeorological conditions (precipitation and evapotranspiration), with the highest levels noted during the winter season and lowest in the summer. The wettest place with the most stable water conditions was in location w1. Median groundwater level in this location in the period September 27, 2013-September 23, 2014 (when groundwater level was recorded at all sites) was 3 cm, and from the end of December 2013 until mid-June 2014 water stagnated about 10 cm above the ground. The groundwater level decreased with increasing distance from the lake. At location w2, median groundwater level was 62 cm, and at w3, it was 90 cm below the surface. During the summer months in 2014, groundwater level dropped to 144 cm at location w2 and 178 cm at location w3. At location w1, it was 52 cm below the surface.

Radiocarbon age of the soils
From profile 1, the transition zone between the gyttja and peat was dated to approximately 1,573 cal. CE, then through profiles 2, 3, and 4 the transitional age was determined as follows: 1,074 cal. CE, 751 cal. BCE, and 850 cal. BCE. The last documented phase of the Lake Rakutowskie overgrowth was dated to about 5,856 cal. BCE (profile 5). Radiocarbon ages of organic material from profiles 6-8 were rejuvenated and thus rejected due to the overly small thickness of topsoil which did not exceed 20 cm, its total transformation into mursh, and the numerous roots of live vegetation which caused the sample to be contaminated (Table 1).

Morphology and typology of the soils
The immediate vicinity of Lake Rakutowskie, which was part of the former lake bottom, was mainly covered with soils developed from organic materials (peat and/or mursh) in the topsoil and lacustrine sediments (gyttja and/or lacustrine chalk) occurred at different depths in the subsoil (Fig. 3). The characteristics of the investigated sediments changed from being biogenic deposits consisting of detritus calcareous gyttja, detritus gyttja, and lacustrine chalk on the bottom to being alder swamp forest peat and reed peat in the upper part. Among the eight soil profiles, five soil reference groups were distinguished in the studied transect: Hypereutric Hemic Histosol (profiles 2 and 3), Hypereutric Hemic Histosol (Limnic) (profile 1), Hypereutric Sapric Histosol (Drainic, Limnic) (profiles 4 and 5), and Hypereutric Regosol (Drainic, Humic) (profiles 6-8) (IUSS Working Group WRB 2015). The soil morphology of each profile is presented in Fig. 3. Additionally, lake sediments which were permanently submerged by water not deeper than 200 cm were hypothetically classified as Subaquatic Histosols (Limnic) and Subaquatic Histosols.

Physical properties
Groundwater level lowered and soil moisture decreased with the distance from the lake, which reflects the percentage of water-filled pore space (Table 2). Peats in the studied soils can be classified as moderately and significantly decomposed ( Table 2).
The percentage of organic soil substrates in the morphology of profiles 1-3 varied from 40 to 100%, while in profiles 4-8, this ranged from 20 to a maximum of 40%. Additionally, the percentage of organic horizons with a low degree of decomposition (R1 and R2) decreased with the distance from the lake, and this also coincided with a decrease in the thickness of organic horizons. Different degree of organic material decomposition in profiles 1 and 5 depends on distance from the lake, despite of morphological similarity of these two soils. The soils with the most decomposed peat mass (profiles 4 to 8), with an R3 degree of decomposition in the whole profile or with features of peat mineralization, were located at the furthest distance from the lake (Fig. 3). In these profiles, peat layers were totally or almost totally transformed into mursh with a mineral fraction content of about 30%. For comparison, in profiles 1 and 2, the mineral fraction content in topsoil did not exceed 22%. In the first three profiles (1-3), peat decomposition was also observed; however, significant transformations (R3 degree of peat decomposition) occurred, maximally to a depth of 30 cm.
Transformation of organic material due to hydrological changes is reflected in soil bulk density and porosity (Table 2). Organic material from profiles 1 to 5 was still classified as peat, however, with a different degree of decomposition. In profiles 1 to 3, a decrease in bulk density with depth was observed, in contrast to the other profiles. An opposite trend occurred in the case of total porosity. A boundary between moderately and significantly transformed soils can be set between profiles 3 and 4, where there was a noticeable change in the soil physical properties ( Table 2). The bulk density of superficial H horizons ranged from 0.21 to 0.24 g· cm −3 in profiles with hemic material (He horizons), 0.35 g· cm −3 in sapric material (Ha horizons) to 0.37-0.47 g·cm −3 in Hypereutric Regosol (Drainic, Humic ) profiles. The average percentage of the mineral fraction in surface H horizons is closely related to the degree of peat mineralization and ranges from 22.9% in profiles 1-3 to 30.7% in the most transformed profiles (4-8).

Chemical properties
The peat materials in the investigated soils had slightly acidic to alkaline reactions (pH from 6.2 to 7.6). Lacustrine sediments were characterized by similar pH values (from 5.9 to 8.1). However, in profiles 1 to 3, differences in soil pH between the topsoil and the subsoil were relatively aligned. Clear  pH differences between peat and lake sediments occur in profiles 4-8. Changes in the TOC and TN content were indicated in topsoil with different degrees of organic material mineralization. In profiles 1-3, TOC contents ranged from 235.91 to 539.49 g kg −1 and TN from 17.01 to 38.44 g kg −1 , with a relatively increasing trend from the top to the bottom part of the profiles (Fig. 4). In profiles 4-8, this dependency was reversed. Gyttja of the least transformed profiles (1 and 2) showed a very high content of TOC (426.60 g kg −1 on average) and TN (30.11 g kg −1 on average), whereas moderately and strongly dehydrated profiles 4-8, on average, had up to 72.70 g kg −1 (TOC) and 4.66 g kg −1 (TN).
Increased Fe content was observed throughout the profiles of the most dehydrated soils (5-8), especially in the layer 10-20 cm below the surface (Fig. 5), while increased Mg content occurred mainly below 20 cm. At a depth below 20 cm, there was a rapid decrease in Fe content. Lacustrine sediments localized in the bottom of profiles 4-8 were much poorer in specified elements than in profiles 1 and 2. In profiles 1-3, Fe, Ca, and Mg contents were higher in topsoil than in the subsoil, in contrast to the more transformed soils (profiles 4-8) in which that trend was reversed. Al content was relatively similar in all profiles, with the highest values in surface horizons (4.40-7.49 g kg −1 ); however, in profiles 1 and 2, values below 20 cm decreased by 3 g kg −1 on average, while in other profiles, these differences reached 6 g·kg −1 . Similar regularity could be noted for K content, which achieved the highest values to a depth of about 20 cm (1.14-2.79 g kg −1 ). Below this depth, K contents decreased to 0.05-1.01 g kg −1 . In profiles 1 and 2, the difference between the maximum and minimum value was 1.15 g kg −1 on average, while in dehydrated profiles (4-8), this increased to 2.48 g kg −1 . In these profiles, a decrease in the K content with depth was clearly marked. The highest P contents were observed to the depth of 20 cm in drained profiles (2.18-2.80 g kg −1 ). In all profiles, below 20 cm, there is a sudden decrease in P content (to 0.06-0.84 g kg −1 ). In profiles 1-3, there were also smaller differences in P contents between surface and deeper horizons. TOC/TN ratio ranged between 12.1 and 20.0. The greatest differences were visible from the depth of 20 cm to the bottom part. In profiles 2 and 3, the range was wider and reached maximum values  (over 18) in the middle and bottom parts. A narrower range of TOC/TN values (13-15) was characteristic for profile 1 located in the immediate vicinity of the lake's current surface. The TOC/P and TN/P profile distribution patterns were fairly similar in each profile, while an increase in the value of these indicators was observed below the depth of 30 cm. Surface horizons of drained soils were characterized by values on average 60 (TOC/P) and 5 (TN/P) units smaller than in other profiles (Fig. 4). On the basis of all the studied profiles, profile 2 was untypical: the average value of the TOC/P index was 559 (with values from 279 to 522 in other profiles), and TN/ P-35 (with values of 19-30 for other profiles).

Mineral composition of the soils
The topsoil of the studied soils were rich in organic matter which resulted in high background and noises in the XRD patterns of these samples (Fig. 6). The predominance of organic matter was also confirmed by results of TOC analysis (Fig. 4). Moreover, organic materials contained admixtures of quartz and calcite (Fig. 6).
Calcite was the major mineral in the lacustrine sediments occurring in subsoil of the studied soils (Fig. 6). Moreover, minor amounts of aragonite were determined in several samples in the upper part of lacustrine chalk layers in contact with organic material layers (e.g., profiles 6 and 8). Aragonite occurred in those horizons where shells of molluscs were found. There was a weak broad peak 0.604 nmin the XRD patterns of samples from the subsoil. This peak most likely belonged to bassanite. Trace amounts of goethite and wavellite were found in the subsoil of profile 8 (Fig. 6).

Phases of soil cover development as an effect of lake water level changes
Radiocarbon ages of lake sediments collected to the depth of 785 cm in Lake Rakutowskie presented by Harasimiuk et al. (2010) encompassed dates between 5590 ± 60 cal. BP (at the depth of about 200 cm) and 11,730 ± 100 cal. BP in the bottom part of the core, which indicates the existence of a reservoir at least from the Allerød. Dates obtained from the samples collected from the bottom part of a 735 cm deep core from Lake Żłoby were older (to 14,530 BP; mean 15,792 cal. BCE) and were determined in the fine sand layer which indicates an accumulation process in the Late Glacial period. These data correspond to those for other lakes from central Poland (Manikowska 1992;Więckowski 1993;Tobolski 1998). The current lake area may vary between 200 and 400 ha depending on the period of the year (Harasimiuk et al. 2010). The oldest Histosols, located at a distance of 1200 m from the lake within the study area, were dated to 5856 cal. BCE which indicates that the water level of Lake Rakutowskie was already in decline in the Atlantic period and that it gradually lowered in the subsequent decades. Due to the occurrence of lake sediments at further distances from the studied transect, it is highly probable that the water level decrease took place much earlier. In the past, Lake Rakutowskie and neighboring Lake Żłoby were probably part of a broad basin with deeper water parts within Lake Żłoby.
The results obtained based on the sequence of the investigated profiles (Figs. 2 and 3) clearly show the transformations of organic soils (Histosols) along with disappearance of a lake. The phases of soil-forming processes in the immediate vicinity of Lake Rakutowskie were closely related to water level fluctuations. Shallow lake sediments such as those of Lake Rakutowskie may be considered as soils provided they were permanently submerged by water not deeper than 200 cm (IUSS Working Group WRB 2015). In the studied area, this condition could probably be met, especially on the outskirts of the lake. For this reason, we assumed that all studied soils before dehydration could primarily be hypothetically classified as Subaquatic Histosols (Limnic) and Subaquatic Histosols (IUSS Working Group WRB 2015). This stage (underwater soil phase-P1) can be treated as the first phase of soil cover development in this area. According to, e.g., Demas et al. (1996), Demas and Rabenhorst (1999), and Erich et al. (2010), aquatic substrates (sediments) in shallow lakes can be transformed into soils as a result of pedogenic processes. The second phase (P2) after the water table lowering and terrestrialization was characterized by full soil saturation with water without traces of peat degradation. This stage was observed in profiles 1 and 2. The next stage of soil development (P3) was represented by profile 3, where as a result of significant water table decrease, the first signs of peat degradation were observed in topsoil. Phase four (P4) was characterized by significant dehydration and clearly marked degradation of peat, which was observed in profiles 4 and 5. Permanent drainage of the area caused a phase of deep peat degradation as a result of the advanced mursh formation process (P5) which is represented by profiles 6-8. The presence of limnic material in Regosols profiles may contribute to the inclusion of Limnic as additional supplementary qualifier in relation also to Regosols which is missing in the current version of the WRB classification.

Soil physical and chemical properties modified as a result of changes in hydrological conditions
Gradual lake drainage led to a decrease in its range and consequently to peatland formation and then to intensification of peat decomposition. Peatland functioning is closely related mainly to hydrological conditions (Strack 2008), which in cases of its radical change can have a strong impact on organic soils. The first noticeable effects of drainage are changes in soil morphology, structure, and physical properties. Mursh is the final product of intensive and long-term dehydration. This process, especially under conditions of constant peat soil drainage, leads to structural changes in soil profiles, gradual delamination of the peat layer and, as a result, its increasing permeability. Differences in soil morphology modifications depend on the distance from the current lake surface. Groundwater level lowering in soils enhances the progressive vegetation succession. Peat dehydration causes its progressive degradation and mursh formation by the intensive mineralization of organic material, which affects the content of elements  Figs. 4 and 5). Such a pattern has also been observed by other researchers (Berglund 1996;Chmieleski et al. 2004;Gonet et al. 2010;Mendyk et al. 2014;Markiewicz et al. 2015;Glina et al. 2016a, b). Degraded peatlands with desiccated soils have a damaged ecosystem which cannot function properly. In profiles located at the shortest distance from the present lake surface, water level ranged from 20 (1) to 40 cm (3) below soil surface, while in other profiles this was deeper than 80 cm. Deterioration of soil physical properties as a result of desiccation is reflected in a increase of the bulk density and decrease of soil porosity and also a higher percentage of mineral fraction in primarily organic sediments, similar to other studies (e.g., Hobbs 1986;Myślińska 2003). In profiles 4-8, murshic soils have been formed with degraded top parts of peat. The best preserved morphology of deep Histosols was observed in profile 3, where peat material was preserved with lacustrine sediments below 100 cm. The peat dehydration process causes a significant deterioration of the physical properties of these soils (e.g., Renger et al. 2002;Kechavarzi et al. 2010;Kiryluk 2017). In the study area, this relationship increases with the distance from the contemporary lake range. The progressive decrease in groundwater level caused a significant acceleration of peat decomposition, which is particularly evident in profiles 4-8. In these profiles, a significant degree of decomposition covers the entire layer of peat material. In other profiles (1-3), in which groundwater occurred at a depth of 20-40 cm, the intensive peat decomposition process was observed only within topsoil.
Peatlands constitute a major terrestrial carbon stock. The problem of organic soil dehydration is closely related to the loss of organic carbon, soil usability, and development resulting in an imbalance in peatland ecosystems (Gonet and Markiewicz 2007). Changes to organic matter in soils occur with varied intensity, expressed as a TOC/TN ratio. The studied profiles were characterized by relatively similar TOC/TN ratios; however, higher values (17.0-20.0) in the bottom part of profile 3 and the middle part of profile 2 could be an indicator of organic matter mineralization slowdown. All surface horizons were subjected to enhanced mineralization processes which is evidenced by lower TOC/TN values (12.1-13.5). Under the influence of the mursh forming process, the content of TN increases by 30%, while the content of TOC decreases (Piaścik and Gotkiewicz 2004). Soil pH in the studied transect was closely related to characteristics of peat and lacustrine sediments. In view of these reactions, the studied profiles can be divided into two groups: (i) soils with relatively homogenous pH across the whole profile and with a marked influence of organic soils on lacustrine sediment reactions (profiles 1-3), and (ii) dehydrated soils with clear reaction changes at the border of organic and lacustrine sediments and marked lacustrine sediment impacts on pH of overdried soils (profiles 4-8). The presence of Ca in organic horizons and thus high pH values, could be correlated with the capillary ascension of water and pedogenic calcite crystallization in organic materials. The large quantities of calcium ions present in solution in the peat layer could also be a result of organic matter binding (Lucas 1982;Okruszko 1976b).
The two or three times higher contents of iron in murshic horizons than in peat may be caused by organic matter mineralization enhancement or oxidation processes intensification (Mocek et al. 2007). Differences, especially in Fe content, between each profile resulted from the rate of the drying process, which first occurs in surface horizons and currently covers the greater part of the Rakutowskie catchment. In overdried Histosols, iron has an increased ability to form complex compounds with an organic substance, which is confirmed by other researchers (Henrot and Wieder 1990;Norrstrom 1995). Iron in the mursh material comes from the peat mineralization process and from groundwater, and it precipitates in the form of oxyhydroxides. A high Fe content in murshic horizons is an indication of increased oxidation and reduction processes and also of intensive mineralization and biological accumulation (Puustrjarvi 1952;Sapek and Sapek 1987). The zone of oxidation in overdried soils can be located at different depths depending on the type of water regime.
The mineral composition of the studied soils is different in the subsoil (a layer of lacustrine deposits) and the topsoil (a layer of organic material). The mineral composition of the subsoil of the studied soils is typical for similar formations in temperate climatic zones. Calcite is the major mineral phase occurring in the lacustrine chalk and certain gyttja types (Wyrwicki 2001;Raukas and Tavast 2002;Maćkowska 2008) present in the subsoil of the studied profiles. Bassanite has previously been found in lacustrine chalk (Maćkowska 2011). This mineral is likely to be an effect of the reaction of weathered pyrite with Ca carbonates occurring in lacustrine chalk (Maćkowska 2011). The lacustrine chalk layer contained numerous shells and therefore contained aragonite. That mineral is one of the major components of mollusc shells (Medaković et al. 2003). Particularly high contents of aragonite occurred typically in the upper part of lacustrine chalk layer in the studied soil profiles (e.g., profiles 6 and 8) (Fig.  6). As mentioned above, the calcite genesis in organic horizons is probably related to the capillary ascension of water and calcite crystallization in organic materials. Quartz present in organic layers of the studied soils is most likely an allochthonous phase.

Impact of human and land use changes on soil transformations
Changes of the lake range caused mainly by anthropogenic factors led to rapid overgrowth of the former lake surface by reed beds, willow clumps, etc., and as a result, soil formation and their progressive overdrying. This process was also observed by other research (Ralska-Jasiewiczowa and van Geel 1992;Smith et al. 2005;Kinder et al. 2009;Dietze et al. 2016;Skowron and Jaworski 2017;Słowiński et al. 2019). Based on the radiocarbon dating results, we can observe the gradual process of the Lake Rakutowskie disappearance (Table 1). Reduction of the lake surface within 75 years from 329.3 to 101.8 ha has accelerated significantly in the last century, and this has an undeniable connection with land use changes in the lake catchment, drainage, and fertilization. The results indicate a significant acceleration of the disappearance of the lake after the implementation of drainage treatments. This process before the melioration was introduced was slow and gradual. Soil morphology furthest away from the contemporary lake surface indicates that peat was degraded completely. The drainage of the area by humans has undoubtedly accelerated this process.
In the early phase of the lakes functioning, they were mainly dependent on precipitation, river water inflow, evaporation, and the variable relationship between natural erosion and accumulation processes leading to gradual lake shallowing (Tobolski 2000;Kramkowski et al. 2014;Ott et al. 2017). Analysis of archival cartographic documentation of the immediate vicinity of Lake Rakutowskie indicates the increased impact of human activity as early as the beginning of the twentieth century (Kramkowski et al. 2014;Bartczak et al. 2019). Initially, this work consisted of the regulation of watercourses, and then by 1930, a network of drainage ditches was established. Further regulatory work on the watercourses led to a significant reduction of the lake area from 329.3 ha (in 1930) to 226.8 ha by 1945. Since the 1960s, the lake area was still rapidly changing. The next phase was the overgrowth of watercourses and more intensive use of adjacent meadows for agricultural purposes. Together with the lowering of the lake's water level, new bottom fragments emerged, which have gradually overgrown. This process led to a lake area decrease to 101.8 ha by 2009 (Kramkowski et al. 2014). On the basis of observations of the lake's range over the last 100 years, a significant acceleration of the lake's rate of decline is clearly visible in relation to the period before the intensification of anthropic pressure. Analogous results have also been obtained during research on other lakes (e.g., Mendyk et al. 2014Mendyk et al. , 2016Markiewicz et al. 2015). The obtained C 14 dates indicate a significant acceleration of the lake's disappearance process over the last 300 years. However, in relation to the results presented in Fig. 2, the fastest decay rate is in the last 70 years. Similar processes are observed in other European countries (e.g., Coops and Hosper 2002;Battarbee and Bennion 2011;Ventelä et al. 2015).
Our studies have shown that soil sequence analysis in the shoreline of disappearing lakes may be an indicator of environmental changes and the basis for reconstruction of former lake ranges. Soil surveys are a valuable complement to other forms of research (hydrological, ecological, etc.).

Conclusions
As a result of more than 12, 000 years of evolution, Lake Rakutowskie has been subjected to lake water level changes. In the past decades, the rate of Lake Rakutowskie disappearance has accelerated due to human activity causing irreversible ecosystem changes. The whole area of the former Lake Rakutowskie bottom is covered with soils originating from lacustrine sediments and gyttja covered with organic materials (peat and mursh) characterized by varied thicknesses and degrees of decomposition. However, a significant acceleration of this process is undoubtedly related to human activity in the immediate catchment (land use changes, melioration, and eutrophication associated with agriculture). Progressive soil dehydration has taken place since 1945, when the lake area was reduced by 30%, which led to vegetation succession on exposed parts of the former lake bottom. According to soil morphology, five phases of soil development were specified, from the stage of subaqueous soils (Subaquatic Histosols (Limnic) and Subaquatic Histosols) still observed in the shallow parts of the lake, to the 5th phase represented by Hypereutric Regosol (Drainic, Humic ) profiles, which began at latest during the Atlantic period (5856 cal BCE (6970 ± 70 14 C BP)). The study stand is characterized by low variability of soil morphology, however differing in peat thickness and its decomposition degree (composed of peat, mursh, and lacustrine sediments) which indicates the relative homogeneity of the conditions of their development affected by the former lake range. Different degrees of soil dehydration influenced soil morphology (changes of peat decomposition degree), and soil physical (decrease in total porosity and increase in bulk density, higher percentage of mineral fraction content, etc.) and chemical properties. In the analysis of soil property changes, the distance from an existing lake plays an important role. The most intensive transformations were characteristic for TOC, TN, Fe, Al, Mg, Ca, and P contents in surface horizons. The calcite content results from the past/present impact of lacustrine sediments. Despite the organic matter, quartz and calcite predominate in the mineral composition of topsoil. In the subsoil, aragonite, basanite, and also trace amounts of wavellite and goethite occurred. It can be predicted that the area of dehydrated ecosystems in the vicinity of Lake Rakutowskie will gradually increase due to the strong impact of drainage and changes in lake trophic status as a result of the supply of artificial fertilizers from the surrounding arable fields causing a further decline in the lake's water level. The results can be extrapolated to similar postglacial areas in temperate climatic zones. Additionally, these studies showed that term Limnic should be added as supplementary qualifier in relation to Regosols in the WRB classification.