Next Article in Journal
Extraction and Preference Ordering of Multireservoir Water Supply Rules in Dry Years
Next Article in Special Issue
Effects of the “Run-of-River” Hydro Scheme on Macroinvertebrate Communities and Habitat Conditions in a Mountain River of Northeastern China
Previous Article in Journal
Implications of Texture and Erodibility for Sediment Retention in Receiving Basins of Coastal Louisiana Diversions
Previous Article in Special Issue
Roles of N:P Ratios on Trophic Structures and Ecological Stream Health in Lotic Ecosystems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Distribution of Benthic Macroinvertebrate Assemblages in Relation to Environmental Variables in Korean Nationwide Streams

1
Department of Environmental Science, Konkuk University, Seoul 143-701, Korea
2
The National Institute of Environmental Research, Incheon 404-170, Korea
3
Department of Life and Nanopharmaceutical Sciences and Department of Biology, Kyung Hee University, Seoul 130-701, Korea
4
Department of Biology, Kyonggi University, Suwon 443-760, Korea
5
Department of Environmental Health Science, Konkuk University, Seoul 143-701, Korea
*
Author to whom correspondence should be addressed.
Water 2016, 8(1), 27; https://doi.org/10.3390/w8010027
Submission received: 23 October 2015 / Revised: 13 January 2016 / Accepted: 15 January 2016 / Published: 20 January 2016

Abstract

:
Conserving and enhancing freshwater biodiversity are global issues to ensure ecosystem integrity and sustainability. To meet this, it is critical to understand how the biological assemblages are determined by environmental gradients in different spatial scales. Nevertheless, information on their large-scale environmental relationships remains scarce in Korea. We aimed to understand nationwide spatial distribution patterns of benthic macroinvertebrates and important environmental factors affecting their distribution in 388 streams and rivers across Korea. A total of 340 taxa, belonging to 113 families in 23 orders of five phyla, were identified. Assemblage composition in most Korean streams included a few predominant colonizers and a majority of rare taxa. Cluster analysis based on benthic macroinvertebrates classified a total of 720 sampling sites into five clusters according to the pollution levels from fast-flowing less polluted streams with low electrical conductivity to moderately or severely polluted streams with high electrical conductivity and slow water velocity. Canonical correspondence analysis revealed that altitude, water velocity and streambed composition were the most important determinants, rather than watershed and water chemistry variables, for explaining the variation in macroinvertebrate assemblage patterns. The results provide basic information for establishing the conservation and restoration strategies of macroinvertebrate biodiversity against anthropogenic disturbances and developing more confident bio-assessment tools for diagnosing stream ecosystem integrity.

1. Introduction

Freshwater ecosystems occupying a very tiny fraction of the Earth’s surface support remarkable biodiversity and abundance of benthic macroinvertebrates [1,2]. However, growing human pressures have substantially deteriorated the overall ecological integrity and induced a biodiversity crisis. This problem now increasingly becomes a global issue, and thereby calls attention to establish relevant conservation strategies and wise management practices to maintain sustainable freshwater environments.
Although Asian riverine ecosystems contain remarkable taxonomic diversity as well as high levels of endangerment and endemism, studies on biodiversity and ecology of these ecosystems have been poor [1,3]. Year-to-year variations in seasonality and its relationship to monsoonal climate, particularly in temperate Asian regions, can be major drivers that profoundly affect the hydrologic regime and geomorphology in stream environments, consequently determining the distribution and abundance of macroinvertebrates [4]. For example, a wet season low and a dry season high are expected for periodic seasonal patterns in abundance, depending on the frequency and intensity of summer monsoon rainfall. However, habitat destruction and water quality degradation in Asian river systems have become more dramatically epidemic than any other continents due to the fast growing human population and demands for economic development [5,6]. Because failure of conservation efforts has mainly resulted from scant evidence on river ecology, the priority over scientific and practical challenges to overcome this circumstance must be put on establishing species inventories and ecological information based on the accumulation of region-specific case studies, particularly at watershed or national scales [4,7].
Macroinvertebrate diversity and abundance are significant community attributes that are controlled by a variety of mechanisms at different spatial scales. A number of studies have documented how macroinvertebrate assemblages respond to environmental variables and which variables best explain their distribution and abundance. Some studies showed good relationships among macroinvertebrate assemblages, chemical variables [8], and the organic energy base [9], whereas habitat-related physical factors were widely demonstrated as primary contributors such as substrate composition [10], flow and current velocity [11], elevation and stream size [12] and temperature [13]. Vegetation, geology and human land use are also important for their spatial distribution [14,15]. While there remains a large number of documents focusing on such environmental relationships at finer spatial scales, studies on large-scale spatial patterns have been relatively small (e.g., at national scale [16,17] and at continental scale [18,19]).
Studies on the spatial distribution patterns of macroinvertebrate assemblages based on their environmental relationships are crucial. Firstly, those studies provide fundamental information for the conservation and restoration of biodiversity against anthropogenic disturbances. Secondly, knowledge of the species response to environmental gradients is important to separate the effects of pollution from the effects of natural variables on community structures [20]. Thirdly, it is possible to develop more confident bio-assessment tools for diagnosing stream ecosystem integrity because macroinvertebrate responds sensitively to environmental alterations that lead to changes in composition and community structures [21,22]. Thus, we performed a synoptic study of the large scale spatial distribution of benthic macroinvertebrates in relation to environmental variables. Specifically, the objectives of this study were: (i) to characterize the spatial distribution and assemblage structures of macroinvertebrates; (ii) to identify environmental distinction of Korean stream ecosystems based on their assemblages; (iii) to determine major environmental variables that affect their distribution; and (iv) to provide methodological considerations to improve biomonitoring programs, particularly focusing on sampling procedures. A relevant scientific basis would then be established with these data for developing sustainable management, conservation practices, and a reliable biomonitoring program. The results of this study could also contribute to a better understanding of the spatial distribution of freshwater macroinvertebrates in poorly explored Asian regions.

2. Materials and Methods

2.1. Study Area

Samples were collected at 720 sampling sites from 388 streams and rivers on a nationwide scale in South Korea. South Korea is located between 37°00′ N latitude and 127°30′ E longitude, encompassing the southern half of the Korean Peninsula with an entire area of approximately 100,033 km2. The annual precipitation is 1308 mm, but there is substantial variation among seasons [23]: severe flooding events during the summer monsoon period, and base flow or even drought conditions in the other seasons.
The entire stream system throughout South Korea comprises five watersheds, including five large rivers, their tributaries, and other small independent streams (Figure 1). Land use types in each watershed generally display well-managed forests upstream, and agricultural and urban development from middle to lower regions. Most Korean streams have suffered from a variety of human activities that alter the physicochemical stream environments, particularly caused by channel modification and eutrophication [24].
Among 720 sampling sites, the largest number of sites belonged to the Han River watershed (HRW) (n = 320), followed by the Nakdong River watershed (NRW) (n = 130), the Geum River watershed (GRW) (n = 130), the Youngsan River watershed (YRW) (n = 76), and the Seomjin River watershed (SRW) (n = 64). It was assumed that such a nationwide survey covered the majority of stream types in Korea to understand how lotic macroinvertebrates are spatially distributed in relation to environmental factors. Field sampling was conducted during spring (May 2009) under base flow conditions because the highest benthic macroinvertebrate diversity was expected at that time. More information of the five major river watersheds in Korea can be found in Hwang et al. [25].
Figure 1. Geographic locations of the study sites. Five macroinvertebrate-based site groups were classified as G1a (n = 126, black circles), G1b (n = 199, white circles), G2a (n = 106, white triangles), G2b (n = 118, gray triangles), and G2c (n = 171, black triangles) based on Sørenson distance measure cluster analysis. Dark lines indicate the five large rivers and the light lines display their tributaries and small independent streams.
Figure 1. Geographic locations of the study sites. Five macroinvertebrate-based site groups were classified as G1a (n = 126, black circles), G1b (n = 199, white circles), G2a (n = 106, white triangles), G2b (n = 118, gray triangles), and G2c (n = 171, black triangles) based on Sørenson distance measure cluster analysis. Dark lines indicate the five large rivers and the light lines display their tributaries and small independent streams.
Water 08 00027 g001

2.2. Measurement of Environmental Variables

Environmental variables were measured both in the field and in the laboratory to define their effects on benthic macroinvertebrate assemblages. Three categories of environmental variables were considered such as watershed-related regional variables, physical in-stream properties, and water quality elements.
Altitude and land-use types were considered as regional variables for watershed characteristics. The altitude at each sampling site was obtained using digital elevation data (Openmate Inc., Seoul, Korea). The proportion of prevalent land use types was determined for each sampling site using a topographic map (1:50,000).
Physical in-stream properties were measured during field surveys: (i) stream width for the distance from bank to bank at a transect representative of the stream channel; (ii) wetted stream width using a range finder (model LRM-1500M, Newcon Optik Inc., Toronto, ON, Canada); (iii) water depth of the vertical distance from the water surface to stream bottom; (iv) current velocity at riffles or gliding runs using a current meter (3000-LX, Swoffer Instruments, Inc., Tukwila, WA, USA), or calculated by the Craig method [26,27]; (v) percent substrate composition visually estimated as fine (<2 mm) and coarse-sized particles (≥2 mm).
Water quality variables such as pH, water temperature, dissolved oxygen (DO), turbidity, and electrical conductivity (EC) were measured with a multi-probe portable meter (e.g., YSI 6920, YSI Inc., Yellow Springs, OH, USA or Horiba U-22XD, Kyoto, Japan) at the center of each sampling stretch. Two liters of water samples were collected in sterilized plastic bottles at each site, kept in a container with ice, and transported to the laboratory. Laboratory measurements were conducted to determine the biochemical oxygen demand (BOD), total nitrogen (TN), and total phosphorus (TP) following Standard methods [28].

2.3. Sampling of Benthic Macroinvertebrate

A Surber sampler (30 cm × 30 cm, 1 mm mesh) was employed to collect benthic macroinvertebrates. Samplings were conducted at fast-flowing riffle or gliding run habitats (in the case when suitable riffles were not available) within 100 m. Quantitative samples were taken from three randomly selected riffles at each site, pooled together in a 500 mL plastic bottle with 80% ethanol, and labeled. The sampling device and procedures followed the guidelines of the National Aquatic Ecological Monitoring Program (NAEMP), Korea [27].
After field sampling, all organisms were hand-sorted from detritus and inorganic material, and stored in 70% ethanol. Subsampling was permitted only for dominant taxonomic groups (e.g., Oligochaeta, Ephemerellidae, Chironomidae, and Hydropsychidae) with large numbers of specimens available in each sample. Each macroinvertebrate specimen was identified to the lowest possible taxonomic level (usually genus or species). However, several taxa with limited systematic information (e.g., Coleoptera and Diptera) were identified only to the family level. All individuals were counted and converted to individuals/m2.

2.4. Data Analysis

Differences in environmental variables and macroinvertebrate assemblages were identified among the five major river watersheds by the Kruskal–Wallis test, followed by Dunn’s nonparametric multiple comparison test if they are significantly different (p < 0.05). Spearman’s rank correlation analysis was used to indicate the relationships between the most dominant taxa and environmental variables. A cluster analysis was conducted to classify the benthic macroinvertebrate communities using the flexible beta method (beta = −0.25) with the Sørenson distance measure [29]. Samples were classified into clusters based on similarities in the community composition. The multi-response permutation procedure (MRPP, [30]), which is a non-parametric method to test for differences in assemblage structure among a priori defined groups [31], was conducted to evaluate differences among clusters. Differences in environmental variables among clusters were evaluated using the Kruskal–Wallis test and Dunn’s nonparametric multiple comparison test. The Kruskal–Wallis test, Dunn’s test and Spearman correlation analysis were performed with STATISTICA software (StatSoft, Inc., version 7, Tulsa, OK, USA).
Indicator species analysis (IndVal, [32]) was performed to evaluate potential indicator species in each cluster defined in advance. The indicator value of a species is the product of its relative abundance and frequency (× 100), ranging from 0 (no indication) to 100 (perfect indication) [33]. A perfect indicator of a particular group should be faithful and exclusive to that group, never occurring in other groups [29]. A species in a cluster that had an indicator value greater than in any other cluster was defined as good indicator species for that cluster in this study. A Monte Carlo method was performed to test the significance of indicator values.
Canonical correspondence analysis (CCA) was used to relate macroinvertebrate assemblages to environmental variables and to identify which environmental variables could best differentiate among the clusters [34]. Monte-Carlo simulations were carried out to verify whether variables exerted a significant effect (p < 0.05) on macroinvertebrate distributions. For the cluster analysis and CCA, those taxa with less than 0.2% of total abundance were excluded to minimize the effects of rare taxa. Taxa abundance data were transformed to log (x + 1) in both analyses to down-weight the effects of dominant taxa. Square-root-transformation was used for the environmental parameters expressed in percentages (e.g., type of land use and substrate composition), and the other variables were transformed to log (x + 1) except pH. After these processes all environmental variables were rescaled in the range of 0 and 1 based on the minimum-maximum range normalization [25]. Spearman correlation coefficients between scores and environmental variables were calculated to assist interpretation of changes in community profile using STATISTICA software (StatSoft, Inc., version 7, Tulsa, OK, USA). Cluster analysis, IndVal, MRPP and CCA were conducted using PC-ORD software (ver. 4.25, MjM Software Design, Gleneden Beach, OR, USA) [34].

3. Results

3.1. Environmental Characteristics

Despite a large variation, physicochemical factors were significantly different among watersheds (Table 1). On average, the altitude, percent forest, and water velocity of the HRW and SRW were the higher than the other watersheds. The SRW streambeds were very heterogeneous with the highest percentage of coarse-sized particles, whereas those of the YRW contained a large amount of fine sediment. High altitude, fast water velocity and substrate complexity in both the HRW and SRW indicated good in-stream habitat conditions and potential to support high biodiversity and abundance of macroinvertebrates (Table 1a). The other three watersheds (NRW, GRW and YRW) were characterized by low altitude and a high proportion of agricultural land use.
Water quality parameters varied more prominently than physical variables among river watersheds (Table 1b). The average concentrations of DO, BOD, TN, TP, and turbidity were generally lower in the SRW than those in the other watersheds. EC varied widely depending on the sampling location, with the highest values near estuaries. Nutrient-related factors were particularly high in the GRW, YRW, and HRW.

3.2. Benthic Macroinvertebrate Assemblages

A total of 340 taxa, belonging to 113 families in 23 orders of five phyla, were identified in the five major river watersheds during the survey. Most of these were aquatic insects (272 species) including 144 Ephemeroptera Plecoptera Trichoptera (EPT; 62 mayflies, 24 stoneflies and 58 caddisflies) taxa and 35 Dipteran species. Taxa richness from the 720 sampling sites ranged from 0 to 49 with a mean of 14.4 (± 9.1) with the highest diversity at HRW and the lowest in YRW (Table 2a). In contrast, homogeneous streambeds and nutrient enrichment attributed to relatively low Shannon diversity index and high dominance index in both the NRW and YRW. The proportion of EPT taxa richness ranged from 40.1% (GRW) to 54.8% (SRW). One-way ANOVA revealed that taxa abundance significantly differed among river watersheds at p < 0.01, whereas there were no great differences for the relative abundance of each taxonomic group among watersheds (Table 2b).
Table 1. Descriptive statistics of (a) regional and physical instream variables and (b) chemical variables in the five major river watersheds: the Han River (HRW), Nakdong River (NRW), Geum River (GRW), Youngsan River (YRW) and Seomjin River Watershed (SRW). Top and bottom lines of each variable indicate the average with standard deviation (in parenthesis) and the range, respectively. Kruskal–Wallis test (K–W) was performed with the environmental variables to compare the differences of each variable among river watersheds. The same small letters indicate no significant difference based on Dunn’s multiple comparison tests. Abbreviations: DO, dissolved oxygen; BOD, biochemical oxygen demand; EC, electrical conductivity; TN, total nitrogen; TP, total phosphorus.
Table 1. Descriptive statistics of (a) regional and physical instream variables and (b) chemical variables in the five major river watersheds: the Han River (HRW), Nakdong River (NRW), Geum River (GRW), Youngsan River (YRW) and Seomjin River Watershed (SRW). Top and bottom lines of each variable indicate the average with standard deviation (in parenthesis) and the range, respectively. Kruskal–Wallis test (K–W) was performed with the environmental variables to compare the differences of each variable among river watersheds. The same small letters indicate no significant difference based on Dunn’s multiple comparison tests. Abbreviations: DO, dissolved oxygen; BOD, biochemical oxygen demand; EC, electrical conductivity; TN, total nitrogen; TP, total phosphorus.
VariableHRW (n = 320)NRW (n = 130)GRW (n = 130)YRW (n = 76)SRW (n = 64)Total (n = 720)p (K–W)
(a) Regional and physical instream variables
Altitude (m)147.5 (151.6) d89.6 (119.1) b,c57.7 (61.2) a,b32.5 (35.1) a118.1 (79.3) c,d106.1 (126.2)0.000
1.0–721.01.0–629.00.0–278.00.0–211.01.0–335.00.0–721.0
% Urban31.7 (32.0) b23.0 (30.8) a,b30.9 (35.8) b31.4 (30.2) b14.6 (12.4) a28.4 (31.5)0.000
0–1000–1000–1000–900–800–100
% Agriculture24.7 (27.5) a42.5 (31.0) b,c46.6 (38.9) c45.4 (30.7) c35.9 (20.7) b35.0 (31.7)0.000
0–1000–1000–1000–900–800–100
% Forest35.8 (33.2) b31.4 (27.8) b17.9 (30.1) a22.6 (24.3) a48.8 (24.6) c31.5 (31.3)0.000
0–1000–1000–1000–1000–950–100
Water velocity (cm/s)53.8 (27.9) c14.7 (14.0) a38.1 (34.5) b17.7 (13.1) a41.9 (26.0) b39.0 (30.4)0.000
0.0–140.00.0–67.00.0–137.70.4–47.51.4–98.20.0–140.0
% Fine particles31.1 (27.9) b36.5 (36.8) b,c43.6 (34.5) c54.7 (31.3) d10.8 (22.9) a35.0 (32.6)0.000
0.0–100.00.0–100.00.0–100.00.0–100.00.0–100.00.0–100.0
% Coarse particles68.6 (28.1) c63.5 (36.8) b,c56.4 (34.5) b45.3 (31.3) a89.2 (22.9) d64.8 (32.7)0.000
0.0–100.00.0–100.00.0–100.00.0–100.00.0–100.00.0–100.0
(b) Chemical variables
pH8.1 (0.8) b,c8.0 (0.8) a,b8.3 (0.9) c8.0 (0.7) a,b7.8 (0.7) a8.1 (0.8)0.003
6.5–10.16.2–10.67.0–11.16.7–10.16.7–9.66.2–11.1
DO (mg/L)9.78 (2.47) a10.81 (2.36) b10.86 (3.36) b10.06 (1.69) a9.57 (1.32) a10.17 (2.54)0.000
2.42–16.102.55–17.342.74–17.866.19–15.407.11–12.692.42–17.86
BOD (mg/L)3.1 (3.7) b,c1.9 (1.3) a3.7 (2.0) c3.5 (3.0) b,c2.7 (1.7) b3.0 (2.9)0.000
0.3–37.50.4–10.40.8–9.10.3–13.30.3–12.30.3–37.5
EC (μS/cm)299.8 (337.6) a1358.7 (6067.1) b404.8 (349.8) a270.7 (402.6) a1284.6 (5823.4) b594.4 (3141.6)0.000
10.3–2729.019.7–44000.086.1–2780.028.5–3082.032.6–33360.010.3–44000.0
TN (mg/L)3.15 (2.97) b2.10 (1.66) a3.60 (2.42) b3.36 (3.59) b2.32 (1.39) a2.99 (2.69)0.000
0.38–23.780.32–11.300.29–14.270.69–27.710.44–5.770.29–27.71
TP (mg/L)0.15 (0.41) b,c0.07 (0.15) a,b0.14 (0.15) b,c0.21 (0.28) c0.05 (0.06) a0.13 (0.30)0.000
0.00–5.590.00–0.910.01–1.010.01–1.660.00–0.320.00–5.59
Turbidity (NTU)9.3 (16.6) b8.6 (7.4) b13.9 (22.3) b7.8 (45.9) b2.1 (3.3) a9.2 (21.2)0.000
0.0–152.00.0–34.20.4–182.40.0–400.00.0–16.90.0–400.0
Table 2. Assemblage attributes (a) and average abundance (individuals·m−2) (b) of benthic macroinvertebrates with standard deviation (in parenthesis) in five major river watersheds: the Han River (HRW), Nakdong River (NRW), Geum River (GRW), Youngsan River (YRW), and Seomjin River Watershed (SRW). Relative abundance (RA) indicated as an average of total density for each taxonomic group. Kruskal–Wallis test (K–W) was performed for each taxonomic group. The same small letters indicate no significant difference based on Dunn’s multiple comparison test.
Table 2. Assemblage attributes (a) and average abundance (individuals·m−2) (b) of benthic macroinvertebrates with standard deviation (in parenthesis) in five major river watersheds: the Han River (HRW), Nakdong River (NRW), Geum River (GRW), Youngsan River (YRW), and Seomjin River Watershed (SRW). Relative abundance (RA) indicated as an average of total density for each taxonomic group. Kruskal–Wallis test (K–W) was performed for each taxonomic group. The same small letters indicate no significant difference based on Dunn’s multiple comparison test.
Biological AttributesHRW (n = 320)NRW (n = 130)GRW (n = 130)YRW (n = 76)SRW (n = 64)Total (n = 720)RAp (K–W)
(a) Assemblage attributes
Taxa richness15.7 (9.4) a11.7 (7.3) b15.5 (10.8) a10.6 (6.3) b15.1 (8.0) a14.4 (9.1)-0.000
EPT richness10.3 (8.0) c5.3 (5.5) a7.6 (8.2) b4.5 (4.7) a8.5 (5.4) b8.1 (7.5)-0.000
Taxa abundance2867.5 (9791.9) b864.9 (1035.2) a3944.4 (5964.7) b947.0 (1054.1) a612.0 (548.2) a2297.1 (7121.0)-0.000
EPT abundance1066.9 (1702.9) b303.2 (726.5) a1552.5 (2441.0) c316.4 (545.0) a278.6 (315.8) a867.4 (1647.2)-0.000
Dominance index0.65 (0.21) b0.70 (0.20) b,c0.71 (0.21) b,c0.73 (0.20) c0.52 (0.23) a0.67 (0.22)-0.000
Shannon diversity index2.33 (1.00) b2.08 (0.91) b2.09 (1.02) b1.79 (0.94) a2.63 (1.09) c2.21 (1.01)-0.000
(b) Taxa abundance of higher taxonomic group
Non-InsectaPlatyhelminthes18.7 (78.8)28.7 (196.7)25.5 (102.4)9.9 (49.2)17.2 (37.8)20.7 (109.6)0.010.447
Nematomorpha0.2 (0.9)0.1 (0.6)0.1 (0.5)0.1 (0.8)0.0 (0.0)0.1 (0.7)0.000.082
Mollusca10.7 (39.4) a44.4 (151.5) c30.4 (54.3) b27.4 (49.5) b21.2 (30.9) a,b23.0 (76.5)0.010.000
Annelida603.5 (6826.2) b39.9 (90.8) a548.3 (2997.6) b21.0 (61.4) a68.1 (91.7) a382.7 (4731.0)0.200.000
Crustacea1.9 (12.2) a42.9 (214.3) a41.2 (243.9) a109.8 (560.6) b1.4 (7.0) a27.7 (230.6)0.010.000
InsectaEphemeroptera692.6 (1132.9) b316.6 (594.6) a565.1 (1015.1) b199.0 (308.7) a377.3 (389.0) a521.5 (936.0)0.270.000
Odonata3.9 (12.3) a4.2 (12.4) a14.5 (37.1) b11.9 (29.2) b0.7 (2.2) a6.4 (21.3)0.000.000
Plecoptera6.6 (21.6) b5.6 (29.6) b4.2 (21.5) a,b0.1 (0.7) a0.7 (4.2) a4.8 (21.3)0.000.000
Hemiptera1.1 (9.3) a38.7 (207.4) b25.0 (129.1) a,b25.5 (213.5) a,b0.0 (0.3) a14.7 (125.7)0.010.000
Megaloptera2.6 (10.4) b1.2 (6.0) a,b1.6 (6.0) a,b0.5 (2.3) a1.0 (9.8) a,b1.8 (8.4)0.000.000
Coleoptera15.3 (97.6) a28.6 (80.6) a80.9 (200.8) c8.9 (31.5) a56.9 (94.9) b32.6 (119.3)0.020.000
Diptera449.3 (780.2) b309.8 (452.1) b806.2 (1853.3) c344.1 (712.1) b105.0 (164.4) a446.9 (1009.5)0.230.000
Trichoptera549.2 (1029.0) b306.6 (771.8) a845.3 (1564.4) c114.5 (296.2) a126.4 (295.1) a475.4 (1044.8)0.240.000
Lepidoptera0.0 (0.1) a-0.1 (0.9) b0.0 (0.6) a,b-0.0 (0.4)0.000.000
Neuroptera---0.0 (0.3)-0.0 (0.1)0.000.076
Korean stream ecosystems were characterized by a few predominant colonizers and a majority of rare taxa in macroinvertebrate assemblage composition. The most abundant and widespread taxa were a worm (Limnodrilus gotoi Hatai) and midge larvae (Chironomidae spp.) with their relative abundance of approximately 50% of total density throughout the whole river watershed. Other dominant species were mayflies (i.e., Baetis fuscatus (Linnaeus), Epeorus pellucidus (Brodsky) and Uracanthella punctisetae (Matsumura)) and netspinning caddisflies (Cheumatopsyche brevilineata Iwata, Hydropsyche valvata Martynov and Hydropsyche kozhantschikovi Martynov) among the different stream and river systems, most of which were dominant in somewhat nutrient-rich habitats at middle or lower streams. However, over 50% (195 taxa) of the fauna was present with low occurring frequency (less than 2% of all sites) and 90% with low abundance (less than 0.2%).
The environmental relationships of the dominant taxa were stronger with the physical variables (i.e., altitude, water velocity and streambed conditions) than with the chemical variables, among which water velocity was the most significant parameter. Consequently, significant positive relationships existed in most dominant taxa with peak abundance at a moderate velocity of 50–100 cm/s, particularly for E. pellucidus (r = 0.410, p < 0.001) and C. brevilineata (r = 0.435, p < 0.001) (Figure 2). Baetiella tuberculata (Kazlauskas), U. punctisetae and H. valvata tended to occur in their highest densities at the fast velocity (120–140 cm/s). No clear tendency was observed for Chironomidae spp.
Figure 2. Distribution patterns of dominant macroinvertebrates along with water velocity. Data are given as arithmetic means with standard deviation. Spearman’s rank correlation coefficients between water velocity and abundance of each dominant species are included. (a) Uracanthella punctisetae; (b) Epeorus pellucidus; (c) Baetiella tuberculata; (d) Chironomini sp.; (e) Cheumatopsyche brevilineata; (f) Hydropsyche valvata.
Figure 2. Distribution patterns of dominant macroinvertebrates along with water velocity. Data are given as arithmetic means with standard deviation. Spearman’s rank correlation coefficients between water velocity and abundance of each dominant species are included. (a) Uracanthella punctisetae; (b) Epeorus pellucidus; (c) Baetiella tuberculata; (d) Chironomini sp.; (e) Cheumatopsyche brevilineata; (f) Hydropsyche valvata.
Water 08 00027 g002

3.3. Macroinvertebrate-Based Site Classification

The cluster analysis, based on the similarity in the benthic macroinvertebrate composition, largely classified the 720 sampling sites into two clusters and subsequently sub-clustered them into five groups. As a result, 126 sampling sites were included in Group 1a; 199 sites in Groups 1b; 106 sites in Group 2a; 118 sites in Group 2b; and 171 sites in Group 2c (Figure 1). MRPP validated these five groups with significant differences (A = 0.484, p < 0.001). The differences in environmental variables among the clusters are shown in Figure 3. The most important indicator taxa for each cluster are shown with their indicator values in Table 3.
Figure 3. Distribution of selected environmental variables among the five clusters of the 720 sampling sites which were identified by cluster analysis with Sørenson distance measure. Box represents the 25th and 75th percentiles and whiskers indicate the 5th and 95th percentiles with standard deviations (error bar). The mean (horizontal dotted line) and median (horizontal solid line) are shown in each box. Different small letters indicate significant difference based on a Dunn’s multiple comparison test at p < 0.05.
Figure 3. Distribution of selected environmental variables among the five clusters of the 720 sampling sites which were identified by cluster analysis with Sørenson distance measure. Box represents the 25th and 75th percentiles and whiskers indicate the 5th and 95th percentiles with standard deviations (error bar). The mean (horizontal dotted line) and median (horizontal solid line) are shown in each box. Different small letters indicate significant difference based on a Dunn’s multiple comparison test at p < 0.05.
Water 08 00027 g003
Table 3. Indicator values (%) for the most important species (p < 0.05) in each cluster group. Monte Carlo tests (999 permutations) were used to assess the significance of each species as an indicator for the respective group (G1a–G2c). In total 50 species whose contribution to total density was higher than 0.2% were analyzed. Less important species were not shown in this table.
Table 3. Indicator values (%) for the most important species (p < 0.05) in each cluster group. Monte Carlo tests (999 permutations) were used to assess the significance of each species as an indicator for the respective group (G1a–G2c). In total 50 species whose contribution to total density was higher than 0.2% were analyzed. Less important species were not shown in this table.
TaxaCluster Groupp
G1aG1bG2aG2bG2c
Rhyacophila nigrocephala Iwata47160000.001
Epeorus nipponicus (Uéno)4610000.001
Glossosoma KUa4440000.001
Drunella aculea (Allen)4410000.001
Hydropsyche orientalis Martynov4071100.001
Uracanthella punctisetae (Matsumura)21486600.001
Hydropsyche valvata Martynov6462900.001
Cheumatopsyche brevilineata Iwata154151000.001
Hydropsyche kozhantschikovi Martynov214031600.001
Psychomyia sp.0321500.001
Ephemera orientalis McLachlan22041520.001
Ecdyonurus levis (Navás)52840100.001
Ecdyonurus joernensis Bengtsson01225100.001
Mataeopsephus KUa01623000.001
Asellus sp.0211360.001
Hirudo nipponia Whitman11022840.001
Chironomini sp.17927210.001
Limnodrilus gotoi Hatai512182390.001
Micronecta sedula Horváth0000200.001
Physa acuta Draparnaud23212160.001
Micronecta sp.005050.007
Total number of significant indicator species1221833-
Each cluster was clearly differentiated according to the differences of instream physicochemical conditions and geographical location of sampling sites. The cluster analysis discriminated less polluted streams with low EC and fast flowing water (Group 1) from moderately or severely polluted steams with high EC and slow water velocity (Group 2). Group 1a (G1a), congregating in the HRW, consisted specifically of mountainous upper streams although a tenth of this group was scattered over the other watersheds except the YRW. This group consisted of oligotrophic streams with distinguishing features of the highest altitude, the lowest BOD and nutrient concentrations, and the fastest water velocity. Additionally, the catchment was predominantly comprised of forested area. The best indicator species for this group were characterized as high sensitivity against organic pollution (e.g., Rhyacophila nigrocephala Iwata, Epeorus nipponicus (Uéno), Drunella aculea (Allen) and Hydropsyche orientalis Martynov). Group 1b (G1b) contained the largest number of sampling sites and was widely distributed throughout all river watersheds but mostly encompassing the agricultural and forested catchment. G1a and G1b closely resembled chemical environments, but steam sites belonging to G1b displayed mesotrophic condition with slightly higher BOD and nutrients than those of G1a. There existed the largest number of indicator species in G1b, among which U. punctisetae showed the highest indicator value.
G2a mostly included sites located in middle reaches of large rivers and their tributaries, particularly in the NRW and SRW. Although Group 2a (G2a) was similar in overall environmental characteristics to G1b, G2a was characterized with slightly lower altitude, higher EC and a lower water velocity when compared with G1b. Sampling sites suffering from poor water quality with organic degradation and nutrient enrichment were confined to Group 2b with the highest BOD, TN, and TP concentrations. These sites were influenced by the highest degree of agriculture and/or urbanization. The significant indicator species of this group were Hirudo nipponia Whitman, Chironomini sp. and L. gotoi, indicating high trophy and saprobity [35]. Water quality conditions in Group 2c (G2c) were as poor as those of G2b with the highest EC. Sites in this group were characterized with the lowest water velocity and the largest proportion of fine particles because G2c was gathered by streams adjacent to estuaries, large rivers, and dammed streams. The only three weak indicator species appeared in G2c.

3.4. Environmental Variables Affecting Macroinvertebrate Distributions

A CCA was performed to understand how macroinvertebrates were distributed along environmental gradients (Figure 4). Total variability explained in the species data was 15.9% (Table 4). The eigenvalues of the first CCA axis (0.281) and the second CCA axis (0.101) were significant (p < 0.01; 99 Monte Carlo permutation test). All three CCA results were significant based on a Monte Carlo permutation test (p < 0.01).
Figure 4. Ordination plots constructed by canonical correspondence analysis for all sampling sites (n = 720). The plots present the ordination of sampling sites (a) and macroinvertebrates (b) with relative contributions of environmental variables (c) for the first two axes. Only species contributing >0.2% to total abundance were included. The summary of ordination results is presented in Table 4 and the codes for each taxon are shown in Appendix Table A1.
Figure 4. Ordination plots constructed by canonical correspondence analysis for all sampling sites (n = 720). The plots present the ordination of sampling sites (a) and macroinvertebrates (b) with relative contributions of environmental variables (c) for the first two axes. Only species contributing >0.2% to total abundance were included. The summary of ordination results is presented in Table 4 and the codes for each taxon are shown in Appendix Table A1.
Water 08 00027 g004
Table 4. Spearman’s rank correlation coefficients and probabilities of environmental variables and the CCA axes (** p < 0.01, * p < 0.05) (a), and summary of CCA results (n = 720) (b). All axes were significant based on Monte Carlo permutation procedures.
Table 4. Spearman’s rank correlation coefficients and probabilities of environmental variables and the CCA axes (** p < 0.01, * p < 0.05) (a), and summary of CCA results (n = 720) (b). All axes were significant based on Monte Carlo permutation procedures.
VariablesAxis 1Axis 2Axis 3
(a) Correlation coefficients
Altitude0.793 **0.236 **0.013
Stream order−0.105 **−0.658 **0.099 **
% Urban−0.218 **0.134 **−0.317 **
% Agriculture−0.199 **−0.273 **0.133 **
% Forest0.504 **0.0550.266 **
Water velocity0.662 **−0.054−0.510 **
% Fine particles−0.744 **0.325 **−0.076 *
% Coarse particles0.747 **−0.327 **0.075 *
pH0.011−0.236 **0.007
DO0.304 **−0.467 **−0.166 **
BOD−0.550 **0.164 **−0.370 **
EC−0.196 **−0.278 **−0.116 **
TN−0.463 **0.266 **−0.387 **
TP−0.373 **0.229 **−0.314 **
(b) Summary of CCA results
Eigenvalue0.2810.1010.073
% variance explained in taxa data9.83.52.5
Cumulative % variance explained9.813.315.9
p value0.0100.0100.010
Total variance2.869--
The five clusters were well separated in the CCA ordination plot (Figure 4a). Both the G1a and G1b groups with good physicochemical environments were positioned on the right side of the ordination plot, and G2b and G2c were on the opposite side. However, the G2a group straddled both sides. Eight environmental variables (i.e., BOD, TN, TP, altitude, water velocity, % forest, % fine particles, and % coarse particles) had significant correlations with the first axis, among which altitude was the most significant contributor (r = 0.793, p < 0.01), followed by % coarse particles (r = 0.747, p < 0.01) and % fine particles (r = −0.744, p < 0.01) (Table 4). All physical factors except for % fine particles were in the opposite direction from the chemical variables (Figure 4c). This result suggests that the decline in altitude reflected deterioration in water quality, accompanied by increases in organic material and nutrients. The heterogeneity of the macroinvertebrate habitats also decreased with the longitudinal gradient toward downstream. The second axis was negatively related with stream order (r = −0.658, p < 0.01) and DO (r = −0.467, p < 0.01). High positive scores with the first axis were denoted for the rhithronic (i.e., pertaining to the headwaters) and intolerant species (e.g., Drunella aculea (Allen), E. nipponicus, Glossosoma sp., R. nigrocephala Iwata and H. orientalis), whereas negative scores were observed for the potamic (i.e., pertaining to rivers) and tolerant taxa (e.g., Micronecta sedula Horváth, Physa acuta Draparnaud, Labiobaetis atrebatinus (Eaton) and Asellus hilgendorfii Bovalius) (Figure 4b).

4. Discussion

The Asian monsoon region is a global biodiversity hotspot suffering from increasing anthropogenic disturbances, but aquatic biodiversity and ecosystem integrity remain poorly explored [3]. This is the same situation in Korea, but recent establishment of the National Aquatic Ecological Monitoring Program [36] opened a new era to assess ecosystem health and biodiversity in Korea. Our work presented in this paper takes advantage of the opportunity of such a nationwide scale of survey. Although the five river watersheds in Korea exhibited differences in environmental conditions and macroinvertebrate taxa abundance, overall taxonomic composition at each watershed displayed little difference (Table 2). Instead, a considerable spatial variation in benthic macroinvertebrate assemblages accounted for the combination of both ultimate (e.g., altitude and the degree of land use) and proximate factors (e.g., flow, stream bed substrate, BOD, EC, and nutrients). This finding was in line with previous studies of benthic diatoms in Korea [25] and macroinvertebrates in other Asian countries [12,37], suggesting the importance of various multi-scale factors in structuring macroinvertebrate assemblages [38]. Our results provide basal information for the sustainable management and conservation practices of stream ecosystems.

4.1. Macroinvertebrate Taxonomic Composition

The macroinvertebrate assemblages in Korean stream ecosystems generally bear resemblance to those in tropical and other temperate streams at higher taxonomic levels [35,39]. However, temperate Korean streams were relatively rich in Ephemeroptera, Diptera, and Trichoptera when compared with tropical streams for the higher biodiversity of Gastropoda, Decapoda and other insect orders such as Odonata and Hemiptera [39,40] (Table 2). The taxonomic composition also included a large number of rhithronic fauna that prefer stony substrates. For example, Baetidae, Heptageniidae, and Ephemerellidae mainly dominated the Ephemeroptera, and Rhyacophilidae and Hydropsychidae composed the Trichoptera in this study.
We confirmed a total of 340 macroinvertebrate taxa in this study. Of the macroinvertebrate taxa, chironomid midge larvae and a small minnow mayfly (B. fuscatus) were extensively encountered throughout Korean stream environments with the highest occurring frequency with 94% and 76% of total sampling sites, respectively. Temperate Korean streams were also characterized with a few predominant colonizers and a majority of rare taxa in macroinvertebrate taxonomic composition. Moreover, only six cosmopolitan taxa were noticeably comprised of >60% of all macroinvertebrate samples, which were L. gotoi, Chironomidae spp., U. punctisetae, C. brevilineata, H. valvata and H. kozhantschikovi. These dominant taxa were found to be significant indicators of mesotrophic or polytrophic streams in our study (Table 3), consistent with a result that most Korean streams and rivers were degraded in both chemical and biological status [36]. On the other hand, most Korean streams were occupied by a great number of rare macroinvertebrate species with a small distribution range and/or low abundance, as was also demonstrated in other studies [41,42].

4.2. Environmental Relationships with Macroinvertebrate Distribution

Benthic habitats are complex, and a variety of environmental variables acting at multiple spatial scales regulate the composition and distribution patterns of stream macroinvertebrate assemblages in an exclusive or synergistic fashion [38,43]. We revealed that the variables associated with altitude and in-stream habitats best accounted for the largest amount of variability in our macroinvertebrate data set, supporting the more important determinants of local environmental factors than broad or regional parameters [16,44]. The importance and role of local environmental variables have also been highlighted in other aquatic communities: aquatic macrophytes [45], freshwater phytoplankton [46], benthic diatoms [44], intertidal macroinvertebrates [47], and fish [48]. However, a comprehensive understanding of multispatial scales is needed because of significant correlations between macrohabitat and microhabitat characteristics, depending on the relative size of the area studied [3].
In our study, the variability among the macroinvertebrate-based stream groups was more prominently explained by the altitudinal gradients together with streambed composition and water velocity than chemical variables (Table 3). First, the most widely accepted theory related to altitudinal changes is the river continuum concept (RCC), which displays structural and functional responses to the longitudinal gradient [49]. Altitude has also been well documented in other studies as a main descriptor to determine macroinvertebrate richness as well as other environmental variables such as temperature, hydrology, food availability, streambed condition, and water chemistry [3,50]. However, we revealed poor relationships of altitude with water chemistry in contrast with significant associations with physical variables and biological attributes such as taxa richness, EPT richness, EPT abundance, and the Shannon diversity index. This may be due to the fact that over half of the studied streams corresponded to lowland streams below 100 m a.s.l., which were characterized by moderate or slightly poor water quality with severe variations in BOD, TN, TP, and Shannon diversity-based saprobity [36,51] (Table 1). Harding et al. [52] demonstrated that increasing human land use intensity along a river continuum caused water quality degradation, consequently leading to changes in taxonomic composition from intolerant EPT dominated to tolerant taxa-dominated assemblages despite a gradual increase in taxa abundance by a few dominant species.
Riffle habitats are commonly characterized by shallow water depth, oxygenating fast-flowing water, and stony beds, and exhibit higher taxa richness and abundance than that at pools or habitats with fine sediment, as in our case [3,42,53]. Such hydraulic conditions are critical determinants for the distribution and species composition of benthic organisms [16,54] and are the driving forces for evolution of their morphologies and life history [11]. Rheophilous (i.e., having an affinity for running waters) and limnophilous taxa were clearly discriminated based on the correlation analysis (Figure 2) and CCA (Figure 4) results. For example, well-known rhithronic species in two Ephemeropteran (Ephemerellidae and Heptageniidae) and one Trichopteran (Hydropsychidae) families mostly displayed high preference to the mid- or fast current conditions, which corresponds to the intolerant scraper or filtering collector groups abundant in upper and middle stream reaches [49]. Thus, biological monitoring and assessment programs using rheophilic macroinvertebrates would benefit from their ecological characteristics and their indication of good environmental quality (e.g., [42]) although seasonal fluctuations in hydraulic parameters by the Asian monsoon more critically cause a catastrophic drift and washout of benthic organisms.
Streambed composition is one of the most important factors to directly influence richness and abundance of macroinvertebrates on local scales [12,16]; thus, close correlations would be expected with variables related to longitudinal changes in stream ecosystems. Boulders and cobbles are typically the major structural elements in steep gradient upper streams, whereas sand and smaller sediments predominate in the lower reaches [55]. There have been a large number of studies on the macroinvertebrate-substrate relationship, most of which have revealed that macroinvertebrate diversity and density increase with higher heterogeneity due to the available stable and diverse microhabitats (e.g., [10,43]). Therefore, factors determining macroinvertebrate communities at local scales obviously put a priority on the streambed conditions rather than water quality or other physical variables. We found that the mainstreams of the Nakdong River, which contained about 60% fine particles, retained the lowest taxa richness (8.2 on average, n = 18) and abundance (747 individuals/m2) among five major rivers despite its good water quality condition (mean BOD, 1.9 and TN, 1.8 mg/L) with the exception of the severely polluted Youngsan River (8.9 taxa richness and 7.3 BOD, n = 13). These results indicate that homogeneous streambeds with greater fine particles support lower diversity and abundance even when streams maintain good water quality [56]. Additionally, organic pollution could transform the coarse substrate into an organic rich soft bottom, altering the community structure from dominated by diverse and intolerant species to communities predominated by a few tolerant species [3,20], as in the case of the Youngsan River.

4.3. Considerations to Improve Macroinvertebrate Biomonitoring Programs

Not until 2006 did Korea adopt biological water quality criteria and the concept of ecological integrity in the water quality program [25,36]. Since then, Korean government has led a biological survey of streams and rivers (i.e., NAEMP) every year to assess the current biological status of stream and river ecosystems, and to develop a strategy for the restoration and management of disturbed systems [25]. As a part of a nationwide survey benthic macroinvertebrates are also monitored based on the guidelines and assessment tools [27]. Notwithstanding their suitability for convenient and rapid bio-assessment, there remain debatable issues as to whether the methods for sampling and treating macroinvertebrates are effective to provide a reliable and accurate indication of the macroinvertebrate fauna throughout the country.
The Korean national biomonitoring program presents field sampling in a cost-effective and time-saving manner for rapid bio-assessment. To satisfy this purpose benthic macroinvertebrates are optimized to be collected at riffle and/or gliding run habitats using a Surber sampler with three replicates [27]. However, such sampling methods probably underestimate overall biodiversity in stream ecosystems, considering the whole stream environment. First, single-habitat sampling possibly produces incomplete taxa lists and includes no target habitats at certain sites despite its advantage that the influences of both water and habitat quality on macroinvertebrates are not confounded by instream habitat variation [57]. In this regard, a multihabitat approach would be more profitable for the estimation of taxonomic diversity due to its consistent application across stream types especially at large-scaled survey, comprehensive taxa lists, and effective assessment of ecological conditions [58]. Second, the Surber sampler is one of the most commonly used quantitative tools in lotic systems and provides high-precision information on the abundance and composition of macroinvertebrate assemblages [21,59]. This method is, on the other hand, usually more appropriate for riffle habitats of shallow streams, presumably underestimating overall biodiversity in a region. Recent studies on the comparison between sampling devices suggested that artificial substrates (e.g., leaf-bags) would be used as complementary tools due to their discriminative taxonomic composition [60,61,62]. Third, a majority of biomonitoring programs widely adopt three to five replicates for lotic ecosystems because one way to reduce monitoring costs is to decrease the sample size [21,63,64]. The previous studies also showed that such a small number of replicates rarely influenced ecological health assessment, particularly for the indices applying sensitivity/tolerance taxa, as the same for NAEMP [65,66,67]. Nevertheless, small sample sizes may influence the values of biological measures and the representativeness for a real benthic community based on the asymptotic relationship of the number of taxa with both the sampling area and the collected individuals [68]. Finally, sieve mesh size could also affect the accurate estimates of taxonomic diversity and environmental quality assessment [67]. Finer mesh sizes more accurately represent macroinvertebrate assemblages than coarser mesh sizes, whereas they need more efforts to handle specimens. On the contrary, coarser mesh sizes undervalue the original density and taxa richness of macroinvertebrate assemblages by passing small individuals through the sieves and result in higher diversity [69]. Considering the tradeoff between cost and effectiveness, many rapid bio-assessments determine 0.5 mm mesh size as a reasonable choice [21,64,70]. Therefore, future researches are required for identifying optimal sampling effort comprehensively considering cost-effectiveness, easy applicability, and well representative of resident macroinvertebrate assemblages.

5. Conclusions

Large-scale knowledge on environmental relationships of aquatic communities is crucial to conserve freshwater biodiversity and sustain ecological integrity. Korean stream macroinvertebrate assemblages were determined not only by regional and physical instream variables but also by pollution-related parameters at nation-wide scale. Macroinvertebrate-based site classification evidently provided an environmental characterization for different types of Korean streams, indicating that benthic macroinvertebrates are a valuable biomonitoring material. The results of this study provide important information and a bridge for further work such as the assemblage—specific responses to environmental disturbance. Our results also contribute to establishing effective management practices, implementing conservation measures, and developing more reliable biological monitoring tools for sustainable freshwater ecosystems.

Acknowledgments

This study was performed under the project of “National Aquatic Ecosystem Health Survey and Assessment” in Korea, and was supported by the Ministry of Environment and the National Institute of Environmental Research, Korea. The authors are grateful to survey members involved in the project. The authors also thank the anonymous reviewers for their help in improving the scientific content of the manuscript.

Author Contributions

Yung-Chul Jun developed the concept of the study, performed data analysis, and wrote the manuscript. Nan-Young Kim assisted data processing and figure construction. Sang-Hun Kim conducted quality assurance of the physico-chemical and biological data. Young-Seuk Park and Dong-Soo Kong assisted study design and statistical analysis. Soon-Jin Hwang supervised the research, and assisted data interpretation and manuscript preparation. All the authors contributed to the review of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Codes for macroinvertebrates contributing to more than 0.2% of total abundance for canonical correspondence analysis in Figure 4.
Table A1. Codes for macroinvertebrates contributing to more than 0.2% of total abundance for canonical correspondence analysis in Figure 4.
CodeTaxon
Phylum Platyhelminthes
DugSpDugesia sp.
Phylum Mollusca
SemLiSemisulcospira libertina (Gould)
PhyAcPhysa acuta Draparnaud
Phylum Annelida
LimGoLimnodrilus gotoi Hatai
HirNiHirudo nipponia Whitman
Phylum Arthropoda
Class Crustacea
AseSpAsellus sp.
GamSpGammarus sp.
Class Insecta
Order Ephemeroptera
BaeTuBaetiella tuberculata (Kazlauskas)
BaeFuBaetis fuscatus (Linnaeus)
BaeUrBaetis ursinus Kazlauskas
LabAtLabiobaetis atrebatinus (Eaton)
NigBaNigrobaetis bacillus (Kluge)
EcdBaEcdyonurus bajkovae Kluge
EcdJoEcdyonurus joernensis Bengtsson
EcdKiEcdyonurus kibunensis Imanishi
EcdLeEcdyonurus levis (Navás)
EpeNiEpeorus nipponicus (Uéno)
EpeLaEpeorus latifolium (Uéno)
EpePeEpeorus pellucidus (Brodsky)
ChoAlChoroterpes altioculus Kluge
ParJaParaleptophlebia japonica (Matsumura)
PotFoPotamanthus formosus Eaton
RhoCoRhoenanthus coreanus (Yoon and Bae)
DruAcDrunella aculea (Allen)
EphOrEphemera orientalis McLachlan
SerSeSerratella setigera (Bajkova)
UraPuUracanthella punctisetae (Matsumura)
CaeNiCaenis nishinoae Malzacher
Order Hemiptera
MicSeMicronecta sedula Horváth
MicSpMicronecta sp.
Order Coleoptera
ElmSpElmidae sp.
EubKaEubrianax KUa
MatKaMataeopsephus KUa
PseKaPsephenoides KUa
Order Diptera
AntKaAntocha KUa
CulSpCulex sp.
ChiSpChironomidae spp. (non-red type)
ChiReChironomini spp. (red-type)
Order Trichoptera
RhyNiRhyacophila nigrocephala Iwata
HydKaHydroptila KUa
GloKaGlossosoma KUa
CheBrCheumatopsyche brevilineata Iwata
CheKaCheumatopsyche KUa
CheKbCheumatopsyche KUb
HydKoHydropsyche kozhantschikovi Martynov
HydKbHydropsyche KUb
HydOrHydropsyche orientalis Martynov
HydVaHydropsyche valvata Martynov
MacRaMacrostemum radiatum McLachlan
PsySpPsychomyia sp.

References

  1. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.I.; Knowler, D.J.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.H.; Soto, D.; Stiassny, M.L.J.; et al. Freshwater biodiversity: Importance, threats, status and conservation challenges. Biol. Rev. 2006, 81, 163–182. [Google Scholar] [CrossRef] [PubMed]
  2. Strayer, D.L.; Dudgeon, D. Freshwater biodiversity conservation: Recent progress and future challenges. J. N. Am. Benthol. Soc. 2010, 29, 344–358. [Google Scholar] [CrossRef]
  3. Li, F.; Chung, N.; Bae, M.-J.; Kwon, Y.-S.; Park, Y.-S. Relationships between stream macroinvertebrates and environmental variables at multiple spatial scales. Freshwater Biol. 2012, 57, 2107–2124. [Google Scholar] [CrossRef]
  4. Dudgeon, D. The ecology of tropical Asian rivers and streams in relation to biodiversity conservation. Annu. Rev. Ecol. Syst. 2000, 31, 239–263. [Google Scholar] [CrossRef]
  5. De Silva, S.S.; Abery, N.W.; Nguyen, T.T.T. Endemic freshwater finfish of Asia: Distribution and conservation status. Divers. Distrib. 2007, 13, 172–184. [Google Scholar] [CrossRef]
  6. Vörösmarty, C.J.; McIntyre, P.B.; Gessner, M.O.; Dudgeon, D.; Prusevich, A.; Green, P.; Glidden, S.; Bunn, S.E.; Sullivan, C.A.; Liermann, C.E.; et al. Global threats to human water security and river biodiversity. Nature 2010, 467, 555–561. [Google Scholar] [CrossRef] [PubMed]
  7. Strayer, D.L. Challenges for freshwater invertebrate conservation. J. N. Am. Benthol. Soc. 2006, 25, 271–287. [Google Scholar] [CrossRef]
  8. Buss, D.F.; Baptista, D.F.; Silveira, M.P.; Nessimian, J.L.; Dorvillé, L.F.M. Influence of water chemistry and environmental degradation on macroinvertebrate assemblages in a river basin in south-east Brazil. Hydrobiologia 2002, 481, 125–136. [Google Scholar] [CrossRef]
  9. Bott, T.L.; Brock, J.T.; Dunn, C.S.; Naiman, R.J.; Ovink, R.W.; Petersen, R.C. Benthic community metabolism in four temperate stream systems: An inter-biome comparison and evaluation of the river continuum concept. Hydrobiologia 1985, 220, 109–117. [Google Scholar] [CrossRef]
  10. Merz, J.R.; Ochikubo Chan, L.K. Effects of gravel augmentation on macroinvertebrate assemblages in a regulated California river. River Res. Appl. 2005, 21, 61–74. [Google Scholar] [CrossRef]
  11. Nelson, S.M.; Lieberman, D.M. The influence of flow and other environmental factors on benthic invertebrates in the Sacramento River, USA. Hydrobiologia 2002, 489, 117–129. [Google Scholar] [CrossRef]
  12. Jiang, X.M.; Xiong, J.; Qiu, J.W.; Wu, J.M.; Wang, J.W.; Xie, Z.C. Structure of macroinvertebrate communities in relation to environmental variables in a subtropical Asian river system. Int. Rev. Hydrobiol. 2010, 95, 42–57. [Google Scholar] [CrossRef]
  13. Vought, L.B.M.; Kullberg, A.; Petersen, R.C. Effect of riparian structure, temperature and channel morphometry on detritus processing in channelized and natural woodland streams in southern Sweden. Aquat. Conserv. Mar. Freshw. Ecosyst. 1998, 8, 273–285. [Google Scholar] [CrossRef]
  14. Townsend, C.R.; Hildrew, A.G.; Francis, J. Community structure in some southern English streams: The influence of physicochemical factors. Freshw. Biol. 1983, 13, 521–544. [Google Scholar] [CrossRef]
  15. Richards, C.; Haro, R.J.; Johnson, L.B.; Host, G.E. Catchment and reach-scale properties as indicator of macroinvertebrate species traits. Freshw. Biol. 1997, 37, 219–230. [Google Scholar] [CrossRef]
  16. Sandin, L. Benthic macroinvertebrates in Swedish streams: Community structure, taxon richness, and environmental relations. Ecography 2003, 26, 269–282. [Google Scholar] [CrossRef]
  17. Leps, M.; Tonkin, J.D.; Dahm, V.; Haase, P.; Sundermann, A. Disentangling environmental drivers of benthic invertebrate assemblages: The role of spatial scale and riverscape heterogeneity in a multiple stressor environment. Sci. Total Environ. 2015, 536, 546–556. [Google Scholar] [CrossRef] [PubMed]
  18. Heino, J. Biodiversity of aquatic insects: Spatial gradients and environmental correlates of assemblage-level measures at large scales. Freshw. Rev. 2009, 2, 1–29. [Google Scholar] [CrossRef]
  19. Shah, D.N.; Tonkin, J.D.; Haase, P.; Jähnig, S.C. Latitudinal patterns and large-scale environmental determinants of stream insect richness across Europe. Limnologica 2015, 55, 33–43. [Google Scholar] [CrossRef]
  20. Rossaro, B.; Pietrangelo, A. Macroinvertebrate distribution in streams: A comparison of CA ordination with biotic indices. Hydrobiologia 1993, 263, 109–118. [Google Scholar] [CrossRef]
  21. Rosenberg, D.M.; Resh, V.H. Freshwater Biomonitoring and Benthic Macroinvertebrates; Chapman and Hall: New York, NY, USA, 1993. [Google Scholar]
  22. Charvet, S.; Statzner, B.; Usseglio-Polatera, P.; Dumont, B. Traits of benthic macroinvertebrates in semi-natural French streams: An initial application to biomonitoring in Europe. Freshw. Biol. 2000, 43, 277–296. [Google Scholar] [CrossRef]
  23. Korea Meteorological Administration (KMA). The Meteorological Yearbook in Korea; KMA: Seoul, Korea, 2010. (in Korean) [Google Scholar]
  24. Jeong, K.S.; Hong, D.G.; Byeon, M.S.; Jeong, J.C.; Kim, H.G.; Kim, D.K.; Joo, G.J. Stream modification patterns in a river basin: Field survey and self-organizing map (SOM) application. Ecol. Inform. 2010, 5, 293–303. [Google Scholar] [CrossRef]
  25. Hwang, S.-J.; Kim, J.-Y.; Yoon, S.-A.; Kim, B.-H.; Park, M.-H.; You, K.-A.; Lee, H.-Y.; Kim, H.-S.; Kim, Y.-J.; Lee, J.; et al. Distribution of benthic diatoms in Korean rivers and streams in relation to environmental variables. Ann. Limnol. Int. J. Limnol. 2011, 47, S15–S33. [Google Scholar] [CrossRef]
  26. Craig, D.A. Some of what you should know about water. Bull. N. Am. Benthol. Soc. 1987, 4, 178–182. [Google Scholar]
  27. The Ministry of Environment/National Institute of Environmental Research (MOE/NIER). Survey and Evaluation of Aquatic Ecosystem Health in Korea; MOE/NIER: Incheon, Korea, 2012. [Google Scholar]
  28. American Public Health Association (APHA). Standard Methods for the Examination of Water and Wastewater, 21st ed.; APHA: Washington, DC, USA, 2001. [Google Scholar]
  29. McCune, B.; Grace, J.B. Analysis of Ecological Communities; MjM Software Design: Gleneden Beach, OR, USA, 2002. [Google Scholar]
  30. Mielke, P.W.; Berry, K.J.; Johnson, E.S. Multiresponse permutation procedures for a priori classifications. Commun. Stat. 1976, A5, 1409–1424. [Google Scholar] [CrossRef]
  31. Zimmerman, G.M.; Goetz, H.; Mielke, P.W. Use of an improved statistical method for group comparisons to study effects of prairie fire. Ecology 1985, 66, 606–611. [Google Scholar] [CrossRef]
  32. Dufrene, M.; Legendre, P. Species assemblages and indicator species: The need for flexible asymmetrical approach. Ecol. Monogr. 1997, 67, 345–366. [Google Scholar] [CrossRef]
  33. Petersen, W.T.; Keister, J.E. Interannual variability in copepod community composition at a coastal station in the northern California Current: A multivariate approach. Deep Sea Res. 2003, 50, 2499–2517. [Google Scholar] [CrossRef]
  34. McCune, B.; Mefford, M.J. Multivariate Analysis of Ecological Data (Version 4.25); MjM Software: Gleneden Beach, OR, USA, 1999. [Google Scholar]
  35. Wang, X.; Cai, Q.; Tang, T.; Yang, S.; Li, F. Spatial distribution of benthic macroinvertebrates in the Erhai basin of southwestern China. J. Freshw. Ecol. 2012, 27, 89–96. [Google Scholar] [CrossRef]
  36. Lee, S.-W.; Hwang, S.-J.; Lee, J.-K.; Jung, D.-I.; Park, Y.-J.; Kim, J.-T. Overview and application of the National Aquatic Ecological Monitoring Program (NAEMP) in Korea. Ann. Limnol. Int. J. Limnol. 2011, 47, S3–S14. [Google Scholar] [CrossRef]
  37. Hoang, T.H.; Lock, K.; Chi Dang, K.; de Pauw, N.; Goethals, P.L.M. Spatial and temporal patterns of macroinvertebrate communities in the Du River Basin in Northern Vietnam. J. Freshw. Ecol. 2010, 25, 637–647. [Google Scholar] [CrossRef]
  38. Bae, M.-J.; Kwon, Y.; Hwang, S.-J.; Con, T.-S.; Yang, H.-J.; Kwak, I.-S.; Park, J.-H.; Ham, S.-A.; Park, Y.-S. Relationships between three major stream assemblages and their environmental factors in multiple spatial scales. Ann. Limnol. Int. J. Limnol. 2011, 47, S91–S105. [Google Scholar] [CrossRef]
  39. Jacobsen, D.; Cressa, C.; Mathooko, J.M.; Dudgeon, D. Macroinvertebrates: Composition, life histories and production. In Aquatic Ecosystems: Tropical Stream Ecology; Dudgeon, D., Ed.; Elsevier Science: London, UK, 2008; pp. 65–105. [Google Scholar]
  40. Hoang, D.H.; Bae, Y.J. Aquatic insect diversity in a tropical Vietnamese stream in comparison with that in a temperate Korean stream. Limnology 2006, 7, 45–55. [Google Scholar] [CrossRef]
  41. Cao, Y.; Williams, D.D.; Williams, N.E. How important rare species in aquatic community ecology and bioassessment? Limnol. Oceanogr. 1998, 43, 1403–1409. [Google Scholar] [CrossRef]
  42. Nijboer, R.C.; Schmidt-Kloiber, A. The effect of excluding taxa with low abundances or taxa with small distribution ranges on ecological assessment. Hydrobiologia 2004, 516, 347–363. [Google Scholar] [CrossRef]
  43. Braccia, A.; Voshell, J.R. Environmental factors accounting for benthic macroinvertebrate assemblage structure at the sample scale in streams subjected to a gradient of cattle grazing. Hydrobiologia 2006, 573, 55–73. [Google Scholar] [CrossRef]
  44. Soininen, J.; Paavola, R.; Muotka, T. Benthic diatom communities in boreal streams: Community structure in relation to environmental and spatial gradients. Ecography 2004, 27, 330–342. [Google Scholar] [CrossRef]
  45. Mikulyuk, A.; Sharma, S.; Egeren, S.V.; Erdmann, E.; Nault, M.E.; Hauxwell, J. The relative role of environmental, spatial, and land-use patterns in explaining aquatic macrophyte community composition. Can. J. Fish. Aquat. Sci. 2011, 68, 1778–1789. [Google Scholar] [CrossRef]
  46. Stomp, M.; Huisman, J.; Mittelbach, G.G.; Litchman, E.; Klausmeier, C.A. Large-scale biodiversity patterns in freshwater phytoplankton. Ecology 2011, 92, 2096–2107. [Google Scholar] [CrossRef] [PubMed]
  47. Rodrigues, A.M.; Quintino, V.; Sampaio, L.; Freitas, R.; Neves, R. Benthic biodiversity patterns in Ria de Aveiro, Western Portugal: Environmental-biological relationships. Estuar. Coast. Shelf Sci. 2011, 95, 338–348. [Google Scholar] [CrossRef]
  48. Waite, I.R.; Carpenter, K.D. Associations among fish assemblage structure and environmental variables in Willamette Basin streams, Oregon. Trans. Am. Fish. Soc. 2000, 129, 754–770. [Google Scholar] [CrossRef]
  49. Vannote, R.L.; Minshall, G.W.; Cummins, K.W.; Sedell, J.R.; Cushing, C.E. The river continuum concept. Can. J. Fish. Aquat. Sci. 1980, 37, 130–137. [Google Scholar] [CrossRef]
  50. Beauchard, O.; Gagneur, J.; Brosse, S. Macroinvertebrate richness patterns in North African streams. J. Biogeogr. 2003, 30, 1821–1833. [Google Scholar] [CrossRef]
  51. Pan, Y.; Hill, B.H.; Husby, P.; Hall, R.K.; Kaufmann, P.R. Relationships between environmental variables and benthic diatom assemblages in California Central Valley streams (USA). Hydrobiologia 2006, 561, 119–130. [Google Scholar] [CrossRef]
  52. Harding, J.S.; Young, R.G.; Hayes, J.W.; Shearer, K.A.; Stark, J.D. Changes in agricultural intensity and river health along ariver continuum. Freshw. Biol. 1999, 42, 345–357. [Google Scholar] [CrossRef]
  53. Boyero, L.; Bosch, J. The effect of riffle-scale environmental variability on macroinvertebrate assemblages in a tropical stream. Hydrobiologia 2004, 524, 125–132. [Google Scholar] [CrossRef]
  54. Potapova, M.G.; Charles, D.F. Benthic diatoms in USA rivers: Distributions along spatial and environmental gradients. J. Biogeogr. 2002, 29, 167–187. [Google Scholar] [CrossRef]
  55. Halwas, K.; Church, M. Channel units in small, high gradient streams on Vancouver Island, British Columbia. Geomorphology 2002, 43, 243–256. [Google Scholar] [CrossRef]
  56. Jun, Y.-C.; Kim, N.-Y.; Kwon, S.-J.; Han, S.-C.; Hwang, I.-C.; Park, J.-H.; Won, D.-H.; Byun, M.-S.; Kong, H.-Y.; Lee, J.-E.; et al. Effects of land use on benthic macroinvertebrate communities: Comparison of two mountain streams in Korea. Ann. Limnol. Int. J. Limnol. 2011, 47, S35–S49. [Google Scholar] [CrossRef]
  57. Parsons, M.; Norris, R.H. The effect of habitat-specific sampling on biological assessment of water quality using a predictive model. Freshw. Biol. 1996, 36, 416–434. [Google Scholar] [CrossRef]
  58. Gerth, W.J.; Herlihy, A.T. Effect of sampling different habitat types in regional macroinvertebrate bioassessment surveys. J. N. Am. Benthol. Soc. 2006, 25, 501–512. [Google Scholar] [CrossRef]
  59. Gilles, C.L.; Hose, G.C.; Turak, E. What do qualitative rapid assessment collections of macroinvertebrates represent? A comparison with extensive quantitative sampling. Environ. Monit. Assess. 2009, 149, 99–112. [Google Scholar] [CrossRef] [PubMed]
  60. Pashkevich, A.; Pavluk, T.; de Vaate, A.B. Efficiency of a standardized artificial substrate for biological monitoring of river water quality. Environ. Monit. Assess. 1996, 40, 143–156. [Google Scholar] [CrossRef] [PubMed]
  61. Saliu, J.K.; Ovuorie, U.R. The artificial substrate preference of invertebrates in Ogbe Creek, Lagos, Nigeria. Life Sci. J. 2007, 4, 77–81. [Google Scholar]
  62. Di Sabatino, A.; Cristiano, G.; Pinna, M.; Lombardo, P.; Miccoli, F.P.; Marini, G.; Vignini, P.; Cicolani, B. Structure, functional organization and biological traits of macroinvertebrate assemblages from leaf-bags and benthic samples in a third-order stream of Central Apennines (Italy). Ecol. Indic. 2014, 46, 84–91. [Google Scholar] [CrossRef]
  63. Vlek, H.E.; Šporka, F.; Krno, I. Influence of macroinvertebrate sample size on bioassessment of streams. Hydrobiologia 2006, 566, 523–542. [Google Scholar] [CrossRef]
  64. Carter, J.L.; Resh, V.H. After site selection and before data analysis: Sampling, sorting, and laboratory procedures used in stream benthic macroinvertebrate monitoring programs by USA state agencies. J. N. Am. Benthol. Soc. 2001, 20, 658–682. [Google Scholar] [CrossRef]
  65. Metzeling, L.; Miller, J. Evaluation of the sample size used for the rapid bioassessment of rivers using macroinvertebrates. Hydrobiologia 2001, 444, 159–170. [Google Scholar] [CrossRef]
  66. Kim, A.R.; Oh, M.W.; Kong, D.S. The influence of sample size on environmental assessment using benthic macroinvertebrates. J. Korean. Soc. Water Environ. 2013, 29, 790–798. [Google Scholar]
  67. Pinna, M.; Marini, G.; Mancinelli, G.; Basset, A. Influence of sampling effort on ecological descriptors and indicators in perturbed and unperturbed conditions: A study case using benthic macroinvertebrates in Mediterranean transitional waters. Ecol. Indic. 2014, 37, 27–39. [Google Scholar] [CrossRef]
  68. Vinson, M.R.; Hawkins, C.P. Effects of sampling area and subsampling procedure on comparisons of taxa richness among streams. J. N. Am. Benthol. Soc. 1996, 15, 392–399. [Google Scholar] [CrossRef]
  69. Marini, G.; Pinna, M.; Basset, A.; Mancinelli, G. Estimation of benthic macroinvertebrate taxonomic diversity: Testing the role of sampling effort in a Mediterranean transitional water ecosystem. Transit. Waters Bull. 2013, 7, 28–40. [Google Scholar]
  70. Buss, D.F.; Borges, E.L. Application of rapid bioassessment protocols (RBP) for benthic macroinvertebrates in Brazil: Comparison between sampling techniques and mesh sizes. Neotrop. Entomol. 2008, 37, 288–295. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Jun, Y.-C.; Kim, N.-Y.; Kim, S.-H.; Park, Y.-S.; Kong, D.-S.; Hwang, S.-J. Spatial Distribution of Benthic Macroinvertebrate Assemblages in Relation to Environmental Variables in Korean Nationwide Streams. Water 2016, 8, 27. https://doi.org/10.3390/w8010027

AMA Style

Jun Y-C, Kim N-Y, Kim S-H, Park Y-S, Kong D-S, Hwang S-J. Spatial Distribution of Benthic Macroinvertebrate Assemblages in Relation to Environmental Variables in Korean Nationwide Streams. Water. 2016; 8(1):27. https://doi.org/10.3390/w8010027

Chicago/Turabian Style

Jun, Yung-Chul, Nan-Young Kim, Sang-Hun Kim, Young-Seuk Park, Dong-Soo Kong, and Soon-Jin Hwang. 2016. "Spatial Distribution of Benthic Macroinvertebrate Assemblages in Relation to Environmental Variables in Korean Nationwide Streams" Water 8, no. 1: 27. https://doi.org/10.3390/w8010027

APA Style

Jun, Y. -C., Kim, N. -Y., Kim, S. -H., Park, Y. -S., Kong, D. -S., & Hwang, S. -J. (2016). Spatial Distribution of Benthic Macroinvertebrate Assemblages in Relation to Environmental Variables in Korean Nationwide Streams. Water, 8(1), 27. https://doi.org/10.3390/w8010027

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