Post Mountain Pine Beetle Impact Assessment of Ts’il?os Provincial Park

 

 

March 2008

 

 

 

 

 

 

Submitted by

 

 

Amelia Needoba and Bruce Blackwell

B.A. Blackwell & Associates Ltd.

3087 Hoskins Rd.

North Vancouver, BC

 

 

 

Submitted to

 

 

 

Tom Hughes

Area Supervisor

South Chilcotin, Cariboo Region

Parks & Protected Areas

Ministry of Environment


Executive Summary

This project was initiated based on a coarse scale mountain pine beetle (MPB) impact assessment of 88 BC Parks in 2005 (B.A. Blackwell & Associates Ltd. 2005). Ts’il?os Provincial Park was ranked as having moderate probability of impact from attack but extreme consequences of impact from attack. This result primarily reflects the mixed species composition and diverse ecosystems within the Park (lower probability of impact due to degree of pine composition) and the many values for which the Park is managed (higher consequence of impact). The results from the multi-Park analysis, in combination with Park planning priorities, led to the initiation of this more detailed, Park level mountain pine beetle impact assessment.

 

The objectives of this impact assessment were to:

  • Identify impacts on Park values resulting from mountain pine beetle attack;
  • Assess impacts over time and space;
  • Provide a reasonable assessment of impact given limited data, time and resources;
  • Provide results that can be used for future planning; and,
  • Develop recommendations for mitigation strategies and budgeting at the Park level.

 

The methods employed to meet these objectives utilized existing MPB pine kill projection modelling data, various inventory data available for the Park and expert opinion to derive predictions of impact. The project was broken down into three phases:

 

  1. Analysis of primary impact;
  2. Analysis of secondary impact; and,
  3. Final impact overlay and recommended actions.

 

Primary impact was defined as the direct impact of mountain pine beetle caused tree mortality and was used to identify where pine mortality resulted in changes in openness or in pine composition. These two primary impacts were used to define the ‘area of primary impact’ within which secondary impacts could be assessed.

 

Secondary impact was defined as the impact on park values resulting from one of two activity alternatives in response to primary impacts: 1) do nothing; and, 2) remove hazard trees. A stepped impact matrix was used to scope Park values and potential impacts on those values for inclusion in the analysis. Significance thresholds were set for each predicted impact based on available literature and expert opinion. Data for each value was loaded into GIS in order to spatially define the area of each value and, ultimately, the area of impact for each value.

 

The significance of each of the values was considered within the Park boundary only as this is where BC Parks has jurisdiction. Therefore, this assessment does not reflect the significance of these values outside the Park where different management objectives apply. This assessment is coarse scale given the available data and its results are not field tested. In all cases, predicted impacts are simplified representations of complex effects. Predictions are limited due to available data, the current state of knowledge and the time frame of the assessment. This assessment is considered a suitable tool for broad strategic planning but is not appropriate for planning at the operational level. In addition, it is possible that values at risk not captured in this assessment may exist at the site level. Verification of actual or potential impacts by qualified professionals is required in the field prior to any activities being planned.

 

The final impact overlay provided a spatial definition of impact within the Park based on identified values judged to be impacted by ‘Do nothing’ or ‘Remove Hazard Trees’. The results of the spatial overlay defined where different actions were recommended based on the secondary impacts that occurred in each location. The area totals were as follows:

 

  • Hazard Tree Removal: 456 ha area of hazard tree removals reflected the significant predicted impact on public safety.
  • Mitigate:108 ha area of mitigation reflected where the impact on red-listed species, invasive plant species, visual quality, special features, fuel hazard and cultural features was predicted to be significant.
  • Monitor: 54,333 ha area of monitoring reflected where impact on structural stage, peak flows, terrain stability, surficial erosion and sedimentation, stream temperature and shading and riparian habitat was predicted to be significant.
  • Monitor Uncertainty: 50,728 ha area of monitoring uncertainty reflected the area of all values within the park for which significance could not be determined.
  • Total Park Area: 235,881 ha

 

Overall, results indicated that impacts are concentrated around Chilko Lake and the valleys radiating out from it (Long Valley, Franklyn Arm, Yohetta Creek, Rainbow Creek, Tchaikazan River and Edmond Creek). Those areas requiring tree removal and/or mitigation are focussed around the trails and public use areas. The large areas of “monitor uncertainty” and “monitor” reflect knowledge and data gaps that exist both in terms of data available for the Park and expected impacts of MPB on Park values.

 

The recommendations made based on the results of the analysis were:

 

1)      Increase the frequency of hazard tree assessments in areas where public safety impacts have been identified as significant or uncertain and concern has been confirmed on-site by a qualified professional.

 

2)      Where hazard tree removals are planned to occur consider the following:

a)      In areas with high pine mortality, consider removing a percentage of dead pine trees (even if not technically assessed as hazardous) to achieve an incremental reduction in dead pine volume. The percentage could be varied to maintain the aesthetic and environmental integrity of the area but should reduce the likelihood of a future spike in resource requirements for hazard tree removals.

b)     If peak flows or erosion/sediment delivery impacts are predicted to be significant and concern is confirmed on-site by a qualified professional, consider doubling riparian reserve zone widths and leaving slash on the ground unless it amounts to more than 5kg/m2 over an area of 100 m2.

c)      If terrain stability impacts are predicted to be significant and concern is confirmed on-site by a qualified professional, consider consulting a Professional Geologist to determine whether special removal methods are required.

d)     Design tree removals to ensure the protection of traditional use values. The current joint management process for the Park should be adequate to ensure the protection of traditional use values during tree removals.

e)      Modify tree removal methods/boundaries if invasive species, red-listed species or cultural features are suspected to occur on site following assessment by a qualified professional.

f)       If fuel hazard impacts are predicted to be significant and concern is confirmed on-site by a qualified professional, ensure that retained CWD is bucked and generally left lying flat to the ground (not elevated). Where slash exceeds 5kg/m2 over an area of 100 m2, consider piling and burning to reduce fuel load.

g)      Undertake tree removals in winter where possible to minimise potential hydrological, terrain stability and soil disturbance impacts.

h)     Design tree removals to minimize windthrow hazard by determining the susceptibility of residual tree species, historic windthrow patterns, topographic features and soil types and conditions. The MoFR procedure (Windthrow Field Card FS 712-2) could provide a suitable assessment method.

i)        If hazard tree removals in a riparian reserve zone are necessary, consult a fisheries biologist regarding potential mitigation strategies.

j)        If the riparian area consists of primarily dead pine (>60%) consult a fisheries biologist to determine whether rehabilitation of the riparian area is desirable.

k)     Retain species other than lodgepole pine unless they are identified as a hazard.

l)        Minimize damage to retained trees and understory vegetation.

m)   Leave slash greater than 15 cm in diameter on the ground unless it amounts to more than 5kg/m2 over an area of 100m2.

n)     Follow the Riparian Management Area Guidebook (FPC 1995) where tree removals occur near streams, rivers, lakes and wetlands.

o)      Avoid any mechanical or other disturbance in or within 20 m of S3 and S4 streams.

p)     Create tall stumps or stubs (3 to 5 m tall) as cavity sites where they do not create a worker safety or wildfire hazard.

q)     Retain CWD on site in a way that mimics its natural distribution of randomness and connectivity. Retain 2 to 5 pieces of greater than 12 cm dbh per square metre.

r)       Plan the timing of treatments to avoid disrupting applicable wildlife using the habitat. For example, avoid treatments in late winter/heavy snow years if ungulates are utilizing the site because the canopy is providing snow interception, or in spring/early summer during the nesting period.

s)       Stubs should be restricted to trees greater than 30 cm in dbh or where cavities already exist.

t)       When re-forestation is required, re-plant as soon as possible following treatment.

 

3)      Given the temporary nature of MPB related impacts, the frequency and intensity of any monitoring that takes place within the Park should be determined based on Park management priorities and no specific schedule or standard is recommended here. The values for which monitoring was recommended include peak flow, surficial erosion, stream temperature and shading, terrain stability, structural stage distribution and riparian habitat. In all cases where uncertainty was identified, the recommendation is monitoring to address uncertainty. The exception to that statement is areas of uncertainty regarding public safety; the current annual parks wildlife danger tree assessment process is considered adequate to address uncertainty in these areas and should be maintained. If assessments identify dead pine in these areas, then the public safety impact is considered ‘significant’ and recommendations 1 and 2 apply.

 

4)      Consider working towards filling knowledge gaps that inhibit the ability to plan with certainty in the Park. Terrestrial Ecosystem Mapping in the park would greatly improve the ability for planners to undertake habitat mapping and identify sensitive ecosystems. In addition, data on traditional use, cultural heritage values, invasive species and red-listed species should be compiled for the Park to improve planning.


Acknowledgements

The authors would like to thank the following individuals for their contributions to the project:

 

The Gwanajegwaghten (Ts’il?os advisory group), Tom Hughes, Rob Dolighan and members of the Ministry of Environment provided valuable information for the development of the values and significance thresholds for the analysis.

 

Rob Scagel implemented the projected pine kill model on the Park’s vegetation inventory data.

 

Avram Sandor, Tamaki Kano, Ben Andrew, Philipp Walter, Colin Mahony and Bob Green of B.A. Blackwell & Associates Ltd. provided assistance with GIS, report editing, peer review, graphics and the compilation of background data.


Table of Contents

1     Introduction. 1

1.1       Project Objectives. 1

1.2       Study area overview.. 2

2     Methodology. 4

2.1       Primary Impact 4

2.1.1     Change in Pine Composition. 5

2.1.2     Change in openness. 6

2.1.3     Area of Primary Impact 6

2.2       Secondary Impacts. 7

2.2.1     Scoping. 7

2.2.2     Significance. 11

Peak Flows. 12

Sediment Delivery/Surficial Erosion. 13

Stream Temperature and Shading. 14

Terrain Stability. 15

Old Forest 16

Structural Stage Distribution. 17

Riparian. 18

Mule Deer Habitat 19

Grizzly Bear Habitat 20

Black Bear Habitat 21

Bighorn Sheep Habitat 22

Moose Habitat 23

Mountain Goat Habitat 24

Fuel Hazard. 25

Public Safety. 26

Visual Quality. 27

2.3       Final Impact Overlay and Recommended Actions. 28

3     Results and Discussion. 29

3.1       Primary Impact 29

3.1.1     Change in Pine Composition and Openness. 30

3.1.2     Overall Primary Impact 31

3.2       Secondary Impact 35

3.2.1     Peak Flow.. 35

3.2.2     Surficial Erosion/ Sediment Delivery. 38

3.2.3     Stream Temperature and Shading. 40

3.2.4     Terrain Stability. 43

3.2.5     Old Forest 45

3.2.6     Structural Stage Distribution. 45

3.2.7     Riparian Habitat 48

3.2.8     Wildlife. 51

3.2.9     Fuel Hazard. 58

3.2.10       Public Safety. 61

3.2.11       Visual Quality. 64

3.3       Overlay of Secondary Impacts and Recommended Management Responses. 66

4     Assumptions, Limitations and Future Considerations. 68

5     Recommended Actions. 69

6     Budget estimates. 71

7     References. 72

 

List of Figures

Figure 1. Summary of methodology. 4

Figure 2. Average cumulative kill percentages for projected pine mortality within the Park. 30

Figure 3. Average % change in pine composition and % change in openness across the Park over time. 31

Figure 4. Representation of fire hazard succession following mountain pine beetle attack. 59

List of Tables

Table 1. Classification of pine composition. 6

Table 2. Classification of forest openness. 6

Table 3. Ts’il?os Provincial Park stepped impact matrix. 9

Table 4. Percent of Park area with watersheds potentially impacted by peak flow increases. 36

Table 5. Macro Stream Reaches with <20% forest cover over >10% of their length. 40

Table 6. Distribution of early-mid seral forest within the Park that is not modelled to meet Biodiversity Guidebook targets post-MPB attack. 45

Table 7. Riparian management areas where <60% of the RMA has >20% forest cover post-MPB attack. 48

Table 8. Summary of fuel hazard impact resulting from MPB attack in the Park. 58

Table 9. Estimated area where public safety is predicted to be impacted by MPB attack. 61

Table 10. Estimated length of named trails potentially impacted by MPB. 61

 


List of Maps

Map 1. Biogeoclimatic Ecosystem Classification within Ts’il?os Provincial Park. 3

Map 2. Change in pine composition from 1999 to 2024 (only considers stands with a pine component). 32

Map 3. Change in openness over time from 1999 to 2024 (only considers stands with a pine component). 33

Map 4. Overall primary impact area. 34

Map 5. Peak flows predicted to be potentially significant in watersheds due to an equivalent clear cut area of >35%. 37

Map 6. Areas where surficial erosion is predicted to be potentially significant due to a substantial reduction in forest cover on steep slopes. 39

Map 7. Fish bearing streams where stream temperature and shading is predicted to be potentially impacted by MPB attack based on substantial reductions in forest cover. 42

Map 8. Areas within the Park where terrain stability may potentially be impacted by substantial reductions in forest cover. 44

Map 9. Biogeoclimatic Ecosystem Classification zone and Natural Disturbance Type where distribution of early-mid seral forest is significant based on the Biodiversity Guidebook. 47

Map 10.  Sections of riparian habitat predicted to be potentially impacted by MPB attack. 50

Map 11. Areas where moderate-high potential mule deer habitat overlaps with primary MPB impact area. 52

Map 12. Areas where moderate-high potential grizzly bear habitat overlaps with primary MPB impact area. 53

Map 13. Areas where moderate-high potential black bear habitat overlaps with primary MPB impact area. 54

Map 14. Areas where moderate-high potential bighorn sheep habitat overlaps with primary MPB impact area. 55

Map 15. Areas where moderate-high potential moose habitat overlaps with primary MPB impact area. 56

Map 16. Areas where moderate-high potential mountain goat habitat overlaps with primary MPB impact area. 57

Map 17. Areas where fuel hazard is predicted to be potentially significant based on MPB impact. 60

Map 18 Areas where public safety is predicted to potentially be impacted by MPB attack. 63

Map 19. Areas where visual quality is predicted to be potentially impacted by MPB attack. 65

Map 20. Final impact overlay delineating areas for recommended actions. 67


1         Introduction

BC Parks has a mandate to manage for ecosystem processes. Park policy states:

Natural ecosystem processes affecting vegetation including fire, insects, disease, weather (i.e., wind, avalanches, etc.), herbivory by wildlife, and tree mortality due to age, are recognized as natural occurrences shaping vegetation. Ecosystems will be managed to maintain ecological processes in as natural a state as possible.

Numerous insects and diseases are active within Ts’il?os Provincial Park (hereinafter referred to as ‘the Park’). These forest health agents, when naturally occurring, contribute to the maintenance of biodiversity within the Park. However, when endemic pest populations reach outbreak or epidemic levels the scale of disturbance increases and cumulative effects can result in negative impacts on numerous forest and non-forest values.

 

The mountain pine beetle (Denroctonus ponderosae Hopk.) (MPB) is responsible for the largest recorded natural disturbance to occur both in British Columbia and its parks and protected areas system (Gawalko 2003). Mountain pine beetle is an endemic species in British Columbia (BC) but it is thought that human related activities have increased forest homogeneity in terms of age and species composition, and increased the density of vegetative cover. These changes have resulted in forest landscapes that are more susceptible to beetle attack (Samman and Logan 2000). If human activity has altered the landscape in a way that has facilitated the current MPB epidemic, then it may be appropriate to develop management strategies aimed at reducing its impacts (Samman and Logan 2000).

 

A coarse scale mountain pine beetle impact assessment was completed for 88 BC Parks in 2005 (B.A. Blackwell & Associates Ltd. 2005). The project quantified current and future impacts of the beetle on forested landscapes and related those impacts to the values that each Park was managed for. This was a coarse scale analysis that prioritized the 88 parks based on their assessed risk of mountain pine beetle impact. Ts’il?os Provincial Park was ranked as having high probability and extreme consequence of impact from attack. This result primarily reflects the high percentage of pine in the park and the many values for which the Park is managed. The results from the multi-Park analysis, in combination with Park planning priorities, led to the initiation of this more detailed, Park level mountain pine beetle impact assessment.

1.1      Project Objectives

The objectives of this impact assessment are to:

1)      Identify impacts on Park values resulting from mountain pine beetle attack;

2)      Assess impacts over time and space;

3)      Provide a reasonable assessment of impact given limited data, time and resources;

4)      Provide results that can be used for future planning; and,

5)      Develop recommendations for mitigation strategies and budgeting at the Park level.

1.2      Study area overview

Ts’il?os Provincial Park (235,881 ha) is a Class A Park and is recognized and managed for its wildlife, their habitat and wilderness recreation opportunities. The primary management goals are the protection of the ecosystems and conservation of biodiversity in the Park. The Park is also of great importance to the Xeni Gwet’in First Nation, which plays an important role in the management of their traditional home and the protection of culturally important sites. An extension to the Park termed ‘Goal 2’ is proposed for the north of the Park and is likely to be formalized in the near future. However, this area was not included in this analysis due to time and data constraints.

 

The Park encompasses mountains, lakes, rivers and glaciers and is distinguished by Chilko Lake, a high elevation lake, and the Tchaikazan Glacier which forms the headwaters of the Tchaikazan River. It has populations of large mammals including mountain goat, black bear, grizzly bear, moose, mule deer and California bighorn sheep.

 

Ts’il?os Park is located in the transitional area between the Coast Mountains and the Chilcotin Plateau and is an important part of the Central Chilcotin Ranges Ecosection. Five biogeoclimatic zones and 13 subzones are represented in the Park. The zones include: Coastal Western Hemlock, Interior Douglas-fir, Montane Spruce, Englemann Spruce-Subalpine Fir and Alpine Tundra (Map 1). Because of the diversity of climate and topography, the Park is biologically diverse with high species richness, functioning ecosystems, and good connectivity between habitat types as fragmentation by industrial use is low.

 

Based on Provincial vegetation resource inventory (VRI) data, approximately 29% (68,799 ha) of the park is vegetated, and of that area 61,700 ha is treed. The tree species that occur within the park, in approximate order of abundance, are Pinus contorta (lodgepole pine), Abies spp. (Balsam), Pinus albicaulis (whitebark pine), Pseudotsuga menziesii (Douglas-fir), Picea spp. (Spruce), Populus tremuloides (trembling aspen), Tsuga spp. (hemlock) and Populus balsamifera ssp. Trichocarpa (black cottonwood).

 

Map 1. Biogeoclimatic Ecosystem Classification within Ts’il?os Provincial Park.


2         Methodology

The project was broken down into three phases:

 

1.      Analysis of primary impact;

2.      Analysis of secondary impact; and,

3.      Final impact overlay and recommended actions.

 

Primary impact was defined as the direct impact of mountain pine beetle caused tree mortality.

 

Secondary impact was defined as the impact on park values resulting from one of two activity alternatives in response to primary impacts: 1) do nothing; and, 2) remove hazard trees.

 

Figure 1 summarizes the methodology of the impact assessment.

2.1      Primary Impact

The primary impacts considered were:

  • Change in pine composition; and
  • Change in openness.

 

Figure 1. Summary of methodology

Initially, change in forest composition from conifer dominated to mixed or deciduous dominated forest was also considered; however, a review of the results indicated that very few areas of the Park were assessed as being affected by this type of change and the outcome of the primary impact layer was not judged to be sensitive to the inclusion of this layer. Therefore, it was assumed that the primary area of impact due to MPB attack was adequately described by change in pine composition and change in openness.

 

The primary impact layers (change in pine composition and change in openness) were developed using lodgepole pine mortality projections from year 3 of the provincial-level projection of the current mountain pine beetle outbreak[1] (Eng et al. 2006) and overlaying these with forest attributes extracted from provincial Vegetation Resource Inventory (VRI) data[2]. The year 3 results were updated with the 2005 Provincial Aerial Overview of Forest Health[3]. Year 4 results updated with 2006 aerial overview data have been published but did not include the data projected into the future and the projected data could not be obtained in time to complete this analysis. However, any discrepancies between the 2007 aerial overview data and the projection data were reviewed and incorporated into the analysis as areas of uncertainty and there are only slight differences between the year 3 and year 4 results (Walton et al. 2007).

 

Cumulative beetle kill scenarios were analyzed for years 1999 to 2024. Each forest cover polygon was assigned an averaged cumulative kill percentage and classified according to the following forest attributes for each beetle kill scenario: 1) pine composition (stand pine composition); and, 2) openness (as measured by changes in canopy cover).

 

Three beetle kill scenarios were used to assign impact ratings to each of the primary impacts. These were: 1) Pre-attack stand (stand status pre-attack from VRI); 2) 2024 projected cumulative kill (projected percentage of pine removed from the stand); and, 3) 100% pine kill (assuming 100% pine removed from the stand). The change between pre-attack stands (1) and 2024 projected stands (2) was used to derive ratings of high, moderate and low primary impact. Uncertainty was accounted for by incorporating 100% pine kill (3) into the layer as an additional rating. For example, if there was no change between the pre-attack and 2024 model results but the 100% pine kill results indicated a change, then the polygon was assigned to ‘uncertainty’.

 

The purpose of the ‘uncertainty’ rating was to capture areas where primary impacts may occur but are not modelled to occur. In addition, model results for beetle kill were averaged based on their spatial location relative to VRI polygons, therefore, beetle kill may potentially be over- or under-estimated within each polygon. The incorporation of uncertainty using 100% pine kill provides an estimate of uncertainty where impact maybe under-estimated and therefore enables a conservative overall estimate of impact. Given the lack of data on the impact of MPB attack in BC Parks, a conservative estimate (i.e., one that incorporates the worst case scenario of 100% pine kill as uncertainty) was judged to be an appropriate approach given the values at risk.

 

Time was incorporated into the analysis by determining the year in which the greatest change was modelled to occur and assigning it to past (1999-2007), present (2008-2015) or future (2016-2024). Broad time categories were used to allow for some assessment of impact over time for planning purposes, while acknowledging that there is uncertainty in projections into the future and the timing of the change may not be accurate to the year.

2.1.1      Change in Pine Composition

The pine composition attribute was a measure of the residual pine composition following beetle attack. It captured the amount of pine that remained in the stand. The pine composition classification is summarized in Table 1.

 

Table 1. Classification of pine composition.

Class

Pine component of stand

No pine

< 10%

Minor pine

>/= 10% and < 20%

Mixed pine

>/= 20% and < 50%

Leading pine

>/= 50%

 

The ratings assigned based on the change from pre-attack pine composition to 2024 pine composition were as follows:

 

  • High: Leading pine class to minor or no pine class
  • Moderate: Leading pine class to mixed pine class / Mixed pine class to no pine class
  • Low: Mixed pine class to minor pine class / Minor pine class to no pine class

 

Uncertainty was assigned when 100% pine kill indicated a change that was different from the 2024 result.

2.1.2      Change in openness

The change in openness attribute was a measure of the change in canopy cover following beetle attack. The openness classification is summarized in Table 2.

 

Table 2. Classification of forest openness.

Class

Canopy Cover

Non-forested

0%

Open

< 20%

Partial

20% < 60%

Closed

>/= 60%

 

The ratings assigned based on the change from pre-attack pine composition to 2024 pine composition were as follows:

 

  • High: Closed class to Open or Non-forested class
  • Moderate: Closed class to Partial class/ Partial class to Open or Non-forested class
  • Low: Open class to Non-forested class

 

Uncertainty was assigned when 100% pine kill indicated a change that was different from the 2024 result.

2.1.3      Area of Primary Impact

The change in pine composition and change in openness ratings were combined spatially to derive an overall assessment of primary impact. Uncertainty was assigned when it was carried through from both primary impact values.

 

Numeric values were assigned to High (3), Moderate (2) and Low (1) ratings for both ‘change in openness’ and ‘change in pine composition’. The values were then summed to derive ‘primary impact’. Primary impact value was defined as follows:

 

  • Low = 1 or 2
  • Moderate = 3 or 4
  • High = 5 or 6
  • Uncertainty = uncertainty exists in both primary impacts


The creation of the overall primary impact layer was based on the assumption that portions of the Park that have had or will experience the greatest change in openness and/or pine composition due to pine mortality are representative of the areas in the Park most impacted by MPB attack. The area of primary impact defines the spatial area of impact and the uncertainty in the area of impact.

2.2      Secondary Impacts

2.2.1      Scoping

Scoping is defined as the process used to identify the issues and impacts that are likely to be important in an impact assessment. In order to determine which secondary impacts may occur, it was necessary to identify the potential activities that would take place in the Park to address primary impacts. Given the BC Parks mandate to maintain natural ecosystem processes, it was determined that two activity alternatives were appropriate:

 

  1. Do nothing (the default state); or,
  2. Remove hazard trees (where public safety is at risk).

 

Either one of these activities could then be assessed based on its impact on values at risk to determine whether impacts were ‘significant’ and therefore whether mitigative actions were required.

 

Values at risk were identified either as Valued Ecosystem Components (VECs) or Valued Social Components (VSCs) that were known or suspected of being sensitive to the primary impacts. VEC/VSCs were initially identified by the project team and then vetted through consultation with BC Parks. A stepped matrix approach was used to aid the determination of how VEC/VSCs could potentially be impacted by changes in openness and pine composition (Table 3).

 

An important assumption that underlies the assessment of secondary impacts is that any tree removal activities will not result in any new area of permanent access or conversion to non-forest (access and fragmentation have not been considered as impacts). Any tree removal that takes place will only be to remove dead pine and hazard trees, no other tree species will be modified.

 

The VEC/VSCs were selected based on Park management values and available data. Initially, an extensive list of potential impacts was prepared. Then, based on available data and expert opinion, the list was reduced down to a manageable number (Table 3) of values that either directly or indirectly represented potential impacts on the Park. The criteria used to select VECs/VSCs were as follows:

  1. The value reflects Park management objectives;
  2. The value is likely to provide a direct or indirect measure of impact by MPB attack or subsequent tree removals;
  3. Adequate data exists to define the value within the Park;
  4. Adequate data or knowledge exists to make a prediction on the magnitude of change in the value and the significance of the impact.

 

Values that did not meet these criteria were not included in the analysis but were considered during scoping are: soil disturbance, noise generation, specific fish species, invasive species, red-listed species, traditional use, cultural heritage and Park visitor use.

 

Traditional Use values and Cultural Heritage values are known to occur within the Park. However, no data could be obtained within the time frame of the analysis. Invasive alien plant species are also known to occur in the park but spatial data was inadequate to draw conclusions on potential distribution. It is also expected that red-listed plants and animals occur within the park but no spatial data exists to confirm their presence/absence.


Table 3. Ts’il?os Provincial Park stepped impact matrix; the X denotes whether or not an impact is predicted as a result of either action.

Values affected by MPB

Type of Action

Condition

 

Take no action

Remove dead/ hazard trees

Initial

Change

Predicted Impact

Valued Ecosystem Components:

 

 

 

 

 

Peak Flow

X

X

Reduced forest cover

Snow interception decreases

Peak flows are increased

Sediment delivery/surficial erosion

X

X

Reduced forest cover

Increased bare ground exposure to precipitation

Surficial erosion/sediment delivery is increased

Stream temperature and shading

X

X

Reduced forest cover

Increased ground and channel exposure to radiation

Stream temperatures are increased

Terrain Stability

X

X

Death/removal of trees

Roots decay

Terrain is less stable

Old Forest

X

X

Death/removal of old growth trees

Forest structure is altered

Less abundance of old seral stages

Structural stage distribution by BEC zone and NDT

X

X

Death/removal of trees

Structural stage distribution is altered

Increased abundance of early-mid seral stages

Riparian

X

X

Removal of trees

Disturbance of riparian habitat

Degraded riparian habitat

Mule Deer Habitat

X

X

Death/removal of trees in habitat

Forest structure and species composition is changed

Habitat is altered

 

Grizzly Bear Habitat

X

X

Death/removal of trees in habitat

Forest structure and species composition is changed

Habitat is altered

Black Bear Habitat

X

X

Death/removal of trees in habitat

Forest structure and species composition is changed

Habitat is altered

Bighorn Sheep Habitat

X

X

Death/removal of trees in habitat

Forest structure and species composition is changed

Habitat is altered

Moose Habitat

X

X

Death/removal of trees in habitat

Forest structure and species composition is changed

Habitat is altered

Mountain Goat Habitat

X

X

Death/removal of trees in habitat

Forest structure and species composition is changed

Habitat is altered

Values affected by MPB

Type of Action

Condition

 

Take no action

Remove dead/ hazard trees

Initial

Change

Predicted Impact

Valued Social Components:

 

 

 

 

 

Fuel hazard

X

X

Standing dead trees

Red attack/trees fall or are felled and left on site

Fuel hazard is increased

Public safety

X

 

Standing dead trees

Increased probability of tree failure

Trees fall

Visual quality

 

X

Removal of trees

Landscape is modified

Visual quality is reduced


2.2.2      Significance

Determining which changes were considered impacts and the magnitude of the impact required the prediction of significance. The indicators listed in Table 3 were used to predict significance and were based on available data, existing research and expert opinion. The following classifications were used:

 

  1. No Impact: If there is no overlap between the activity and the value in space or time the activity is assessed as having no impact.
  2. Significant Impact: A secondary impact that has the potential to negatively affect the value according to modelled primary impact results.
  3. Uncertain Impact: A secondary impact that is not predicted to potentially occur based on modelled primary impact results, but that could occur based on the 100% pine kill scenario.

 

The following tables outline how significance was determined for the impact on each value, the references used to determine the significance thresholds and the assumptions and limitations of the significance thresholds. Significance was then mapped for each impact within the Park.

 

In all cases, values that met significance thresholds pre-attack were not considered ‘significant’ post-attack because the impact was not significant due to MPB attack.


Peak Flows

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Peak Flows

x

x

Reduced forest cover

Snow interception decreases

Peak flows are increased

Equivalent Clear-cut Area > 35% of total watershed

GIS Layer Name: Peak_Flows

Data Sources used to Derive Layer:

Community watersheds (LRDW), WSA 3rd order and greater watersheds (LRDW), VRI data (LRDW), cumulative pine kill (Eng et al. 2006).

Significance Thresholds:

ECA > 35% of watershed in 2024 (Significant)

ECA > 35% of watershed assuming 100% pine kill (Uncertain)

Other:

ECA = Area of opening x (1 – hydrological recovery)

Divide ECA by area of entire watershed and multiply by 100 to achieve %

Assumptions and Limitations:

Removal of overstory canopy increases peak flow.

Stand volume reduced by < 30% = 100% hydrological recovery, Volume reduced by >/= 30% </= 60% = 50% hydrological recovery, Volume reduced > 60% = 0% hydrological recovery.

Elevation >800 m therefore ECA weighted x 1.

ECA is an appropriate calculation for beetle killed areas; residual forest structure in MPB killed stands may well be different from that of forest harvesting.

No other potential or recent disturbance within the watershed has been considered.

3rd order and greater watersheds are an appropriate scale at which to assess the ECA; each was assessed as an entire watershed or sub-basin, not broken down by sub-basins or, alternatively, aggregated into watersheds (e.g., 3rd order watersheds were not calculated as sub-basins of a 4th order watershed but as discrete units).

This assessment was not completed by a professional hydrologist and was not part of a formal Watershed Assessment Procedure.

References: Watershed Assessment Procedure Guidebook (http://www.for.gov.bc.ca/tasb/legsregs/fpc/FPCGUIDE/wap/WAPGdbk-Web.pdf)

 


Sediment Delivery/Surficial Erosion

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Increased risk of sediment delivery/surficial erosion

X

X

Reduced forest cover

Increased bare ground exposure to precipitation

Surficial erosion/sediment delivery is increased

Slope is >60% and canopy cover is reduced from >60% (closed) to < 20% (open)

GIS Layer Name: Surficial_Erosion_Sediment_Delivery

Data Sources used to Derive Layer:

TRIM derived slope raster (LRDW), WSA 3rd order and greater watersheds (LRDW), community watersheds (LRDW), VRI data (LRDW), cumulative pine kill (Eng et al. 2006).

Significance Thresholds:

Slope >60% and forest cover reduced from >60% (Closed) to < 20% (Open) by 2024 (Significant)

Slope >60% and forest cover reduced from >60% (Closed) to < 20% (Open) assuming 100% pine kill (Uncertain)

Assumptions and Limitations:

Areas with > 60% slope are likely to be at risk from erosion.

Areas with dense pine cover killed by mountain pine beetle are likely to have bare ground exposed following beetle kill until ground cover has had time to establish; this is likely to be highly variable in reality.

No advanced regeneration in stands.

Areas where heavy pine removal occurs will have bare ground exposed due to mechanical tree removal.

Lack of appropriate scale surficial geology data requires reliance on >60% slope as a crude estimate of erosion hazard; some soil types have high erosion hazard at lower slope gradients.

References: Community Watershed Guidebook (http://www.for.gov.bc.ca/tasb/legsregs/FPC/fpcguide/WATRSHED/watertoc.htm); Mapping and Assessing Terrain Stability Guidebook (http://www.for.gov.bc.ca/hfd/library/ffip/BCMoF1999_C.pdf)


Stream Temperature and Shading

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Stream temperature and shading

X

X

Reduced forest cover

Increased ground and channel exposure to radiation

Stream temperatures are increased

Canopy cover is <20% over >10% of the stream length

GIS Layer Name: Stream_Temperature_Shading

Data Sources used to Derive Layer:

Macro stream reaches (LRDW), FISS fish occurrences (LRDW), VRI data (LRDW), cumulative pine kill (Eng et al. 2006)

Significance Thresholds:

Canopy cover < 20% (open) over 10% of the stream length by 2024 (significant)

Canopy cover < 20% (open) over 10% of the stream length assuming 100% pine kill (uncertain)

Assumptions and Limitations:

Macro stream reaches with gradients of less than 12% could support fish (sections with gradients greater than that cannot).

Forest cover loss within the VRI polygon through which the stream passes is an appropriate assessment criterion for ‘significantly’ reduced forest cover and may be adequate to increase stream temperature.

Pine may have been identified in the VRI polygon (source of forest cover estimates) but may not be present within the riparian corridor due to the scale of mapping.

Increased stream temperature will not necessarily have a negative impact in all cases.

References: Riparian Management Area Guidebook (http://www.for.gov.bc.ca/tasb/legsregs/fpc/fpcguide/riparian/rip-toc.htm)

 


Terrain Stability

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Terrain Stability

X

X

Death/removal of trees

Roots decay

Terrain is less stable

Slopes >60% and change in openness from >20% to </= 20%

GIS Layer Name: Terrain_Stability

Data Sources used to Derive Layer:

TRIM derived slope raster, VRI data (LRDW), cumulative pine kill (Eng et al. 2006)

Significance Thresholds:

Slope >60% (significant) where change in openness is moderate or high by 2024 (significant)

Slope >60% (significant) where change in openness is moderate or high assuming 100% pine kill (uncertain)

Assumptions and Limitations:

Slopes >60% have a higher risk of landslide following timber harvesting or substantial reduction in forest cover.

No advanced regeneration in stands.

Inadequate data to produce terrain stability class mapping therefore crude measure of >60% is used to indicate terrain stability.

Will likely over-estimate areas where terrain stability impact is significant.

References: Mapping and Assessing Terrain Stability Guidebook (http://www.for.gov.bc.ca/hfd/library/ffip/BCMoF1999_C.pdf)


Old Forest

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Old Forest

X

X

Death/removal of old growth trees

Forest structure is altered

Less abundance of old seral stages

CWH/ESSF NDT2: Area of old forest <13%

MS NDT3: Area of old forest <21%

IDF NDT4: Area of old forest <19%

IMA NDT5: Area of mature and old forest <85%

GIS Layer Name: N/A

Data Sources used to Derive Layer:

VRI data (LRDW), BEC data (LRDW), cumulative pine kill (Eng et al. 2006)

Significance Thresholds:

When retention of total stand is =/>70% old forest structure is retained, when retention is < 70% old seral structure is lost

CWH/ESSF NDT2: Area of old forest <13%

MS NDT3: Area of old forest <21%

IDF NDT4: Area of old forest <19%

IMA NDT5: Area of mature and old forest <85%

Assumptions and Limitations:

Forests in old seral stage as defined by the Biodiversity Guidebook and the Field Manual for Describing Terrestrial Ecosystems are representative of old forest structure.

BEC is intersected with VRI to assign BEC and Natural Disturbance Type (NDT) to polygons: There may be inaccuracies in the line work leading to some polygons being assigned to inappropriate BEC and NDT.

In reality, old seral structure may be retained in some stands even though less than 70% of the stand survives.

All forest that could shift to old forest by 2024 was included in the post-attack analysis to account for old forest recruitment over time.

Other disturbances that could alter the distribution of old forest within the Park are not accounted for in this analysis.

References: Biodiversity Guidebook (http://www.for.gov.bc.ca/tasb/legsregs/fpc/fpcguide/biodiv/biotoc.htm), Field Manual for Describing Terrestrial Ecosystems (http://www.for.gov.bc.ca/hfd/pubs/Docs/Lmh/Lmh25.htm)

Structural Stage Distribution

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Structural stage distribution by BEC zone and NDT

X

X

Death/removal of trees

Structural stage distribution is altered

Increased abundance of early-mid seral stages

CWH/ESSF NDT2: Area of early and mid seral forest > 27%

MS NDT3: Area of early and mid seral forest > 35%

IDF NDT4: Area of early and mid seral forest > 23%

IMA NDT5: Area of early and mid seral forest > 15%

GIS Layer Name: Structural_Stage

Data Sources used to Derive Layer:

VRI data (LRDW), BEC data (LRDW), cumulative pine kill (Eng et al. 2006)

Significance Thresholds:

By 2024: When retention of total stand is <70% early and mid seral forest structures occur

CWH/ESSF NDT2: Area of early and mid seral forest > 27%

MS NDT3: Area of early and mid seral forest > 35%

IDF NDT4: Area of early and mid seral forest > 23%

IMA NDT5: Area of early and mid seral forest > 15%

Assumptions and Limitations:

Forests in early/mid stage as defined by the Biodiversity Guidebook and the Field Manual for Describing Terrestrial Ecosystems are representative of early/mid seral forest structure.

In reality, old and mature seral structure may be retained in some stands even though less than 70% of the stand survives.

BEC is intersected with VRI to assign BEC and Natural Disturbance Type (NDT) to polygons: There may be inaccuracies in the line work leading to some polygons being assigned to inappropriate BEC and NDT.

Other disturbances that could alter the distribution of early/mid seral stages within the Park are not accounted for in this analysis.

References: Biodiversity Guidebook (http://www.for.gov.bc.ca/tasb/legsregs/fpc/fpcguide/biodiv/biotoc.htm), Field Manual for Describing Terrestrial Ecosystems (http://www.for.gov.bc.ca/hfd/pubs/Docs/Lmh/Lmh25.htm)

 


Riparian

 

Type of Action

 

 

 

 

 

1. Take no action

2. Remove dead/ hazard trees

 

 

 

 

Values affected by MPB

 

Condition

Significance Threshold

Initial

Changes

Predicted Impact

 

Riparian

X

X

Removal of trees

Disturbance of riparian habitat

Degraded riparian habitat

<60% of the Riparian Management Area has > 20% canopy cover

GIS Layer Name: Riparian

Data Sources used to Derive Layer: Stream macro reaches (LRDW), TRIM lakes (LRDW), TRIM wetlands (LRDW), VRI data (LRDW)