Drivers of Macrophyte and Diatom Diversity in a Shallow Hypertrophic Lake

: We studied macrophyte and diatom assemblages and a range of environmental factors in the large hypertrophic Deht á ˇr ﬁshpond (Southern Bohemia, Czech Republic) over the course of several growing seasons. The spatial diversity of the environment was considered when collecting diatoms and water samples in three distinct parts of the ﬁshpond, where automatic sensor stations continually measuring basic factors were established. Macrophytes were mapped in 30 segments of the ﬁshpond littoral altogether. High species richness and spatiotemporal variability were found in assemblages of these groups of autotrophs. Water level ﬂuctuations, caused by the interaction of ﬁsh farming management and climatic extremes, were identiﬁed as one of the most important factors shaping the structure and species composition of diatom and macrophyte assemblages. The distance of the sampling sites from large inﬂows reﬂected well the spatial variability within the ﬁshpond, with important differences in duration of bottom drainage and exposure to disturbances in different parts of the ﬁshpond. Disturbances caused by intensive wave action are most probably a crucial factor allowing the coexistence of species with different nutrient requirements under the hypertrophic conditions of the Deht á ˇr ﬁshpond. Due to a range of variables tested and climatic extremes encountered, our study may be considered as a basis for predictive model constructions in similar hypertrophic water bodies under a progressing climate change.


Introduction
Shallow lakes are often considered unique habitats of rich aquatic and wetland biota [1][2][3]. Increasing attention is paid particularly to natural oligotrophic and mesotrophic lakes, which started to be rare due to overall eutrophication, or, in some regions, acidification [4,5]. On the other hand, heavily eutrophicated water bodies usually become the focus of scientists as the objectives of restoration projects [6,7]. Many are considered sites where a substantial part of biodiversity has already been lost and where mitigation measures are needed in order to turn back this unfavourable trend [1,8]. Many fishponds of these more than 20,000 artificial water bodies occurring in the Czech Republic, a country that lacks natural lakes, belong among aquatic habitats with the status of high eutrophy to hypertrophy [9]. Since their construction, mainly in the Middle Ages [10], most of these fishponds have served for semi-intensive carp breeding and recently to a certain extent also for other Our earlier studies [16,55] showed that Dehtář is a strongly eutrophic to slightly hypertrophic fishpond (according to the OECD classification of 1982 [65], adjusted for fishpond conditions [66]). Similarly to the majority of Czech fishponds, it is a carp pond, i.e., serving mainly common carp (Cyprinus carpio) breeding and to a lesser extent also other fish species with similar habitat requirements as the carp (more details are included in Figures S1-S15). Dehtář is a so-called main fishpond, i.e., the production of marketable fish (fish of consumable size) is its farming purpose. Like the majority of large fishponds in the Czech Republic, Dehtář is managed in bi-annual management cycles. In the first year of the management cycle, the fishpond is stocked with young fish; as the whole fish biomass is, due to a smaller size of fish individuals, much lower than in the second year of the management cycle, the pressure of fish stock on the aquatic environment, particularly on its vegetation, is lower as well [19,20]. In the second year of the management cycle, the fish biomass strongly increases; the whole management cycle ends in October/November by harvesting the fish stock. Before the harvesting, a substantial part of the fishpond is drained (Figure S1), and although it is flooded again after the fish harvesting, the water level is kept lower during the first year of the management cycle than it is in the second year. Therefore, the values of other fishpond parameters (the total area, total volume, and average depth) in the first year should also be lower than in the second year and increase only slowly. Our earlier studies [16,55] showed that Dehtář is a strongly eutrophic to slightly hypertrophic fishpond (according to the OECD classification of 1982 [65], adjusted for fishpond conditions [66]). Similarly to the majority of Czech fishponds, it is a carp pond, i.e., serving mainly common carp (Cyprinus carpio) breeding and to a lesser extent also other fish species with similar habitat requirements as the carp (more details are included in Figures S1-S15). Dehtář is a so-called main fishpond, i.e., the production of marketable fish (fish of consumable size) is its farming purpose. Like the majority of large fishponds in the Czech Republic, Dehtář is managed in bi-annual management cycles. In the first year of the management cycle, the fishpond is stocked with young fish; as the whole fish biomass is, due to a smaller size of fish individuals, much lower than in the second year of the management cycle, the pressure of fish stock on the aquatic environment, particularly on its vegetation, is lower as well [19,20]. In the second year of the management cycle, the fish biomass strongly increases; the whole management cycle ends in October/November by harvesting the fish stock. Before the harvesting, a substantial part of the fishpond is drained ( Figure S1), and although it is flooded again after the fish harvesting, the water level is kept lower during the first year of the management cycle than it is in the second year. Therefore, the values of other fishpond parameters (the total area, total volume, and average depth) in the first year should also be lower than in the second year and increase only slowly.
Strips of the fishpond bottom are usually exposed in the first year of the management cycle, offering a suitable habitat to a range of wetland vascular plant species and their  (Figures S2-S12). Large areas of shallow water support the growth of some amphibious and aquatic plant taxa. During our research period, however, the low water level from the first year of the management cycle in 2015 continued until May 2016, i.e., the second year of the management cycle ( Figures S11 and S12). This was a consequence of extremely dry weather in 2015 and in the spring of 2016, enhanced by extraordinarily high temperatures in the summer of 2015 (for basic climatic characteristics of the Dehtář fishpond from the study period compared with the climate of the whole region, see Table A1). This dry period ended in heavy rains in the summer of 2016 ( [17]; Table A1, Figure S13), which elevated the water level above its common management value. Moreover, 2014, also the 2nd year of the management cycle, was also characterised by a lower water level due to previous dry years ( Figure S11). According to our information from fish farmers [67], a lack of water was always quite frequent in the Dehtář fishpond; for this reason, the exposed bottom used to be sown by clover and other crops in the 1960s-1970s in order to get forage for domestic animals and to increase the fishpond fertility (see also [68]).
Semi-intensive fishpond management practiced on the Dehtář fishpond (and the majority of fishponds throughout the Czech Republic) includes fertilisation with manure, supplementary feeding of fish stock with cereals, and occasional liming of some fishpond parts [68]. Another source of nutrients is the surrounding agricultural landscape and the water brought by the inflows, particularly the large ones, the Babický stream, the Kamenný stream, and the Dehtářský stream.
As already reported elsewhere (e.g., [16,38,39,55]), hypertrophy combined with the relatively high average depth of the Dehtář fishpond contributes to strong fluctuations of oxygen concentrations and other water parameters. Another important factor is strong winds, mainly from the northwest, causing intensive wave action and associated mechanical disturbances in the littoral zone ( Figure S14). Thus, despite the cyanobacterial blooms appearing in the last years ( Figure S15), Dehtář is popular as a recreational fishpond, particularly among windsurfers [9,55].
Dehtář is a fishpond with several bays separated by natural barriers such as peninsulas and an island. Its surroundings include large areas of arable land and cattle pastures, and small patches of wet to semi-dry meadows, woodlands, and settlement areas. As the natural barriers and diversity of the surrounding landscape may cause differences in some environmental parameters as well as in the biota between the different fishpond parts, for the aim of this study, we divided the fishpond into a total of 30 segments (Figure 1), expecting more or less homogeneous environmental conditions within each of them. The length of the segments varied according to local topography of the fishpond littoral and the surrounding landscape, and their area (=vegetated zone) was also dependent on water level fluctuations determining the extent of exposed pond bottom and shallow water zones (see Table S1 for more details).
Additionally, three parts of the Dehtář fishpond with presumably the largest differences in environmental conditions were selected for the research of the physico-chemical parameters of water and epiphytic diatom assemblages. All these sites have been equipped with automatic, sensor-based stations designed for high-frequency monitoring of several basic environmental parameters; a range of additional parameters was recorded at regular intervals at the stations (for detailed description, see the Methods). The sites are, namely, (1) the fishpond dam-FD, (2) the bay of the Dehtářský stream-DSB, and (3) the bay of the Babický stream-BSB ( Figure 1). Epiphytic diatom assemblages were studied in the vegetation zones as close as possible to the three monitored sites, about 200 m from the stations in the Dehtářský stream bay (DSB) and Babický stream bay (BSB) and about 600-650 m from the FD station. The sampling site at FD was situated in a narrow bay at the northeastern corner of the fishpond, naturally preserved against waves (segment 2 near the border with segment 3, and segments 3 and 4; Figures S2 and S3). In contrast, BSB and DSB sampling sites were in broad bays, open and fully exposed to wave action. While the fishpond bottom in the DSB and especially in the BSB exhibits moderate slope and thus remains non-flooded for a relatively long time after the fish harvesting (Figures S11 and S12), Water 2021, 13, 1569 6 of 39 the site near the fishpond dam is characterised by a relatively steep slope and very fast water level increase.
Since the 1930s, the Dehtář fishpond has attracted the interest of botanists (e.g., [69,70]). It is also an important ornithological locality, protected as a part of the Natura 2000 network [40,71]. The earlier research of Dehtář was, however, not systematic. Thus, only individual records on plant and animal species are reported in the literature.

Physico-Chemical Parameters of Water
In 2014, three monitoring sensor-based stations were established in three different parts of the Dehtář fishpond (for the location and description of these sites see Section 2.1. Study Site and Figure 1). Each station was equipped with several sensors. The following parameters were measured: water temperature and concentration of dissolved oxygen (M4016, Fiedler AMS) at all stations, global radiation, rainfall, speed and direction of wind (M4016-A-G3, Fiedler AMS) at the BSB station, and photosynthetically active radiation (PhAR; LI-COR sensors) at the FD station (further details on the stations and parameters measured are included in Figures S16-S22). At the end of 2015, the station at FD was further equipped with a sensor for water level detection. After 2016 (end of the targeted research project and of the substantial part of funding), only the station at FD was kept, and some sensors from the other stations were moved there.
In addition to the measurements provided by the sensor system, basic physicochemical parameters (temperature, concentration and saturation of oxygen, pH, conductivity, and water transparency) were measured by the multiparametric probe YSI 6600 V2, and samples for water chemistry were taken at bi-weekly intervals near the stations (some of the parameters were taken only at FD).
All samples were analysed for total (TP) and soluble reactive phosphorus (SRP); total (TN), ammonia (N-NH4), and nitrate (N-NO3) nitrogen; total (TC), total organic (TOC), and dissolved organic carbon (DOC); dissolved ions of sodium (Na), magnesium (Mg), calcium (Ca), and potassium (K); and total iron (Fe) according to standard methods certified and intercalibrated in the laboratories of the Vltava River Authority, state enterprise. Moreover, the concentration of chlorophyll-a (Chl-a), and alkalinity (AT) were also analysed in epilimnetic samples (integrated samples from 0-1 m water layer). Water transparency (WT) was measured simultaneously. An overview of the physico-chemical characteristics obtained for the growing seasons (i.e., IV-VIII) of 2014-2016 and 2019 is given in Table S2.

Other Environmental Factors, Fishpond Management Data
A range of environmental factors, related to vascular plants, was recorded for each segment directly in the field. These were, e.g., type of substrate, depth of mud, occurrence of stones, or organic detritus (see Table S1 for the complete list of factors). Additionally, the length of segments, their area in individual years, the proportion of the landscape types surrounding each of the segments, and other characteristics (Table S1) were measured in the national internet mapping application (www.mapy.cz, accessed on 15 November 2020). Fishpond management data on the fish stock composition and biomass, amounts of manure, lime, and supplementary feeding applied were provided by the Hluboká nad Vltavou Fish Farm Cz, s.r.o. However, these data were used only for interpretations of some parameters studied, as they usually could not be linked to particular segments of the fishpond (see Table S1 for more details).

Vascular Plant Data Recording
Vascular plant and charophyte species were mapped in the field separately for each of the 30 segments using a three degree scale: 1-rare (up to ten individuals in a segment), 2-frequent throughout the segment but not dominant, and 3-dominant, forming its own stands. The data were collected throughout the growing season of a particular year in monthly intervals (IV-IX 2014, III-IX 2015, III-IX 2016). Segments 1-19 were mapped  during each visit, while segments 20-30, mostly difficult to access, were visited twice  a year in 2014 and 2015 and only once in 2016. In 2019 (after the end of the research  project), only segments 1-19 were mapped during three sampling events in May, June, and August. Taxonomically critical taxa and/or difficult to identify taxa were collected into a herbarium and revised by experts (see the Acknowledgements). The herbarium vouchers are deposited in the herbaria PR and BRNU (the abbreviations follow Thiers [72]).

Diatom Sampling and Processing
The sampling was performed during three one-day sampling events (16 June 2015, 23 June 2016, and 19 June 2019). Macrophyte species, water depth, distance of a sample from the bank, and disturbance category were noted for each of the samples directly in the field (for more details see Table S3 with the list of samples).
Epiphyton was collected with the whole substrate, which means that the plants were cut with scissors above the bottom and rinsed in water to wash away nonepiphytic species. The samples for light microscopy analysis (LM) were placed into polyethylene bags, kept in a refrigerator, and transported to the laboratory. For Environmental Scanning Electron Microscopy (ESEM), up to 10 cm long sections of stalks and leaves of Phalaris arundinacea and Phragmites australis in 2016 and Limosella aquatica in 2019 were cut, rinsed in water, attached to a polystyrene matrix in order to keep diatom assemblages in the state best corresponding to natural conditions, placed into plastic containers containing water from the locality, and transported to the laboratory in a cooling thermo box.
In order to get better insight into the structure of epiphytic diatom assemblages collected on sites with contrasting environmental conditions and/or different macrophyte taxa, observation using ESEM was employed. Sample observation and photography in 2016 took place at the Institute of Scientific Instruments of the Czech Academy of Sciences in Brno. The Low Temperature Method for sample stabilisation in ESEM was used. This method was already successfully applied to diatom assemblages [84]. Observations were performed on an ESEM QUANTA 650 FEG at a temperature of -20 • C. Sample observation and photography in 2019 were conducted at the Institute of Botany in Brno on a table scanning electron microscope Phenom ProX using a temperature-controlled sample holder at a temperature of -8 • C.

Data Processing and Analyses
Data from the automatic sensor stations, recorded at 10 minute intervals, were downloaded from the hosting server and stored regularly. They were checked for possible inconsistencies caused by station failures (e.g., during extreme storms) or by removal of biofilms from some of the sensors. The potentially erroneous data were removed. In this paper, we analysed data on the following factors: air temperature, precipitation, wind speed, water temperature, oxygen concentration, and photosyntetically active radiation (PhAR), while the data from the sensors for global radiation and water level fluctuation data recording were used only in descriptions of site conditions. After the prescreening, only the data measured by the sensors at the depth of 0.3 m were selected and further elaborated. We used absolute values of wind speed, air temperature, oxygen concentrations, water temperature, and PhAR and displayed them as line graphs in order to compare the course and fluctuations of these parameters between years. We also prepared short graphs with a detailed visualisation of wind speed, oxygen concentrations, and PhAR at the time of diatom sampling, as a several-day course of these factors may potentially impact fast-changing diatom communities [85]. Additionally, for air temperature and precipitation, we calculated monthly means and sums, respectively. For the months that were missing or incomplete in our own data, we used publicly available data from the nearest weather station of the Czech Hydrometeorological Institute (CHMI) inČeské Budějovice-Rožnov [86]. In order to highlight climatic extremes during the study period, we also used CHMI data [86], where long-term means for the whole South-Bohemian region are available. To get better insight into the frequency of wave action in individual years, we classified each wind speed record in one of the classes defined according to the Beaufort wind force scale (https://cs.wikipedia.org/wiki/Beaufortova_stupnice, accessed on 1 March 2021) and calculated the proportion of winds of each class in each year (Table 1). Table 1. Frequency of selected wind events based on the data on wind speed from an automatic sensor station (10 min measurement interval). Absolute counts of each event in each year, along with their share in total numbers of records, are given. All other data were computerised and prepared for statistical analyses (vascular plants and charophytes, diatoms, environmental and management-related data from the field, data from the maps, and average monthly chlorophyl-a concentrations as a proxy of nutrient concentrations and water transparency) and/or transformed into the tables for direct presentations (all physico-chemical water parameters, vascular plant, charophyte, and diatom abundances and frequencies). Vascular plants were classified according to their origin [87] and threat status (Grulich [88], with the updates by Grulich & Chobot [89]). Additionally, the classification of vascular plants and charophytes into functional groups was performed according to Francová et al. [19,20], with an extension to the species not presented in the cited papers (Table S4). An adjusted version of this classification was also used in macrophytes collected for diatom analyses (Table S3). In order to get a better overview of the terrestrial species imported to the fishpond, the ecological classification was performed in parallel, based on our own field experience and species ecological characteristics summarised within the PLADIAS database (www.pladias.cz, accessed on 20 November 2020; [90]). The nomenclature of vascular plants follows Danihelka et al. [91], with the exception of Spergularia kurkae, where we accepted the concept of Kúr et al. [33]. The nomenclature, threat status, and origin of charophytes are based on the charophyte overview by Caisová and Gąbka [92]. All the classification categories for each species are given in Table S4.

Wind
Based on the diatom relation to the host plant surface (according to Fránková et al. [55]), five groups, according to diatom functional types (FT1-FT5) were used: FT1-planktonic taxa represented by centric diatoms; FT2-typically periphytic taxa adhering to the surface directly by a mucous film or with a mucilaginous stalk; FT3-facultatively periphytic araphid, passively moving diatoms able to attach; FT4-facultatively periphytic taxa with raphe with fibulae, actively moving; and FT-epipelic taxa with raphe, actively moving, mainly symmetrical biraphid pennate diatoms (see Table S5 for more details). The affiliation of diatom species to functional types is given in Tables S4 and S5. Planktonic and epipelic taxa were present in the diatom samples even though the collected macrophytes were washed in order to remove non-epiphytes. It was decided to keep and classify the whole diatom assemblage in order to get a realistic view of the community.
The patterns in species composition of both diatoms and macrophytes were investigated by non-metric Multidimensional scaling (NMDS). The data in each dataset were standardized by dividing by margin total, and the metaMDS function from the vegan package [93]) was used. Correlations between the ordination axes and selected environmental and management factors were calculated using the envfit function with 999 permutations. Additionally, the impact of the factors "year of the management cycle" and "organic detritus" for vascular plants and "year of sampling" for diatoms was tested by PER-MANOVA [94] with 999 permutations. The ordinations of both datasets were plotted, with these factors passively projected to illustrate the associations with gradients in species compositions. For diatom assemblages, selected environmental factors were also plotted as a smooth surface by using the ordisurf function. All the analyses were conducted in the R program [95].
Constrained ordination (RDA) was used to evaluate the affinity of diatoms to macrophytes. Here, individual macrophyte species and life history trait-and leaf morphologybased functional groups of macrophytes were used as factors to explain patterns in the ordination, with the year of sampling used as a covariable to filter out interannual variability. To evaluate the affinity of specific diatom functional groups, ordinations were also done for each of the FT2 and FT5 species subsets. p values were calculated based on 999 permutations.

Environmental Parameters
Our data confirmed a high trophy of the aquatic environment in the Dehtář fishpond during the study period (Table S2). The values of total phosphorus (TP), chlorophyll-a concentrations, and water transparency (WT) correspond to hypertrophy (the thresholds between hypertrophy and eutrophy according to the OECD 1982 classification [65] are: TP 0.75 mg/L, chl-a 100 µg/L, and WT 0.5 m; the thresholds adjusted for fishpond conditions [66] are: TP 0.2 mg/L, chl-a 30 µg/L, and WT 0.5 m). Associated criteria such as the values of total nitrogen (TN) are also consistent with this classification. There were only small differences in the analysed parameters between the three sites. BSB showed higher average concentrations of several nutrients (e.g., TP and total Fe) and usually lower water transparency than DSB and FD; however, the highest measured chlorophyll-a concentrations were detected at FD, probably due to the cyanobacterial mats concentrated at this site because of wave action. Most of the trophic parameters show a gradual increase during the study period (Table S2).
The photosynthetically active radiation (PhAR), oxygen concentration, and water temperature measured continually showed large fluctuations (besides the diurnal ones, which are typical of these parameters) in some periods ( Figures S16-S20). The course of oxygen concentration followed the same trends at all three sites; however, the maxima were the lowest at FD and the highest at DSB ( Figure S18). Low oxygen concentrations were also common, particularly at BSB and FD; however, at the beginning of August 2016, serious oxygen depletion was detected at all three sites ( Figure S18). Shortly before and during the diatom sampling in June 2015 and 2016, the oxygen concentrations fell into the range of (3.5)5-10(15) mg/L, with a marked difference between FD and BSB in 2016 ( Figure S19d,e). Water temperature was rather high during the summers, usually ranging around (15)20-25(30) • C ( Figure S20). Although the temperature courses followed the same trends at all three sites, some differences were identified. BSB exhibited the highest temperature fluctuations and most frequent maxima at around 30 • C or more. On the other hand, the DSB temperature curve showed the smallest fluctuations and a particularly low occurrence of high temperatures. There also was a difference between individual years, with 2015 and 2019 having the warmest and 2016 the coldest water in summer ( Figure S20).
The data from the meteorological sensors point out high precipitation amounts in the summer of 2014, extremely low precipitation during the whole growing season in 2015 and in the spring of 2016, and a subsequent precipitation extreme in July 2016 (Table A1). Similarly, the mean monthly temperature and temperature course during the growing seasons exhibit extremes, particularly in the summer of 2015, when day temperatures often reached values around 35 • C ( Figure S21). In 2014 and 2016, summer temperatures only rarely reached values 30 • C or more. In an open fishpond space, fully exposed to the sun, air temperature was, however, probably much higher in all years (data not available). The large day/night summer temperature fluctuations are remarkable, whereas the night minima in some periods fell below 10 • C ( Figure S21). The data on wind speed, available for the whole study period, show high frequency of speed above 5.5 m/s (=the winds raising small waves) in 2014 and 2015 and its lower frequency in 2016, but the year 2015 had a higher proportion of situations without any winds. Nevertheless, all these years may be considered windy, with the year 2015 also having a higher number of strong winds (categories > 8.0 and > 13.9 m/s) with the potential to raise large waves (Table 1, Figure 2, Figure S22). In contrast, the year 2019 had a high proportion of situations without any wind, and winds with speed higher than 5.5 m/s were rather rare.  Our evaluation of factors with a potentially important impact on vascular plants, particularly the emergent species, showed that the extent of exposed bottom zone was substantially higher in 2015 than in the other years, particularly in the summer of 2016.
This fact is directly reflected in the size of individual mapping segments (=vegetated zone; Table S1), which was several times larger in 2015 than in the other years, particularly in 2014 and 2016. These fluctuations are related to changes in the basic parameters of the first and the second year of the fishpond management cycle, particularly to the water depth and extent of the flooded area (Table S1). However, there were also differences between particular years of the same phase of the management cycle, e.g., 2015 and 2019 (both first years of the management cycle). In 2015, exposed muddy substrata including the relatively deep muds, important for wetland annuals, were more often available. However, the Dehtář fishpond is principally poor in very deep muddy sediments in the vegetated zone, and most of the muddy zones do not reach more than 10 cm in depth. Organic detritus regularly occurred in the bays of both streams, but in 2016 it was also frequent elsewhere (remnants of vegetation from 2015). Some segments were under the direct influence of manuring, liming, or supplemental fish feeding with cereals (Table S1).   In total, 28 threatened species have been identified in the Dehtář fishpond habitats during the whole study period. The categories C3 (vulnerable) and C4 (near threatened) were the most numerous, but also the C1 (critically threatened) and C2 (endangered) categories were represented by a few species (Figure 4). Wetland annuals (e.g., Carex bohemica, Coleanthus subtilis, Lindernia procumbens, and Lythrum hyssopifolia) were the most numerous group of the threatened species, followed by wetland perennials (e.g., Bolboschoenus laticarpus, B. yagara, and Leersia oryzoides). Among the wetland annuals, Montia arvensis (Table S4, Figure S9) is one of the rarest plant species of the Czech Republic. Red-listed aquatics (Chara braunii) and amphibious species (Elatine hydropiper) were represented as well. Some terrestrial plants also occurred among the threatened taxa, but they mostly appeared very rarely (e.g., Aphanes arvensis; Table S4).

Vegetation Structure and Species Composition of Vascular Plant and Charophyte Assemblages
Plant assemblages of the Dehtář fishpond follow a moisture gradient of the freshwater littoral. Three basic zones may be distinguished: 1) Zone of permanent reed and tall sedge beds and of willow shrubs, 2) zone of disturbed grasslands (usually wet), temporary reed beds and wetland annuals, 3) zone of amphibious plants and submerged aquatics. The first zone is placed on the transition between the high fishpond bank and the upper littoral. It is rather narrow, usually formed by Phalaris arundinacea, Carex acuta, and Salix cinerea, more rarely by Phragmites australis or Salix triandra (Figures S2, S3, S5 and S6). Depending on the littoral morphology and substrate distribution, the second zone may either include several sub-zones ( Figures S4 and S5), predominated usually by Trifolium spp. (mainly T. hybridum, on dry sites T. arvense) and Alopecurus geniculatus (grasslands), Bolboschoenus maritimus agg. (mainly B. laticarpus; temporary reed beds) and a mixture of tall and low-growing wetland annuals (e.g., Bidens radiatus, Rumex maritimus, and Veronica anagallis-aquatica in muddy places and Juncus bufonius, Myosurus minimus, and Stellaria alsine at sandy sites), or it is formed by a mosaic of above plant assemblages. The second zone occupies many hectares of the fishpond bottom during low water levels in the

Vegetation Structure and Species Composition of Vascular Plant and Charophyte Assemblages
Plant assemblages of the Dehtář fishpond follow a moisture gradient of the freshwater littoral. Three basic zones may be distinguished: 1) Zone of permanent reed and tall sedge beds and of willow shrubs, 2) zone of disturbed grasslands (usually wet), temporary reed beds and wetland annuals, 3) zone of amphibious plants and submerged aquatics. The first zone is placed on the transition between the high fishpond bank and the upper littoral. It is rather narrow, usually formed by Phalaris arundinacea, Carex acuta, and Salix cinerea, more rarely by Phragmites australis or Salix triandra (Figures S2, S3, S5 and S6). Depending on the littoral morphology and substrate distribution, the second zone may either include several sub-zones ( Figures S4 and S5), predominated usually by Trifolium spp. (mainly T. hybridum, on dry sites T. arvense) and Alopecurus geniculatus (grasslands), Bolboschoenus maritimus agg. (mainly B. laticarpus; temporary reed beds) and a mixture of tall and low-growing wetland annuals (e.g., Bidens radiatus, Rumex maritimus, and Veronica anagallis-aquatica in muddy places and Juncus bufonius, Myosurus minimus, and Stellaria alsine at sandy sites), or it is formed by a mosaic of above plant assemblages. The second zone occupies many hectares of the fishpond bottom during low water levels in the growing season (Figures S4, S11 and S12). The third zone usually includes only a narrow strip of vegetated shallow water, predominated by an amphibious species Elatine hydropiper ( Figure S8) and an aquatic species Zannichellia palustris.
While the extent of the first zone was more or less unchanged during the study period, particularly the second but also the third zone varied greatly in size between the years or even in different periods of a single year (Figures S2, S3, S11-S13). For instance, the second zone was best developed during the extremely dry year 2015 which corresponds to the large extent of bottom exposure (Table S1, Figures S4 and S12). On the other hand, the third zone disappeared at the end of the summer of 2015, as it was drained. In contrast, in 2016 the second zone only occurred under a low water level in early spring (low water level continued from 2015) and later disappeared due to extremely high precipitation (Table A1) and increased water levels. The third zone with aquatics, which has a summer phenology, did not appear at all in 2016.
The species which occurred as dominants in the vegetation usually also showed very high frequency within the 30 segments. For instance, the species Phalaris arundinacea with the frequency 100%, followed by, e.g., Alopecurus aequalis, Bidens radiatus, Carex acuta, Juncus bufonius, Myosurus minimus ( Figure S10), Trifolium hybridum, Rumex maritimus, and Veronica anagallis-aquatica with the frequency 90-100% in the mapped segments, were also the most common dominants (Table S4). Still, some of the species were an exception to this rule: Phragmites australis occurred in only 40% of segments, nearly always forming stands, although they were not very large. Among the species occurring with high frequency (90-100%) but not forming stands were, e.g., common wetland species Lythrum salicaria, Oenanthe aquatica, and Persicaria hydropiper, or Tripleurospemum inodorum as a ruderal. Only 50 plant species appeared as dominants in at least one segment. Within the segments the frequency of most of them was at least 50%, although there were some exceptions, e.g., terrestrial species dominating only locally at the dam, such as Geranium pusillum and Lamium purpureum.
The frequency of many species fluctuated highly between the years (Table S4). While in 2019 many species were not recorded due to lower research intensity, the data from 2014-2016 reflect the presence or absence of some of the vegetation zones in which the particular species usually grow.

Ecological and Functional Groups of Vascular Plants and Charophytes
Only about one-third of the 281 vascular plant and charophyte taxa were classified within some of the groups of wetland and aquatic species (Figure 5a,b). The rest of the species pool was represented by terrestrial plants with an optimum in various types of habitats, mainly in grasslands and in ruderal and arable weed vegetation (Figure 5a,b). cipitation (Table A1) and increased water levels. The third zone with aquatics, which has a summer phenology, did not appear at all in 2016.
The species which occurred as dominants in the vegetation usually also showed very high frequency within the 30 segments. For instance, the species Phalaris arundinacea with the frequency 100%, followed by, e.g., Alopecurus aequalis, Bidens radiatus, Carex acuta, Juncus bufonius, Myosurus minimus ( Figure S10), Trifolium hybridum, Rumex maritimus, and Veronica anagallis-aquatica with the frequency 90-100% in the mapped segments, were also the most common dominants (Table S4). Still, some of the species were an exception to this rule: Phragmites australis occurred in only 40% of segments, nearly always forming stands, although they were not very large. Among the species occurring with high frequency (90-100%) but not forming stands were, e.g., common wetland species Lythrum salicaria, Oenanthe aquatica, and Persicaria hydropiper, or Tripleurospemum inodorum as a ruderal. Only 50 plant species appeared as dominants in at least one segment. Within the segments the frequency of most of them was at least 50%, although there were some exceptions, e.g., terrestrial species dominating only locally at the dam, such as Geranium pusillum and Lamium purpureum.
The frequency of many species fluctuated highly between the years (Table S4). While in 2019 many species were not recorded due to lower research intensity, the data from 2014-2016 reflect the presence or absence of some of the vegetation zones in which the particular species usually grow.

Ecological and Functional Groups of Vascular Plants and Charophytes
Only about one-third of the 281 vascular plant and charophyte taxa were classified within some of the groups of wetland and aquatic species (Figure 5a,b). The rest of the species pool was represented by terrestrial plants with an optimum in various types of habitats, mainly in grasslands and in ruderal and arable weed vegetation (Figure 5a Table S4. Among the wetland plants, annual and shortly perennial species of exposed pond bottoms formed the most numerous group, followed by reed bed species (i.e., helophytes, Figure 5a,b). On the other hand, the number of aquatic species (floating-leaved, freefloating, and submerged plants) was very low. The representation of particular groups varied between the years. This holds true particularly for small wetland annuals, amphibious, and submerged aquatic plants, which were most frequent in 2014 and 2015, some of them also in 2019 (the years with exposed bottom and shallow water available throughout the growing season, see Tables S1 and S4) but rare or even missing in 2016, i.e., the year with high water levels in summer (Table S4, Figure 5a,b). The wetland perennials growing in lower littoral parts (e.g., Bolboschoenus maritimus agg.), ruderal and arable weeds, and trees and shrubs (species occurring only as seedlings) also exhibit a strong decline in species numbers in 2016 (Table S4). On the other hand, grassland species, e.g., Holcus lanatus, Lychnis flos-cuculi, and Potentilla anserina profited from a long dry phase in the spring of The year 2019 shows apparently lower proportion of terrestrial species due to the fact that the research focused on the shallow water and exposed bottom zones (Figure 5a,b).

Main Drivers of Species Composition in Vascular Plant and Charophyte Assemblages
Most of the factors tested have been significantly correlated with the NMDS ordination, where the first ordination axis corresponds to the gradient of flood and exposure duration, while the second axis exhibits the relation to the depth of muddy sediment and the distance from large inflows ( Figure 6). The extent of summer exposure and the distance from large inflows have been identified as the factors with the highest impact on the species composition of vascular plants and charophytes in the plant assemblages of the Dehtář fishpond. A less important relationship was detected for mud depth and, with a weak significance, for a number of small inflows. Research intensity has a negative relationship with the second ordination axis, but the significance (p = 0.067) was below the commonly used threshold. The two remaining factors, i.e., the surrounding landscape with prevailing arable land and shrub share in the buffer zone did not show any relationship to the species composition of the plant assemblages. also in 2019 (the years with exposed bottom and shallow water available throughout the growing season, see Tables S1 and S4) but rare or even missing in 2016, i.e., the year with high water levels in summer (Table S4, Figure 5a,b). The wetland perennials growing in lower littoral parts (e.g., Bolboschoenus maritimus agg.), ruderal and arable weeds, and trees and shrubs (species occurring only as seedlings) also exhibit a strong decline in species numbers in 2016 (Table S4). On the other hand, grassland species, e.g., Holcus lanatus, Lychnis flos-cuculi, and Potentilla anserina profited from a long dry phase in the spring of 2016 and partly in 2014. The year 2019 shows apparently lower proportion of terrestrial species due to the fact that the research focused on the shallow water and exposed bottom zones (Figure 5a,b).

Main Drivers of Species Composition in Vascular Plant and Charophyte Assemblages
Most of the factors tested have been significantly correlated with the NMDS ordination, where the first ordination axis corresponds to the gradient of flood and exposure duration, while the second axis exhibits the relation to the depth of muddy sediment and the distance from large inflows ( Figure 6). The extent of summer exposure and the distance from large inflows have been identified as the factors with the highest impact on the species composition of vascular plants and charophytes in the plant assemblages of the Dehtář fishpond. A less important relationship was detected for mud depth and, with a weak significance, for a number of small inflows. Research intensity has a negative relationship with the second ordination axis, but the significance (p=0.067) was below the commonly used threshold. The two remaining factors, i.e., the surrounding landscape with prevailing arable land and shrub share in the buffer zone did not show any relationship to the species composition of the plant assemblages. . NMDS ordination diagram of species composition of vascular plant and charophyte assemblages of the Dehtář fishpond with passively projected selected environmental and management factors. Significant relationship of the factors to species composition is indicated based on 999 permutations. ***-highly significant (p ≤ 0.001), *-weakly significant (p ≤ 0.05). The correlations were calculated using the envfit function from the vegan package (Oksanen [93]). Explanations of abbreviated names of the selected factors: summer exposure (=extent of bottom exposure in summer), research intensity (=number of visits in the given segment and year), shrub share (% share of shrubs in the buffer zone between open water and the surrounding landscape), arable (=surrounding landscape with prevailing arable land). For more detail, see Table S1.
In the upper left part of the NMDS ordination diagram ( Figure 6) the species of exposed substrates, particularly wetland annuals (e.g., Cyperus fuscus, Isolepis setacea, Persicaria lapathifolia) and amphibious species (e.g., Eleocharis acicularis) are concentrated, with some of them showing relationship to a higher depth of mud (e.g., Veronica anagallis-aquatica and Peplis portula). On the opposite side of the gradient, mainly grassland species (e.g., Ranunculus repens, Plantago lanceolata) and other terrestrial plants (e.g., Galium aparine) are displayed. The gradient of the distance from large inflows shows species composition of plant assemblages near the fishpond dam that occurred on dry sandy substrates. These sites were occupied mainly by weeds of arable land (e.g., Apera spica-venti, Thlaspi arvense), ruderal weeds (e.g., Descurainia sophia, Artemisia vulgaris), and annual ephemeral herbs of dry disturbed habitats (e.g., Arabidopsis thaliana, Valerianella locusta). All these terrestrial species have been recorded, particularly in the years with high research intensity. The opposite part of this gradient shows the vicinity of large inflows, with Scirpus radicans growing directly in one of them.
The "spider" diagram of the NMDS ordination showing the influence of the year of the management cycle indicates that the samples (species lists collected separately for each segment and year) are well separated between year 1 and 2 (Figure 7a), and the difference in their species composition is significant (PERMANOVA, F = 8.1855, p = 0.005).
The "spider" diagram of the NMDS ordination showing the influence of the year of the management cycle indicates that the samples (species lists collected separately for each segment and year) are well separated between year 1 and 2 (Figure 7a), and the difference in their species composition is significant (PERMANOVA, F = 8.1855, p = 0.005).
The figure also clearly demonstrates the predominance of annual species in the species spectrum of year 1, and, by contrast, of perennial species in the samples of year 2, including wetland species of the upper littoral zone. Figure 7b shows the "spider" NMDS ordination diagram with the relationship between organic detritus and species composition. Visually, there are rather large overlaps in the species composition of the segments with and without organic detritus. Nevertheless, the difference in species composition between the contrast sites was significant (PERMANOVA, F=6.2197, p = 0.012). Wetland annuals such as Cyperus fuscus, Isolepis setacea, and seedlings of Typha sp. are clearly related to the segments without substantial amounts of organic detritus. On the other hand, nitrophytes such as Galium aparine and Urtica dioica show the preference for segments with organic detritus deposits.  Table S1.

Diatom Species and Functional Group Richness
Diatom species richness varied greatly among the samples, particularly in 2015 (min. 14, max. 46 species) and 2016 (min. 13, max. 48 species); in 2019, the differences were smaller (min. 23, max. 35) (Figure 8a-c). The same applied even to samples from the same site (e.g., FD 2015 or BSB 2016) and, in some cases, for the same macrophyte species. Nevertheless, the tendency to support higher or lower diatom species richness was visible in several macrophyte species. Macrophytes with high diatom species richness, reaching the  Table S1.
The figure also clearly demonstrates the predominance of annual species in the species spectrum of year 1, and, by contrast, of perennial species in the samples of year 2, including wetland species of the upper littoral zone. Figure 7b shows the "spider" NMDS ordination diagram with the relationship between organic detritus and species composition. Visually, there are rather large overlaps in the species composition of the segments with and without organic detritus. Nevertheless, the difference in species composition between the contrast sites was significant (PERMANOVA, F = 6.2197, p = 0.012). Wetland annuals such as Cyperus fuscus, Isolepis setacea, and seedlings of Typha sp. are clearly related to the segments without substantial amounts of organic detritus. On the other hand, nitrophytes such as Galium aparine and Urtica dioica show the preference for segments with organic detritus deposits.

Diatom Species and Functional Group Richness
Diatom species richness varied greatly among the samples, particularly in 2015 (min. 14, max. 46 species) and 2016 (min. 13, max. 48 species); in 2019, the differences were smaller (min. 23, max. 35) (Figure 8a-c). The same applied even to samples from the same site (e.g., FD 2015 or BSB 2016) and, in some cases, for the same macrophyte species. Nevertheless, the tendency to support higher or lower diatom species richness was visible in several macrophyte species. Macrophytes with high diatom species richness, reaching the maxima above 40 and the minima higher than 25 diatom species per sample, were, e.g., Elatine hydropiper (n = 5, min. 28  All five functional groups of diatoms (i.e., FT1-FT5) were represented in most of the samples (Figure 8a-c). The proportion of particular groups in the species composition was rather stable across the samples, with low species numbers at FT1 and FT4, middle high at FT2 and FT3, and the highest at FT5 (Figure 8a-c). However, FT5 exhibited the largest fluctuations in species numbers among the samples. The extraordinary species-rich samples had all five functional groups, with especially high numbers at FT5, and in 2016 also FT2 and FT3 (Figure 8a Table S3. Host plant species is listed at each of the samples. All five functional groups of diatoms (i.e., FT1-FT5) were represented in most of the samples (Figure 8a-c). The proportion of particular groups in the species composition was rather stable across the samples, with low species numbers at FT1 and FT4, middle high at FT2 and FT3, and the highest at FT5 (Figure 8a-c). However, FT5 exhibited the largest fluctuations in species numbers among the samples. The extraordinary species-rich samples had all five functional groups, with especially high numbers at FT5, and in 2016 also FT2 and FT3 (Figure 8a

Quantitative Representation of Diatom Species and Their Functional Groups
Nitzschia palea (FT4 functional group) was the most abundant diatom species in our samples in terms of its frequency and quantity. It occurred in 56 of the total of 58 samples (frequency 96.6%) with 2,986 individuals altogether (Table S5). Its highest concentrations were detected in 2016 when it was found in all the samples and strongly dominated in some of them (Table S6). Puncticulata balatonis, a planktonic diatom (FT1), occurred in 54 samples altogether (frequency 93%) but in low quantities (low numbers of individuals of FT1 species are typical for epiphytic diatom samples). Nine additional taxa from various  All five functional groups of diatoms (i.e., FT1-FT5) were represented in most of the samples (Figure 8a-c). The proportion of particular groups in the species composition was rather stable across the samples, with low species numbers at FT1 and FT4, middle high at FT2 and FT3, and the highest at FT5 (Figure 8a

Quantitative Representation of Diatom Species and Their Functional Groups
Nitzschia palea (FT4 functional group) was the most abundant diatom species in our samples in terms of its frequency and quantity. It occurred in 56 of the total of 58 samples (frequency 96.6%) with 2,986 individuals altogether (Table S5). Its highest concentrations were detected in 2016 when it was found in all the samples and strongly dominated in some of them (Table S6). Puncticulata balatonis, a planktonic diatom (FT1), occurred in 54 samples altogether (frequency 93%) but in low quantities (low numbers of individuals of FT1 species are typical for epiphytic diatom samples). Nine additional taxa from various  Table S3 for more detail).

Quantitative Representation of Diatom Species and Their Functional Groups
Nitzschia palea (FT4 functional group) was the most abundant diatom species in our samples in terms of its frequency and quantity. It occurred in 56 of the total of 58 samples (frequency 96.6%) with 2,986 individuals altogether (Table S5). Its highest concentrations were detected in 2016 when it was found in all the samples and strongly dominated in some of them (Table S6). Puncticulata balatonis, a planktonic diatom (FT1), occurred in 54 samples altogether (frequency 93%) but in low quantities (low numbers of individuals of FT1 species are typical for epiphytic diatom samples). Nine additional taxa from various functional groups had a frequency of between 75% and 90% within the whole dataset (=highly frequent taxa). Among them, e.g., Navicula capitatoradiata and Nitzschia archibaldii also occurred in high quantities, strongly predominating the diatom assemblages of some of the samples. Most of the highly frequent and frequent species (occurring with the frequency between 50% and 75% within the dataset) have shown a strong affinity to a particular year and/or site (Table S6; see also Section 3.3.3). Pseudostaurosira elliptica, P. brevistriata, and Staurosirella leptostauron var. dubia in 2015; Achnanthidium saprophilum (probably on ESEM photo- Figure 14c), Cocconeis placentula, and Navicula reichardtiana in 2016, and Gomphonema parvulum (Figure 14d) in 2016 and 2019 occurred with a high frequency and, at the same time, in a high number of individuals. The frequent species also included those found in a small number of individuals only, e.g., Staurosira binodis, Navicula gregaria, and Placoneis clementiodes, concentrated into the year 2015, or Craticula subminuscula and Ulnaria acus, typical for 2016. The species that had, despite low frequencies, relatively high numbers of individuals were exceptional: Lemnicola hungarica and Eunotia bilunaris (both FT2) occurred only in eight and ten samples, respectively, but especially the former sometimes produced very dense populations (Tables S5 and S6).
The proportion of the functional groups of diatoms strongly varied between the years and sites. In 2015, the groups FT3 and FT5 predominated in the samples, with the former being more frequent at BSB and DSB and the latter in the samples from FD (Figure 11a). In 2016, only FD and BSB were sampled. While typically periphytic diatoms (FT2 group) predominated the samples from FD, the FT4 taxa occurred with the highest quantities in most samples from BSB (Figure 11b). In 2019, when two microsites at FD were sampled, the share of the functional groups among the samples was rather balanced. While the samples from the more muddy site were predominated by the FT2 and FT5 groups (Figure 11c, samples F1-19 to F4-19), the highest share was achieved by FT2 and FT3 in the samples from the sandy site exposed to wave action (Figure 11c, FN5-19 to FN8-19).
The predominance of a single or a few diatom taxa is strongly reflected in the share of their functional species groups within the samples (Figure 11a

Impact of Environmental Factors on Structure and Species Composition of Diatom Assemblages
All the tested factors have been significantly correlated with the NMDS ordination. The first ordination axis corresponds to the gradients of flood and exposure duration, mud depth, and, in the opposite direction, chlorophyll-a concentrations during the growing season. The second axis corresponds to the gradient of the distance from large inflows where the largest distance corresponds to the proximity of the fishpond dam (FD). The other factors correlated with this gradient, i.e., shrub share in the buffer zone, surrounding landscape formed by arable fields, and number of small inflows, are clearly associated with the FD sampling site ( Figure 12).
Water 2021, 13, x FOR PEER REVIEW 19 o functional groups had a frequency of between 75% and 90% within the whole dat (=highly frequent taxa). Among them, e.g., Navicula capitatoradiata and Nitzschia archib also occurred in high quantities, strongly predominating the diatom assemblages of so of the samples. Most of the highly frequent and frequent species (occurring with the quency between 50% and 75% within the dataset) have shown a strong affinity to a ticular year and/or site (Table S6; (Table S5, S6). The proportion of the functional groups of diatoms strongly varied between the y and sites. In 2015, the groups FT3 and FT5 predominated in the samples, with the for being more frequent at BSB and DSB and the latter in the samples from FD (Figure 1 In 2016, only FD and BSB were sampled. While typically periphytic diatoms (FT2 gro predominated the samples from FD, the FT4 taxa occurred with the highest quantitie most samples from BSB (Figure 11b). In 2019, when two microsites at FD were samp the share of the functional groups among the samples was rather balanced. While samples from the more muddy site were predominated by the FT2 and FT5 groups (Fig  11c, samples F1-19 to F4-19), the highest share was achieved by FT2 and FT3 in the s ples from the sandy site exposed to wave action (Figure 11c, FN5-19 to FN8-19).
The predominance of a single or a few diatom taxa is strongly reflected in the sh of their functional species groups within the samples (Figure 11a

Impact of Environmental Factors on Structure and Species Composition of Diatom Assemblages
All the tested factors have been significantly correlated with the NMDS ordinatio The first ordination axis corresponds to the gradients of flood and exposure duratio mud depth, and, in the opposite direction, chlorophyll-a concentrations during the grow ing season. The second axis corresponds to the gradient of the distance from large inflow where the largest distance corresponds to the proximity of the fishpond dam (FD). Th other factors correlated with this gradient, i.e., shrub share in the buffer zone, surroundin landscape formed by arable fields, and number of small inflows, are clearly associate with the FD sampling site (Figure 12).
The species composition was significantly different between the years of samplin (here a calendar year) (PERMANOVA, F = 34.638, p = 0.001). The differences are als clearly demonstrated in the NMDS ordination diagram (Figure 12). For the year 2015, r markable for its large area of exposed pond bottom and high depth of mud (vegetatio zones 2 and 3, see part 3.2.2), FT3 diatoms such as Staurosirella leptostauron and Pseudosta rosira elliptica were the most typical, followed by some Navicula species (FT5 group), e.g N. cryptocephala and N. trivialis ( Figure 12, Table S6).  Table S1.
In contrast, the year 2016 was characterized by high water levels, and therefore only the vegetation zones 1 and 2 (deeply flooded) could be sampled. FT4 diatoms of the genus Nitzschia, particularly N. archibaldii and N. palea, as well as the FT5 species such as Pinnularia marchia and Ulnaria ulna differentiated the sites distant from the fishpond dam at BSB from the sites distant from the large inflows at FD. Samples from FD 2016 were predominated by periphytic diatoms (FT2), e.g., Achnanthidium minutissimum, A. saprophilum, Coc-  Table S1. The species composition was significantly different between the years of sampling (here a calendar year) (PERMANOVA, F = 34.638, p = 0.001). The differences are also clearly demonstrated in the NMDS ordination diagram (Figure 12). For the year 2015, remarkable for its large area of exposed pond bottom and high depth of mud (vegetation zones 2 and 3, see Section 3.2.2), FT3 diatoms such as Staurosirella leptostauron and Pseudostaurosira elliptica were the most typical, followed by some Navicula species (FT5 group), e.g., N. cryptocephala and N. trivialis (Figure 12, Table S6).
In contrast, the year 2016 was characterized by high water levels, and therefore only the vegetation zones 1 and 2 (deeply flooded) could be sampled. FT4 diatoms of the genus Nitzschia, particularly N. archibaldii and N. palea, as well as the FT5 species such as Pinnularia marchia and Ulnaria ulna differentiated the sites distant from the fishpond dam at BSB from the sites distant from the large inflows at FD. Samples from FD 2016 were predominated by periphytic diatoms (FT2), e.g., Achnanthidium minutissimum, A. saprophilum, Cocconeis placentula, and Lemnicola hungarica. Increased chlorophyll-a concentrations probably modulated the species composition at FD in 2016. The FT2 species Gomphonema parvulum and FT5 species Navicula capitatoradiata were present in the samples from each year and site, but they were particularly frequent at FD in 2015 (higher share of N. capitatoradiata) and 2019 (higher share of G. parvulum).
The relationship of diatom assemblages to three important factors, not considered above, is shown in Figure 13: distance from the fishpond bank (a), water depth (b), and exposure to disturbances caused by wave action or streaming (c). These NMDS ordination diagrams show not only the species distribution patterns between the three sites and sampling years, but also the differentiation or overlaps between the samples.
The samples from the three sites collected in 2015 exhibit a high level of similarity in species composition, all being situated a long distance from the bank and in very shallow, disturbed parts of the lower littoral (see Table S3). The two sites sampled in 2016, FD and BSB, exhibit a high level of differentiation between the samples: while the disturbance intensity and distance from the bank were similar for both sites, water depth was markedly higher at FD. At the same time, the samples from both sampling sites in 2016 were clearly differentiated from those from 2015 by a substantially higher water depth, and shorter distance from the bank and disturbance intensity.  Tables S3 and S6).
On a small scale, a particular spatial position of the sampled macrophyte individuals and the distances between them probably also play a role. For instance, diatom assemblages of Glyceria maxima, Persicaria hydropiper, Veronica anagallis-aquatica, and Ranunculus sceleratus (samples B17-16 to B20-16) growing at BSB in the distance of up to about 1 m from one another, were highly similar in their species composition, species richness, and proportions of functional groups of diatoms within the species pools (Figure 8b, Table S6). The quantitative share of particular diatom species and functional groups varied, however, even at these spatially close samples. Similarly, the samples of Phragmites australis and Phalaris arundinacea (B12-16 and B13-16) collected at the same site in a sandy littoral at BSB (segment 6-see Figure 1 and Table S3) strongly differed in the share of functional groups of diatoms (Figure 11b). These differences were also well visible on the microphotographs of plant surfaces with diatom assemblages acquired using ESEM (Figure 14a,b). The overall density of diatom assemblages also strongly varied, as shown by a set of ESEM photographs from the sites with contrasting conditions (e.g., exposure to wave-mediated disturbances) or the different macrophytes from one and the same site and year (Figure 14a-e).  Table S3. 14a-e).

Relationships between Macrophytes and Epiphytic Diatom Assemblages
We did not find any significant relationship between diatom species composition and particular macrophyte species either in the dataset of all macrophyte samples and all diatom functional groups (Permutation test for RDA; p = 0.72) or in the test of all macrophyte samples, and FT2 (p = 0.565) or FT5 (p = 0.825) diatom groups, respectively.

Relationships between Macrophytes and Epiphytic Diatom Assemblages
We did not find any significant relationship between diatom species composition and particular macrophyte species either in the dataset of all macrophyte samples and all diatom functional groups (Permutation test for RDA; p = 0.72) or in the test of all macrophyte samples, and FT2 (p = 0.565) or FT5 (p = 0.825) diatom groups, respectively.
When we used the ecology and life history traits-based functional groups of macrophytes instead of vascular plant and charophyte species, we detected a highly significant response of all diatom species (p = 0.001) as well as of diatom species selected only from the functional groups FT2 (p = 0.002) and FT5 (p = 0.001). The use of leaf morphology-based functional groups of macrophytes did not reveal a significant relationship to species composition within the two functional groups of diatoms (p = 0.116 for FT2, p = 0.127 for FT5). The relationship between all diatom species and morphology-based macrophyte groups was only weakly significant (p = 0.048). The effect of the macrophyte species groups on diatom species richness is shown in Figure 15. In the graph of diatom species richness in the groups of samples sorted according to ecology and trait-based classification (Figure 15a), a trend of higher species richness in amphibious and submerged plants is visible. On the other hand, floating-leaved species exhibit low species richness; however, they were present only in 2016. When the morphological classification of macrophytes was used, no clear relationship between the species richness and a particular macrophyte group was found (Figure 15b). The relationship between all diatom species and morphology-based macrophyte groups was only weakly significant (p = 0.048). The effect of the macrophyte species groups on diatom species richness is shown in Figure 15. In the graph of diatom species richness in the groups of samples sorted according to ecology and trait-based classification (Figure 15a), a trend of higher species richness in amphibious and submerged plants is visible. On the other hand, floating-leaved species exhibit low species richness; however, they were present only in 2016. When the morphological classification of macrophytes was used, no clear relationship between the species richness and a particular macrophyte group was found (Figure 15b).  Table S3.

Joint Characteristics of Plant and Epiphytic Diatom Assemblages
Vascular plants (including a single charophyte species) and epiphytic diatoms in the Dehtář fishpond exhibit similarities in several general trends. These are, in particular: (1) the species composition, as well as the proportion of functional and ecological groups, fluctuates significantly between individual years and sites. (2) A substantially higher spe-  Table S3.

Joint Characteristics of Plant and Epiphytic Diatom Assemblages
Vascular plants (including a single charophyte species) and epiphytic diatoms in the Dehtář fishpond exhibit similarities in several general trends. These are, in particular: (1) the species composition, as well as the proportion of functional and ecological groups, fluctuates significantly between individual years and sites. (2) A substantially higher species richness in 2015 (first year of the fishpond management cycle with extreme summer drought and the largest extent of fishpond bottom exposure) as compared to 2016 (second year of the management cycle with a very dry spring but a very wet summer, causing extremely high water levels): We suppose that these differences are mainly related to the differences in the extent of fishpond bottom exposure, which is typically much larger in the first than in the second year of the management cycle; however, irregularities may also occur due to climatic extremes. Our incomplete data from 2014 (second year of the management cycle but with water levels lower than typical for this period) and 2019 (first year of the management cycle) at least partly support this presumption. For macrophytes, similar results have been obtained also by Francová et al. ([19,20]). (3) The proportion of taxa with high tolerance or preference of water level fluctuations and/or drought tolerance: In particular, there is a high share of non-wetland species among vascular plants, a high frequency of wetland annual species preferring sandy substrata and the occurrence of aerotolerant diatoms, and, in contrast, there is the absence of species preferring a stable water column. The overall species composition and share of functional and ecological species groups of vascular plants and diatoms correspond to a temporary water body (e.g., [18,[96][97][98][99]. (4) The occurrence of taxa with various nutrient demands at one site or even within a single stand/sample: Particularly in the first year of the management cycle, the joint occurrence of species preferring oligotrophic to mesotrophic habitats with species typical of highly eutrophic habitats was rather common, although the former usually did not exhibit high abundances in case of diatoms. The final trend was (5) an exceptionally high overall species richness in comparison to the studies from similar temperate aquatic habitats (lakes, fishponds, water reservoirs (see [52][53][54]60] for diatoms and [18][19][20]26,36,100] for macrophytes).

Nutritional Gradient
The plant and diatom species and functional diversity of the Dehtář fishpond is very high and complex. We analysed a number of environment-and management-related factors in order to disentangle their influence on the studied groups of autotrophs. Basic nutrients, pH, and conductivity are considered to be important characteristics driving patterns in the distribution of aquatic plants and algae [23,36,49,57,[101][102][103][104][105]. This, however, holds true for the level of separated localities, particularly on a larger spatial scale [3,21,23,50,52], rather than within a single water body. Although Dehtář is a large fishpond and its bays are somewhat separated from the main part, which is connected to the dam, spatial differences in water chemistry are, in our opinion, too small to be interpreted as an important source of variability in autotrophic communities occurring at the different sites. As the system is hypertrophic, we do not suppose any limitations by nutrients. Diatoms could be limited by silicon (Si) availability; however, we have no data regarding this nutrient. Yet, as waters in the study region do not suffer from its lack [106], we suppose it is not a limiting element here. Moreover, silicon is usually positively correlated with iron concentrations [107], which were high in some years. Aquatic plants found in the Dehtář fishpond are able to grow in waters with various trophy levels [108]. For macrophytes rooted in the bottom, the sediment is the main source of nutrients; muddy fishpond sediments are generally nutrient-rich [16].

Light Limitation
Light limitation might be crucial for autotrophs in hypertrophic water bodies [47]. Therefore, we involved monthly averages of chlorophyll-a concentrations from each site into the analyses of epiphytic diatom assemblages; this factor corresponds well with nutrient amounts and is negatively correlated with water transparency. Indeed, we found a clear relationship between chlorophyll-a concentrations and the diatom species composition. This result is in agreement with the findings of Poulíčková et al. [105], where, e.g., Lemnicola hungarica was shown as an important indicator of high trophy and chlorophyll-a concentrations. These characteristics also depict the high nutrient requirements of the host plants of Lemnicola hungarica, in our case Lemna gibba and Spirodela polyrhiza [109].
A generally low light availability is also illustrated by water transparency (Table S2), which often did not exceed 30 cm at our study sites in the summer months, and PhAR, reaching only negligible values in the depth of 1.5 m; this situation is common in fishponds with adult carp [19]. Although in the depth of 0.3 m, PhAR reached rather high values in some periods, these values did not occur continually throughout the daytime. Thus, light could be an important factor limiting species richness and biomass of submerged aquatic and amphibious plants in the Dehtář fishpond, particularly during the second year of the management cycle, when the water level was high. Although in the first year of the management cycle the areas of shallow water, potentially suitable for submerged vegetation, are fairly large in various parts of the fishpond, the instability of daylight in the water column may be of high importance. Nevertheless, most submerged aquatics and amphiphytes, regularly growing in carp ponds, are rather tolerant of low water transparency and are also known to be able to survive periods of long-term low light availability and high fish stock pressure (i.e., the second year of the management cycle) in the sediment seed bank [19,20,30]. Our data show that these species nearly disappeared from the Dehtář fishpond in 2016 (Table S4, Figure 5) but were present again in 2019. Therefore, the very low species richness of aquatic plants (particularly the submerged ones) and amphiphytes in the Dehtář fishpond is surprising. Only several species of the above functional groups occurred in Dehtář, and only two of them, Elatine hydropiper and Zannichellia palustris, were rather frequent, forming species-poor stands in shallow waters along the shoreline. Species common in the fishponds in the study region such as Myriophyllum spicatum and Stuckenia pectinata, which are tolerant of turbid water and high nutrient amounts, did not occur in the Dehtář fishpond at all, and others such as Potamogeton crispus and P. pusillus were extremely rare here [19,110]. We conclude that light availability is responsible for the interannual population fluctuations of submerged aquatics and amphiphytes, but it is not the decisive factor limiting their overall species richness. It is likely that the propagules of many aquatics are dispersed into the Dehtář fishpond from the surrounding water bodies via waterfowl [111], but site conditions do not enable their establishment. We suppose that the reason for the low species richness of aquatic species could be similar to that identified by Jupp and Spence [42]. They described macrophyte limitation by wave action in a lake shallow water zone, while phytoplankton dominance and associated light shortage restricted macrophyte growth in a deeper water zone protected against the impact of waves.

Moisture and Temperature Fluctuations
Another important factor in wetlands is drought. The strong water level decrease in the Dehtář fishpond in 2015 led to the whole vegetated zone finally being out of the water in the late summer; under such conditions, only amphiphytes survived, while the aquatics died. However, the ability of fishpond macrophytes to survive as seeds or spores in the soil seed bank is usually not eliminated by substrate desiccation [30], as mechanisms of ecological filtering in regularly drained fishponds probably selected species with relevant traits decades or even centuries ago. Nevertheless, in very dry years, the fast desiccation of muddy substrates may be an important factor eliminating not only aquatics and amphiphytes, but also some wetland annuals with high moisture demands (e.g., Coleanthus subtilis, Eleocharis ovata, and Peplis portula; [18]). In contrast, droughttolerant wetland annuals such as Cerastium dubium, Montia arvensis, Myosurus minimus, and Spergularia kurkae, colonising large areas of exposed sandy shores, seem to profit from dry years, although a prolonged dry period may increase the competition by perennials, particularly of clonal grassland species. However, the former three wetland annuals with very early development (March and April) usually make use of a narrow time window between the early beginning of the growing season and the progressive development of grasslands in May. S. kurkae is capable of long-term persistence in a soil seed bank [33].
Our analyses revealed that the fishpond summer drainage is one of the most important factors shaping species composition and increasing species richness of wetland and nonwetland vascular plants (see also [19,20]), but also of epiphytic diatoms. We expected that epiphytic diatom communities would be impoverished due to the environmental stress caused by highly unstable moisture during summers and winters (see also [112,113]). However, permanent and temporary wetland habitats usually differ in their diatom species composition, with the taxa dependent on a stable water column or moisture occurring in the former and missing in the latter wetland types [62,114]. Some of the diatom species found in our samples with very low frequency and/or abundance are considered by Evans [114] as species preferring permanent water bodies (e.g., Melosira varians, Nitzschia acicularis, and Gyrosigma acuminatum) and could be introduced from elsewhere by stream water (see also [52]). In this sense, we did not find analogies between diatoms and aquatic plants: the hydrophytes sensitive to even very short-term water drawdown apparently do not occur in fishponds due to regular drainage. Even the aquatics such as some Potamogeton species or Zannichellia palustris are supported by a water level drawdown in some of their development phases [115].
While there is a lot of literature on the ecology of vascular plants and charophytes (e.g., [28][29][30]99,108,[115][116][117]), this is only partly true for diatoms; unresolved taxonomy is one of the reasons why details on the ecology of some diatom species may highly differ between papers [52,105]. An FT2 diatom species Achnanthidium minutissimum, recently regarded as species aggregate, is a typical example: while many authors consider it as an indicator of good water quality and, at the same time, tolerant of stress caused by water level fluctuations (e.g., [62,112,118]), we found its highest density in 2016, when the nutrient concentrations were particularly high and the water level was, after the bottom re-flooding, more stable than in 2015 and 2019. Recent studies prove the existence of several ecologically distinct microspecies [105,119], which makes A. minutissimum agg. an apparently euryvalent species. We propose that the mechanisms that allow the survival of epiphytic diatoms under the conditions of a temporary water unavailability could be one of the crucial attributes of distinct microspecies within the species aggregates. At the same time, these mechanisms should be seen as important determinants shaping the species composition in temporary wetlands [114,120]. Our data suggest that the diatom taxa such as Achnanthidium minutissimum agg., A. saprophillum, and Gomphonema parvulum agg. are able to tolerate high water level fluctuations but prefer a long phase of flooding (see also [98,121]). Their highest densities were identified at FD, and the only sample with abundant G. parvulum agg. from BSB came directly from the Babický brook mouth, which holds water even in the driest years. All these taxa belong to the most frequent dominants of epiphytic diatom assemblages [52,53,105,122,123] where they probably are highly competitive, as we have also shown on ESEM photographs (Figure 14c,d). According to some authors, high water temperatures and nutrient amounts are important for G. parvulum agg. [52,124,125], which corresponds to our water chemistry and temperature data (Table S2, Figure S20j).
At the FD site, diatoms were better indicators of habitat conditions than macrophytes, because sites with quickly increasing water levels are only suitable for a limited number of macrophytic species. Nevertheless, the narrow bay at FD is one of the few sites in the littoral zone of the Dehtář fishpond with deep sapropelic sediments and high substrate moisture, predominated by nutrient-and moisture-demanding wetland annuals and perennials (e.g., Bidens spp., Bolboschoenus laticarpus, Persicaria spp., Oenanthe aquatica; Figure S3a) in late summer. This site is also exceptional for its lack of drought-tolerant small wetland annuals and aerotolerant diatoms. The BSB sampling site had the highest share of species with a more frequent occurrence in temporary than in permanent waters, e.g., Eunotia bilunaris, Nitzschia archibaldii, N. palea (present there in the highest abundances), and Navicula reichardtiana, and species classified as aerotolerant, i.e., highly tolerant of a long phase without flooding, e.g., Hantzschia abundans, Pinnularia borealis, P. marchia, and P. obscura [58,96,114,120,122]. Particularly the species of the latter group occurred at sites close to the bank. Considering vascular plants, a mosaic of moisture demanding and moisture tolerant species occurred at BSB, as the sites with rather deep mud (segments 7, 8, and 9) and sandy and stony fishpond margins (segments 6, 10, and 11) are situated side-by-side there. This part of the fishpond is very shallow and is the last to be overflooded. Therefore, it is an important refuge for small wetland annuals such as Coleanthus subtilis ( Figure S7) even in wet years when suitable sites at FD and DSB are flooded early.
The temperature course reasonably modulated the overall species richness of vascular plants of the Dehtář fishpond during the study period. While diurnal temperature fluctuations are crucial for the germination and early development of a summer annual Coleanthus subtilis, winter annuals such as Cerastium dubium, Montia arvensis, and Myosurus minimus have profited from mild winters in the last decade [116] (see also [86] and Table A1 for climatic data). The frequency of Cyperus fuscus and Lythrum hyssopifolia markedly increased in 2019, the year with extraordinary high summer temperatures but also sufficient precipitation (Table A1 and Table S4). Both of these species are relatively thermophilous and late-germinating [116], but the former has recently been spreading in colder regions such as southern Bohemia [99]. Lindernia procumbens (Table S4) appeared only in the summer of 2015, probably as a consequence of extremely hot and dry weather. This highly thermophilous species otherwise mostly survives in the soil seed bank at many of its sites in Central Europe [110,117].

Omitted Role of Wave Action in Fishponds
Wave action is an important factor in large water bodies including fishponds [126]. Among all the large fishponds in the Czech Republic, there are probably only few in which winds causing wave action occur regularly [9]. The size as well as overall volume, water body shape, and position in the landscape are all important predictors of wave action [43,48]. Indirectly, the waves act through the transport of fine sediment particles into deeper and/or calm parts of the water body, thus forming pure sandy or stony beaches in the upper littoral zone [126]. This effect is particularly important for plant species with specific edaphic requirements [48]. The direct impact of wave action on macrophytes proceeds through the uprooting and fragmenting of submerged species, disturbance of the root system of helophytes, burial of plant parts with sand or organic detritus, and decrease in water transparency [42,127]. Macrophytes able to form creeping morphotypes are at an advantage. In the Dehtář fishpond, Zannichellia palustris, Eleocharis acicularis, and Elatine hydropiper possess the necessary morphological traits and occupy places under a strong wave action [116,128]. Nevertheless, at some sites or in some years, wave-mediated disturbances are probably too strong for any macrophyte species. The substantially lower incidence of strong winds in 2019 (as compared to 2015) probably contributed to the colonisation of the sites that were previously without any submerged vegetation by Zannichellia palustris.
For diatoms, the published information on wave mediated disturbance is scarce, but the FT2 species attached to the substrate with mucilaginous stalks, e.g., Achnanthidium spp. and Gompohonema spp., are reported as more disturbance-sensitive than the adnate FT2 species or the species from other functional groups [98]. This would be in accordance with our results, although it seems to be difficult to separate the influence of wave action and water level fluctuations/water depth in the above species. The year 2016, with the highest and most stable water level, and high Achnanthidium spp. and medium high Gomphonema parvulum agg. abundances, was rather windy (Table 1, Figure 2 and Figure S22). Our ESEM photographs (Figure 14a,c) of diatom assemblages in reed (Phragmites australis) samples from two sites at which the level of exposure to wave actions was different suggest, however, that strong wave-mediated disturbances may have a rather devastating effect on the FT2 stalked taxa. It is also likely that the influence of wave action on the littoral plant and algal communities is different under different water levels [126].
Among vascular plants, even Phragmites australis, which is one of the most common helophytes in fishponds [20], forms stand only at those sites of the Dehtář fishpond that are at least partly protected against wave action. Nevertheless, its overall cover is negligible in comparison to other large fishponds without intensive wave-mediated disturbances (see, e.g., [12,13]). Similarly, Typha latifolia, very common even in fishponds with high nutrient levels [20,129], occurred in the Dehtář fishpond only as young plants during the exposure of muddy sediments. Although Gaberščik et al. [130] have demonstrated the negative impact of water level fluctuations on Phragmites australis in an intermittent lake, studies from fishponds show the positive effect of regular summer drainage on the reed bed regeneration [12], which is why we do not consider this factor to be responsible for the very limited extent of reed beds in the Dehtář fishpond. Lemnaceae, known for being very sensitive to streaming and wave action [131], occurred in small quantities only, usually within reed bed stands or in small inflows.
To a certain extent, wave action probably decreases the impact of high nutrient amounts and enables side-by-side coexistence of vascular plant and diatom species with contrasting nutrient demands, increasing the overall plant and diatom species diversity. According to Morris et al. [48], disturbances caused by wave action might even support rare and threatened plants. Also in our case, species that are the rarest in the Czech Republic grow preferentially on shores under strong wave action. In the case of Montia arensis, the Dehtář fishpond is one of the four sites in the Czech Republic and a single fishpond locality where this species occurs. Spergularia kurkae is more widespread [110], but the population from Dehtář is probably one of the richest in the country. Cerastium dubium is considered to be secondary in Southern Bohemia, colonising exclusively the sandy shores of large fishponds there [132].
Among diatoms, Staurosirella leptostauron var. dubia, Pseudostaurosira brevistriata, and P. elliptica are considered as species preferring nutrient poor habitats [55,122], while, e.g., Nitzschia palea, Fistulifera saprophila, and Planothidium frequentissimum are classified as the species of highly eutrophic to hypertrophic water bodies [58,122,133]. Most species on our list, however, occur between the two extremes, e.g., Cocconeis placentula, Navicula capitatoradiata, and Nitzschia archibaldii [122]. While the species with lower nutrient demands seem to be also highly tolerant of wave action and low water level (Figure 13b,c), those preferring highly eutrophic habitats occurred more frequently in places with fine sediment and organic detritus sedimentation and in the absence of wave-mediated disturbances. It is clearly visible not only when comparing findings from individual years and sites with varying disturbance intensities, but also when comparing two sites at FD (2019), differing in disturbance levels and also in quantities of, e.g., Pseudostaurosira spp. on the one hand and Planothidium frequentissimum on the other. Surprisingly, stalked FT2 diatom Gomphonema parvulum agg., predominating most of the 2019 samples, occurred in higher quantities on the more disturbed shore; however, its abundance varied between samples. Species from all the above groups occurred frequently in a single sample, but the samples dominated by species with low nutrient demand/higher disturbance tolerance were generally richer in species and had more balanced species abundances than the samples dominated by species with high nutrient demands/lower disturbance tolerance (Table S6, Figure 8). Analogically, the non-disturbed reed beds, e.g., stands of Phalaris arundinacea (zone 1), a strong competitor preferring nutrient-rich soils [134], were extremely species-poor (data not displayed). The phenomenon of mechanical disturbances and their positive influence on vascular plant species diversity in eutrophic wetlands has been already discussed (e.g., [18,20,68]).

A Real Biennial Fluctuation or a Trend of Hypertrophy Progression?
Potužák et al. [16] published substantially lower values of TN, TP, and chlorophyll-a concentrations from the Dehtář fishpond for the growing seasons of 2011 and 2012 than was detected during our research. The years 2016 and 2019 in particular show a great increase of trophy levels, potentially associated with the mineralisation of nutrients in the fishpond sediment during bottom exposure in dry summers. There is a high nutrient input related to the fish farming management, the surrounding agricultural landscape, and an outflow from a small sewage plant treating waste water from the surrounding settlements, but these factors have been constant over at least two decades. We suggest that sediment desiccation accompanied by nutrient mineralisation and the subsequent nutrient release into a water column may accelerate the influence of the above factors. Such a situation occurred at the Dehtář fishpond after the long dry phase in 2015 and the spring of 2016, followed by a water level increase in the summer of 2016. However, also a continual release of phosphorus from the submerged sediment due to high water temperatures is likely, as summer water temperatures measured by automatic sensor stations often reached 30 • C (see, e.g., [135,136] for the description of these processes).
There are no data available for an earlier period (2011-2012) on epiphytic diatom and macrophyte assemblages of the Dehtář fishpond. Therefore, we cannot be sure if the supposed continual increase in trophic levels in the Dehtář fishpond already caused sub-recent associated change in species composition of biotic communities or not. Only scarce data on selected vascular plant species of exposed pond bottoms were published several decades ago [69,70]. Thanks to these papers we know that, e.g., Juncus tenageia and Pseudognaphalium luteoalbum, the wetland annuals of nutrient-poor habitats [46,116] not found during our research, grew in the Dehtář fishpond in the past, but we have no information on what exactly preceded the disappearance of these species. On the one hand, potential hypertrophy-related changes may be very fast, particularly in algal assemblages [85]. On the other hand, frequent wind-mediated disturbances might have reduced the relative importance of nutrients in the hypertrophic ecosystem. Moreover, a parametes that could seem to be a consequence of recent global changes, namely a rather frequent prolonged summer drainage due to summer drought, has, according to our information, been occurring in the Dehtář fishpond for at least several decades [67]. Nevertheless, manuring and fish feeding based on the real needs of this aquacultural system would not only reduce the nutrient load but would also save money for the fish farm, as a large part of these inputs are not involved in the food chain and thus cannot be reflected in the fish production [37,68].

Diatoms and Macrophytes-Do Direct Relationships Exist?
Most published papers confirm the so-called "neutral substrate hypothesis", i.e., the absence of diatom assemblages strictly specific to a particular macrophyte species (plant substrate; [137]); our own early study from the Dehtář fishpond brought similar results [55]. On the other hand, we cannot say that there are no relationships between epiphytic diatom assemblages and their host plants. We proved that the ecological and life-history traits of the host plant species are an important factor influencing species richness, structure, and species composition of associated diatom assemblages. In some cases, the whole vegetation zone probably functions as a refuge for a certain epiphytic diatom species pool. This fact is reflected in the high similarity in diatom assemblages of the samples from different macrophyte species collected within a few square meters in a single zone. The importance of vegetation zonation for diatom assemblages is also mentioned by, e.g., Letáková et al. [60].
On the other hand, sometimes even the samples collected close to each other within a single zone showed very important differences in epiphytic diatom assemblages. For instance, our analyses revealed substantial differences between the samples from Phragmites australis and Phalaris arundinacea collected during the same sampling event at BSB. Similarly, diatom assemblages recorded in the samples from Spirodela polyrhiza and Lemna gibba differed markedly from those from reed bed species growing in the same stand: while Lemnaceae were densely colonised by a specialised diatom Lemnicola hungarica, only rarely occurring elsewhere [138], reed bed species samples were mainly predominated by Achnanthidium spp. and Cocconeis placentula. It is supposed that plant macro-morphology (i.e., stem branching, leave dissection) and surface microstructure (e.g., trichomes, spikes) have a selective effect on diatom assemblages, where e.g., species of a certain cell size, shape, or type of attachment may better cope with conditions offered by particular macrophytes [53,60].
The filtering of particular diatom size categories by macrophytes occurred also in some of our samples: while small (up to 12.5 µm) and middle-sized (up to 37 µm) diatoms such as Achnanthidium spp. predominated at FD in 2016, large species (>40-45 µm) such as Pinnularia spp. and Ulnaria spp. were common at BSB in 2016. In 2015, there was a mixture of medium-sized (e.g., Staurosira leptostauron) and large diatoms (e.g., Craticula cuspidata) at all the study sites. While the large and medium-sized diatoms seemed to be related particularly to small growing wetland annuals, shortly perennial species, and submerged plants, the small-sized species were more common on large wetland perennials (with the exception of Lemnicola hungarica on Lemnaceae). As among the small-sized species, the attached periphytic taxa dominated, one of the possible explanations might be that the periphytes need to attach to a stable structure in order to be able to persist in conditions of strong wave action. Remarkably, in 2019 when hydrological conditions at the time of diatom sampling were similar to those in 2015, but the year was markedly less windy than 2015 (Table 1, Figure 2 and Figure S22b,d), small-sized species such as Gomphonema parvulum occurred in high densities in FD samples from various types of macrophytes, mainly the non-reed bed species (Figure 14d, Table S6). G. parvulum is considered to be a species highly sensitive to mechanical disturbances due to its long stalks [139], and thus the course of the weather in 2019 with less frequent wave action could be the decisive factor in colonising the potentially unstable substrate such as Zannichellia palustris. On the other hand, epipelic species of the Navicula genus (FT5) were among the large-sized diatoms. Their higher species richness and/or quantities in the samples of low-growing macrophytes may be related to their occurrence on the surface of muddy sediments and to the active movement between the sediment and the plants [140]. This strategy would enable the diatoms to profit from their host plants by, e.g., gaining better access to light (the diatom may move upwards on the plant).

Importance of Research Intensity
After the end of the research grant in 2016, our data sampling in the Dehtář fishpond was limited. Therefore, the data on vascular plants after 2017 are rather incomplete, as we preferentially recorded taxa in the second and third zones, because many of them are of high conservation value. Despite this incompleteness, we included the data from 2019 when we also collected a limited number of diatom samples. These data serve partly for comparison with data from the research period of 2014-2016, but they particularly highlight the necessity of more intensive research in large water bodies. Although due to higher substrate moisture, the conditions in 2019 were even more suitable to plant development than in the very dry year 2015, during three short visits in 2019, we recorded less than 1/2 of the vascular plant and charophyte species listed for 2015 (Table S4). We excluded some difficult-to-access segments at the Dehtářský stream inflow from our 2019 research; however, this part of the Dehtář fishpond hosts only a few species that only grow there. Thus, it seems that a substantially higher time investment would be necessary to obtain data fully comparable with the research grant period.
For diatoms, the number of samples is crucial. There was a similarly significant decrease in the total species number in 2019 as compared to 2015 and especially to 2016 (the highest amounts of samples). However, the most important reduction of information was related to the large-scale distribution patterns of diatoms, because all 2019 samples came only from FD. Another possible source of bias in diatom data is related to differences in the zones sampled in particular years: while in 2015 and 2019, all the samples were collected on the transition between the second and the third zone, in 2016, the first and the second zones were sampled. Therefore, somewhat different sets of diatom substrates were sampled each year. The sampling of the same plant species at each site and each year would allow us to formulate more general conclusions; however, it was not realistic due to the differences in macrophytic species composition between the sites and years and because not all of the plant species occurring at particular sites were always available growing in water. Thus, some rare macrophyte species could be collected only once. In some cases, our sample was probably the first one of the given macrophyte species (e.g., Veronica scutellata, Alopecurus geniculatus) that was ever analysed for epiphytic diatoms.

Conclusions
Macrophytes and diatoms exhibit similar responses to some of the basic, usually locally acting environmental factors, such as water depth, water level fluctuations, nutrient amounts, water transparency, type of bottom sediment, occurrence of organic detritus, morphology of littoral zone, and mechanical disturbance intensity. Locally as well as regionally, important environmental and management-related factors such as the yearly course of weather conditions, type of management (e.g., fish farming), and landscape structure usually act in close synergy.
Our study of the Dehtář fishpond offers a complex view of multiple factors and their interactions within a large hypertrophic water body. For instance, wind-mediated wave action is related to specific weather conditions, landscape structure, and fishpond morphology, but the controlled manipulation of the water level may largely modify its influence on macrophyte and diatom species in particular littoral habitats. Similarly, the timing and intensity of water level fluctuations in a water body are basically determined by weather conditions, the hydrology of a local catchment basin, and the type of surrounding landscape. In fishponds specifically, however, the overall hydrological conditions may also be greatly modified by an intended drainage before fish harvesting and the continual filling of a fishpond after restocking. Due to the synergy between hydrological, climatic, and management factors, the functioning of the littoral zones of these fishponds resembles temporary wetlands, as we have also shown with our complex data. Our results also reveal other stress factors, such as the eutrophication and mechanical disturbances with a joint effect on water transparency and, consecutively, on macrophyte and microalgal assemblages. Yet, the Dehtář fishpond is an important hotspot of wetland biodiversity and a refuge of threatened taxa.
The study period was characteristic of various extremes in climate and related local factors such as water level and substrate moisture fluctuations, with associated responses in plant and diatom assemblages. As such extremes are supposed to be more frequent in the future due to the progressing global climate change, our study from the Dehtář fishpond may be considered as the basis for predictive model constructions in similar hypertrophic water bodies, including for proposals of mitigation measures preventing accelerated habitat deterioration.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/w13111569/s1: Figures S1-S15: Photographs documenting the environment, selected vascular plant species and their assemblages, and fish farming practices in the Dehtář fishpond, Figures S16-S22: Line graphs documenting a course of selected environmental factors in the Dehtář fishpond, Table S1: Summary of the factors used in NMDS ordinations of macrophyte and diatom species assemblages; Table S2: Basic physico-chemical parameters of water; Table S3: List of epiphytic diatom samples; Table S4: List of vascular plant species and charophytes along with their frequencies, quantities and classification to species groups; Table S5: Frequency of occurrence and cumulative numbers of diatom taxa in the samples from the Dehtář fishpond; Table S6: List of diatom taxa found in particular samples of epiphyton of the Dehtář fishpond.