Abstract Introduction Methods Results Discussion Conclusion Acknowledgements References Appendix 1 About the Author(s) Nabil Majdi Department of Animal Ecology, Faculty of Biology, University of Bielefeld, Bielefeld, Germany Lizaan de Necker Unit for Environmental Sciences and Management, North-West University, Potchefstroom, South Africa South African Institute for Aquatic Biodiversity, Makhanda, South Africa Hendrika Fourie Unit for Environmental Sciences and Management, North-West University, Potchefstroom, South Africa Ané Loggenberg Unit for Environmental Sciences and Management, North-West University, Potchefstroom, South Africa Edward C. Netherlands Unit for Environmental Sciences and Management, North-West University, Potchefstroom, South Africa Jasmine Bunte-Tschikin Department of Animal Ecology, Faculty of Biology, University of Bielefeld, Bielefeld, Germany Walter Traunspurger Department of Animal Ecology, Faculty of Biology, University of Bielefeld, Bielefeld, Germany Gerhard C. du Preez Unit for Environmental Sciences and Management, North-West University, Potchefstroom, South Africa Citation Majdi, N., de Necker, L., Fourie, H., Loggenberg, A., Netherlands, E.C., Bunte-Tschikin, J. et al., 2022, ‘Diversity and distribution of benthic invertebrates dwelling rivers of the Kruger National Park, South Africa’, Koedoe 64(1), a1702. https://doi.org/10.4102/koedoe.v64i1.1702 Original Research Diversity and distribution of benthic invertebrates dwelling rivers of the Kruger National Park, South Africa Nabil Majdi, Lizaan de Necker, Hendrika Fourie, Ané Loggenberg, Edward C. Netherlands, Jasmine Bunte-Tschikin, Walter Traunspurger, Gerhard C. du Preez Received: 04 Dec. 2021; Accepted: 29 Apr. 2022; Published: 27 June 2022 Copyright: © 2022. The Author(s). Licensee: AOSIS. This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Meiobenthos (or meiofauna) are microscopic invertebrates that inhabit biofilms and interstitial spaces in rivers. They are diverse and extremely abundant, and they perform essential ecological functions by linking microbial production to higher trophic levels (e.g. macrobenthic invertebrates and fishes). However, meiobenthic communities remain poorly studied in Africa. Here, we sampled meio- and macrobenthic invertebrate communities associated with biofilms and sediments across an upstream–downstream gradient along the Olifants, Sabie and Crocodile rivers flowing through the Kruger National Park (KNP). We expected to link differences in community structure to environmental gradients as those rivers show different degrees of anthropogenic stress as they enter the park. Both meio- and macrobenthic communities differed across rivers and also structured along an upstream–downstream gradient. The upstream sites, which were the closest to the park borders, consistently showed a lower diversity in all three rivers. There, the invasive snail Tarebia granifera strongly dominated (making up 73% – 87% of the macrobenthos), crowding hard substrates, while concomitantly the abundances of biofilm-dwelling meiobenthos like nematodes and rotifers were substantially reduced. Nevertheless, the diversity and evenness of communities then tended to increase as water flowed downstream through the park, suggesting a beneficial effect of protected river reaches on benthic invertebrate diversity. However, for the Crocodile River, which makes up the southern border of the park, this trend was less conspicuous, suggesting that this river may experience the greatest threats. More generally, benthic invertebrate communities were driven by the concentrations of phosphates, sulphates, ammonium and organic matter and by substrate characteristics. Conservation implications: Meiobenthic organisms are very abundant in KNP rivers and react to environmental gradients; thus, they should be more considered for bio-monitoring or conservation of comprehensive assemblages of animals. Interestingly, protected reaches tended to show a reduced dominance of the invasive T. granifera and a higher diversity of benthic invertebrates. Keywords: biodiversity; aquatic ecology; community structure; meiofauna; macrofauna; invasive species; pollution. Introduction Rivers represent only a minor part of Earth’s water resources, but they play a key role for the connectivity of ecosystems and the maintenance of biological diversity at a global scale (Dudgeon et al. 2006). However, most river ecosystems have been degraded by human activities throughout history (e.g. fisheries, commercial routes, irrigation, production of energy, receptacles for wastewater and introduction of invasive species), resulting in a current situation where rivers are one of the most threatened ecosystems worldwide (Malmqvist & Rundle 2002). As an example, South African rivers are intensively exploited to ensure water security for urban areas, agriculture and industry, where water use often exceeds natural water availability (O’Keeffe 1989). Those rivers further serve as outlets for urban wastewater, excess nutrients and xenobiotics from agriculture, as well as metals draining from mine exploitation or industrial activities (e.g. Jackson et al. 2013; Rimayi et al. 2018; Riddell et al. 2019). Nel et al. (2007) analysed the health status of South African rivers and found that 84% of river reaches (and especially those located in the largest, permanently flowing river systems) were degraded to the point they could be considered as highly threatened. Nature conservation regions including relatively large portions of river catchments can help mitigate pollution inputs and restore water quality (Nel et al. 2007). Although some xenobiotics like organochlorine pesticides can still drift over long distances, their bio-accumulation potential in aquatic organisms of conservation areas is correlated with their proximity to pollution sources (Wolmarans et al. 2021). This means that river reaches further downstream from pollution sources should show a better ecological status. To reduce threats and develop proactive management practices, one requires an in-depth understanding of the response of natural communities and ecosystem processes to water quality. However, in the case of South African rivers, management of water resources is often decided in a context of incomplete knowledge of the response of the different biological communities present in the field, leading to uncertainty in decision-making (see, e.g., Roux, Kleynhans & Thirion 1999a; Roux et al. 1999b). For example, management decisions concerning the conservation actions applicable in the large protected area of the Kruger National Park (KNP) have been traditionally based on the assessment of abiotic rather than biotic factors (Solomon et al. 1999). This is based on the assumption that resultant biotic patterns are likely to be correlated with abiotic components. However, to measure, protect and restore natural river integrity in the KNP, it would be more relevant to include multiple biological indicators of water quality (Rogers & Biggs 1999). Pollutants not only alter the chemistry of the water column but also persist in sediments (Gerber et al. 2015) and within tissues of organisms (Gerber et al. 2016; Seymore 2014; Wolmarans et al. 2021). The consequences of altered environmental conditions have been found to translate into conspicuous modifications of the structure of bacteria, diatom, fish and insect communities in KNP water (Farrell et al. 2019; Rasifudi et al. 2018; Riddell et al. 2019; Shikwambana et al. 2021). Another cause of stress to the native aquatic communities of KNP is the accidental introductions followed by problematic blooms of invasive pathogenic micro-organisms, fishes, crayfishes and molluscs species (e.g. Jones et al. 2017; Macdonald 1988; Petersen et al. 2017). While the response of vertebrates to abiotic and biotic threats is conspicuous, the response of benthic invertebrates (living hidden in biofilm matrices and aquatic sediments) is mostly overlooked when it comes to assess ecosystem health status. Among those benthic invertebrates, one may further make a distinction between the macrobenthos, such as aquatic insects, snails and leeches that are large enough to be visible and usually are retained on 500 µm meshes, and the meiobenthos, such as nematodes, rotifers, copepods and tardigrades that are so tiny that they are invisible and usually pass through 500 µm-meshes but are retained on 20 µm meshes (Ptatscheck, Gehner & Traunspurger 2020a). Meiobenthic invertebrates are little studied although they show complex behaviours and extraordinary physiologies that allow them to colonise most, if not all, benthic habitats (Brüchner-Hüttemann, Ptatscheck & Traunspurger 2020; Rebecchi, Boschetti & Nelson 2020) where they can reach remarkable abundances (between 105–106 ind. m−2 on any submerged substrate) (Majdi, Schmid-Araya & Traunspurger 2020; Traunspurger, Wilden & Majdi 2020). Meiobenthos play an important role in riverine food webs as intermediaries: grazing on microbes (e.g. bacteria, protozoans and micro-algae) and serving as prey for a variety of macro-invertebrates and fish juveniles (Majdi & Traunspurger 2015; Ptatscheck et al. 2020b; Schmid-Araya et al. 2002). Some studies have highlighted the potential of meiobenthic invertebrates to reflect aquatic ecosystem health in South Africa: for instance, Gyedu-Abadio et al. (1999) observed that the density and diversity of nematode assemblages were affected by the concentration of metals in a subtidal portion of the Swartkops estuary. Further experimentation in the laboratory showed that nematode genera like Axonolaimus, Theristus and Paramonhystera were tolerant to metal pollution (Gyedu-Ababio & Baird 2006). In Europe, a nematode ‘species at risk’ index has been developed to monitor river sediment pollution (Höss et al. 2011). This index adds value to routinely used macrobenthic-based indices because nematodes are ubiquitous and experience the effects of pollutants during their whole-life cycle within the sediment (Brüchner-Hüttemann et al. 2021). Whether such meiobenthos-based indices could be tested, developed and added to the bio-indication toolbox outside Europe is a question of general interest for ecologists and stakeholders, but lacks support from an insufficient collection of field data so far. The specific objective of this study is to provide fresh insights into the diversity and distribution of benthic invertebrates (meio- and macrobenthos) in three large rivers (Olifants, Sabie and Crocodile) flowing through an emblematic protected area of austral Africa, the KNP. We sampled the main benthic habitats in those rivers (i.e. biofilms associated with hard substrates and surface sediment) to examine the structure of resident communities. The three rivers were expected to show a gradient of environmental stress as they enter the park: Olifants and Crocodile rivers experience different types of pressures (mostly related to mining practices for Olifants and to agricultural practices for Crocodile), while Sabie is quite pristine in comparison. Thus, we expected that the structure of benthic invertebrate communities could reflect gradients of pollution experienced by the three rivers. Finally, an upstream–downstream gradient was also expected to be manifest through a beneficial effect of preserved areas of the park on the taxonomic richness of benthic assemblages. Methods Study sites The landscape in the southern part of the KNP consists of plains showing a gentle eastward slope and being drained by three major river systems (Crocodile, Sabie and Olifants rivers; Figure 1). The Crocodile (catchment area: 10 455 km2) and Sabie rivers (6252 km2) start off in the Drakensberg region. The Olifants River has the largest catchment area (54 434 km2) and starts in the Highveld region (Muller & Villet 2004; Venter & Bristow 1986). It generally shows the highest annual run-off followed by the Crocodile and Sabie rivers (Muller & Villet 2004) but, at the time of sampling, all rivers had similar, relatively low discharges (3.82 m3 s−1, 3.87 m3 s−1 and 4.11 m3 s−1 in Olifants, Crocodile and Sabie rivers, respectively). This was the result of a steady decrease in flow over 1 month prior to sampling (Figure 1-A1). FIGURE 1: Map of the southern part of the Kruger National Park showing residential and agricultural areas, the different ecozones and river systems with the location of sampling sites. The three rivers drain residential and agricultural areas and then flow eastwards through well-preserved ecozones in the park, except for the Crocodile River that delineates the southern border of the park, one of its bank draining the large agricultural zone of the Lowveld (Figure 1). The Crocodile River is considered highly threatened, draining untreated industrial and household wastewater from large cities such as Johannesburg and Tshwane, its water being increasingly used for irrigation and receiving agricultural effluents. The Olifants River drains the industrial fertiliser complex of Phalaborwa as well as phosphate rock mining facilities (foskorite and pyroxenite) before entering the park and is considered one of the most polluted rivers in South Africa (Gerber et al. 2015 and references cited therein). In contrast, the Sabie River is recognised as a remarkable hotspot of aquatic biodiversity and one of the most pristine river systems of South Africa (Riddell et al. 2019). Three sites were selected in each river system, with the help of South African National Park officials and KNP game rangers. Each site was chosen according to its accessibility, the availability of riverine microhabitats (hard and soft substrates) and its coherence regarding our aim to examine the effects of an environmental gradient through the park. As far as possible, we selected river reaches that showed the most similar hydro-morphologies and without riparian canopy cover (e.g. Figure 2-A1). Upstream sites were located close (< 500 m) to the park’s border, and to ensure substantial spatial coverage, each sampling site was located 10 km – 60 km apart from each other (Figure 1). Sampling Samples were collected over a 4-day sampling campaign in May 2019 during austral autumn, a period of low flow conditions (Figure 1-A1). To homogenise sampling, at each site, a 25-m reach was selected along which five 1 m2 shallow (< 30 cm water depth) plots located approximately 5 m apart were set (e.g. Figure 3-A1). Each plot comprised both hard substrates (in the form of large bedrock or cobble stones) and soft substrates (mostly in the form of fine, sandy sediment). To avoid trampling disturbance, plots were always sampled from downstream to upstream (Figure 3-A1). Samples were placed at 4°C in the dark upon collection and further preserved with fixatives or frozen until further laboratory analyses could take place. At each plot, 500 mL of water was collected for nutrient analyses; and water temperature, pH, conductivity, dissolved inorganic salts, total dissolved solids and turbidity were measured in situ with a portable probe (HI9813-5, Hanna Instruments Inc., Bedfordview, South Africa). Sediments were sampled to measure total organic carbon by pushing one 2-cm diameter Perspex corer (Figure 1-A1) through 5 cm of the upper sediment layer at one location chosen randomly within each plot. Sediments were sampled for meiobenthos by pushing two 2-cm diameter Perspex corers through 5 cm of the upper sediment layer at two different locations within each plot. The two sediment cores collected were merged in the same tube and the resulting sample was fixed in a final concentration of 4% buffered formaldehyde. A calibrated brush-sampler (Figure 3-A1; Peters et al. 2005) was used for underwater sampling of 3.14 cm2 areas of biofilm (and its associated meiobenthos) growing on hard substrates such as rocks and cobbles. This process was repeated three times at different random locations within each plot, using a 20-µm-aperture nylon mesh sieve each time to concentrate the collected subsamples. The contents of the sieve were then poured into a tube and the resulting total sample (representing an area of 9.42 cm2) was fixed in 4% buffered formaldehyde. After meiobenthos sampling, a ‘kick-sampling’ procedure was used to obtain a semi-quantitative sample representative of the macrobenthic community dwelling in each plot. Stones, macrophytes and superficial sediments were thoroughly agitated by hand and foot all over the plot for 3 min. The resulting suspended organisms were collected in a sturdy hand net (500 µm meshes) held downstream of the plot. The content of the net was poured in a jar and preserved in 70% ethanol. Sample processing Water chemistry variables were measured in the laboratory using a Merck Spectroquant Pharo 300 Spectrophotometer and appropriate test kits following standard protocols (Merck KGaA 2014). Variables were measured on defrosted, unfiltered water samples and included Cl- (further referred to as chloride) (Merck protocol #14897), P-PO43-(phosphate) (#14848), SO42- (sulphate) (#14791), N-NH4+ (ammonium) (#14752), N-NO3- (nitrate) (#14773). Subsamples of frozen sediments were defrosted, dried in an oven at 60°C for 96 h, weighted, burned to ashes at 600°C for 6 h, then weighted once again to measure the percentage of total organic content (TOC) as a ratio between ash-free dry weight and dry weight of the sediment. To extract meiobenthos from sediment samples, a density-centrifugation procedure was used in the laboratory involving the flotation of organic particles on Ludox HS-40 (specific gravity set at 1.14), following the method of Pfannkuche and Thiel (1988). The supernatant, containing the meiobenthic invertebrates, was poured through 20-µm meshes, preserved in 4% formaldehyde and stained with a few drops of Rose Bengal. The pellet (i.e. inorganic sediment particles from which the meiobenthos had been extracted) was further passed through a series of stacked sieves (1000 µm, 500 µm, 250 µm, 125 µm, 63 µm and 32-µm-aperture meshes). Each sediment size fraction retained on the respective sieve (further coined F1000, F500, F250, F125, F63 and F32) was dried in an oven at 60°C for 96 h and weighted to estimate the sediment particle size distribution of each sample. Meiobenthic invertebrates were counted and identified under a 40–80x magnification stereomicroscope (Nikon SMZ1500). Morpho-taxonomic features were used to classify the meiobenthic individuals into coarse taxonomic groups: nematodes, rotifers, gastrotrichs, copepods and their Nauplii larvae, ostracods, oligochaetes, tardigrades, mites, molluscs (mostly veliger and pediveliger larval stages) and larval stages of plecopters, ephemeropters and dipters (detailing Chironomidae and Ceratopogonidae families). The weight of sediment from which meiobenthos had been extracted was used to express abundances per gram sediment dry weight. Meiobenthic organisms dwelling on hard substrates were not extracted using the Ludox procedure and they were directly counted and identified as described above, except that abundances were expressed per area of biofilm. Individuals of macrobenthic invertebrates were counted in each sample under a 1–5x magnification dissection microscope (Nikon C-LEDS) and identified to lowest possible taxonomic level based on morphological characteristics and regional identification key for aquatic macro-invertebrates. Total abundances were expressed as number of individuals per 1 m2 plot. However, relative abundances were used for community structure analyses to avoid potential misinterpretations because of the semi-quantitative nature of the ‘kick-sampling’ procedure used to sample (Table 1-A1). Data analysis Statistical analyses were performed in R (v 4.0.3.; R Development Core Team 2020) using packages ‘vegan’ and ‘indicspecies’ (De Cáceres & Legendre 2009). Effects of river identity (Olifants, Crocodile or Sabie) and reach location (upstream, midstream or downstream) on abiotic parameters, faunal abundances and macrobenthic diversity indices were tested using a two-way analysis of variance (ANOVA) method followed by a post-hoc Tukey’s HSD test for pairwise comparisons. The univariate data were checked for normality and homoscedasticity using Shapiro–Wilk test and Levene’s test, respectively. If data were not normally distributed and homoscedastic, the non-parametric Kruskal–Wallis rank sum test followed by a multiple comparison test was used instead. For macrobenthos data, effects of river identity and reach were tested on Shannon’s diversity index (H’) and Pielou’s evenness (J’). Multivariate effects of river and reach on the benthic community structure were tested using permutational multivariate analysis of variance using distance matrices (PERMANOVA: adonis function with 9999 permutations). The assumption of multivariate homogeneity of variances was first checked using the PERMDISP2 procedure (betadisper function), which is a multivariate analogous to Levene’s test. The function rankindex was also used before tests to rank the performance of different dissimilarity indices. In our case, Bray–Curtis distance was the most relevant dissimilarity index. Non-metric multidimensional scaling (NMDS) based on Bray–Curtis distances (metaMDS function) was used to ordinate samples and species scores. To highlight the most important environmental factors describing the variability in community structure of macrobenthos and meiobenthos, envfit function (that fits vectors of continuous environmental variables) was used. The significance of environmental variables (vectors) fitted onto the ordination was assessed using a goodness-of-fit (R2) statistics and empirical p-values based on a permutation test (N = 999 permutations). The direction of vector showed the direction of the gradient, and the length of the arrow was proportional to the correlation between the variable score and the ordination space. Only significant vectors were further displayed on NMDS plot using ‘p.max=0.05’ as a filter argument when plotting vectors on the ordination space. The response of individual taxon abundances to river identity and reach location was examined using multi-level pattern analysis after De Cáceres and Legendre (2009), which is an extended method of indicator species analysis. Multi-level pattern analysis (multipatt function) is a permutation test examining which taxa are significantly responsible for differences among a group of samples. Multi-level pattern analysis was chosen over other species indicator methods because it is less sensitive to the weight of over-dominant species, and it is not based on the relative contribution to differences (De Cáceres & Legendre 2009). Ethical considerations This article followed all ethical standards for research without direct contact with human or animal subjects. Results Water and sediment environment Overall, water pH was slightly alkaline but did not differ significantly across rivers or reaches (Table 1). However, conductivity was significantly higher in the Olifants River (on average 663 µS cm−1) in comparison to Crocodile (450 µS cm−1) and Sabie (104 µS cm−1) rivers (Kruskal–Wallis test: χ2 = 38.7, p < 0.001). The concentrations of sulphate and chloride were likely the dominant contributors to variance in conductivity values as strong positive correlations were observed between conductivity and sulphate and chloride concentrations (Figure 4-A1). Turbidity was significantly lower in Sabie in comparison to Crocodile and Olifants rivers (ANOVA, F2,36 = 31, p < 0.001), and in the latter two more turbid rivers, turbidity was found to be significantly higher at upstream sites (ANOVA, F2,36 = 9.8, p < 0.001). Nitrate and ammonium concentrations did not show significant associations with river or with reach, although they were found to be positively correlated with phosphate concentration (Figure 4-A1). Phosphate was significantly more concentrated in the Crocodile River (ANOVA, F2,36 = 37.1, p < 0.001) and interestingly its concentration tended to increase in downstream areas in all rivers (ANOVA, F2,36 = 187.7, p < 0.001). TABLE 1: Water and sediment environmental parameters measured at three sites along three rivers flowing through the Kruger National Park. The sediment of the Olifants River showed a significantly greater proportion of organic matter in comparison to the other rivers (Kruskal–Wallis test: χ2 = 8.2, p = 0.015). In terms of granulometry, relatively coarse sediment dominated (fractions > 500 µm represented on average 63.5% of the sediment grain size distribution). However, it is worth noting that the downstream site on Olifants River showed unusually low proportions of the largest grain size categories (Table 1). The proportion of coarser sediment was negatively correlated with the proportion of finer grain size categories (Figure 4-A1). The coarsest sediment fraction (> 1000 µm) showed a significant river*reach interaction (ANOVA, F4,36 = 6.2, p < 0.001), essentially confirming that Sabie’s riverbed showed coarser sediment as we moved downstream. Abundance and diversity of benthic invertebrates Macrobenthos was at least three orders of magnitude less abundant than the meiobenthos: on average, 130 macrobenthic individuals per m2 against 332 000 ind. m−2 for biofilm-dwelling meiobenthos (Table 2). Macrobenthos abundance appeared significantly lower in the Sabie River (ANOVA, F2,36 = 8.4, p < 0.001). The macrobenthic community consisted of 82 taxa from 53 families or coarser taxonomic groups. Fifty-seven taxa could be further resolved to genus or species level (Table 1-A1). The taxon dominating the macrobenthic community was Tarebia granifera, making up 36.95% of the macrobenthic community on average. However, T. granifera made up 73% – 87% of individuals in upstream reaches (Figure 2a). As a result, reach position had a significant effect on the relative abundance of T. granifera (ANOVA, F2,36 = 101, p < 0.001). The larval stage of the ephemeropteran genus Caenis sp. was the second most common macrobenthic taxon, representing 16.8% of the assemblage on average and being especially abundant in mid- and downstream reaches of the Sabie and Crocodile rivers (Figure 2a, Table 1-A1). Other macrobenthic taxa collectively represented 46.2% of the assemblage on average (Figure 2a), with only 14 taxa making up individually more than 1% of the assemblage (Tricorythus sp., Corbicula fluminalis, Paragomphus sp., Acanthiops sp., Leptophlebiidae, Cheumatopsyche sp., Notogomphus sp., Hydroglyphus sp., Laccocoris sp., Pseudagrion sp., Simulium sp., Chironominae Anisops sp., Elminae larvae; see Table 1-A1). On average, 9.2 different macrobenthic taxa were identified from each sample (Table 2); however, taxon richness was significantly lower in upstream reaches in comparison to mid- or downstream reaches (Table 2, ANOVA, F2,36 = 12.7, p < 0.001). Shannon’s diversity index (H’) was also lower in upstream reaches (Table 2) but was further found to be significantly lower in the Crocodile River as well (ANOVA, F2,36 = 11, p < 0.001). Pielou’s equitability index followed a similar pattern as H’ with communities generally more equitable in mid- and downstream reaches (ANOVA, F2,36 > 11.2, p < 0.001), except for the downstream reach of the Crocodile River, which was uneven and strongly dominated (87%) by Caenis sp. (Table 2, Figure 2a). FIGURE 2: Composition of the benthic invertebrate community in three rivers flowing through the Kruger National Park, South Africa. (a) Mean (N = 5) relative composition of the macrobenthic community at each sampling site. Only the first 51 most abundant taxa out of 82 have assigned coloured legend. Taxa abbreviations are detailed in Table 1-A1. (b) Mean (N = 5) abundance and composition of the meiobenthic community in the sediment. (c) Mean (N = 5) abundance and composition of the meiobenthic community in biofilms. Examples of sampling site. TABLE 2: Abundance, diversity and taxonomic composition of benthic invertebrate communities in three rivers flowing through the Kruger National Park (detailed taxonomic composition for macrobenthic invertebrates is presented in Table 1-A1). Meiobenthos abundance in sediment did not show any significant trend (Table 2). However, biofilm-dwelling meiobenthic organisms were significantly less abundant in upstream reaches (ANOVA, F2,36 = 9.3, p < 0.001). Nematodes and rotifers were the numerically dominant representatives of the meiobenthos in all rivers and reaches (Figure 2b and Figure 2c), but while nematodes made up on average 21% of the meiobenthic assemblage both in the sediment and in biofilms, rotifers were relatively more abundant in biofilms (45.7% of individuals) than in sediments (31.9%; Table 2). Structural responses of the macrobenthic community to environmental gradients The structure of the macrobenthic community was significantly affected by river identity (PERMANOVA, F2,42 = 29, R2 = 0.12, p = 0.002). Multivariate homoscedasticity was met without data transformation (PERMDISP2, p = 0.72); however, the effects of reach location on macrobenthic community structure could not be tested as multivariate homoscedasticity of variances was not met, even after reasonable data transformation techniques. This was caused by a largely heterogeneous distribution of distances to centroids among reaches: samples from upstream reaches had a very skewed community structure because of the large dominance of T. granifera (Figure 2a), resulting in a very low dispersion of sample scores around group centroids (N = 15, average distance to centroid: 0.15 for upstream samples), while for mid- and downstream sections, average distances of samples to centroids were 0.51 and 0.54, respectively. Phosphate was found to be the first driver of macrobenthic community structure (R2 = 0.29, p < 0.001), followed by variables associated with sediment granulometry: F125 (R2 = 0.24, p = 0.002), F1000 (R2 = 0.21, p = 0.009), F250 (R2 = 0.20, p = 0.01) and F63 (R2 = 0.15, p = 0.034). Sulphate was also a significant driver (R2 = 0.14, p = 0.036) although its direction was opposed to that of phosphate (Figure 3). FIGURE 3: Non-metric multidimensional scaling ordination of the macrobenthic community. Ellipses show the 50% confidence regions for the locations of group of samples’ centroids. ‘Spider webs’ show distances of samples to group centroids they belong. Species scores are directly ordinated on the biplot as text abbreviations (in Table 1-A1). Stars highlight taxa significantly associated with one or two river systems after multi-level pattern analysis. Environmental fitting was used to determine significant effects of environmental factors: the direction of vector showed the direction of the gradient, and the length of the arrow was proportional to correlation strength. Only significant environmental variables were displayed. Other abbreviations include NMDS, non-metric multidimensional scaling; PO4, concentration of phosphate in the water; SO4, concentration of sulphate in the water; F1000, fraction of sediment particles > 1000 µm; F250, fraction of sediment particles of size 250 µm – 500 µm; F125, fraction of sediment particles of size 125 µm – 250 µm; F63, fraction of sediment particles of size 63 µm – 125 µm. The multi-level pattern analysis identified 14 taxa significantly associated with one river group. Those groups are highlighted with stars on the NMDS biplot (Figure 3). Namely, only macrobenthic nematodes were found to be significantly associated with the Crocodile River. The ephemeropteran larvae Tricorythus sp. and Elassoneuria sp., the trichopteran larvae Cheumatopsyche sp. and Macrostemum sp., the dipteran larvae of Simulium sp., the dragonfly larvae of Pseudagrion sp. and coleopteran adults of Berosus sp. were significantly associated with the Olifants River. The ephemeropteran larvae of Machadorythus maculatus, the hemipteran Laccocoris sp., the dipteran larvae of subfamily Chironominae, the dragonfly larvae of Lestinogomphus sp. and Notogomphus sp., and the coleopteran larvae of subfamily Larainae were significantly associated with the Sabie River. The clam Corbicula fluminalis was associated with both Crocodile and Olifants rivers. Dragonfly larvae Coenagrionidae and Paragomphus sp., and ephemeropteran larvae of the Leptophlebiidae family were associated with both Crocodile and Sabie rivers (Figure 3). Structural responses of the meiobenthic community to environmental gradients The structure of the meiobenthic community dwelling sediments was significantly affected by river identity (PERMANOVA, F2,36 = 2.3, R2 = 0.08, p = 0.003), multivariate homoscedasticity being met without data transformation (PERMDISP2river, p = 0.21; PERMDISP2reach, p = 0.4). It was significantly affected by the interaction river*reach (PERMANOVA, F4,36 = 2.3, R2 = 0.17, p < 0.001). This was caused by the high nematode and rotifer abundances found in the downstream reach of Olifants River (Figure 2b). Sediment TOC was found to be the first driver of sediment-dwelling meiobenthic community structure (R2 = 0.38, p < 0.001; Figure 4), followed by variables associated with fine sediment granulometry: F63 (R2 = 0.34, p < 0.001), F32 (R2 = 0.28, p < 0.001). Water quality variables such as total dissolved solids (TDS) and conductivity (Cond) were also significant drivers (R2 = 0.18, p = 0.012, and R2 = 0.17, p = 0.017, respectively). The multi-level pattern analysis identified only one taxon: nematodes, as being significantly associated with the Olifants and Sabie rivers (Figures 2b, Figure 4 and Table 2). When performed to examine the associations with a specific river reach, the multi-level pattern analysis identified two taxa, namely, the larvae of Ceratopogonidae and copepods being associated with both mid- and downstream reaches (Table 2). FIGURE 4: Non-metric multidimensional scaling ordination sediment-dwelling meiobenthic community. Ordination is based on untransformed absolute abundances (ind. gSedDM−1). The structure of the meiobenthic community dwelling biofilms growing on hard substrates was significantly affected by river identity (PERMANOVA, F2,36 = 4.9, R2 = 0.15, p < 0.001), river reach (PERMANOVA, F2,36 = 5.1, R2 = 0.16, p < 0.001) and by the river*reach interaction (PERMANOVA, F4,36 = 1.8, R2 = 0.11, p = 0.023). Multivariate homoscedasticity was achieved without data transformation (PERMDISP2river, p = 0.79; PERMDISP2reach, p = 0.65). This may have been caused by the relatively coherent pattern of lower abundances in upstream reaches coupled with relatively higher abundances of chironomids in the Sabie River (Figure 2b and Table 2). Concentration of ammonium in the water was found to be the strongest driver of biofilm-dwelling meiobenthic community structure (R2 = 0.17, p =0.02; Figure 5), followed by TDS and conductivity (R2 = 0.145, p = 0.035 and R2 = 0.146, p = 0.032, respectively). The multi-level pattern analysis identified three taxa significantly associated with the Sabie River, namely, oligochaetes, tardigrades and gastrotrichs; while one taxon (Chironomids) was significantly associated with both Olifants and Sabie rivers (Figure 5). When performed to examine associations with a specific river reach, the multi-level pattern analysis identified four taxa being associated with both mid- and downstream reaches, namely, ostracods and larval stages of Ceratopogonidae, Chironomidae and Ephemeropterans. FIGURE 5: Non-metric multidimensional scaling ordination of the biofilm-dwelling meiobenthic community. Ordination is based on untransformed abundances (ind. cm−2). Discussion Here we studied for the first time the entire assemblage of benthic invertebrates dwelling hard and soft substrates of three major rivers flowing through the KNP, South Africa. Overall, we observed a great diversity of benthic invertebrates of various sizes and morphologies belonging to a total of seven major animal groups: Arthropoda, Mollusca, Nematoda, Rotifera, Gastrotricha, Tardigrada and Annelida. We found the highest diversity of families in the Sabie River (42 families), followed by the Olifants (36 families) and Crocodile River (25 families). These results generally confirmed a better biodiversity status for the Sabie River as was already observed in previous studies focusing on macrobenthic invertebrates dwelling rivers of the KNP. For example, Muller and Villet (2004) identified 58, 51 and 46 families in the Sabie, Crocodile and Olifants rivers, respectively. The Crocodile and Olifants rivers are known to experience substantial anthropogenic pressures in comparison to the Sabie River (Riddell et al. 2019). We found different structures for both meio- and macrobenthic communities across the rivers, confirming our hypothesis that benthic communities would react to environmental gradients across rivers. Crocodile River’s water had more phosphates and showed the largest alterations of benthic communities, suggesting that this system faced acute stress. Olifants River’s water showed more sulphates and a higher conductivity and turbidity, and the sediment texture was finer housing higher amounts of organic material. But as Olifants ran through the KNP, diversity and evenness tended to improve. Sabie showed lower concentrations of nutrients and the most diverse assemblages. However, we found that upstream sites located at the park border were consistently less diverse, supporting that the richness of benthic assemblages benefits from the course of the river through protected areas of the park. The Crocodile River, making the southern border of the park, did not show such improvement presumably because it receives effluents from intensive agriculture plots all along its southern catchment outside the KNP (Van der Laan, Van Antwerpen & Bristow 2012). One of the most conspicuous effects of the upstream-downstream gradient was the alteration of community structure because of the strong dominance of T. granifera at all upstream sites. This snail species originates from Asia and was first reported from South Africa in 1999 in northern KwaZulu-Natal (Appleton & Nadasan 2002). It has further rapidly spread into the Mpumalanga province, KNP and Swaziland and will probably continue its expansion to northern sub-tropical lowlands in Zimbabwe and Mozambique (Appleton, Forbes & Demetriades 2009). Fifteen years ago, Wolmarans and De kock (2006) sampled T. granifera from the Sabie and Crocodile rivers, but did not find it in Olifants. In our study we found substantial numbers of T. granifera in all river sites, including all river sites along the Olifants. T. granifera can attain high densities in invaded areas, thus impacting other benthic fauna and flora (Appleton et al. 2009; Jones et al. 2017). In invaded estuarine reaches, it has been shown that T. granifera mostly exploits algal biofilms and detritus (Miranda, Perissinotto & Appleton 2011). Thus, their grazing activity on hard substrates might particularly affect resource availability for biofilm-dwelling meiobenthos, as was evidenced empirically in the case of other aquatic snails dampening biofilm-associated assemblages (Peters, Hillebrand & Traunspurger 2007). Here, we found a significant reduction of biofilm-dwelling meiobenthos in upstream reaches, reduction that could have been caused by a more severe biofilm grazing activity by T. granifera (we observed that ostracods and small instars of insect larvae were rare in the locations strongly dominated by T. granifera). However, stronger evidences are needed to show direct causality so we recommend experimental evaluation of the direct top-down or indirect bottom-up effect of the invasive T. granifera on local biofilm-dwelling meiobenthos. Interestingly, our results showed that T. granifera did not over-dominate other downstream reaches in the KNP. A possible rationale is that the more diverse communities in mid- and downstream reaches could have reduced the invasion success of T. granifera through increased competition, predation, or parasitism (Miranda et al. 2016). This would also need further experimental evaluation. Nevertheless, our results contribute to support that the success of invasive species is all the more increased when indigenous fauna is impoverished because of stress like habitat deterioration (Van der Velde et al. 2002). The positive association of sulphates with the Olifants River is not surprising as this river has historically suffered from salt enrichment because of sulphates that enter the river from effluent discharge which is one of the main contributors to the high salinity measurements in this river (Goetsch & Palmer 1997; Riddell et al. 2019). Furthermore, we found a positive association of phosphates in the Crocodile River as this system receives a great deal of organic effluents through fertilisers from agricultural activities taking place along its catchment outside the KNP (Van der Laan et al. 2012). Similar results were obtained by Riddell et al. (2019) who attributed impacts of orthophosphates on rivers in KNP to agricultural activities taking place immediately adjacent to the park. In our study, increased sulphate and phosphate levels, along with other organic and inorganic variables emerged as significant drivers of the meio- and macrobenthic community structure. Biota such as nematodes, the water scavenger beetles Berosus spp., the Hydropsychid caddisfly Cheumatopsyche spp., Caenis spp. mayflies, dipteran Simulidae and Chironominae are adapted to tolerate a wide range of water quality parameters, such as organic enrichment and high salinity (Erasmus et al. 2021; Malherbe, Wepener & Van Vuren 2010), which is likely why they were associated with the Olifants and Crocodile rivers. Habitat likely also played a role in the distribution of macrobenthos in the rivers as most of the biota specifically associated with either the Olifants or the Sabie river also have particular habitat and water flow preferences. For instance, the mayflies Elassoneuria sp. and Tricorythus sp. and the caddisflies Macrostemum sp. were significantly associated with the Olifants River presumably because of their habitat preference for rocky surfaces in fast-flowing water (Gerber & Gabriel 2002), a habitat that was more represented at sampling sites in this river. Similarly, the Sabie River habitats consisted of dense reed vegetation (Figure 2-A1), and most of the macrobenthos that were significantly associated with the Sabie River have affinities with such dense vegetation patches, particularly Laccocoris, or show affinities with the sandy to muddy substrate in slower flowing water at the edges of streams, particularly the gomphids Lestinogomphus sp. and Notogomphus sp. (Gerber & Gabriel 2002). Although they are omnipresent in limnetic habitats, meiobenthic communities are poorly considered in ecology and conservation studies (Majdi et al. 2020). This knowledge gap is even more preoccupying for tropical and subtropical regions because studies about freshwater meiobenthos ecology and taxonomy have traditionally focused on temperate regions from the northern hemisphere (e.g. Fontaneto et al. 2012; Traunspurger et al. 2020; Zullini 2014). Here we found an abundant meiobenthic community dominated by nematodes and rotifers, and the taxonomic composition and abundance values reported in the present study are quite comparable to meiobenthic communities found in streams and rivers in Europe (e.g. Brüchner-Hüttemann et al. 2020; Majdi et al. 2012; Majdi, Threis & Traunspurger 2017). Some meiobenthic groups have been previously investigated in the KNP: for example, 30 years ago, Botha and Heyns (1991, 1992, 1993) identified a total of 33 nematode species (including species new to science) in sediment samples collected from Crocodile and Olifants rivers. Here we also observed an abundant nematode community, particularly in fine sediments of the Olifants River. In South African estuaries, Pillay and Perissinotto (2009) and Nozais, Perissinotto and Tita (2005) have observed that the abundance and community structure of the meiobenthos respond to river inflow and drought conditions (mostly through a decrease in physiologically sensitive taxa under harsher conditions). Here we found that biofilm- and sediment-dwelling meiobenthos exhibited different responses: nutrient and water quality proxies affected the biofilm-dwelling organisms, while sediment-dwelling meiofauna was also shaped by sediment texture and the availability of benthic organic material. Overall, the observed patterns are in agreement with previous studies showing that freshwater meiobenthic communities can become relevant bio-indicators of sediment pollution, water eutrophication or sediment texture (e.g. Haegerbaeumer et al. 2017; Schenk et al. 2020; Traunspurger et al. 2020). However, further studies should detail the response of dominant meiobenthic phyla (e.g. nematodes or rotifers) to the genus or species level in order to identify sensitive and tolerant species and develop subtler bio-indication tools. Conclusion The KNP is known worldwide for its unique and emblematic terrestrial megafauna but less so for its remarkable freshwater fauna of benthic invertebrates. In this study we observed that the diversity of those invertebrates increased as we moved away from park borders, highlighting a potential beneficial effect of protected river reaches. Furthermore, in protected reaches, the assemblage was more even, less affected by the domination of the invasive snail T. granifera. The effects of this species on native assemblages should be further assessed, our results suggesting that T. granifera could affect the structure of biofilm-dwelling communities. Other drivers (like the concentration of nutrients in water and sediment granulometry) also showed significant effects on the distribution of benthic invertebrates. We recommend that the interplay between biotic and abiotic drivers should be further studied for a more comprehensive management and conservation of aquatic resources. Finally, our results stress that the Crocodile River showed the poorest and most unbalanced communities all the way. This suggests that further environment-damaging project on the Crocodile River (such as current applications for large-scale coal-mining on the southern boundary of the park) could have catastrophic impacts on an already stressed river system. Acknowledgements The authors are grateful to South African National park game guards and officials for providing kind support during sampling and stay at Phalaborwa Research House and Skukuza Research Camp. They thank Robin Petersen, Mike Silberbauer, Harold van Niekerk, Jacques Venter, Marié du Toit and Stefanie Gehner for their helpful technical support and expertise. Competing interests The authors declare that they have no financial or personal relationships that may have inappropriately influenced them in writing this article. Authors’ contributions All authors contributed to the conception and design of the study. Material preparation, data collection and analysis were performed by N.M., A.L., E.C.N., J.B.-T. and G.C.d.P. The first draft of the manuscript was written by N.M. and all authors commented on previous versions of the manuscript. The final manuscript has been read and approved by all the authors. Funding information The authors declare that no funds, grants or other support were received during the preparation of this manuscript. Data availability The data sets generated during the current study are available from the corresponding author upon reasonable request. Disclaimer The views and opinions expressed in this article are those of the authors and do not necessarily reflect the official policy or position of any affiliated agency of the authors. The research complies with ethical standards and sample collection, and the experimental procedures were approved in strict accordance with institutional guidelines. Sampling in the park was authorised under strict guidance and after permissions were accorded by SANParks officials. Sampling only concerned invertebrates, and reasonable sample sizes were collected. References Appleton, C.C., Forbes, A.T. & Demetriades, N.T., 2009, ‘The occurrence, bionomics and potential impacts of the invasive freshwater snail Tarebia granifera (Lamarck, 1822) (Gastropoda: Thiaridae) in South Africa’, Zoologische Mededelingen 83, 525–536. Appleton, C.C. & Nadasan, D.S., 2002, ‘First report of Tarebia granifera (Lamarck, 1822) (Gastropoda: Thiaridae) from Africa’, Journal of Molluscan Studies 68(4), 399–402. https://doi.org/10.1093/mollus/68.4.399 Botha, A. & Heyns, J., 1991, ‘Dorylaimoidea (Nematoda) from rivers in the Kruger National Park’, KOEDOE 34(2), 1–24. https://doi.org/10.4102/koedoe.v34i2.420 Botha, A. & Heyns, J., 1992, ‘Freshwater nematodes of the genera Thornenema and Mesodorylaimus from the Kruger National Park with a diagnostic species compendium for South African species of the genus Mesodorylaimus (Nematoda: Dorylaimida)’, KOEDOE 35(1), 25–42. https://doi.org/10.4102/koedoe.v35i1.387 Botha, A. & Heyns, J., 1993, ‘Species of the genera Oxydirus, Dorylaimellus (Axodorylaimellus), Laimydorus and Rhabdolaimus from rivers in the Kruger National Park (Nematoda: Dorylaimida and Araeolaimida)’, KOEDOE 36(1), 49–60. https://doi.org/10.4102/koedoe.v36i1.361 Brüchner-Hüttemann, H., Höss, S., Ptatscheck, C., Brinke, M., Schenk, J. & Traunspurger, W., 2021, ‘Added value of the NemaSPEAR [%]-index to routinely used macrofauna-based indices for assessing the quality of freshwater sediments’, Ecological Indicators 121, 107015. https://doi.org/10.1016/j.ecolind.2020.107015 Brüchner-Hüttemann, H., Ptatscheck, C. & Traunspurger, W., 2020, ‘Meiofauna in stream habitats: Temporal dynamics of abundance, biomass and secondary production in different substrate microhabitats in a first-order stream’, Aquatic Ecology 54, 1079–1095. https://doi.org/10.1007/s10452-020-09795-5 De Cáceres, M. & Legendre, P., 2009, ‘Associations between species and groups of sites: Indices and statistical inference’, Ecology 90(12), 3566–3574. https://doi.org/10.1890/08-1823.1 Dudgeon, D., Arthington, A.H., Gessner, M.O., Kawabata, Z.I., Knowler, D.J., Lévèque, C. et al., 2006, ‘Freshwater biodiversity: Importance, threats, status and conservation challenges’, Biological Reviews 81(2), 163–182. https://doi.org/10.1017/S1464793105006950 Erasmus, J.H., Lorenz, A.W., Zimmermann, S., Wepener, V., Sures, B., Smit, N.J. et al., 2021, ‘A diversity and functional approach to evaluate the macroinvertebrate responses to multiple stressors in a small subtropical austral river’, Ecological Indicators 131, 108206. https://doi.org/10.1016/j.ecolind.2021.108206 Farrell, M.J., Govender, D., Hajibabaei, M., Van Der Bank, M. & Davies, T.J., 2019, ‘Bacterial diversity in the waterholes of the Kruger National Park: An eDNA metabarcoding approach’, Genome 62(3), 229–242. https://doi.org/10.1139/gen-2018-0064 Fontaneto, D., Barbosa, A.M., Segers, H. & Pautasso, M., 2012, ‘The ‘rotiferologist’effect and other global correlates of species richness in monogonont rotifers’, Ecography 35, 174–182. https://doi.org/10.1111/j.1600-0587.2011.06850.x Gerber, A. & Gabriel, M.J.M., 2002, Aquatic invertebrates of South African rivers field guide, Institute for Water Quality Studies Department of Water Affairs and Forestry, Pretoria, South Africa. Gerber, R., Smit, N.J., Van Vuren, J.H., Nakayama, S.M., Yohannes, Y.B., Ikenaka, Y. et al., 2015, ‘Application of a Sediment Quality Index for the assessment and monitoring of metals and organochlorines in a premier conservation area’, Environmental Science and Pollution Research 22, 19971–19989. https://doi.org/10.1007/s11356-015-5206-z Gerber, R., Smit, N.J., Van Vuren, J.H., Nakayama, S.M., Yohannes, Y.B., Ikenaka, Y. et al., 2016, ‘Bioaccumulation and human health risk assessment of DDT and other organochlorine pesticides in an apex aquatic predator from a premier conservation area’, Science of the Total Environment 550, 522–533. https://doi.org/10.1016/j.scitotenv.2016.01.129 Goetsch, P.-A. & Palmer, C.G., 1997, ‘Salinity tolerances of selected macroinvertebrates of the Sabie River, Kruger National Park, South Africa’, Archives of Environmental Contamination and Toxicology 32, 32–41. https://doi.org/10.1007/s002449900152 Gyedu-Ababio, T.K. & Baird, D., 2006, ‘Response of meiofauna and nematode communities to increased levels of contaminants in a laboratory microcosm experiment’, Ecotoxicology and Environmental Safety 63(3), 443–450. https://doi.org/10.1016/j.ecoenv.2005.01.010 Gyedu-Ababio, T.K., Furstenberg, J.P., Baird, D. & Vanreusel, A., 1999, ‘Nematodes as indicators of pollution: A case study from the Swartkops River system, South Africa’, Hydrobiologia 397, 155–169. https://doi.org/10.1023/A:1003617825985 Haegerbaeumer, A., Höss, S., Ristau, K., Claus, E., Heininger, P. & Traunspurger, W., 2017, ‘The use of meiofauna in freshwater sediment assessments: Structural and functional responses of meiobenthic communities to metal and organics contamination’, Ecological Indicators 78, 512–525. https://doi.org/10.1016/j.ecolind.2017.03.048 Höss, S., Claus, E., Von der Ohe, P.C., Brinke, M., Güde, H., Heininger, P. et al., 2011, ‘Nematode species at risk – A metric to assess pollution in soft sediments of freshwaters’, Environment International 37(5), 940–949. https://doi.org/10.1016/j.envint.2011.03.013 Jackson, V.A., Paulse, A.N., Odendaal, J.P. & Khan, W., 2013, ‘Identification of point sources of metal pollution in the Berg River, Western Cape, South Africa’, Water, Air, & Soil Pollution 224, 1–8. https://doi.org/10.1007/s11270-013-1477-5 Jones, R.W., Hill, J.M., Coetzee, J.A., Avery, T.S., Weyl, O.L.F. & Hill, M.P., 2017, ‘The abundance of an invasive freshwater snail Tarebia granifera (Lamarck, 1822) in the Nseleni River, South Africa’, African Journal of Aquatic Science 42(1), 75–81. https://doi.org/10.2989/16085914.2017.1298984 Macdonald, I.A.W., 1988, ‘The history, impacts and control of introduced species in the Kruger National Park, South Africa’, Transactions of the Royal Society of South Africa 46(4), 251–276. https://doi.org/10.1080/00359198809520128 Majdi, N., Mialet, B., Boyer, S., Tackx, M., Leflaive, J., Boulêtreau, S. et al., 2012, ‘The relationship between epilithic biofilm stability and its associated meiofauna under two patterns of flood disturbance’, Freshwater Science 31(1), 38–50. https://doi.org/10.1899/11-073.1 Majdi, N., Schmid-Araya, J.M. & Traunspurger, W., 2020, ‘Preface: Patterns and processes of meiofauna in freshwater ecosystems’, Hydrobiologia 847, 2587–2595. https://doi.org/10.1007/s10750-020-04301-2 Majdi, N., Threis, I. & Traunspurger, W., 2017, ‘It’s the little things that count: Meiofaunal density and production in the sediment of two headwater streams’, Limnology and Oceanography 62(1), 151–163. https://doi.org/10.1002/lno.10382 Majdi, N. & Traunspurger, W., 2015, ‘Free-living nematodes in the freshwater food web: A review’, Journal of Nematology 47, 28–44. Malherbe, W., Wepener, V. & Van Vuren, J.H.J., 2010, ‘Anthropogenic spatial and temporal changes in the aquatic macro invertebrate assemblages of the lower Mvoti River, KwaZulu-Natal, South Africa’, African Journal of Aquatic Science 35(1), 13–20. https://doi.org/10.2989/16085914.2010.466638 Malmqvist, B. & Rundle, S., 2002, ‘Threats to the running water ecosystems of the world’, Environmental Conservation 29(2), 134–153. https://doi.org/10.1017/S0376892902000097 Miranda, N.A., Measey, G.J., Peer, N., Raw, J.L., Perissinotto, R. & Appleton, C.C., 2016, ‘Shell crushing resistance of alien and native thiarid gastropods to predatory crabs in South Africa’, Aquatic Invasions 11, 303–311. https://doi.org/10.3391/ai.2016.11.3.08 Miranda, N.A., Perissinotto, R. & Appleton, C.C., 2011, ‘Feeding dynamics of the invasive gastropod Tarebia granifera in coastal and estuarine lakes of northern KwaZulu-Natal, South Africa’, Estuarine, Coastal and Shelf Science 91(3), 442–449. https://doi.org/10.1016/j.ecss.2010.11.007 Muller, W.J. & Villet, M.H., 2004, Similarities and differences between rivers of the Kruger National Park, WRC Report, Water Research Commission, Grahamstown, South Africa. Nel, J.L., Roux, D.J., Maree, G., Kleynhans, C.J., Moolman, J., Reyers, B. et al., 2007, ‘Rivers in peril inside and outside protected areas: A systematic approach to conservation assessment of river ecosystems’, Diversity and Distributions 13(3), 341–352. https://doi.org/10.1111/j.1472-4642.2007.00308.x Nozais, C., Perissinotto, R. & Tita, G., 2005, ‘Seasonal dynamics of meiofauna in a South African temporarily open/closed estuary (Mdloti Estuary, Indian Ocean)’, Estuarine, Coastal and Shelf Science 62(1–2), 325–338. https://doi.org/10.1016/j.ecss.2004.09.020 O’Keeffe, J.H., 1989, ‘Conserving rivers in southern Africa’, Biological Conservation 49(4), 255–274. https://doi.org/10.1016/0006-3207(89)90047-5 Peters, L., Hillebrand, H. & Traunspurger, W., 2007, ‘Spatial variation of grazer effects on epilithic meiofauna and algae’, Journal of the North American Benthological Society 26(1), 78–91. https://doi.org/10.1899/0887-3593(2007)26[78:SVOGEO]2.0.CO;2 Peters, L., Scheifhacken, N., Kahlert, M. & Rothhaupt, K.O., 2005, ‘Note: An efficient in situ method for sampling periphyton in lakes and streams’, Archiv für Hydrobiologie 163(1), 133–141. https://doi.org/10.1127/0003-9136/2005/0163-0133 Petersen, R.M., Hoffman, A.C., Kotze, P. & Marr, S.M., 2017, ‘First record of the invasive Australian redclaw crayfish Cherax quadricarinatus (von Martens, 1868) in the Crocodile River, Kruger National Park, South Africa’, KOEDOE 59(1), 1–3. https://doi.org/10.4102/koedoe.v59i1.1435 Pfannkuche, O. & Thiel, H., 1988, ‘Sample processing’, in R.P. Higgins & H. Thiel (eds.), Introduction to the study of meiofauna, pp. 134–145, Smithsonian Institution Press, Washington, DC. Pillay, D. & Perissinotto, R., 2009, ‘Community structure of epibenthic meiofauna in the St. Lucia Estuarine Lake (South Africa) during a drought phase’, Estuarine, Coastal and Shelf Science 81(1), 94–104. https://doi.org/10.1016/j.ecss.2008.10.014 Ptatscheck, C., Gehner, S. & Traunspurger, W., 2020a, ‘Should we redefine meiofaunal organisms? The impact of mesh size on collection of meiofauna with special regard to nematodes’, Aquatic Ecology 54, 1135–1143. https://doi.org/10.1007/s10452-020-09798-2 Ptatscheck, C., Brüchner-Hüttemann, H., Kreuzinger-Janik, B., Weber, S. & Traunspurger, W., 2020b, ‘Are meiofauna a standard meal for macroinvertebrates and juvenile fish?’, Hydrobiologia 847, 2755–2778. https://doi.org/10.1007/s10750-020-04189-y Rasifudi, L., Addo-Bediako, A., Bal, K. & Swemmer, T.M., 2018, ‘Distribution of benthic macroinvertebrates in the Selati River of the Olifants River system, South Africa’, African Entomology 26(2), 398–406. https://doi.org/10.4001/003.026.0398 R Development Core Team, 2020, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, viewed n.d., from http://R-project.org. Rebecchi, L., Boschetti, C. & Nelson, D.R., 2020, ‘Extreme-tolerance in meiofaunal organisms: A case study with tardigrades, rotifers and nematodes’, Hydrobiologia 847, 2779. https://doi.org/10.1007/s10750-019-04144-6 Riddell, E.S., Govender, D., Botha, J., Sithole, H., Petersen, R.M. & Shikwambana, P., 2019, ‘Pollution impacts on the aquatic ecosystems of the Kruger National Park, South Africa’, Scientific African 6, e00195. https://doi.org/10.1016/j.sciaf.2019.e00195 Rimayi, C., Odusanya, D., Weiss, J.M., De Boer, J. & Chimuka, L., 2018, ‘Contaminants of emerging concern in the Hartbeespoort Dam catchment and the uMngeni River estuary 2016 pollution incident, South Africa’, Science of the Total Environment 627, 1008–1017. https://doi.org/10.1016/j.scitotenv.2018.01.263 Rogers, K. & Biggs, H., 1999, ‘Integrating indicators, endpoints and value systems in strategic management of the rivers of the Kruger National Park’, Freshwater Biology 41(2), 439–451. https://doi.org/10.1046/j.1365-2427.1999.00441.x Roux, D.J., Kleynhans, C.J. & Thirion, C., 1999a, ‘Biological monitoring and assessment of rivers as a basis for identifying and prioritising river management options’, Water Science and Technology 39(10–11), 207–210. https://doi.org/10.2166/wst.1999.0657 Roux, D.J., Kleynhans, C.J., Thirion, C., Hill, L., Engelbrecht, J.S., Deacon, A.R. et al., 1999b, ‘Adaptive assessment and management of riverine ecosystems: The Crocodile/Elands River case study’, Water SA 25, 501–511. Schenk, J., Höss, S., Brinke, M., Kleinbölting, N., Brüchner-Hüttemann, H. & Traunspurger, W., 2020, ‘Nematodes as bioindicators of polluted sediments using metabarcoding and microscopic taxonomy’, Environment International 143, 105922. https://doi.org/10.1016/j.envint.2020.105922 Schmid-Araya, J.M., Hildrew, A.G., Robertson, A., Schmid, P.E. & Winterbottom, J., 2002, ‘The importance of meiofauna in food webs: Evidence from an acid stream’, Ecology 83(5), 1271–1285. https://doi.org/10.1890/0012-9658(2002)083[1271:TIOMIF]2.0.CO;2 Seymore, T., 2014, ‘Bioaccumulation of Metals in Barbus Marequensis from the Olifants River, Kruger National Park and Lethal Levels of Manganese to Juvenile Oreochromis Mossambicus’, PhD thesis, University of Johannesburg. Shikwambana, P., Taylor, J.C., Govender, D. & Botha, J., 2021, ‘Diatom responses to river water quality in the Kruger National Park, South Africa’, Bothalia-African Biodiversity & Conservation 51(1), 1–11. https://doi.org/10.38201/btha.abc.v51.i1.5 Solomon, M., Zambatis, N., Biggs, H.C. & Mare, N., 1999, ‘Comparison of classifications commonly used as templates for management, scientific and GIS work in the Kruger National Park’, KOEDOE 42(2), 131–142. https://doi.org/10.4102/koedoe.v42i2.238 Traunspurger, W., Wilden, B. & Majdi, N., 2020, ‘An overview of meiofaunal and nematode distribution patterns in lake ecosystems differing in their trophic state’, Hydrobiologia 847, 2665–2679. https://doi.org/10.1007/s10750-019-04092-1 Van der Laan, M., Van Antwerpen, R. & Bristow, K.L., 2012, ‘River water quality in the northern sugarcane-producing regions of South Africa and implications for irrigation: A scoping study’, Water SA 38(1), 87–96. https://doi.org/10.4314/wsa.v38i1.11 Van der Velde, G., Nagelkerken, I., Rajagopal, S. & De Vaate, A.B., 2002, ‘Invasions by alien species in inland freshwater bodies in Western Europe: The Rhine Delta’, in E. Leppäkoski, S. Gollasch & S. Olenin (eds.), Invasive aquatic species of Europe. Distribution, impacts and management, pp. 360–372, Springer, Dordrecht. Venter, F.J. & Bristow, J.W., 1986, ‘An account of the geomorphology and drainage of the Kruger National Park’, KOEDOE 29(1), 117–124. https://doi.org/10.4102/koedoe.v29i1.524 Wolmarans, N.J., Bervoets, L., Gerber, R., Yohannes, Y.B., Nakayama, S.M., Ikenaka, Y. et al., 2021, ‘Bioaccumulation of DDT and other organochlorine pesticides in amphibians from two conservation areas within malaria risk regions of South Africa’, Chemosphere 274, 129956. https://doi.org/10.1016/j.chemosphere.2021.129956 Zullini, A., 2014, ‘Is a biogeography of freshwater nematodes possible?’, Nematology 16(1), 1–18. https://doi.org/10.1163/15685411-00002779 Appendix 1 FIGURE 1-A1: Mean daily discharge of Olifants, Sabie and Crocodile rivers before sampling over the period January–May 2019. FIGURE 2-A1: Upstream sampling sites and surrounding landscape. FIGURE 3-A1: (a) material used to sample meiobenthos: brush-sampler to scrape stones and cobbles underwater, perspex corer to sample sediment, sieve, funnel and tubes. (b) Localisation of sampling plots along a reach of the Crocodile River. The downstream plot (#1) was sampled first FIGURE 4-A1: Correlation plot showing upper triangular matrix of environmental variables measured in the water and sediment of three South African rivers. Correlation coefficients are displayed and clustered. Non-significant correlations are not shown. TABLE 1-A1: Relative taxonomic composition of the macrobenthic invertebrate community dwelling the riverbed of three rivers flowing through the Kruger National Park.