Biodiversity Indicators Show Climate Change Will Alter Vegetation in Parks and Protected Areas

While multifaceted, a chief aim when designating parks and protected areas is to support the preservation of biological diversity, in part, through representing and conserving the full range of landscape conditions observed throughout a representative area. Parks and protected areas are, however, typically developed using a static interpretation of current biodiversity and landscape conditions. The observed and potential climate change impacts to biodiversity have created a need to also contemplate how parks and protected areas will respond to climate change and how these areas will represent the future range of landscape conditions. To assess change in biodiversity, broad-scale ecosystem information can be sourced from indirect remotely sensed indicators. Quantifying biodiversity through indirect indicators allows characterization of inter-relationships between climate and biodiversity. Such characterizations support the assessment of possible implications of climatic change, as the indicators can be generated using modeled forecasts of future climatic conditions. In this paper we model and map impacts of climate change on British Columbia’s parks and protected areas by quantifying change in a number of remotely sensed indicators of biodiversity. These indicators are based on the measured amount of incoming solar energy used by vegetation and map the overall annual energy utilization, variability (seasonality), and latent or baseline energy. We compare current conditions represented by parks and protected areas, to those forecasted in the year 2065. Our results indicate that parks and protected areas are forecasted to become more productive and less seasonal, due to increased vegetation productivity in higher elevation environments. While increased vegetation productivity may be beneficial for biodiversity overall, these changes will be particularly problematic for sensitive and specialist species. Future gaps in vegetation conditions protected by parks and protected areas are observed in the eastern edge of the Rocky Mountains and the central interior region of British Columbia. Protected areas along the Coast Mountains, Vancouver Island highlands, and the Rocky Mountains show the greatest levels of change in the biodiversity indicators, including decreasing seasonality, with the Mountain Hemlock ecozone most at risk. Examples of large parks that are predicted to experience rapid change in vegetation characteristics include Strathcona, Garabaldi, and Kitlope. Our maps of future spatial distributions of indirect biodiversity indicators fill a gap in information products available for adaptive parks management and provide an opportunity for dialogue and further research on the use of future scenarios of landscape conditions in conservation planning.

Abstract: While multifaceted, a chief aim when designating parks and protected areas is to support the preservation of biological diversity, in part, through representing and conserving the full range of landscape conditions observed throughout a representative area.Parks and protected areas are, however, typically developed using a static interpretation of current biodiversity and landscape conditions.The observed and potential climate change impacts to biodiversity have created a need to also contemplate how parks and protected areas will respond to climate change and how these areas will represent the future range of landscape conditions.To assess change in biodiversity, broad-scale ecosystem information can be sourced from indirect remotely sensed indicators.Quantifying biodiversity through indirect indicators allows characterization of inter-relationships between climate and biodiversity.Such characterizations support the assessment of possible implications of climatic change, as the indicators can be generated using modeled forecasts of future climatic conditions.In this paper we model and map impacts of climate change on British Columbia's parks and protected areas by quantifying change in a number of remotely sensed indicators of biodiversity.These indicators are based on the measured amount of incoming solar energy used by vegetation and map the overall annual energy OPEN ACCESS

Introduction
Parks and protected area systems play a key role in conserving biological diversity.Many international and national conservation targets call for between 10-12% protection of every ecosystem or habitat type [1] with a new target of 17% recently set by Aichi Biodiversity Target [2].The conservation targets have been criticized for their over-simplified and arbitrary values placed on biological elements yet applauded for their intention in achieving a more resilient and representative global protected areas system.As one might expect, opinion in areas to be protected often differ between those making recommendations based upon ecological principles and those responsible for land management and the often monetary implications.Regardless, in order for protection to be effective, all habitat types must have representation, resiliency, and redundancy within a protected area system [3], and while the allocation of land for conservation is often opportunistic [4] conservation networks should be designed to consider the spatial distribution of biodiversity, ecosystem threats [5], habitat configuration, environmental representation, and species at risk [6].
An assumption of protected area planning is that designated areas offer long term environmental conservation.However, planning does not typically account for anticipated landscape and ecological change that may occur over time [7][8][9].Climate change is understood to threaten global ecosystem functions and has been shown to potentially cause rapid shifts to species' geographic range, diversity, and abundance [10][11][12].As a result, protected areas will rapidly diminish in their intended purpose of preserving natural conditions and representing a particular suite of landscape conditions, indicating a need for a more adaptive planning process and prioritization of appropriate policy and/or management responses [13].Previous research has predicted that 37-48% of Canada's protected areas could experience a conversion in ecological classification type due to climate change [8] requiring information to support the prioritization, revision, or expansion of protected areas most expected to be impacted by climate change is required to inform proactive planning and management [14][15][16].
Climate change has been identified as a threat to the preservation of biological diversity where it is altering environmental conditions at a rate at which species cannot effectively adapt in situ or migrate [17,18] resulting in different potential impacts to geographic ranges [19].Species and ecosystem responses to climate change will thus be diverse and complex.Given climate is the main driving and limiting force determining productivity, measures of productivity offer insights into the amount of species richness and abundance a terrestrial habitat can sustain [12,[20][21][22].
Flexible, broad scale, monitoring opportunities such as afforded by remote sensing provide quantitative information for mapping ecosystem productivity dynamics over large areas that are difficult to measure and assess on the ground [23,24].Satellite based measures of vegetation productivity for example, have been utilized extensively as indirect indicators of biodiversity [21,[25][26][27] and show promise in detecting conditions of sensitive ecosystems and parks ecosystem representation [28].To date, satellite measures of productivity have been correlated with plant canopy cover [29,30] and overall productivity [23] in historical [31][32][33] and current contexts [34,35].Relationships between productivity indicators and biodiversity have also been used to assess future spatial distributions of species [36] and biodiversity [26,37].Time series of remotely sensed productivity components have been linked with climate data and applied to forecast future indicators of biodiversity [37][38][39][40][41].For example, [26] predicted future productivity in Washington and Oregon by linking climate data to metrics of fractional Photosynthetically Active Radiation (fPAR) (a measure for the amount of energy absorbed by the plant canopy).Other studies have used remotely sensed productivity indicators to predict climate alterations to ecosystem dynamics such as [42] and [43].The findings from these studies showed that climate change is often predicted to cause warmer and wetter conditions that are likely to increase forest productivity and, hence, biological diversity [25].
Our goal in this paper is to examine potential alterations to vegetation conditions in protected areas due to climate change for the over 94 million hectare province of British Columbia, Canada.To meet our research goal, forecasted changes in biodiversity metrics are quantified to (1) characterize how the frequency of biodiversity indicators in parks and protected areas will behave under a changing climate in order identify possible future gaps in the types of vegetation conditions protected by parks; (2) evaluate potential changes to biodiversity indicators by biogeoclimatic zones in order to identify broad species assemblages and ecosystem types that are most threatened by climate change; and (3) rank protected areas by mean level of potential change as a mechanism for identifying parks that may experience the most change.To complete this research we initially modeled and mapped future indirect indicators of biodiversity based on a historical archive of remotely sensed and climate data and future climate change scenarios [44].We compared biodiversity conditions represented within the parks and protected areas system at present and expected in the future.By forecasting changes in vegetation conditions province-wide and within protected areas, we can inform strategies for long-term preservation of biodiversity.

Study Area
British Columbia is located in the Pacific Northwest of North America from 49°00'N to 60°00' N and spans 944,735 km 2 .The province is highly diverse with temperate rainforests, deserts, and tundra.British Columbia also has highly variable topography and a wide-range of climates from mild coastal maritime to highly seasonal northern-interior continental conditions [45].Urban and agriculture centers are concentrated in the south juxtaposed with biologically rich landscapes [46].The rural areas of British Columbia are relatively untouched with most disturbances traced to forest management activities, ongoing infestation by mountain pine beetle (Dendroctonus ponderosae), regionalized drought conditions, and wildfires [47,48].Protected areas as classified by the International Union for Conservation and Nature (IUCN) as categories I-VI, include ecological reserves, conservancies, provincial and national parks.Using this definition British Columbia has 13.7% land area protected [49] within 1008 parks.

The Dynamic Habitat Index
The Dynamic Habitat Index (DHI) is an example of an indirect indicator of biodiversity derived from remotely sensed imagery and is based on annual trends in the amount of photosynthetically active energy absorbed by the plant canopy (400 nm-700 nm) (fPAR) [48].Shown to represent annual ecosystem vegetation dynamics, fPAR has been associated with biodiversity including butterfly and avian species richness [50,51].Using high temporal resolution satellite sensors, such as AVHRR and MODIS (AVHRR: Advanced Very High Resolution Radiometer; MODIS: MODerate-resolution Imaging Spectrometer.),fPAR can potentially be mapped twice daily and can be used to map intraannual productivity conditions [52].
The DHI is a combination of three annual fPAR metrics: cumulative greenness (greenness), coefficient of variation (seasonality), and minimum cover.High greenness values indicate landscapes of year-round productivity with stable and temperate climates.Higher greenness equates to more productive habitat [53,54] and greater species diversity [55,56].Therefore, protected areas with possible increases in greenness may see an increase in both species abundance and/or richness and existing species may become more dynamic.Examples of high greenness habitats include productive coniferous forests at lower elevations.Reduction in greenness, as observed in some dry interior habitats, may indicate that the environment will become less productive [50,57].For instance, bunchgrass may become more prolific as forested (higher productivity) environments are being extirpated from the region.
High seasonality values indicate habitats that are often snow covered for at least some portion of the year (interior or northern locations), or that occur at higher elevations, with vegetation growth mostly occurring in the warmer summer months.Examples of high seasonality environments include high elevation alpine habitats and deciduous forests with significant winter snowfall.Decreasing seasonality indicates less variation in productivity between winter and summer growing seasons.Decreased seasonality is a key factor for species succession and movement, especially at higher elevation environments because it can limit nutrient cycles and biomass accumulation and therefore restrict species movement up elevation gradients [58].
High minimum cover is indicative of habitats with increased productivity, often densely vegetated and with limited snow cover [48].An increase to minimum cover can result in a greater forecasted baseline year-round productivity levels and therefore equates to more year round sustainable food supplies and biomass [50].High minimum cover environments are often associated with riparian areas and productive coniferous forests along the coast and central interior regions of British Columbia.
In previous work, a historical archive of DHI [31] derived from AVHRR data was compared to historical climate data [44], and models developed to generate potential future spatial distributions of each DHI metric [44].Historical productivity metrics were generated using archival remote sensing information as processed by [31].Annual (1 km × 1 km resolution) DHI data from 1987-2007 were linked, using recursive partitioning (Rpart) regression trees [59], to 10 corresponding climate and landscape variables that consisted of mean annual temperature, mean annual precipitation, temperature difference, growing degree days, number of frost free days, mean coldest month temperature, Julian date on which frost free period begins, climate moisture deficit, evapotranspiration, and mean summer precipitation.
Observed climate-DHI relationships were then projected to 2065 using established climate change scenarios.To do so, three climate scenarios from the Canadian Centre for Climate Modelling and Analysis (CCCma) derived from Climate WNA (Western North America) version 4.6.[60] were used.The future climate models were projected to a thirty year average 2050-2080, hereafter referred to as 2065, and ranged from B1 (AR4-R1), representing the least extreme scenario, A1 (AR4-R1), representing the business-as-usual scenario where atmospheric CO 2 to continues to increase unabated, and A2 (AR4-R4), representing the most extreme scenario.To simplify the comparison between present day and 2065 results, the climate model A1 (CCCma A1-AR4-R1) was chosen because it was positioned in the middle of the projected precipitation and temperature ranges for mainstream climate models [61].Monserud et al. [62] utilized three general circulation model climate outputs (CGCM2, UK Hadley Centre and the Max Plank Institute) and found all three predicted similar trends of changing climate over the next 100 years across Alberta.Mote et al. [63] compared the ability of 10 global climate models to track recently recorded trends in temperature and precipitation in the Pacific Northwest and reported that the Canadian model was consistent in predicting a relatively rapid rate in warming and increases in precipitation when compared to actual observations.

Protected Area and Ecological Classification Data
Parks data were provided by the government of British Columbia data server last updated in July 2012 [64].Zones from British Columbia's province-wide Biogeoclimatic Ecosystem Classification (BEC) system (Version 8) developed by the BC Ministry of Forests were used to provide more species specific information about protected areas and what species assemblages might be most impacted.The BEC system categorizes ecosystems hierarchically according to vegetation, soil, topography, and climate classifications to create homogenous units of similar physical and ecosystem characteristics [65].At the most broad scale, 16 zones provide broad spatial species assemblage information for the entire province [65].

Methods
In order to categorize potential change in the DHI characteristics of parks through time, DHI values were classified into nine categories using natural breaks [44].The changing area (or number of pixels) in each DHI class, between the present and future (modelled year 2065) DHI maps, provided an indication on how well a current park's characteristics may be represented by parks in the future.
Gaps in protection are observed when DHI characteristics are currently well represented by protected areas, however not represented in the future.Gaps indicate vegetation conditions that may become minimally or no longer protected by the parks system as climate change impacts ecological processes.Several methodological steps were used to identify gaps.First, the current and future DHI values for protected areas were extracted from each DHI element (cumulative greenness, seasonality, and minimum cover).Second, each of the nine categories per DHI component (cumulative greenness, seasonality, or minimum cover) was mapped by park for current and future conditions.Lastly, changes in each DHI component, by category, were assessed by scenario to identify the locations of gaps.Categories less than 1000 km 2 in area were excluded from the analysis.
At the ecozone level to quantify spatial variation in predicted broad-scale habitat change, image differencing was used to measure regional change in indicators by BEC zones from present day to 2065.Greenness, seasonality, and minimum cover metrics were each subtracted by pixel to produce a forecasted change map.Change was summarized by ranking the mean change of BEC zones.
Park size was identified as a potential factor in ranking change, as smaller parks had lower sample sizes and thus could be entirely located within regions of high or low change.To account for the influence of park size, parks were divided into three categories for ranking, large (>1000 km 2 ), medium-sized (50-1000 km 2 ), and small (<50 km 2 ) based on natural breaks.Parks were ranked by size and change in DHI categories to highlight the most threatened parks and characterize the nature of modelled change.

Results
Presently, all DHI classes are unequally represented across the total province and within parks (Figure 1 a,b,c).Forecasted change indicates that in the future the DHI value categories will potentially be more evenly represented, both generally and within parks.Currently, parks and protected areas overrepressent low greenness (Figure 1a), high seasonality (Figure 1b) and low minimum cover (Figure 1c).By 2065 there is potential for increased greenness and reduced seasonality within parks.Gaps in cumulative greenness protection occurred in the mid-low values to mid-high values and were located along the eastern and central regions of the province.Gaps in seasonality, only, were more geographically limited compared to greenness only gaps which occur in moderate greenness locations from the Rocky Mountains to the central interior.Minimum cover only gaps were observed in mid-low values stretching the latitudinal length of the province with clusters situated in the eastern edge of the Coastal Mountains and valleys of the Rocky Mountains.Figure 2 shows 81 parks that are identified as having gaps in all three DHI variables, and have a total area of 5,105 km 2 .
The 10 top ranking large, medium-sized, and small parks, according to mean level of forecasted change to their greenness, seasonality, and minimum cover, are shown in Figures 3-5.Impacts for individual parks were spatially variable with some parks having consistently high degrees of change for all three indirect biodiversity indicators.Protected areas along the Coast Mountains were consistently ranked among those with the greatest levels of change, followed by the Rocky Mountains and the northern central interior.Protected areas with the least amount of forecasted change were located on low lying coastal zones and leeward interior parks.Protected Areas with all forecasted to have decreases in seasonality and increases in greenness and minimum cover.Of these, Strathcona Provincial Park is the most sensitive to change with mean greenness increasing 173.1% greater than the provincial average, seasonality 242.7% above average levels, and minimum cover increasing 498.3% above average levels.Other large parks were impacted by only one indicator, such as the Stitkine River Provincial Park, which was forecasted to have a large decrease in seasonality while attaining low to average levels of change in greenness and minimum cover.
Medium-sized parks (50-1000 km 2 ) varied greatly in their rankings between DHI categories; however, many are located along the western mountainous region of British Columbia.Parks consistently predicted to experience high levels of change included Tetrahedron Provincial Park, Fiordland Conservancy, and Ecstall-Spoksuut Conservancy, all located along the Coast Mountains.For example, Tetrahedron Provincial Park by comparison to mean medium-sized parks was forecasted to have greenness increase 230.4% above average, seasonality increase 301.2% above average, and minimum cover increase 405.4% above average.In the southcentral interior, some parks had negative forecasted values.For example, Kalamalka Lake Protected Area, which was forecasted to have a decrease in greenness of 158.5% below the mean for medium-sized parks, is indicative of the complex interrelationships observed in the drier interior environments.
Like medium-sized parks, small parks (<50 km 2 ) with considerable projected change were located in the Coast Mountains of British Columbia.The rankings for high changing small area parks were often not consistent across DHI categories resulting in a wide variety of protected areas identified as high risk.Examples include Mount Elliot Ecological Reserve, Tsitika Mountain Ecological Reserve, and Williams Creek Ecological Reserve.The Mount Elliot Ecological Reserve (3.3 km 2 ) on northern Vancouver Island had a predicted decrease in the seasonality value of 734.8% above average, greenness value increase of 869.3% above average, and minimum cover value increase 490.3% above small park average.Small parks are often located strategically to protect a particular set of values.For instance, Milligan Hills Provincial Park is located in the northeast Boreal region and protects woodland caribou habitat.Milligan Hills had a projected decrease in greenness metric of 47.5% below the provincial mean and an increase in minimum cover of 141.2% above average indicating possible alterations to vegetation and the spatial distribution of forging resources.
At the BEC zone level, projected changes in ecosystem biodiversity were represented by ranked mean change in DHI component (Table 1).The Mountain Hemlock zone had a 332.0%increase in greenness greater than the provincial average, followed by Coastal Western Hemlock ecozone (213.2%above average).Interior Douglas-fir zone was observed to have a decrease in mean greenness of 141.7%, Ponderosa Pine zone decreased 110.4%, and Bunchgrass zone decreased 148.2%; however, they did not change considerably in seasonality or minimum cover.BEC zones with the least amount of change were observed in geographically restricted zones or areas of limited area such as the Coastal Douglas-fir, Bunchgrass, and Ponderosa Pine zones.

Discussion
Our research provides maps of future, potential spatial distributions of indirect indicators of biodiversity, which are informative and quantitative input to adaptive management.Both the Canadian Biodiversity Strategy [16] and British Columbia's Climate Change Adaptation Strategy [66] state that improved data and monitoring are essential to aid adaptive management and ensure ecological resilience.Although climate change is perceived to be impacting Canadian parks, there is limited information on what may occur in the future [67].Mapping potential change to biodiversity supports landscape-scale strategic conservation frameworks that are designed to aid adaptive management [67].
Provincially, our model showed an increase in greenness and minimum cover and a decrease in seasonality.Generally, provincial trends were mirrored by the parks system.The current system favours higher elevation environments that have lower productivity (greenness and minimum cover) and higher seasonality [4], and due to high rates of change forecasted to occur by 2065 in mid-to high-elevation landscapes, increased greenness and reduced seasonality are predicted.The forecasted changes to DHI varied and, as a result, there was also variability to how well the parks represent the range of landscape conditions.In some respects, the findings appear positive.Given that parks presently over-represent low greenness values, while high greenness values are under-represented, the forecasted shift has potential to create an increasingly balanced representation of greenness landscapes within protected areas.Greater greenness may mean more biomass and enhanced productivity [50] and lower seasonality habitats can sustain greater biodiversity [58].It is worth noting that this representation is province wide, with the regionalized considerations where large changes in conditions are expected remaining important (i.e., localized changes, and possible loss of unusual conditions, niches, or particular species).Some species will not benefit as well (in relation to the general trends) based on predicted scenarios of future DHI distributions.Consider, for instance, the northern woodland caribou (Rangifer tarandus-caribou) that inhabit the large Omineca Provincial Park (1359.5 km 2 ), which was forecast to have a high rate of change with large increases in greenness, minimum cover, and decreases in seasonality.Omineca Park is an important calving location in the spring for the blue-listed caribou species.For food, shelter, and reduction of predation risk, caribou favour alpine and mid-elevation environments in Omineca Park [68] which are among the highest changing types of environments.As a result, increased vegetation productivity in the alpine will have negative impacts on the caribou in Omineca.Likewise small southern parks that host the most endangered species [46] are less able to adapt to changes in their habitats and a greening effect may have a negative impact on their survivability [25].Other examples of negative impacts include the landscape becoming more homogenous, with loss of heterogeneity, potentially resulting in a reduction in overall biodiversity.Further, rapidly shifting environmental conditions (i.e., DHI) have many unknown consequences since individual species' resilience to change and degradation thresholds are difficult to assess [69].Not all DHI values are better represented in the future.For instance, reduced representation occurred in locations with mid-to low-level minimum greenness values, which are found throughout most of the interior (i.e., Tweedsmuir and Wells Grey Provincial Parks) and in areas associated with lower productivity forests.
Meeting biodiversity related mandates requires representative protection of the landscape [67].By identifying under-represented biological characteristics in the existing protected areas system, we identified gaps in protection [70], which in turn are predicted to be most at risk to climate change.The Rocky Mountains and central interior were both identified as locations where gaps in all three DHI characteristics were evident.Gaps had DHI characteristics of high greenness and mid-level minimum cover and seasonality, and were located in the following BEC zones: Interior Cedar Hemlock, Sub-Boreal Spruce, and Montane Spruce.Other gaps were observed at Wells Grey Provincial Park and appear to be located in the valleys dominated by Interior Cedar Hemlock zones up in elevation to the Engelmann Spruce-Subalpine Fir zones.The Wells Gray Provincial Park Master Plan states the need to preserve park values of wildlife and ecosystem health [71] and our results indicate further assessment of climate change impacts to wildlife and ecosystem health will be beneficial.
In contrast, General regions such as the Boreal (in the northeast) and Okanagan (in the south region) demonstrated decreasing predicted greenness.Although the interior landscape did not rank amongst the greatest changing region, it did present complex combinations of negative greenness, decreased seasonality, and increased minimum cover which indicates dry forests associated with this unique DHI characteristic may shift northward.The southern interior region has a number of parks that are projected to observe a decrease in greenness.The South Okanagan Grasslands Protected Area, for example, has a projected decrease in greenness.Decreased greenness in the grasslands region equates to an expansion of dry grassland habitats and less favourable conditions for the megafauna species like the dry forests that have a higher greenness DHI characteristic.The Grasslands Protected Area hosts 15 endangered species such as the red listed Lyall's mariposa lily (Calochortus lyallii) and the flat-topped broomrape (Orobanche corymbosa ssp.Mutabilis) [72].Many of the endangered plants that thrive in grasslands may benefit from expansion of their habitat due to diminished climate limiters that once restricted migration into new regions.Changes to climate may also impact species negatively in the south, such as the bird species that require dry forests for their habitat.Decreased greenness indicates that the dry forests within the Grasslands Protected Area may decline causing red-listed species such as the Brewer's sparrow (Spizella breweri breweri) and the sage thrasher (Oreoscoptes montanus) to reduce in population.To protect these bird species in the future, parks must either expand into territories where habitats are projected to be more favourable or park conservation efforts should implement strategies to preserve species abundance and diversity within park boundaries.
At the ecosystem level, our models indicate the Mountain Hemlock BEC zone will see the greatest change by 2065.Mountain Hemlock BEC zone (407-1428 m elevation) is located in the transition zone of Coastal Western Hemlock and higher alpine habitats.Mountainous temperate regions have been found to be highly sensitive landscapes to climate change worldwide [13].High alpine forests and BEC zones such as Mountain Hemlock and Englemann Spruce-Subalpine Fir are sensitive to slight changes in climate since their distribution is predominately determined by climate and are pressed between highly competitive ecosystems [73,74].Protected areas closer to the coast are forecasted to have little change in biodiversity indicators; however, the Coastal Western Hemlock BEC zone borders on the highest changing ecosystems (Mountain Hemlock) and has a high mean rate of change.Low elevation parks in the Coastal Douglas-fir BEC zone were forecasted to have minimal change with only increases in minimum cover, while inland Coastal Hemlock BEC zone parks changed considerably.Interior BEC zones (i.e., Bunchgrass, Interior Douglas-fir, and Ponderosa Pine) were found to expand due to their forecasted decreases in greenness.Decreased greenness signifies that drier ecosystems will dominate over wetter forests of the interior such as those in Montane Spruce and Sub-Boreal Spruce BEC zones.Grasses and drier forests will spread throughout much of the southern interior of the province.Our findings are corroborated by BEC zone forecasting completed by Wang et al. [75] where warmer and wetter future conditions will lead to species groupings to shift northward and to higher elevations.[75].
Small protected areas (>50 km 2 ) are commonly designed to protect specific species or unique habitats and often protect the most threatened regions [6].The most extreme range in forecasted DHI results were found for small parks, they are often surrounded by regions of strong change.Williams Creek Ecological Reserve (7km 2 ) in the northwest of the province represents one of the highest forecasted increases in greenness because it is located along the Coast Mountains, is situated at an elevation gradient of 490-1,735m, and includes fast changing ecosystem zones such as Mountain Hemlock and Coastal Mountain-heather Alpine [76].Increased greenness may result in a push of dominant ecosystem species up elevation gradients, allowing coastal hemlock to overtake traditional Mountain Hemlock zones and Mountain Hemlock trees can advance into traditional alpine environments.
Achieving equal representation of landscape conditions within parks can be difficult, especially for biologically rich environments because of conflicting interests such as resource extraction, urbanization, and agriculture practices [77].In fact, climate change begs the question of whether representativeness of parks should continue to be a goal in parks planning [67].Southern Vancouver Island, for example, has a unique Garry oak (Quercus garryana) ecosystem with the greatest high priority species counts in British Columbia [78], but protection of this ecosystem is limited by urbanization.Evaluation of present protection coverage shows capricious results with some biological characteristics lacking adequate protection.Biases exist, especially in large reserves over 10,000 km 2 , toward protection of more seasonal, higher elevation, lower biodiversity environments [4,79].Critiques of protected area systems identify a bias toward low biodiversity regions that protect barren or glacial environments [4,8,19] rather than southern forests and riparian habitats that host many more species at risk [80].As vegetation on the landscape becomes more homogenous, parks will be more representative of the overall landscape.However, a likely cost is the loss of species vulnerable to change and species with specific habitat requirements [9].
The CCCma climate models have been shown to have consistent results, when compared to 10 global climate models, and are calibrated for Canada [63].However, there are climate models, such as the GFD GCM (Geophysical Fluid Dynamic, Generalized Circulation Model), that predict less precipitation in the future [81] and would forecast different vegetation trends.In the future, it will be beneficial to expand modeling of future DHI spatial distributions using a wider range of future climate scenarios and even updated climate scenarios.These mapped trends in future DHI are therefore not viewed as the final solution, rather highlight possible scenarios of the future distribution of indirect indicators of biodiversity.One of the strengths of this research is that it can be used to highlight locations where large amounts of change are possible and future benefits will be obtained by integrating this research with other results to build consensus models [82].
Our method of assessing change to individual parks and ecosystems using remotely sensed biodiversity indicators (i.e., fPAR metrics) provides a broad understanding of the vectors of change.By using indirect indicators of biodiversity we avoid the problems and assumptions associated with species specific niche-based modelling [83].It is unrealistic to represent complex ecological processes such as adaptation, tolerance, and migration [84,85] over a broad scale.Indirect indicators of biodiversity provide a broad lens for identifying where and how biodiversity may change and should be considered as a tool for strategic-level planning and resource prioritization.

Conclusions
Our remote sensing-based approach, modelling indirect indicators of biodiversity, provides information on possible changes to vegetation support for biodiversity in parks, which is essential for adaptive management.We have identified protected areas and ecosystem assemblages that are most threatened by climate change, such as Strathcona Provincial Park and the Mountain Hemlock Ecozone.Coast Mountains, Vancouver Island, and southern regions of the province are of particular importance because of their projected large change and the need to preserve the region's high concentrations of endangered species.Overall, trends in biodiversity indicators suggest future parks and protected areas may better represent landscape conditions.Climate change will lead to increased vegetation productivity and reduced seasonality in many alpine parks, which may also improve the biases currently observed in how well parks represent the range of landscape conditions.However, the cost of more representative parks may likely cause negative impacts to sensitive habitats and species due to their inability to adapt to changes in climate and ecosystem functions [9].
In the future, landscapes with high greenness and mid-to low-level seasonality and minimum cover will be poorly represented by parks.These are located in the valleys of the Rocky Mountains and central interior.Planning new protected areas should be focused around identified gap regions which are associated with Interior Cedar Hemlock and Sub-Boreal Spruce BEC zones.
With repeated coverage and spatially continuous data, our research, linking remote sensing-based biodiversity indicators to historical climate data and future climate models, responds to repeated calls for better information on future park conditions to aid adaptive management.The next phase is to better understand how to integrate scenarios of future landscape conditions into adaptive management to meet long term conservation goals.

Figure 1 .Figure 2 .
Figure 1.(a) Changes in greenness representation showing present day, B1, A1, and A2 2065 climate scenarios.Overall changes indicate a better representation for higher greenness due to climate change.(b) Changes in seasonality representation showing present day, B1, A1, and A2 2065 climate scenarios.Overall changes indicate a better representation for lower seasonality due to climate change.(c) Changes in minimum cover representation showing present day, B1, A1, and A2 2065 climate scenarios.Overall changes indicate a better representation for higher minimum cover due to climate change.

Figure 3 .
Figure 3. Ranking of highest changing greenness in large, medium, and small protected areas.

Figure 4 .
Figure 4. Ranking of highest changing seasonality in large, medium, and small protected areas.

Figure 5 .
Figure 5. Ranking of highest changing minimum cover in large, medium, and small protected areas.

Table 1 .
Ranking changes to the Dynamic Habitat Index within each Biogeoclimatic Ecosystem Classification (BEC) zone.