Using a Geospatial Model to Relate Fluvial Geomorphology to Macroinvertebrate Habitat in a Prairie River—part 1: Genus-level Relationships with Geomorphic Typologies

Modern river ecosystems undergo constant stress from disturbances such as bank stabilization, channelization, dams, and municipal, agricultural, and industrial water use. As these anthropogenic water requirements persist, more efficient methods of characterizing river reaches are essential. Benthic macroinvertebrates are helpful when evaluating fluvial health, because they are often the first group to react to contaminants that can then be transferred through them to other trophic levels. Hence, the purpose of this research is to use a geospatial model to differentiate instream macroinvertebrate habitats, and determine if the model is a viable method for stream evaluation. Through the use of ArcGIS and digital elevation models, the fluvial geomorphology of the Qu'Appelle River in Saskatchewan (SK) was assessed. Four geomorphological characteristics of the river were isolated (sinuosity, slope, fractal dimension, and stream width) and clustered through Principle Component Analysis (PCA), yielding sets of river reaches with similar geomorphological characteristics, called typologies. These typologies were mapped to form a geospatial model of the river. Macroinvertebrate data were aligned to the locations of the typologies, revealing several relationships with the fluvial geomorphology. A Kruskal-Wallis analysis and post hoc pairwise multiple comparisons were completed with the macroinvertebrate data to pinpoint significant genera, as related to the geospatial model.


Introduction
The study of aquatic macroinvertebrates within fluvial systems can reveal information about water quality [1,2], interactions with other trophic levels [3,4], as well as anthropogenic influences on river health [5,6].Traditional methods used to examine the fluvial benthic community rely on fieldwork and data collection; as such, the methods can be time consuming, labour intensive, and expensive.More efficient methods have emerged for studying macroinvertebrates, but these are often tailored to the interests of specific studies.This could be a result of the sheer vastness of fluvial ecosystems and macroinvertebrates.The dynamic nature of river ecosystems and the highly diverse macroinvertebrate community commonly found in fluvial ecosystems make it difficult to study them as a whole.Therefore, studies often choose to organize macroinvertebrate data collection and analysis by metahabitat (i.e., riffle, pool, run) or instream habitats that are classified using physical traits.Methods used in studies related to analyzing physical habitat include the IndVal method for communities [7][8][9], analyses using the species-based habitat suitability criteria (HSC) that are subsequently incorporated into habitat simulation models [10][11][12][13], the fairly new random forests nonparametric method [9,14], and many more.These methods relate collected macroinvertebrate data back to the river metahabitats for correlation.Other studies focus on certain groups of macroinvertebrates, such as the more sensitive Ephemeroptera, Plecoptera, and Trichoptera (EPT) orders [2,4], or macroinvertebrates found in niche stream habitats (e.g., rock substrates).Methods for stone substrate macroinvertebrate studies began using labour-intensive field collection, but modern methods, such as surface area estimation, have since been developed [15,16].Methods used to understand macroinvertebrates within large fluvial systems that do not rely on field data have not been widely developed.
Creating more effective and efficient methods of studying macroinvertebrates is a worthy goal, because a better understanding of macroinvertebrates leads to increased knowledge of the watershed as a whole.Macroinvertebrates are the food source for many trophic levels, are often the first to respond to ecosystem degradation, and can be the point of contaminant uptake [17,18].In the case of contaminants, macroinvertebrates can pass harmful substances to those that feed on them, like fish and birds, creating a pattern of bioaccumulation.Contaminants can be linked to geomorphology as well; fluvial geomorphology directly influences instream flow and subsequently the distribution of sediment, creating the many instream habitats macroinvertebrates live in [19,20].Consequently, the sediment type and distribution (driven by geomorphology) directly relates to many dangerous substances that macroinvertebrates may ingest.Clay, for instance, has such a small grain size that it is negatively charged, making bonds with the positively charged metals much easier [19].The metals derived from industrial processes, such as mining, increase aquatic metal concentrations (i.e., mercury, cadmium, copper, zinc, lead) as they settle and bond to the clay in the river [19,21].Macroinvertebrates living in this clay can ingest these metals, leading to bioaccumulation in later trophic levels.In the case of the more toxic metals, like mercury, bioaccumulation can reach terrestrial consumers, including humans [21].Understanding fluvial geomorphology and instream macroinvertebrate habitats is vital to further understanding dangerous bioaccumulation and contaminant relations.
The River Continuum Concept (RCC) exemplifies the notion of viewing rivers as predictable systems assuming that macroinvertebrate habitat is organized by linear location [22].The RCC has been disputed for not applying to real and complex fluvial environments.Instead of concepts like the RCC, spatial dynamics and geomorphological processes at a local-level are widely accepted as the driving factors behind instream habitats and communities [23,24].Furthermore, fluvial geomorphology has long been associated with instream biodiversity, meaning that river macroinvertebrates are influenced by river structure [4,25,26].Aquatic macroinvertebrate habitat selection is influenced by local features, which result from extensive geomorphological and hydrological processes [20,27].Through classifying reaches of the river with similar geomorphological characteristics (typologies), rather than by metacommunities, patterns in macroinvertebrate habitat can be pursued.This information could then be applied to rivers from which data are difficult to collect, due to their large spatial extent [28].Therefore, creating a model of geomorphological characteristics within a river is a novel method for identifying macroinvertebrate habitat more broadly.Unlike the linear RCC, the use of geomorphic typologies acknowledges the fact that diverse instream habitats, resulting from channel morphology interacting with dynamic hydrological processes, are the reality of modern rivers [29].Four geomorphological characteristics (sinuosity, stream width, fractal dimension, and slope) have been focused on, in the creation of the geospatial model for this study.These characteristics are all affected by anthropogenic influences in modern rivers.Understandably, channelization changes the geomorphology of a river substantially.For the purposes of flood control, land drainage, or erosion mitigation, the channel may be straightened, realigned, or enlarged [30].Other modifications include dredging, building embankments, bank stabilization, and, of course, dams [30].These drastic anthropogenic alterations are examples of how the geomorphological characteristics studied here are certain to change with increased societal use of rivers.In this study, valuable macroinvertebrate We hypothesize that geomorphic typologies within the river will reveal significant relationships between geomorphology and the benthic macroinvertebrate community in a fluvial system.The study concentrates on macroinvertebrates at the genus level; the data set was tested to determine which genera contain significant relationships to geomorphological typologies.These genera were then the focus of further study.

Study Site
The Qu'Appelle River flows in the Saskatchewan plains region of Canada, which is home to sweeping river valleys, wide lake plains, ground moraines, and large spillways prone to flooding [31].The Qu'Appelle is a river and lake network; it begins at the Qu'Appelle Dam at Lake Diefenbaker and flows toward the Manitoba border (Figure 1).The basin's dynamic geomorphology provides diverse habitats for many species, including more than 30 rare and endangered animal and plant species.Among these species are bigmouth buffalo fish, loggerhead shrike, and smooth arid goosefoot, all relying on the varied landscape and water within the valley [32].The Qu'Appelle River Watershed provides water to several municipalities, including two cities (Moose-Jaw and Regina), whose combined population could reach half a million by the year 2050 [33].Besides municipal water use, the system accommodates many local interests including recreational, agricultural, industrial, and environmental requirements [34].Water quality and supply are fundamental for the Qu'Appelle system and its residents, as it boasts a diverse fish assemblage and active recreational fishery.The Qu'Appelle River Watershed area also hosts a large amount of Saskatchewan potash production, for which water is essential; demand for water will continue to rise in this region as potash mining and other industrial practices expand [33,35].Channelization has occurred from the outlet at the Qu'Appelle River Dam to about 32 km downstream; the remaining river stretches' freely-developed and meandering expanses are prone to flooding, erosion, sloughing, silting, as well as water quality and quantity issues [31,33].Heavy reliance on water by people and industry within the watershed necessitates sustainable and effective monitoring and protection efforts for the Qu'Appelle River.
increased societal use of rivers.In this study, valuable macroinvertebrate data were analyzed in relation to a geomorphological model of the Qu'Appelle River, Saskatchewan.We hypothesize that geomorphic typologies within the river will reveal significant relationships between geomorphology and the benthic macroinvertebrate community in a fluvial system.The study concentrates on macroinvertebrates at the genus level; the data set was tested to determine which genera contain significant relationships to geomorphological typologies.These genera were then the focus of further study.

Study Site
The Qu'Appelle River flows in the Saskatchewan plains region of Canada, which is home to sweeping river valleys, wide lake plains, ground moraines, and large spillways prone to flooding [31].The Qu'Appelle is a river and lake network; it begins at the Qu'Appelle Dam at Lake Diefenbaker and flows toward the Manitoba border (Figure 1).The basin's dynamic geomorphology provides diverse habitats for many species, including more than 30 rare and endangered animal and plant species.Among these species are bigmouth buffalo fish, loggerhead shrike, and smooth arid goosefoot, all relying on the varied landscape and water within the valley [32].The Qu'Appelle River Watershed provides water to several municipalities, including two cities (Moose-Jaw and Regina), whose combined population could reach half a million by the year 2050 [33].Besides municipal water use, the system accommodates many local interests including recreational, agricultural, industrial, and environmental requirements [34].Water quality and supply are fundamental for the Qu'Appelle system and its residents, as it boasts a diverse fish assemblage and active recreational fishery.The Qu'Appelle River Watershed area also hosts a large amount of Saskatchewan potash production, for which water is essential; demand for water will continue to rise in this region as potash mining and other industrial practices expand [33,35].Channelization has occurred from the outlet at the Qu'Appelle River Dam to about 32 km downstream; the remaining river stretches' freely-developed and meandering expanses are prone to flooding, erosion, sloughing, silting, as well as water quality and quantity issues [31,33].Heavy reliance on water by people and industry within the watershed necessitates sustainable and effective monitoring and protection efforts for the Qu'Appelle River.

Geospatial Factors
In order to characterize the geospatial nature of the river, four characteristics were selected as points of comparison; these were stream width, sinuosity, fractal dimension, and slope.These characteristics were chosen because they provide a good representation of the watershed geomorphology.Other studies wishing to use the geospatial model delineated in this study could include a variety of parameters, such as average flow, depth, or water temperature to create a different river profile.First, the centerline of the river was delineated through the use of river polygons, and points placed along the line every 50 m.At each point, transects were placed, crossing the polygon from bank to bank (Figure 2).The 1:50,000 digital elevation model (DEM) data (Source: Department of Natural Resources Canada) provided elevation at each bank point using Geographic Information System (GIS) software, and its average elevation linked to the point at the centerline.Finally, at each centerline point, the four geospatial characteristics were extracted [28].Slope, stream width, sinuosity, and fractal dimension were calculated and extracted at each centerline point delineated.Slope is simply the level of vertical drop in river bed elevation at each 50 m point [30].Stream width is the length of the transect at each centerline point (Figure 2).The minimum stream width recorded at the centerline points was 6.8 m, the maximum was 2226.9 m, and the average was 168.7 m.The commercial software package Mathcad ® v.15 (PTC, Inc., Needham, MA, USA) [36] was used in calculating both the sinuosity and fractal dimension of the river reaches.For the sinuosity, a marginally larger scale was used than for the stream width and slope; 40 adjacent points were used (20 upstream and 20 downstream of each centerline point).Sinuosity, or how much a river meanders, is a significant fluvial characteristic because it impacts how water flows, how sediment is distributed, and ultimately how local instream habitats form [37]. Sinuosity is also defined as the interchange between the valley slope and the river slope [38], or between the length of the valley (shortest possible distance) and the length of the river.In addition to different definitions of sinuosity, various terms are used to describe the spectrum of sinuosity values (S).These include S = 1 for a linear (straight) reach, and in succession of increasing sinuosity, elongated, oscillating, tortuous, and lastly meandering with S > 2 [39,40].Rivers with high sinuosity also have increased bed material of a smaller grain size (i.e., silt, clay) [41].Where channels meander, augmented bank erosion transpires, thus providing a higher load of easily transferred (finematerials) sediment to the stream [19], with S > 1.5 indicating a riverbed containing over 92% fine textured material [38].Within the Qu'Appelle River, 63% has a sinuosity of S >1.5, suggesting that the riverbed has a share of both coarse and fine materials, corresponding to the till blanket that comprises much of the surficial geology of the river basin.Sinuosity, in our study, is defined as the length of a river reach divided by the shortest distance between the reach end points (a straight path).A sinuosity score of one represents a straight channel, and >1 is a meandering one [37,42].In order to characterize the geospatial nature of the river, four characteristics were selected as points of comparison; these were stream width, sinuosity, fractal dimension, and slope.These characteristics were chosen because they provide a good representation of the watershed geomorphology.Other studies wishing to use the geospatial model delineated in this study could include a variety of parameters, such as average flow, depth, or water temperature to create a different river profile.First, the centerline of the river was delineated through the use of river polygons, and points placed along the line every 50 m.At each point, transects were placed, crossing the polygon from bank to bank (Figure 2).The 1:50,000 digital elevation model (DEM) data (Source: Department of Natural Resources Canada) provided elevation at each bank point using Geographic Information System (GIS) software, and its average elevation linked to the point at the centerline.Finally, at each centerline point, the four geospatial characteristics were extracted [28].Slope, stream width, sinuosity, and fractal dimension were calculated and extracted at each centerline point delineated.Slope is simply the level of vertical drop in river bed elevation at each 50 m point [30].Stream width is the length of the transect at each centerline point (Figure 2).The minimum stream width recorded at the centerline points was 6.8 m, the maximum was 2226.9 m, and the average was 168.7 m.The commercial software package Mathcad ® v.15 (PTC, Inc., Needham, MA, USA) [36] was used in calculating both the sinuosity and fractal dimension of the river reaches.For the sinuosity, a marginally larger scale was used than for the stream width and slope; 40 adjacent points were used (20 upstream and 20 downstream of each centerline point).Sinuosity, or how much a river meanders, is a significant fluvial characteristic because it impacts how water flows, how sediment is distributed, and ultimately how local instream habitats form [37]. Sinuosity is also defined as the interchange between the valley slope and the river slope [38], or between the length of the valley (shortest possible distance) and the length of the river.In addition to different definitions of sinuosity, various terms are used to describe the spectrum of sinuosity values (S).These include S = 1 for a linear (straight) reach, and in succession of increasing sinuosity, elongated, oscillating, tortuous, and lastly meandering with S > 2 [39,40].Rivers with high sinuosity also have increased bed material of a smaller grain size (i.e., silt, clay) [41].Where channels meander, augmented bank erosion transpires, thus providing a higher load of easily transferred (fine materials) sediment to the stream [19], with S > 1.5 indicating a riverbed containing over 92% fine textured material [38].Within the Qu'Appelle River, 63% has a sinuosity of S >1.5, suggesting that the riverbed has a share of both coarse and fine materials, corresponding to the till blanket that comprises much of the surficial geology of the river basin.Sinuosity, in our study, is defined as the length of a river reach divided by the shortest distance between the reach end points (a straight path).A sinuosity score of one represents a straight channel, and >1 is a meandering one [37,42].Analogous to sinuosity, but at a greater scale, fractal dimension represents how much Euclidean space a river reach fills.Different than sinuosity, fractal dimension quantifies form and indicates both the number and amplitude of meander changes within a channel length, at a larger order of magnitude [42][43][44].As a tool, fractal dimension is helpful for understanding patterns within fluvial systems at a more expansive, geological scale [28,42].Fractal dimension is uniquely suited to the complex and irregular nature of a river.It uses a minimum value of one and a Analogous to sinuosity, but at a greater scale, fractal dimension represents how much Euclidean space a river reach fills.Different than sinuosity, fractal dimension quantifies form and indicates both the number and amplitude of meander changes within a channel length, at a larger order of Water 2016, 8, 42 5 of 17 magnitude [42][43][44].As a tool, fractal dimension is helpful for understanding patterns within fluvial systems at a more expansive, geological scale [28,42].Fractal dimension is uniquely suited to the complex and irregular nature of a river.It uses a minimum value of one and a maximum of two; in linear river systems, the higher the value is, the more space-filling it is within the drainage area [43].To calculate fractal dimension, the number of centerline points contained within a 40 km ˆ40 km window moving along the river was used.These scales were selected, as they were the values where the peak variation was calculated.Once these geospatial characteristics were found, a Principal Component Analysis (PCA) was performed with the statistical package R 2.15.2 [45] to pinpoint groups (typologies) of similar characteristics [28].The PCA makes a large dataset more manageable by identifying patterns within it and reducing its dimensionality [46].To improve normality of the dataset while considering the high incidence of zero and minute slope values, the dataset was log 10 + 1 transformed prior to PCA.The explained variance of principal components (PC) one through four was 49%, 24%, 14% and 13%, respectively.Only the first two principal components (accounting for 73% cumulative variance) were used to derive typologies as, following the Kaiser-Guttman rule, eigenvalues for components 3 and 4 were far less than 1 (0.57 and 0.52, respectively).At each centerline point, binary values were assigned to the scores associated with PCs one and two, 1 for positive PC scores and 0 for negative PC scores.Binary values were then summed, resulting in four unique typologies.Each of the four typologies was designated a colour, which was then matched to the corresponding centerline points (Figure 3).
Water 2016, 8, 42 5 of 17 maximum of two; in linear river systems, the higher the value is, the more space-filling it is within the drainage area [43].To calculate fractal dimension, the number of centerline points contained within a 40 km × 40 km window moving along the river was used.These scales were selected, as they were the values where the peak variation was calculated.Once these geospatial characteristics were found, a Principal Component Analysis (PCA) was performed with the statistical package R 2.15.2 [45] to pinpoint groups (typologies) of similar characteristics [28].The PCA makes a large dataset more manageable by identifying patterns within it and reducing its dimensionality [46].To improve normality of the dataset while considering the high incidence of zero and minute slope values, the dataset was log10 + 1 transformed prior to PCA.The explained variance of principal components (PC) one through four was 49%, 24%, 14% and 13%, respectively.Only the first two principal components (accounting for 73% cumulative variance) were used to derive typologies as, following the Kaiser-Guttman rule, eigenvalues for components 3 and 4 were far less than 1 (0.57 and 0.52, respectively).At each centerline point, binary values were assigned to the scores associated with PCs one and two, 1 for positive PC scores and 0 for negative PC scores.Binary values were then summed, resulting in four unique typologies.Each of the four typologies was designated a colour, which was then matched to the corresponding centerline points (Figure 3).

Macroinvertebrate Sampling
Benthic macroinvertebrate collection took place between 2006 and 2009 along the Qu'Appelle River, Saskatchewan.The macroinvertebrate collection and compilation of the database used in this study were completed before the delineation of the geospatial model, and separate from this study.Optimally, the typologies would have been delineated prior to sampling, so that sites could be selected with a more even distribution amongst typologies.The sampling method used was the travelling kick and sweep (TK & S) technique, a standardized sampling method that permits the comparison of communities within similar habitats between sites.Sampling sites were selected to accommodate for road access, at a pre-determined minimum river distance of 10 km between sites.During this process, a bank-to-bank transect is sampled at five different positions and combined into one composite sample.Since the TK & S method involves sampling along transects, the subsequent data set spans an array of habitats found within a reach.A D-frame net (500 µm mesh),

Macroinvertebrate Sampling
Benthic macroinvertebrate collection took place between 2006 and 2009 along the Qu'Appelle River, Saskatchewan.The macroinvertebrate collection and compilation of the database used in this study were completed before the delineation of the geospatial model, and separate from this study.Optimally, the typologies would have been delineated prior to sampling, so that sites could be selected with a more even distribution amongst typologies.The sampling method used was the travelling kick and sweep (TK & S) technique, a standardized sampling method that permits the comparison of communities within similar habitats between sites.Sampling sites were selected to accommodate for road access, at a pre-determined minimum river distance of 10 km between sites.During this process, a bank-to-bank transect is sampled at five different positions and combined into one composite sample.Since the TK & S method involves sampling along transects, the subsequent data set spans an array of habitats found within a reach.A D-frame net (500 µm mesh), positioned downstream from the person collecting the sample, was used in the present study.As the net is held, the area (~30 cm ˆ30 cm) is kicked to ~5 cm depth for a 10 s time period at each of the five positions along all transects.All the sampled sites were similar prairie reaches that were straight, slow flowing, and with no riffles.Substrate was generally soft silt, with less frequent and small patches of cobble and vegetation on bank positions.Ultimately, our objective was to employ a standardized method of collecting benthic macroinvertebrates in comparable habitats across all sites, to ensure the dominant communities for major habitat features would be represented.Although the TK & S method does not guarantee that all taxa present at a sampling site are collected, we feel the most dominant assemblages and habitat features have been adequately characterized.The collected macroinvertebrate samples were saved in 80% ethanol and subsequently moved to the laboratory, in which 7ˆmagnification was used to separate the macroinvertebrate samples from any organic material.Identification of the macroinvertebrates using keys for North America [17] and Western Canada, when available [47][48][49], was then completed.The Water Security Agency of Saskatchewan Invertebrate Voucher Collection (Saskatoon, SK, Canada), and the Royal Saskatchewan Museum (Regina, SK, Canada) hold the voucher series from this sample.Taxa occurrence records from the study were submitted to the Saskatchewan Conservation Data Centre with the Ministry of Environment.

Statistical Analysis
Using the R 3.1.2statistical software [45], a Kruskal-Wallis analysis of variance on ranks was performed on the macroinvertebrate genera, to determine whether their distribution yielded statistical differences between the 4 typologies.In the event of a significant outcome (p value ď 0.05), post hoc pairwise multiple comparisons were completed.A one-tailed Dunn test with Bonferroni p value correction for multiple paired tests was used, and graphical analyses were accomplished through use of the statistical program R 3.1.2[45].The data set was narrowed down for the analysis, to include only genera with ě30 individuals sampled and collected in >9% (3+) of sampling sites.This test can be used to understand whether the mean ranks are the same or different in all the groups, when dealing with an unequal sampling distribution among typologies.

Results and Discussion
The PCA scores (colour coded by typology) and variable vectors are plotted in Figure 4, with principal component one on the x-axis, and component two on the y-axis.This biplot is used to surmise the overall relationships between the four geomorphological variables (sinuosity, slope, fractal dimension, channel width), as well as their relationships within the distinct typologies (Figure 4).
Principal components one and two account for 73% of the total variation found in the dataset.Normalized values of sinuosity, slope, fractal dimension, and channel width were also charted in density plots.The plots were used to observe the qualitative contribution of these four geomorphological variables to each individual geomorphic typology based on absolute means (Figure 5).
The subsequent relationships established for all four typologies to sinuosity, slope, fractal dimension, and channel width are summarized in Table 1.Typology 0 is positively related to both slope and fractal dimension and negatively related to width.Typology 1 is positively related to slope, yet negatively related to all the other variables.Typology 2 is positively related to both sinuosity and fractal dimension, but negatively related to width.Typology 3 is negatively related to all variables, except for width.

Typology Sinuosity Slope Fractal Dimension Width
River stretches designated as Typology 0 are therefore characterized by more narrow regions of high slope and fractal dimension.Typology 1 is associated with narrow areas of high slope, low fractal dimension, and sinuosity.The features related with Typology 2 are highly sinuous regions with high fractal dimension values, yet are still narrow.Finally, Typology 3 is associated with wide reaches that are less sinuous, with low slope and fractal dimension values.

Typology Sinuosity Slope Fractal Dimension Width
River stretches designated as Typology 0 are therefore characterized by more narrow regions of high slope and fractal dimension.Typology 1 is associated with narrow areas of high slope, low fractal dimension, and sinuosity.The features related with Typology 2 are highly sinuous regions with high fractal dimension values, yet are still narrow.Finally, Typology 3 is associated with wide reaches that are less sinuous, with low slope and fractal dimension values.

Typology
Sinuosity Slope Fractal Dimension Width River stretches designated as Typology 0 are therefore characterized by more narrow regions of high slope and fractal dimension.Typology 1 is associated with narrow areas of high slope, low fractal dimension, and sinuosity.The features related with Typology 2 are highly sinuous regions with high fractal dimension values, yet are still narrow.Finally, Typology 3 is associated with wide reaches that are less sinuous, with low slope and fractal dimension values.

Macroinvertebrate Survey
Macroinvertebrate sampling occurred at 35 sites along the Qu'Appelle River, resulting in a database of 115,696 individuals, which were identified to 128 groups at the genus level.It should be noted that some of these 128 taxonomic groups found at the genus level were identified only to class (e.g., Oligochaeta), order (e.g., Amphipoda), or family (e.g., Heptageniidae), in cases where further taxonomic breakdown could not be determined.Sampling occurred prior to delineation of the geospatial model, so an imbalanced 14 of the macroinvertebrate sampling sites correspond with Typology 0 (blue).In total, these sights contain 14,764 individuals, an average of 1055 per site, with the highest taxa richness of 85. Typology 1 (red) also contains a large portion of sampling sites (11/35), and contained 75 taxa at the genus level.This typology also has the highest macroinvertebrate count (74% of the whole population) at 85,270 individuals, averaging 7752 per site.With only five macroinvertebrate sampling sites, Typology 2 (yellow) has the lowest population sum (2124), averaging 425 per site.It makes up merely 2% of the entire macroinvertebrate survey, and has a taxa richness of 49.Typology 3 (green) is represented by five sites as well, and contains a count of 66 different taxa.This typology has a total of 13,538 individual macroinvertebrates that average 2708 individuals per site.

Statistical Analysis of Macroinvertebrate Survey
The 128 genera identified in the macroinvertebrate survey were narrowed down to 50 genera for use in the statistical analyses (Kruskal-Wallis test; Table 2), as the remaining 78 genera were found to be less common.The parameters to designate rare genera were that they occurred at <3 sampling sites and/or they had a sum of ď30 total individuals collected in the macroinvertebrate survey.The Dunn's post hoc test provides adjusted p values for each typology comparison of the six genera with significant Kruskal-Wallis test results.If the adjusted p value is ď0.05, we assume that the ranking of values is significantly different between typologies; therefore, we can infer that the genus likely favours certain typologies over others.The outcomes of the Kruskal-Wallis test (Table 2) reveal significant differences in genus distribution among typologies for six of the 50 tested genera.These were Callicorixa, Dubiraphia, Haliplus, Heptagenia, Heptageniidae (identified to family), and Sigara.Subsequently, post hoc pairwise multiple comparisons were carried out using Dunn's Method on these six genera.A Bonferroni p value correction for multiple paired analyses pinpoints the typologies that are different from each other (Table 3).After the post hoc test, it was found that four species have significant differences between typologies (Figure 6).Typology 3 and Typology 1 are significantly different for two genera, and Typology 2 and Typology 3 are significantly different for two genera.Typology 0 has the highest total counts of all four Dubiraphia, Heptagenia, Heptageniidae, and Sigara genera.Typology 3 has the second highest total for both Dubiraphia and Sigara.Typology 2 was the second highest for Heptageniidae and Heptagenia (Figure 6).After the post hoc test, it was found that four species have significant differences between typologies (Figure 6).Typology 3 and Typology 1 are significantly different for two genera, and Typology 2 and Typology 3 are significantly different for two genera.Typology 0 has the highest total counts of all four Dubiraphia, Heptagenia, Heptageniidae, and Sigara genera.Typology 3 has the second highest total for both Dubiraphia and Sigara.Typology 2 was the second highest for Heptageniidae and Heptagenia (Figure 6).3, as related to the typologies.

Ephemeroptera, Plecoptera, Trichoptera (EPT)
The mayfly (Ephemeroptera), stonefly (Plecoptera) and caddisfly (Trichoptera) orders (EPT) are commonly used together to indicate the overall water health at a site.The groups are often sensitive to water condition and have specific habitat preferences [3,50,51].The presence of a healthy EPT population indicates a healthy and well-oxygenated ecosystem, as opposed to one dominated by macroinvertebrates more tolerant to pollution (e.g., Chironomidae midges), which may indicate poor water quality [17,52].The total EPT individuals, percent of the total population they represent (% EPT), and EPT taxa richness were calculated for each typology (Figure 7).Typology 0 (blue) contains the largest total EPT population, followed closely by Typology 1 (red), while Typology 2 (yellow) contains the smallest sum.However, Typology 2 has the highest % EPT at 52%, followed by Typology 0 (35%), and Typology 3 (18%), with Typology 1 only containing 6% EPT (Figure 7).Typology 0 has the highest EPT richness (39) of its total 85 genera.Of the 75 genera found in Typology 1, only 25 are EPT.Typology 2 has an EPT taxa richness of 18, and 49 genera overall.Typology 3 has 22 EPT taxa, and a richness of 66 genera total.3, as related to the typologies.

Ephemeroptera, Plecoptera, Trichoptera (EPT)
The mayfly (Ephemeroptera), stonefly (Plecoptera) and caddisfly (Trichoptera) orders (EPT) are commonly used together to indicate the overall water health at a site.The groups are often sensitive to water condition and have specific habitat preferences [3,50,51].The presence of a healthy EPT population indicates a healthy and well-oxygenated ecosystem, as opposed to one dominated by macroinvertebrates more tolerant to pollution (e.g., Chironomidae midges), which may indicate poor water quality [17,52].The total EPT individuals, percent of the total population they represent (% EPT), and EPT taxa richness were calculated for each typology (Figure 7).Typology 0 (blue) contains the largest total EPT population, followed closely by Typology 1 (red), while Typology 2 (yellow) contains the smallest sum.However, Typology 2 has the highest % EPT at 52%, followed by Typology 0 (35%), and Typology 3 (18%), with Typology 1 only containing 6% EPT (Figure 7).Typology 0 has the highest EPT richness (39) of its total 85 genera.Of the 75 genera found in Typology 1, only 25 are EPT.Typology 2 has an EPT taxa richness of 18, and 49 genera overall.Typology 3 has 22 EPT taxa, and a richness of 66 genera total.
Water 2016, 8, 42 10 of 17 After the post hoc test, it was found that four species have significant differences between typologies (Figure 6).Typology 3 and Typology 1 are significantly different for two genera, and Typology 2 and Typology 3 are significantly different for two genera.Typology 0 has the highest total counts of all four Dubiraphia, Heptagenia, Heptageniidae, and Sigara genera.Typology 3 has the second highest total for both Dubiraphia and Sigara.Typology 2 was the second highest for Heptageniidae and Heptagenia (Figure 6).3, as related to the typologies.

Ephemeroptera, Plecoptera, Trichoptera (EPT)
The mayfly (Ephemeroptera), stonefly (Plecoptera) and caddisfly (Trichoptera) orders (EPT) are commonly used together to indicate the overall water health at a site.The groups are often sensitive to water condition and have specific habitat preferences [3,50,51].The presence of a healthy EPT population indicates a healthy and well-oxygenated ecosystem, as opposed to one dominated by macroinvertebrates more tolerant to pollution (e.g., Chironomidae midges), which may indicate poor water quality [17,52].The total EPT individuals, percent of the total population they represent (% EPT), and EPT taxa richness were calculated for each typology (Figure 7).Typology 0 (blue) contains the largest total EPT population, followed closely by Typology 1 (red), while Typology 2 (yellow) contains the smallest sum.However, Typology 2 has the highest % EPT at 52%, followed by Typology 0 (35%), and Typology 3 (18%), with Typology 1 only containing 6% EPT (Figure 7).Typology 0 has the highest EPT richness (39) of its total 85 genera.Of the 75 genera found in Typology 1, only 25 are EPT.Typology 2 has an EPT taxa richness of 18, and 49 genera overall.Typology 3 has 22 EPT taxa, and a richness of 66 genera total.

Geomorphological Typologies
Six of the common genera tested were found to have significant differences in distribution between typologies, signifying the predilection of several macroinvertebrates for specific geomorphic typologies in this river.These six significantly different genera are the water boatman Callicorixa, the riffle beetle Dubiraphia, the water beetle Haliplus, the water boatman Sigara, the flathead mayfly only identified to the family Heptageniidae, and those Heptageniidae further identified to the genus Heptagenia (Table 2).Four of these genera (excluding the Callicorixa and Haliplus) ranked as significant, according to Dunn's post hoc test (Table 3).

Typology 0 (Blue)
The blue typology occurs throughout the system but dominates in a few areas, most notably a large expanse south of Esterhazy, SK, stretching to just before the Manitoba border (Typology 0; Figure 3).Typology 0 is characterized as being relatively narrow, with high bank slope and fractal dimension, and no relationship with sinuosity (Table 1).Narrower reaches such as these tend to have more vegetation falling into the water, creating an ideal habitat for herbivores and detrivores [19].Although this typology consists primarily of the functional feeding group collector gatherers (46%), when compared to the others Typology 0 contains the most detrivores, housing 49% of them (Figure 8).The predominance of collector gatherers in this typology could be attributed to the high abundance of Chironomidae (midges) and Amphipod Hyalella, accounting for 19% and 15% of the population respectively (Figure 9).It should be noted that Typology 0 contains the most EPT, highest EPT richness, and the second highest % EPT of the sample (Figure 7).

Geomorphological Typologies
Six of the common genera tested were found to have significant differences in distribution between typologies, signifying the predilection of several macroinvertebrates for specific geomorphic typologies in this river.These six significantly different genera are the water boatman Callicorixa, the riffle beetle Dubiraphia, the water beetle Haliplus, the water boatman Sigara, the flathead mayfly only identified to the family Heptageniidae, and those Heptageniidae further identified to the genus Heptagenia (Table 2).Four of these genera (excluding the Callicorixa and Haliplus) ranked as significant, according to Dunn's post hoc test (Table 3).

Typology 0 (Blue)
The blue typology occurs throughout the system but dominates in a few areas, most notably a large expanse south of Esterhazy, SK, stretching to just before the Manitoba border (Typology 0; Figure 3).Typology 0 is characterized as being relatively narrow, with high bank slope and fractal dimension, and no relationship with sinuosity (Table 1).Narrower reaches such as these tend to have more vegetation falling into the water, creating an ideal habitat for herbivores and detrivores [19].Although this typology consists primarily of the functional feeding group collector gatherers (46%), when compared to the others Typology 0 contains the most detrivores, housing 49% of them (Figure 8).The predominance of collector gatherers in this typology could be attributed to the high abundance of Chironomidae (midges) and Amphipod Hyalella, accounting for 19% and 15% of the population respectively (Figure 9).It should be noted that Typology 0 contains the most EPT, highest EPT richness, and the second highest % EPT of the sample (Figure 7).The population distribution is highest in Typology 0 for all four of the genera with significant links to typologies according to the statistical analysis; thus, all four demonstrate a preference for Typology 0 (Figure 6).Of the four genera, one is a water boatman, one is a riffle beetle, and two are Heptageniidae mayflies.The mayflies (Heptageniidae; Heptagenia) are nearly exclusively found in Typology 0, followed by Typology 2 at a distant second (Figure 6).Sometimes known as flat-headed mayflies because of adaptive general flattening and streamlining of the body, they cling to substrate in lotic streams [17,19].All four of the genera showing preference to Typology 0 succeed in lotic environments.The mayfly Heptageniidae belongs to the scraper functional feeding group, which makes up 19% of the population in Typology 0. They thrive in habitats where sunlight can reach, scraping biofilm and periphyton as their food source [19,50].The highest concentration of the two The population distribution is highest in Typology 0 for all four of the genera with significant links to typologies according to the statistical analysis; thus, all four demonstrate a preference for Typology 0 (Figure 6).Of the four genera, one is a water boatman, one is a riffle beetle, and two are Heptageniidae mayflies.The mayflies (Heptageniidae; Heptagenia) are nearly exclusively found in Typology 0, followed by Typology 2 at a distant second (Figure 6).Sometimes known as flat-headed mayflies because of adaptive general flattening and streamlining of the body, they cling to substrate in lotic streams [17,19].All four of the genera showing preference to Typology 0 succeed in lotic environments.The mayfly Heptageniidae belongs to the scraper functional feeding group, which makes up 19% of the population in Typology 0. They thrive in habitats where sunlight can reach, scraping biofilm and periphyton as their food source [19,50].The highest concentration of the two other significant genera (Dubiraphia; Sigara) is also found in Typology 0, followed by Typology 3. The water boatman Sigara are herbivores and swimmers; they prefer areas with ample vegetation and use their modified mouthparts to pierce plant tissue [21].The narrow reaches of Typology 0 are apt to provide access to vegetation for these herbivores.The riffle beetles Dubiraphia are collector gatherers, who cling on to substrate in lotic environments, similar to the Heptageniidae mayflies, collecting organic matter for food.
Water 2016, 8, 42 12 of 17 other significant genera (Dubiraphia; Sigara) is also found in Typology 0, followed by Typology 3. The water boatman Sigara are herbivores and swimmers; they prefer areas with ample vegetation and use their modified mouthparts to pierce plant tissue [21].The narrow reaches of Typology 0 are apt to provide access to vegetation for these herbivores.The riffle beetles Dubiraphia are collector gatherers, who cling on to substrate in lotic environments, similar to the Heptageniidae mayflies, collecting organic matter for food.

Typology 1 (Red)
Typology 1 dominates first at the culmination of the river until Buffalo Pound Lake (Figure 3; 1.a), then again from near Kennel, SK to the fishing lakes around Fort Qu'Appelle (Figure 3; 1.b), and again for a short stretch at the outlet of Katepwa Lake (Figure 3; 1.c).Geomorphological characteristics associated with Typology 1 are high bank slopes that are narrow with low sinuosity and low fractal dimension (Table 1).The channelized areas at the onset of the river are dominated by the red typology.It has the lowest percent of the sensitive EPT individuals, at only 6% of the total population, yet houses the highest total population at 74% of the sample (Figure 7).Correspondingly, this typology has a high taxa richness of 75, and merely 25 of them are EPT genera.Further examination of the expansive Typology 1 population at the genus level reveals that it is dominated by the amphipod Hyalella.This amphipod genus represents 43% of the total population held in Typology 1.The Hyalella here mark the highest portion represented by one genus in any of the typologies (Figure 9).Following Hyalella, by population size, in the red typology is the black fly larva Simulium (22%), and those only identified to the Amphipoda order (16%).Simulium can be found in almost any lotic water body, as they have many mechanisms to adapt to selective pressures [53,54].Amphipods are very tolerant as well; Hyalella in particular has been found in extremely different climates, from Guatemala to the Northwest Territories [55].Out of all four typologies, Typology 1 appears to have the most imbalanced distribution and least complex community, dominated by the tolerant amphipods and black fly larvae (Figure 9).Typology 1 also consists mostly of the collector gatherer functional feeding group (69%), followed by filterers (24%), and scrapers (4%), with very few representatives of the other groups relative to its large population (Figure 8).As compared with the other typologies, Typology 1 contains the majority of the collector gatherers (81%) and filterers (84%), as well as the most predators (44%), omnivores (37%), and scrapers (35%) of any typology.This is not surprising due to the large population held within this typology, as compared to the other three.The second highest EPT population is found here; nevertheless, it has the lowest % EPT

Typology 1 (Red)
Typology 1 dominates first at the culmination of the river until Buffalo Pound Lake (Figure 3; 1.a), then again from near Kennel, SK to the fishing lakes around Fort Qu'Appelle (Figure 3; 1.b), and again for a short stretch at the outlet of Katepwa Lake (Figure 3; 1.c).Geomorphological characteristics associated with Typology 1 are high bank slopes that are narrow with low sinuosity and low fractal dimension (Table 1).The channelized areas at the onset of the river are dominated by the red typology.It has the lowest percent of the sensitive EPT individuals, at only 6% of the total population, yet houses the highest total population at 74% of the sample (Figure 7).Correspondingly, this typology has a high taxa richness of 75, and merely 25 of them are EPT genera.Further examination of the expansive Typology 1 population at the genus level reveals that it is dominated by the amphipod Hyalella.This amphipod genus represents 43% of the total population held in Typology 1.The Hyalella here mark the highest portion represented by one genus in any of the typologies (Figure 9).Following Hyalella, by population size, in the red typology is the black fly larva Simulium (22%), and those only identified to the Amphipoda order (16%).Simulium can be found in almost any lotic water body, as they have many mechanisms to adapt to selective pressures [53,54].Amphipods are very tolerant as well; Hyalella in particular has been found in extremely different climates, from Guatemala to the Northwest Territories [55].Out of all four typologies, Typology 1 appears to have the most imbalanced distribution and least complex community, dominated by the tolerant amphipods and black fly larvae (Figure 9).Typology 1 also consists mostly of the collector gatherer functional feeding group (69%), followed by filterers (24%), and scrapers (4%), with very few representatives of the other groups relative to its large population (Figure 8).As compared with the other typologies, Typology 1 contains the majority of the collector gatherers (81%) and filterers (84%), as well as the most predators (44%), omnivores (37%), and scrapers (35%) of any typology.This is not surprising due to the large population held within this typology, as compared to the other three.The second highest EPT population is found here; nevertheless, it has the lowest % EPT (Figure 7).It also contains few of the four genera found to be statistically linked to typology choice, indicating these specialized genera prefer geomorphic regions different to that of Typology 1 (Figure 6).In combining all of these indices, Typology 1 seems to house the least healthy macroinvertebrate population.It is dominated by less sensitive genera, which could result from nearby anthropogenic influences on the geomorphology.

Typology 2 (Yellow)
Typology 2 dominates the most in two separate stretches, first south of Kenlis, SK and again at Hyde, SK (Figure 3; 2.a; 2.b) It has high sinuosity and fractal dimension, with a negative relation to width, and no relation to slope (Figure 3; Table 1).Although Typology 2 houses the smallest population at merely 2% of the total population (and only five sample sites), it has the highest % of EPT individuals.The intolerant EPT orders represent a notable 52% of the Typology 2 population, indicating healthy water and habitat quality (Figure 7).At the genus level, the yellow typology consists mostly of Stenacron (mayfly; 15%) and Hexagenia (mayfly; 9%), followed by 7% each of Brachycercus (caddisfly), Chironomidae (midge), Baetidae (mayfly), and Probythinella (snail).With a fairly proportional distribution of the remaining genera, Typology 2 appears to have the most even taxa makeup of all four typologies; no single genus takes obvious dominance over the others, as is the case with Hyalella in Typology 1 (Figure 9).It contains predominantly collector gatherers (39%), scrapers (29%), and filterers (16%).Within its small population, Typology 2 has 25% of the total omnivores (Figure 8).

Typology 3 (Green)
The green typology is characterized by large stream width (Table 1), and correspondingly consists mostly of the large water bodies (Buffalo Pound Lake, Pasqua Lake, Echo Lake, and Mission Lake) as well as areas scattered throughout the system (Figure 3; 3.a; 3.b).Typology 3 contains 12% of the total population within its five sampling sites, 18% of which are EPT.When examining Typology 3 at the genus level, it contains mostly the amphipod Hyalella (37%), followed by a fairly proportional distribution with Hydropsyche (caddisfly) and Physa (snail) at 5%, and Valvata (snail), Chironomidae (midge), Corixidae (water boatmen), Baetidae (mayfly), and Oligochaeta (worms) at 4% each.Besides a large amount of Hyalella, this typology has one of the more even taxa distributions, second to Typology 2 (Figure 9).In regards to functional feeding groups, the green typology has mostly collector gatherers (44%), followed by scrapers (23%), and herbivores (13%; Figure 8).When compared to the other typologies, Typology 3 contains the most collector gatherer scrapers (74%), shredders (46%), and herbivores (37%); it also has the least omnivores (4%).Two of the four genera that were revealed to be significantly linked to the typologies (Dubiraphia; Sigara) shared population distribution attributes, with the highest amount of both found in Typology 0, followed by Typology 3 (Figure 6).The water boatman Sigara has very similar numbers within both typologies, indicating it prefers both Typology 0 and 3 (446 individuals were collected in Typology 0, and a close 406 in Typology 3).As herbivores, Sigara need access to plants, and as swimmers, they would do well in the wide streams found in Typology 3. The riffle beetles Dubiraphia more heavily favour Typology 0 (74 individuals), with Typology 3 holding second place with 35 individuals.

Conclusions
A need for rapid and visual methods of river assessment has been identified, especially those capable of reflecting and adjusting to local stream variability [56].Within current systems of understanding fluvial macroinvertebrates, there is a lack of comprehensive watershed approaches.The geospatial model here fills the gap in understanding river typologies and benthic macroinvertebrate distributions.The use of typologies in this study resulted in the successful classification the Qu'Appelle River by similarities in geomorphology and provided a useful model to be used as a reference for these geomorphic characteristics.Geospatial models have the potential to be applied to larger rivers, and to be a predictive method for macroinvertebrate habitat.The geospatial model here was further used as a means to reveal relationships between macroinvertebrate distribution and geomorphology.The effectiveness of geomorphology as a potential means of identifying macroinvertebrate habitat in a more efficient manner has been proven in this way.Certain genera in the Qu'Appelle River are significantly linked to specific regions with similar geomorphic characteristics, showing potential for the use of this geospatial model as a means of pinpointing where certain macroinvertebrates can be found.This research has been an important step in revealing significant relationships between geomorphology and benthic macroinvertebrate community in a fluvial system.We recognize that additional data would help the analysis and complement the results discussed in this study.Whereas in this study, the geospatial data was extracted at a very fine resolution of 50 m intervals along the river, faunistic sampling was carried out only at certain stations.Future studies will include not only geomorphological variables but also biological and physical variables extracted using different methods.A water quality model of the Qu'Appelle River is currently being developed, which will provide water temperatures and dissolved oxygen concentrations at 50 m intervals along the river.Bathymetric surveys of the river are also planned, which will provide water depths and substrate types at the same fine-scale resolution.
The model would be valuable to use prior to, or in lieu of field data collection.Through gaining understanding into the basic geomorphic makeup of the river before sampling, choosing sampling sites can be a more informed process.Since geomorphology drives sediment and contaminant transport, geomorphic typology maps could be a valuable tool to strategically locate new monitoring stations.Macroinvertebrates are established as indicators of environmental health, because they often respond more quickly than other groups, like fish, to ecosystem degradation and contaminants [17,18].Understanding how geomorphology influences macroinvertebrates can be applied further up the food chain as well, as energy and contaminants move bottom up.Geomorphic structure and instream habitats are key limiting factors for freshwater fish, but it can be hard to pinpoint habitat preferences throughout their life stages [57].Macroinvertebrates are important food sources for fish, and are also heavily influenced by geomorphic variables; incorporating macroinvertebrates may help to refine fish studies [2,57].The geomorphic model used in this study is valuable for the sustainable preservation and management of rivers, which requires an understanding of the geomorphological factors that drive species distributions.

Figure 1 .
Figure 1.Location of the Qu'Appelle River, the river centerline, and the 35 macroinvertebrate sampling sites.

Figure 1 .
Figure 1.Location of the Qu'Appelle River, the river centerline, and the 35 macroinvertebrate sampling sites.

Figure 2 .
Figure 2. Delineation of a river polygon, centerline, points placed on the line every 50 m, and transects through each point that cross the polygon to both banks.

Figure 2 .
Figure 2. Delineation of a river polygon, centerline, points placed on the line every 50 m, and transects through each point that cross the polygon to both banks.

Figure 3 .
Figure 3. Delineated Typologies on the Qu'Appelle River centerline, key reaches labeled, site pictures corresponding to the key reaches.

Figure 3 .
Figure 3. Delineated Typologies on the Qu'Appelle River centerline, key reaches labeled, site pictures corresponding to the key reaches.

Figure 5 .
Figure 5. Density plot of the normalized values of channel sinuosity, slope, fractal dimension and width exemplified in Typology 3.

Figure 5 .
Figure 5. Density plot of the normalized values of channel sinuosity, slope, fractal dimension and width exemplified in Typology 3.

Figure 6 .
Figure 6.Bar charts depicting population distribution for the four genera found significant in Table3, as related to the typologies.

Figure 6 .
Figure 6.Bar charts depicting population distribution for the four genera found significant in Table3, as related to the typologies.

Figure 6 .
Figure 6.Bar charts depicting population distribution for the four genera found significant in Table3, as related to the typologies.

Figure 8 .
Figure 8. Proportional bar charts depicting (A) functional feeding groups contained in each typology; and (B) functional feeding group distribution amongst typologies.

Figure 8 .
Figure 8. Proportional bar charts depicting (A) functional feeding groups contained in each typology; and (B) functional feeding group distribution amongst typologies.

Figure 9 .
Figure 9. Population breakdown in each typology, with the most abundant genera (≥5% of the population) named in the legend.

Figure 9 .
Figure 9. Population breakdown in each typology, with the most abundant genera (ě5% of the population) named in the legend.
were analyzed in relation to a geomorphological model of the Qu'Appelle River, Saskatchewan. data

Table 2 .
Kruskal-Wallis Test results for the 50 common genera, with H values, degrees of freedom (df), and p values.
Note: * Significant results are denoted with an asterisk.

Table 3 .
Post hoc Dunn test results for each typology comparison for each genus that had a significant Kruskal-Wallis result, including the Genus Identifier, Z-Value (ZV), p Value (PV), Adjusted PV (APV).
Note: * Significant results are denoted with an asterisk.