Next Article in Journal
Assessing Spatial Variability and Trends of Droughts in Eastern Algeria Using SPI, RDI, PDSI, and MedPDSI—A Novel Drought Index Using the FAO56 Evapotranspiration Method
Previous Article in Journal
Detection of Pesticides in Water through an Electronic Tongue and Data Processing Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benthic Macroinvertebrate Communities in Wadeable Rivers and Streams of Lao PDR as a Useful Tool for Biomonitoring Water Quality: A Multimetric Index Approach

by
Jutamas Sripanya
1,
Chanda Vongsombath
2,
Viengkhone Vannachak
3,
Kaewpawika Rattanachan
4,
Chutima Hanjavanit
1,
Wuttipong Mahakham
1,5 and
Narumon Sangpradub
1,*
1
Department of Biology, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Environmental Science, Faculty of Environmental Sciences, National University of Laos, Vientiane 7322, Laos
3
Department of Biology, Faculty of Science, National University of Laos, Vientiane 7322, Laos
4
Department of National Parks, Wildlife and Plant Conservation, Bangkok 10900, Thailand
5
Applied Taxonomic Research Center, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
*
Author to whom correspondence should be addressed.
Water 2023, 15(4), 625; https://doi.org/10.3390/w15040625
Submission received: 6 December 2022 / Revised: 1 February 2023 / Accepted: 2 February 2023 / Published: 5 February 2023
(This article belongs to the Section Biodiversity and Functionality of Aquatic Ecosystems)

Abstract

:
Lao PDR, a landlocked country in the lower Mekong River basin of Southeast Asia, has been considered a global biodiversity hotspot with a high level of biological endemism. In recent years, urban development and industrialization have affected the water quality of freshwater ecosystems in Lao PDR. However, the assessment of water quality in the country is primarily focused on a physicochemical method, while the application of a multimetric index (MMI) approach using benthic macroinvertebrates for biomonitoring in rivers and streams has not been established. MMI, based on benthic macroinvertebrates, is a biomonitoring tool that considers the effects of multiple anthropogenic impacts on benthic macroinvertebrate metrics associated with their biological attributes (e.g., taxa richness, composition, pollution tolerance, habits, and functional feeding) and aggregates individual metrics into a single value for assessing the water quality and health conditions of aquatic ecosystems. Here, we developed an MMI based on macroinvertebrate communities collected during 2016–2018 from 10 localities of streams and wadeable rivers in Lao PDR. Of the 54 potential metrics tested, 35 candidate macroinvertebrate metrics representing richness, composition, trophic structure, habit, and tolerance to pollution were selected, while 19 metrics were excluded. Of the 35-candidate metrics, a total of 11 core metrics (Total taxa, EPT taxa, Ephemeroptera taxa, %Diptera, %Plecoptera, %Tolerant, Beck’s biotic index, %Intolerant, Filterers taxa, %Sprawlers, and %Burrowers) were finally selected for the development of MMI based on their sensitivity, redundancy, and easy-to-apply tool for the biomonitoring program. These metrics can be used to distinguish the reference (seven sites) from stressed conditions (seven sites). In addition, the final MMI scores classified 40 sampling sites into four classes of water quality, including excellent (25%), good (10%), fair (60%), and poor (5%), which the conventional physicochemical method could not clearly distinguish. The Lao MMI developed in this study is an effective tool for evaluating the water conditions of sites affected by human activities, particularly agricultural areas, and, thus, is appropriate for use in future studies for assessing the ecological conditions of rivers and streams in the Mekong region.

1. Introduction

The Lao People’s Democratic Republic, or Lao PDR, is a landlocked country located in mainland Southeast Asia. Boarded by China, Vietnam, Cambodia, Thailand, and Myanmar, the country encompasses an area of 236,800 square kilometers and has a population of approximately 6.9 million [1]. As a global biodiversity hotspot with high levels of biological endemism, Lao PDR boasts diverse aquatic ecosystems, including lakes, ponds, rivers, and streams, which are home to a wide range of fish, amphibians, mollusk, and small aquatic invertebrate species, as well as aquatic macrophytes and algae [2,3,4]. These freshwater ecosystems play a vital role in supporting food and livelihood security for Lao PDR’s fisheries [5]. Thus, the aquatic ecosystems of Lao PDR are of great scientific interest and economic importance. However, in recent decades, Lao PDR’s aquatic ecosystems have been severely threatened by rapid population growth, urban development, industrialization, land use changes, dam construction, mining excavation, and the removal of riparian vegetation [6]. According to estimates from the last decade, about 35% of liquid effluent discharged into inland surface waters has been treated, and the quantity is unknown [7]. Furthermore, high concentrations of nutrients (NO3 and P), nitrogenous matter, and total suspended solids (TSS) have been observed in the Mekong River at Vientiane City in Lao PDR [8]. Additionally, the lower Mekong countries, including Lao PDR, are increasingly concerned about the environmental and agricultural damage caused by hydropower development and dam construction [6]. Taken together, these anthropogenic activities have resulted in the deterioration of water resources, which, in turn, has had a negative impact on the aquatic environment and biodiversity in Lao PDR.
To achieve the sustainable development and management of water resources, it is crucial to understand the relationship between human activity and the health of aquatic ecosystems, as well as to develop reliable tools for monitoring water quality. Historically, water quality has been primarily evaluated using physicochemical methods, which analyze various physical and chemical parameters of water against the established standards [9,10]. However, these methods are insufficient in assessing the consequences of certain human disturbances, as they only reflect the water quality at the time and place of sampling and do not directly measure the biological response to pollution [11,12]. Therefore, physicochemical methods cannot detect biogeochemical changes that occur in aquatic ecosystems over an extended period of time [13]. On the other hand, biomonitoring methods can reflect the consequences of changes in water quality and the biological response to human disturbance and pollution, both in the past and present [12,14,15,16,17,18]. Therefore, biomonitoring can provide a more accurate assessment of the true health of aquatic ecosystems, making it a valuable tool when used in conjunction with physicochemical approaches.
Benthic macroinvertebrates are a commonly used group of animals for biomonitoring in freshwater ecosystems [14,19]. These organisms are taxonomically and functionally diverse, with each group having specific environmental needs and ecological preferences. As a result, they show varying responses to changes in the ecological conditions of their habitat and water quality [17]. Environmental factors such as dissolved oxygen (DO), biological oxygen demand (BOD), nitrate nitrogen content, phosphorus compounds, water velocity, and depth, as well as the quality and quantity of available habitats, have been found to directly and indirectly, affect the diversity, composition, and distribution of benthic macroinvertebrates [15,20,21,22,23,24]. For example, in freshwater habitats with low DO concentrations, tolerant taxa such as Mollusca and Oligochaeta are often prominent [25]. Habitat heterogeneity in streams and rivers can also influence the functional composition and diversity of benthic macroinvertebrates. Research has shown that high levels of habitat heterogeneity promote faunal diversity, particularly among benthic macroinvertebrates [15]. Substrate and sediment grain size also play an important role in regulating the composition, distribution, diversity, and geographic preference of benthic macroinvertebrates. A small sediment grain size in sandy habitats affects the spatial distribution and density of benthic macroinvertebrate taxa [26,27]. Additionally, the decrease in the percentage of pebbles can lead to low habitat diversity for benthic fauna, resulting in a decline in the richness and density of benthic macroinvertebrates [28]. Therefore, changes in substrate conditions resulting from urbanization increased anthropogenic and agricultural activity, and the reduction in riparian vegetation can significantly threaten benthic macroinvertebrates and other aquatic organisms in streams and rivers [29]. As such, benthic macroinvertebrates can be used as indicators of various types of anthropogenic disturbance [30,31].
The use of biological indices for biomonitoring the water quality in freshwater ecosystems has been established for a significant amount of time [11,13,14,15,16,17]. However, there are limitations when using a single index for water quality monitoring [32]. While evaluation based on species occurrence can detect changes in water quality and pollution, these indices are relatively subjective and require knowledge of ecological life strategies and intensive sampling. As a result, a single index is often not able to accurately reflect the overall status of an aquatic ecosystem under various anthropogenic activities [17]. To overcome these limitations, a more integrated approach known as the multimetric index (MMI) has been introduced. This approach has gained increasing attention worldwide for its ability to integrate different environmental drivers that are impacted by human and natural activities over an extended period of time [11,12,14,15,16,17,18,19,30,32,33,34,35,36,37,38,39]. The MMI approach based on benthic macroinvertebrates has become a popular tool for biomonitoring programs in the United States [15,30,34] and European countries [32,35,40], as well as for monitoring the ecological health of rivers and streams in Asian countries such as China, Korea, Malaysia, Thailand, and Vietnam [17,18,37,41,42,43,44]. The MMI approach is flexible and can easily be adapted by adding or removing biological metrics or refining the index threshold values. However, it is important to exercise caution when using MMIs in different ecoregions or environmental gradients in specific geographical regions, which may have different reference conditions, anthropogenic pressures, and regional species assemblages [19,37].
In Lao PDR, water quality assessment is primarily focused on physicochemical and microbial analysis [45,46], while the use of biomonitoring is still relatively new [47]. Additionally, there is a lack of research on the application of the MMI approach using benthic macroinvertebrates for the biomonitoring of rivers and streams in Lao PDR. Furthermore, there is limited knowledge of the ecology and biodiversity of rivers and streams in this country, making the Mekong ecoregion in Lao PDR particularly important for assessing its current aquatic ecological conditions. The main objective of this study is to develop an MMI for biomonitoring wadeable rivers and streams in Lao PDR. To the best of our knowledge, this study is the first attempt to create a “Lao MMI”: an easy-to-apply and cost-effective tool for monitoring and evaluating the ecological condition of rivers and streams in Lao PDR.

2. Materials and Methods

2.1. Study Site

In the present study, 10 localities situated in three main ecoregions, the Lower Lancang, Khorat Plateau, and Kratie-Stung Treng, were used as sampling areas (Figure 1, Table 1). Reference sites were chosen based on the minimally disturbed areas, while the stressed sites were primarily assessed based on the presence of anthropogenic activities close to the studied areas. Samplings were conducted four times during the cool (December 2016 and November 2017) and hot seasons (April 2017 and 2018).

2.2. Measurement of Environmental Variables

Habitat assessment was taken at each site following the method of the USEPA [15]. Ten parameters of physical habitats were assessed, including the epifaunal substrate/available cover, embeddedness, velocity/depth regime, sediment deposition, channel flow status, channel alteration, frequency of riffles (or bends), bank stability, vegetative protection, riparian vegetative zone width. Total habitat scores (THS) were calculated and ranged from 0−200 points. THS values <130 points were considered as the stressed physical habitat condition [48].
Physico-chemical characteristics of water (15 parameters) were measured along with the collection of benthic macroinvertebrates at each of the sampling sites. Air and water temperatures (°C) were measured using a thermometer, while dissolved oxygen (DO, mg/L) was determined using a dissolved oxygen meter (YSI model 550A, Yellow Spring Instrument Co., Inc., Yellow Spring, MP, USA). Electrical conductivity (EC, µS/cm), total dissolved solids (TDS, mg/L), and pH were collected using a portable multi-probe meter (Hanna model HI 98129, Hanna Instruments, Woonsocket, Rhode Island, USA). In addition, 2 L water samples from each site were collected in sterilized plastic bottles and transported in an ice box (approximately 20 °C) to the laboratory for the standard measurement of orthophosphate (PO43−, mg/L), nitrate nitrogen (NO3-N), and chlorophyll a (µg/L) using previous methods [18,49]. For the five-day biochemical oxygen demand (BOD5) analysis, water samples were stored in dark-colored glass bottles at room temperature. The initial BOD5 concentration was measured with a dissolved oxygen meter (YSI model 550A). After an incubation period of 5 days at 20 °C in the dark, the water sample was checked again for DO content. The BOD5 was calculated from the difference in oxygen content between the start and end of the measurement [49]. Suspended solids (SS, mg/L) were determined using the photometric method [50]. Turbidity (NTU) was measured with a turbidimeter (Hach model 2100N, HACH Company, Loveland, CO, USA). Moreover, water channel width (m) and depth (cm) were measured using tape and a steel rod, respectively, while water velocity (m/s) was recorded using a pygmy water current meter (Genuine Gurley® current meter model D625 digital pygmy meter with Model 1100 Flow Velocity Indicator, Gurley Precision Instruments, New York, NY, USA).

2.3. Benthic Macroinvertebrate Sampling

Benthic macroinvertebrates were sampled from each site using the multihabitat approach by a D-frame net (30 × 30 cm2, 450 μm mesh size). A total of 20 kicks were collected from all habitats within a 100 m stretch and combined into a single sample. Animal samples were immediately fixed in 95% ethanol for observation in the laboratory. Then, they were sorted, subsampled, and counted for 300 ± 60 individuals [18]. Benthic macroinvertebrates were identified to the finest possible taxonomic level (usually family or genus) by compound and stereo microscope using the available identification keys [51,52]. All individuals were counted following identification. In addition, they were classified into a habit, functional feeding groups, and tolerance values according to previous reports [52,53,54,55,56,57,58,59]. Some taxa with unknown tolerance values were assigned based on their distribution in land use.

2.4. Site Classification

Based on the preliminary investigation, most of the study localities had quite a similar sediment composition except for Nam Houng downstream (HL), Nam Khan (NK), and Nam Hinboun downstream (BL), which demonstrated sites of human disturbances with less gradient substrate types. While Nam Houng upstream (HU), Nam Khod (KD), Nam Song (NS), Nam Thang (NT), Nam Hinboun upstream (BU), Xe Katam (KT), and Xe Namnoy (NN) were sites with fewer human stressors and were defined as having characteristics of dense riparian vegetation areas with high gradient substrate types (Figure 1, Table 1).
The reference sites used for developing a multimetric index were selected in the areas that were as minimally impaired as possible based on the site classification criteria [37,44]. Sample sites were classified into reference, intermediate, and stressed sites based on the following 14 environmental parameters: dissolved oxygen (DO), pH, electrical conductivity (EC), nitrate nitrogen (NO3-N), percent of land use (%land use), epifaunal substrate score (EpifauSub score), velocity score, sediment deposition score (SedDep score), bank stability score (BankSta score), vegetative protection score (VegPro score), percent of total habitat score (%THS), dam present, point source present, and the percentage of benthic macroinvertebrate taxa (%BMC taxa) [34,37,44,48,60]. A sampling site was classified as a reference site when all 14 criteria reached a standard or good level; sampling sites were classified as stressed sites if only 1 criterion failed with poor water quality or a high degree of development pressure [44]. The study site was designated as an intermediate site, as it did not fulfill all 14 criteria used to define a reference site and also did not meet any of the 14 criteria used to identify a stressed site.

2.5. Data Analysis

2.5.1. Habitat Score, Environmental Variables and Benthic Macroinvertebrate Composition

The mean ± SD of habitat scores, environmental parameters, and benthic macroinvertebrate communities were calculated for reference, intermediate, and stressed site groups. According to data distribution, one-way ANOVA or the Kruskal–Wallis test was used to determine significant differences in the total habitat score, environmental parameters, and benthic macroinvertebrate communities among the three groups (reference, intermediate, and stressed sites). In addition, if two groups were tested, the Man-Whitney U test was used to determine their statistical difference. All statistical analyses were performed using the IBM SPSS Statistics version 23 software [61]. Canonical correspondence analysis (CCA) was used to find relationships between environmental variables and benthic macroinvertebrate communities using PC-ORD software (version 5) [62].

2.5.2. MMI Development

Similarities in the benthic macroinvertebrate faunal composition between ecoregions were calculated with the ANOSIM by Past (Paleontological Statistics) version 4.03 [63]. An R-value close to one indicated a high difference, and close to zero indicated a low difference in the faunal composition between ecoregions. Independent sample T-test analysis was used to test significant differences in the number of benthic macroinvertebrate taxa between seasons and years. The benthic macroinvertebrate data from each sampling site was imported into the Ecological Data Application System (EDAS) version 3.3 for data management and the calculation of metrics [64].
The development of the macroinvertebrate-based MMI for Lao PDR was adapted from the procedure guided by the USEPA [15]. Fifty-four candidate metrics covering five categories of benthic macroinvertebrates, i.e., richness, composition, tolerance/intolerance, functional feeding group, and habit, were screened and examined according to the previous report [44]. The suitability of metrics for discriminating the reference and stressed conditions was assessed according to previous reports [37,44]. The comparison of the metric value between reference and stressed sites was performed with the aim of selecting the most appropriate metrics for Lao streams and wadeable rivers. The selection of potential metrics from the candidate metrics was based on two applications: sensitivity tests and redundancy tests [19,37].
In the first step of the sensitivity test, the capacity of each metric to differentiate the minimally disturbed sites (reference sites) from the most disturbed sites (stressed sites) was tested and visualized by box and whisker plot analysis. The metrics which showed a clear discriminatory power between the reference and stressed sites, i.e., no overlap of interquartile ranges, were selected [37]. The discriminatory power score of each metric was assessed according to the degree of overlap for medians and interquartile ranges [34]. Metrics with discriminatory power score 3 were mainly selected as potential metrics. If discriminatory power 2 was selected, a high (>70%) discrimination efficiency (DE) was also considered. DE was calculated as the percentage of stressed sites that scored below the 25th percentile of the reference sites for decreasing metrics with disturbance and above the 75th percentile for increasing metrics with disturbance [19,65]. Based on the discriminatory power or sensitivity score, 35 metrics were retained. These metrics were then selected based on their DE value [65].
In the second step, redundancy between the selected metrics was tested and evaluated using Pearson’s correlation analysis [37]. Pairs of metrics with r > 0.85 were considered redundant [17,41]. Redundant metrics were then removed, and at least one metric per category of benthic macroinvertebrates was retained [66,67]. Apart from the two testing steps, i.e., sensitivity and redundancy tests, ecological importance and the wider applicability of each metric were also used as criteria for selecting redundant metrics into the final core metrics. Spearman′s correlation was used to test a correlation among the core metrics, environmental variables, and grain sizes before MMI integration.
The final MMI was constructed using an integration of selected core metrics representing five categories of benthic macroinvertebrates, as mentioned above. Prior to integration, selected metrics were standardized by using the minimum value, the lower quartile (25%), and upper quartile (75%), and the maximum value of each metric dataset [44]. The scoring system of 1, 3, or 5 is representative of the respective threshold value (the 25th percentile and the 75th percentile) of each metric component, according to its response to environmental degradation [37,41]. For the metrics where numbers were expected to decrease with increasing disturbance or pollution, a score of 5 was awarded if the metric value of reference sites was above or equal to the lower quartile (25%), and a score of 3 and 1 was given for the bisected metric value of reference sites with less than the lower quartile. On the other hand, metrics that were expected to increase with increasing disturbance or pollution were awarded a score of 5 if the metric value was equal to or lower than the upper quartile (75%) of the reference sites and a score of 3 and 1 was assigned by the bisected metric value of reference sites with more than the upper quartile [34]. Thus, the appropriate quartile was used as a threshold depending on the type of response to degradation. A score of 5 indicated that the sample was part of the reference groups, a score of 3 indicated an intermediate condition, and a score of 1 indicated the highest deviation from the expected numbers for the reference sites [68].
After aggregating the individual core metric scores into a final MMI value for each sampling site, the index range was then divided into five classes: Class A (Excellent), which indicated that the site was comparable in condition to the reference biological conditions; Class B (Good) reflecting slight disturbance; Class C (Fair) indicating moderate disturbance; Class D (Poor) indicating highly disturbed biological integrity; and Class E (Impaired) reflecting severely disturbed biological integrity.
The sensitivity of Lao MMI was determined by assessing whether there was clear discrimination between the classified site group, i.e., the reference and stressed sites. The samples collected in 2016–2018 (n = 40) were used to validate Lao MMI for its appropriateness in assessing the water quality of wadeable rivers and streams in Lao PDR. Box-and-whisker plots were performed to distinguish the MMI value of the reference and stressed sites. Spearman′s correlation was used to test a correlation between the Lao MMI and some environmental variables.

3. Results

3.1. Site Classification and Environmental Characteristics

Of the 40 sampling sites from 10 localities, the site classification screening procedure identified seven reference sites (HU-C16, HU-C17, NT-C16, NT-C17, KT-C16, KT-H18, and NN-C16) and seven stressed sites (HL-C16, KN-C16, KH-N17, BL-H17, BL-C17, KH-H18, BL-H18). In the reference sites identified, most of them were mainly distributed in forest areas that were minimally impacted by human activities (Supplementary Table S1). On the contrary, the stressed sites were the areas characterized by high levels of agricultural land uses indicating anthropogenic activities (Supplementary Table S1).
Based on the data of physicochemical parameters and physical and benthic macroinvertebrate characteristics, as shown in Supplementary Table S1, the level of impairment of the studied sites was assessed and categorized into three groups, i.e., reference, intermediate, and stressed sites (Table 2). For statistical analysis, one-way ANOVA was used to test significant differences in the THS (Table 2) and total taxa (Table 3) due to the fact that both data had a normal distribution. While the Kruskal–Wallis test was used to determine significant differences in other data, including environmental variables (Table 2), the number of benthic macroinvertebrate individuals (Table 4), and the number of taxa and individuals in each order and data (Table 3 and Table 4) were not normally distributed. The Man-Whitney U test was used to determine statistical differences in the number of taxa and individuals in the order of Isopoda, Orthoptera, Unionoida, and Veneroida (Table 3 and Table 4), due to these orders only having datasets and the distribution of data was not normal.
The water of these site groups was moderately oxygenated with a slightly neutral pH. However, the Kruskal–Wallis test revealed that intermediate and stressed site groups showed higher average values of the EC, TDS, and NO3-N than those of the reference sites (p < 0.05) (Table 2). For the THS, the value of the reference site was significantly higher than those of stressed sites, but there was no significant difference in the total habitat score between the reference and the intermediate sites (Table 2).

3.2. Structure of Benthic Macroinvertebrate Community

A total of 13,651 benthic macroinvertebrate individuals in the samples were sorted and identified. They were classified into 3 phyla, 17 orders (Table 3), 102 families, and 241 taxa (Supplementary Table S2). Among them, insects (Arthropoda) were the most diverse group, with 220 taxa belonging to 11 orders and 86 families (91.29 %), followed by Mollusca (6.22 %) and Annelida (0.41 %) (Supplementary Table S2). The dominant groups of insects with the highest taxa richness were Ephemeroptera (18.62%), followed by Coleoptera (16.18%), Trichoptera (14.94 %), and Odonata (14.94%). Based on the relative total abundant and occurrence taxa (Supplementary Table S2), the Chironomidae (Diptera) showed the highest percentages (7.985% and 100%), followed by Choroterpes (Dilatognathus) (7.457%, 75%), Cheumatopsyche (5.487%, 87.5%), Potamyia (5.260%, 55%), and Baetis (3.780%, 85%), respectively. It was noted that the Chironomidae showed higher abundance in the stressed sites (Table 4). Some representative families and genera of these benthic macroinvertebrates are shown in Figure 2.
The highest number of taxa (taxa richness) of benthic macroinvertebrates was found in Nam Thang (128 taxa) and the lowest in Nam Hinboun downstream (77 taxa). The total taxa richness, Ephemeroptera taxa richness, Plecoptera taxa richness, and Trichoptera taxa richness were significantly different between stream conditions (p < 0.05). It was observed that these four taxa richness values were higher in reference sites than those of the stressed sites, while intermediate sites had moderate values of the taxa richness between the reference and stressed sites (Table 3). Moreover, the relative abundance of Plecoptera was significantly decreased in the stressed sites, while the relative abundance of Diptera and Chironomids was considerably higher at the stressed sites (Table 4).
The CCA ordination analysis showed the differences in benthic macroinvertebrate fauna between the three sites as explained by THS, water velocity (VEL), water channel width (WD), EC, and TDS. CCA axis 1 was primarily related to the gradient of THS and VEL, while axis 2 was related to EC, TDS, and WD. All sites in the reference sites were positively correlated to high THS and VEL. The reference sites were also clearly separated from the stressed sites, which were correlated to high EC and TDS. Intolerant taxa (tolerance values 1–3), including EPT, were prominently on the positive side of the CCA axis 1, which was correlated to high THS and VEL. The moderately tolerant taxa (tolerance values 7–8), e.g., some Ephemeroptera, Diptera, Odonata, Mollusca, and tolerant taxa (tolerance values 9–10), such as Chironomidae, Culicidae, and Oligochaeta were abundant on the negative side of the CCA axis 1, which was related to high EC and TDS (Figure 3). In addition, the ordination results confirmed the same trend as the previous results and strongly supported the highly significant relationships between physicochemical variables and benthic macroinvertebrate communities (r = 0.813, p = 0.001).

3.3. Metric Sensitivity and Selection

Before the metric sensitivity test, ANOSIM was used to test the differences in the fauna composition among different ecoregions. The results from this test indicated no significant difference in fauna composition among the three ecoregions (r = 0.1106, p > 0.05). Moreover, an independent sample T-test analysis revealed that the number of benthic macroinvertebrate taxa observed in different seasons (hot and cold) (p = 0.793) and years (2016–2018) (p = 0.651) were not significantly different. Then, box-and-whisker plot analysis was used to test if a metric was sensitive enough, i.e., if it could be used to discriminate between the reference and stressed sites. Of the 54 metrics initially evaluated, 35 metrics were considered sensitive showing a significant difference between the reference and stressed sites (Supplementary Figures S1–S5). The remaining 35 metrics represented five categories of benthic macroinvertebrates, i.e., richness (seven metrics), composition (eight metrics), tolerance (seven metrics), trophic structure (nine metrics), and habit (four metrics). In the sensitivity test, the discrimination power score (sensitivity score), together with a higher DE value (>70%), was used to select potential metrics to integrate into MMI. Metrics with a discriminatory power score of three and the highest DE value were mainly selected as potential metrics. If discriminatory power two was selected, a high (>70%) discrimination efficiency (DE) was also considered. In addition, a redundancy test was used to select the potential metrics. According to Pearson’s correlation analysis, metrics with a high redundancy showed high correlation values (r > 0.85) with each other, and those metrics represent how a similar category of the macroinvertebrate community was excluded from the candidate metrics. For example, in the tolerance/intolerance metric, the intolerant taxa metric was redundant with Beck’s biotic index metric (r = 0.863, p < 0.01); thus, the intolerant taxa metric was eliminated. Another example is the richness metric that Trichoptera taxa and EPTC taxa metrics were discarded due to redundancy with the total taxa metric (r = 0.897, 0.979, respectively, p < 0.05) (Supplementary Table S3). Apart from sensitivity and redundancy tests, ecological importance and the wider applicability of each metric were also used as criteria for selecting redundant metrics into final core metrics. Moreover, in order to obtain MMI covering five categories of benthic macroinvertebrates, at least one metric per category was retained. Regarding these criteria, a total of 11 core metrics that met the sensitivity and redundancy criteria with characteristics of simplicity, cost-benefit, and common use for biomonitoring were finally selected to integrate into the MMI. The 11-core metrics of the Lao MMI consisted of three richness metrics (Total taxa, EPT taxa, Ephemeroptera taxa), two composition metrics (%Diptera and %Pleoptera), three tolerance metrics (%Tolerant, Beck’s Biotic Index, and %Intolerant), one functional feeding group metric (Filterers taxa), and two habit metrics (%Sprawlers and %Burrowers) (Table 5). All of the 11-core metrics clearly distinguished reference and stressed sites (Figure 4). None of them showed partial or considerable interquartile overlaps. A non-parametric analysis utilizing Spearman’s rank correlation was conducted to examine the relationship between environmental variables and eleven metrics related to benthic macroinvertebrates. The results, presented in Supplementary Table S4, indicate that a number of the selected metrics exhibited a significant correlation with physicochemical parameters and sediment grain size. Specifically, metrics related to the total taxa (richness category), %Plecoptera (composition category), %Intolerant (tolerant/intolerant category), and Filterer taxa (functional feeding group category) were found to have a negative correlation with physical parameters such as EC and TDS. Additionally, the Filterer taxa and %Burrowers (habit category) metrics were found to have a negative correlation with the chemical parameter, i.e., nitrates. Furthermore, these metrics also demonstrated either a positive or negative correlation with the sediment grain size. For example, the total taxa, %intolerant, and filterer taxa metrics were positively correlated with boulder and pebble, but negatively correlated with sand and muck-mud, with the exception of the total taxa, which had a negative correlation with only muck-mud. Conversely, the %Burrower metric was found to be positively correlated with gravel, sand, and muck-mud but negatively correlated with boulder and cobble.

3.4. Development of Lao MMI

After selecting 11 metrics for integration into the Lao MMI, the value for the appropriate quartile of each of the 11 selected metrics at the reference sites was used as a threshold for segregating the maximum possible score from the lower score (Table 6). Using the metric scores in Table 6, the MMI was calculated by aggregating the scores of each of the 11 core metrics. The range of possible scores for the Lao MMI was determined by the minimum and maximum scores from <13 to ≥53. The range of the Lao MMI was then quadrisected in order to establish the five classes of water quality assessment; excellent condition (score ≥ 53), which pertained to the desired reference biological condition with a low degree of alteration in biological integrity; good (49–52) with good water quality and slightly disturbed biological integrity; fair (25–48) with fair quality of water and moderately disturbed biological integrity; poor (13–24), with poor quality of water and highly disturbed biological integrity, and impaired condition (score < 13) with very poor water quality and severe impairment of biological integrity (Table 7).

3.5. Validation of the Multimetric Index

In the assessment of the water quality of wadeable rivers and streams in Lao PDR, the developed Lao MMI was tested at the 14 sampling sites (seven reference and seven stressed sites) and was used to develop the index and another 26 sampling sites. Interestingly, the Lao MMI was able to discriminate the reference and stressed site groups (Table 8, Figure 5). Among the forty sampling sites tested, ten (25%) showed waters of excellent quality, four (10%) presented good quality, twenty-four (60%) had fair quality, and two (5%) had poor quality. On the other hand, the conventional physicochemical method established by the Department of Pollution Control of Lao PDR could not clearly distinguish the water quality of 40 sampling sites into different classes: excellent (one site), good (thirty-eight sites) and fair (one site) conditions (Table 8). Therefore, the Lao MMI developed in this study was sensitive enough to distinguish not only the reference and stressed sites but also intermediate environmental conditions. The Lao MMI developed in this study also showed a correlation with some environmental variables using Spearman’s correlation analysis. For example, THS was positively correlated with the index (r = 0.610, p < 0.001). In contrast, there were negative correlations with EC (r = −0.649, p < 0.001), TDS (r = −0.647, p < 0.001), and NO3-N (r = −0.348, p < 0.05) (Supplementary Table S5). As increasing EC, TDS, and NO3-N values were related to the sites closest to agricultural areas (Supplementary Table S1), this suggested that the developed MMI had good responsiveness to anthropogenic disturbance.

4. Discussion

Urbanization, industrialization, and agricultural land use have led to the widespread degradation of water quality in Southeast Asian countries, including the Lao People’s Democratic Republic (Lao PDR) [6,8]. The current methods of water quality assessment in Lao PDR primarily focus on physicochemical parameters, with microbial analysis being used as an optional component [69]. The use of benthic macroinvertebrates for biomonitoring water quality, which is commonly used in other regions, has not yet been implemented in Lao PDR. A review of the literature suggests that freshwater biomonitoring in Lao PDR has recently begun through the Mekong River Commission project [47]. Therefore, in this study, a benthic macroinvertebrate-based multimetric index (MMI) was developed to assess the water quality of some wadeable rivers and streams in Lao PDR.
The selection of appropriate reference sites is a fundamental aspect of developing an MMI for biomonitoring and evaluating water quality [37,70]. Ideally, these reference sites should be located in areas that are non-impaired or minimally impaired [70,71]. However, identifying such sites can be challenging, particularly in the lower Mekong region, where many aquatic ecosystems have been disturbed by human activities. To address this issue, the least disturbed areas from anthropogenic activities in each ecoregion were chosen as reference sites in this study. The data on benthic macroinvertebrates, together with the physicochemical properties of waters and habitats from seven reference and seven stressed sites of wadeable rivers and streams, were used to generate the Lao MMI (Supplementary Table S1). A total of 54 candidate metrics from five categories of benthic macroinvertebrates (richness, composition, tolerance/intolerance to pollution, habit, and functional feeding group) were initially considered for integration into the MMI, but some were eliminated due to a lack of effective discrimination between the reference and stressed sites. The final MMI included 11 metrics covering five categories of benthic macroinvertebrates, which were previously used for MMI development in various climatic regions [14,15,17,18,30,32,34,35,37,40,41,42,43,44,72,73,74,75] and were selected based on criteria such as cost-benefit, wide applicability, simplicity, and ease of calculation [11,17,19,37].
The taxa richness metric, a measure of the number of different species present, has been shown to be a reliable indicator of anthropogenic disturbance and exhibits good responsiveness to changes in water quality [15,37]. In this study, the total taxa (number of taxa), which are representative of the richness metric, was sufficiently sensitive to distinguishing reference and stressed sites. Previous research has demonstrated that high species richness is indicative of physical habitat diversity, good water quality, high availability of food resources, undisturbed conditions, and overall ecosystem health [14,15,37]. Furthermore, studies have shown that taxa richness, particularly that of sensitive taxa, decreases substantially in areas impacted by pollution, as pollution-intolerant species are lost [76,77,78]. In the current study, a decrease in taxa richness was observed in stressed sites, such as the Nam Khan, which may be attributed to moderate increases in nutrient availability from nearby agricultural activities. Previous studies have also reported the use of total taxa, the number of EPT taxa, and the number of Ephemoroptera taxa metrics in developing MMIs [37,41,42,43,44,71,79]. These metrics have been found to be excellent indicators in river and stream biomonitoring due to their sensitivity to perturbation and the fact that crucial environmental factors affect their spatial distribution [37].
The EPT (Ephemeroptera, Plecoptera, and Trichoptera) richness metric, which comprises sensitive taxa, is widely acknowledged as a significant bioindicator for determining the water quality of aquatic ecosystems [11,76,80]. Studies have demonstrated that human activities near rivers can have negative effects on the abundance and diversity of EPT taxa [80]. Since they like clean, cold running water, Plecoptera (stoneflies) are among the EPT that is particularly sensitive to environmental changes, specifically changes in water DO and temperature [81,82]. Due to their limited ability to move, stoneflies will inevitably respond to changes in the quality of the water [80,82]. Therefore, the existence of stoneflies suggests that environmental conditions in habitats are within their tolerance range [80]. Ephemeroptera, or mayflies, are another significant component of the EPT and are typically prevalent in highland streams with rather high DO concentrations. With rising habitat disturbance, pollution, and conductive water, the abundance of mayflies drastically decreases [83,84]. Thus, Ephemeroptera can be potentially used as a bioindicator of anthropogenic disturbance, except for some groups, e.g., Baetidae and Caenidae, which are capable of tolerating a broad range of anthropogenic disturbances [19,37,83]. Trichoptera (caddisflies), the final member of the EPT group, and some cased caddisflies construct their cases from leaves and other debris found in rivers and streams [85]. The habitats that are close to forests are likely to obtain available leaf materials from the surrounding trees; hence Trichopterans can use the leaf litter for their case construction and as a food source because they are shredders and collector-gatherers [85,86]. Given the ecological characteristics and functions of the EPT group, the EPT richness metric has become widely used in many biomonitoring programs due to its sensitivity to anthropogenic disturbance [19,80,83].
The composition metric is an essential core metric commonly used in conjunction with the richness metric for MMI development. The high degree of composition and the richness of certain groups of benthic macroinvertebrates, such as EPT taxa, can reflect minimally disturbed areas that provide heterogeneous habitats for the diverse niche partitioning of macroinvertebrate communities and support habitats for a diverse array of macroinvertebrates [86]. In contrast, rivers and streams that have been disturbed due to anthropogenic activities, such as agricultural land use and urbanization, tend to have a low composition, and richness of sensitive taxa, such as EPT [19,80,83]. Given that they are highly sensitive to human disturbance, composition and richness metrics are two core essential components for developing MMIs that have been previously used in the biomonitoring of river and stream ecosystems [11,15,37,41,86]. Additionally, certain macroinvertebrate groups can be tolerant (e.g., Diptera) or intolerant (e.g., EPT group) to pollution and human disturbance [15,37,86]. Diptera, including its diverse Chironomidae family and other tolerant taxa such as Oligochaetes (Annelida) and mollusk Bithynidae, are frequently predominant in freshwater ecosystems that are rich in increasing nutrient concentration and depleting dissolved oxygen concentration [19,52,86,87]. Due to the possession of hemoglobin in their body, Chironomidae and Oligochaetes are able to tolerate and survive in disturbed areas with depleted oxygen concentration because oxygen molecules can be trapped by hemoglobin molecules [52,86,88,89]. Both groups of macroinvertebrates are, therefore, useful indicators for assessing organic pollution in river and stream ecosystems and are important for the development of ecological integrity MMI and other biomonitoring tools globally [19,86,89].
In this study, two habit metrics, %Burrowers, and %Sprawlers, were selected for incorporation into the established MMI. Habit metrics take into account the mode of locomotion and position of benthic macroinvertebrates in aquatic environments [55]. The percentage of burrowers, primarily composed of chironomid larvae (Diptera), was found to be positively correlated with EC, TDS, nitrate, gravel, and muck mud. Similarly, the percentage of sprawlers, primarily composed of caenid mayfly nymphs, was found to be positively correlated with nitrate, gravel, and sand. Furthermore, EC, TDS, and nitrate values were found to be significantly higher in stressed sites, where the substrate grain size was dominated by fine sediment. This change in the substrate in stressed sites, may be caused by agricultural activities, leading to the deposition of fine sediment due to the instability of bank sides. Previous research also suggests that the tight packing of sand grains and muck mud reduces the trapping of detritus, thus limiting the availability of oxygen concentrations [90]. Thus, the occurrence of fine sediments due to agricultural activity becomes a preferred substratum for dipteran larvae. Similarly, many caenids are able to inhabit the surface of very fine (silty) substrata, with nymphs having the first abdominal gill modified to operculate gill cover and protect the underlying gills from siltation [55].
Regardless of their taxonomic affiliations, many species that obtain food in comparable ways were categorized in the same functional feeding category [91]. In the current study, the functional feeding group metric was included in the MMI as it provided insight into the degradation of ecosystem functioning by measuring the relative abundance and proportion of different functional groups in response to disturbances that affect food availability. It has been found that certain functional feeding groups, such as shredders and scrapers, may be more sensitive to environmental changes, while others, such as gatherers and filterers, may be more tolerant of pollution [34,92]. Many of the filterers found in this study belonged to the family Hydropsychidae (Trichoptera) (Figure 2), which are able to tolerate moderately polluted water [52]. Although Hydropsychidae members are moderately tolerant, they tended to decrease in stressed sites in this study, as previously observed in stream biomonitoring studies in the United States [93] and Thailand [94]. It is important to note that there are limitations to the use of functional feeding groups as it can be difficult to properly assign some taxa to functional feeding groups as species can change their feeding mode during the developmental stage. However, the use of filterer taxa in this study is reliable as most of them are trichopterans in the Hydropsychidae family, which have morpho-behavioral characteristics that are easily classified as filterers [95]. To ensure the proper assignment of functional feeding groups, guidelines from various works in the literature for categorizing the group based on life stage are recommended [15,95].
The assessment of water quality based solely on physicochemical parameters established by the Department of Pollution Control of Lao PDR did not effectively distinguish between the reference and stressed sites. However, the Lao MMI, which is based on benthic macroinvertebrate data, was able to differentiate between reference and stressed sites with a high level of discrimination (Table 8), indicating its sensitivity to environmental changes. Additionally, eleven metrics from the MMI were found to be significantly correlated with environmental variables and sediment grain size (see Supplementary Table S4), suggesting that the MMI was responsive to anthropogenic disturbances that could alter the distribution and abundance of benthic macroinvertebrates in rivers and streams. This is supported by previous research, which has shown that changes in the substrate condition due to anthropogenic and agricultural activity can affect the distribution and abundance of benthic macroinvertebrates [29,96]. For instance, the richness and density of benthic macroinvertebrates clearly decreased in sandy habitats, leading to a reduction in food sources for fish [26].
Prior to being employed in biomonitoring, the development of MMIs typically needs to undergo a stability test for an appropriate amount of time in order to prevent errors or misinterpretations about the health of river and stream ecosystems caused by the seasonal variability of the index [15,17,41]. In some cases, studies have revealed that monsoonal rains can cause seasonal differences in the assemblage structure of macroinvertebrates [97], but other studies have found that seasonality did not affect the ecological status by using MMI in tropical Southeast countries such as Vietnam [17], Thailand [44], and Malaysia [41]. The results of this study are consistent with these reports, and the 11-core metrics are, therefore, sufficiently qualified as biomonitoring components for assessing water quality and ecological integrity in wadeable rivers and streams in Lao PDR.
Rapid biomonitoring methods, which use coarse taxonomy to save time and money, are widely accepted for their practicality [78]. It is acknowledged that the process of identifying benthic organisms at a detailed level is a costly and time-consuming endeavor, requiring a comprehensive library of keys, reference collections, specialized knowledge, techniques, and equipment [78,98]. For detailed identification at the species level, taxonomic knowledge may be incomplete as keys and descriptions of specific life stages or sexes for many species are not readily available or easily accessible. Additionally, the use of morphology-based taxonomies may be limited by intraspecific variability and the difficulty of distinguishing between closely related species [78]. Therefore, the use of a higher taxonomic resolution leads to more accurate and informative results with low cost [78], although the use of a genus or species level could be useful to detect specific pollution [11]. In this study, the Lao MMI was designed for use in a biomonitoring program in the Mekong region, including the Lao PDR, where available keys for the identification of benthic macroinvertebrates at the phylum, class, order, family, and genus levels were established [51]. The advantage of the developed Lao MMI is that it is suitable for ecologists who are not familiar with benthic macroinvertebrate taxonomy, and it saves time and cost as the Lao MMI utilizes coarse taxonomy for sample diagnosis. However, it should be noted that the MMI developed here has some limitations. For example, it can only be used to monitor freshwater ecosystems in the explored ecoregions, particularly the Lower Lancang, Khorat Plateau, and Kratie-Stung Treng, as MMIs developed from different ecoregions or environmental gradients in specific geographical regions must be used with caution due to the differences in reference conditions, anthropogenic pressures, and regional species assemblages [19,37]. Additionally, the application of MMI in various field areas is recommended to ensure the stability of the MMI in the biomonitoring program.

5. Conclusions

In the current study, benthic macroinvertebrate-based MMI has been developed to assess the water quality of streams and wadeable rivers in Lao PDR. A total of 40 samples were collected across the country, and 11 core metrics, including Total taxa, EPT taxa, Ephemeroptera taxa, % Diptera, %Plecoptera, %Tolerant, Beck’s biotic index, %Intolerant, Filterers taxa, %Sprawlers, and %Burrowers, were selected for incorporation into the final index. The developed Lao MMI was found to be sensitive enough to differentiate between minimal human disturbance (reference sites) and areas of anthropogenic impact (stressed sites). Additionally, the Lao MMI had an advantage over conventional physicochemical methods as it was able to clearly classify the water quality of 40 sampling sites into four classes, while the physicochemical method was unable to do so. Furthermore, the Lao MMI is cost-effective and easy to use for biomonitoring as it primarily utilizes family and EPT genera levels for benthic macroinvertebrate diagnosis. To the best of our knowledge, this research represents the first attempt to utilize a biomonitoring method based on multimeric data generated from benthic macroinvertebrates for monitoring the water quality of rivers and streams in Lao PDR. Further research in other areas of rivers and streams is needed to validate the stability of the Lao MMI in biomonitoring programs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15040625/s1, Table S1: Site classification based on environmental parameter and benthic macroinvertebrate da-ta; Table S2: Data on relative abundance, occurrence, tolerance values (TolVal), functional feeding groups (FFGs), and habits of benthic macroinvertebrates observed from 40 sampling sites; Table S3: Pearson’s correlation analysis between candidate metrics of benthic macoinvertebrates in the reference sites; Table S4: Spearman’s correlation coefficients among the eleven core metrics, envi-ronmental variables and grian sizes; Table S5: Spearman′s correlation between final index scores of Lao MMI and some environmental variables; Figure S1: Box and whisker plots of benthic ma-croinvertebrate metrics to discriminate between reference and stressed sites. A-G richness catego-ry: A, total taxa; B, EPT taxa; C, Ephemeroptera taxa; D, Plecoptera taxa; E, Trichoptera taxa; F, EPTC taxa; G, Coleoptera taxa: H-O composition category: H, %EPT; Figure S2: Box and whisker plots of benthic macroinvertebrate metrics to discriminate between reference and stressed sites. H-O composition category: I, %Ephemeroptera; J, Margalef’s index; K, %Odonata; L, %Chirono-midae; M, %Diptera; N, %Plecoptera; O, %Trichoptera: P-V tolerance/intolerance category: P, In-tolerant taxa; Figure S3: Box and whisker plots of benthic macroinvertebrate metrics to discrimi-nate between reference and stressed sites. P-V tolerance/intolerance category: Q, %Tolerant; R, %Dominant taxon; S, Beck’s Biotic Index; T, Simpson Index; U, Hilsenhof’s Biotic Index; V, %In-tolerant: W-AE functional feeding group category: W, %Filterers; X, %Scrapers; Figure S4: Box and whisker plots of benthic macroinvertebrate metrics to discriminate between reference and stressed sites. W-AE functional feeding group category: Y, %Collectors; Z, Collector taxa; AA, Fil-terer taxa; AB, Predator taxa; AC, Scraper taxa; AD, %Shredders; AE, Shredder taxa: AF-AI habit category: AF, Clinger taxa; Figure S5: Box and whisker plots of benthic macroinvertebrate metrics to discriminate between reference and stressed sites. AF-AI habit category: AG, %Clingers; AH, %Sprawlers; AI, %Burrowers.

Author Contributions

Conceptualization, N.S.; methodology, C.H., C.V., J.S., K.R., N.S., and V.V.; validation, J.S., N.S., and W.M.; formal analysis, J.S. and N.S.; investigation, C.H., C.V., J.S., N.S., V.V., and W.M.; resources, N.S.; data curation, J.S. and N.S.; writing—original draft preparation, J.S., N.S., and W.M.; writing—review and editing, C.H., C.V., K.R., N.S., V.V., and W.M.; visualization, J.S., N.S., and W.M.; supervision, N.S.; funding acquisition, N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Center of Excellence on Biodiversity, Thailand, grant number BDC-PG1-159008.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Department of Natural Resources and Environment of Luangprabang, Xaignabouly, Vientiane, Khammouan and Xekong Provinces for allowing us to access the studied areas.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. The World Factbook. Available online: https://www.cia.gov/the-world-factbook/ (accessed on 19 November 2022).
  2. Japan International Cooperation Agency. Profile on Environmental and Social Considerations in Lao P.D.R.; Japan International Cooperation Agency: Tokyo, Japan, 2013; pp. 1–333. [Google Scholar]
  3. MoNRE; IUCN. Fifth National Report to the United Nations Convention on Biological Diversity; Department of Forest and Natural Resources Management, Ministry of Natural Resource and Environment: Vientiane, Laos, 2016; pp. 1–91. [Google Scholar]
  4. Singkham, P. An Introduction to the Fisheries of Lao PDR; Mekong River Commission: Phnom Penh, Cambodia, 2013; pp. 1–56. [Google Scholar]
  5. UNDP. Available online: https://www.la.undp.org/ (accessed on 7 January 2023).
  6. Soutullo, J. The Mekong River: Geopolitics over Development, Hydropower and the Environment; European Union: Brussels, Belgium, 2019; pp. 1–55. [Google Scholar]
  7. Mekong River Commission Lao National Mekong Committee. Annual Report on Water Quality Data Assessment Lao PDR; Water Quality Monitoring Network in Mekong River Basin, Natural Resources and Environment Institute, Ministry of Natural Resources and Environment: Vientiane, Laos, 2016; pp. 1–46. [Google Scholar]
  8. Chea, R.; Grenouillet, G.; Lek, S. Evidence of Water Quality Degradation in Lower Mekong Basin Revealed by Self-Organizing Map. PLoS ONE 2016, 11, e0145527. [Google Scholar] [CrossRef] [PubMed]
  9. Sharma, R.; Kumar, R.; Satapathy, S.C.; Al-Ansari, N.; Singh, K.K.; Mahapatra, R.P.; Agarwal, A.K.; Le, H.V.; Pham, B.T. Analysis of Water Pollution Using Different Physicochemical Parameters: A Study of Yamuna River. Front. Environ. Sci. 2020, 8, 581591. [Google Scholar] [CrossRef]
  10. Braga, F.H.R.; Dutra, M.L.S.; Lima, N.S.; Silva, G.M.; Miranda, R.C.M.; Firmo, W.C.A.; Moura, A.R.L.; Monteiro, A.S.; Silva, L.C.N.; Silva, D.F.; et al. Study of the Influence of Physicochemical Parameters on the Water Quality Index (WQI) in the Maranhão Amazon, Brazil. Water 2022, 14, 1546. [Google Scholar] [CrossRef]
  11. Tampo, L.; Kaboré, I.; Alhassan, E.H.; Ouéda, A.; Bawa, L.M.; Djaneye-Boundjou, G. Benthic Macroinvertebrates as Ecological Indicators: Their Sensitivity to the Water Quality and Human Disturbances in a Tropical River. Front. Water 2021, 3, 662765. [Google Scholar] [CrossRef]
  12. Abdelkarim, M.S. Biomonitoring and bioassessment of running water quality in developing countries: A case study from Egypt. Egypt. J. Aquat. Res. 2020, 46, 313–324. [Google Scholar] [CrossRef]
  13. Mezgebu, A. A review on freshwater biomonitoring with benthic invertebrates in Ethiopia. Environ. Sustain. Indic. 2022, 14, 100174. [Google Scholar] [CrossRef]
  14. Morse, J.C.; Bae, Y.J.; Munkhjargal, G.; Sangpradub, N.; Tanida, K.; Vshivkova, T.S.; Wang, B.; Yang, L.; Yule, C.M. Freshwater Biomonitoring with Macroinvertebrates in East Asia. Front. Ecol. Environ. 2007, 5, 33–42. [Google Scholar] [CrossRef]
  15. Barbour, M.T.; Gerritsen, J.; Snyder, B.D.; Stribling, J.B. Rapid Bioassessment Protocols for Use in Streams and Wadeable Rivers: Periphyton, Benthic Macroinvertebrates and Fish, 2nd ed.; U.S. Environmental Protection Agency, Office of Water: Washington, DC, USA, 1999; pp. 1–339.
  16. Karr, J.R.; Allan, J.D.; Benke, A.C. River Conservation in the United States and Canada. Global Perspectives on River Conservation: Science, Policy, Practice; J. Wiley: Chichester, UK, 2000; pp. 3–39. [Google Scholar]
  17. Nguyen, H.H.; Everaert, G.; Gabriels, W.; Hoang, T.H.; Goethals, P.L.M. A multimetric macroinvertebrate index for assessing the waterquality of the Cau river basin in Vietnam. Limnologica 2014, 45, 16–23. [Google Scholar] [CrossRef]
  18. Boonsoong, B.; Sangpradub, N.; Barbour, M.T. Development of rapid bioassessment approaches using benthic macroinvertebrates for Thai streams. Environ. Monit. Assess. 2009, 155, 129–147. [Google Scholar] [CrossRef]
  19. Lakew, A.; Moog, O. A multimetric index based on benthic macroinvertebrates for assessing the ecological status of streams and rivers in central and southeast highlands of Ethiopia. Hydrobiologia 2015, 751, 229–242. [Google Scholar] [CrossRef]
  20. Rosenberg, D.M.; Resh, V.H. (Eds.) Introduction to freshwater biomonitoring and benthic macroinvertebrates. In Freshwater Biomonitoring and Benthic Macroinvertebrate; Chapmam and Hall: New York, NY, USA, 1993; pp. 1–10. [Google Scholar]
  21. Koudenoukpo, C.Z.; Chikou, A.; Toko, I.I.; Togouet, S.H.Z.; Tchakonté, S.; Hazoume, R.; Piscart, C. Diversity of aquatic macroinvertebrates in relationship with the environmental factors of a lotic ecosystem in tropical region: The Sô river in South-East of Benin (West Africa). J. Entomol. Zool. Stud. 2017, 5, 1–10. [Google Scholar]
  22. Souilmi, F.; Ghedda, K.; Fahde, A.; EL Fihri, F.Z.; Tahraoui, S.; Elasri, F.; Malki, M. Taxonomic diversity of benthic macroinvertebrates along the Oum Er Rbia River (Morocco): Implications for water quality bio-monitoring using indicator species. West Afr. J. Appl. Ecol. 2019, 27, 137–149. [Google Scholar]
  23. Brysiewicz, A.; Czerniejewski, P.; Dabrowski, J.; Formicki, K. Characterisation of Benthic Macroinvertebrate Communities in Small Watercourses of the European Central Plains Ecoregion and the Effect of Different Environmental Factors. Animals 2022, 12, 606. [Google Scholar] [CrossRef] [PubMed]
  24. Sripanya, J.; Rattanawilai, K.; Vongsombath, C.; Vannachak, V.; Hanjavanit, C.; Sangpradub, N. Benthic Macroinvertebrates and Trichoptera Adults for Bioassessment Approach in Streams and Wadeable Riversin Lao People’s Democratic Republic. Trop. Nat. Hist. 2022, 22, 12–24. [Google Scholar]
  25. Souto, R.M.G.; Facure, K.G.; Pavanin, L.A.; Jacobucci, G.B. Influence of environmental factors on benthic macroinvertebrate communities of urban streams in Vereda habitats, Central Brazil. Acta Limnol. Bras. 2011, 23, 293–306. [Google Scholar] [CrossRef]
  26. Leitner, P.; Hauer, C.; Ofenböck, T.; Pletterbauer, F.; Schmidt-Kloiber, A.; Graf, W. Fine sediment deposition affects biodiversity and density of benthic macroinvertebrates: A case study in the freshwater pearl mussel river Waldaist (Upper Austria). Limnologica 2015, 50, 54–57. [Google Scholar] [CrossRef]
  27. Pacioglu, O.; Duţu, F.; Pavel, A.B.; Duţu, L.T. The influence of hydrology and sediment grain-size on the spatial distribution of macroinvertebrate communities in two submerged dunes from the Danube Delta (Romania). Limnetica 2022, 41, 85–100. [Google Scholar] [CrossRef]
  28. Rezende, R.S.; Santos, A.M.; Henke-Oliveira, C.; Gonçalves, J.F., Jr. Effects of spatial and environmental factors on benthic a macroinvertebrate community. Zoologia 2014, 31, 426–434. [Google Scholar] [CrossRef]
  29. Zweig, L.D.; Rabeni, C.F. Biomonitoring for deposited sediment using benthic invertebrates: A test on 4 Missouri streams. J. N. Am. Benthol. Soc. 2001, 20, 643–657. [Google Scholar] [CrossRef]
  30. Shull, D.R.; Smith, Z.M.; Selckmann, G.M. Development of a benthic macroinvertebrate multimetric index for large semiwadeable rivers in the Mid-Atlantic region of the USA. Environ. Monit. Assess. 2019, 191, 22. [Google Scholar] [CrossRef]
  31. Indicators: Benthic Macroinvertebrates. Available online: https://www.epa.gov/national-aquatic-resource-surveys/indicators-benthic-macroinvertebrates (accessed on 19 November 2022).
  32. Gabriels, K.; Lock, N.; De Pauw, N.; Goethals, P.L.M. Multimetric macroinvertebrate index Flanders (MMIF) for biological assessment of rivers and lakes in Flanders (Belgium). Limnologica 2010, 40, 199–207. [Google Scholar] [CrossRef]
  33. Mereta, S.T.; Boets, P.; Meester, L.D.; Goethals, P.L.M. Development of a multimetric index based on benthic macroinvertebrates for the assessment of natural wetlands in Southwest Ethiopia. Ecol. Indic. 2013, 29, 510–521. [Google Scholar] [CrossRef]
  34. Barbour, M.T.; Gerritsen, J.; Griffith, G.E.; Frydenborg, R.; McCarron, E.; White, J.S.; Bastian, M.L. A framework for biological criteria for Florida streams using benthic macroinvertebrates. J. N. Am. Benthol. Soc. 1996, 15, 185–211. [Google Scholar] [CrossRef]
  35. Vlek, H.E.; Verdonschot, P.F.M.; Nijboer, R.C. Towards a multimetric index for the assessment of Dutch streams using benthic macroinvertebrates. Hydrobiologia 2004, 516, 173–189. [Google Scholar] [CrossRef]
  36. Hering, D.; Moog, O.; Sandin, L.; Verdonschot, P.F.M. Overview and application of the AQEM assessment system. Hydrobiologia 2004, 516, 1–20. [Google Scholar] [CrossRef]
  37. Jun, Y.C.; Won, D.H.; Lee, S.H.; Kong, D.S.; Hwang, S.J. A Multimetric Benthic Macroinvertebrate Index for the Assessment of Stream Biotic Integrity in Korea. Int. J. Environ. Res. Public Health 2012, 9, 3599–3628. [Google Scholar] [CrossRef]
  38. Karr, J.R. Assessment of Biotic Integrity Using Fish Communities. Fisheries 1981, 6, 21–27. [Google Scholar] [CrossRef]
  39. Vadas, R.L., Jr.; Hughes, R.M.; Bae, Y.J.; Baek, M.J.; Gonzáles, O.C.B.; Callisto, M.; de Carvalho, D.R.; Chen, K.; Ferreira, M.T.; Fierro, P.; et al. Assemblage-based biomonitoring of freshwater ecosystem health via multimetric indices: A critical review and suggestions for improving their applicability. Water Biol. Secur. 2022, 1, 100054. [Google Scholar] [CrossRef]
  40. Haase, P.; Lohse, S.; Pauls, S.; Schindehütte, K.; Sundermann, A.; Rolauffs, P.; Hering, D. Assessing streams in Germany with benthic invertebrates: Development of a practical standardised protocol for macroinvertebrate sampling and sorting. Limnologica 2004, 34, 349–365. [Google Scholar] [CrossRef]
  41. Arman, N.Z.; Salmiati, S.; Said, M.I.M.; Aris, A. Development of macroinvertebrate-based multimetric index and establishment of biocriteria for river health assessment in Malaysia. Ecol. Indic. 2019, 104, 449–458. [Google Scholar] [CrossRef]
  42. Huang, Q.; Gao, J.; Cai, Y.; Yin, H.; Gao, Y.; Zhao, J.; Liu, L.; Huang, J. Development and application of benthic macroinvertebrate-based multimetric indices for the assessment of streams and rivers in the Taihu Basin, China. Ecol. Indic. 2014, 48, 649–659. [Google Scholar] [CrossRef]
  43. Liu, L.; Xu, Z.; Yin, X.; Li, F.; Dou, T. Development of a multimetric index based on benthic macroinvertebrates for the assessment of urban stream health in Jinan City, China. Environ. Monit. Assess. 2017, 189, 205. [Google Scholar] [CrossRef] [PubMed]
  44. Rattanachan, K.; Boonsoong, B.; Getwongsa, P.; Uttaruk, Y.; Sangpradub, N. A benthic macroinvertebrates multimetric index for assessment of the ecological integrity of northeast streams, Thailand. Environ. Asia 2016, 9, 186–194. [Google Scholar]
  45. Mekong River Commission. 2018 Lower Mekong Water Quality Monitoring Report; Mekong River Commission: Vientiane, Laos, 2021; pp. 1–60. [Google Scholar]
  46. Brindha, K.; Pavelic, P.; Sotoukee, T. Environmental assessment of water and soil quality in the Vientiane Plain, Lao PDR. Groundw. Sustain. Dev. 2019, 8, 24–30. [Google Scholar] [CrossRef]
  47. Mekong River Commission. Report on the 2017 Biomonitoring Survey of the Lower Mekong River and Selected Tributaries; Mekong River Commission: Vientiane, Laos, 2019; pp. 1–85. [Google Scholar]
  48. Tetra Tech. A Stream Condition Index for West Virginia Wadeable Streams; Tetra Tech, Inc.: Baltimore, MD, USA, 2000; pp. 1–80. [Google Scholar]
  49. APHA. Standard Methods for the Examination of Water and Wastewater, 20th ed.; APHA: Washington, DC, USA, 2001. [Google Scholar]
  50. Krawczyk, D.; Gonglewski, N. Determining Suspended Solids Using a Spectrophotometer. Sewage. Ind. Waste 1959, 31, 1159–1164. [Google Scholar]
  51. Sangpradub, N.; Boonsoong, B. Identification of Freshwater Invertebrates of the Mekong River and Tributaries; Mekong River Commission: Vientiane, Laos, 2006; pp. 1–274. [Google Scholar]
  52. Morse, J.C.; Yang, L.; Tian, L. Aquatic Insects of China Useful for Monitoring Water Quality; Hohai University Press: Nanjing, China, 1994; pp. 1–570. [Google Scholar]
  53. Wiggins, G.B. Larvae of the North American Caddisfly Genera (Trichoptera); University of Toronto Press: Toronto, ON, Canada, 1996; pp. 1–457. [Google Scholar]
  54. Boonsoong, B. Rapid Bioassessment of Thai Streams with Benthic Macroinvertebrates. Ph.D. Thesis, Khon Kean University, Khon Kean, Thailand, 2007. [Google Scholar]
  55. Merritt, R.T.; Cummins, K.W.; Berg, M.B. An Introduction to the Aquatic Insects of North America, 5th ed.; Kendall Hunt Publishing Company: Dubuque, IA, USA, 2008; pp. 1–1158. [Google Scholar]
  56. Getwongsa, P. Development of Biotic Index for Rapid Bioassessment in Mekong 2 Basin. Ph.D. Thesis, Khon Kean University, Khon Kean, Thailand, 2010. [Google Scholar]
  57. Uttaruk, P. The Development of Ecological Indicators for Streams: Using Benthic Macroinvertebrates to Rapid Bioassessment for Wadeable Streams in Chi and Mun Headwater, the Northeast Thailand. Ph.D. Thesis, Khon Kean University, Khon Kaen, Thailand, 2011. [Google Scholar]
  58. Rattanachan, K. Benthic Macroinvertebrate Multimetric Index for Rapid Bioassessment of Streams in Thailand. Ph.D. Thesis, Khon Kean University, Khon Kaen, Thailand, 2017. [Google Scholar]
  59. Somnark, R. Community Structure of Freshwater Benthic Macroinvertebrates, Phenology of Some Species of Aquatic Insects and Leaf Litter Breakdown in Nam Nao National Park. Ph.D. Thesis, Khon Kean University, Khon Kaen, Thailand, 2018. [Google Scholar]
  60. Pollution Control Department of Thailand. Water Quality Criteria and Standards in Thailand; Pollution Control Department, Ministry of Natural Resources and Environment: Bangkok, Thailand, 1997. [Google Scholar]
  61. IBM Corp. IBM SPSS Statistics for Windows. Version 23.0; IBM Corp.: Armonk, NY, USA, 2015. [Google Scholar]
  62. McCune, B.; Mefford, M.J. PC-ORD. Mutivariate Analysis of Ecological Data. Version 5; MjM Software: Gleneden Beach, OR, USA, 2006. [Google Scholar]
  63. Hammer, O. PAST—PAlaeontological STatistics, Version 4.03. Reference Manual; Natural History Museum, University of Oslo: Oslo, Norway, 2020. [Google Scholar]
  64. Tetra Tech. Ecological Data Application System (EDAS): A User’s Manual; Tetra Tech, Inc.: Baltimore, MD, USA, 2000. [Google Scholar]
  65. Stribling, J.B.; Jessup, B.K.; Gerritsen, J. Development of Biological and Physical Habitat Criteria for Wyoming Streams and Their Use in the TMDL Process; Tetra Tech, Inc.: Baltimore, MD, USA, 2000. [Google Scholar]
  66. Stoddard, J.L.; Herlihy, A.T.; Peck, D.V.; Hughes, R.M.; Whittier, T.R.; Tarquinio, E. A process for creating multimetric indices for large-scale aquatic surveys. J. N. Am. Benthol. Soc. 2008, 27, 878–891. [Google Scholar] [CrossRef]
  67. Wu, H.; Lu, K.; Lu, X.; Xue, Z. A Macroinvertebrate Multimetric Index for the Bioassessment of Wetlands Adjacent to Agriculture Fields in the Sanjiang Plain, China. Chin. Geogr. Sci. 2019, 29, 974–984. [Google Scholar] [CrossRef]
  68. Ferreira, W.R.; Paiva, L.T.; Callisto, M. Development of a benthic multimetric index for biomonitoring of a neotropical watershed. Braz. J. Biol. 2011, 71, 15–25. [Google Scholar] [CrossRef]
  69. Department of Pollution Control. Standard of National Environment; Ministry of Natural Resources and Environment: Vientiane, Laos, 2017; pp. 15–19. [Google Scholar]
  70. Sánchez-Montoya, M.M.; Vidal-Abarca, M.R.; Puntı, T.; Poquet, J.M.; Prat, N.; Rieradevall, M.; Alba-Tercedor, J.; Zamora-Muñoz, C.; Toro, M.; Robles, S.; et al. Defining criteria to select reference sites in Mediterranean streams. Hydrobiologia 2009, 619, 39–54. [Google Scholar] [CrossRef]
  71. Li, Y.; Li, Y.; Xu, Z.; Li, L. Assessment of the Huntai River in China using a multimetric index based on fish and macroinvertebrate assemblages. J. Freshw. Ecol. 2016, 31, 169–190. [Google Scholar] [CrossRef]
  72. Ofenbӧck, T.; Otto Moog, O.; Gerritsen, J.; Barbour, M. A stressor specific multimetric approach for monitoring running waters in Austria using benthic macro-invertebrates. Hydrobiologia 2004, 516, 251–268. [Google Scholar] [CrossRef]
  73. Stoddard, J.L.; Larsen, D.P.; Hawkins, C.P.; Johnson, R.K.; Norris, R.H. Setting expectations for the ecological condition of streams: The concept of reference condition. Ecol. Appl. 2006, 16, 1267–1276. [Google Scholar] [CrossRef] [PubMed]
  74. Baptista, D.F.; Buss, D.F.; Egler, M.; Giovanelli, A.; Silveira, M.P.; Nessimian, J.L. A multimetric index based on benthic macroinvertebrates for evaluation of Atlantic Forest streams at Rio de Janeiro State, Brazil. Hydrobiologia 2007, 575, 83–94. [Google Scholar] [CrossRef]
  75. Hu, X.; Zuo, D.; Xu, Z.; Huang, Z.; Liu, B.; Han, Y.; Bi, Y. Response of macroinvertebrate community to water quality factors and aquatic ecosystem health assessment in a typical river in Beijing, China. Environ. Res. 2022, 212, 113474. [Google Scholar] [CrossRef] [PubMed]
  76. Karr, J.R.; Chu, E.W. Restoring Life in Running Waters: Better Biological Monitoring; Island Press: Washington, DC, USA, 1999. [Google Scholar]
  77. Clements, W.H. Metal Tolerance and Predator-Prey Interactions in Benthic Macroinvertebrate Stream Communities. Ecol. Appl. 1999, 9, 1073–1084. [Google Scholar]
  78. Jones, F.C. Taxonomic sufficiency: The influence of taxonomic resolution on freshwater bioassessments using benthic macroinvertebrates. Environ. Rev. 2008, 16, 45–69. [Google Scholar] [CrossRef]
  79. Helson, J.E.; Williams, D.D. Development of a macroinvertebrate multimetric index for the assessment of low-land streams in the neotropics. Ecol. Indic. 2013, 29, 167–178. [Google Scholar] [CrossRef]
  80. Suhaila, A.H.; Che Salmah, M.R. Application of aquatic insects (Ephemeroptera, Plecoptera and Trichoptera) in water quality assessment of Malaysian headwater. Trop. Life Sci. Res. 2017, 28, 143–162. [Google Scholar]
  81. Bowles, D.E.; Cribbs, J.T.; Hinsey, J.A.; Williams, J.M. Aquatic Invertebrate Monitoring at Wilson’s Creek National Battlefield: 1996–2016 Status Report; National Park Service: Fort Collins, CO, USA, 2020; pp. 1–21. [Google Scholar]
  82. Lee, D.Y.; Lee, D.S.; Hwang, S.J.; Lee, K.L.; Park, Y.S. Distribution patterns and vulnerability of stoneflies (Plecoptera: Insecta) in South Korean streams with conservation perspectives. Glob. Ecol. Conserv. 2022, 34, e02030. [Google Scholar] [CrossRef]
  83. Alemneh, T.; Ambelu, A.; Zaitchik, B.F.; Bahrndorff, S.; Mereta, S.T.; Pertoldi, C. A macroinvertebrate multi-metric index for Ethiopian highland streams. Hydrobiologia 2019, 843, 125–141. [Google Scholar] [CrossRef]
  84. Arimoro, F.O.; Muller, W.J. Mayfly (Insecta: Ephemeroptera) community structure as an indicator of the ecological status of a stream in the Niger Delta area of Nigeria. Environ. Monit. Assess. 2010, 166, 581–594. [Google Scholar] [CrossRef] [PubMed]
  85. Morse, J.C.; Frandsen, P.B.; Graf, W.; Thomas, J.A. Diversity and Ecosystem Services of Trichoptera. Insects 2019, 10, 125. [Google Scholar] [CrossRef] [PubMed]
  86. Edegbene, A.O.; Akamagwuna, F.C.; Odume, O.N.; Arimoro, F.O.; Edegbene Ovie, T.T.; Akumabor, E.C.; Ogidiaka, E.; Kaine, E.A.; Nwaka, K.H. A Macroinvertebrate-Based Multimetric Index for Assessing Ecological Condition of Forested Stream Sites Draining Nigerian Urbanizing Landscapes. Sustainability 2022, 14, 11289. [Google Scholar] [CrossRef]
  87. White, J.C.; Hill, M.J.; Bickerton, M.A.; Wood, P.J. Macroinvertebrate taxonomic and trait compositions within lotic habitats affected by river restoration practices. Environ. Manag. 2017, 60, 513–525. [Google Scholar] [CrossRef] [PubMed]
  88. Arimoro, F.O.; Odume, O.N.; Uhunoma, S.I.; Edegbene, A.O. Anthropogenic impact on water chemistry and benthic macroinvertebrate associated changes in a southern Nigeria stream. Environ. Monit. Assess. 2015, 187, 14. [Google Scholar] [CrossRef]
  89. Odume, O.N.; Palmer, C.G.; Arimoro, F.O.; Mensah, P.K. Chironomid assemblage structure and morphological response to pollution in an effluent-impacted river, Eastern Cape, South Africa. Ecol. Indic. 2016, 67, 391–402. [Google Scholar] [CrossRef]
  90. Harrison, E.T.; Norris, R.H.; Wilkinson, S.N. The impact of fine sediment accumulation on benthic macroinvertebrates: Implications for river management. In Proceedings of the 5th Australian Stream Management Conference. Australian Rivers: Making a Difference, Charles Sturt University, New South Wales, Australia, 21–25 May 2007; Wilson, A.L., Dehaan, R.L., Watts, R.J., Page, K.J., Bowmer, K.H., Curtis, A., Eds.; Charles Sturt University, Institute for Land, Water and Society: New South Wales, Australia, 2007; pp. 139–144. [Google Scholar]
  91. Park, Y.S.; Lek, S.; Chon, T.S.; Verdonschot, P.F.M. Evaluation of Environmental Factors to Determine the Distribution of Functional Feeding Groups of Benthic Macroinvertebrates Using an Artificial Neural Network. J. Ecol. Field Biol. 2008, 31, 233–241. [Google Scholar] [CrossRef] [Green Version]
  92. Addo-Bediako, A. Spatial distribution patterns of benthic macroinvertebrate functional feeding groups in two rivers of the olifants river system, South Africa. J. Freshw. Ecol. 2021, 36, 97–109. [Google Scholar] [CrossRef]
  93. Purcell, A.H.; Bressler, D.W.; Paul, M.J.; Barbour, M.T.; Rankin, E.T.; Carter, J.L.; Resh, V.H. Assessment Tools for Urban Catchments: Developing Biological Indicators Based on Benthic Macroinvertebrates. JAWRA 2009, 45, 306–319. [Google Scholar] [CrossRef]
  94. Uttaruk, P.; Voharndee, P.; Jusanit, P.; Bunyaadunyakit, P.; Suriya, S.; Jaengjaithum, T.; Nittaya, K.; Sangpradub, N. Development of biotic index based on rapid bioassessment approaches using benthic macroinvertebrates for Chi and Mun headwater streams, Northeast Thailand. AES Bioflux 2011, 3, 29–43. [Google Scholar]
  95. Ramírez, A.; Gutiérrez-Fonseca, P.E. Functional feeding groups of aquatic insect families in Latin America: A critical analysis and review of existing literature. Rev. Biol. Trop. 2014, 62, 155–167. [Google Scholar] [CrossRef] [PubMed]
  96. Gething, K.J.; Ripley, M.C.; Mathers, K.L.; Chadd, R.P.; Wood, P.J. The influence of substrate type on macroinvertebrate assemblages within agricultural drainage ditches. Hydrobiologia 2020, 847, 4273–4284. [Google Scholar] [CrossRef]
  97. Leung, A.S.L.; Li, A.O.Y.; Dudgeon, D. Scales of spatiotemporal variation in macroinvertebrate assemblage structure in monsoonal streams: The importance of season. Freshw. Biol. 2012, 57, 218–231. [Google Scholar] [CrossRef]
  98. Marshall, J.C.; Steward, A.L.; Harch, B.D. Taxonomic resolution and quantification of freshwater macroinvertebrate samples from an Australian dryland river: The benefits and costs of using species abundance data. Hydrobiologia 2006, 572, 171–194. [Google Scholar] [CrossRef]
Figure 1. Map showing the distribution of the 10 study localities (wadeable rivers and streams) in Lao PDR. Photographs of 10 study localities. (A) Nam Houng upstream (HU), (B) Nam Houng downstream (HL), (C) Nam Khod (KD), (D) Nam Khan (KN), (E) Nam Song (NS), (F) Nam Thang (NT), (G) Nam Hinboun upstream (BU), (H) Nam Hinboun downstream (BL), (I) Xe Katam (KT), (J) Xe Namnoy (NN).
Figure 1. Map showing the distribution of the 10 study localities (wadeable rivers and streams) in Lao PDR. Photographs of 10 study localities. (A) Nam Houng upstream (HU), (B) Nam Houng downstream (HL), (C) Nam Khod (KD), (D) Nam Khan (KN), (E) Nam Song (NS), (F) Nam Thang (NT), (G) Nam Hinboun upstream (BU), (H) Nam Hinboun downstream (BL), (I) Xe Katam (KT), (J) Xe Namnoy (NN).
Water 15 00625 g001
Figure 2. Benthic macroinvertebrates found in the study localities with their tolerance values. A-F Diptera: (A) Athericidae, Suragina, tolerance value three; (B) Blephariceridae, Blepharicera, tolerance value one; (C) Ceratopogonidae, Bezzia, tolerance value seven; (D) Chironomidae, tolerance value nine; (E) Culicidae, tolerance value nine; (F) Dixidae, tolerance value three: (GL) Ephemeroptera: (G) Caenidae, Caenis, tolerance value seven; (H) Caenidae, Clypeocaenis, tolerance value seven; (I) Heptageniidae, Thalerosphyrus lamuriensis, tolerance value two; (J) Potamanthidae, Potamanthus formosus, tolerance value two; (K) Prosopistomatidae, Prosopistoma, tolerance value three; (L) Vietnamellidae, Vietnamella thani, tolerance value one: (MQ) Plecoptera: (M) Nemouridae, Indonemoura, tolerance value one; (N) Peltoperlidae, Cryptoperla, tolerance value one; (O) Perlidae, Etrocorema, tolerance value three; (P) Perlidae, Phanoperla, tolerance value two; (Q) Perlidae, Togoperla, tolerance value three: (RV) Trichoptera: (R) Calamoceratidae, Anisocentropus, tolerance value one; (S) Calamoceratidae, Ganonema, tolerance value one; (T) Glossosomatidae, Agapetus, tolerance value one; (U) Hydropsychidae, Cheumatopsyche, tolerance value seven; (V) Stenopsychidae, Stenopsyche, tolerance value three: (WX) Mesogastropoda: (W) Planorbidae, Indoplanorbis exustus, tolerance value eight; (X) Thiaridae, Thiara scabra, tolerance value eight: (Y) Veneroida, Corbiculidae, Corbicula, tolerance value eight. Scale bars: 2 mm.
Figure 2. Benthic macroinvertebrates found in the study localities with their tolerance values. A-F Diptera: (A) Athericidae, Suragina, tolerance value three; (B) Blephariceridae, Blepharicera, tolerance value one; (C) Ceratopogonidae, Bezzia, tolerance value seven; (D) Chironomidae, tolerance value nine; (E) Culicidae, tolerance value nine; (F) Dixidae, tolerance value three: (GL) Ephemeroptera: (G) Caenidae, Caenis, tolerance value seven; (H) Caenidae, Clypeocaenis, tolerance value seven; (I) Heptageniidae, Thalerosphyrus lamuriensis, tolerance value two; (J) Potamanthidae, Potamanthus formosus, tolerance value two; (K) Prosopistomatidae, Prosopistoma, tolerance value three; (L) Vietnamellidae, Vietnamella thani, tolerance value one: (MQ) Plecoptera: (M) Nemouridae, Indonemoura, tolerance value one; (N) Peltoperlidae, Cryptoperla, tolerance value one; (O) Perlidae, Etrocorema, tolerance value three; (P) Perlidae, Phanoperla, tolerance value two; (Q) Perlidae, Togoperla, tolerance value three: (RV) Trichoptera: (R) Calamoceratidae, Anisocentropus, tolerance value one; (S) Calamoceratidae, Ganonema, tolerance value one; (T) Glossosomatidae, Agapetus, tolerance value one; (U) Hydropsychidae, Cheumatopsyche, tolerance value seven; (V) Stenopsychidae, Stenopsyche, tolerance value three: (WX) Mesogastropoda: (W) Planorbidae, Indoplanorbis exustus, tolerance value eight; (X) Thiaridae, Thiara scabra, tolerance value eight: (Y) Veneroida, Corbiculidae, Corbicula, tolerance value eight. Scale bars: 2 mm.
Water 15 00625 g002
Figure 3. CCA ordination analysis showing the association between environmental variables and benthic macroinvertebrate abundance in the reference, intermediate, and stressed sites. Black color denotes intolerant taxa; yellow denotes moderate and tolerant taxa.
Figure 3. CCA ordination analysis showing the association between environmental variables and benthic macroinvertebrate abundance in the reference, intermediate, and stressed sites. Black color denotes intolerant taxa; yellow denotes moderate and tolerant taxa.
Water 15 00625 g003
Figure 4. Box-and-whisker plots for each of the 11 core metrics sensitive to discrimination between reference and stressed sites for MMI development. (A) Total taxa, (B) EPT taxa, (C) Ephemeroptera taxa, (D) %Diptera, (E) %Plecoptera, (F) %Tolerant, (G) BeckBI, (H) %Intolerant, (I) Filterer taxa, (J) %Sprawlers, (K) %Burrowers.
Figure 4. Box-and-whisker plots for each of the 11 core metrics sensitive to discrimination between reference and stressed sites for MMI development. (A) Total taxa, (B) EPT taxa, (C) Ephemeroptera taxa, (D) %Diptera, (E) %Plecoptera, (F) %Tolerant, (G) BeckBI, (H) %Intolerant, (I) Filterer taxa, (J) %Sprawlers, (K) %Burrowers.
Water 15 00625 g004
Figure 5. Box-and-whisker plot of the final index scores for Lao MMI used to discriminate between reference, intermediate, and stressed sites.
Figure 5. Box-and-whisker plot of the final index scores for Lao MMI used to discriminate between reference, intermediate, and stressed sites.
Water 15 00625 g005
Table 1. General characteristics and geographical position of 10 localities of wadeable rivers and streams conduced in this study in Lao PDR.
Table 1. General characteristics and geographical position of 10 localities of wadeable rivers and streams conduced in this study in Lao PDR.
EcoregionStudy LocalitiesType of Lotic EcosystemsProvinceCoordinatesAltitude (masl)ActivitySubstrate Types (%)
The Lower Lancang1. Nam Houng upstream (HU)Wadeable riverXaignabouly19°18′30.68″ N
101°43′24.94″ E
294Forest areas, small agriculture patchesCobble (60%), pebble (10%), gravel (20%), sand (10%)
2. Nam Houng downstream (HL)Wadeable riverXaignabouly19°15′38.89″ N
101°42′46.01″ E
283Residential areas, agriculture areasCobble (70%), gravel (10%), sand (15%), detritus (5%)
3. Nam Khod (KD)StreamLuangprabang19°44′14.12″ N
102°8′52.63″ E
375Forest areas, small agriculture patchesBoulder (20%), cobble (30%), pebble (15%), gravel (20%), sand (5%), detritus (10%)
4. Nam Khan (NK)StreamLuangprabang19°43′55.58″ N
102°9′24.31″ E
332Agriculture areas, intensive erosionCobble (50%), detritus (10%), muck-mud (40%)
Khorat Plateau5. Nam Song (NS)Wadeable riverVientiane19°6′13.07″ N
102°30′3.16″ E
315Forest areas, small agriculture patchesBoulder (10%), cobble (65%), pebble (15%), gravel (5%), sand (5%)
6. Nam Thang (NT)StreamVientiane19°6′15.12″ N
102°29′54.98″ E
324Forest areas, small agriculture patchesCobble (60%), pebble (15%), gravel (10%), sand (10%), detritus (5%)
7. Nam Hinboun upstream (BU)Wadeable riverKhammouan17°57′19.47″ N
104°45′28.73″ E
164Preserve areasCobble (15%), pebble (5%), gravel (60%), sand (20%)
8. Nam Hinboun downstream (BL)Wadeable riverKhammouan17°57′33.58″ N
104°43′32.96″ E
156Agriculture areaGravel (20%), sand (40%), muck-mud (40%)
Kratie-Stung Treng9. Xe Katam (KT)Wadeable riverChampasak15°7′48.93″ N
106°40′09.92″ E
257Forest areas, small agriculture patchesBoulder (60%), cobble (30%), gravel (5%), sand (5%)
10. Xe Namnoy (NN)Wadeable riverXekong15°13′37.99″ N
106°44′45.30″ E
132Forest areas, small agriculture patchesBoulder (50%), cobble (30%), pebble (10%), gravel (5%), sand (5%)
Table 2. Environmental variables of reference, intermediate, and stressed sites. Data are expressed as Mean ± SD.
Table 2. Environmental variables of reference, intermediate, and stressed sites. Data are expressed as Mean ± SD.
ParametersReference Group
(n = 7)
Intermediate Group
(n = 26)
Stressed Group
(n = 7)
p-Value
Air temperature (°C)24.64 ± 4.5327.12 ± 4.8630.64 ± 4.190.152
Water temperature (°C)22.59 ± 2.1624.73 ± 2.6925.91 ± 3.500.103
Water channel width (m)19.85 ± 11.9417.62 ± 10.2014.18 ± 11.040.592
Water depth (cm)21.35 ± 7.3130.82 ± 13.9933.43 ± 14.020.182
Water velocity (m/s)0.49 ± 0.180.43 ± 0.180.32 ± 0.210.249
Turbidity (NTU)2.83 ± 1.842.19 ± 1.932.56 ± 2.430.772
Suspended solids (mg/L)5.38 ± 2.966.19 ± 3.986.05 ± 2.160.743
Electrical conductivity (µS/cm)113.33 ± 65.56 a256.40 ± 124.29 b379.48 ± 101.02 b0.001 *
Total dissolved solids (mg/L)58.31 ± 34.23 a132.09 ± 64.03 b194.57 ± 53.45 b0.001 *
Dissolved oxygen (mg/L)6.30 ± 0.166.63 ± 0.756.62 ± 1.780.325
pH7.74 ± 0.487.91 ± 0.367.85 ± 0.460.665
Nitrate nitrogen (mg/L)0.20 ± 0.05 a0.30 ± 0.12 b0.32 ± 0.10 b0.021 *
Orthophosphate (mg/L)0.17 ± 0.060.19 ± 0.100.24 ± 0.130.647
Biochemical oxygen demand (mg/L)1.45 ± 0.681.14 ± 0.811.11 ± 0.840.573
Chlorophyll a (μg/L)0.83 ± 0.200.97 ± 0.431.15 ± 0.820.786
THS152.57 ± 14.33 a137.27 ± 17.40 a80.71 ± 21.78 b0.000 *
Within each row of means different letters (a, b) indicate significant differences at p < 0.05 according to pairwise comparisons of the Kruskal–Wallis test and LSD multiple comparisons of one-way ANOVA. * Significant difference at p < 0.05.
Table 3. Numbers of benthic macroinvertebrate taxa found in reference, intermediate, and stressed sites. Data are expressed as Mean ± SD.
Table 3. Numbers of benthic macroinvertebrate taxa found in reference, intermediate, and stressed sites. Data are expressed as Mean ± SD.
TaxaReference GroupIntermediate GroupStressed Groupp-Value
ANNELIDA
Oligochaeta0.14 ± 0.380.08 ± 0.270.29 ± 0.490.338
ARTHROPODA
Decapoda0.43 ± 0.530.96 ± 0.961.14 ± 0.690.247
Isopoda0.29 ± 0.490.04 ± 0.2000.330
Coleoptera8.14 ± 4.986.85 ± 2.944.14 ± 1.680.085
Collembola00.04 ± 0.200-
Diptera3.71 ± 1.114.92 ± 2.334.57 ± 0.980.395
Ephemeroptera14.29 ± 1.80 a13.15 ± 2.71 ab10.29 ± 3.20 b0.038 *
Hemiptera4.71 ± 2.434.96 ± 2.223.00 ± 2.240.109
Lepidoptera0.71 ± 0.761.15 ± 1.050.57 ± 0.980.298
Megaloptera0.71 ± 0.490.50 ± 0.510.43 ± 0.530.522
Odonata5.14 ± 3.394.96 ± 2.244.00 ± 2.310.593
Orthoptera0.29 ± 0.490.15 ± 0.3700.620
Plecoptera2.14 ± 0.90 a1.92 ± 1.74 ab0.71 ± 0.49 b0.029 *
Trichoptera10.57 ± 2.94 a7.85 ± 2.91 ab4.86 ± 2.48 b0.006 *
MOLLUSCA
Mesogastropoda0.57 ± 0.531.19 ± 0.942.00 ± 1.410.072
Unionoida00.04 ± 0.200.29 ± 0.490.330
Veneroida00.31 ± 0.470.71 ± 0.490.109
Total taxa51.86 ± 15.20 a49.08 ± 7.56 a37.00 ± 9.00 b0.009 *
Within each row of means different letters (a, b) indicate significant differences at p < 0.05 according to pairwise comparisons of the Kruskal–Wallis test and LSD multiple comparisons of one-way ANOVA. * Significant difference at p < 0.05.
Table 4. Relative abundance of benthic macroinvertebrates found in reference, intermediate, and stressed sites. Data are expressed as Mean ± SD.
Table 4. Relative abundance of benthic macroinvertebrates found in reference, intermediate, and stressed sites. Data are expressed as Mean ± SD.
TaxaReference GroupIntermediate GroupStressed Groupp-Value
ANNELIDA
Oligochaeta0.29 ± 0.760.12 ± 0.431.71 ± 4.110.320
ARTHROPODA
Decapoda0.43 ± 0.533.00 ± 4.855.86 ± 8.650.089
Isopoda0.43 ± 0.790.12 ± 0.5900.352
Coleoptera35.29 ± 35.0034.46 ± 27.8725.00 ± 35.420.639
Collembola00.04 ± 0.200-
Diptera32.00 ± 32.34 a38.50 ± 23.03 a78.57 ± 61.95 b0.045 *
Ephemeroptera126.57 ± 62.66122.92 ± 51.9096.57 ± 48.590.539
Hemiptera19.14 ± 15.7020.73 ± 13.5318.57 ± 17.930.801
Lepidoptera5.00 ± 10.652.23 ± 2.691.14 ± 2.040.489
Megaloptera2.00 ± 2.161.42 ± 2.270.43 ± 0.530.311
Odonata12.71 ± 8.6017.69 ± 15.2926.43 ± 23.960.352
Orthoptera0.29 ± 0.490.23 ± 0.6500.651
Plecoptera11.71 ± 11.60 a13.81 ± 20.54 a1.00 ± 0.82 b0.020 *
Trichoptera84.29 ± 28.1580.15 ± 51.4462.43 ± 70.960.450
MOLLUSCA
Mesogastropoda10.71 ± 25.747.38 ± 9.749.43 ± 12.660.374
Unionoida00.08 ± 0.391.86 ± 3.290.308
Veneroida00.73 ± 1.284.00 ± 4.550.060
Total individuals340.86 ± 19.99343.62 ± 16.11333.00 ± 17.300.168
Chironomidae17.57 ± 21.83 a20.19 ± 10.68 b63.14 ± 68.04 c0.005 *
Within each row of means different letters (a, b, c) indicate significant differences at p < 0.05 according to pairwise comparisons of the Kruskal Wallis test. * Significant difference at p < 0.05.
Table 5. Thirty-five candidate metrics selected for the development of MMI and their predicted response to disturbance or pollution, discrimination power score, and DE value. Metrics that met or did not meet the test criteria are expressed as symbols (/) and (-), respectively.
Table 5. Thirty-five candidate metrics selected for the development of MMI and their predicted response to disturbance or pollution, discrimination power score, and DE value. Metrics that met or did not meet the test criteria are expressed as symbols (/) and (-), respectively.
MetricExpected Response of Metrics to PollutionDiscrimination Power ScoreDE ValueMetric Selection
Richness category
1. Total taxaDecrease271.43/
2. EPT taxaDecrease3100.00/
3. Ephemeroptera taxaDecrease371.43/
4. Plecoptera taxaDecrease3100.00-
5. Trichoptera taxaDecrease385.71Redundant
6. Coleoptera taxaDecrease142.86-
7. EPTC taxaDecrease3100.00Redundant
Composition category
8. %EPTDecrease157.14-
9. %EphemeropteraDecrease142.86-
10. Margalef’s indexDecrease271.43-
11. %OdonataIncrease257.14-
12. %ChironomidaeIncrease257.14-
13. %DipteraIncrease371.43/
14. %PlecopteraDecrease3100.00/
15. %TrichopteraDecrease257.14-
Tolerance value category
16. Intolerant taxaDecrease3100.00Redundant
17. %TolerantIncrease385.71/
18. %Dominant taxonIncrease371.43-
19. Beck’s Biotic IndexDecrease3100.00/
20. Simpson IndexIncrease271.43-
21. Hilsenhof’s Biotic IndexIncrease3100.00Redundant
22. %IntolerantDecrease3100.00/
Functional feeding group category
23. %FilterersDecrease157.14-
24. %ScrapersDecrease371.43-
25. %CollectorsIncrease257.14-
26. Collector taxaDecrease157.14-
27. Filterer taxaDecrease3100.00/
28. Predator taxaDecrease242.86-
29. Scraper taxaDecrease357.14-
30. %ShreddersIncrease142.86-
31. Shredder taxaDecrease371.43-
Habit category
32. Clinger taxaDecrease3100.00Redundant
33. %ClingersDecrease271.43-
34. %SprawlersIncrease385.71/
35. %BurrowersIncrease385.71/
Table 6. Score thresholds for each of the 11 core metrics integrated in Lao MMI.
Table 6. Score thresholds for each of the 11 core metrics integrated in Lao MMI.
MetricExpected Response to PollutionStatistic Value of Reference SitesScoring Criteria
Min.25th50th75thMax.531
Richness category
Total taxaDecrease3737405463≥3736–19<19
EPT taxaDecrease2424253032≥2423–12<12
Ephemeroptera taxaDecrease1113151616≥1312–7<7
Composition category
%DipteraIncrease0.61.96.314.426.1≤14.414.5–21.6>21.6
%PlecopteraDecrease0.30.62.35.211.1≥0.60.5–0.3<0.3
Tolerance value category
%TolerantIncrease2.34.49.62525.8≤2525.1–37.5>37.5
Beck’s Biotic IndexDecrease1214192325≥1413–7<7
%IntolerantDecrease3.1619.923.524.4≥65.9–3.0<3.0
Functional feeding category
Filterers taxaDecrease6681111≥65.0–3.0<3
Habit category
%SprawlersIncrease3.94.96.7825.3≤8.08.1–12.0>12.0
%BurrowersIncrease0.91.310.519.722≤19.719.8–28.5>28.5
Table 7. Proposed water quality classes and index score for the studied wadeable rivers and streams in Lao PDR.
Table 7. Proposed water quality classes and index score for the studied wadeable rivers and streams in Lao PDR.
Stream Quality ClassPercentileIndex Score
Excellent≥75th53
Good≥25th49–52
Fair<25th25–48
Poor-13–24
Impaired-<13
Table 8. Water quality of 40 sampling sites (wadeable rivers and streams) in Lao PDR assessed by using surface water quality standards of Lao PDR (physicochemical method) and the developed Lao MMI.
Table 8. Water quality of 40 sampling sites (wadeable rivers and streams) in Lao PDR assessed by using surface water quality standards of Lao PDR (physicochemical method) and the developed Lao MMI.
Sampling SiteConditionThe Surface Water Quality Standard of Lao PDRLao MMI
Stream ClassWater QualityIndex ScoreWater Quality
Cool season 2016
HUC16Reference2Good49Good
HLC16Stressed2Good33Fair
KDC16Intermediate2Good47Fair
KNC16Stressed2Good31Fair
NSC16Intermediate2Good53Excellent
NTC16Reference2Good55Excellent
BUC16Intermediate2Good45Fair
BLC16Intermediate2Good31Fair
KTC16Reference2Good53Excellent
NNC16Reference2Good53Excellent
Hot season 2017
HUH17Intermediate2Good35Fair
HLH17Intermediate2Good39Fair
KDH17Intermediate2Good53Excellent
KNH17Stressed2Good37Fair
NSH17Intermediate1Excellent43Fair
NTH17Intermediate2Good53Excellent
BUH17Intermediate2Good27Fair
BLH17Stressed2Good17Poor
KTH17Intermediate2Good53Excellent
NNH17Intermediate2Good53Excellent
Cool season 2017
HUC17Reference2Good43Fair
HLC17Intermediate2Good53Excellent
KDC17Intermediate2Good45Fair
KNC17Intermediate2Good39Fair
NSC17Intermediate2Good45Fair
NTC17Reference2Good51Good
BUC17Intermediate2Good33Fair
BLC17Stressed2Good37Fair
KTC17Intermediate2Good53Excellent
NNC17Intermediate2Good51Good
Hot season 2018
HUH18Intermediate2Good39Fair
HLH18Intermediate2Good41Fair
KDH18Intermediate2Good43Fair
KNH18Stressed3Fair35Fair
NSH18Intermediate2Good47Fair
NTH18Intermediate2Good45Fair
BUH18Intermediate2Good47Fair
BLH18Stressed2Good19Poor
KTH18Reference2Good51Good
NNH18Intermediate2Good47Fair
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sripanya, J.; Vongsombath, C.; Vannachak, V.; Rattanachan, K.; Hanjavanit, C.; Mahakham, W.; Sangpradub, N. Benthic Macroinvertebrate Communities in Wadeable Rivers and Streams of Lao PDR as a Useful Tool for Biomonitoring Water Quality: A Multimetric Index Approach. Water 2023, 15, 625. https://doi.org/10.3390/w15040625

AMA Style

Sripanya J, Vongsombath C, Vannachak V, Rattanachan K, Hanjavanit C, Mahakham W, Sangpradub N. Benthic Macroinvertebrate Communities in Wadeable Rivers and Streams of Lao PDR as a Useful Tool for Biomonitoring Water Quality: A Multimetric Index Approach. Water. 2023; 15(4):625. https://doi.org/10.3390/w15040625

Chicago/Turabian Style

Sripanya, Jutamas, Chanda Vongsombath, Viengkhone Vannachak, Kaewpawika Rattanachan, Chutima Hanjavanit, Wuttipong Mahakham, and Narumon Sangpradub. 2023. "Benthic Macroinvertebrate Communities in Wadeable Rivers and Streams of Lao PDR as a Useful Tool for Biomonitoring Water Quality: A Multimetric Index Approach" Water 15, no. 4: 625. https://doi.org/10.3390/w15040625

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop