Alluvial Floodplain Classification and Organization in Low-Relief Glacially Conditioned River Catchments
by
Roger Thomas James Phillips
A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Department of Geography University of Toronto
© Copyright by Roger Phillips, 2014 Alluvial Floodplain Classification and Organization in Low-Relief Glacially Conditioned River Catchments
Roger Thomas James Phillips
Doctor of Philosophy
Department of Geography University of Toronto
2014 Abstract
The imprint of late Pleistocene glaciation on river systems is an essential theme in Canadian geomorphology. Existing ideas about glacial legacy effects tend to focus on mountainous environments, which are different from the low‐relief physiography of the Laurentian Great
Lakes region. This study investigates river landforms in southern Ontario to develop an improved conceptual fluvial landscape model, reflecting glacial legacy and post‐glacial fluvial processes. The analysis is based on an original dataset, including basic channel observations from over 500 field sites, alluvial floodplain properties from 109 sites, and published physiographic mapping from digital sources.
Glacial signatures are evident in river profiles extracted from a digital elevation model (DEM) for 22 river catchments in southern Ontario. Stream power and slope–area analysis stratify river slopes by glacial landform types and demonstrate significant differences between rivers incised into glacial moraines versus plains. A stream length–gradient index provides a relative measure of how river profiles are oversteepened or understeepened by glacial landforms relative to a theoretical graded profile.
ii Four first‐order alluvial floodplain classifications are presented using k‐means clustering analysis. Predictive variables are explored using PCA and discriminant analysis, producing two principal components: (1) stream power‐resistance and (2) floodplain sedimentology (or floodplain sand equivalent, FSE). The river classifications from this study are consistent with previous literature, but special consideration is required for the inherited sources of cobble and sand materials.
General glacial–fluvial landform relationships can broadly be clustered into: 1) topographic and sedimentological glacial legacy effects; 2) landforms resulting from isostatic and lake baselevel change, and 3) superimposed patterns of Holocene fluvial sediment supply. The adapted fluvial landscape model and conceptual framework presented for rivers in low‐relief glacially conditioned landscapes have the potential to enhance interdisciplinary and applied science in the areas of biogeochemistry, geoarcheology, conservation ecology, and environmental water resources management.
iii Acknowledgments
It doesn’t hurt to start with your PhD supervisory committee. Joe Desloges, Sarah Finkelstein,
Bill Gough, Brian Branfireun, and Nick Eyles, thank you greatly for all your academic and intellectual advice over the years. Joe, you have been an extraordinary mentor in science and teaching. I marvel at your work ethic and dedication to students. Sarah, your enthusiasm and thoughtfulness in teaching are inspiring; and thanks to you I will never forget counting diatoms under the microscope for hours in PGB. Bill, you may not have started as a geographer, but I very much admire your passion for geography now. Brian, our pre‐comps discussions of representative elementary areas were compelling and formative. Nick, our tour of southern
Ontario geological sites in 2008 was invaluable to my learning.
I am grateful to André Roy for his positive and constructive external review of this dissertation and for our insightful discussions during the thesis defense. Thanks also to Mike Church for reading parts of this thesis, and specifically for his review and comments on Chapter 3 resulting in considerable enhancements to the final thesis.
Thanks to the many field assistants over the period of 2010–2013, namely Stephanie Mah,
Beata Opalinska, and Joyce Arabian; as well as James Thayer, Jennifer Henshaw and Tina Hui who also donated some of their own field data. Thanks also to the many public agency staff from OMNR and Conservation Authorities who responded to our information requests, including Brynn Upsdell, Kari Jean, and Ross Wilson (ABCA); Peter Dragunas and Tony Difazio
(CCCA); Shannon Wood and David Pybus (SVCA); Glenn Switzer (NVCA); Muriel Andreae and
Rick Batterson (SCRCA); Joe Gordon and Brian Widner (KCCA); and Kent Todd (OMNR).
iv This research was supported by the Natural Sciences and Engineering Research Council of
Canada (NSERC‐CGS graduate scholarship) and research funding from the University of Toronto to J.R. Desloges.
Again to my patient PhD advisor Joe Desloges, my sincerest gratitude for your guidance and support over the last 6 years. As the saying goes, or it goes without saying, I could not have done this without you. The thesis might be a bit long, or I guess in the words of Mike Church, a bit overdone. So thank you for accepting a few literary indulgences.
Oh well, my mother Doris never minded burnt toast. But she also said I could be rather verbose. I suppose she would be proud of me no matter what. Thanks mom.
To my late father Reginald Phillips (1927–2013) I owe my perseverance for learning.
To my wife Kate, thank you for your loving commitment, for your editorial reviews, and for remembering that Snuffbox has an “n” in the name.
Dedication For Kate, Amelie, and Lachlan. Yes Amelie, I am done my work now. Can you find the three little book mice?
v Table of Contents
Abstract ...... ii Acknowledgments ...... iv Table of Contents ...... vi List of Tables ...... viii List of Figures ...... ix List of Appendices ...... xi
Chapter 1 Introduction ...... 1 Research Statement ...... 1 1.1 Glacial Legacy Effects on Fluvial Systems ...... 4 1.2 Fluvial Process and Landform Interactions ...... 9 1.3 The Laurentian Great Lakes Region ...... 17 1.3.1 Modern human impacts ...... 21 1.4 Summary of Study Approach ...... 21 1.4.1 Statement of authorship and publication status ...... 23
Chapter 2 Glacially conditioned specific stream powers in low‐relief river catchments of the southern Laurentian Great Lakes ...... 24 2.1 Introduction ...... 25 2.2 Theoretical Background ...... 26 2.2.1 The graded river concept ...... 26 2.2.2 Specific stream power approach ...... 30 2.3 Regional Setting ...... 33 2.4 Materials and Methods ...... 36 2.5 Specific Stream Power Inputs ...... 41 2.5.1 Discharge regime models ...... 41 2.5.2 Bankfull width regime models ...... 44 2.5.3 DEM longitudinal profile extraction and slope generalization ...... 47 2.6 Results and Discussion ...... 51 2.6.1 Specific stream power mapping ...... 51 2.6.2 Profile analysis and SL/K index...... 54 2.6.3 Slope–area analysis ...... 59 2.6.4 Glacial conditioning of stream power ...... 64 2.7 Conclusions ...... 66
Chapter 3 Alluvial floodplain classification by multivariate clustering and discriminant analysis for low‐relief glacially conditioned river catchments ...... 70 3.1 Introduction ...... 71 3.1.1 Interdisciplinary floodplains ...... 71 3.1.2 Genetic floodplains ...... 73 3.1.3 Hyperdimensional floodplains ...... 75 3.2 Study Area ...... 77 3.3 Data Collection and Methods ...... 80 vi 3.3.1 Floodplain dataset ...... 80 3.3.2 Multivariate analysis ...... 83 3.4 Results ...... 85 3.4.1 Field floodplain classifications ...... 85 3.4.2 K‐means clustering analysis ...... 88 3.4.3 Principal component analysis ...... 89 3.4.4 Discriminant analysis ...... 95 3.5 Discussion ...... 100 3.5.1 Floodplain parsimony...... 100 3.5.2 Floodplain alphabet ...... 104 3.6 Conclusions ...... 108
Chapter 4 Glacial legacy effects on the spatial organization of alluvial floodplain types in the Laurentian Great Lakes region ...... 110 4.1 Introduction ...... 111 4.2 Study Area ...... 112 4.2.1 Post‐glacial baselevel changes ...... 113 4.3 Data and Methods ...... 118 4.3.1 Field and GIS data ...... 119 4.3.2 Isostatic paleo‐DEM and river incision ...... 122 4.3.3 Radiocarbon and OSL ...... 123 4.4 Results and Discussion ...... 125 4.4.1 Isostatic and baselevel variations ...... 125 4.4.2 Reach classification and river profiles ...... 130 4.4.3 Post‐glacial fluvial adjustments ...... 142 4.5 Conclusions ...... 146
Chapter 5 Conclusions ...... 148 Introduction: “The Beginning of the End” ...... 148 5.1 Thesis Questions ...... 149 5.1.1 Thesis question #1: Glacial signatures and river profiles ...... 149 5.1.2 Thesis question #2: Alluvial floodplain classifications ...... 151 5.1.3 Thesis question #3: Glacial legacy and landform organization ...... 154 5.2 Adapted Fluvial Landscape Model ...... 156 5.3 Research Significance ...... 159 5.3.1 Interdisciplinary significance ...... 160 5.3.2 Contribution to applied geoscience ...... 163
References ...... 164
vii List of Tables
Table 1.1: Comparison of fluvial depositional units and floodplain building process ...... 13 Table 1.2: Summary of floodplain classifications by Nanson and Croke (1992) ...... 14 Table 1.3: Summary of authorship and publication status ...... 23
Table 2.1: Area–discharge regime models for southern Ontario ...... 43 Table 2.2: Area–width regime models (bankfull width) for southern Ontario ...... 46 Table 2.3: Profile analysis results for Hack and Flint equations ...... 55
Table 3.1: Statistical transformation and normality tests of 20 reach variables ...... 80 Table 3.2: KMC analysis results based on 9 selected floodplain variables ...... 89 Table 3.3: PCA eigenvalues of the correlation matrix for 12 variables ...... 93 Table 3.4: PCA squared cosines of the 12 variables for the first four principal components ...... 93 Table 3.5: KMC‐DA confusion matrix for jackknife cross‐validation of floodplain types ...... 97 Table 3.6: Select examples of floodplain sites from 12‐variable DA jackknife cross‐validation .. 98 Table 3.7: KMC‐DA percent correct classification matches from jackknife cross‐validation ...... 99 Table 3.8: Three‐variable classification domains for primary floodplain archetypes ...... 104 Table 3.9: Comparison of floodplain classifications from this study with previous classifications ...... 105
Table 4.1: Number of field sites per catchment...... 120 Table 4.2: Summary of floodplain classifications for southern Ontario rivers...... 121 Table 4.3: Radiocarbon samples from southern Ontario alluvial river floodplains for this study ...... 124 Table 4.4: Summary of general landform relationships for glacially conditioned catchments . 143
Table 5.1: Summary of floodplain classifications for southern Ontario rivers ...... 153 Table 5.2: Glacial–fluvial landform relationships for southern Ontario river catchments ...... 155
viii List of Figures
Figure 1.1: Laurentian Great Lakes watershed boundary and study area of southern Ontario .... 3 Figure 1.2: Fluvial landscape models for fluvial landforms and sedimentary cascades ...... 7 Figure 1.3: Aerial images of river channel patterns from Western Canada and southern Ontario ...... 12 Figure 1.4: Generalized water level curves for the three eastern Great Lakes ...... 18 Figure 1.5: Summary of paleo‐hydrography with emphasis on eastern Great Lakes ...... 19
Figure 2.1: Watershed for the Laurentian Great Lakes of North America ...... 34 Figure 2.2: Map of study area showing major river drainages selected in southern Ontario ..... 36 Figure 2.3: Flow chart of GIS and spreadsheet analysis for DEM slope generalization and stream power mapping ...... 38 Figure 2.4: Area–discharge regime model statistical regression of empirical data ...... 43 Figure 2.5: Area–width regime model statistical regression of empirical data ...... 45 Figure 2.6: Example DEM profile from the Ausable River ...... 48 Figure 2.7: Stream power mapping of select major river drainages in southern Ontario ...... 52 Figure 2.8: Histogram of specific stream power results in southern Ontario ...... 54 Figure 2.9: Longitudinal profiles for select rivers in southern Ontario ...... 57 Figure 2.10: Slope–area plots of 146 river reaches in southern Ontario ...... 60 Figure 2.11: Slope–area plots with reach bed material and reach planform classifications ...... 63 Figure 2.12: Generalized patterns of bedload transport based on modeled stream power ...... 65
Figure 3.1: Study area of southern Ontario ...... 77 Figure 3.2: Schematic cross‐sections and mapping of alluvial floodplain types ...... 86 Figure 3.3: Photographs of channels representing four interpreted floodplain types ...... 87 Figure 3.4: PCA correlation circles with variable projections ...... 90 Figure 3.5: PCA ordinations for floodplain types based on first three principal components ..... 95 Figure 3.6: First two factors (F1 vs. F2) and correlation circle for DA test for 12‐variables ...... 97 Figure 3.7: Group‐split cross‐validation DA analysis for 3‐variable model ...... 100 Figure 3.8: Floodplain classifications for low‐relief glacially conditioned river catchments ..... 102 Figure 3.9: Four floodplain classifications with specific stream power versus FSE ...... 103
ix
Figure 4.1: Study area of southern Ontario ...... 113 Figure 4.2: Post‐glacial water‐level curve for Lake Huron basin ...... 117 Figure 4.3: Isoline maps and histograms of isostatic tilting for each study catchment ...... 126 Figure 4.4: Post‐glacial baselevel highstands in Lake Huron basin ...... 128 Figure 4.5: Reach classification for the Ausable River and select tributaries ...... 131 Figure 4.6: Reach classification for the Saugeen River and select tributaries ...... 135 Figure 4.7: Reach classification for the Nottawasaga River and select tributaries ...... 139 Figure 4.8: Radiocarbon results from alluvial floodplain deposits for this study ...... 145
Figure 5.1: Glacial legacy signatures of moraines and plains for rivers in southern Ontario .... 150 Figure 5.2: Adapted fluvial landscape model for low‐relief glacially conditioned catchments . 157
x List of Appendices
Appendix A: Raw Floodplain Dataset and Example Photographs……………………………………………186
Appendix B: Optically Stimulated Luminescence Dating Potential of Quartz for Holocene Alluvial Deposits in the Southern Laurentian Great Lakes Glaciated Region……………………………………….194
Appendix C: Interdisciplinary Conference Presentations 2009–2012……………………………………..205
Appendix C.1: Snuffbox and the Three Bars: Investigating Geomorphological Approaches to Assess the Distribution of Freshwater Mussel Species at Risk in the Lower Ausable River (Latornell 2012)…………………………………………………………………………………………………205
Appendix C.2: The Geomorphology‐Ecology Balance of Designing Headwater Channels (Natural Channel Systems 2010)…………………………………………………………………………………211
Appendix C.3: “Well my heart’s in the Highlands”: Historical Patterns of Specific Stream Power in Urbanized Highland Creek (CAG‐ONT 2012)………………………………………………..216
Appendix C.4: Tightening the River Meander‐Belt: Application of a Topographic Erodible Corridor Concept Using DEM Raster Analysis – A Case Study of Highland Creek, Ontario (AGU 2009)…………………………………………………………………………………………………………………221
xi
Chapter 1 Introduction
1 Research Statement
Rivers are idealized. They are natural phenomena, but are idealized in that our understanding
must balance simplicity with complexity — generalization with nuance. How well do prevailing
scientific theories represent the environmental complexity and geomorphological diversity of
rivers? The glacial legacy of Quaternary ice ages can be a considerable source of discontinuity
imposed on fluvial systems relative to theories of equilibrium and self‐organization. While the
inscription of glaciation is prominent in much of the Canadian landscape, the low‐relief
topography of the southern Laurentian Great Lakes region is a relatively untested environment
with respect to glacial legacy effects on the spatial organization of fluvial processes and river landforms. As a result, fundamental concepts in fluvial geomorphology must often be highly contextualized.
In such environments, applied geoscience in fluvial geomorphology and river engineering often cannot rely on classic theories which idealize river systems. For example, understanding catchment‐scale patterns of sediment supply and transport is essential in stream restoration, channel design, and erosion hazard assessment projects. Similarly, understanding the spatial organization of fluvial processes is important in interdisciplinary ecological research on aquatic
and riparian habitat, now with growing interest for the sake of sensitive and endangered species. In a low‐relief glacially conditioned landscape not only are steep bedrock‐hillslope and mass wasting processes absent in the headwaters, but idealized trends in channel slope and sediment grain‐size are frequently interrupted by inherited glacial topography and sediments, disrupting systematic patterns of fluvial sediment transport.
1
The broad purpose of this thesis is to investigate the nature of fluvial systems and alluvial floodplains in southern Ontario conditioned to varying degrees by glacial landforms, glacial sediments, and post‐glacial landscape history. Fluvial‐floodplain systems of the Laurentian
Great Lakes region, and southern Ontario in particular (Figure 1.1), have developed through a
complex history of environmental changes which are expressed at a variety of spatial and
temporal scales. In paleoenvironmental reconstructions, and potential insights gained with
respect to future environmental change responses, documentation of this complexity is the
necessary point of scientific explanation and prediction. However, in a traditional view of
scientific parsimony, theoretical expectations of geomorphic process‐landform relationships are
also anticipated to yield some explanation of the spatial distribution of contemporary (or late
Holocene) river landforms and genetic floodplains (Nanson and Croke, 1992).
With respect to the low‐relief glacially conditioned river catchments of the southern Laurentian
Great Lakes region, the following thesis questions have been proposed and refined:
1. To what degree are glacial signatures imposed on the longitudinal profiles and stream power distributions of river systems in this low‐relief landscape?
2. Are natural groupings of channel morphology, alluvial landforms, and fluvial process domains in this environment consistent with previous river and floodplain classifications from other environments, and in view of scientific parsimony which environmental variables are most effective at explaining and predicting distinct classes?
3. How are the spatial arrangements of river landforms and process domains spatially organized in the context of glacial landforms and post‐glacial landscape histories, and how have the fluvial systems responded to the glacial legacy over the Holocene?
2
Figure 1.1: Late Pleistocene glaciation of North America, with the limit of the Laurentide Ice Sheet at 18 ka BP as per data from Dyke et al. (2003). Inset map of modern Laurentian Great Lakes watershed boundary and study area of southern Ontario (SO), with Niagara escarpment (NE), Oak Ridges moraine (ORM), and Algonquin Arch, Precambrian basement bedrock (AA).
The rationale for this research arises from the need to adapt existing ideas about glacial legacy effects to the low‐relief physiography and river systems of the Laurentian Great Lakes region — an environment which has received little systematic attention with respect to fluvial landscape
evolution. This research is also important to inform applied geoscience in fluvial
geomorphology, which would benefit from developing a consistent framework to explain
spatial patterns of river and stream morphology in this context.
The theoretical foundation of this thesis is outlined in the following sections of Chapter 1,
starting with a review of some established concepts in our understanding of glacial legacy
effects on fluvial systems. Fundamental theories of fluvial process and landform interactions
are then discussed with respect to river morphology, channel patterns, floodplain
classifications, and fluvial process domains. A brief review of the post‐glacial history of the
Great Lakes region is then provided to set the stage for an investigation of rivers and alluvial floodplains in southern Ontario. The overall study approach is then summarized at the end of
Chapter 1 to layout the general structure of the thesis. 3
1.1 Glacial Legacy Effects on Fluvial Systems
A recurring discourse in fluvial geomorphology is the effect of past glaciations on river systems.
A notable contribution is that of Church and Ryder (1972) who propose the term paraglacial to describe how fluvial processes are conditioned by glaciation. The paraglacial concept emphasizes the increased sediment supply available to fluvial sediment transport processes
following glaciation and consequently the time period required to effectively exhaust glacial sources of debris. More generally, Slaymaker (2009) discusses this in terms of fluvial responses in “disturbed” landscapes and in terms of how glaciated landscapes respond to non‐glacial
conditions. What emerges from this perspective is that most post‐glacial fluvial adjustments
remain incomplete since the last Pleistocene glaciation (Eyles and Kocsis, 1989; Brardinoni and
Hassan, 2006; Collins and Montgomery, 2011; McCleary et al., 2011; Snyder et al., 2013).
A practical question then is what constitutes an effective framework for measuring glacial
legacy effects? A promising answer seems to hinge on equilibrium theory in geomorphology
(Thorn and Welford, 1994), and specifically relies on the idealized concept of graded river
profiles with systematic increases in drainage area and discharge matched by systematic
decreases in channel slope and bed material grain size. Longitudinal variations in river
morphology or channel pattern can then be expected in a downstream succession of fluvial
process‐landform relationships. In terms of equilibrium responses, river profiles and morphologies may in theory adjust to the prevailing physiographic conditions of topography, geology, and climate. While it is recognized that equilibrium states in geomorphic systems are scale‐dependent in time and space (Schumm and Lichty, 1965), the assumption is that an equilibrium river profile, often envisioned with a smooth concave‐up geometric form, provides
4
a conceptual benchmark from which to gauge external factors which complicate fluvial landscape evolution.
Particularly, since the work of Hack (1957), the equilibrium graded river profile has been demonstrated as a useful measure to investigate landscape diversity, specifically with respect to tectonic, lithological, sedimentological, and base level controls on long‐term river profile
adjustments. The focus being on deviations from a theoretically graded equilibrium‐profile as
evidence of underlying non‐fluvial influences in the landscape. This conceptual approach has
also been applied to glacial legacy effects on fluvial systems (Fonstad, 2003; Brardinoni and
Hassan, 2006; Collins and Montgomery, 2011; McCleary et al., 2011). It is expected that systematic downstream trends in fluvial systems may be interrupted by the lingering effects of
glacial landforms, including both erosional and depositional features. Brardinoni and Hassan
(2006) refer to these features as glacial signatures within stream profiles. A similar concept which has been eluded to by others is the idea of fluvial discontinuities (such as in sediment flux), specifically for glacially conditioned environments (Snyder et al., 2013) and more generally (Burchsted et al., 2014). A complication of this approach is that topographic signatures and sedimentological discontinuities in fluvial systems may be the result of any number of complex environmental factors, of which glacial legacy effects may be just one. As
such, the spatial scale and strength of other potential controlling factors (e.g., tectonic uplift,
structural faulting, mass wasting, geologic variability and base level controls) must be carefully
considered for different environmental contexts.
A number of practical measurement and analysis methods have been previously tested with
respect to glacial legacy signatures on river profiles. Based on deviations from expected
5
equilibrium stream profiles, some notable methods include stream power maps (Fonstad,
2003); slope‐drainage area plots (Brardinoni and Hassan, 2006; Collins and Montgomery, 2011); gradient‐index profiles (or the SL/K index; McCleary et al., 2011); and grain‐size prediction
models (Snyder et al., 2013). Effectively, each of these methods detect slope‐profile anomalies
relative to a theoretical equilibrium slope profile. Each of these methods will be discussed
further in Chapter 2.
Glacial signatures in river profiles and discontinuities in fluvial processes then provide a basis to
assess the downstream spatial organization of channel morphology. Reach‐scale process
domains (Montgomery, 1999) provide a convenient spatial unit by which to delineate fluvial
landforms within the context of relic glacial landform effects. For mountainous catchments in
British Columbia (Canada), Brardinoni and Hassan (2006) present data showing how the
topographic signatures of glacial macro‐forms (e.g., cirques, hanging‐valleys) complicate the
sequence of reach‐scale channel morphology (e.g., cascade, step‐pool, riffle‐pool) due to
imposed variations in fluvial processes and patterns of hillslope mass transfer (e.g., colluvial‐
alluvial coupling and decoupling). Investigating the transition between mountainous
headwaters and lowland alluvial rivers in Washington state (USA), Collins and Montgomery
(2011) document the “valley‐scale” imprint of late Pleistocene glacial features on the
organization of alluvial river landforms (e.g., channel patterns, sinuosity, channel width). These
two examples from British Columbia and Washington demonstrate how morphological
sequences of fluvial landforms can be organized by larger‐scale glacial landscape features.
Thus, the incomplete Holocene response of fluvial systems to the topographic and
sedimentological legacy of late Pleistocene glaciation reflects a long‐term disequilibrium in
6
modern catchment‐scale erosion and sediment transport patterns (Church and Slaymaker,
1989). The downstream succession of fluvial landforms can also be discussed in terms of
sediment cascades (Burt and Allison, 2010) defined as the transfer of sediments through the fluvial landscape, including their source, transport, temporary storage (e.g., in alluvial landforms), and ultimate sink. In terms of paraglacial landscape responses, much of the discussion of sediment cascades tends to focus on a mountains–to–marine succession (Figure
1.2A) of sedimentary land systems (Ballantyne, 2002). However over the continental scale of
Pleistocene glaciation in North America (Figure 1.1), less emphasis has been placed on the vast
low‐relief areas of the continental craton where sediment cascades can be considered largely truncated with respect to the ultimate end‐member sources and sinks of sediment. Specifically within the Laurentian Great Lakes region, the extent and complex architecture of glacial drift together with the large inland lake system leads to a truncated moraines–to–lake‐mud succession of sedimentary cascades and fluvial landforms (Figure 1.2B).
Figure 1.2: Fluvial landscape models for fluvial landforms and sedimentary cascades. (A) Classic mountains–to–marine fluvial landform succession (based on a similar diagram of the fluvial system from Charlton, 2009). (B) Truncated moraines–to–lake‐mud fluvial landform succession, this study (dashed curves represent theoretical equilibrium graded profile surface for mountainous catchments).
7
As a general terminology, glacial conditioning of fluvial systems represents the lasting
impression of late Pleistocene glaciation, including the paraglacial legacy of sediment
availability and the inherited glacial signatures on river profiles relative to a theoretical self‐
organized equilibrium state for fluvial landscapes. Ferguson (1981) describes this condition as a
persistent passive disequilibrium. More specifically, the balance between equilibrium‐
disequilibrium timescales for fluvial response is largely dependent on spatial scale (Schumm
and Lichty, 1965). This idea has been acknowledged in the paraglacial literature where it is expected that small catchments may adjust rapidly to post‐glacial conditions, while at the scale
of large river basins and regional landscapes the paraglacial timescale may potentially span
entire interglacial periods (Church and Ryder, 1972; Slaymaker, 2009; Ballantyne, 2002).
However, it also follows that reach‐scale fluvial processes can exhibit quasi‐equilibrium states
within the overall persistent disequilibrium of much slower river profile adjustments. Quasi‐
equilibrium states of rivers were considered on a theoretical basis by Langbein and Leopold
(1964) where they rationalize the competing tendencies of rivers to adjust rapidly within reaches and to respond slowly in river profiles over geologic time. This idea is further demonstrated for post‐glacial fluvial landscapes by Brardinoni and Hassan (2006) and Collins
and Montgomery (2011) where reach‐scale channel morphologies are organized by non‐fluvial
topographic signatures and glacial features at the catchment scale. In other words, fluvial
process‐landform equilibriums (or quasi‐equilibriums) can in theory be nested within a long‐ term disequilibrium state of a river profile. Reach‐scale equilibriums provide a basis to consider
fluvial process and landform interactions in terms of previous literature on river morphology, channel patterns, floodplain classifications, and fluvial process domains.
8
1.2 Fluvial Process and Landform Interactions
The morphology of a natural river channel is the product of interactions between the variable
flow of water within the channel and the potentially erodible boundary materials, with
consideration of upstream supplies of sediment. Fluvial processes then make‐up the
instantaneous physical mechanics of fluid forces in terms of sediment entrainment, transport,
and deposition; as well as the cumulative patterns of sediment movement (and selective transport) which are responsible for molding and building the alluvial landforms of the river channel. As such, fluvial landforms may be considered the product of a suite of fluvial processes. However, what is less clear is the degree to which these process‐landform
interactions can be considered in equilibrium (i.e., in terms of stationary environmental
controls) balanced with the relative imprint of historical contingency (i.e., importance of
preceding landforms and system trajectories; see Section 1.1 discussion of glacial legacy and equilibrium). These concepts and questions are of course fundamental in geomorphology.
River Morphology and Channel Patterns
A common starting point for a discussion of river morphology is the meandering‐braiding
paradigm of Leopold and Wolman (1954) who discriminated river planforms based on only
discharge (Q) and channel slope (S). Although the classic Q‐S planform discrimination‐analysis
provides a compelling theoretical explanation of river channel patterns based on a threshold
stream power (where stream power, Ω, is proportional to QS), it is understood that the
erodibility of the channel boundary, and namely sediment grain sizes (or calibre), is another
essential factor required to constrain channel pattern predictions (e.g., Ferguson, 1987). It is also recognized that multiple‐channel planforms (or multi‐thread) are not all strictly subject to
braiding processes of dynamic bar formation (i.e., below the bankfull elevation). Specifically, 9
the term anabranching as per Nanson and Knighton (1996) evolved from descriptions of
anastomosing rivers as multiple‐channel planforms where channels are separated by semi‐ stable islands which are elevated above the active channel and are flooded less frequently than braid‐bars (e.g., island surfaces are close to or above the bankfull elevation).
The broadly used term anabranching has also been appropriate to describe other transitional types of multi‐thread channel patterns. For example, wandering gravel‐bed rivers of Western
Canada (Church, 1983; Desloges and Church, 1989) exhibit transitional properties of braiding, meandering, and anabranching as paraglacial sediment supplies are gradually being exhausted
since the last glaciation (e.g., Church and Slaymaker, 1989; see discussion Section 1.1).
Although formerly encompassing many types of anabranching, anastomosing rivers are now
generally considered most typical of extremely low‐gradient rivers in aggrading environments such as estuary deltas, mountainous trenches with cross‐valley alluvial fans (Tabata and Hickin,
2003); and areas of isostatic rebound (Smith, 1983) for example.
The over‐simplification of the original Leopold and Wolman (1954) meandering‐braiding threshold has also been addressed with further consideration of boundary erodibility. A notable example is that the braiding threshold in Q‐S for sand beds is considerably lower
compared to gravel‐bed rivers (e.g., Simpson and Smith, 2001). Further, bank strength has
been shown to be an important factor in determining channel stability and associated
adjustments in channel width (w), with consequential changes in specific stream power (ω ~
Ωw‐1) and thus sediment transport potential (Ferguson 1987; Eaton et al., 2004). While only some studies directly consider bank strength, recent treatments of the meandering‐braiding threshold are more sophisticated, with notable examples being from the work of van den Berg
10
(1995) and Kleinhans and van den Berg (2011) based on specific stream power (ω) and bed
material grain size (D). However, Kleinhans and van den Berg (2011) are careful not to
overstate the discrimination of planform types suggesting that their equations are indicators of transitions (‘transitionators’) between distinguishable bar patterns (i.e., scroll‐bars versus
chute‐bars), rather than hard thresholds between meandering and braiding planforms.
Thus from the planform view, environmental factors which govern fluvial processes, such as
stream power (Q‐S) and sediment calibre (D), have been shown to explain variations in channel
pattern. While such process‐landform associations may have some relevance in the context of southern Ontario rivers (Figure 1.3), the low‐relief landscape of the Laurentian Great Lakes
region does considerably constrain patterns of stream power and sediment supply. Effectively,
multiple‐channel planforms are rare and single‐thread rivers dominate the landscape. As such,
investigation of river morphologies in such environments requires more in‐depth consideration
of the alluvial landforms and their sedimentological characteristics.
Alluvial Floodplain Classification
River planform types (or channel patterns) have provided a basis to interpret alluvial floodplain
landforms, with the seminal article by Nanson and Croke (1992) summarizing a framework for
genetic floodplain classifications. However, interpretation of modern and ancient river
deposits also has a long history in geomorphology and geology in terms of fluvial facies models
(e.g., Miall, 1985; Hickin, 1993). Church (2006) highlights the interconnection between channel
morphology and bed material sediment transport by classifying rivers based on Shields number
and relative roughness. These parameters summarize stream competence and the relative
roles of bedload and suspended sediment load which contribute to building alluvial floodplains.
11
Figure 1.3: Google Earth™ aerial images of river channel patterns from Western Canada (left) and southern Ontario (right) based on planform continuum of Church (1992). A) Braiding Sunwapta River, Alberta (Goff and Ashmore, 1994). B) Braiding Lower Maitland River, southern Ontario (43°45'17" N., 81°42'20" W.). C) Wandering gravel‐bed Bella Coola River, British Columbia (Church, 1983; Desloges and Church, 1989). D) Anabranching Grand River, southern Ontario (Croil, 2002). E) Meandering Milk River, Alberta (Simpson and Smith, 2001). F) Meandering Thames River, southern Ontario (Stewart and Desloges, 2013).
12
Miall (1985) presents a detailed framework for interpreting fluvial deposits in terms of building
three‐dimensional stratigraphical models using a suite of what he terms architectural‐elements.
Nanson and Croke (1992) on the other hand offer a simplified list of six alluvial floodplain
formation processes, including three main and three secondary processes. Table 1.1 summarizes the eight architectural‐elements of Miall (1985) and the six floodplain formation processes of Nanson and Croke (1992). While there is overlap in the first four definitions of fluvial deposits, the remaining depositional architectural‐elements of Miall (1985) hint at the diversity of depositional units relative to the more generalized floodplain formation processes
of Nanson and Croke (1992).
The genetic floodplain classifications of Nanson and Croke (1992) are basically constrained by
two factors: 1) specific stream power (ω), which summarizes the available forces of sediment transport; and 2) cohesive versus non‐cohesive boundary materials, which summarizes in very
general terms the erodibility of the alluvial boundary in terms of sediment calibre and bank strength. A condensed version of the Nanson and Croke (1992) classification scheme is provided in Table 1.2.
Table 1.1: Comparison of fluvial depositional units and floodplain building process from key sources in the literature.
Miall (1985) Nanson and Croke (1992) * Architectural‐Elements Floodplain Formation Processes 1. Lateral accretion deposits = 1. Lateral point‐bar accretion 2. Overbank fines = 2. Overbank vertical‐accretion 3. Gravel bars and bedforms ~ 3. Braid‐channel accretion main 4. Channels ~ 4. Abandoned‐channel accretion secondary 5. Sandy bedforms 5. Oblique accretion 6. Laminated sand sheets 6. Counterpoint accretion 7. Foreset macro‐forms 8. Sediment gravity flows
* Most recently revised in Miall (2010).
13
Table 1.2: Summary of floodplain classifications by Nanson and Croke (1992).
Classification Stream Power Description Examples (ω, Wm‐2) Class A > 300 Disequilibrium floodplains which Single‐thread, often confined, High Energy erode in response to extreme events, dominantly vertical accretion types, Non‐Cohesive typically located in steep headwater include catastrophic cut and fill areas where channel migration is floodplains prevented by valley confinement
Class B 10–300 Equilibrium floodplains formed by Braiding, wandering, and laterally Medium Energy regular flow events in relatively migrating meandering rivers ‐2 Non‐Cohesive unconfined valleys Note: meandering, ω =10–60 Wm
Class C < 10 Floodplains formed by regular flow Single‐thread to anastomosing Low Energy events along laterally stable single‐ Note: anastomosing is a low‐energy Cohesive thread or anastomosing low‐gradient form of anabranching (Nanson and channels Knighton, 1996)
Given that the low‐relief landscape of the Laurentian Great Lakes region is anticipated to constrain patterns of stream power and sediment supply, it is expected that high energy Class A floodplains of Nanson and Croke (1992) will not typically be found in southern Ontario (Note: a few exceptions may be found for localized stream reaches along the Niagara escarpment).
Further, the discharge variability, stream power, and sediment supply characteristics of braiding
and wandering rivers (two sub‐orders of Nanson and Croke (1992) Class B floodplains) do not
typically occur within the study area, limiting the potential development of these planforms.
Other than a few known examples of local anabranching (e.g., Croil, 2002), this generally leaves single‐channel (or single‐thread) rivers as the dominant floodplain types in southern Ontario, including laterally active meandering rivers (Nanson and Croke (1992) Class B, suborder B3) and more stable low‐gradient single‐thread channels with cohesive boundaries (Nanson and Croke
(1992) Class C, suborder C1).
14
In terms of floodplain formation processes (or fluvial architectural‐elements), typical Nanson
and Croke Class B3 meandering river floodplains are primarily built of some combination of lateral point‐bar accretion and overbank vertical‐accretion (Table 1.1). For Nanson and Croke
Class C1 floodplains, overbank vertical accretion is the dominant process. Both Class B3 and C1
floodplains of Nanson and Croke (1992) may also contain minor contributions from oblique
accretion (a sort of steep‐angle cohesive lateral accretion) and abandoned‐channel accretion.
However, in the context of glacial legacy effects in southern Ontario, it is unclear to what
degree these floodplain formation processes hold‐up and it is unknown if the inherited
sediments and topography of the landscape have produced any distinct departures from the
Nanson and Croke (1992) classification framework.
Fluvial Process Domains
A reach of river has long been considered an appropriate spatial scale to organize river
landforms as coherent landscape features derived from distinct assemblages of fluvial
processes (Leopold and Wolman, 1957). However, particularly in the context of environmental heterogeneity the concept of river reaches is juxtaposed with the perspective that fluvial systems are multi‐dimensional physical and biological continuums varying gradually in space and time (Montgomery and Buffington, 1997). While these two perspectives may be at odds, each retains its own value in interpreting fluvial landscapes.
As effectively codified in the work of Montgomery (1999), stream reaches have been rationalized not only as a convenient spatial scale to describe river landforms, but also as
conspicuous spatial units characterized by different suites of dominant geomorphic processes
(i.e., fluvial process domains). As such reaches are intended to represent relatively
15
homogeneous units within the reach‐patchwork of a fluvial drainage network. In abstract statistical terms, it is expected that the morphological variability between reaches will be greater than the variability within reaches. Reaches as defined by Montgomery (1999) are
typically at least 10–20 channel widths in length and may generally be delineated by variations
in topography, geology, vegetation, and hydrology (e.g., tributary confluences).
Montgomery (1999) further suggests that because reaches tend to have similar geology and topography (or slopes), that reaches provide a basis to define what he calls lithotopo units.
Slope–drainage area plots for reaches can then provide a first‐order approach to spatially
discriminate between distinct fluvial process domains, which is essentially a means of stratifying reach‐scale morphology based on stream power (cf. Flores et al., 2006). Montgomery (1999) even advocates that process domains provide an effective framework to interpret fluvial responses to environmental complexity and landscape disturbance (e.g., glacial legacy effects,
Collins and Montgomery, 2011).
The reach framework allows for inspection of the spatial arrangement and structure of
individual reach morphologies within a drainage network mosaic of process domains. It follows
that exploring the spatial organization of fluvial process‐landform groups across a landscape such as southern Ontario, and in the context of complex glacial legacy effects, should provide
insights into patterns of post‐glacial fluvial adjustment. Indeed, the late Pleistocene and
Holocene history of fluvial systems in the Laurentian Great Lakes region is complicated by a
variety of environmental changes following glacial retreat.
16
1.3 The Laurentian Great Lakes Region
The Laurentian Great Lakes are a significant landscape inscription of northern hemisphere glaciations. The story of the Great Lakes has been an exercise in paleoenvironmental reconstruction by scientists for more than a century. Since the late Cenozoic initiation of northern hemisphere glaciations (Raymo, 1994), ice has subsequently scoured large areas of the North American continent. In particular, the Great Lake basins have been carved into bedrock strata of the Paleozoic Michigan sedimentary basin, where weaker strata have been preferentially scoured and resistant bedrock dominates existing islands and peninsulas (Larsen and Schaetzl, 2001). Reworking the landscape of previous Pleistocene glaciations, the melting
Laurentide Ice Sheet of the last glacial period has left behind widespread evidence of the pro‐
glacial and post‐glacial lake systems which evolved in its wake (Larsen and Schaetzl, 2001).
The evidence of paleo‐Great Lakes evolution is complicated by dynamic interactions between changing meltwater inputs, changing topography, and changing climate (i.e., precipitation,
temperature, wind) (Teller, 1995). Areas previously scoured by glacial ice lobes provided some
physical constraints for basins of the proto‐Great Lakes, however initial impoundment during
deglaciation was also accentuated by glacial isostatic depression of the earth’s crust (Peltier,
1994). The processes of deglaciation beginning approximately 15 ka BP were not characterized
by a continuous retreat of ice, but rather by fluctuations in the ice margin superimposed on a
millennium scale melting trend (Teller, 1995; Larsen and Schaetzl, 2001). The evolution of the
Great Lakes can be generalized into four broad lake phases, similar to those discussed by Lewis
et al. (2008). These four phases include: glacial lakes; extreme lakes; Nipissing lakes; and
modern lakes (Figure 1.4). The paleo‐lakes with glacial retreat are also illustrated in Figure 1.5.
17
Paleo‐Lake Names 1. Lake Maumee 2. Lake Ypsilanti 3. Lake Whittlesey 4. Lake Iroquois 5. Early Lake Algonquin – Kirkfield 6. Early Lake Erie 7. Early Lake Ontario 8. Lake Algonquin – Main 9. Lake Stanley/Hough 10. Mattawa highstands 11. Nipissing lakes 12. Modern lakes
Figure 1.4: Generalized water level curves for the three eastern Great Lakes (Huron, Erie, Ontario) relative to modern lake elevations (Adapted from Anderson and Lewis, 1985; Barnett, 1985; Eschman and Karrow, 1985; Pengelly et al., 1997; Lewis et al., 2007), including generalized lake phases (Adapted from Lewis et al., 2008) and paleo‐lake names.
During the glacial lake phase (~14–10 ka BP), pro‐glacial lakes were formed within the basins by
ice‐damning along the margins of the ice sheet, with initial drainage south to the Mississippi
basin (Teller, 1995). These drainage connections likely represented important entry routes to
the early Great Lakes for flora and fauna from the south (Lewis et al., 2008).
During the transition from glacial to post‐glacial conditions (~12–7 ka BP), the lake levels in each basin were adjusting to changing topography (i.e., glacial isostatic rebound) and changing climates, while still receiving indirect pulses of meltwater from glacial lakes to the north (Lewis et al., 2008). This extreme lake phase was most significant in the Lake Huron basin which was characterized by dramatic fluctuations in water levels (Lewis et al., 2008). In comparison, the
Lake Erie and Ontario basins were dominated by low water levels during this phase, with slowly
rising lake levels associated with adjustment to changing paleoclimate and differential isostatic
rebound of their outlet sills (Lewis et al., 2008).
18
Figure 1.5: Summary of paleo‐hydrography with emphasis on eastern Great Lakes (Adapted from data by Dyke et al., 2003). Chronology in uncalibrated radiocarbon years before present.
19
The rising Nipissing lake phase (~7–4 ka BP) followed complete deglaciation of the Laurentide
Ice Sheet with continued adjustments to paleoclimate and differential isostatic rebound (Lewis
et al., 2008). Effectively, inferred increases in precipitation and continued uplift in northeastern
areas of the Great Lakes region resulted in rising lake levels and drainage pattern adjustments
(Lewis et al., 2008). In particular, the lower lakes (i.e., Erie and Ontario) began to receive water
from the upper lakes (i.e., Huron, Michigan, and Superior) by drainage through the Port Huron
outlet starting at about 6 ka BP (Lewis et al., 2008). Water levels then slowly lowered to
modern lake levels (~3–0 ka BP) with incision at the Port Huron outlet and a slight reduction in
climate precipitation in the late Holocene (Lewis et al., 2008).
Late Pleistocene deglaciation in the Great Lakes region exposed glaciogenic sediments to renewed fluvial processes of erosion, sedimentation, and alluvial floodplain development over
the Holocene period. Post‐glacial drainage patterns were inherited from physiographic features of glacial erosion and deposition, as well as from antecedent bedrock controls including
structural features (e.g., Algonquin Arch, Figure 1.1) and neotectonic joints (Eyles, 2002). As discussed above in Section 1.1, glacial legacy effects on fluvial systems generally include
paraglacial influences on sediment availability as well as inherited topographic signatures of
glacial landforms in river profiles. Consequently, reach‐scale fluvial processes and landform
relationships may be spatially organized by larger scale glacial features in the landscape.
However, other paleoenvironmental factors have also influenced post‐glacial fluvial
adjustments over the Holocene, including lake level history, paleoclimate variations, and
differential isostatic rebound. Thus for catchments of the Laurentian Great Lakes region, the
spatial arrangement of river landforms is complicated by diverse paleoenvironmental histories,
with many confounding effects potentially expressed at several spatial and temporal scales.
Still, documenting this spatial complexity relative to theoretical concepts of equilibrium and 20
self‐organization in fluvial geomorphology is expected to provide insights into a more generalized framework for understanding fluvial process‐landform associations in low‐relief glacially conditioned landscapes.
1.3.1 Modern human impacts
Consideration of modern human impacts is acknowledged as a necessary component in the study of river landforms in southern Ontario. While glacial features in the landscape are
expected to impose physiographic controls on southern Ontario rivers, the impacts of European settlement starting in the 1700s and 1800s are associated with widespread deforestation, extensive agricultural development, and local river engineering projects. Post‐settlement
alluvium eroded from cleared land surfaces, grade changes from mill and flood‐control dams,
and channel erosion from runoff in urbanized areas have resulted in modern river adjustments.
Even so, it is expected that late Holocene alluvial records can largely be distinguished from modern impacts in river channel and floodplain landforms. Ultimately, modern human impacts will impart additional variability in a regional dataset of river landforms. While these impacts
may confound glacial influences, the persistence of glacial features can still be tested based on
interpretation of pre‐settlement river landforms and alluvial records.
1.4 Summary of Study Approach
The discussion and literature review presented in Chapter 1 provide the theoretical foundation
to investigate fluvial process‐landform relationships in the context of complex glacial legacy
effects for the low‐relief landscape of southern Ontario. The purpose of this thesis is to
investigate the nature of fluvial systems and alluvial floodplains conditioned to varying degrees
by the inherited glacial landforms and sediments; as well as by other post‐glacial environmental
changes. To engage the proposed thesis questions, this dissertation is structured in a series of
21
three papers (i.e., thesis‐by‐papers), ultimately building to synthesize the study conclusions in
the final chapter. As such, each paper (Chapters 2, 3, and 4) will begin with its own abstract and introduction.
Chapter 2 investigates the glacial signatures contained within river profiles extracted from a digital elevation model of southern Ontario. Anomalies in channel slope and stream power are assessed relative to a theoretical graded profile and are compared to glacial landforms using stream power mapping, slope‐area analysis, and a stream length‐gradient index (i.e., an SL/K index).
Chapter 3 tests alluvial floodplain classifications to characterize the morphological and sedimentological variations of single‐channel planforms in southern Ontario. The classifications are intended to represent distinct morphological groups with differing fluvial process domains in the landscape.
Chapter 4 examines the spatial arrangement of morphological groups to interpret post‐ glacial fluvial adjustments in the context of glacial landforms and post‐glacial landscape histories for select catchments and drainage networks within the study area.
The concluding discussion in Chapter 5 presents a theoretical framework for understanding the
spatial organization of fluvial process‐landform relationships in low‐relief glacially conditioned
environments such as southern Ontario, as adapted from the theoretical foundations derived in
other environments and presented in previous literature. It is intended that this conceptual framework will provide an improved approach to explain the geomorphological diversity of rivers and alluvial floodplain landforms in complex glacially inherited landscapes, particularly those with modest relief. The significance of the research will also be discussed, particularly in
terms of its relevance to applied geoscience in fluvial geomorphology and to interdisciplinary research.
22
1.4.1 Statement of authorship and publication status
This thesis is submitted in conformity with the requirements for the degree of Doctor of Philosophy by the author Roger T.J. Phillips. The introductory and concluding chapters 1 and 5 (respectively) have been written by the author for the purposes of this dissertation. The core chapters 2, 3, and 4 make up the substantive research articles of this thesis. As such, this research has been completed in collaboration with, and with support from, the author’s supervising professor Joseph R. Desloges. The authorship contributions and publication status of each research paper (as of June 2014) are outline in Table 1.3 below.
Table 1.3: Summary of authorship and publication status. Research Paper #1 (Chapter 2)
Glacially conditioned specific stream powers in low relief river catchments of the southern Laurentian Great Lakes Authorship Contribution: R.T.J. Phillips (90%): primary author of thesis; project development; field and lab work; data collection and analysis. J.R. Desloges (10%): supervision of project development (theory and methodology); funding and editorial support. Journal: Geomorphology Publication Status: Published February 1, 2014 (Volume 206, pg. 271–287). Available online: October 19, 2013.
Research Paper #2 (Chapter 3) Alluvial floodplain classification by multivariate clustering and discriminant analysis for low‐relief glacially conditioned river catchments Authorship Contribution: R.T.J. Phillips (90%): primary author of thesis; project development; field and lab work; data collection and analysis. J.R. Desloges (10%): supervision of project development (theory and methodology); funding and editorial support. Journal: Earth Surface Processes and Landforms Publication Status: Submitted March 7, 2014 (Comments received May 28, 2014). Revisions pending.
Research Paper #3 (Chapter 4) Glacial legacy effects on the spatial organization of alluvial floodplain types in the Laurentian Great Lakes region Authorship Contribution: R.T.J. Phillips (90%): primary author of thesis; project development; field and lab work; data collection and analysis. J.R. Desloges (10%): supervision of project development (theory and methodology); funding and editorial support. Journal: TBD Publication Status: In preparation.
23
Chapter 2 Glacially conditioned specific stream powers in low-relief river catchments of the southern Laurentian Great Lakes
2 Chapter 2
Abstract: Fluvial systems of the southern Laurentian Great Lakes region are carved into a complex glacial landscape shaped by continental ice and meltwater of the late Pleistocene. These glacially conditioned river catchments are typically small with drainage areas < 104 km2. A 10‐m digital elevation model (DEM) is used to map the spatial distribution of stream gradient for 22 major river catchments of peninsular southern Ontario, which drain to base levels in the lower Great Lakes (Huron, St. Clair, Erie, and Ontario). Raw data from the DEM show stream gradients that exhibit multiscale variance from real and from artifact sources. Based on a vertical slice and multiple‐pass moving‐window averaging approach, slope data are generalised to the river reach scale (1–2 km) as a representative spatial scale for fluvial processes operating over Holocene timescales. Models of specific stream power are then compared with glacial landform and surface geology mapping. Inherited glacial signatures in river slope appear as deviations in a stream length‐gradient index (SL/K index), where river reaches are frequently oversteepened or understeepened. Based on a slope–area analysis, and complementary to theories of channel pattern discrimination, constant stream power curves (with power‐law exponent of ‐0.4) provide a first‐order approach to stratify river reaches in terms of glacial conditioning and expected planform morphologies. However, multiple‐channel planform types are rare and localized in southern Ontario, indicating that oversteepened reaches with high stream powers may often be moderated by (1) sediment calibre, with cobble‐beds from inherited glacial sediments; and/or (2) relative bank strength, with limited channel widening particularly in gravel and sand‐bed channels. Further discrimination of glacially conditioned fluvial process domains will ultimately require consideration of alluvial floodplain characteristics in addition to general observations of river morphology and channel pattern.
24
2.1 Introduction
Fluvial systems of the southern Laurentian Great Lakes region are geologically young, carved
into a landscape of complex glacial drift and subtle topography. A product of Quaternary
glaciations, the Great Lakes were formed by glacial erosion of Precambrian bedrock of the
Canadian Shield and overlying Paleozoic sedimentary strata of the Michigan basin. At the
southern margins, late Pleistocene deglaciation deposited an extensive landscape of glacial
landforms and sediments, including moraines and glaciolacustrine sequences, in the wake of
the melting Laurentide Ice Sheet (Larson and Schaetzl, 2001).
Glacial conditioning of stream longitudinal profiles has been documented by recent studies for
mountainous landscapes in terms of the topographic reorganization of fluvial processes, as well
as the coupling (or decoupling) of fluvial with hillslope and mass wasting processes (Fonstad,
2003; Brardinoni and Hassan, 2006; McCleary et al., 2011). The study of glacially conditioned
fluvial catchments within the geologic context of past continental ice sheets in nonmountainous
regions has received less attention, and any existing research has yet to be synthesized to a
landscape scale. This is particularly true of the southern Laurentian Great Lakes region where
the topography is subtle and the sedimentary architecture of the glacial palimpsest is vast and
complex.
This paper investigates glacial signatures within river profiles of southern Ontario to extend
previous research on glacial conditioning from mountainous landscapes to the low‐relief
catchments of the southern Laurentian Great Lakes. From river longitudinal profiles, the spatial
properties of stream power are mapped within the context of glacial landforms and sediments.
Previous approaches to evaluate environmental controls on river slope, and consequently
25
channel morphology, are tested with respect to the concept of stream power, in terms of
channel patterns (e.g., meandering vs. braiding; Leopold and Wolman, 1957) and fluvial process
domains (Montgomery, 1999). To evaluate glacial signatures, some theory of fluvial profile
evolution must be acknowledged, traditionally defined as the graded river profile.
2.2 Theoretical Background
2.2.1 The graded river concept
The idealized concept of graded rivers is a compelling paradigm with a long history in the
discipline of fluvial geomorphology (see Chorley, 2000, for historical retrospective). While a
graded stream is by tradition envisioned as a geometrical manifestation, a smooth longitudinal
profile concave to the sky, its scientific establishment as a product of physical channel
processes is largely founded on Mackin’s (1948) definition that it represents the slope of transportation — the hydraulic gradient of the river is adjusted to transport the sediment
supplied to it. Thus changes to the quantity and calibre of sediment supply should cause
adjustments in slope produced by channel aggradation or degradation to compensate for the change, giving the concept of grade a strong connection to theories of geomorphic equilibrium and disequilibrium (e.g., Schumm and Lichty, 1965; Thorn and Welford, 1994).
In general terms, the graded river state trends downstream with systematic increases in discharge and decreases in channel slope (and sediment size), so it has also been considered to
be closely associated with spatial distributions of stream power (i.e., the product of discharge
and slope) (Knighton, 1999). Considerable research relies on the concept of grade or at least
on the theoretical expectation of a smooth concave‐up profile, most often with the assumption
that it can be mathematically represented as an exponential curve (Hack, 1957, 1973; Seeber
26
and Gornitz, 1983; Sinha and Parker, 1996; Morris and Williams, 1997; Knighton, 1999; Smith et
al., 2000; Fonstad, 2003; Jain et al., 2006; Goldrick and Bishop, 2007; Barker et al., 2009; Pérez‐
Peña et al., 2009; Gonga‐Saholiariliva et al., 2011; McCleary et al., 2011).
Since the work of Hack (1957, 1973) and Flint (1974), mathematical representations of an ideal concave‐up profile regularly rely on empirical relationships of slope scaled to either river distance (L) or drainage area (Ad). Hack’s (1957) equation for channel slope (S) can be
summarized by
S = k Ln (1)
where k and n are empirically derived constants, and assuming that bed material size is
constant. Hack (1957) found that the case of n = ‐1 provides a useful graded profile index (or
gradient index) for many of the Appalachian rivers he observed. This produces a simplified
version of Eq. (1), where k is simply equal to the product of SL. Assuming that n = ‐1, the
integrated version of Eq. (1) is of the form
H = C ‐ k ln L (2) where H is the channel bed elevation, and k and C are empirically derived constants. This is a
straight line semilog relationship between channel elevation (linear) and channel distance
(logarithmic), which represents an idealized graded profile, and has been branded by some
authors as the Hack profile (Pérez‐Peña et al., 2009; McCleary et al., 2011). For comparison of
river profiles of different lengths, normalization of the semilog profile using the graded river
gradient (K) has shown to be useful for revealing profile anomalies relative to a theoretical SL/K
index (Seeber and Gornitz, 1983; Pérez‐Peña et al., 2009; McCleary et al., 2011). The
normalization factor K being the SL index, but calculated for the entire river profile as
27
(h h ) (3) K = s f ln Lt
where hs is the elevation of the drainage divide, hf is the elevation of the river outlet, and Lt is
the total length of the entire river.
The primary contribution of Flint (1974) was the equation relating channel slope (S) to drainage
area (Ad) based on an empirically derived power‐law:
–θ S = ks Ad (4)
where ks is known as the steepness index and θ as the concavity index (Whipple, 2004; Gonga‐
Saholiariliva et al., 2011). Values of the concavity index (θ) tend to vary between 0.4 and 1,
with the average often considered to be in the range of 0.6 (Flint, 1974; Whipple, 2004; Gonga‐
Saholiariliva et al., 2011). Assuming that discharge (Q) and drainage area (Ad) increase at roughly an equivalent rate, the lower range value of θ ≈ 0.4 from the Flint equation produces a similar relationship to that proposed by Leopold and Wolman (1957) for the slope‐discharge threshold between meandering and braiding channel patterns:
S* = 0.0125 Q –0.44 (5)
where S* is the meandering–braiding threshold slope (Leopold and Wolman, 1957).
Criticism of the graded river concept tends to focus on the questionable universality of the
smooth concave‐up longitudinal profile as the ultimate state of river evolution (e.g., Phillips and
Lutz, 2008; Phillips, 2011). It is essentially a concern with treating the graded river state as the
normative fluvial condition, in the sense that the geometrically concave‐up profile is not
necessarily the most typical (i.e., exceptions may be more common than the rule) or in terms of
a graded slope profile representing an ideal equilibrium steady state for the entire length of a
28
fluvial system (but concepts of equilibrium are scale‐dependent in space and time; e.g.,
Schumm and Lichty, 1965; Hickin, 1983; Phillips, 2011). So, at best the geometric graded river
profile represents a deep property of fluvial systems (Smith et al., 2000), or at worst it
represents an arbitrary benchmark (Phillips et al., 2010), from which expected variations in
channel slope can be gauged against a conceptual standard.
The graded river concept from Mackin (1948), and his slope of transportation, seems most in
tune with the ideas of geomorphic equilibrium, whereby slope is a dependent variable and is easily adjusted to spatial and temporal changes in sediment supply and calibre. However, this
idea becomes troublesome if coupled with the geometric notion of smooth concave‐up profiles, particularly when valley fills, sediment inputs, and bed material sizes are spatially variable and
do not always trend downstream in a systematic way. Profile irregularities that may be defined
as knickpoints, knickzones (Whipple, 2004; Phillips and Lutz, 2008; Phillips et al., 2010; Gonga‐
Saholiariliva et al., 2011), or other convex features embedded within a long profile may not be
in disequilibrium if viewed at the reach scale. At reach scales, slope may be considered an
independent variable over periods of years, centuries, and perhaps even millennia (Schumm
and Lichty, 1965; Hickin, 1983). On the other hand, over geologic timescales of the Holocene or
Quaternary, profile evolution and channel slope adjustments may tend toward some ideal form, regardless of whether or not it can ever or will ever be achieved.
The concept of the graded river profile, insofar as it represents an exponential concave‐up
Hack‐type profile, has been demonstrated as a useful model for interpreting landscape diversity
(Hack, 1973), particularly with respect to tectonic, lithological, sedimentological, glacial legacy, and base level controls on fluvial profile evolution. The assumption being that deviations from
29
the theoretical benchmark profile represent interruptions to a graded condition (idealized by systematic trends in slope, discharge, and sediment size). Thus, the long profile may be used to tease out clues of the underlying environmental controls that complicate fluvial landscape evolution.
2.2.2 Specific stream power approach
In the most general physical terms, power is the rate at which energy is used, or the rate at which work is performed, expressed in units of watts (W). The concept of stream power thus is an expression of the potential for flowing water to perform geomorphic work, specifically in terms of sediment transport rates. As formulated by Bagnold (1966), the potential energy of
water flowing downslope with gravity can provide an expression of the available stream power
in fluvial channels:
= QS (6)
where is the specific weight of water (9792 kg∙m∙m‐3∙s‐2 at 20°C; Yang, 1996), Q is the
discharge (m3s‐1), and S is the channel slope. The total stream power () per unit length of
stream (in units of Wm‐1), alternatively referred to as cross‐sectional stream power (Rhoads,
1987), can also be expressed per unit bed area if divided by the channel width (w):
QS = (7) w
where ω is the specific stream power in units of Wm‐2. Expressed as the average cross‐ sectional stream power per unit width of the channel (or stream power per unit bed area, m2),
specific stream power may also be referred to as the mean stream power (Rhoads, 1987).
30
The introduction of the stream power concept by Bagnold (1966) was explicitly to engage with matters of sediment transport; and in particular, bedload transport was the initial concern.
Mostly proceeding from the work of Shields (1936), an enormous amount of research by engineers and geomorphologists in the last century has focused on the incipient motion of bed material in terms of critical shear stress related to the forces of flowing water acting on the bed
(for recent reviews see: Buffington and Montgomery, 1997; Ferguson, 2005, 2012; Church,
2006; Lamb et al., 2008; Recking, 2009; Parker et al., 2011). However, as reasoned by Eaton and Church (2011), sediment transport rates (or the volume of sediment transported with time) may be more closely related to the bulk transfer of momentum from the fluid to the bed material, which is a stream power phenomenon.
A growing number of recent studies have updated Bagnold’s (1966, 1980) stream power approach, with the notion of critical stream power for bedload sediment transport founded in the parameter of specific stream power (Ferguson 2005, 2012; Petit et al., 2005; Eaton and
Church, 2011; Parker et al., 2011). Parker et al. (2011) proposed a dimensionless form of critical stream power ( c) as it relates to grain size (Di):