UBC Wiki - User contributions [en]

Course:GEOB406/LynnCr

2011-10-18T00:11:43Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]

Here is a photograph of a laminated silt/clay deposit overlain by a coarse bouldery deposit [http://wiki.ubc.ca/File:LynnCrUpper.jpg found adjacent to the trail]

Here is an example of typical sediment deposits in the [http://wiki.ubc.ca/File:UpperLynn.jpg lower part] of the channel.

To get a sense of the range of flows that are common at Lynn Creek, [http://wiki.ubc.ca/File:LynnLowFlow.jpg here] is a photo from the bridge taken at low flow, and [http://wiki.ubc.ca/File:Teaching122.jpg here] is one taken at high flow.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-18T00:10:33Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]

Here is an example of typical sediment deposits in the [http://wiki.ubc.ca/File:UpperLynn.jpg lower part] of the channel.

To get a sense of the range of flows that are common at Lynn Creek, [http://wiki.ubc.ca/File:LynnLowFlow.jpg here] is a photo from the bridge taken at low flow, and [http://wiki.ubc.ca/File:Teaching122.jpg here] is one taken at high flow.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

File:LynnCrUpper.jpg

2011-10-18T00:10:08Z

Bretteaton: A typical view of the deposits in the upper part of the Lynn Creek watershed.

== Summary ==
A typical view of the deposits in the upper part of the Lynn Creek watershed.
== Copyright status: ==
copyright
== Source: ==
Brett Eaton

Course:GEOB406/LynnCr

2011-10-18T00:03:18Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]

Here is an example of typical sediment deposits in the [http://wiki.ubc.ca/File:UpperLynn.jpg lower part] of the channel

To get a sense of the range of flows that are common at Lynn Creek, [http://wiki.ubc.ca/File:LynnLowFlow.jpg here] is a photo from the bridge taken at low flow, and [http://wiki.ubc.ca/File:Teaching122.jpg here] is one taken at high flow.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

File:LynnLowFlow.jpg

2011-10-18T00:03:05Z

Bretteaton: Here is a picture of Lynn Creek from the bridge at low flow

== Summary ==
Here is a picture of Lynn Creek from the bridge at low flow
== Copyright status: ==
copyright
== Source: ==
Brett Eaton

File:Teaching122.jpg

2011-10-18T00:00:38Z

Bretteaton: Lynn Creek from the bridge at flood

== Summary ==
Lynn Creek from the bridge at flood
== Copyright status: ==
copyright
== Source: ==
Brett Eaton

Course:GEOB406/LynnCr

2011-10-17T20:02:20Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]

Here is an example of typical sediment deposits in the [http://wiki.ubc.ca/File:UpperLynn.jpg lower part] of the channel and the [upper part] of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T20:01:55Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]

----

Here is an example of typical sediment deposits in the [http://wiki.ubc.ca/File:UpperLynn.jpg lower part] of the channel and the [upper part] of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

File:UpperLynn.jpg

2011-10-17T20:01:05Z

Bretteaton: Typical channel conditions in Lynn Creek downstream of the bridge, showing stone steps and clusters.

== Summary ==
Typical channel conditions in Lynn Creek downstream of the bridge, showing stone steps and clusters.
== Copyright status: ==
Copyright
== Source: ==
Brett Eaton

Course:GEOB406/LynnCr

2011-10-17T19:59:46Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]

----

Here is an example of typical sediment deposits in the [lower part] of the channel and the [upper part] of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:58:27Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a [[http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg basal till exposure]]

----

Here is an example of typical sediment deposits in the lower part of the channel and the upper part of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:57:52Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg |thumb|alt=Lynn Creek Till]]

----

Here is an example of typical sediment deposits in the lower part of the channel and the upper part of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:57:22Z

Bretteaton: /* Hazard Assessment of Lynn Creek Watershed */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[http:// File:Lynn_Cr_RB_Till.jpg |thumb|alt=Lynn Creek Till]]

----

Here is an example of typical sediment deposits in the lower part of the channel and the upper part of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:56:01Z

Bretteaton: /* Hazard Assessment of Lynn Creek Watershed */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[File:Lynn_Cr_RB_Till.jpg |thumb|alt=Lynn Creek Till]]

----

Here is an example of typical sediment deposits in the lower part of the channel and the upper part of the channel.

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:55:06Z

Bretteaton: /* Hazard Assessment of Lynn Creek Watershed */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[File:Lynn_Cr_RB_Till.jpg |thumb|alt=Lynn Creek Till]]
---

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:53:12Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[File:Lynn_Cr_RB_Till.jpg |thumb|alt=Lynn Creek Till]]

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:52:14Z

Bretteaton: /* Typical Sediment Deposits */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[File:Lynn_Cr_RB_Till.jpg]]

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

Course:GEOB406/LynnCr

2011-10-17T19:51:24Z

Bretteaton: /* Hazard Assessment of Lynn Creek Watershed */

=Hazard Assessment of Lynn Creek Watershed=

Sections have been created for each of the four sites to be assessed. Please add a new section for each additional new topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Typical Sediment Deposits===
There are a range of sediment deposits in the drainage basin. Here is an example of a basal till exposure [[http://wiki.ubc.ca/File:Lynn_Cr_RB_Till.jpg#file]]

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

File:Lynn Cr RB Till.jpg

2011-10-17T19:48:52Z

Bretteaton: This photo was taken downstream of the bridge at the Lynn Headwaters Park entrance, and shows a good exposure of basal till on the right bank of the river.

== Summary ==
This photo was taken downstream of the bridge at the Lynn Headwaters Park entrance, and shows a good exposure of basal till on the right bank of the river.
== Copyright status: ==
Copyright
== Source: ==
Brett Eaton

Course:GEOB406

2011-10-17T19:43:31Z

Bretteaton:

{{Infobox_New_Course

|title=Watershed Geomorphology

|picture=Image:wiki.png

|subject code=GEOB

|course number=406

|section number=101

|instructor=Brett Eaton

|instructor 2=

|instructor 3=

|instructor 4=

|instructor 5=

|email=brett.eaton@ubc.ca

|office=Geography 232

|office hours=

|schedule=

|classroom=

}}


==Course Description==

This is the WIKI site for GEOB 406, a course designed to introduce undergraduate students to the geomorphic processes and hazards active in mountainous regions. It is a hands-on course requiring students to work with real data on laboratory-style assignments.

This WIKI is intended for use by the students in resolving technical difficulties by putting questions out there for others to read and hopefully answer, posting solutions that have been found, and general discussion of the course material.

==Pages for Discussing the Assignments==
The following pages have been created for students to refer to, edit and comment upon.
#Assessing hazards to a proposed development at [[Course:GEOB406/LynnCr | (Lynn Creek)]]

==Useful Links==
It was suggested that a catalogue of useful web links would be a valuable asset. Please add any useful links in the space below.

Course:GEOB406/LynnCr

2011-10-17T19:42:13Z

Bretteaton: /* Hazard Assessment of Lynn Creek Watershed */

Course:GEOB406/LynnCr

2011-10-17T19:40:49Z

Bretteaton: Created page with "=Hazard Assessment of Lynn Creek Watershed= Please add a new section for each topic, and add your contributions to the appropriate topic area. The assignment is found [http://w..."

=Hazard Assessment of Lynn Creek Watershed=

Please add a new section for each topic, and add your contributions to the appropriate topic area.

The assignment is found [http://www.geog.ubc.ca/courses/geob406/Field_Trip_G406.pdf here].

The development proposal is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_extent.pdf here].

An overlay showing the aerial photograph and topographic map is found [http://www.geog.ubc.ca/courses/geob406/images/Lynn_Cr_AP_detail.pdf here].

===Lower Road: Site 1===
This is the section of the proposed road starting at the bridge and extending about half way to the bridge location
----
comments/questions go here
----

===Upper Road: Site 2===
This is the section proposed road extending up to the proposed bridge location
----
comments/questions go here
----

===Proposed Bridge: Site 3===
The location is just upstream of the vegetated island at the mouth of Kennedy Creek.

----
comments/questions go here
----

===Proposed Development: Site 4===
Located on the fan north of Kennedy Creek.

----
comments/questions

----

A tool identifying leap years may be useful:

http://www.advancedmcode.org/isleap.html

Course:GEOB406

2011-10-17T19:31:14Z

Bretteaton: /* Pages for Discussing the Assignments */

{{Infobox_New_Course

|title=Watershed Geomorphology

|picture=Image:wiki.png

|subject code=GEOB

|course number=406

|section number=101

|instructor=Brett Eaton

|instructor 2=

|instructor 3=

|instructor 4=

|instructor 5=

|email=brett.eaton@ubc.ca

|office=Geography 232

|office hours=

|schedule=

|classroom=

}}


==Course Description==

This is the WIKI site for GEOB 406, a course designed to introduce undergraduate students to the geomorphic processes and hazards active in mountainous regions. It is a hands-on course requiring students to work with real data on laboratory-style assignments.

This WIKI is intended for use by the students in resolving technical difficulties by putting questions out there for others to read and hopefully answer, posting solutions that have been found, and general discussion of the course material.

==Pages for Discussing the Assignments==
The following pages have been created for students to refer to, edit and comment upon.
#Analyzing data from the Water Survey of Canada [[Course:GEOB406/LynnCr | (Lynn Creek Hazards)]]

==Useful Links==
It was suggested that a catalogue of useful web links would be a valuable asset.

http://www.math.iitb.ac.in/~pde08/matlab1.pdf

Course:GEOB406

2011-10-17T19:29:11Z

Bretteaton:

{{Infobox_New_Course

|title=Watershed Geomorphology

|picture=Image:wiki.png

|subject code=GEOB

|course number=406

|section number=101

|instructor=Brett Eaton

|instructor 2=

|instructor 3=

|instructor 4=

|instructor 5=

|email=brett.eaton@ubc.ca

|office=Geography 232

|office hours=

|schedule=

|classroom=

}}


==Course Description==

This is the WIKI site for GEOB 406, a course designed to introduce undergraduate students to the geomorphic processes and hazards active in mountainous regions. It is a hands-on course requiring students to work with real data on laboratory-style assignments.

This WIKI is intended for use by the students in resolving technical difficulties by putting questions out there for others to read and hopefully answer, posting solutions that have been found, and general discussion of the course material.

==Pages for Discussing the Assignments==
The following pages have been created for students to refer to, edit and comment upon. One page will be created for each assignment that we will attempt.
#Analyzing data from the Water Survey of Canada [[Course:GEOB406/Assign1 | (Week 1)]]
#Analyzing data from the Water Survey of Canada [[Course:GEOB503/Week_2 | (Week 2)]]

==Useful Links==
It was suggested that a catalogue of useful web links would be a valuable asset.

http://www.math.iitb.ac.in/~pde08/matlab1.pdf

Course:GEOB406

2011-10-17T19:24:16Z

Bretteaton:

Course:GEOB406

2011-10-17T19:19:54Z

Bretteaton:

Course:GEOB406

2011-10-17T19:19:28Z

Bretteaton:

Course:GEOB406

2011-10-17T19:18:59Z

Bretteaton:

Course:GEOB406

2011-10-17T19:12:19Z

Bretteaton:

Course:GEOB406

2011-10-17T19:11:39Z

Bretteaton:

Course:GEOB406

2011-10-17T19:11:06Z

Bretteaton: This page is intended for the use of students in GEOB 406. They should use this space to post questions and discuss ideas related to the course assignments.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:18:15Z

Bretteaton:

==Introduction to Open Channel Hydraulics==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.

==Flow resistance==
It is the friction between the flowing water and the bed and banks of the stream that controls how fast water can flow at a given stream gradient. The depth of the flow matters, too, since the further a molecule of water is from the bed or banks, the less friction it feels and the faster it can move. Since friction takes place everywhere along the wetted perimeter, ''P'', it is most accurate to think of flow depths in terms of the hydraulic radius, ''R'', rather than the mean hydraulic depth, ''d''. The most common means for expressing the relation between velocity, depth, slope and flow resistance is the ''Manning equation'':

<math>
u = \frac {R^{2/3} \cdot S^{1/2}}{n}
</math>

Here, ''n'' is an empirically determined flow resistance parameter called ''Manning's n''. The value of ''n'' increases as the mean particle diameter of sediment on the bed of the river increases (''D50), and it also varies with the channel morphology.

For natural channels that have steep, boulder-covered beds and that carry relatively small discharges, ''n'' can be as high as 0.06 to 0.07. In natural streams with lower gradients, large discharges and beds composed of sand or fine gravel, ''n'' will range from 0.025 to 0.035.

==Fluid shear forces==
===Vertical velocity distribution and shear velocity===
So far, we have addressed the means by with the average velocity is defined and estimated. However, at at single location across the stream, velocity varies systematically with distance from the bed. The variation of velocity with distance from the bed is called a ''velocity profile''. The water flows most slowly near the bed, and increases with distance from the boundary until it reaches the ''mean free velocity''. The region where the velocity changes continuously is called the ''boundary layer'', and this is where shear stresses are generated. In many streams where the flows are relatively shallow, the boundary layer extends right to the surface of the water, and the flow never reaches the mean free velocity.

Within the boundary layer, the variation in velocity follows the ''log law of the wall'':

<math>
u(z)=\frac {u^*}{\kappa} \ln{\frac {z}{z_o}}
</math>

where ''κ'' is von Karman's constant (''κ ≈'' 0.4), ''z'' is distance measured upwards from the bed surface, ''u*'' is the ''shear velocity'', which is related to the diffusion of momentum from the flow towards the bed, and ''zo'' is the height above the bed at which ''u = 0'' m/s. The term ''zo'' is related to the size of the sediment found on the bed and is called the ''roughness height''.

The force that acts on the bed depends on how quickly momentum is being diffused towards the bed. In a flow, forces are transmitted been molecules that slide over each other but which exhibit some molecular attraction or ''stickiness''. Thus the shear force of gravity acting to pull water down the slope is transmitted from molecule to molecule until it is ultimately transferred to the bed of the stream. The size of the force acting on the bed depends on the molecular viscosity of the fluid (i.e. the stickiness) and how fast the molecules are moving past each other (which is essentially the shear velocity, ''u*''). The shear force per unit area of the bed is given by:

<math>
\tau = \rho \cdot u^{*\,2}
</math>

Therefore, by measuring the velocity profile at a given point, we can use the data to estimate the shear velocity, ''u*''. Then, we can estimate the shear stress, ''τ'', from the estimate of the shear velocity. Im this way, we could in theory estimate ''τ'' at a large number of places within the channel for a given flow, and then map the distribution of shear forces acting on the bed of the channel. This is seldom possible, given the difficulty of collecting enough data before flows change. While this kind of detailed information of the forces acting on the channel bed is highly desirable, fluvial geomorphologists often have to settle for estimates of the average shear stress acting on the channel boundary, which are more easily obtained.

===Average boundary shear stress===
If we think of a a length of stream channel down which water is flowing at a constant rate, we can immediately make the following observations using the derrivation of the ''infinite slope equation'' as a guide:

* the force of gravity acting on the is given by <math>F_g=m \cdot g \cdot \sin{\beta} = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}</math>
* this force must be opposed by an equal and opposite force if the mass of water is not accelerating (which we have said it is not)
* the opposing force is the shear stress generated at the interface between the flowing water and the bed, which is equivalent the the average shear stress times the surface area over which shear stress is being generated, so <math>F_s = \tau \cdot P \cdot L</math>

If we set <math>F_g=F_s</math> and then rearrange to isolate ''τ'', we get:

<math>
F_s=F_g
</math>

<math>
\tau \cdot P \cdot L = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}
</math>

<math>
\tau = \frac{A}{P} \cdot \rho \cdot g \cdot \sin{\beta}
</math>

Now if we recall that ''R = A / P'' and ''γ = ρg'' we can further simplify the equation. Also, for low gradients, we know that ''sin(β)'' ≈ ''S'', and since river gradients are usually less than 0.02 m/m, we can write:

<math>
\tau = \gamma \cdot R \cdot S
</math>

remembering, of course, that this is only appropriate for estimating the average shear stress acting along the boundary, and that it does not tell us anything about the distribution of shear stresses across the channel.

==Flow indices==
===Reynolds number===
Most of the potential energy expended in a stream does not end up transporting sediment downstream. Most of it is turned into turbulence within the flow. Turbulent flows only occur if the mean flow velocity (''u'') is high enough, the flow is deep enough and the viscosity of the flow is low enough. The ''Reynolds number'' (''Re'') is a dimensionless ratio that is related to the degree of turbulence in a fluid flow. The equation is:

<math>
Re=\frac{\rho \cdot v \cdot d}{\mu}
</math>

where μ is the fluid viscosity and ''d'' is the flow depth. When ''Re'' < 500, the flow is purely laminar and there is no turbulence. When ''Re'' > 2000 the flow is fully turbulent. In laminar flow, momentum diffuses through the flow as a result of just the molecular viscosity, while in turbulent flow, eddies carry high and low energy packets upwards and downwards in the flow, adding bulk flow momentum transfers to those associated with viscosity.

===Froude number===
In steep, fast flowing streams things called ''hydraulic jumps'' occur. At the jump, the flow transitions from a shallow flow with a very high velocity to a deeper flow with a low velocity. Both upstream and downstream of the jump the ''specific discharge'', ''q'', remains constant so that (''ud'')upstream = (''ud'')downstream. For the upstream section, the velocity is faster than the speed of a ''gravity wave'' moving through the flow (<math>u_g = \sqrt{gd}</math>). A gravity wave is generated by, for example, throwing a stone into a pond, wherein the gravity wave takes the form of the ripples migrating away from the point-of-impact. In fluids, information about the presence of obstacles protruding into the flow (such as a bridge pier or a large boulder) is transmitted through the flow as a gravity wave, so for the upstream section, where the stream velocity is greater than the speed of a gravity wave, the flow does not see obstacles downstream until it crashes into it. This kind of flow is called ''supercritical flow''. Downstream of the jump, the flow is deeper and slower, and the stream velocity is always less than the speed of a gravity wave. In this case, the flowing water feels the presence of an obstacle before it actually makes contact with it. This is called ''subcritical flow''. The equation that is used to determine if the flow is supercritcal or subcrtical is:

<math>
Fr=\frac {v}{\sqrt{gd}}
</math>

where ''Fr'' is the Froude number. For supercritical flow, ''Fr'' > 1, and for subcritical flow, ''Fr'' < 1 (''Fr'' =1 is associated with ''critical'' flow, where the stream velocity and wave speed are exactly equal). In most streams, the flow is usually subcritical.

Interestingly, the occurrence of hydraulic jumps is one phenomenon that can be easily predicted from hydraulic theory. Consider the potential sources of energy for an element of water in a stream: it possesses potential energy (related to its elevation, ''z''), energy due to the pressure force of the overlying water (related to the water depth, ''d'') and due to the motion of the particle (related to ''u''). The total energy for a 1 m by 1 m by 1 m cube of water, moving at velocity ''u'' is:

<math>
E_{total} = \rho \cdot g \cdot z + \rho \cdot g \cdot d + \frac {\rho \cdot v^2}{2}
</math>

If we divide the entire expression by the unit weight of water (''ρg''), then each component can be expressed as a length scale (called `heads', in m):

<math>
H=z+d+\frac{u^2}{g}
</math>

where ''H'' is called the ''total head'', ''z'' is called the elevation head, ''d'' is the pressure head, and <math>u^2/2g</math> is the velocity head.

In the vicinity of a hydraulic jump, we can assume that the change in elevation is relatively small, and so we can ignore the potential energy term. The remaining terms are referred to as the ''specific energy'' (<math>E = d + u^2/2g</math>, in m) of the flow. Over short distances, ''E'' is effectively constant. Now, remember that <math>q=Q/W=d \cdot u</math>. Since ''q'' is also constant (as long as the width of the stream does not change), we can express ''E'' as a function of one variable (either ''u'' or ''d'') and a couple of constants (''g'' and ''q''). In the equation below, we have substituted ''u=q/d'' into the specific energy equation:

<math>
E=d+\frac{q^2}{2gd^2}
</math>

<math>
E-d =\frac{q^2}{2gd^2}
</math>

<math>
(E-d)\cdot d^2 =\frac{q^2}{2g}
</math>

Since this equation is a cubic equation, there are three solutions for a given value of ''E'', two of which are `real' solutions with positive ''d'' values, and the other is a mathematical artifact (''d'' < 0). One of the real solutions corresponds to supercritical flow, having a small ''d'' value and a high ''u'' value (but the same ''q'' and same ''E''). The other solution has a large ''d'' value and a smaller ''u'' (''q'' and ''E'' again are constant). A hydraulic jump is what happens when one hops from one solution (i.e. the supercritical one) to the other. It does not happen when transitioning from the subcritical flow to the supercritical flow. The reason a jump occurs when transitioning from one solution to the other but not when transitioning in the opposite direction as to do with conservation of momentum. In fact, in order to conserve momentum and specific energy at the same time it is necessary that energy be dissipated at the hydraulic jump as turbulence, so in fact ''E'' is smaller downstream of a jump than upstream of it, even if momentum is conserved.

===Reynolds particle number===
Another important flow index expresses the flow conditions near the bed to the size of grains found on the bed. The ''Reynolds particle number'' (''Re*'') is a dimensionless ratio with the same form as that for predicting the Reynolds number, but with different characteristic velocities and length scales. Referring to the equation for the Reynolds number and thinking about the typical velocity profile in the stream, it is clear that since both ''u'' and ''d'' get smaller as we approach the bed, ''Re'' also gets smaller, and since the flow becomes ''laminar'' for ''Re'' < 500, there should be a very thin layer of laminar-flowing fluid at the boundary: this is called the ''laminar sub-layer'' (δ, in m). The Reynolds particle number considers how fast the velocity is changing near the bed, as measured by the shear velocity, $u*$, which replaces ''u'' in the equation, and it considers only the flow behavior around the particles making up the bed, so the particle diameter (''D'', in m) replaces the flow depth, ''d''. The equation is:

<math>
Re^*= \frac {\rho u^* D}{\mu}
</math>

When ''Re*'' > 70, then ''u*'' is large enough (and ''μ'' is small enough) that the laminar sub-layer is so thin that particles with diameter ''D'' protrude above it, into the turbulent part of the boundary layer, where they are exposed to shear forces due to both molecule-to-molecule transfers of momentum and bulk momentum transfers due to eddies that introduce a vertical component of the flow (i.e. ''D'' > ''δ''). These conditions are referred to as ''hydraulically rough'' flows. When ''Re*'' < 10, then particles of diameter ''D'' are entirely enveloped within the laminar sub-layer, where they are protected from turbulent eddies (i.e. ''D'' < ''δ''). These conditions are referred to as ''hydraulically smooth'' flows. As we will see, ''Re*'' plays and important role in determining the relative ease with which a particle can be ''entrained''.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:17:17Z

Bretteaton: /* Flow indices */

==Stream channel hydraulics & sediment transport==

==Introduction==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.

==Flow resistance==
It is the friction between the flowing water and the bed and banks of the stream that controls how fast water can flow at a given stream gradient. The depth of the flow matters, too, since the further a molecule of water is from the bed or banks, the less friction it feels and the faster it can move. Since friction takes place everywhere along the wetted perimeter, ''P'', it is most accurate to think of flow depths in terms of the hydraulic radius, ''R'', rather than the mean hydraulic depth, ''d''. The most common means for expressing the relation between velocity, depth, slope and flow resistance is the ''Manning equation'':

<math>
u = \frac {R^{2/3} \cdot S^{1/2}}{n}
</math>

Here, ''n'' is an empirically determined flow resistance parameter called ''Manning's n''. The value of ''n'' increases as the mean particle diameter of sediment on the bed of the river increases (''D50), and it also varies with the channel morphology.

For natural channels that have steep, boulder-covered beds and that carry relatively small discharges, ''n'' can be as high as 0.06 to 0.07. In natural streams with lower gradients, large discharges and beds composed of sand or fine gravel, ''n'' will range from 0.025 to 0.035.

==Fluid shear forces==
===Vertical velocity distribution and shear velocity===
So far, we have addressed the means by with the average velocity is defined and estimated. However, at at single location across the stream, velocity varies systematically with distance from the bed. The variation of velocity with distance from the bed is called a ''velocity profile''. The water flows most slowly near the bed, and increases with distance from the boundary until it reaches the ''mean free velocity''. The region where the velocity changes continuously is called the ''boundary layer'', and this is where shear stresses are generated. In many streams where the flows are relatively shallow, the boundary layer extends right to the surface of the water, and the flow never reaches the mean free velocity.

Within the boundary layer, the variation in velocity follows the ''log law of the wall'':

<math>
u(z)=\frac {u^*}{\kappa} \ln{\frac {z}{z_o}}
</math>

where ''κ'' is von Karman's constant (''κ ≈'' 0.4), ''z'' is distance measured upwards from the bed surface, ''u*'' is the ''shear velocity'', which is related to the diffusion of momentum from the flow towards the bed, and ''zo'' is the height above the bed at which ''u = 0'' m/s. The term ''zo'' is related to the size of the sediment found on the bed and is called the ''roughness height''.

The force that acts on the bed depends on how quickly momentum is being diffused towards the bed. In a flow, forces are transmitted been molecules that slide over each other but which exhibit some molecular attraction or ''stickiness''. Thus the shear force of gravity acting to pull water down the slope is transmitted from molecule to molecule until it is ultimately transferred to the bed of the stream. The size of the force acting on the bed depends on the molecular viscosity of the fluid (i.e. the stickiness) and how fast the molecules are moving past each other (which is essentially the shear velocity, ''u*''). The shear force per unit area of the bed is given by:

<math>
\tau = \rho \cdot u^{*\,2}
</math>

Therefore, by measuring the velocity profile at a given point, we can use the data to estimate the shear velocity, ''u*''. Then, we can estimate the shear stress, ''τ'', from the estimate of the shear velocity. Im this way, we could in theory estimate ''τ'' at a large number of places within the channel for a given flow, and then map the distribution of shear forces acting on the bed of the channel. This is seldom possible, given the difficulty of collecting enough data before flows change. While this kind of detailed information of the forces acting on the channel bed is highly desirable, fluvial geomorphologists often have to settle for estimates of the average shear stress acting on the channel boundary, which are more easily obtained.

===Average boundary shear stress===
If we think of a a length of stream channel down which water is flowing at a constant rate, we can immediately make the following observations using the derrivation of the ''infinite slope equation'' as a guide:

* the force of gravity acting on the is given by <math>F_g=m \cdot g \cdot \sin{\beta} = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}</math>
* this force must be opposed by an equal and opposite force if the mass of water is not accelerating (which we have said it is not)
* the opposing force is the shear stress generated at the interface between the flowing water and the bed, which is equivalent the the average shear stress times the surface area over which shear stress is being generated, so <math>F_s = \tau \cdot P \cdot L</math>

If we set <math>F_g=F_s</math> and then rearrange to isolate ''τ'', we get:

<math>
F_s=F_g
</math>

<math>
\tau \cdot P \cdot L = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}
</math>

<math>
\tau = \frac{A}{P} \cdot \rho \cdot g \cdot \sin{\beta}
</math>

Now if we recall that ''R = A / P'' and ''γ = ρg'' we can further simplify the equation. Also, for low gradients, we know that ''sin(β)'' ≈ ''S'', and since river gradients are usually less than 0.02 m/m, we can write:

<math>
\tau = \gamma \cdot R \cdot S
</math>

remembering, of course, that this is only appropriate for estimating the average shear stress acting along the boundary, and that it does not tell us anything about the distribution of shear stresses across the channel.

==Flow indices==
===Reynolds number===
Most of the potential energy expended in a stream does not end up transporting sediment downstream. Most of it is turned into turbulence within the flow. Turbulent flows only occur if the mean flow velocity (''u'') is high enough, the flow is deep enough and the viscosity of the flow is low enough. The ''Reynolds number'' (''Re'') is a dimensionless ratio that is related to the degree of turbulence in a fluid flow. The equation is:

<math>
Re=\frac{\rho \cdot v \cdot d}{\mu}
</math>

where μ is the fluid viscosity and ''d'' is the flow depth. When ''Re'' < 500, the flow is purely laminar and there is no turbulence. When ''Re'' > 2000 the flow is fully turbulent. In laminar flow, momentum diffuses through the flow as a result of just the molecular viscosity, while in turbulent flow, eddies carry high and low energy packets upwards and downwards in the flow, adding bulk flow momentum transfers to those associated with viscosity.

===Froude number===
In steep, fast flowing streams things called ''hydraulic jumps'' occur. At the jump, the flow transitions from a shallow flow with a very high velocity to a deeper flow with a low velocity. Both upstream and downstream of the jump the ''specific discharge'', ''q'', remains constant so that (''ud'')upstream = (''ud'')downstream. For the upstream section, the velocity is faster than the speed of a ''gravity wave'' moving through the flow (<math>u_g = \sqrt{gd}</math>). A gravity wave is generated by, for example, throwing a stone into a pond, wherein the gravity wave takes the form of the ripples migrating away from the point-of-impact. In fluids, information about the presence of obstacles protruding into the flow (such as a bridge pier or a large boulder) is transmitted through the flow as a gravity wave, so for the upstream section, where the stream velocity is greater than the speed of a gravity wave, the flow does not see obstacles downstream until it crashes into it. This kind of flow is called ''supercritical flow''. Downstream of the jump, the flow is deeper and slower, and the stream velocity is always less than the speed of a gravity wave. In this case, the flowing water feels the presence of an obstacle before it actually makes contact with it. This is called ''subcritical flow''. The equation that is used to determine if the flow is supercritcal or subcrtical is:

<math>
Fr=\frac {v}{\sqrt{gd}}
</math>

where ''Fr'' is the Froude number. For supercritical flow, ''Fr'' > 1, and for subcritical flow, ''Fr'' < 1 (''Fr'' =1 is associated with ''critical'' flow, where the stream velocity and wave speed are exactly equal). In most streams, the flow is usually subcritical.

Interestingly, the occurrence of hydraulic jumps is one phenomenon that can be easily predicted from hydraulic theory. Consider the potential sources of energy for an element of water in a stream: it possesses potential energy (related to its elevation, ''z''), energy due to the pressure force of the overlying water (related to the water depth, ''d'') and due to the motion of the particle (related to ''u''). The total energy for a 1 m by 1 m by 1 m cube of water, moving at velocity ''u'' is:

<math>
E_{total} = \rho \cdot g \cdot z + \rho \cdot g \cdot d + \frac {\rho \cdot v^2}{2}
</math>

If we divide the entire expression by the unit weight of water (''ρg''), then each component can be expressed as a length scale (called `heads', in m):

<math>
H=z+d+\frac{u^2}{g}
</math>

where ''H'' is called the ''total head'', ''z'' is called the elevation head, ''d'' is the pressure head, and <math>u^2/2g</math> is the velocity head.

In the vicinity of a hydraulic jump, we can assume that the change in elevation is relatively small, and so we can ignore the potential energy term. The remaining terms are referred to as the ''specific energy'' (<math>E = d + u^2/2g</math>, in m) of the flow. Over short distances, ''E'' is effectively constant. Now, remember that <math>q=Q/W=d \cdot u</math>. Since ''q'' is also constant (as long as the width of the stream does not change), we can express ''E'' as a function of one variable (either ''u'' or ''d'') and a couple of constants (''g'' and ''q''). In the equation below, we have substituted ''u=q/d'' into the specific energy equation:

<math>
E=d+\frac{q^2}{2gd^2}
</math>

<math>
E-d =\frac{q^2}{2gd^2}
</math>

<math>
(E-d)\cdot d^2 =\frac{q^2}{2g}
</math>

Since this equation is a cubic equation, there are three solutions for a given value of ''E'', two of which are `real' solutions with positive ''d'' values, and the other is a mathematical artifact (''d'' < 0). One of the real solutions corresponds to supercritical flow, having a small ''d'' value and a high ''u'' value (but the same ''q'' and same ''E''). The other solution has a large ''d'' value and a smaller ''u'' (''q'' and ''E'' again are constant). A hydraulic jump is what happens when one hops from one solution (i.e. the supercritical one) to the other. It does not happen when transitioning from the subcritical flow to the supercritical flow. The reason a jump occurs when transitioning from one solution to the other but not when transitioning in the opposite direction as to do with conservation of momentum. In fact, in order to conserve momentum and specific energy at the same time it is necessary that energy be dissipated at the hydraulic jump as turbulence, so in fact ''E'' is smaller downstream of a jump than upstream of it, even if momentum is conserved.

===Reynolds particle number===
Another important flow index expresses the flow conditions near the bed to the size of grains found on the bed. The ''Reynolds particle number'' (''Re*'') is a dimensionless ratio with the same form as that for predicting the Reynolds number, but with different characteristic velocities and length scales. Referring to the equation for the Reynolds number and thinking about the typical velocity profile in the stream, it is clear that since both ''u'' and ''d'' get smaller as we approach the bed, ''Re'' also gets smaller, and since the flow becomes ''laminar'' for ''Re'' < 500, there should be a very thin layer of laminar-flowing fluid at the boundary: this is called the ''laminar sub-layer'' (δ, in m). The Reynolds particle number considers how fast the velocity is changing near the bed, as measured by the shear velocity, $u*$, which replaces ''u'' in the equation, and it considers only the flow behavior around the particles making up the bed, so the particle diameter (''D'', in m) replaces the flow depth, ''d''. The equation is:

<math>
Re^*= \frac {\rho u^* D}{\mu}
</math>

When ''Re*'' > 70, then ''u*'' is large enough (and ''μ'' is small enough) that the laminar sub-layer is so thin that particles with diameter ''D'' protrude above it, into the turbulent part of the boundary layer, where they are exposed to shear forces due to both molecule-to-molecule transfers of momentum and bulk momentum transfers due to eddies that introduce a vertical component of the flow (i.e. ''D'' > ''δ''). These conditions are referred to as ''hydraulically rough'' flows. When ''Re*'' < 10, then particles of diameter ''D'' are entirely enveloped within the laminar sub-layer, where they are protected from turbulent eddies (i.e. ''D'' < ''δ''). These conditions are referred to as ''hydraulically smooth'' flows. As we will see, ''Re*'' plays and important role in determining the relative ease with which a particle can be ''entrained''.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:10:50Z

Bretteaton: /* Froude number */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:09:58Z

Bretteaton: /* Froude number */

==Stream channel hydraulics & sediment transport==

==Introduction==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.

==Flow resistance==
It is the friction between the flowing water and the bed and banks of the stream that controls how fast water can flow at a given stream gradient. The depth of the flow matters, too, since the further a molecule of water is from the bed or banks, the less friction it feels and the faster it can move. Since friction takes place everywhere along the wetted perimeter, ''P'', it is most accurate to think of flow depths in terms of the hydraulic radius, ''R'', rather than the mean hydraulic depth, ''d''. The most common means for expressing the relation between velocity, depth, slope and flow resistance is the ''Manning equation'':

<math>
u = \frac {R^{2/3} \cdot S^{1/2}}{n}
</math>

Here, ''n'' is an empirically determined flow resistance parameter called ''Manning's n''. The value of ''n'' increases as the mean particle diameter of sediment on the bed of the river increases (''D50), and it also varies with the channel morphology.

For natural channels that have steep, boulder-covered beds and that carry relatively small discharges, ''n'' can be as high as 0.06 to 0.07. In natural streams with lower gradients, large discharges and beds composed of sand or fine gravel, ''n'' will range from 0.025 to 0.035.

==Fluid shear forces==
===Vertical velocity distribution and shear velocity===
So far, we have addressed the means by with the average velocity is defined and estimated. However, at at single location across the stream, velocity varies systematically with distance from the bed. The variation of velocity with distance from the bed is called a ''velocity profile''. The water flows most slowly near the bed, and increases with distance from the boundary until it reaches the ''mean free velocity''. The region where the velocity changes continuously is called the ''boundary layer'', and this is where shear stresses are generated. In many streams where the flows are relatively shallow, the boundary layer extends right to the surface of the water, and the flow never reaches the mean free velocity.

Within the boundary layer, the variation in velocity follows the ''log law of the wall'':

<math>
u(z)=\frac {u^*}{\kappa} \ln{\frac {z}{z_o}}
</math>

where ''κ'' is von Karman's constant (''κ ≈'' 0.4), ''z'' is distance measured upwards from the bed surface, ''u*'' is the ''shear velocity'', which is related to the diffusion of momentum from the flow towards the bed, and ''zo'' is the height above the bed at which ''u = 0'' m/s. The term ''zo'' is related to the size of the sediment found on the bed and is called the ''roughness height''.

The force that acts on the bed depends on how quickly momentum is being diffused towards the bed. In a flow, forces are transmitted been molecules that slide over each other but which exhibit some molecular attraction or ''stickiness''. Thus the shear force of gravity acting to pull water down the slope is transmitted from molecule to molecule until it is ultimately transferred to the bed of the stream. The size of the force acting on the bed depends on the molecular viscosity of the fluid (i.e. the stickiness) and how fast the molecules are moving past each other (which is essentially the shear velocity, ''u*''). The shear force per unit area of the bed is given by:

<math>
\tau = \rho \cdot u^{*\,2}
</math>

Therefore, by measuring the velocity profile at a given point, we can use the data to estimate the shear velocity, ''u*''. Then, we can estimate the shear stress, ''τ'', from the estimate of the shear velocity. Im this way, we could in theory estimate ''τ'' at a large number of places within the channel for a given flow, and then map the distribution of shear forces acting on the bed of the channel. This is seldom possible, given the difficulty of collecting enough data before flows change. While this kind of detailed information of the forces acting on the channel bed is highly desirable, fluvial geomorphologists often have to settle for estimates of the average shear stress acting on the channel boundary, which are more easily obtained.

===Average boundary shear stress===
If we think of a a length of stream channel down which water is flowing at a constant rate, we can immediately make the following observations using the derrivation of the ''infinite slope equation'' as a guide:

* the force of gravity acting on the is given by <math>F_g=m \cdot g \cdot \sin{\beta} = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}</math>
* this force must be opposed by an equal and opposite force if the mass of water is not accelerating (which we have said it is not)
* the opposing force is the shear stress generated at the interface between the flowing water and the bed, which is equivalent the the average shear stress times the surface area over which shear stress is being generated, so <math>F_s = \tau \cdot P \cdot L</math>

If we set <math>F_g=F_s</math> and then rearrange to isolate ''τ'', we get:

<math>
F_s=F_g
</math>

<math>
\tau \cdot P \cdot L = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}
</math>

<math>
\tau = \frac{A}{P} \cdot \rho \cdot g \cdot \sin{\beta}
</math>

Now if we recall that ''R = A / P'' and ''γ = ρg'' we can further simplify the equation. Also, for low gradients, we know that ''sin(β)'' ≈ ''S'', and since river gradients are usually less than 0.02 m/m, we can write:

<math>
\tau = \gamma \cdot R \cdot S
</math>

remembering, of course, that this is only appropriate for estimating the average shear stress acting along the boundary, and that it does not tell us anything about the distribution of shear stresses across the channel.

==Flow indices==
===Reynolds number===
Most of the potential energy expended in a stream does not end up transporting sediment downstream. Most of it is turned into turbulence within the flow. Turbulent flows only occur if the mean flow velocity (''u'') is high enough, the flow is deep enough and the viscosity of the flow is low enough. The ''Reynolds number'' (''Re'') is a dimensionless ratio that is related to the degree of turbulence in a fluid flow. The equation is:

<math>
Re=\frac{\rho \cdot v \cdot d}{\mu}
</math>

where μ is the fluid viscosity and ''d'' is the flow depth. When ''Re'' < 500, the flow is purely laminar and there is no turbulence. When ''Re'' > 2000 the flow is fully turbulent. In laminar flow, momentum diffuses through the flow as a result of just the molecular viscosity, while in turbulent flow, eddies carry high and low energy packets upwards and downwards in the flow, adding bulk flow momentum transfers to those associated with viscosity.

===Froude number===
In steep, fast flowing streams things called ''hydraulic jumps'' occur. At the jump, the flow transitions from a shallow flow with a very high velocity to a deeper flow with a low velocity. Both upstream and downstream of the jump the ''specific discharge'', ''q'', remains constant so that (''ud'')upstream = (''ud'')downstream. For the upstream section, the velocity is faster than the speed of a ''gravity wave'' moving through the flow (<math>u_g = \sqrt{gd}</math>). A gravity wave is generated by, for example, throwing a stone into a pond, wherein the gravity wave takes the form of the ripples migrating away from the point-of-impact. In fluids, information about the presence of obstacles protruding into the flow (such as a bridge pier or a large boulder) is transmitted through the flow as a gravity wave, so for the upstream section, where the stream velocity is greater than the speed of a gravity wave, the flow does not see obstacles downstream until it crashes into it. This kind of flow is called ''supercritical flow''. Downstream of the jump, the flow is deeper and slower, and the stream velocity is always less than the speed of a gravity wave. In this case, the flowing water feels the presence of an obstacle before it actually makes contact with it. This is called ''subcritical flow''. The equation that is used to determine if the flow is supercritcal or subcrtical is:

<math>
Fr=\frac {v}{\sqrt{gd}}
</math>

where ''Fr'' is the Froude number. For supercritical flow, ''Fr'' > 1, and for subcritical flow, ''Fr'' < 1 (''Fr'' =1 is associated with ''critical'' flow, where the stream velocity and wave speed are exactly equal). In most streams, the flow is usually subcritical.

Interestingly, the occurrence of hydraulic jumps is one phenomenon that can be easily predicted from hydraulic theory. Consider the potential sources of energy for an element of water in a stream: it possesses potential energy (related to its elevation, ''z''), energy due to the pressure force of the overlying water (related to the water depth, ''d'') and due to the motion of the particle (related to ''u''). The total energy for a 1 m by 1 m by 1 m cube of water, moving at velocity ''u'' is:

<math>
E_{total} = \rho \cdot g \cdot z + \rho \cdot g \cdot d + \frac {\rho \cdot v^2}{2}
</math>

If we divide the entire expression by the unit weight of water (''ρg''), then each component can be expressed as a length scale (called `heads', in m):

<math>
H=z+d+\frac{u^2}{g}
</math>

where ''H'' is called the ''total head'', ''z'' is called the elevation head, ''d'' is the pressure head, and <math>u^2/2g</math> is the velocity head.

In the vicinity of a hydraulic jump, we can assume that the change in elevation is relatively small, and so we can ignore the potential energy term. The remaining terms are referred to as the ''specific energy'' (<math>E = d + u^2/2g</math>, in m) of the flow. Over short distances, ''E'' is effectively constant. Now, remember that <math>q=Q/W=d \cdot u</math>. Since ''q'' is also constant (as long as the width of the stream does not change), we can express ''E'' as a function of one variable (either $u$ or $d$) and a couple of constants ($g$ and $q$). In the equation below, we have substituted $u=q/d$ into the specific energy equation:

<math>
E=d+\frac{q^2}{2gd^2}
</math>

<math>
E-d =\frac{q^2}{2gd^2}
</math>

<math>
(E-d)\cdot d^2 =\frac{q^2}{2g}
</math>

Since this equation is a cubic equation, there are three solutions for a given value of ''E'', two of which are `real' solutions with positive ''d'' values, and the other is a mathematical artifact (''d'' < 0). One of the real solutions corresponds to supercritical flow, having a small ''d'' value and a high ''u'' value (but the same ''q'' and same ''E''). The other solution has a large ''d'' value and a smaller ''u'' (''q'' and ''E'' again are constant). A hydraulic jump is what happens when one hops from one solution (i.e. the supercritical one) to the other. It does not happen when transitioning from the subcritical flow to the supercritical flow. The reason a jump occurs when transitioning from one solution to the other but not when transitioning in the opposite direction as to do with conservation of momentum. In fact, in order to conserve momentum and specific energy at the same time it is necessary that energy be dissipated at the hydraulic jump as turbulence, so in fact ''E'' is smaller downstream of a jump than upstream of it, even if momentum is conserved.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:09:41Z

Bretteaton: /* Froude number */

==Stream channel hydraulics & sediment transport==

==Introduction==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.

==Flow resistance==
It is the friction between the flowing water and the bed and banks of the stream that controls how fast water can flow at a given stream gradient. The depth of the flow matters, too, since the further a molecule of water is from the bed or banks, the less friction it feels and the faster it can move. Since friction takes place everywhere along the wetted perimeter, ''P'', it is most accurate to think of flow depths in terms of the hydraulic radius, ''R'', rather than the mean hydraulic depth, ''d''. The most common means for expressing the relation between velocity, depth, slope and flow resistance is the ''Manning equation'':

<math>
u = \frac {R^{2/3} \cdot S^{1/2}}{n}
</math>

Here, ''n'' is an empirically determined flow resistance parameter called ''Manning's n''. The value of ''n'' increases as the mean particle diameter of sediment on the bed of the river increases (''D50), and it also varies with the channel morphology.

For natural channels that have steep, boulder-covered beds and that carry relatively small discharges, ''n'' can be as high as 0.06 to 0.07. In natural streams with lower gradients, large discharges and beds composed of sand or fine gravel, ''n'' will range from 0.025 to 0.035.

==Fluid shear forces==
===Vertical velocity distribution and shear velocity===
So far, we have addressed the means by with the average velocity is defined and estimated. However, at at single location across the stream, velocity varies systematically with distance from the bed. The variation of velocity with distance from the bed is called a ''velocity profile''. The water flows most slowly near the bed, and increases with distance from the boundary until it reaches the ''mean free velocity''. The region where the velocity changes continuously is called the ''boundary layer'', and this is where shear stresses are generated. In many streams where the flows are relatively shallow, the boundary layer extends right to the surface of the water, and the flow never reaches the mean free velocity.

Within the boundary layer, the variation in velocity follows the ''log law of the wall'':

<math>
u(z)=\frac {u^*}{\kappa} \ln{\frac {z}{z_o}}
</math>

where ''κ'' is von Karman's constant (''κ ≈'' 0.4), ''z'' is distance measured upwards from the bed surface, ''u*'' is the ''shear velocity'', which is related to the diffusion of momentum from the flow towards the bed, and ''zo'' is the height above the bed at which ''u = 0'' m/s. The term ''zo'' is related to the size of the sediment found on the bed and is called the ''roughness height''.

The force that acts on the bed depends on how quickly momentum is being diffused towards the bed. In a flow, forces are transmitted been molecules that slide over each other but which exhibit some molecular attraction or ''stickiness''. Thus the shear force of gravity acting to pull water down the slope is transmitted from molecule to molecule until it is ultimately transferred to the bed of the stream. The size of the force acting on the bed depends on the molecular viscosity of the fluid (i.e. the stickiness) and how fast the molecules are moving past each other (which is essentially the shear velocity, ''u*''). The shear force per unit area of the bed is given by:

<math>
\tau = \rho \cdot u^{*\,2}
</math>

Therefore, by measuring the velocity profile at a given point, we can use the data to estimate the shear velocity, ''u*''. Then, we can estimate the shear stress, ''τ'', from the estimate of the shear velocity. Im this way, we could in theory estimate ''τ'' at a large number of places within the channel for a given flow, and then map the distribution of shear forces acting on the bed of the channel. This is seldom possible, given the difficulty of collecting enough data before flows change. While this kind of detailed information of the forces acting on the channel bed is highly desirable, fluvial geomorphologists often have to settle for estimates of the average shear stress acting on the channel boundary, which are more easily obtained.

===Average boundary shear stress===
If we think of a a length of stream channel down which water is flowing at a constant rate, we can immediately make the following observations using the derrivation of the ''infinite slope equation'' as a guide:

* the force of gravity acting on the is given by <math>F_g=m \cdot g \cdot \sin{\beta} = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}</math>
* this force must be opposed by an equal and opposite force if the mass of water is not accelerating (which we have said it is not)
* the opposing force is the shear stress generated at the interface between the flowing water and the bed, which is equivalent the the average shear stress times the surface area over which shear stress is being generated, so <math>F_s = \tau \cdot P \cdot L</math>

If we set <math>F_g=F_s</math> and then rearrange to isolate ''τ'', we get:

<math>
F_s=F_g
</math>

<math>
\tau \cdot P \cdot L = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}
</math>

<math>
\tau = \frac{A}{P} \cdot \rho \cdot g \cdot \sin{\beta}
</math>

Now if we recall that ''R = A / P'' and ''γ = ρg'' we can further simplify the equation. Also, for low gradients, we know that ''sin(β)'' ≈ ''S'', and since river gradients are usually less than 0.02 m/m, we can write:

<math>
\tau = \gamma \cdot R \cdot S
</math>

remembering, of course, that this is only appropriate for estimating the average shear stress acting along the boundary, and that it does not tell us anything about the distribution of shear stresses across the channel.

==Flow indices==
===Reynolds number===
Most of the potential energy expended in a stream does not end up transporting sediment downstream. Most of it is turned into turbulence within the flow. Turbulent flows only occur if the mean flow velocity (''u'') is high enough, the flow is deep enough and the viscosity of the flow is low enough. The ''Reynolds number'' (''Re'') is a dimensionless ratio that is related to the degree of turbulence in a fluid flow. The equation is:

<math>
Re=\frac{\rho \cdot v \cdot d}{\mu}
</math>

where μ is the fluid viscosity and ''d'' is the flow depth. When ''Re'' < 500, the flow is purely laminar and there is no turbulence. When ''Re'' > 2000 the flow is fully turbulent. In laminar flow, momentum diffuses through the flow as a result of just the molecular viscosity, while in turbulent flow, eddies carry high and low energy packets upwards and downwards in the flow, adding bulk flow momentum transfers to those associated with viscosity.

===Froude number===
In steep, fast flowing streams things called ''hydraulic jumps'' occur. At the jump, the flow transitions from a shallow flow with a very high velocity to a deeper flow with a low velocity. Both upstream and downstream of the jump the ''specific discharge'', ''q'', remains constant so that (''ud'')upstream = (''ud'')downstream. For the upstream section, the velocity is faster than the speed of a ''gravity wave'' moving through the flow (<math>u_g = \sqrt{gd}</math>). A gravity wave is generated by, for example, throwing a stone into a pond, wherein the gravity wave takes the form of the ripples migrating away from the point-of-impact. In fluids, information about the presence of obstacles protruding into the flow (such as a bridge pier or a large boulder) is transmitted through the flow as a gravity wave, so for the upstream section, where the stream velocity is greater than the speed of a gravity wave, the flow does not see obstacles downstream until it crashes into it. This kind of flow is called ''supercritical flow''. Downstream of the jump, the flow is deeper and slower, and the stream velocity is always less than the speed of a gravity wave. In this case, the flowing water feels the presence of an obstacle before it actually makes contact with it. This is called ''subcritical flow''. The equation that is used to determine if the flow is supercritcal or subcrtical is:

<math>
Fr=\frac {v}{\sqrt{gd}}
</math>

where ''Fr'' is the Froude number. For supercritical flow, ''Fr'' > 1, and for subcritical flow, ''Fr'' < 1 (''Fr'' =1 is associated with ''critical'' flow, where the stream velocity and wave speed are exactly equal). In most streams, the flow is usually subcritical.

Interestingly, the occurrence of hydraulic jumps is one phenomenon that can be easily predicted from hydraulic theory. Consider the potential sources of energy for an element of water in a stream: it possesses potential energy (related to its elevation, ''z''), energy due to the pressure force of the overlying water (related to the water depth, ''d'') and due to the motion of the particle (related to ''u''). The total energy for a 1 m by 1 m by 1 m cube of water, moving at velocity ''u'' is:

<math>
E_{total} = \rho \cdot g \cdot z + \rho \cdot g \cdot d + \frac {\rho \cdot v^2}{2}
</math>

If we divide the entire expression by the unit weight of water (''ρg''), then each component can be expressed as a length scale (called `heads', in m):

<math>
H=z+d+\frac{u^2}{g}
</math>

where ''H'' is called the ''total head'', ''z'' is called the elevation head, ''d'' is the pressure head, and <math>u^2/2g</math> is the velocity head.

In the vicinity of a hydraulic jump, we can assume that the change in elevation is relatively small, and so we can ignore the potential energy term. The remaining terms are referred to as the ''specific energy'' (<math>E = d + u^2/2g</math>, in m) of the flow. Over short distances, ''E'' is effectively constant. Now, remember that <math>q=Q/W=d \cdot u</math>. Since ''q'' is also constant (as long as the width of the stream does not change), we can express ''E'' as a function of one variable (either $u$ or $d$) and a couple of constants ($g$ and $q$). In the equation below, we have substituted $u=q/d$ into the specific energy equation:
<math>
E=d+\frac{q^2}{2gd^2}
</math>

<math>
E-d =\frac{q^2}{2gd^2}
</math>

<math>
(E-d)\cdot d^2 =\frac{q^2}{2g}
</math>

Since this equation is a cubic equation, there are three solutions for a given value of ''E'', two of which are `real' solutions with positive ''d'' values, and the other is a mathematical artifact (''d'' < 0). One of the real solutions corresponds to supercritical flow, having a small ''d'' value and a high ''u'' value (but the same ''q'' and same ''E''). The other solution has a large ''d'' value and a smaller ''u'' (''q'' and ''E'' again are constant). A hydraulic jump is what happens when one hops from one solution (i.e. the supercritical one) to the other. It does not happen when transitioning from the subcritical flow to the supercritical flow. The reason a jump occurs when transitioning from one solution to the other but not when transitioning in the opposite direction as to do with conservation of momentum. In fact, in order to conserve momentum and specific energy at the same time it is necessary that energy be dissipated at the hydraulic jump as turbulence, so in fact ''E'' is smaller downstream of a jump than upstream of it, even if momentum is conserved.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:09:22Z

Bretteaton: /* Froude number */

==Stream channel hydraulics & sediment transport==

==Introduction==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.

==Flow resistance==
It is the friction between the flowing water and the bed and banks of the stream that controls how fast water can flow at a given stream gradient. The depth of the flow matters, too, since the further a molecule of water is from the bed or banks, the less friction it feels and the faster it can move. Since friction takes place everywhere along the wetted perimeter, ''P'', it is most accurate to think of flow depths in terms of the hydraulic radius, ''R'', rather than the mean hydraulic depth, ''d''. The most common means for expressing the relation between velocity, depth, slope and flow resistance is the ''Manning equation'':

<math>
u = \frac {R^{2/3} \cdot S^{1/2}}{n}
</math>

Here, ''n'' is an empirically determined flow resistance parameter called ''Manning's n''. The value of ''n'' increases as the mean particle diameter of sediment on the bed of the river increases (''D50), and it also varies with the channel morphology.

For natural channels that have steep, boulder-covered beds and that carry relatively small discharges, ''n'' can be as high as 0.06 to 0.07. In natural streams with lower gradients, large discharges and beds composed of sand or fine gravel, ''n'' will range from 0.025 to 0.035.

==Fluid shear forces==
===Vertical velocity distribution and shear velocity===
So far, we have addressed the means by with the average velocity is defined and estimated. However, at at single location across the stream, velocity varies systematically with distance from the bed. The variation of velocity with distance from the bed is called a ''velocity profile''. The water flows most slowly near the bed, and increases with distance from the boundary until it reaches the ''mean free velocity''. The region where the velocity changes continuously is called the ''boundary layer'', and this is where shear stresses are generated. In many streams where the flows are relatively shallow, the boundary layer extends right to the surface of the water, and the flow never reaches the mean free velocity.

Within the boundary layer, the variation in velocity follows the ''log law of the wall'':

<math>
u(z)=\frac {u^*}{\kappa} \ln{\frac {z}{z_o}}
</math>

where ''κ'' is von Karman's constant (''κ ≈'' 0.4), ''z'' is distance measured upwards from the bed surface, ''u*'' is the ''shear velocity'', which is related to the diffusion of momentum from the flow towards the bed, and ''zo'' is the height above the bed at which ''u = 0'' m/s. The term ''zo'' is related to the size of the sediment found on the bed and is called the ''roughness height''.

The force that acts on the bed depends on how quickly momentum is being diffused towards the bed. In a flow, forces are transmitted been molecules that slide over each other but which exhibit some molecular attraction or ''stickiness''. Thus the shear force of gravity acting to pull water down the slope is transmitted from molecule to molecule until it is ultimately transferred to the bed of the stream. The size of the force acting on the bed depends on the molecular viscosity of the fluid (i.e. the stickiness) and how fast the molecules are moving past each other (which is essentially the shear velocity, ''u*''). The shear force per unit area of the bed is given by:

<math>
\tau = \rho \cdot u^{*\,2}
</math>

Therefore, by measuring the velocity profile at a given point, we can use the data to estimate the shear velocity, ''u*''. Then, we can estimate the shear stress, ''τ'', from the estimate of the shear velocity. Im this way, we could in theory estimate ''τ'' at a large number of places within the channel for a given flow, and then map the distribution of shear forces acting on the bed of the channel. This is seldom possible, given the difficulty of collecting enough data before flows change. While this kind of detailed information of the forces acting on the channel bed is highly desirable, fluvial geomorphologists often have to settle for estimates of the average shear stress acting on the channel boundary, which are more easily obtained.

===Average boundary shear stress===
If we think of a a length of stream channel down which water is flowing at a constant rate, we can immediately make the following observations using the derrivation of the ''infinite slope equation'' as a guide:

* the force of gravity acting on the is given by <math>F_g=m \cdot g \cdot \sin{\beta} = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}</math>
* this force must be opposed by an equal and opposite force if the mass of water is not accelerating (which we have said it is not)
* the opposing force is the shear stress generated at the interface between the flowing water and the bed, which is equivalent the the average shear stress times the surface area over which shear stress is being generated, so <math>F_s = \tau \cdot P \cdot L</math>

If we set <math>F_g=F_s</math> and then rearrange to isolate ''τ'', we get:

<math>
F_s=F_g
</math>

<math>
\tau \cdot P \cdot L = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}
</math>

<math>
\tau = \frac{A}{P} \cdot \rho \cdot g \cdot \sin{\beta}
</math>

Now if we recall that ''R = A / P'' and ''γ = ρg'' we can further simplify the equation. Also, for low gradients, we know that ''sin(β)'' ≈ ''S'', and since river gradients are usually less than 0.02 m/m, we can write:

<math>
\tau = \gamma \cdot R \cdot S
</math>

remembering, of course, that this is only appropriate for estimating the average shear stress acting along the boundary, and that it does not tell us anything about the distribution of shear stresses across the channel.

==Flow indices==
===Reynolds number===
Most of the potential energy expended in a stream does not end up transporting sediment downstream. Most of it is turned into turbulence within the flow. Turbulent flows only occur if the mean flow velocity (''u'') is high enough, the flow is deep enough and the viscosity of the flow is low enough. The ''Reynolds number'' (''Re'') is a dimensionless ratio that is related to the degree of turbulence in a fluid flow. The equation is:

<math>
Re=\frac{\rho \cdot v \cdot d}{\mu}
</math>

where μ is the fluid viscosity and ''d'' is the flow depth. When ''Re'' < 500, the flow is purely laminar and there is no turbulence. When ''Re'' > 2000 the flow is fully turbulent. In laminar flow, momentum diffuses through the flow as a result of just the molecular viscosity, while in turbulent flow, eddies carry high and low energy packets upwards and downwards in the flow, adding bulk flow momentum transfers to those associated with viscosity.

===Froude number===
In steep, fast flowing streams things called ''hydraulic jumps'' occur. At the jump, the flow transitions from a shallow flow with a very high velocity to a deeper flow with a low velocity. Both upstream and downstream of the jump the ''specific discharge'', ''q'', remains constant so that (''ud'')upstream = (''ud'')downstream. For the upstream section, the velocity is faster than the speed of a ''gravity wave'' moving through the flow (<math>u_g = \sqrt{gd}</math>). A gravity wave is generated by, for example, throwing a stone into a pond, wherein the gravity wave takes the form of the ripples migrating away from the point-of-impact. In fluids, information about the presence of obstacles protruding into the flow (such as a bridge pier or a large boulder) is transmitted through the flow as a gravity wave, so for the upstream section, where the stream velocity is greater than the speed of a gravity wave, the flow does not see obstacles downstream until it crashes into it. This kind of flow is called ''supercritical flow''. Downstream of the jump, the flow is deeper and slower, and the stream velocity is always less than the speed of a gravity wave. In this case, the flowing water feels the presence of an obstacle before it actually makes contact with it. This is called ''subcritical flow''. The equation that is used to determine if the flow is supercritcal or subcrtical is:

<math>
Fr=\frac {v}{\sqrt{gd}}
</math>

where ''Fr'' is the Froude number. For supercritical flow, ''Fr'' > 1, and for subcritical flow, ''Fr'' < 1 (''Fr'' =1 is associated with ''critical'' flow, where the stream velocity and wave speed are exactly equal). In most streams, the flow is usually subcritical.

Interestingly, the occurrence of hydraulic jumps is one phenomenon that can be easily predicted from hydraulic theory. Consider the potential sources of energy for an element of water in a stream: it possesses potential energy (related to its elevation, ''z''), energy due to the pressure force of the overlying water (related to the water depth, ''d'') and due to the motion of the particle (related to ''u''). The total energy for a 1 m by 1 m by 1 m cube of water, moving at velocity ''u'' is:

<math>
E_{total} = \rho \cdot g \cdot z + \rho \cdot g \cdot d + \frac {\rho \cdot v^2}{2}
</math>

If we divide the entire expression by the unit weight of water (''ρg''), then each component can be expressed as a length scale (called `heads', in m):

<math>
H=z+d+\frac{u^2}{g}
</math>

where ''H'' is called the ''total head'', ''z'' is called the elevation head, ''d'' is the pressure head, and <math>u^2/2g</math> is the velocity head.

In the vicinity of a hydraulic jump, we can assume that the change in elevation is relatively small, and so we can ignore the potential energy term. The remaining terms are referred to as the ''specific energy'' (<math>E = d + u^2/2g</math>, in m) of the flow. Over short distances, ''E'' is effectively constant. Now, remember that <math>q=Q/W=d \cdot u</math>. Since ''q'' is also constant (as long as the width of the stream does not change), we can express ''E'' as a function of one variable (either $u$ or $d$) and a couple of constants ($g$ and $q$). In the equation below, we have substituted $u=q/d$ into the specific energy equation:
<math>
E=d+\frac{q^2}{2gd^2}
</math>
<math>
E-d =\frac{q^2}{2gd^2}
</math>
<math>
(E-d)\cdot d^2 =\frac{q^2}{2g}
</math>

Since this equation is a cubic equation, there are three solutions for a given value of ''E'', two of which are `real' solutions with positive ''d'' values, and the other is a mathematical artifact (''d'' < 0). One of the real solutions corresponds to supercritical flow, having a small ''d'' value and a high ''u'' value (but the same ''q'' and same ''E''). The other solution has a large ''d'' value and a smaller ''u'' (''q'' and ''E'' again are constant). A hydraulic jump is what happens when one hops from one solution (i.e. the supercritical one) to the other. It does not happen when transitioning from the subcritical flow to the supercritical flow. The reason a jump occurs when transitioning from one solution to the other but not when transitioning in the opposite direction as to do with conservation of momentum. In fact, in order to conserve momentum and specific energy at the same time it is necessary that energy be dissipated at the hydraulic jump as turbulence, so in fact ''E'' is smaller downstream of a jump than upstream of it, even if momentum is conserved.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-14T00:08:56Z

Bretteaton: /* Flow indices */

==Stream channel hydraulics & sediment transport==

==Introduction==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.

==Flow resistance==
It is the friction between the flowing water and the bed and banks of the stream that controls how fast water can flow at a given stream gradient. The depth of the flow matters, too, since the further a molecule of water is from the bed or banks, the less friction it feels and the faster it can move. Since friction takes place everywhere along the wetted perimeter, ''P'', it is most accurate to think of flow depths in terms of the hydraulic radius, ''R'', rather than the mean hydraulic depth, ''d''. The most common means for expressing the relation between velocity, depth, slope and flow resistance is the ''Manning equation'':

<math>
u = \frac {R^{2/3} \cdot S^{1/2}}{n}
</math>

Here, ''n'' is an empirically determined flow resistance parameter called ''Manning's n''. The value of ''n'' increases as the mean particle diameter of sediment on the bed of the river increases (''D50), and it also varies with the channel morphology.

For natural channels that have steep, boulder-covered beds and that carry relatively small discharges, ''n'' can be as high as 0.06 to 0.07. In natural streams with lower gradients, large discharges and beds composed of sand or fine gravel, ''n'' will range from 0.025 to 0.035.

==Fluid shear forces==
===Vertical velocity distribution and shear velocity===
So far, we have addressed the means by with the average velocity is defined and estimated. However, at at single location across the stream, velocity varies systematically with distance from the bed. The variation of velocity with distance from the bed is called a ''velocity profile''. The water flows most slowly near the bed, and increases with distance from the boundary until it reaches the ''mean free velocity''. The region where the velocity changes continuously is called the ''boundary layer'', and this is where shear stresses are generated. In many streams where the flows are relatively shallow, the boundary layer extends right to the surface of the water, and the flow never reaches the mean free velocity.

Within the boundary layer, the variation in velocity follows the ''log law of the wall'':

<math>
u(z)=\frac {u^*}{\kappa} \ln{\frac {z}{z_o}}
</math>

where ''κ'' is von Karman's constant (''κ ≈'' 0.4), ''z'' is distance measured upwards from the bed surface, ''u*'' is the ''shear velocity'', which is related to the diffusion of momentum from the flow towards the bed, and ''zo'' is the height above the bed at which ''u = 0'' m/s. The term ''zo'' is related to the size of the sediment found on the bed and is called the ''roughness height''.

The force that acts on the bed depends on how quickly momentum is being diffused towards the bed. In a flow, forces are transmitted been molecules that slide over each other but which exhibit some molecular attraction or ''stickiness''. Thus the shear force of gravity acting to pull water down the slope is transmitted from molecule to molecule until it is ultimately transferred to the bed of the stream. The size of the force acting on the bed depends on the molecular viscosity of the fluid (i.e. the stickiness) and how fast the molecules are moving past each other (which is essentially the shear velocity, ''u*''). The shear force per unit area of the bed is given by:

<math>
\tau = \rho \cdot u^{*\,2}
</math>

Therefore, by measuring the velocity profile at a given point, we can use the data to estimate the shear velocity, ''u*''. Then, we can estimate the shear stress, ''τ'', from the estimate of the shear velocity. Im this way, we could in theory estimate ''τ'' at a large number of places within the channel for a given flow, and then map the distribution of shear forces acting on the bed of the channel. This is seldom possible, given the difficulty of collecting enough data before flows change. While this kind of detailed information of the forces acting on the channel bed is highly desirable, fluvial geomorphologists often have to settle for estimates of the average shear stress acting on the channel boundary, which are more easily obtained.

===Average boundary shear stress===
If we think of a a length of stream channel down which water is flowing at a constant rate, we can immediately make the following observations using the derrivation of the ''infinite slope equation'' as a guide:

* the force of gravity acting on the is given by <math>F_g=m \cdot g \cdot \sin{\beta} = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}</math>
* this force must be opposed by an equal and opposite force if the mass of water is not accelerating (which we have said it is not)
* the opposing force is the shear stress generated at the interface between the flowing water and the bed, which is equivalent the the average shear stress times the surface area over which shear stress is being generated, so <math>F_s = \tau \cdot P \cdot L</math>

If we set <math>F_g=F_s</math> and then rearrange to isolate ''τ'', we get:

<math>
F_s=F_g
</math>

<math>
\tau \cdot P \cdot L = A \cdot L \cdot \rho \cdot g \cdot \sin{\beta}
</math>

<math>
\tau = \frac{A}{P} \cdot \rho \cdot g \cdot \sin{\beta}
</math>

Now if we recall that ''R = A / P'' and ''γ = ρg'' we can further simplify the equation. Also, for low gradients, we know that ''sin(β)'' ≈ ''S'', and since river gradients are usually less than 0.02 m/m, we can write:

<math>
\tau = \gamma \cdot R \cdot S
</math>

remembering, of course, that this is only appropriate for estimating the average shear stress acting along the boundary, and that it does not tell us anything about the distribution of shear stresses across the channel.

==Flow indices==
===Reynolds number===
Most of the potential energy expended in a stream does not end up transporting sediment downstream. Most of it is turned into turbulence within the flow. Turbulent flows only occur if the mean flow velocity (''u'') is high enough, the flow is deep enough and the viscosity of the flow is low enough. The ''Reynolds number'' (''Re'') is a dimensionless ratio that is related to the degree of turbulence in a fluid flow. The equation is:

<math>
Re=\frac{\rho \cdot v \cdot d}{\mu}
</math>

where μ is the fluid viscosity and ''d'' is the flow depth. When ''Re'' < 500, the flow is purely laminar and there is no turbulence. When ''Re'' > 2000 the flow is fully turbulent. In laminar flow, momentum diffuses through the flow as a result of just the molecular viscosity, while in turbulent flow, eddies carry high and low energy packets upwards and downwards in the flow, adding bulk flow momentum transfers to those associated with viscosity.

===Froude number===
In steep, fast flowing streams things called ''hydraulic jumps'' occur. At the jump, the flow transitions from a shallow flow with a very high velocity to a deeper flow with a low velocity. Both upstream and downstream of the jump the ''specific discharge'', ''q'', remains constant so that (''ud'')upstream = (''ud'')downstream. For the upstream section, the velocity is faster than the speed of a ''gravity wave'' moving through the flow (<math>u_g = \sqrt{gd}</math>). A gravity wave is generated by, for example, throwing a stone into a pond, wherein the gravity wave takes the form of the ripples migrating away from the point-of-impact. In fluids, information about the presence of obstacles protruding into the flow (such as a bridge pier or a large boulder) is transmitted through the flow as a gravity wave, so for the upstream section, where the stream velocity is greater than the speed of a gravity wave, the flow does not see obstacles downstream until it crashes into it. This kind of flow is called ''supercritical flow''. Downstream of the jump, the flow is deeper and slower, and the stream velocity is always less than the speed of a gravity wave. In this case, the flowing water feels the presence of an obstacle before it actually makes contact with it. This is called ''subcritical flow''. The equation that is used to determine if the flow is supercritcal or subcrtical is:

<math>
Fr=\frac {v}{\sqrt{gd}}
</math>

where ''Fr'' is the Froude number. For supercritical flow, ''Fr'' > 1, and for subcritical flow, ''Fr'' < 1 (''Fr'' =1 is associated with ''critical'' flow, where the stream velocity and wave speed are exactly equal). In most streams, the flow is usually subcritical.

Interestingly, the occurrence of hydraulic jumps is one phenomenon that can be easily predicted from hydraulic theory. Consider the potential sources of energy for an element of water in a stream: it possesses potential energy (related to its elevation, ''z''), energy due to the pressure force of the overlying water (related to the water depth, ''d'') and due to the motion of the particle (related to ''u''). The total energy for a 1 m by 1 m by 1 m cube of water, moving at velocity ''u'' is:

<math>
E_{total} = \rho \cdot g \cdot z + \rho \cdot g \cdot d + \frac {\rho \cdot v^2}{2}
</math>

If we divide the entire expression by the unit weight of water (''ρg''), then each component can be expressed as a length scale (called `heads', in m):

<math>
H=z+d+\frac{u^2}{g}
</math>

where ''H'' is called the ''total head'', ''z'' is called the elevation head, ''d'' is the pressure head, and <math>u^2/2g</math> is the velocity head.

In the vicinity of a hydraulic jump, we can assume that the change in elevation is relatively small, and so we can ignore the potential energy term. The remaining terms are referred to as the ''specific energy'' (<math>E = d + u^2/2g</math>, in m) of the flow. Over short distances, ''E'' is effectively constant. Now, remember that <math>q=Q/W=d \cdot u</math>. Since ''q'' is also constant (as long as the width of the stream does not change), we can express ''E'' as a function of one variable (either $u$ or $d$) and a couple of constants ($g$ and $q$). In the equation below, we have substituted $u=q/d$ into the specific energy equation:
<math>
E&=d+\frac{q^2}{2gd^2}
</math>
<math>
E-d &=\frac{q^2}{2gd^2}
</math>
<math>
(E-d)\cdot d^2 &=\frac{q^2}{2g}
</math>

Since this equation is a cubic equation, there are three solutions for a given value of ''E'', two of which are `real' solutions with positive ''d'' values, and the other is a mathematical artifact (''d'' < 0). One of the real solutions corresponds to supercritical flow, having a small ''d'' value and a high ''u'' value (but the same ''q'' and same ''E''). The other solution has a large ''d'' value and a smaller ''u'' (''q'' and ''E'' again are constant). A hydraulic jump is what happens when one hops from one solution (i.e. the supercritical one) to the other. It does not happen when transitioning from the subcritical flow to the supercritical flow. The reason a jump occurs when transitioning from one solution to the other but not when transitioning in the opposite direction as to do with conservation of momentum. In fact, in order to conserve momentum and specific energy at the same time it is necessary that energy be dissipated at the hydraulic jump as turbulence, so in fact ''E'' is smaller downstream of a jump than upstream of it, even if momentum is conserved.

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:58:52Z

Bretteaton:

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:56:39Z

Bretteaton: /* {Average boundary shear stress */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:56:18Z

Bretteaton: /* {Average boundary shear stress */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:55:51Z

Bretteaton:

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:49:01Z

Bretteaton: /* Vertical velocity distribution and shear velocity */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:48:38Z

Bretteaton: /* Vertical velocity distribution and shear velocity */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:48:09Z

Bretteaton: /* Vertical velocity distribution and shear velocity */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:47:36Z

Bretteaton: /* Vertical velocity distribution and shear velocity */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:45:46Z

Bretteaton:

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:40:57Z

Bretteaton:

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:38:17Z

Bretteaton: /* Stream channel hydraulics & sediment transport */

Course:GEOBrefGuide/Geomorphology/Hydraulics

2010-12-13T23:37:01Z

Bretteaton: /* Unit stream power */

=Stream channel hydraulics & sediment transport=

==Introduction==
In order to consider fluvial landforms from a process point-of-view, it is first necessary to consider the laws governing the flow of water. These laws are typically covered under the broad heading of ''open channel hydraulics'', and they have been determined through theoretical and empirical means, mostly by hydraulic engineers. There are a number of key concepts, including:

* total stream power (''Ω'')
* unit stream power (''ω'')
* shear stress (''τ'') and shear velocity (''u*''), which are basically the same thing
* critical shear stress (''τc'')
* mean velocity (''u'')
* Froude number (''Fr'')
* Reynolds number (''Re'')
* Reynolds particle number (''Re*'')

We will define and explore all of these terms, as well as the equations used for estimating them. They are the basic building blocks used to understand fluvial geomorphology.

==Continuity==
Before getting into the key concepts, it is useful to define the relation between the geometric variables and stream discharge. The most basic relation is that which relates the mean velocity of the stream (''u'', in m/s) to the total discharge being carried by the stream (''Q'', in m 3/s) using the cross sectional area of the flow (''A'', in m2).

<math>
Q = A \cdot u
</math>

The cross sectional area can be further subdivided into a horizontal and a vertical component, one of which is measured, the other being estimated. If we divide ''A'' by the width of the water surface (''W'', in m), we get an estimate of the average water depth, referred to as the ''mean hydraulic depth'' (''d'', in m). Substituting these variables into the equation above gives:

<math>
Q = W \cdot d \cdot u
</math>

which is referred to as the ''continuity equation''. This is the key equation that constrains many others, since it is really stating that the width, depth, and/or velocity with increase/decrease as discharge increases/decreases.

There is another possible way of subdividing the cross section area of flow: if we divide ''A'' by the total length of the channel bed and banks, starting at one water's edge and moving across the stream to the other (which is called the ''wetted perimeter'', ''P'', in m), we get a quantity that is similar to but slightly less than ''d''. This quantity is called the ''hydraulic radius'' (''R'', in m). We can also use these variables to define a continuity equation:

<math>
Q = P \cdot R \cdot u
</math>

but this version is seldom if ever used by fluvial geomorphologists.

To summarize, the dimensions of the water-filled part of the channel change with discharge; the mean velocity is given by the ratio ''u=Q/A''; and the cross sectional area can be related to either the channel width and depth (''W, d'') or to wetted perimeter and hydraulic radius (''P,R'').

==Stream Power==
The loss of potential energy per unit length of time (''ΔEp'' / Δ''t'') determines the rate at which geomorphic work can be done. This is referred to as ''stream power''. In some cases, geomorphologists use the total stream power that a river has (called the ''total stream power'') and in others it is more useful to consider the stream power acting on a 1 m by 1 m area of the bed of the stream (called the ''unit stream power'').

===Total stream power===
Recall that the loss of potential energy per unit length of stream channel is proportional to the change in elevation.

<math>
\frac{\Delta E_p}{\Delta x} = m \cdot g \cdot \frac{\Delta z}{\Delta x}
</math>

where the mass, ''m'', per unit length of channel depends on the width and average depth of the river (''W'', ''d''), as well as the density of water (ρ). Thus, the equation becomes:

<math>
\frac{\Delta E_p}{\Delta x} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x}
</math>

where ''γ = ρg''. To turn ''ΔEp / Δ x'' into a loss per unit rate of time, we simply need to multiply both sides of the equation by ''Δx / Δt'', as follows:

<math>
\frac{\Delta E_p}{\Delta x}\cdot \frac{\Delta x}{\Delta t} = W \cdot d \cdot \gamma \cdot \frac{\Delta z}{\Delta x} \cdot \frac{\Delta x}{\Delta t}
</math>

On the left-hand side, the ''Δx'' terms cancel, leaving ''ΔEp / Δt'', which is referred to as the ''total stream power'', and represented using the symbol ''Ω'' (measured in W/m). On the right-hand side, we realize that the term ''Δz / Δx'' is simply the slope or gradient of the stream channel, which is represented by the symbol ''S''. The term ''Δx / Δt'' is the average velocity for the flowing water, represented by ''u''. Making these simplifications and substitutions, we get:

<math>
\Omega = W \cdot d \cdot u \cdot \gamma \cdot S
</math>

Referring back to the continuity equation above, it is clear that the geometric terms (''Wdu'') simply represent the discharge carried by the stream. Thus, ''total stream power'' can be most easily predicted knowing nothing more than the gradient of the river (''S'') and the discharge the river happens to be carrying (''Q''):

<math>
\Omega = \gamma \cdot Q \cdot S
</math>

===Unit stream power===
The stream power acting on a 1m by 1 m area of the channel bed is simply given by dividing ''Ω'' by the width of the stream. This is called the ''unit stream power'' (''ω'', in W/m2).

<math>
\omega = \frac{\Omega}{W}
</math>

<math>
\omega= \frac {\gamma \cdot Q \cdot S}{W}
</math>

<math>
\omega = \gamma \cdot q \cdot S
</math>

where ''q'' is called the ''specific discharge'' and can be estimated at any point from the local depth (''d'') and velocity (''u'').

<math>
q=d \cdot u
</math>

So, if we know only the width of a channel and the discharge, we can estimate only an average value for ''ω'' that represents the average rate of potential energy loss per unit area of the channel bed. But, if we have measured the local depths and velocities during, for example, a stream gauging exercise, we can estimate ''ω'' at any point across the channel where ''d'' and ''u'' have both been measured.