From oxbow to mire: Chironomidae and Cladocera as habitat palaeoindicators

Chironomidae and Cladocera are useful palaeoindicators to describe long-term dynamics in biodiversity of freshwater aquatic communities. However, palaeoecological studies relying on these aquatic invertebrates have generally focused on post-glacial lakes, with relatively few applications to floodplain water bodies such as palaeo-oxbow lakes. In this study, zone-level (stratigraphic) trends in Chironomidae and Cladocera assemblage composition were analysed from a small palaeo-oxbow lake in a Central European river valley of Poland. Summer air temperature and hydrological conditions were found to be the main drivers of Cladocera and Chironomidae populations, causing a cascade effect in assemblage composition. At the local habitat scale, there were autecological differences in both taxonomic groups reflecting different responses to abiotic and biotic factors. Specifically, Cladocera responded distinctly to plant vegetation and fish predation during the Late Glacial lake phase, whereas Chironomidae were more sensitive to paludification processes and disappeared in the Northgrippian dry phase of the mire. Both groups indicated fluvial activity of the River Grabia, although the response was different. Also, species richness was higher in the more taxonomically diverse Chironomidae compared to Cladocera. As both groups responded to the same factors in a different way, it is recommended that they should be used in parallel in palaeoecological research studies.

Whilst palaeoecological studies relying on Chironomidae and Cladocera have generally focused on post-glacial lakes, relatively few investigations based on these aquatic invertebrates have dealt with floodplain water bodies such as palaeo-oxbow lakes. This is a notable knowledge gap, as paludification is a common process that in Europe has intensified since the Greenlandian (lower age of the Holocene), with river valleys in Central Europe undergoing transformation from the Late Glacial (end of Pleistocene) to the Meghalayan age (uppermost Holocene) and with valley stagnant waters turning from lakes in the Younger Dryas to mires in the Northgrippian (middle age of the Holocene) (Vandenberghe, 2003;Starkel et al., 2013).
The present study investigates the responses of Chironomidae and Cladocera in assemblage composition and abundance to regional and habitat-related environmental factors in a palaeo-habitat that turned from oxbow to mire. It is anticipated that the outcomes of the present study will contribute to the use of Chironomidae and Cladocera as proxies in palaeoenvironmental and palaeoclimatic reconstructions based on data from sediments of lakes, flood plains and valley mires. A novel analytical approach is also implemented that can be applied more widely in palaeolimnological studies.

Study area
The Pawłowa site (51°30 0 19 00 N; 19°19 0 59 00 E) is located in the Łódź region of Central Poland (Fig. 1a) at an elevation of 185 m a.s.l. in the middle part of the course of the River Grabia (Fig. 1b) and in the marginal part of a valley floor that occupies a slight depression within the flood plain (Fig. 1c). The catchment of the Pawłowa mire is about 7.7 ha and is presently covered for the major part by an alder swamp and meadows. An area surrounding the Grabia valley consists of glacial and glaciofluvial deposits that accumulated during the Saalian glaciations, especially of the Warthanian Cold Stage (the last one to cover the Łódź region: Klatkowa, 1984), and fluvial sediments (Late Glacial and Holocene alluvium) that fill the main part of the valley. General descriptions of the geological and geomorphological characteristics of the study site, including lithology and chronology of sediments, are provided by Pawłowski et al. (2016a) (Fig. 2).
Currently, the Łódź region has a temperate climate influenced by westerly airflow from the Atlantic and by eastern Eurasian high-pressure cells (Woś, 1999). Mean annual temperature in the study area is 7.7°C, mean temperature of the warmest month (July) 18°C and of the coldest month (January) -3.3°C. Average annual precipitation is 590 mm, ranging from 438 to 937 mm (Kłysik, 2001).

Sample collection and processing
A 450-cm sediment core from Pawłowa consisting of gyttja and various types of peat was taken in the southern part of the mire (Fig. 1c, d), where the thickest organic deposits were found, using a manual Russian Corer with a diameter of 5 cm and a length of 50 cm. The chronology of the core was determined based on AMS 14 C dating. Pollen biostratigraphy was based on the results of radiocarbon dating analysis of sediment deposits and on the sequence of eight local pollen assemblage zone points from the Oldest Dryas to the Subatlantic (upper Meghalayan) periods (Fig. 2). Chironomidae subfossils were analysed at 8 cm resolution. Preparation methods for larval head capsules (HCs) followed Brooks et al. (2007), with sediments passed through a 63 lm sieve. As HC concentrations in the upper core sequence were low, kerosene flotation was used (after Rolland & Larocque, 2007). Taxon identification followed mainly keys by Klink & Moller Pillot (2003) and Brooks et al.  Pawłowski et al., 2016a) (2007), and ecological preferences of the identified taxa were based on Brooks et al. (2007), Vallenduuk & Moller Pillot (2007) and Moller Pillot (2009Pillot ( , 2013. [The reference collection is deposited at the Department of Invertebrate Zoology and Hydrobiology, University of Lodz.] Cladocera subfossils were analysed at 4 cm resolution. One cm 3 of fresh sediment was prepared in the laboratory following standard procedures (Frey, 1986), with each sample boiled in a 10% KOH solution for 20 min and stirred using a magnetic stirrer to remove humic matter, and later treated with HCl to eliminate carbonates. The residue was then washed and sieved (50 lm) and diluted in 10 cm 3 of distilled water. One tenth of a millilitre of solution was then used for each microscope slide. The extracted remains were identified after Szeroczyńska & Sarmaja-Korjonen (2007), Van Damme & Dumont (2008), Van Damme et al. (2010 and Faustova et al. (2011). All skeletal elements (i.e. head shield, shell, post-abdomen, claw, ephippium) in each sample were counted and used to compute the number of individuals per cm 3 of fresh sediment. Classification of Cladocera habitat preferences followed Flössner (1972,2000), Bjerring et al. (2009) and Błędzki & Rybak (2016).

Statistical analysis
The number of zones (and sub-zones) in the biostratigraphic sequence was determined after Bennett (1996), using optimal sum-of-squares partitioning (Birks & Gordon, 1985;Birks, 1986) using the ZONE software (Lotter and Juggins, 1991), with testing for statistical significance based on the broken-stick model (MacArthur, 1957) using the BSTICK software (unpublished software). To identify overall trends in assemblage composition, Detrended Correspondence Analysis (DCA) was performed on percentage squareroot transformed data using CANOCO 4.5 (ter Braak & Smilauer, 2002), with detrending by segments, nonlinear re-scaling, and down-weighting of rare species. Stratigraphic diagrams were created using C2 (Juggins, 2007). Rarefaction with ES(n) = 15 for Chironomidae and with ES(n) = 20 for Cladocera was computed using PRIMER 6 software (Clarke & Gorley, 2006).
Trends in the relative abundance and composition of Chironomidae and Cladocera guilds (i.e. habitat Fig. 2 Lithology, chronostratigraphy and representative plant and animal macrofossils and geochemical elements from the Pawłowa site (after Pawłowski et al., 2016a). Lithology: 1 = fluvial sands; 2 = gyttja; 3 = mid-decomposed peat; 4 = well-decomposed peat indicators) and taxa were further analysed by Minmax Autocorrelation Factor Analysis (MAFA: Zuur et al., 2007) to provide a more detailed description at the guild and (representative) taxon level. MAFA is a statistical method for estimating trends in multivariate (time) series related to principal component analysis with applications in ecology (Vilizzi, 2012;Whiterod et al., 2015;Shojaei et al., 2016;Doray et al., 2018) but not yet in palaeoecological studies. In MAFA, axes have maximum autocorrelation with time lag k. The first MAF axis represents the main trend in the data associated with the highest autocorrelation at lag 1, the second MAF axis at lag 2, and so forth, in decreasing order of importance. Cross-correlations (or canonical correlations) between variables (i.e. guilds and taxa) and trends (as indicated by the MAF axes) are then computed and tested for significance. Given that in MAFA the number of 'points' (in the present case, the zones and sub-zones under study) is constrained to be larger than the number of variables, a maximum of five guilds or taxa (i.e. the most abundant) could be analysed for both Chironomidae and Cladocera. Following square-root transformation (Chironomidae) and double-square-root transformation (Cladocera) and centring of the relative abundance (count) data, the first two MAF axes were estimated. Canonical correlations of the individual variables (i.e. guilds and taxa within guilds) with the two MAF axes of variation were then tested for significance (a = 0.05). Analyses were carried out in Brodgar v. 2.5.7 (http://www. brodgar.com).

Chironomidae
Sampling material of the Pawłowa sequence contained altogether 4460 HCs from 73 taxa of Chironomidae identified at least to the genus level, 619 HCs of unidentified Tanytarsini (excluding Corynocera ambigua Zetterstedt, 1837), 166 of Ceratopogonidae, and three of Simuliidae. Chironomidae stratigraphy was divided into six zones (Fig. 3).
Ch1 (450-408 cm): Oxbow establishment -Coldadapted taxa were present, with Corynocera ambigua representing 60% of the community at the top of this zone. Notably, at 444 cm a singular HC of Lipiniella spp. Shilova, 1961 (a taxon usually typical of large, sandy rivers) was found.
Ch3 (344-256 cm): Second well-developed lake phase -The relative abundance of Corynocera ambigua declined at the top of this zone from 70% to about 10% of the total and was replaced by a range of warm-adapted eutrophic taxa again mainly represented by Chironomus plumosus-type and Dicrotendipes nervosus-type.
Ch6 (152-0 cm): Dry terrestrial conditions -The number of specimens was overall very low. Besides Ceratopogonidae, a few amphibious taxa were identified.
Cl1 (  Cl2 (344-232 cm): Final stage of oxbow-lake development and transformation to mire -This zone was divided into two sub-zones. The first sub-zone Cl2a (344-300 cm), was characterised by an increase (especially at the beginning) in the pelagic pelagic Cladocera Eubosmina spp., which represented [ 40% of the specimens. Macrophyte/sediment-associated taxa including Alona affinis, Chydorus sphaericus and Coronatella rectangula were also quite abundant, and the macrophyte-associated Graptoleberis testudinaria (Fischer, 1848) and Sida crystalina (O. F. Müller, 1776) appeared for the first time, whereas sediment-associated taxa were still rare. The second sub-zone Cl2b (300-232 cm), was characterised by the highest frequency and diversity of Cladocera in the whole core, which exceeded more than 13 taxa and 7000 specimens per cm 3 . The dominant littoral macrophyte/sediment-associated taxa included Chydorus sphaericus and Coronatella rectangula, though littoral sediment-associated taxa and pelagic forms were still present, and sediment-associated Leydigia acanthocercoides (Fischer, 1854) occurred for the first time. At the end of this sub-zone, macrophyteassociated Alona guttata Sars, 1862 also occurred for the first time.

Trends in assemblage composition and relative abundance
For Chironomidae, the first DCA axis explained 17.9% of the cumulative percentage variance in the data (Fig. 3). After an increase from Ch1 to Ch2, values of the first axis of variation were high and did not change from Ch2 until the end of Ch3. At Ch4, values of the first axis decreased through to the middle of Ch6, and then sharply increased towards surface level. For Cladocera, the first DCA axis explained 30.2% of the cumulative percentage variance in the species data (Fig. 4). Along this axis, the highest values were achieved at Cl1 and Cl2, whereas at Cl2b values started to decline and were lowest at Cl3, and then increased moderately at Cl4. The first axis displayed peaks from the Oldest Dryas to the early Younger Dryas, reflecting the oxbow-lake development.
Species richness derived from rarefaction revealed different trends for Chironomidae and Cladocera. Chironomidae showed high rarefaction at Ch1, followed by a decrease to lower values at the Ch2/Ch3 boundary, and then by an increase at the Ch3/Ch4 boundary, which was followed by relatively high values at Ch4 (Fig. 3). Afterwards, at Ch5 and Ch6 the relative abundance was too small to count rarefaction for Chironomidae. Whereas, for Cladocera rarefaction could be computed for all four zones due to their higher relative abundance (Fig. 4). The highest values were identified at the bottom part of the sequence, increasing from Cl1 to Cl2b, and then decreasing at the Cl3/Cl4a boundary to reach the lowest value (0) at Cl4b and finally to increase slightly again.
For Chironomidae (noting that rheophiles were excluded from MAFA due to low relative abundances), there was a strong correlation with the first axis of variation of the oligotrophic, mesotrophic and eutrophic guilds (statistically significant for the last two) due to an overall decrease in relative abundance starting from Ch4 through to Ch6 relative to Ch1-Ch3 (Fig. 5). At the guild level, among the most abundant oligotrophic and cold-adapted morphotypes Chironomus anthracinus-type (sensu Brooks et al. 2007 after Chironomus (Chironomus) anthracinus Zetterstedt, 1860) and Corynocera ambigua were strongly correlated with the first axis (the former taxon also significantly) due to a sharp decrease in relative abundance starting from Ch4. And a similar trend was revealed for the most abundant mesotrophic and eutrophic taxa, though significantly only for the mesotrophic Microtendipes pedellus-type (sensu  (De Geer, 1776)). Conversely, no significant and clear trends were detected across the phytophile, rheophile and eurytopic guilds (Fig. 6). However, at the guild level the phytophile Glyptotendipes severini-type (sensu Brooks et al. 2007 after Glyptotendipes subg. Caulochironomus Heyn, 1993) (significantly), Ablabesmyia Johannsen, 1905 and Paratanytarsus penicillatus-type were highly correlated with the first axis due to a peak in relative abundance between Ch3 and Ch4, whereas Monopelopia tenuicalcar (Kieffer, 1918) was significantly correlated with the second axis as a result of a peak at Ch5. For Cladocera (Fig. 7), all guilds were highly correlated with the first MAFA axis of variation, and significantly so for the macrophyte associated and the pelagic guilds due to an overall decreasing trend starting from Cl3. This trend was overall reflected at the guild level, except for the sediment-associated, which were mainly represented by Pleuroxus trigonellus (O. F. Müller, 1785) and Pleuroxus uncinatus that significantly decreased in relative abundance only at Cl4. Finally, the macrophyte-associated guild was mainly represented by Acroperus harpae and Alonella nana (Baird, 1850), the macrophyte/sediment-associated by Alona affinis and Coronatella rectangula, and the pelagic by Eubosmina spp. Seligo, 1900 and Daphnia pulex-group.

Chironomidae
Mean summer temperature and water trophic status are known to be main drivers of Chironomidae community composition and relative abundance in temperate climates (Brodersen & Quinlan, 2006;Brooks, 2006), though some taxa have specific preferences for certain mesohabitats, hence eluding general classification (Brooks et al., 2007). In the present study, the overall change in Chironomidae communities reflected the climate-related Late Glacial-Holocene transition from oxbow to mire environment of Pawłowa.
Overall, the taxonomic composition of Chironomidae communities suggests that during the Late Glacial Pawłowa was a mesotrophic water body, with all trophic groups of taxa still finding appropriate habitat conditions in the different zones of the lake. Although the trophic status of Pawłowa lake in the Late Glacial was strongly related to temperature, it also depended on climate-related variables such as denudation processes and flooding. In this respect, in-wash of nutrients and inorganic particles to Pawłowa lake has been confirmed by geochemical, palaeobotanical and Cladocera proxies (Pawłowski et al., 2016a), and similar processes have been reconstructed from the Younger Dryas sequence of Lake Ś wierczyna in the Grabia valley (Pawłowski et al., 2015a), where mean annual temperature was lower during the Late Glacial compared to the Holocene. This explains why the first Chironomidae generation in Pawłowa lake has been found to include cold-adapted taxa despite high summer temperatures (Brodersen & Lindegaard, 1999;Kołaczek et al., 2018b). Whereas, low annual temperatures prevented fast eutrophication of the small Pawłowa pond, when nutrients could only be assimilated by the biota during the short and warm summers, with the result that Pawłowa lake maintained rich benthic assemblages during the Late Glacial. Paludification then started at the end of the Younger Dryas and occurred in parallel with eutrophication but caused a decline in most species from the interstadial generation. Palaeobotanical proxies indicate high trophic status and slightly alkaline conditions in Pawłowa lake during the early Younger Dryas, although during its second part and at the onset of the Holocene the lake transitioned to a nutrient-poor, acidified water body (Pawłowski et al., 2016a).
Amongst the groups of taxa associated with specific mesohabitats, phytophile taxa revealed a clear pattern in assemblage composition in the present study (Fig. 7). Plant macrofossil analysis (Pawłowski et al., 2016a) indicates the development of rich aquatic vegetation (i.e. Chara spp. Linnaeus, 1753 and Potamogeton spp. Linnaeus, 1753) during the Younger Dryas and of amphibious plant vegetation (Carex spp. Linnaeus, 1753) during the Younger Dryas-Holocene transition (Fig. 2). Since the Greenlandian, the lake transformed into a mire with ample Bryales communities (Fig. 2), which could maintain abundant Chironomidae populations in the valley wetlands (Brusven et al., 1990). In this regard, Monopelopia tenuicalcar has often been recorded in peat sections of palaeolakes of the region (e.g. Płóciennik et al., 2011;Pawłowski et al., 2015a;Słowiński et al., 2016), and this is known to be a species characteristic of shallow, acid and overgrown peat pools elsewhere in Europe (Vallenduuk & Moller Pillot, 2007).
Rheophile taxa, albeit rare, suggest episodes of connection of the Pawłowa lake to the main River Grabia channel (Howard et al., 2010). Thus, Lipiniella spp. included taxa typical of large rivers with a sandy bottom and low silt/clay content, and the presence of a singular head capsule of Lipiniella spp. in the lowest

Canonical correlations
Pawłowa layers suggests that, during the Oldest Dryas, the River Grabia channel was large, sandy and meandering. The Lipiniella spp. head capsule detected in the present study could belong to a specimen inwashed to Pawłowa lake from the adjacent river during a flooding episode, and a small population of this taxon could have developed in the lake if habitat conditions resembled those of the river channel (Moller Pillot, 2009;Klukowska et al., 2011). The second generation of rheophile guilds appeared in Pawłowa at Younger Dryas-Holocene transition, with Simuliidae clearly associated with small-medium streams (Bass, 1998;Palmer et al., 2002). The parallel presence of a few rheophile species in the same zones indicates higher river activity and a periodical connection of the River Grabia to the Pawłowa mire during the Greenlandian. The most prominent floods were recorded in the Grabia valley during the Younger Dryas, when major flood flows are estimated to have been at least five times larger than present-day ones (Pawłowski et al., 2015a(Pawłowski et al., , 2016b, and higher river activity in the Younger Dryas and the Greenlandian has also been recorded from other places of the study region (Starkel et al., 2013;Kittel et al., 2016). Also, during the Greenlandian the Pawłowa mire experienced phases of substantial floods from the River Grabia, which however had different and less prominent sedimentological impact than the floods recorded from the Younger Dryas layers (Pawłowski et al., 2016a). The 'Preboreal oscillation' cooling then caused increased climate instability and led to a higher rate of extreme events including floods; whereas, the presence of more developed riparian and swamp vegetation might have limited the impact of floods on the oxbow lakes of central European valleys (Kołaczek et al., 2018a). Research in this regard has shown that riparian wetlands freely connected to the main channel are characterised by a higher biodiversity and greater biomass of macrobenthos (Pan et al., 2012), and that reconnection of an oxbow to the main river channel may cause an increase in Chironomidae relative abundance and benthic community diversity (Obolewski & Glińska-Lewczuk, 2011). Therefore, stronger fluvial activity of the River Grabia during the Younger Dryas and the early Greenlandian could have prolonged the presence of more diverse communities, including fully aquatic taxa even if not abundant. Starting from the late Greenlandian, when river activity decreased, aquatic species were replaced solely by semi-terrestrial Chironomidae taxa. Despite the absence of trends in amphibious guilds (likely due to low subfossil densities), the stratigraphic data of the present study clearly revealed that amphibious Chironomidae almost exclusively dominated in the Holocene peat sequence of the Pawłowa mire. This is not surprising, as the Greenlandian oxbow paludified to a semi-terrestrial habitat, and semi-terrestrial assemblages are known to occur generally at low to very low subfossil densities (e.g. Pawłowski et al., 2015a;Kittel et al., 2016), with Chironomidae being more typical of shallow lakes (Whitehouse et al., 2008). Starting from the Northgrippian, seasonal peat desiccation and decomposition negatively influenced the preservation of Chironomidae subfossils, and even if terrestrial taxa were present, their head capsules were partly decomposed during the Northgrippian and Meghalayan.

Cladocera
Changes in Cladocera composition and relative abundance are generally related to water level fluctuations, with the formation and functioning of lakes on flood plains being dependent on the hydrological regime of the associated rivers combined with physical and chemical parameters (Amoros & Jacquet, 1987;Korhola & Rautio, 2001). There are also several environmental variables (e.g. summer water temperature, fish predation, macrophytes, nutrients, and sediment in-wash) that have been identified as important in influencing these aquatic crustacean assemblages (Medley & Havel, 2007;Schöll & Kiss, 2008;Vadadi-Fülöp, 2013;Pawłowski et al., 2015b).
Pelagic and macrophyte-associated taxa had a peak in relative abundance from the Oldest Dryas to the early Younger Dryas after which they declined. The decrease in the relative abundance of these taxa may be related to a systematic water level drop with transition of Pawłowa from lake to mire. The oligo/ mesotrophic conditions during the first stage of Pawłowa lake also would have favoured the abovementioned taxa, which however are characterised by a wide range of habitat preferences (Duigan, 1992). For example, Acroperus harpae and Alonella nana can occur in the rheopelagial as potamoplankton (Robertson, 1990;Błędzki & Rybak, 2016), and Acroperus harpae can act as a 'pioneer' taxon with early immigrants following ice retreat (Hofmann, 2000;Pawłowski, 2011;Błędzki & Rybak, 2016). Therefore, a gradual drop in dominance in both guilds may also indirectly reflect climate changes in the River Grabia basin with increased summer water temperatures.
Amongst the taxa associated with sediments, only two species (i.e. Pleuroxus trigonellus and Pleuroxus uncinatus) revealed a clear pattern in the assemblages of the Pawłowa palaeolake with relatively high dominance from the Oldest Dryas to the early Younger Dryas. The high dominance of these benthic taxa has been correlated to the geochemical characteristics of Pawłowa (Pawłowski et al., 2016a), and especially with an increase in the admixture of mineral matter (SiO 2 ) and presence of detrital components (i.e. K, Al, Na, Ti, and Mg). This can be explained as the result of higher river activity and a periodical connection of the River Grabia to the Pawłowa palaeo-oxbow but also of denudation processes in the catchment. The distinct decrease in sediment-associated taxa starting from the end of the Greenlandian was then a response to the weaker fluvial activity of the River Grabia and the transformation of the oxbow-lake to mire.
Amongst the macrophyte/sediment-associated taxa revealing clear patterns, Alona affinis has highly adaptive traits, including the ability to live in lakes of various depths and trophic state and feed on detritus under extreme conditions (Flössner, 2000). Coronatella rectangula also has high ecological plasticity (Flössner, 2000) and is an important indicator of increased in-lake productivity and warm waters in the Pawłowa palaeolake, where its high relative abundance in the Younger Dryas-Greenlandian transition can be linked to higher trophic state and a gradual increase in temperature. The following sharp decrease in the relative abundance of this species is correlated with the presence of Bryales communities but could also be attributed to low water pH and the transformation of the Pawłowa lake into a semi-terrestrial habitat. Finally, it is also relevant that in zones Cl2a and Cl2b fish remains were collected from the core sediment. Fish generally affect zooplankton communities, as size-selective predation on larger Cladocera could increase the relative abundance of some other planktonic species (Davidson et al., 2007;Feniova et al., 2015).
Autecology and diversity as determinants of fossil records Chironomidae and Cladocera have different life cycles, dispersion mechanisms and resistance to environmental factors such as temperature and hydrology. These autoecological factors were the main drivers of their different responses to the global changes that took place at Pawłowa during the Late Glacial and Holocene. The global factor driving the Grabia valley ecosystems in the long timescale was climate change, which was mainly related to temperature fluctuations and duration of ice cover. However, both Chironomidae and Cladocera also responded to local-scale abiotic factors including water trophic status, hydrology, pH and paludification, and to biotic factors related to aquatic vegetation and fish predation.
Temperature fluctuations were the main driver of freshwater biocenoses during the Late Glacial (Nevalainen et al., 2012;Luoto et al., 2019) and the onset of the Holocene (Kołaczek et al., 2018b). In temperate zones, Chironomidae are regarded as being related directly to air temperature, mostly because of their terrestrial adult stage. In this respect, adult emergence in univoltine species has been found to correlate to higher (summer taxa) or lower (species that emerge only in early spring) day temperatures (Płóciennik et al., 2018). Summer day length influences the rate of Chironomidae larval development, thereby affecting relative abundance and number of generations per year. Chironomidae are fast reproducers and their voltinism may change annually (Brooks, 2006;Armitage et al., 1995). Also, the Chironomidae fauna of glacial lakes in cold climate regions appears to be more sensitive to local continentalism, being composed of profundal, oligotrophic and cold-stenothermic taxa; whereas, floodplain assemblages of the boreal zone are closely related to July air temperatures and reveal stronger habitat influence (Nazarova et al., 2017). Cladocera, albeit generally more sensitive to changes in water trophic status, are also influenced by temperature (Lotter et al., 1997;Korhola, 1999). In this regard, parthenogenetic reproduction in Cladocera is a response to seasonal weather changes, with a quick switch to sexual reproduction and production of resting eggs if temperature conditions become unfavourable (Brendock & De Meester, 2003). As the Pawłowa oxbow was a shallow water body during the Late Glacial, summer temperatures should have acted as main driver in the assemblage composition and relative abundance of both Chironomidae and Cladocera. However, reconstructions of July air temperatures from Pawłowa show different results, with Chironomidae associated to substantially higher (& 2-3°C) temperatures compared to Cladocera (Pawłowski et al., 2016a;Płóciennik et al., unpublished data). Also, climatic reconstructions from Ś wierczyna (see above) have pointed to the importance of assemblage diversity for reliable palaeolimnological reconstructions (Pawłowski et al., 2015b), which however are often confounded by the presence of Corynocera ambigua -a common dominant species in Late Glacial lake sequences of Europe (Brodersen & Lindegaard, 1999).
Another crucial climatic aspect of Pawłowa lake during the Late Glacial-Holocene transition, and one that has been rarely accounted for, is the duration of ice cover, which can affect the phenology of both Chironomidae and Cladocera. Chironomidae seem to be more dependent on mean July temperatures than on level of continentalism, although ice cover influences the survival of their wintering larvae . During winter, these remain in diapause and may suffer from factors indirectly related to low air temperatures and insolation level, including limited availability of food resources and duration of the icefree period (Nazarova et al., 2017). On the other hand, Cladocera produce ephippia, which are highly resistant to freezing and desiccation (Fryer, 1996) and in cold climates their larger numbers (cf. chydorids) are generally found in lake sediments as the relative duration of the asexual reproduction period is limited in the presence of an open-water shorter season (Sarmaja-Korjonen & Seppä, 2007;Frisch et al., 2012;Pawłowski, 2017).
Although climate is a leading factor at northern latitudes, temperature combines with water trophic status and hydrology in affecting both Chironomidae and Cladocera assemblages (Brown, 2007;Schöll & Kiss, 2008;Kotrys et al., 2020). Previous research in the Grabia valley indicates that the valley's peripheral stagnant water bodies (i.e. oxbow lakes) had a much more complex ecological profile compared to postglacial lakes. Chironomidae and Cladocera are reliable indicators of episodic hydrological fluctuations (Pawłowski et al., 2015a(Pawłowski et al., , 2016b, and temporary connections of the floodplain water bodies (such as Pawłowa) with the main channel of the River Grabia had a positive influence on the life condition of these aquatic invertebrates. This is because connectivity with the river channel would favour water exchange in the oxbow, thereby increasing the diversity of plantassociated invertebrate communities. Collector-filtrators (hence, including Cladocera) are known to be dominant during stagnant phases, whereas in presence of water movement collector-gatherers (mostly, Chironomidae) play an important role. Limited water exchange between an oxbow and the main river channel would lead to unfavourable conditions for most macrobenthos and meiofauna species due to a decrease in water aeration and light penetration occurring at the onset of paludification processes (Obolewski, 2011). In this regard, aquatic insect subfossil assemblages (including mainly Chironomidae) have provided useful information regarding habitat characteristics and flow activity in palaeochannels, even in a semi-quantitative way (Gandouin et al., 2006;Howard et al., 2010). For Cladocera, the stability of Pawłowa lake was crucial for dominant macrophyte/sediment-associated taxa. In general, microcrustaceans are more exposed to the effects of floods than benthic macrofauna, with Chironomidae larvae able to migrate into the hyporheic zone. Flooding can also be critical for populations located in oxbows adjacent to the river channel or during exceptional overbank flows (Kołaczek et al., 2018a), and that is the reason why Chironomidae provided more fragmentary information about fluvial activity in the River Grabia. Well-developed microcrustacean potamoplankton is an accidental phenomenon (Illyová et al., 2008;Czerniawski & Domagała, 2012), and overbank flood episodes in Pawłowa lake have been also indirectly supported in the present study by geochemical profiles and sedimentology. During floods, there were episodic increases in sediment associated Cladocera and in fully planktonic species, and this pattern in assemblage composition occurred especially during cold climate events. Additionally, these changes were displayed by the higher rarefaction observed in the present study (Fig. 4). It is therefore argued that the increase in Cladocera species richness and relative abundance can be related to enhanced eutrophication during floods, which are considered one of the most important factors contributing to the development of aquatic biota.
Paludification, which causes pH decrease and terrestrialisation during dry phases, is also known to strongly affect macroinvertebrate communities (Verdonschot, 1992). During the Late Glacial, in Pawłowa lake there were slightly alkaline water conditions (Pawłowski et al., 2016a), similar to several shallow water bodies and mires across Poland, although at the end of the Younger Dryas and onset of the Holocene pH generally decreased (i.e. Kołaczek et al., 2015;Forysiak, 2018). Water pH below 6.0 is known to eliminate many Chironomidae taxa, with non-biting midges replaced by biting ones (i.e. Ceratopogonidae) (Henrikson et al., 1982;Walker et al., 1985;Brooks et al., 2007). Generally, paludification was one of the dominant factors in stagnant waters of mid-European lakes since the Northgrippian, and Chironomidae are regarded as a group of desiccation-resistant organisms that spend dry periods in the soil layer (Batzer & Boix, 2016). Current ecosystem studies indicate that wetland and terrestrial Chironomidae communities may be diverse and abundant (Rosenberg et al., 1988;Frouz, 1999), although species of common terrestrial genera (e.g. Limnophyes spp. or Pseudorthocladius spp. Goetghebuer, 1943) can be usually distinguished only at their adult stage. Cladocera remain far more abundant, although in the present study they were represented by few taxa of which Chydorus sphaericus is the most resistant -a species that appears in abundant populations even in ephemeral pools (Płóciennik et al., 2015a).
A biotic factor that influenced invertebrate communities in the Pawłowa palaeo-oxbow is aquatic vegetation. In the Greenlandian, phytophile Chironomidae taxa may have been predominant (Kołaczek et al., 2018a), but during the Late Glacial this role was overtaken by sediment-associated guilds , and this was also observed in Pawłowa. Whereas, macrophyte-and macrophyte/sediment-associated Cladocera were a dominant component in Late Weichselian assemblages. As Pawłowa was a relatively shallow water body, plants including Characeae and Potamogeton spp. were abundant along the palaeochannel bank zone (Pawłowski et al., 2016a), and abundant vegetation is known to limit resources (e.g. light, nutrients) for phytoplankton and planktonic crustaceans (Kołaczek et al., 2015). A high rate of Ca:Mg has also been observed at Pawłowa (Pawłowski et al., 2016a), causing precipitation of carbonates on plant leaves thereby increasing surface roughness and strongly enhancing periphyton production (Obolewski, 2011). Chironomidae can also benefit from sediment re-suspension overload as a result of denudation processes, as observed during the Bøling-Allerød Interstadial and Younger Dryas (Pawłowski et al., 2016b) -even though Chironomidae have been shown to prefer aquatic vegetation such as Fontinalis spp. Hedwig, 1801 moss patches over inorganic sediment microhabitats (Brusven et al., 1990). In lakes with rich vegetation abundance and diversity, biotic gradients including bentho-pelagic production, fish predation, plant substratum, food availability and competition override the indirect mechanisms associated with nutrient availability Sayer et al., 2010). The quantitative and qualitative structure of macrophyteassociated invertebrate communities is driven, to a great extent, by the hydrological connectivity between an oxbow lake and the parent river channel (Obolewski, 2011), and this is why phytophile and rheophile Chironomidae may appear at the same time during the Younger Dryas-Holocene transition.
Fish predation, a second biotic factor, is considered as equally strong factor as hydrology in shaping Cladocera communities (Pawłowski et al., 2016a). A slight increase in zooplanktivorous fish density in Pawłowa lake during the late Allerød and Younger Dryas resulted in the absence of the larger taxa, especially from the Daphniidae family, paralleled by a relative increase in the relative abundance of species from the Bosminidae family (Figs. 2 and 4). Experimental studies have shown that increasing density and dominance of smaller-sized fish due to water warming conditions would augment predation pressure on large zooplankton species, favouring the development of small-sized ones (Meerhoff et al., 2007). Similarly, a shorter ice-cover period caused by increasing warming conditions in European lakes might be correlated with increasing survival of young fish during winter. images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.