Triggers for the Impoverishment of the Macroinvertebrate Communities in the Human-Impacted Rivers of Two Central European Ecoregions

Human activity triggers negative alternations in river habitats, including changes to the physical and chemical parameters of the water, its hydromorphological features and the introduction and spread of invasive alien species. These modifications are expected to be intensified by climate change. Eight rivers in one of the most urbanised and industrialised regions in Europe, i.e. the Upper Silesian Coal Basin, were surveyed in order to explain the impact of anthropopressure on the distribution of macroinvertebrates. Conductivity, altitude, hydromorphological transformations, hardness, the organic matter content and certain fractions of benthic sediments significantly affect (p < 0.01) the occurrence of macroinvertebrates in Central European rivers. Our results proved that the hydromorphological transformation of watercourses, which is expressed by the relevant indices, is one of the most predictive factors that contribute to the distribution of macroinvertebrates. Anthropogenic inland waters that have been salinised by the discharge of hard coal mine waters create new habitats for brackish and marine species that replace native freshwater species. An increase in salinity causes an impoverishment of macroinvertebrate biodiversity therefore all possible actions should be taken to reduce the anthropogenic salinity of inland waters. Secondary saline rivers may prove to be prescient for climate-induced changes to river macroinvertebrates.


Introduction
Global biodiversity is now at the beginning of the sixth mass extinction of species. Projections indicate that approximately 75% of species may be lost due to this human-induced, global ecological crisis (Ceballos et al. 2015;Payne et al. 2016;Briggs 2017). Ten key factors, the so-called big killers (Maxwell et al. 2016), have been implicated in the accelerating the loss of biodiversity worldwide: overexploitation, agricultural activity, urban development, invasion and disease, pollution, system modification, climate change, human disturbance, transportation and energy production. Almost all of these "big killers" impact species dependent on freshwater environments-environments themselves that are the most threatened ecosystems on earth (Nieto et al. 2017).
In freshwater environments, nearly 10% of all known species are present (Winemiller 2018). This includes approximately one-third of known vertebrates and more than 2% (> 150 000 species) of global invertebrate biodiversity (Strayer 2006;Dudgeon et al. 2007). Changes in the physical and chemical characters of water along with the hydromorphological (e.g. riverbeds, water regime, connection with floodplain) degradation of rivers are direct threats to freshwater invertebrate diversity worldwide. In addition, anthropogenically derived point and diffuse pollution sources impact the diversity and abundance of macroinvertebrates. Anthropogenic salinisation is increasingly recognised as contributing to the degradation of water quality worldwide and represents a growing ecological threat to rivers and streams (Cañedo-Argüelles et al. 2013).
Hydromorphological transformations of most streams and rivers in Europe constitute one of the major threats to the habitats of mussels (Geist 2014). The removal of weeds and dredging of riverbeds, and above all channelisation and dams, alter the characteristics of aquatic ecosystems, disrupt the natural meta-population structure and contribute to local macroinvertebrate extirpations (Aldridge 2000;Cosgrove and Hastie 2001;Geist and Kuehn 2005).
Freshwater environments have long been experiencing the dangers of human activity, but in recent years, this process has been intensifying as urbanisation progress and the human population increases (Amoatey and Baawain 2019; Kumaraswamy et al. 2020). In recent years, anthropogenic salinisation has been considered one of the key factors contributing to the decrease of water quality worldwide (Cañedo-Argüelles et al. 2013). Nowadays, the major contributors to freshwater pollution include wastewater and industrial effluents, metals and metalloids, flame retardants, polycyclic aromatic hydrocarbons, nutrients and persistent organic pollutants, herbicides, pesticides, pharmaceuticals and illicit drugs and endocrine disruptors. In addition, new emerging pollutant sources in freshwater such as microplastics and engineered nanoparticles are being discovered (Amoatey and Baawain 2019). These changes in pollution affect the ecosystem functions of clean water for humanity and also affect aquatic biota (including macroinvertebrates) and may affect their ecological services (Cañedo-Argüelles et al. 2016;Lopes-Lima et al. 2017;Schmeller et al. 2018;Amoatey and Baawain 2019;Berger et al. 2019;Kumaraswamy et al. 2020). For example, heavy metals and the organochlorine insecticide DDT can affect the calcification and formation of mussel shells (Pynnönen 1995;Lopes-Lima et al. 2017). Phosphorous and nitrogen concentrations are higher and are increasing in most European wetlands (Galloway et al. 2008;Douda 2010;Grizzetti et al. 2011). The deposition of fine-sediments and the introduction of nutrients from agricultural run-off are considered major threats to the most endangered freshwater macroinvertebrates such as mussels (Geist and Auerswald 2007;Lopes-Lima et al. 2017). Anthropogenic (secondary) salinisation is responsible for the degradation of water quality on a worldwide scale and results in biological changes in ecosystems, mainly in freshwater biota (Bäthe and Coring 2011;Kang and King 2012;Arle and Wagner 2013;Bąk et al. 2020;Sowa et al. 2020). In addition, dioxins, pharmaceutical compounds and other chemicals can have chronic and acute effects on aquatic macroinvertebrates (Lopes-Lima et al. 2017).
Climate change impacts abiotic factors including temperature and precipitation, which in turn, affect the functioning of aquatic ecosystems including the rate of reproduction and feeding of freshwater organisms (Vörösmarty et al. 2000;Schmeller et al. 2018). Both climate change and the anthropogenic degradation of water ecosystems contribute to the reduction in the quality and quantity of freshwater resources. Together with the increase in human population, they intensify the large-scale phenomenon of the impoverishment of biodiversity of freshwater environments (Vörösmarty et al. 2000;Malmqvist and Rundle 2002;Dudgeon et al. 2007). What is more, because of the interaction between global climate change and the anthropogenic alterations of the hydrologic cycle, secondary salinisation may become more intense in the coming decades (Neubauer and Craft 2009). Anthropogenic salinisation can degrade wetlands, affect a variety of the processes and dynamics of water environments and the ecosystem services that are provided by freshwater (Herbert et al. 2015).
However, the scientific understanding of the mechanisms by which increasing salinity destroys freshwater ecosystems is poor, which makes protecting freshwater habitats difficult (Cañedo-Argüelles et al. 2016). Referring to the above arguments, research on the factors that impact the changes in freshwater biodiversity in the age of climate change and environmental emergency have substantial value in the policy and decisionmaking processes of administrations and other stakeholders about the risk of climate change and the potential impact on human well-being (Lioy and Smith 2013).
Europe is divided into 25 ecoregions on the basis of the fauna living in European inland waters. The European Union Water Framework Directive (EU WFD) (Directive 2000) uses abiotic variables to classify streams and rivers into types. For rivers, the Directive fixed the typology, i.e. "System A" typology is defined by the ecoregions and size based on the catchment area, catchment geology and altitude. Therefore, based on the "System A" typology, the homogenous water bodies for the sampling sites were selected in accordance with the above requirements of the Directive including the ecoregions (Directive 2000).
The objectives of our survey were to analyse the structure of the macroinvertebrate communities in rivers that are located in one of the most industrialised and urbanised regions in Europe and to determine the most predictive environmental factors that affect the structure of macroinvertebrate communities. The physical, chemical and hydromorphological human pressure on the macroinvertebrate communities in various abiotic types of rivers have rarely been studied at the same time including the most salinised river in the world by the discharge of underground coal mine water. We hypothesise that the physical and chemical anthropogenic transformations are more important than the hydromorphological transformations in explaining the distribution of macroinvertebrates in the selected human-impacted Central European rivers.

Study Area and Site Selection
The research was carried out from 2017 to 2018 in Upper Silesia and adjacent areas (Southern Poland), which comprise part of the Central Plains (Ecoregion 14) and the Carpathians (Ecoregion 10) according to the division of the EU WFD (Directive 2000) (Fig. 1). The initial field recognition (site visits) and selection of sampling sites took place in 2016 prior to the fundamental field survey. The degree of human pressure was determined based on preliminary field surveys and the reports on ecological status of rivers obtained from the Voivodeship Inspectorates of Environmental Protection, Katowice and Cracow, Poland.
The research area included both protected areas and the most urbanised and industrialised region, i.e. the Upper Silesian Coal Basin (USCB), which is one of the largest hard coal mining areas in the world. Hard coal has been exploited in the USCB since 1740 (Dulias 2016). About 65 underground coal mines have been discharging the saline mine waters of Carboniferous rocks that have high concentrations of salts (mainly chlorides and sulphates) via the coal mine dewatering systems into the tributaries of the Vistula and Oder rivers, which caused their deterioration by the last decades of the twentieth century. Nowadays, about 35 hard coal mines are still active in the USCB. However, even after hard coal mines are closed, the saline mine waters are still discharged into the surface waters including rivers in order to eliminate the flooding risk to adjacent coal mine workings (Strozik 2017).
Eight rivers of four abiotic types, which are under different anthropogenic pressures, were selected in the catchments of the Vistula and Oder rivers within two ecoregions, i.e. the Central Plains and the Carpathians (Table 1). Two sampling sites (one in the upper course and one in the lower course) were selected for each river depending on the degree of anthropopressure. The sampling site in the upper course of each river (reference sampling site) was characterised by the absence or smaller scale of hydromorphological transformations including the degree of regulation of the river, more natural forms of land use and the lower values of the physical and chemical parameters of the water (e.g. conductivity, salinity or the concentrations of nutrients) compared with the sampling site located in the lower course. A total of 16 sampling sites are located according to the gradient of human pressure (Table 1). The four sampling sites (the upper courses of the Centuria, Wiercica, and Vistula and the lower course of the Korzenica River) are located within the protected areas of nature reserves and Natura 2000, whereas the sampling sites of the Mleczna and Bolina rivers are located within the most urbanised and industrialised parts of the USCB.

Environmental Surveys
Prior to the macroinvertebrate sampling, water samples and bottom sediments were collected from each sampling site. The physical and chemical parameters of the water, i.e. pH, temperature and the concentration of oxygen were measured in the field using CO-401 Elmetron and HI-9811-5 Hanna Instruments portable meters. The concentrations of nutrients, iron, chlorides, sulphates, alkalinity, magnesium, calcium and total hardness were analysed in the laboratory. The concentration of ammonium (Nessler method); nitrites (adaptation of the EPA diazotization method 354.1.); nitrates (cadmium reduction method); phosphates (adaptation of the ascorbic acid method); iron (phenanthroline method); sulphates (turbidimetric method); calcium (adaptation of the oxalate method); magnesium (adaptation of the calmagite method) in the water and hardness (adaptation of the EPA method 130.1.) were determined using Hanna Instruments (USA) reagents and portable photometers. The concentration of chlorides was analysed according to the titrimetric determination with a mercury nitrate solution ( m e r c u r i m e t r i c d e t e rm i n a ti o n a g a in s t 1 . 5diphenylcarbazone) and alkalinity according to the acidimetric titration against phenolphthalein or mixed indicator using Merck reagents.
Electrical conductivity (EC) and total dissolved solids (TDS) were measured in the field using a Multi 3410 WTW meter and the results were converted to the salinity according to Piscart et al. (2005). The width, depth and flow velocity of the rivers were measured in the field (Hauer and Lamberti 2007). The organic matter content (%) in the bottom sediments was determined according to the methods of Myślińska (2001). The grain size composition of the bottom sediments was determined using both the sieve and aerometric methods.
The anthropogenic pressure on the watercourses and adjacent land use was assessed in accordance with the field survey method, i.e. the Hydromorphological Index for Rivers (HIR) and the calculation of two numerical metrics, i.e. the Hydromorphological Diversity Index (WRH) and the Hydromorphological Transformation Index (WPH) (Szoszkiewicz et al. 2017(Szoszkiewicz et al. , 2020. The HIR method, which was prepared mainly on the basis of the British River Habitat Survey (RHS), was adapted to Polish conditions (Szoszkiewicz et al. 2017(Szoszkiewicz et al. , 2020. These studies included both fieldwork and analyses of the Geographic Information Systems data and remote sensing materials. High values of the WRH and HIR indices reflect a large number of the natural characteristics of rivers and the adjacent land use, while a high value of the WPH reflects strong anthropopressure on a river and the adjacent land use. The macroinvertebrate samples were collected using a quadrat frame and hydrobiological net with a 0.4-mm mesh sieve from different microhabitats (bottom sediments and macrophytes) according to the standard  methods. A total of 367 biological samples were collected from different types of microhabitats. The samples were sorted and preserved in 80% ethanol. The collected material was identified to the lowest possible level according to Rozkošný (1980), Kołodziejczyk and Koperski (2000), Eggers and Martens (2001), Tończyk et al. (2013), Piechocki and Wawrzyniak-Wydrowska (2016), Mauch (2017) and Rudolph et al. (2018) and counted under a stereoscopic microscope (OLYMPUS SZ61 SZ2-ILST). The wet weight of the macroinvertebrates was measured on a laboratory balance with an accuracy of 0.001 g after removing excess water by drying the macroinvertebrates on filter paper. The density and biomass of the macroinvertebrates was estimated as the number of individuals per square metre and gramme (g), respectively.

Statistical and Zoocenological Analyses
The zoocenological analysis of the macroinvertebrate communities was based on the Shannon-Wiener index H' (Hauer and Lamberti 2007). The significance of the differences in the median values of the environmental variables, the number of taxa, the density and biomass of macroinvertebrates and the values of the Shannon-Wiener index H' between the rivers was calculated using the Kruskal-Wallis one-way ANOVA and Dunn's multiple comparison post hoc tests. The values of the biological data and environmental variables did not reveal a normal distribution according to the Lilliefors test of normality, which justified the use of a non-parametric test. The statistical analyses were performed using Statistica version 13.1.
Canonical ordination analyses for relating the taxonomic composition of macroinvertebrates to the environmental variables were performed using CANOCO for Windows version 4.5 (Ter Braak and Smilauer 2002). The appropriate type of analysis was selected to analyse the biological data using detrended correspondence analysis (DCA) and the length of the gradient. The gradient length exceeded 3 SD (the standard deviation), therefore a unimodal direct ordination, canonical correspondence analysis (CCA) with a forward selection was used to reduce the large set of environmental variables (a maximum gradient length of 5.161). The following environmental variables were included in the CANOCO analysis: the physical and chemical parameters of the water, altitude, stream gradient, the width, depth and flow velocity of the rivers, the area of catchments, the grain size composition of the bottom sediments, the organic matter content in the bottom sediments and the values of the hydromorphological indices. Macroinvertebrate taxa that occurred at fewer than 10% of the sampling sites were excluded from the statistical analyses following a preliminary exploration of their influence in the initial DCA analysis (McCune and Grace 2002). Both the biological and environmental data were log-transformed. The statistical significance of the relationship between the biological data and the environmental variables was evaluated using the Monte Carlo permutation test (499 permutations) (Ter Braak and Smilauer 2002).

Results
The values of the WRH, which reflect the number and abundance of the natural features of the rivers and the adjacent land use, were relatively high (up to 88.5) for the sampling sites that are located within the protected areas: the nature reserves and Natura 2000 (the upper courses of the Vistula and Wiercica rivers) ( Table 2). The relatively high values of the WPH reflected a high degree of the anthropogenic modification of the habitats and were recorded for the river sections within the city located outside the most urbanised and industrialised parts of the USCB (abiotic type 12: the Dziechcinka and Vistula rivers) ( Table 2). Total hardness, alkalinity, the concentration of nutrients, calcium, magnesium, salinity and the parameters associated with salinity, i.e. EC, TDS, the concentrations of chlorides and sulphates were very high for the sampling sites located in the most urbanised and industrialised parts of the USCB (the Bolina and Mleczna rivers) compared with the sampling site located outside of the USCB and the protected areas (Table 3).
The differences in the median values of most of the physical and chemical parameters of the water and the morphological features between the rivers were significant (p < 0.01) (the Kruskal-Wallis one-way ANOVA and multiple comparison post hoc tests) (Tables 2 and 3 and Appendix Table 4). The Kruskal-Wallis one-way ANOVA and the multiple comparison post hoc tests test revealed statistically significant differences (p < 0.05) in the median number of taxa, density and the median values of the Shannon-Wiener index H' between the rivers (Fig. 2).   Table 3 The physical and chemical parameters of the water and the results of the Kruskal-Wallis one-way ANOVA and multiple comparison post hoc tests (superscript denotes significant differences between the rivers, p < 0.01 The Shannon-Wiener index H' summarises the species richness, evenness and the density of the taxa in samples including rare and not numerous species. Therefore, higher median values of H' were recorded for the rivers with a higher number of the natural features and the adjacent land use, which was reflected by the maximum values of the WRH (less human-impacted rivers) (Fig. 2). A higher median number of taxa and values of the H' index were recorded for the less humanimpacted rivers within the same abiotic types, i.e. for the Centuria River within abiotic type 5, for the Mitręga River within abiotic type 6, for the Vistula River within abiotic type 12 and for the Wiercica River within abiotic type 17 (Fig. 2). The lowest median number of taxa and values of the Shannon-Wiener index H' were obtained for the most salinised rivers, i.e. for the Bolina and Mleczna rivers, even though the maximum density of macroinvertebrates was relatively high in the Mleczna caused by invasive alien species (Fig. 2).
A total 143 macroinvertebrate taxa were recorded in the studied rivers (Appendix Table 5). The number of macroinvertebrate taxa ranged from 40 in the Bolina River to 99 in the Wiercica River. Among them, three invasive alien species were recorded, i.e.  Bo-Bolina River, Ce-Centuria River, Mi-Mitręga River, Ml-Mleczna River, Dz-Dziechcinka River, Vi-Vistula River, Ko-Korzenica River, Wi-Wiercica River the Korzenica River) and Physa acuta (in the Mleczna River). One species of unionid mussel, i.e. Anodonta anatina, which made up a small percentage in the macroinvertebrate communities, only occurred in the rivers of abiotic type 6 (the Mitręga, the upper course of the Mleczna).
Canonical correspondence analysis (CCA) based on the macroinvertebrate data and environmental variables showed that the first two axes explained 17.4% of the variance in the biological data and 71.8% of the variance in the biological data and environment relationships. The results of CCA proved that many environmental factors affect the occurrence of macroinvertebrates simultaneously. In addition to the physical and chemical parameters of the water, the geographical and habitat factors as well as the hydromorphological parameters of rivers were also important. Selection of the environmental variables based on the forward selection (FS) results and the marginal and conditional effects showed that EC, altitude, the values of the WRH, hardness, medium sand (particles of 0.5-0.25 mm), pebble (particles above 20 mm) and the organic matter content in the bottom sediments were the parameters most associated (statistically significant, p = 0.002) with the distribution of the macroinvertebrates (Fig. 3). According to the forward selection, among the environmental factors, EC had the strongest effect on the distribution of the macroinvertebrates including the invasive alien species. The anthropogenic salinisation of the Mleczna and Bolina rivers of up to 33.6 PSU (23 300 mg dm −3 TDS, salinity similar to the North Sea) negatively impacted the macroinvertebrate communities directly. Gastropod Ancylus fluviatilis and Bythinella sp., the crustacean Gammarus fossarum, dipterans (Blephariceridae, Dixidae, Athericidae) and most of ephemeropteran, plecopteran and trichopteran taxa (e.g. Odontocera albicorne, Glossosomatidae, Odontoceridae, Polycentropodidae, Sericostomatidae, Rhyacophilidae, Perlidae, Perlodidae, Nemouridae, Heptageniidae, Leuctridae) were affected by altitude and the natural features of the rivers, which was reflected by the values of the WRH and pebbles (Fig. 3). Gammarus jazdzewskii, Leptophlebiidae, Limnephilidae, Polycentropodidae, Pisidium sp., Sphaerium sp. and Leptoceridae were the taxa positively correlated with an increasing value of the WRH or by both the values of the WRH and medium sand. The abundance of gastropods (e.g. Bithynia tentaculata, Planorbarius corneus, Gyraulus crista and Radix balthica) was associated with the hardness and a lower altitude, whereas the abundance of Oligochaeta, Chironomidae or Tabanidae was associated with the organic matter content in the bottom sediments. The brackish and marine invasive alien species, i.e. Gammarus tigrinus, Potamopyrgus antipodarum and the alien gastropod Physa acuta occurred at the sampling sites with a higher EC. The relationship between the taxonomic composition of the macroinvertebrates and the environmental variables was statistically significant (Monte Carlo test of significance of the first canonical axis: p = 0.002, F-ratio = 44.88; test of significance of all of the canonical axes: p = 0.002, F-ratio = 16.42).

Discussion
This study proved that the physical, chemical and hydromorphological anthropogenic transformations are important in explaining the distribution of macroinvertebrates in human-impacted Central European rivers. Our results, which showed a decrease in macroinvertebrate biodiversity in the most salinised rivers reflected by the lowest median values of the Shannon-Wiener index H' for the macroinvertebrate communities, are consistent with numerous studies that have been carried out in secondary saline rivers worldwide (Battaglia et al. 2005;Piscart et al. 2005;Bäthe and Coring 2011;Arle and Wagner 2013;Cañedo-Argüelles et al. 2013Ladrera et al. 2017). According to Pinder et al. (2005), the threshold value of salinity is 4100 mg dm −3 TDS, above which the biodiversity begins to decrease. Even a small increase in salinity up to about 800 mg dm −3 TDS inhibits macroinvertebrate reproduction and reduces their growth, while at a salinity of approximately 9000 mg dm −3 TDS, their osmoregulatory functions begin to fail (Hart et al. 1991;James et al. 2003). The juvenile forms of these organisms (e.g. larvae such as glochidia) are generally more sensitive to salinity and their occurrence is limited to the maximum ranges of salinity tolerance for adult macroinvertebrates (Hart et al. 1991;Blakeslee et al. 2013;Herbert et al. 2015;Lopes-Lima et al. 2017). Nevertheless, a slight increase in EC to 300-500 μS cm −1 may increase species richness substantially (Kefford et al. 2011). The salinity gradient that was obtained in our survey proves that with increasing EC, salinity-sensitive taxa are replaced by eurytopic and euryhaline taxa including invasive a l i e n s p e c i e s . T h e b r a c k i s h a n d m a r i n e amphipodGammarus tigrinus is an example of such a species which, among others, because of the saltwater pollution of rivers, has had colonisation success in Europe and is still expanding its range (Lewin et al. 2018;Rewicz et al. 2019).
An increase in the concentration of salts in rivers induces physiological stress, and ultimately, can cause major changes in the macroinvertebrate communities and the functions of the ecosystem. For example, water with an EC higher than 1500 μS cm −1 is not suitable for irrigating most crops and extirpates many species of freshwater invertebrates, while water of EC of 2500 μS cm −1 , according to standards established in some countries, is not suitable for drinking or for use in some areas of industry (Cañedo-Argüelles et al. 2013, 2016). Freshwater mussels are sensitive to increasing salinity, whose values higher than 1400 μS cm −1 may negatively affect their metabolic processes and cause mortality (Blakeslee et al. 2013). As a consequence, rivers with an elevated salinity lose their ecosystem services. Our results revealed that highly saline rivers, e.g. the Bolina River and the lower course of the Mleczna River were devoid of mussels. Hence, the anthropogenic salinity of waters has a substantial impact not only on freshwater biota, but also on the goods and services that rivers provide to people, and therefore, they can be associated with high economic costs (Cañedo-Argüelles et al. 2013, 2016. In addition, rivers with the discharge of underground salty mine waters can cause local environmental contamination due to excessive salt loads as well as heavy metals and sometimes radioactive substances (Harat and Grmela 2008).
The anthropogenic salinisation of rivers (especially the strong salinity that is associated with the discharge of salty mine waters, but also with the agents that are used to de-ice roads) and its impact on freshwater biota illustrates the consequences of climate change. The consequences of climate change will not only increase the number of violent phenomena and change the volume and time of precipitations, but will also result in an increase in the average annual temperatures and droughts, which may contribute to an increasing salt concentration in freshwater environments (Nijssen et al. 2001;Fischer et al., 2007;Neubauer and Craft 2009;Herbert et al. 2015;Cañedo-Argüelles et al. 2016;Maxwell et al. 2016;Berger et al. 2019). Therefore, rivers that have been polluted with salty mine waters, e.g. the Bolina, the Gostynia, the Mleczna and the Soldevila rivers (Ladrera et al. 2017;Lewin et al. 2018;) may constitute an appropriate model for research on the impact of climate change on freshwater organisms in Europe. Our results showed that among the rivers that are impacted by the discharge of hard coal mine waters, the Bolina River is currently the most salinised river in the world. Undoubtedly, comprehensive guidelines are required for monitoring, regulating, managing and dealing with the salinity of freshwater (Schuler et al. 2019). However, in order to prevent or limit climate change, a rapid reduction of coal mining is necessary and the use of energy from this source should be discontinued. This is particularly important in a country such as Poland, which was responsible for producing 61.6 million tonnes of hard coal in 2019, which accounted for 95% of the total production in the European Union (Eurostat Statistics Explained 2020). The output and processing of hard coal affects both the anthropogenic salinisation of inland water and climate change. Hard coal, which is a natural resource, is extracted from the underground seams and then to the earth's surface throughout the world. Poland hosts the largest number of coal mines (35) followed by Spain, Germany and Bulgaria (Alves Dias et al. 2018). According to McGinley (2011), annual world coal production is projected to grow to seven billion tonnes by 2030, and China will account for half of the increase. Poland is in the 10th largest world coal consumer. The enormous consumption of billions of tonnes of coal for generating electricity pollutes the Earth's atmosphere with greenhouse gases including SO 2 , CO 2 , NO x and ozone (McGinley 2011). Coal contributes more climate change-inducing greenhouse gases than any other fossil fuel. In the USCB, the extraction of hard coal that is sent to the earth's surface requires that underground mine waters, which carry high concentrations of salts, mainly chlorides and sulphates and also heavy metals and radionuclides, be pumped out. In order for the underground mine waters not to flood the mine workings, they are pumped out to the surface through mine dewatering systems and are initially discharged into central settling ponds and then directed to surface waters (waterbodies) or directly to the rivers (Strozik 2017;Halabowski et al. 2019). Therefore, the rivers of the Vistula and Oder catchments are strongly affected by the highly mineralised coal mine waters originating from the dewatering systems of hard coal mines. This result showed that among the rivers that are impacted by the discharge of hard coal mine waters, the Bolina River is currently the most salinised river in the world. The Bolina River is an example of a saline tributary of the upper Vistula River . Moreover, the increase in the consumption of hard coal means that it has to be extracted from deeper layers of the earth. The deeper the hard coal seams that are exploited, the more saline waters occur in underground layers that are then pumped out to the surface watercourses of the USCB. What is more, even if hard coal mines are closed for economic reasons, pumping out salty mine waters to the surface is continued so that these waters do not flood the mine workings and adjacent areas. Although the output of hard coal has decreased in the USCB, the discharge of saline waters (325 thousand m 3 per day) and the load of salts (355 tonnes per day) from hard coal mines has not followed this trend (Strozik 2017). All mine waters with an average salt concentration of about 10 kg m −3 are discharged into the tributaries of the catchments of the Vistula and Oder including the Bolina and the Mleczna and then affect river biota.
The transformations of flowing waters that are caused by regulating rivers, constructing dams and other hydrotechnical devices and intentionally or accidentally introducing alien species as well as climate change all negatively affect the biota and natural abiotic elements of rivers (Allan and Castillo 2009;Herbert et al. 2015). It also intensifies the effects of salinity. Population growth (which involves an intensification of agriculture) and climate change will certainly intensify this problem (Vörösmarty et al. 2000). Therefore, it has become n e c e s s a r y t o d e t e r m i n e t h e d e g r e e o f t h e hydromorphological changes and to measure them using the indices for determining the quality of river habitats. Our research showed that the natural features of rivers and the adjacent land use, which are reflected by the values of the WRH, were the most important factors that affected the distribution of macroinvertebrates in the rivers. Both the most anthropogenically transformed rivers and eurytopic taxa were negatively associated with the WRH. Several other studies (Hering et al. 2006;Lewin et al. 2013Lewin et al. , 2014Paller et al. 2014;Lewin et al. 2015) have reported similar correlations. Our survey revealed the relatively high values of the WPH reflected a high degree of the anthropogenic modification of the habitat for the river sections of the Dziechcinka and Vistula (flysch streams) within the city located outside the most urbanised and industrialised parts of the USCB. This surprising pattern may be explained by the location of these sections of the rivers in the centre of the Wisła city (famous sports and recreational resort with high volumes of tourist traffic in the Beskid Śląski Mts.). Strong water flow and quick increase in the water level during a rainfall in these sections determined heavily transformation of the riverbeds (regulation, concrete reinforcement of the bank) and land use to protect the city from flooding. Therefore, higher values of the WPH were recorded for these sampling sites than for other sampling sites. However, despite the fact that the substantial impact of hydromorphological transformation (number and abundance of the natural features of the rivers and the adjacent land use) that has been determined using the appropriate indices has been proven, these environmental variables are still not taken into account by researchers that are trying to explain the occurrence of macroinvertebrates in watercourses, especially in human-impacted rivers. In our research, salinity and altitude proved to be more decisive for the occurrence of macroinvertebrates in the rivers. This is justified in the light of the above arguments as an effect of secondary salinisation on freshwater macroinvertebrates, while altitude, which reflects geology, is mentioned as a very important factor that is associated with the structure of the benthic macroinvertebrate communities in rivers (Rawer-Jost et al. 2004;Šporka et al. 2009;Lewin et al. 2013Lewin et al. , 2014. These surveys found a relatively high content of organic matter in bottom sediments at the sampling sites that were located at different distances below dam reservoirs. In addition, some sampling sites were located within agricultural areas. According to the results of the CCA, some macroinvertebrate taxa were affected by the organic matter content in the bottom sediments. This finding confirmed that the occurrence of eurytopic species is positively correlated with a high organic matter content in the bottom sediments (Allan and Castillo 2009;Ward et al. 2002). What is more, our results also showed that the extremely high organic matter content in the bottom sediments consisted mainly of mining waste in the rivers into which there was a discharge of hard coal mine waters. This phenomenon can be explained by the increased sedimentation in saltwater and the flocculation of sediment that is caused by the divalent cations in saline water (e.g. Ca 2+ and Mg 2+ ) because of the aggregation of suspended matter (Sholkovitz 1976;Grace et al. 1997;Meiggs and Taillefert 2011;de Nijs and Pietrzak 2012). This is consistent with our results, because in both the Bolina and Mleczna rivers, high concentrations of calcium, magnesium and total hardness were recorded, and those higher values significantly (p = 0.002) affected the occurrence of eurytopic macroinvertebrate taxa.

Conclusions
Our results showed the crucial role both of the physical, chemical and hydromorphological anthropogenic transformations in explaining the gradient in the distribution of macroinvertebrates in the studied rivers, which are under different degrees of human impact. In contrast to the assumed hypothesis, these results proved that the hydromorphological transformation of human-impacted Central European watercourses, which is expressed by the relevant indices and is often overlooked or underestimated by scientists, is one of the most predictive factors that contribute to the distribution of macroinvertebrates including invasive alien species. The application of the HIR methods and the calculation of two indices, i.e. the WRH and WPH enabled the impact of human pressure to be assessed. Our results also showed that both the WRH and WPH indices are indispensable tools for assessing the degree of human pressure within the catchment area of the human-impacted Central European rivers. The WRH not only reflected the degrees of the natural features of a riverbed and the land use adjacent to a river but also the relationships between the habitat features and the structure of the macroinvertebrate communities. Therefore, in future surveys, we recommend that the Hydromorphological Diversity Index be applied and that the hydromorphological status of rivers based on the HIR methods be considered to be drivers that determine the structure of the macroinvertebrate communities in rivers. Our survey revealed that anthropogenic salinity and the hydromorphological transformations of river habitats reflect a high level of human pressure, which triggers the occurrence and spread of the alien gastropod Physa acuta and brackish and marine invasive alien species including Gammarus tigrinus and Potamopyrgus antipodarum. We can conclude that the anthropogenic salinised inland waters due to the discharge of hard coal mine waters create new habitats for brackish and marine species that replace the native freshwater species. Moreover, these salty habitats can illustrate future changes that are connected with an increase in the human population and climate change and can be used to model the impact of these changes on freshwater biota. In addition, rivers impacted by salty mine waters discharge lose their ecosystem service indirectly (e.g. unionid mussels, which are effective filtrators, are unable to support the process of selfpurification of waters) and therefore such waters are not suitable for human use. Among the rivers that are impacted by the discharge of hard coal mine waters, the Bolina River, which is currently the most salinised river in the world, is an example of such a watercourse. However, all possible actions should be taken to minimise the anthropogenic salinity of inland waters including implementing comprehensive legal regulations.     Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.