TABLE OF CONTENTS

 

1. INTRODUCTION *

2. THE NAZKO RIVER WATERSHED *

2.1 LOCATION *

2.2 PHYSIOGRAPHY *

Figure 1 – The Nazko River IWAP Study Area *

2.3 CLIMATE *

Table 1. Mean maximum and minimum monthly temperature and total precipitation at Nazko, B.C. (52º55’N 123º37’W, 899m), 1974-1977. *

2.4 HYDROLOGY *

2.5 FOREST HARVESTING *

Table 2. Licensed water demand on the Nazko River (after Rood and Hamilton 1995) *

3.0 THE NAZKO WATERSHED ASSESSMENT PROCEDURE *

3.1 PURPOSE OF THE ANALYSIS *

3.2 METHODOLOGY 7

3.2.1 PROCEDURAL MODIFICATIONS 7

3.3 RESULTS OF THE ASSESSMENT *

Table 4. Areal measurements of the study area watersheds by elevation. *

Table 5. Peak Flow index measurements and results *

Table 6. Road inventory and density *

Table 7. Roads adjacent to streams *

Table 8. Riparian Buffer impacts *

Table 9. Mass Wasting hazard measurements *

Table 10. Other landuses *

Table 11. Watershed ownership characteristics *

Table 12. Physical characteristics by watershed *

Table 13. Watershed Report Card. *

3.4 INTERPRETATION OF THE RESULTS 19

3.4.1 ANEKO CREEK (APPENDIX D) *

Table 14. Aneko Creek Watershed Report Card. 20

3.4.2 BROWN CREEK (APPENDIX E) *

Table 15. Brown Creek Watershed Report Card. 20

3.4.3 CANYON CREEK (APPENDIX F) *

Table 16. Canyon Creek Watershed Report Card *

3.4.4 CLISBAKO RIVER (APPENDIX G) *

Table 17. Clisbako Creek Watershed Report Card *

3.4.5 GRUIDAE CREEK (APPENDIX H) *

Table 18. Gruidae Creek Watershed Report Card *

3.4.6 MICHELLE CREEK (APPENDIX I) 22

Table 19. Michelle Creek Watershed Report Card *

3.4.7 REDWATER CREEK (APPENDIX J) *

Table 20. Redwater Creek Watershed Report Card *

3.4.8 ROSS CREEK (APPENDIX K) *

TABLE OF CONTENTS (CONTINUED)

Table 21. Ross Creek Watershed Report Card *

3.4.9 SNAKING RIVER (APPENDIX L) *

Table 22. Snaking River Watershed Report Card *

3.4.10 SUMMIT CREEK (APPENDIX M) *

Table 23. Summit Creek Watershed Report Card *

3.4.11 TAUTRI CREEK (APPENDIX N) *

Table 24. Tautri Creek Watershed Report Card *

3.4.11 UDY CREEK (APPENDIX O) 24

Table 25. Udy Creek Watershed Report Card *

3.4.13 WENTWORTH CREEK (APPENDIX P) *

Table 26. Wentworth Creek watershed report card *

Table 27. Summary of Results 26

Table 28. Summary of Recommendations 27

REFERENCES *

 

 

 

LIST OF APPENDICES

APPENDIX A - Complete List of Nazko River Watershed Assessment Multi-agency Roundtable Members

APPENDIX B - Interior Watershed Assessment Procedure Blank Digital Forms

APPENDIX C - Forest Cover Data Statistics

APPENDIX D - Aneko Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX E - Brown Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX F - Canyon Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX G - Clisbako River Watershed IWAP Data Sheets and Map with Overlays

APPENDIX H - Gruidae Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX I - Michelle Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX J - Redwater Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX K - Ross Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX L - Snaking River Watershed IWAP Data Sheets and Map with Overlays

APPENDIX M - Summit Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX N - Tautri Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX O - Udy Creek Watershed IWAP Data Sheets and Map with Overlays

APPENDIX P - Wentworth Creek Watershed IWAP Data Sheets and Map with Overlays

 

 

1. INTRODUCTION

This document was commissioned by the Ministry of Environment, Lands and Parks, Cariboo Region. It contains the results of the application of the Level I Interior Watershed Assessment Procedure (IWAP) to 13 separate watersheds within the greater Nazko River Watershed. The procedure was undertaken as described in the September 1995 version of the Forest Practices Code of British Columbia’s Interior Watershed Assessment Guidebook. Any deviation from the procedure and methodology outlined in this guidebook is fully detailed in section 3.2.

Integral to the IWAP is the formation of a multi-agency roundtable whose purpose is to inform, direct and focus the scope of the analysis. Members to the roundtable should include representatives of resource agencies or land tenure holders with legal obligations in the area undergoing analysis. Members to the Nazko roundtable included representatives from Ministry of Forests; Ministry of Environment, Lands and Parks; timber licensees; and Department of Fisheries and Oceans. A complete list of roundtable members is given in Appendix A.

As described in the guidebook, the IWAP is designed to be applied to complete watersheds of an area between 5 and 500 km2 (500 to 50 000 ha). At 4150 km2, the Nazko River study area, as defined by the Ministry of Environment, is more than four times the upper limit of 500 km2. Carmanah Resources Ltd. was then directed by the multi-agency roundtable to choose representative watersheds in consultation with Pat Teti, the Regional Hydrologist for the Ministry of Forests in the Cariboo Region, and other roundtable members. Using fisheries values, timber harvesting history, and development plans as the decision criteria, the 13 watersheds were chosen for analysis.

The terms of reference for this assessment differs significantly from those presumed by the Forest Practices Code guidebook. That is, the IWAP was not initiated by the Ministry of Forests as part of the forest planning process. Instead, the Nazko River Watershed Assessment was funded as part of Forest Renewal British Columbia’s Watershed Restoration Program. The purpose of this analysis is to provide information for watershed restoration activities. This information may be used by landuse managers to make more informed decisions regarding further forestry development, fisheries and other land-based issues.

The balance of this report contains information directly relevant to the Nazko River Watershed Assessment. It is not intended to be, nor is it, a comprehensive hydrological analysis of the watershed. For additional information regarding the specific requirements and products of the IWAP, refer to the Forest Practices Code of British Columbia guidebook to the Interior Watershed Assessment Procedure (September 1995).

 

 

2. THE NAZKO RIVER WATERSHED

 

2.1 LOCATION

The Nazko River watershed (Figure 1) is located in the Interior Plateau of central British Columbia. The region is typified by low relief and a gentle, rolling landscape interspersed with generally round hills of 200 to 400 meters in relative elevation. The relief never varies more than 400 meters over the entire study area,. The drainage basin area is 4150 km2 (415 000 ha), and the Nazko River mainstem is 157 km in length.. The centre of the watershed is located at approximately 52º 45’N latitude and 123º 30’W longitude. In order, from north to south, the watersheds chosen for the application of the IWAP Level 1 are Redwater Creek, Udy Creek, Michelle Creek, Snaking River, Clisbako River, Canyon Creek, Wentworth Creek, Tautri Creek, Brown Creek, Aneko Creek, Ross Creek, Gruidae Creek, and Summit Lake.. Settlements within the watershed include the town of Nazko near the confluence of the Nazko River and Redwater Creek, and the Nazko Indian Reserve at the confluence of the mainstem and Michelle Creek.

 

2.2 PHYSIOGRAPHY

The Nazko River watershed lies within the physiographic regions of the Chilcotin and Nechako Pleauteaux. The topography is made up of gently rolling, undissected uplands interspersed with numerous small lakes and wetland areas. The bedrock geology is of volcanic origin, consisting of gently dipping, Miocene/Pliocene olivine basalt flows. The area was heavily glaciated during the Pleistocene. This is reflected by the often thick surficial deposits of mainly unconsolidated glacial deposits. Fluvial and glaciofluvial deposits dominate the mainstems of the Nazko and Clisbako rivers, while the remainder of the basin is overlain by basal and ablationa; moraine deposits of till (Lord and Walmsley, 1988).

The biogeoclimatic zones of the study area fall within the Subboreal Spruce and Engelmann Spruce-subalpine Fir zones. The Chilcotin pine subzone (SBSa) is predominant in the watershed, and is characterized by white spruce and lodgepole pine, and black spruce in wetland areas. The Douglas fir-white spruce subzone (SBSb) follows mainly along the corridor of the lower mainstem, and is dominated by douglas fir, white spruce, lodgepole pine, trembling aspen and paper birch. The wet subzone of the Engelmann spruce-subalpine fir zone (ESSFh) is characterized by subalpine fir and engelmann spruce, and occurs in localized regions of the eastern portion of the drainage basin (Annas and Coupé, 1979).

 

 

Figure 1 – The Nazko River IWAP Study Area

2.3 CLIMATE

The climate of the study area is characteristically dry. Average annual temperature is approximately 2° C, with average summer highs in July-August of 14° C and January winter lows averaging -13° C. Average annual precipitation is approximately 440 mm, with an estimated 40 percent falling as snow. The majority of rainfall occurs in June, July and August. Table 1 shows a summary of climatic data for the hamlet of Nazko, B.C.

 

Table 1. Mean maximum and minimum monthly temperature and total precipitation at Nazko, B.C. (52º55’N 123º37’W, 899m), 1974-1977.

                 
 

Precipitation (mm)

 

Temperature (ºC)

 

Total

StanDev

 

Max.

StanDev

Min.

StanDev

Mean

                 

January

28.66

12.86

 

-2.57

2.20

-16.13

1.89

-9.35

February

13.30

5.80

 

.10

6.00

-14.63

7.25

-7.27

March

22.80

6.00

 

4.00

1.01

-9.87

3.10

-2.93

April

9.97

4.36

 

12.03

1.85

-5.03

1.46

3.50

May

34.23

22.18

 

14.73

.87

-.40

.37

7.16

June

58.90

49.56

 

18.08

1.86

2.73

.67

10.4

July

71.48

26.68

 

20.08

1.94

5.58

1.13

13.19

August

62.08

44.50

 

20.58

2.79

5.30

1.31

12.94

September

31.63

21.15

 

18.48

2.53

1.08

1.00

9.78

October

25.87

16.42

 

10.27

2.06

-2.03

.65

4.11

November

36.83

35.07

 

1.73

1.76

-8.33

1.97

-3.30

December

35.60

6.61

 

-.33

.75

-12.50

2.05

-6.42

                 

Annual

431.35

251.19

           

 

 

2.4 HYDROLOGY

The Water Survey of Canada gauging station, Nazko River above Michelle Creek (No. 08KF001), has a flow record going back to 1965. Mean annual flow over the period of record is 4.61 m3s-1. To estimate the mean annual flow for the Nazko River at the mouth, a conversion of the measured discharge by a ratio of the total watershed area (415 000 ha) over the gauged watershed area (324 000 ha) gives a discharge of 5.9 m3s-1. Average maximum and minimum monthly discharges take place in May and January-February respectively. Summer and winter 7 day low flows calculated by Rood and Hamilton (1995) using records from 1981 to 1990 are 1.27 m3s-1 and 0.96 m3s-1 respectively. Low summer flows can affect fish habitat by increasing water temperatures, stranding juveniles and reducing rearing habitat, while low winter flows can dewater redds which harm incubating eggs by increasing the possibility of freezing (Rood and Hamilton, 1995).

 

2.5 FOREST HARVESTING

Of the 415 000 ha of the Nazko River watershed, 90% of the area is forested. Harvesting in the study area began in the mid-1960’s, the majority taking place over the past 15 years. Slocan Forest Product Ltd. is the largest licensee operating in the study area. From 1982 to 1997, areas logged and areas proposed for logging are estimated at 6.3 percent of the watershed; less than 2 percent of the area was logged before 1982. Within the Nazko watershed, both the Snaking (Rood and Hamilton 1995) and lower Clisbako rivers have been identified as areas of concern with regard to logging.

The major land use activities within the watershed are agriculture and forestry. Some mineral exploration also takes place in the upper Clisbako River sub-basin. Farming is predominant on the mainstem downstream of the Nazko reserve. Surface and groundwater extraction for irrigation, stock watering and domestic use impact the hydrological regime of the watershed mainly in the summer. Furthermore, stream sedimentation via bank erosion caused by cattle may have an effect on the stream. Table 2 summarizes licensed water demand for the Nazko River.

 

Table 2. Licensed water demand on the Nazko River (after Rood and Hamilton 1995)

   

Monthly Demand (Ls-1)

License Type

Total Demand

February

August

September

Domestic

73 m3d-1

0.84

.084

0.84

Irrigation

2,227,429 m3

0

208

86

Waterworks

0 m3d-1

     

Industrial

70 m3d-1

0.82

0.82

0.82

Conservation

0 m3 s-1

     

Mean Stream Flow (Ls-1)

1,270

2,700

2,090

 

 

 

3.0 THE NAZKO WATERSHED ASSESSMENT PROCEDURE

 

3.1 PURPOSE OF THE ANALYSIS

The purpose of the Level 1 Interior Watershed Assessment Procedure is twofold. Firstly it is designed to assess the sensitivity of a watershed to hydrological changes as a result of landuse change or disturbance. Secondly, it assists in quantifying the extent and impacts of current and historical landuse practices, and forestry in particular. Additionally, the Assessment is intended to meet the requirements of future assessments mandated by the Forest Practices Code

The procedure was carried out according to its hierarchical structure which assesses thirteen impact measurements and related anecdotal information. The data was collected from various sources. The two primary sources of relevant data were the Ministry of Forests and the Ministry of Environment. The Ministry of Environment provided the Terrain Resource Inventory Maps (TRIM) which formed the spatial and positional base for the analysis. These digital maps contained all of the landform, contour, hydrological and road information required for the analysis. They are the standard for spatial analysis in British Columbia, and are the most accurate mapping products presently available. The Ministry of Forests provided the digital Forest Cover maps which contained detailed information relating to all aspects of forestry activity. A blank data entry form is located in Appendix B. There are 4 chief areas of analysis: Peak Flow, Surface Erosion, Riparian Buffer, and Landslide measurement.

The purpose of the Peak Flow measurements are to evaluate the potential for changes in both the timing and magnitude of maximum discharge. In a snowmelt dominated hydrological regime, such as that of the Nazko study area, the peak flow occurs during the spring freshet as the winter snowpack melts. The H60 line, or line above which 60% of the watershed area lies, is an attempt at defining the location of the snowpack prior to the initiation of the spring freshet. Once this has been calculated, it is then possible to determine the areal extent of logging and categorize it as being either above or below the H60 line. Forest harvesting can change the shape of the natural hydrograph. When trees are removed water is no longer taken up by the vegetation thus increasing the soil moisture content. The primary effect of this removal is increased peak flows, as water is no longer required by the vegetation or to recharge soil moisture. The increased soil moisture content reduces the ability of the soil to absorb water. As a result, large storm events or rapid snowmelt more often results in the overland flow of water streamward. This process is additionally facilitated by the existence of roads which provide an additional sediment source. This process then delivers sediment to the stream. The water yield from logged areas is also increased as trees are no longer available to either intercept rain or snow, or shade the snow from direct solar radiation and wind during spring. So then, there is both more snow available for melting (ablation) and it melts faster. Taken together, the above affect the channel system by increasing sediment load and peak discharge.

Roads and especially road crossings are significant potential sources of sediment supply to the hydrological system. Roads expose soil to increased erosion by rain and snow, and can often provide pathways to streams, creeks and rivers. Road crossings are of particular concern because they are immediately adjacent to streams, shortening the distance from the sediment source. The Surface Erosion impact indicators are consequently measures of these variables.

The third impact category is Riparian Buffer. The riparian buffer refers to the vegetative region that borders a waterway or body. Riparian areas often contain smaller side channels, overhanging trees and shrubs, and deliver large woody debris to the channel. In short, they provide many of the characteristics necessary for productive fish habitat. Increased sedimentation or logging of the riparian buffer can destroy fish rearing and spawning sites, thus destroying a fishery value of the stream.

The increased soil moisture content sometimes associated with logging can provide additional impetus for increased occurrences of mass movement events. The addition of relatively modest amounts of water to an already saturated slope may result in slope failure and landslides. The last impact category, Landslides, provides a measure of the number of landslides as well as an inventory of steep and therefore potentially unstable slopes (>60%).

 

3.2 METHODOLOGY

3.2.1 PROCEDURAL MODIFICATIONS

The Forest Practices Code Interior Watershed Assessment Procedure Guidebook defines the methodology and procedures of the Interior Watershed Assessment. However, since its publication in 1995, data availability and the application parameters have changed. Most importantly for this assessment are the modification of the terms of reference (see section 1.0) and the delivery agency. The data from this analysis is not in any way part of the Ministry of Forests’ formal planning process. It is, rather, intended for inventory and watershed restoration purposes only. In addition, the delivery agency is not the Ministry of Forests as presumed by the guidebook. Instead, the report was commissioned by the Ministry of Environment, whose digital data delivery requirements differ significantly from those of the MoF. These two factors affected the prescribed procedure in the following manners.

The Ministry of Environment’s Terrain Resource Inventory Maps (TRIM) are now available for almost the entire province. With respect to positional accuracy and landform descriptions, these maps are the highest quality available today. Consequently, rather than changing the TRIM base from the more accurate North American Datum (NAD) 1983 to NAD 27 as suggested in the guidebook, the procedure was reversed. This method has the advantages of maintaining the positional accuracy inherent in the TRIM base, while at the same time producing a series of coverages in the MELP standard form.

A related modification was the preferential use of TRIM data whenever possible. The guidebook suggests that TRIM be used to provide a base with which to augment the Forest Cover data (FC1). In the interest of positional accuracy and appropriate formatting, the Forest Cover information was selectively chosen from the FC1 data to provide information pertinent to the analysis. The positional accuracy of the TRIM base was maintained and any further data compatibility problems were minimized.

 

3.3 RESULTS OF THE ASSESSMENT

This section presents the results of the application of the Interior Watershed Assessment Procedure to the Nazko River watershed. The tabular manner in which the results are presented is intended to allow the reader to easily compare the results of individual watersheds. Most useful are the density measurements, which allow for fair comparison of individual watersheds.

An extremely useful product of the IWAP is the mapping products which are produced as an integral part of the final product. These maps facilitate intuitive comparison of watersheds and visual analysis of the many phenomena which cumulatively result in the hazard impact scores. Small scale maps of the study area watersheds showing roads, rivers, wetlands, lakes, slopes greater than 60 percent, H60 lines and the tree stand heights used in the equivalent to clear-cut area (ECA) calculations are contained in Appendices D through P. With each watershed map is the completed IWAP calculation sheets. These sheets summarize the measurements, indicators, scores and hazard indices in each of the four hazard analysis areas (i.e. Peak Flow, Surface Erosion, Riparian Buffer and Landslides).

The watersheds under consideration vary greatly in relative areal extent. Table 4 gives the area breakdown for each watershed. Of note is the relative consistency of the elevations of the H60 lines, which range over 200 metres only. The purpose of the H60 line is to indicate the likely location of the snowline at spring freshet. It is designed to be applied over all Interior watersheds. Like all generalizations, it compromises specificity for general applicability. Consequently, there are exceptions. The Nazko study area watersheds are exceptions. They are low relief watersheds and the snow melts both earlier and more rapidly. The hydrograph peaks in early May and descends rapidly through June and July achieving baseflow in mid to late August (Rood and Hamilton 1995).

 

 

Table 4. Areal measurements of the study area watersheds by elevation.

Watershed Name

Watershed

Code

H60 elevation (metres)

Area below H60 line (ha)

Area above H60 line (ha)

Total watershed area (ha)

Aneko Creek

170-3243-725

1200

8242

12364

20606

Brown Creek

170-3243-704

1200

5311

7966

13277

Canyon Creek

170-3243-476

1220

2220

3331

5552

Clisbako River

170-3243-384

1240

32675

49013

81688

Gruidae Creek

170-3243

1120

2636

3953

6589

Michelle Creek

170-3243-246

1060

2046

3070

5116

Redwater Creek

170-3243-133

1020

2232

3348

5580

Ross Creek

170-3243-866

1100

4868

7303

12171

Snaking River

170-3243-087

1100

7553

11329

18883

Summit Creek

170-3243-900

1100

1896

2843

4739

Tautri Creek

170-3243-654

1180

26256

39383

65639

Udy Creek

170-3243-087-380

1100

15592

23388

38980

Wentworth Creek

170-3243-477

1220

3887

5831

9718

The Peak Flow index score results from the weighted combination of equivalent clear cut area (ECA). The ECA is calculated to estimate the hydrological recovery of a watershed following harvesting. As such, burn areas , slides and hydro-cuts are considered clear-cuts with 0 percent recovery. Private land of less than 15 percent of the total watershed area is excluded from these calculations. In the ECA calculations, tree height is used as a surrogate to estimate the extent of hydrologic recovery. The lower the tree height, the more hydrological equivalent the area is to a clear cut. The relative importance of the ECA is weighted according to whether it is above or below the H60 line. The ECA above the line is given 50 percent more weight in the summary calculation of the Peak Flow index score. Below, table 5 presents the ECAs for each watershed.

 

 

Table 5. Peak Flow index measurements and results

   

Below H60 line

Above H60 line

Peak

Watershed Name

Watershed

Code

Weighted ECA

(ha)

Weighted ECA

(ha)

flow index

score

Aneko Creek

170-3243-725

460

1392

0.21

Brown Creek

170-3243-704

1020

1965

0.50

Canyon Creek

170-3243-476

136

303

0.18

Clisbako River

170-3243-384

3491

3311

0.17

Gruidae Creek

170-3243

317

255

0.20

Michelle Creek

170-3243-246

254

630

0.50

Redwater Creek

170-3243-133

0

311

0.14

Ross Creek

170-3243-866

238

401

0.15

Snaking River

170-3243-087

794

2105

0.35

Summit Creek

170-3243-900

181

231

0.19

Tautri Creek

170-3243-654

2113

4745

0.25

Udy Creek

170-3243-087-380

2147

3905

0.34

Wentworth Creek

170-3243-477

306

1069

0.33

Road lengths for the watersheds were measured digitally directly from the relevant TRIM maps. In general, these maps do not show the full extent of roads within the watershed. It is often then preferred that the forest cover maps be used for calculations involving roads. However, the positional accuracy of these maps is significantly less than that of the TRIM maps. As this directly affects the hazard calculations, the TRIM data was preferentially used.

 

 

Table 6. Road inventory and density

   

 

 

Road above

H60 line

Road for

entire watershed

Watershed Name

Watershed

Code

Watershed area (ha)

Length

(m)

Density

(km/km2)

Length

(m)

Density

(km/km2)

Aneko Creek

170-3243-725

20606

46872

0.23

83463

0.41

Brown Creek

170-3243-704

13277

56022

0.42

100759

0.76

Canyon Creek

170-3243-476

5552

4450

0.08

5070

0.09

Clisbako River

170-3243-384

81688

113355

0.14

274663

0.34

Gruidae Creek

170-3243

6589

17316

0.26

30038

0.46

Michelle Creek

170-3243-246

5116

37602

0.74

58318

1.14

Redwater Creek

170-3243-133

5580

8643

0.15

19682

0.35

Ross Creek

170-3243-866

12171

24557

0.20

53458

0.44

Snaking River

170-3243-087

18883

66967

0.35

141973

0.75

Summit Creek

170-3243-900

4739

673

0.01

6226

0.13

Tautri Creek

170-3243-654

65639

220010

0.34

366651

0.56

Udy Creek

170-3243-087-380

38980

152882

0.39

296948

0.76

Wentworth Creek

170-3243-477

9718

36853

0.38

43603

0.45

Digital terrain maps are not available for the Nazko study area. It was consequently not possible to use these maps in the delineation of sensitive or erodible soils. For the purpose of the calculations, erodible soils were considered to be all slopes of greater than 60 percent. This is in adherence to the guidebook (pg. 59), but given the relatively flat nature of the watersheds it is not a satisfactory surrogate measurement. There are thick surficial deposits of glacial materials over the study area. Some of these materials are potentially erodible. The spatial union of erodible soils, roads and streams is shown in table 7. With respect to their genesis, the erodible soils density measurements are likely underestimates.

 

Table 7. Roads adjacent to streams

 

Water-

 

Water

Water

shed

Roads on erodible soils

Roads within 100m of stream

Roads within 100m of stream on erodible soils

Density of stream crossings

shed

shed

area

Length

Density

Length

Density

Length

Density

Number

Density

Name

Code

(ha)

(m)

(km/km2)

(m)

(km/km2)

(m)

(km/km2)

(no.)

(km/km2)

Aneko Creek

170-3243-725

20606

0

0

8640

0.04

0

0

12

0.06

Brown Creek

170-3243-704

13277

0

0

5800

0.04

0

0

11

0.08

Canyon Creek

170-3243-476

5552

0

0

146

0.00

0

0

0

0.00

Clisbako River

170-3243-384

81688

0

0

21162

0.03

0

0

39

0.05

Gruidae Creek

170-3243

6589

0

0

5555

0.08

0

0

8

0.12

Michelle Creek

170-3243-246

5116

0

0

3347

0.07

0

0

4

0.08

Redwater Creek

170-3243-133

5580

0

0

5600

0.10

0

0

4

0.07

Ross Creek

170-3243-866

12171

0

0

7065

0.06

0

0

8

0.07

Snaking River

170-3243-087

18883

0

0

13322

0.07

0

0

18

0.10

Summit Creek

170-3243-900

4739

0

0

229

0.00

0

0

1

0.02

Tautri Creek

170-3243-654

65639

0

0

29975

0.05

0

0

47

0.07

Udy Creek

170-3243-087-380

38980

0

0

13555

0.03

0

0

20

0.05

Wentworth Creek

170-3243-477

9718

0

0

1664

0.02

0

0

2

0.02

 

Table 8 compares the length of streams logged to total stream length. These numbers are used to calculate the Riparian Buffer hazard index.

 

Table 8. Riparian Buffer impacts

       

Portion of

Watershed

Watershed

Length of

Total

stream

Name

Code

stream

stream

logged

   

logged (m)

length (m)

(km/km)

         

Aneko Creek

170-3243-725

1084

129766

0.01

Brown Creek

170-3243-704

2082

42026

0.05

Canyon Creek

170-3243-476

151

36378

0.00

Clisbako River

170-3243-384

8781

435760

0.02

Gruidae Creek

170-3243

245

38452

0.01

Michelle Creek

170-3243-246

565

20515

0.03

Redwater Creek

170-3243-133

0

24432

0.00

Ross Creek

170-3243-866

64

54742

0.00

Snaking River

170-3243-087

3594

100163

0.04

Summit Creek

170-3243-900

0

16548

0.00

Tautri Creek

170-3243-654

10868

261162

0.04

Udy Creek

170-3243-087-380

2189

130047

0.02

Wentworth Creek

170-3243-477

522

27569

0.02

The full extent of fish bearing stream through the watersheds is not accurately known. For the purposes of this analysis, then, all non-intermitent streams were assumed to be fish bearing.

In the IWAP, unstable terrain is defined as slopes greater than 60 percent. Which, in the case of this analysis, is the same criterion used to determine the extent of erodible soils. Although the measurements are in some cases the same, they are attributed different impact hazard scores because they are reported in a different context. In other words, 1 km of road on a slope greater than 60 percent is given a different hazard potential than 1 km of road on erodible soil. Table 9 shows the length of roads on unstable terrain in the context of mass movement potential. It also gives the density of streams whose banks have been logged on slopes greater than 60 percent in km per square km. Again, the relatively flat nature of the study area is reflected in these measurements.

 

Table 9. Mass Wasting hazard measurements

     

Landslides

Roads on

Streams whose banks

Watershed

Watershed

Watershed

in

Unstable

have been logged on

Name

Code

area

Watershed

Terrain

slopes >60%

   

(ha)

Number

Density

Length

Density

Length

Density

       

(no./km2)

(m)

(km/km2)

(m)

(km/km2)

Aneko Creek

170-3243-725

20606

0

0.00

0

0.00

0

0.00

Brown Creek

170-3243-704

13277

0

0.00

0

0.00

0

0.00

Canyon Creek

170-3243-476

5552

0

0.00

0

0.00

0

0.00

Clisbako River

170-3243-384

81688

0

0.00

0

0.00

0

0.00

Gruidae Creek

170-3243

6589

0

0.00

0

0.00

0

0.00

Michelle Creek

170-3243-246

5116

0

0.00

0

0.00

0

0.00

Redwater Creek

170-3243-133

5580

0

0.00

0

0.00

0

0.00

Ross Creek

170-3243-866

12171

0

0.00

0

0.00

0

0.00

Snaking River

170-3243-087

18883

0

0.00

0

0.00

0

0.00

Summit Creek

170-3243-900

4739

0

0.00

0

0.00

0

0.00

Tautri Creek

170-3243-654

65639

1

0.00

0

0.00

0

0.00

 

Table 9 (CONTINUED). Mass Wasting hazard measurements

     

Landslides

Roads on

Streams whose banks

Watershed

Watershed

Watershed

in

Unstable

have been logged on

Name

Code

area

Watershed

Terrain

slopes >60%

   

(ha)

Number

Density

Length

Density

Length

Density

       

(no./km2)

(m)

(km/km2)

(m)

(km/km2)

Udy Creek

170-3243-087-380

38980

0

0.00

0

0.00

0

0.00

Wentworth Creek

170-3243-477

9718

0

0.00

0

0.00

0

0.00

Although forestry is the principle landuse activity in the study area, the area is also used for unrelated activities. Among these is agriculture. The area, particular in the northern part of the study area and around the village of Nazko, is used for both crops and cattle. The area is also used for game hunting, and many of the roads are open to recreational use. Some other relevant landuses are found detailed below.

 

Table 10. Other landuses

Watershed

Watershed

Range use close

Mining close

All-terrain vehicles

Name

Code

to streams?

to streams?

close to streams?

   

(yes/no)

(yes/no)

(yes/no)

Aneko Creek

170-3243-725

no

no

no

Brown Creek

170-3243-704

no

no

no

Canyon Creek

170-3243-476

no

no

no

Clisbako River

170-3243-384

no

no

no

Gruidae Creek

170-3243

no

no

no

Michelle Creek

170-3243-246

yes*

no

no

Redwater Creek

170-3243-133

no

no

no

Ross Creek

170-3243-866

no

no

no

Snaking River

170-3243-087

no

no

no

Summit Creek

170-3243-900

no

no

no

Table 10 (CONTINUED) Other landuses

Watershed

Watershed

Range use close

Mining close

All-terrain vehicles

Name

Code

to streams?

to streams?

close to streams?

   

(yes/no)

(yes/no)

(yes/no)

Tautri Creek

170-3243-654

no

no

no

Udy Creek

170-3243-087-380

yes*

no

no

Wentworth Creek

170-3243-477

no

no

no

 

The Nazko study area is almost entirely crown land, most of which is used for timber supply. Table 11 summarizes these ownership properties.

 

 

Table 11. Watershed ownership characteristics

Water-

Water

Water-

Crown

Private

Indeterminate Crown/Private

shed

Shed

shed

land

land

land

Name

Code

area (ha)

area (ha)

%

area (ha)

%

area (ha)

%

Aneko Creek

170-3243-725

20606

20152

97.8

454

2.2

0

0.0

Brown Creek

170-3243-704

13277

13180

99.3

97

0.7

0

0.0

Canyon Creek

170-3243-476

5552

5550

100.0

0

0.0

0

0.0

Clisbako River

170-3243-384

81688

81140

99.3

548

0.7

0

0.0

Gruidae Creek

170-3243

6589

6407

97.2

182

2.8

0

0.0

Michelle Creek

170-3243-246

5116

4875

95.3

21

0.4

0

0.0

Redwater Creek

170-3243-133

5580

5368

96.2

171

3.1

0

0.0

Ross Creek

170-3243-866

12171

11936

98.1

224

1.8

0

0.0

Snaking River

170-3243-087

18883

18397

97.4

457

2.4

0

0.0

Summit Creek

170-3243-900

4739

3190

67.3

0

0.0

0

0.0

Tautri Creek

170-3243-654

65639

55492

84.5

418

0.6

0

0.0

Udy Creek

170-3243-087-380

38980

37565

96.4

914

2.3

0

0.0

Wentworth Creek

170-3243-477

9717

9715

100.0

0

0.0

0

0.0

Table 12 provides physical descriptions of the watersheds, including hydrological zone, bedrock geology and area with unstable slopes. As described earlier, unstable areas are defined as those areas with slopes greater than 60 percent. In the Nazko study area, there is a small percentage of area which meets this criteria due to the relatively flat to rolling nature of the landscape. Appendices D through P contain individual overlays which illustrate the extent of these areas within each watershed. Dominant bedrock geology is numerically coded to a table contained in the IWAP guidebook. Class "9" is described as "volcanic (basalt) - flows/breccias/porphyries/greenstone" (page 79). As such, it best describes the bedrock geology of the Nazko study area. However, the area was glaciated repeatedly during the Pleistocene (at least three times). In some areas, then, there exist deep deposits of surficial materials of glacial origin.

Table 12. Physical characteristics by watershed

     

Area with

Do DFO

Hydro-

Dominant

Are there

Water-

Watershed

Watershed

unstable

and/or

logical

bedrock

any

shed

Code

area

slope

MELP

zone?

geology

glaciers?

name

 

(ha)

(%)

have any

   

(yes/no)

       

concerns?

     

Aneko Creek

170-3243-725

20606

0

no

Central Interior Plateau

9

no

Brown Creek

170-3243-704

13277

0

no

Central Interior Plateau

9

no

Canyon Creek

170-3243-476

5552

0

no

Central Interior Plateau

9

no

Clisbako River

170-3243-384

81688

0

no

Central Interior Plateau

9

no

Gruidae Creek

170-3243

6589

0

no

Central Interior Plateau

9

no

Michelle Creek

170-3243-246

5116

0

no

Central Interior Plateau

9

no

Redwater Creek

170-3243-133

5580

0

no

Central Interior Plateau

9

no

Ross Creek

170-3243-866

12171

0

no

Central Interior Plateau

9

no

Snaking River

170-3243-087

18883

0

no

Central Interior Plateau

9

no

Summit Creek

170-3243-900

4739

0

no

Central Interior Plateau

9

no

Tautri Creek

170-3243-654

65639

0

no

Central Interior Plateau

9

no

Udy Creek

170-3243-087-380

38980

0

no

Central Interior Plateau

9

no

Wentworth Creek

170-3243-477

9718

0

no

Central Interior Plateau

9

no

The watershed report card, table 13, summarizes the hazard index scores in each of the impact categories for each watershed. The index scores were calculated from the data derived from the digital database which was described in detail in sections 3.2.1, 3.2.2 and 3.2.3. The calculations were completed using a combination of the IWAP Excel data entry and calculations spreadsheets, and long-hand checks. The measurements and calculations undertaken as part of this assessment followed closely the procedure described in the guidebook. The hazard impact index numbers displayed below provide only coarse indications regarding watershed hydrological impacts. As a coarse filter, the numbers simply point to potential concerns. It is the interpretation of the impact number which contextualizes the concerns.

 

Table 13. Watershed Report Card.

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Aneko Creek

170-3243-725

0.21

0.12

0.03

0.00

Brown Creek

170-3243-704

0.50

0.18

0.17

0.00

Canyon Creek

170-3243-476

0.18

0.02

0.01

0.00

Clisbako River

170-3243-384

0.17

0.09

0.07

0.00

Gruidae Creek

170-3243

0.20

0.18

0.02

0.00

Michelle Creek

170-3243-246

0.50

0.27

0.09

0.00

Redwater Creek

170-3243-133

0.14

0.18

0.00

0.00

Ross Creek

170-3243-866

0.15

0.15

0.00

0.00

Snaking River

170-3243-087

0.35

0.21

0.12

0.00

Summit Creek

170-3243-900

0.19

0.03

0.00

0.00

Tautri Creek

170-3243-654

0.25

0.15

0.14

0.01

Udy Creek

170-3243-087-380

0.34

0.17

0.06

0.00

Wentworth Creek

170-3243-477

0.33

0.10

0.06

0.00

 

 

3.4 INTERPRETATION OF THE RESULTS

Alone, the hazard impact scores provide a gross indication of potential hydrological impacts. The purpose of this section is to interpret the numbers within the context of the data quality and limitations described.

The following section will review and the interpret the specific assessment results on a watershed by watershed basis. The individual watershed data entry and calculations sheets, together with the watershed maps, should be used to augment the understanding of the results and aid interpretation. The sheets and maps are to be found in Appendices D through P. All of the information referenced in the following watershed sections comes directly from either the sheets or the maps

 

 

3.4.1 ANEKO CREEK (APPENDIX D)

The hazard indices for Aneko Creek indicate no cause for concern. The scores are all low, and there are no obvious data quality problems.

Table 14. Aneko Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Aneko Creek.

170-3243-725

0.21

0.12

0.03

0.00

 

3.4.2 BROWN CREEK (APPENDIX E)

The Peak Flow hazard index for Brown Creek is moderately high. All other scores indicate no cause for concern. The Peak Flow index is the result, in this case, of high ECA, particularly above the H60 line. Because of the rather gentle topography in the watershed, the hazard index may be overestimating the potential dangers of peak flow. However, there are a large number of streams in the watershed, many of which intersect with forest activity, which are marked as intermittent on the TRIM sheets and were consequently removed from the analysis. Many of these features appear non-intermittent on aerial photographs and on NTS maps. Consequently, the Surface Erosion and Riparian Buffer indices may be underestimating the extent of damage, and the potential for impacts in the watershed.

Table 15. Brown Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Brown Creek

170-3243-704

0.50

0.18

0.17

0.00

 

3.4.3 CANYON CREEK (APPENDIX F)

The hazard indices for Canyon Creek indicate no areas of concern. The scores are all low, and there are no obvious data quality problems.

 

 

Table 16. Canyon Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Canyon Creek

170-3243-476

0.18

0.02

0.14

0.00

 

3.4.4 CLISBAKO RIVER (APPENDIX G)

The hazard scores for the Clisbako watershed all indicate no significant impacts or potential impacts. One possible concern, however, is the concentration of logging activity along the mainstem of the Clisbako River. Because the watershed is quite large, the unlogged portions of the watershed act as a statistical buffer, and the density measurements are consequently decreased. Although the overall portion of the watershed which has been logged is quite small, the portion of the mainstem itself which has been logged is cause for potential concern.

 

Table 17. Clisbako Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Clisbako Creek

170-3249-384

0.17

0.09

0.07

0.00

 

3.4.5 GRUIDAE CREEK (APPENDIX H)

The hazard indices in the Gruidae Creek watershed are all low, and no areas of concern are indicated.

 

Table 18. Gruidae Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Gruidae Creek

170-3243

0.20

0.18

0.02

0.00

 

 

 

3.4.6 MICHELLE CREEK (APPENDIX I)

The Peak Flow index for Michelle Creek is moderately high, indicating potential concern. This score is driven largely by the road density above the H60 line, suggesting that it may be overestimating the degree of impact. However, the total road density (1.14 km of road per square km) is sufficiently high as to be of concern itself. A second, lesser, concern, which also contributes to the high Peak Flow score, is the high ECA.

Table 19. Michelle Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Michelle Creek

170-3243-246

0.50

0.27

0.09

0.00

 

3.4.7 REDWATER CREEK (APPENDIX J )

The hazard indices all indicate no cause for concern.

 

Table 20. Redwater Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Redwater Creek

170-3243-133

0.14

0.18

0.00

0.00

 

3.4.8 ROSS CREEK (APPENDIX K)

 

The hazard indices show no areas of concern.

 

Table 21. Ross Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Ross Creek

170-3243-866

0.15

0.15

0.00

0.00

 

 

 

3.4.9 SNAKING RIVER (APPENDIX L)

 

The hazard scores indicate no obvious impacts in the Snaking River watershed. One possible area of concern, though, is the high density of logging activity along the mainstem of the Snaking River. Although the actual density of logging in the watershed is quite low, there may be impacts on the mainstem itself.

Table 22. Snaking River watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Snaking River

170-3243-087

0.35

0.21

0.12

0.00

 

 

3.4.10 SUMMIT CREEK (APPENDIX M)

There are no obvious areas of concern in the Summit Creek watershed. The hazard scores are all low, and there are no apparent data quality problems.

 

Table 23. Summit Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Summit Creek

170-3243-900

0.19

0.03

0.00

0.00

 

 

3.4.11 TAUTRI CREEK (APPENDIX N)

There are no areas of concern indicated by the hazard indices. One potential concern not indicated, though, is the density of stream crossings and density of roads within 100 m of streams. If intermittent streams are included in the analysis — and for the reasons stated in section 3.4.2 there is reason to believe in some cases they ought to be — the number and density of these indicators would increase dramatically. This would consequently increase the Surface Erosion score.

 

Table 24. Tautri Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Tautri Creek

170-3243-654

0.25

0.15

0.14

0.01

 

 

3.4.12 UDY CREEK (APPENDIX O)

 

Although the hazard scores do not indicate any immediate cause for concern in the Udy Creek watershed, two other factors could indicate impacts on the hydrological system. First, the number of intermittent streams indicated in the watershed may be erroneous, and their inclusion in the assessment could significantly affect the results. Second, there is evidence of range use next to the creek which is not accounted for in the Forest Cover data. The possible consequences of range use are increased damage to the riparian buffer, increased sedimentation of the stream, and possible small scale slope failures along the banks.

Table 25. Udy Creek watershed report card

 

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Udy Creek

170-3243-087-380

0.34

0.17

0.06

0.00

 

 

3.4.13 WENTWORTH CREEK (APPENDIX P)

 

There are no evident causes for concern in the Wentworth Creek watershed.

 

Table 26. Wentworth Creek watershed report card

Watershed Name

Watershed

Code

Peak Flow

Surface Erosion

Riparian Buffer

Mass Wasting

Wentworth Creek

170-3243-477

0.33

0.10

0.06

0.00

 

4.0 CONCLUSIONS

4.1 SUMMARY OF RESULTS AND RECOMMENDATIONS

 

The Interior Watershed Assessment Procedure is a general assessment tool designed for application throughout the interior of British Columbia. As a general diagnostic filter, it measures many important hydrological and landuse parametres. The hazard index calculations used to determine the scores are formulated to be used in an average interior watershed. When considered from this point of view, the Nazko watershed is arguably not an average interior watershed. Its geological and geomorphological history conspire to make the IWAP too general to result in a completely accurate assessment. And this should be taken into consideration when interpreting the results and considering future assessments. The flat to rolling topography of the study area renders the H60 line largely insignificant: the low relief means that rather than having a snowpack which lasts further into the spring, it undergoes snowmelt earlier and more rapidly. That is, there is not a significant higher elevation snowpack which remains later into the spring to be affected by clearcutting. Instead, it melts earlier due to the relatively low elevations. Thus the differential ECA weighing for areas above and below the H60 line does not accurately estimate the actual situation.

Another anomaly of the study area is the glacial history. The area was glaciated numerous times. Most recently, it was glaciated during the late Wisconsinan Fraser glaciation. Pre-existing thick deposits of surficial materials and relatively soft volcanic bedrock combined to create a heavily incised landscape which has a drainage density which exceeds that necessary to balance hydrological inputs. This is illustrated by the large number of wetlands within the study area. In the context of the IWAP procedure, the effect is to remove many riparian habitats from the assessment. Consequently, in some cases the damage to the riparian buffer may be underestimated. Table 27 summarizes the findings of the assessment.

 

Table 27. Summary of Results

Watershed Name

Watershed

Code

Primary

Concern

Secondary

Concern

Aneko Creek

170-3243-725

None

None

Brown Creek

170-3243-704

Peak Flow due to logging activity

Possible damage to Riparian Buffer due to possible underestimation of stream length

Canyon Creek

170-3243-476

None

None

Clisbako River

170-3243-384

Possible damage to Riparian Buffer along mainstem

None

Gruidae Creek

170-3243

None

None

Michelle Creek

170-3243-246

Road Density

Peak Flow due to logging activity

Redwater Creek

170-3243-133

None

None

Ross Creek

170-3243-866

None

None

Snaking River

170-3243-087

Possible damage to Riparian Buffer due to underestimation of stream length

None

Summit Creek

170-3243-900

None

None

Tautri Creek

170-3243-654

None

Possible surface erosion due to underestimation of stream length

Udy Creek

170-3243-087-380

Possible damage resulting from range use

None

Wentworth Creek

170-3243-477

None

None

 

Two of the watersheds assessed in this analysis proved to contain potential significant hydrological concerns. As such they are eligible for additional Level 2 (Channel Assessment Procedure). Furthermore, concerns over data quality prompt the recommendation of aerial photograph based CAPs for three additional watersheds. The application of a Sediment Source Survey to two of these watersheds is also recommended.

As a watershed assessment and inventory project, this analysis has shown the forest harvesting and associated activities, such as road building, have potentially affected the hydrological regimes of five of the thirteen watersheds included in this study.

The application of the Channel Assessment Procedure (CAP) to selected watersheds within the study area would further define the location and magnitude of the hydrological impacts. Following this, restoration prescriptions may be appropriate. The decision to apply the CAP should be dependent on both the value of the stream and its prospects for restoration. Table 28 summarizes the recommendations following the application of the Interior Watershed Assessment Procedure to the Nazko River study area.

 

Table 28. Summary of Recommendations

Watershed Name

Watershed

Code

Initial

Recommendations

Aneko Creek

170-3243-725

None

Brown Creek

170-3243-704

Channel Assessment Procedure; Sediment Source Survey

Canyon Creek

170-3243-476

None

Clisbako River

170-3243-384

Aerial Photograph based Channel Assessment Procedure

Gruidae Creek

170-3243

None

Michelle Creek

170-3243-246

Channel Assessment Procedure; Sediment Source Survey

Redwater Creek

170-3243-133

None

Ross Creek

170-3243-866

None

Snaking River

170-3243-087

Aerial Photograph based Channel Assessment Procedure

Summit Creek

170-3243-900

None

Tautri Creek

170-3243-654

None

Udy Creek

170-3243-087-380

Aerial Photograph based Channel Assessment Procedure

Wentworth Creek

170-3243-477

None

 

 

REFERENCES

Annas, R.M. and R. Coupé (eds.) 1979. Biogeoclimatic zone and subzones of the Cariboo Forest Region. Ministry of Forests, Victoria, B.C. 103pp.

Carmanah Research Ltd. 1997. Nazko River Stream and Fish Habitat Assessment. Prepared for the Ministry of Environment, Lands and Parks, Cariboo Region, Williams Lake, B.C.

Lord, T.M. and Walmsley. 1988. Soils of the Nazko Area, British Columbia. Land Resource Institute, Vancouver, B.C. 62pp, 2 maps.

Struik, L.C. 1988. Structural Geology of the Cariboo Gold Mining District, East-central British Columbia. Geological Survey of Canada, Memior 421. 100 pp, 4 maps.

Forest Practices Code of British Columbia. 1995. Interior Watershed Assessment Procedure Guidebook (IWAP). Province of British Columbia. 82pp.

Geological Survey of Canada. 1959. Map no. 12-1959, sheet 93B, bedrock geology, Quesnel, B.C. Department of Energy, Mines and Resources, Ottawa, Ont. (map only)

Geological Survey of Canada. 1960. Map no. 49-1960, sheet 93G, bedrock geology, Quesnel, B.C. Department of Energy, Mines and Resources, Ottawa, Ont. (map only)

Rood, K.M. and R.E. Hamilton. 1995. Hydrology and water use for salmon streams in the West Road Habitat Management Area. British Columbia. Can. Manuscr. Rep. Fish. Aquat. Sci. 2295: 82pp.

Tipper, H.W. 1971a. Multiple Glaciations in Central British Columbia. Canadian Journal of Earth Sciences, 8: 743-752.

Tipper, H.W. 1971b. Glacial Geomorphology and Pleistocene History of Central British Columbia. Geological Survey of Canada, Bulletin 196.