Homehttp://conserveonline.org/Not FoundHTTP/1.1 200 OK Server: Zope/(Zope 2.10.5-final, python 2.4.4, linux2) ZServer/1.1 Plone/3.0.6 Date: Fri, 05 Dec 2008 01:53:57 GMT Content-Length: 2998784 Content-Disposition: attachment; filename="Integration of Conservation Portfolios with SITEsv3.doc" ࡱ> .@ bjbjFF R,,8*|T*(...555{}}}}}}$Rh5..m8..{{J. PUt <**j**5!7555HnM< vnMAutomated Integration of Aquatic and Terrestrial Conservation Areas in Conservation Planning: A New Method Michael Schindel Oregon Chapter The Nature Conservancy Abstract The integration of conservation area designs for aquatic and terrestrial species has been a challenge for planners. The difficulties of crafting a single suitability index which reflects the landscape condition relative to all species, terrestrial and aquatic, and choosing assessment units appropriate for both realms, have made integration especially problematic. Here I will introduce a new technique, vertical integration, which allows planners to analyze aquatic and terrestrial targets simultaneously by using separate layers of assessment units, crafted to match the natural boundaries of the targets being assessed, with suitability indices incorporating impacts specific to those targets. Conservation Areas identified through a vertically integrated solution are efficient, and offer specific information about where to capture each target group individually, and where to capture them together. Examples of the use of this technique in planning will be presented for the Pacific Northwest Coast and the Alaska-Yukon Arctic bioregions. Key words: Planning, Conservation area design, Integration, ANWR Introduction Planning has emerged as one of the fundamental steps in conservation, with the aim of making sure conservation targets are protected with maximum efficiency. As planning has evolved, a wider variety and mix of targets have been brought into planning. Some of the earliest conservation plans focused only on imperiled species. Later plans have focused on all known species, and/or vegetation types. Most recently, aquatic systems and ecological processes have been included in the targets and goals of conservation plans. Whatever the targets may be, some sort of simulated annealing automated site selection algorithm (e.g., SITES, MARXAN, Kirkpatrick et al. 1983, Otten et al. 1989 ) is commonly used to create a map of conservation priority area (Andelman et al. 1999, Ball et al. 2000, Possingham et al. 2000). One of the biggest challenges in planning is how to incorporate aquatic and terrestrial targets into a single suite of conservation areas. Some plans have analyzed terrestrial and aquatic species and systems separately then attempted to merge the results manually. Others have analyzed both target types together in one layer of assessment units (AUs) and allowed the computer to find an optimal solution. A third approach is to merely overlay the outputs of each assessment. All of these approaches have serious shortcomings. Manual integration may be feasible for small areas, but large-scale planning efforts often cover millions of hectares. It is simply impossible to synthesize enough information to ensure reasonable outcomes. Analyzing both aquatic and terrestrial realms with the one-layer approach pushes a large portion of the solution into sub-optimal territory for each set of targets. Site selection algorithms look at the world through the lens of a suitability index which incorporates a combination of factors such as road density, percent land conversion or monetary value. An index crafted for an aquatic species will have little relevance for upland terrestrial systems. Similarly, an index crafted for both realms will tend to mask impacts specific to a single realm. The simple overlay of the independent assessments is perhaps the most robust solution, but often leads to massive conservation area designs, as identifying areas where it makes good sense to work on both aquatic and terrestrial systems at the same time is not an overt criteria. Many opportunities for efficiency will be overlooked. Site selection algorithms require that all species and System information for a planning area be attributed to wall-to-wall coverages of AUs, usually small-scale watersheds or hexagons of several thousand hectares. A computer then examines millions of AU combinations, and chooses the best combination from among them that meet the goals at the smallest cost. This cost is the combination of the sum of the suitability index for all the selected AUs and the sum of penalties for not achieving desired goal levels, combined with the sum of the boundary length, a measure of the outer perimeter of all selected AUs. Boundary length is directly proportional to fragmentation. A conservation area design comprised of many small, isolated patches will have a larger boundary length than one comprised of fewer, large patches. The best output of the site selection algorithm then becomes the departure point for human planners to review, modify and craft a final suite of conservation areas. In this paper I will introduce a technique that allows planners to analyze different target types simultaneously by using multiple layers of AUs crafted to match the natural boundaries of the targets being assessed with suitability indices incorporating impacts specific to those targets. This technique, vertical integration, enables planners to identify conservation areas which capture the best locales for each target group, while simultaneously looking for efficiencies by seeking overlap in areas where multiple target types may be effectively conserved at once. The Mechanics of Vertical Integration One of the focal concerns of conservation area planning is the design of protected areas that are not so fragmented that even though they may contain all of the species of interest now, over time the species are likely to disappear as a result of habitat fragmentation. In order to address this concern, automated assessments utilize the length of the conservation area perimeter to apply a penalty for fragmentation. Groupings of contiguous AUs have a shorter total perimeter, as the edge/area ratio is smaller than in a Conservation Area comprised of isolated AUs.   The ultimate Conservation Area would be circular to maximize area while minimizing boundary. Site selection algorithms utilize a boundary modifier option to modify the clustering in a conservation area design. This works by altering the penalty for fragmentation. As the computer examines possible AU combinations, the tendency to prefer solutions with contiguous groupings of AUs increases as the boundary modifier is increased. In our new analytical technique, vertical integration, the boundary relations between AUs are used to allow the model to recognize that two or more polygons stacked upon each other are also adjacent. In these situations, the model attempts to minimize the length of the total solution boundary by clustering vertically through a stack of AUs. If the boundary modifier is set to 0, the solution will pick the minimum number of AUs from each layer to meet the goals with no regard for adjacency. As the boundary modifier is increased, the importance of clustering, horizontally as well as vertically, is increased. This 3 dimensional approach mimics GIS analysis, though no spatial analysis is involved in the selection algorithms.  A major advantage of vertical integration is that it frees planners from using the same AU polygons for all targets. It is often quite useful to use polygons which more closely match the natural expression of a target type. Aquatic Systems, for example, are often classified as nesting polygons of increasing watershed size. Tributary and headwater drainages (Class 1) nest within small river drainages (Class 2), which in turn nest within large river drainages (Class 3). These classes of watershed can all be represented by polygons depicting their full contributing area. Their nesting is utilized with the vertical analysis so that each polygon is aware of all the polygons contributing to it, or which it contributes to. This larger, landscape scale context is a key advantage of this technique. The selection of the larger watersheds is greatly influenced by their attraction to basins with a greater selected proportion of their constituent tributaries. Techniques which rely on only one layer of AUs will often only select isolated reaches with no regard for their relation to the larger stream network. Multiple AU layers allow specific, relevant information for each target group to be factored into the suitability index. In the one-layer approach, a single suitability value was expected to account for all conditions which may impact any target group. This works well for pristine or heavily degraded AUs with similar degrees of impact to all targets, but fails where impacts are specific to one target group. A fish hatchery, for example, may threaten a wild salmon stock, but present no danger at all to a ridgeline plant species. In the one-layer approach, several AUs can have similar suitability values, but each may be inappropriate for one target group while well suited for conservation of another. With vertical integration, the aquatic suitability index can factor in the hatchery, while the terrestrial index is free to ignore it. The majority of AUs for any target group are interchangeable in that many different combinations of AUs can meet similar proportions of goals at similar costs. Vertical integration attempts to maximize the overlap between layers, allowing the site selection algorithm to actively seek efficiencies while maintaining the discrimination to avoid sites where conditions are unsuitable for a specific target group. The outputs from a vertically integrated solution offer more specific information about the conservation area design. Where does it make sense to capture all targets or to capture target groups individually? Methods Site selection algorithms utilize a file that contains the lengths of shared boundaries between adjacent AUs to determine how to cluster AUs into conservation areas. This file is the key to the proper functioning of vertical integration. Lets examine the simplest vertical integration, two spatially identical AU layers: one for aquatic and another for terrestrial targets. The length of each boundary between all adjacent terrestrial AUs is measured. These relations are then stored in the boundary relations file. The aquatic AUs will then be related to the terrestrial AUs they overlap. In this case the aquatic and terrestrial AUs are spatially identical; the length of their shared boundaries could be measured as the area of the polygons, or set at some synthetic value. We will initially set all of the aquatic-terrestrial boundaries at the mean of the terrestrial to terrestrial boundaries, so the model will generally be as likely to clump upwards through the stack as from side to side within a layer. These relations will also be stored to the boundary relations file. Two components will be part of the complete boundary relations file; the traditional boundary relations between the terrestrial AUs, and the relations of the aquatic AUs to the terrestrial AUs they overlap. An iteration of a site selection algorithm analysis begins with any "locked in" AUs that should be part of any conservation area, and a partial random selection of additional AUs. All selected AUs will then be scored for how well they meet target goals, the total cost of the solution, and total length of boundary. All exposed boundaries of selected AUs are included in the boundary length score. In vertical integration, those exposed boundaries will also include the values relating a selected AU with other non-selected AUs above or below it. For example, we are using 2 layers of AUs stacked in our analysis. If a terrestrial AU and the aquatic unit above it are both selected, there will be no penalty in the vertical plane, while a terrestrial unit selected without any corresponding aquatic AU would accrue a penalty. Similarly, the aquatic AUs would accumulate penalties for the unselected terrestrial AUs beneath them. Solutions which maximize the overlap between AU layers will be favored by the algorithm. However, the algorithm is not forced to select overlapping AUs in all cases. If the costs of an AU are prohibitive, or if the conservation targets in an AU are no longer required to meet goals, the algorithm can choose to forgo its selection even when the unit above or below it has been selected. Another advantage of utilizing the boundary for integration is that the weight of the boundary penalty can easily be changed in the site selection algorithm. If the weighting is set to 0, AUs required to meet the goals at a low cost are selected without regard for adjacency. At low weightings the effects of clustering will begin to be seen. A high weighting will clump so tightly that virtually every selected terrestrial AU will correspond with a selected aquatic unit, and the patch size of the terrestrial conservation areas will increase dramatically. It is important to remember that as the weighting increases more extraneous AUs will be selected merely to reduce the exposed boundary of the conservation area. Iterative runs, with increasing boundary penalty weightings, will allow the planning team to select the level at which clustering is appropriately balanced with the size of the total conservation area. Planning in the Alaska Yukon Arctic The Alaska-Yukon Arctic (AYA) bioregion covers more than 30,000,000 hectares, stretching from the Brooks Range to the North and South Slopes along the Chukchi and Beaufort Seas. The Arctic National Wildlife Refuge (ANWR) and the National Petroleum Reserve-Alaska are included within the planning area. This was the first large scale conservation area plan to use vertical integration as an analytical method. A simple, two-layer approach was used for transparency and ease of trouble shooting. All terrestrial target data was attributed to 5,000 hectare hexagonal AUs. The aquatic systems were maintained in their native polygonal format, watersheds ranging from 40 km to 100,000 km in size. Suitability indices were crafted for the terrestrial and aquatic AU layers, and boundary relations were built between adjacent terrestrial AUs, as well as between overlapping aquatic and terrestrial AUs. Because of the remote nature and large size of this bioregion, comprehensive data was not available for all of the aquatic or terrestrial species across the planning area. The data available for the terrestrial realm included life-history maps for 85 species, including data for each of the 4 caribou herds within the planning area, and species richness for 89 bird species, 29 terrestrial mammals and 438 vascular plants (The Nature Conservancy Alaska, June 2004). 36 Terrestrial Ecological Systems were also used (Jorgenson, M. T., et al. 2003). For the aquatic realm, 36 aquatic systems were mapped across the region, and anadromous fishes (where available) were attributed to the appropriate aquatic AU (The Nature Conservancy Alaska, March 2004). Iterative analysis runs were performed to establish appropriate boundary and suitability ranges for the two AU layers. The first run, for example, had a perfect coincidence between selected aquatic and terrestrial AUs. Looking at the output table containing the components of the objective function, it was revealed that the boundary values were swamping the suitability values in the equation. Multiplying the suitability values for all AUs by a factor of 10 gave them parity with the boundary relations in the objective function. Boundary and suitability values in the objective function were well balanced for all subsequent runs. Three goal-level scenarios were analyzed by the site selection algorithm, 25%, 50% and 75% capture of the representation of all targets (Figure 3). The vertical integration performed well, meeting a very substantial proportion of target goals at reasonable cost, with a large degree of overlap between selected terrestrial and aquatic AUs. However, at the 75% goal level so much of the planning area was selected there was no benefit from the vertical integration; a very high proportion of all AUs were required to meet a 75% goal level for all targets. Because the bioregion is relatively pristine and many of the targets are free to range over large expanses, a nearly infinite number of AU combinations can meet goals. For these reasons, no comparison was performed with a single-layer analysis for the entire bioregion. However, some areas do contain disproportionate amounts of habitat for some species. The Beaufort coastal plain, including the Arctic National Wildlife Refuge (ANWR) and the National Petroleum Reserve-Alaska, is just such a region. A one-layer analysis was performed just for the terrestrial targets within the Beaufort coastal plain, specifically to measure the relative importance of the biological resources of the Teshekpuk Lake portion of the National Petroleum Reserve-Alaska to the wider coastal plain. 11% of the 85 primary species included in this assessment are heavily represented in the Teshekpuk Lake area, and of those, 5 species are very dependant on the area for life-history stages. Black brant, Canada, white-fronted and snow geese exclusively use the Teshekpuk Lake region for molting, a crucial life-history phase. Almost half of all black brant nests are also found there. The local caribou herd only utilizes this area for mosquito relief during the summer months, and over half of their caribou calves are born there (The Nature Conservancy Alaska, August 2004). The Beaufort coastal plain is also very important for many of the aquatic species within the planning area. The majority of deepwater lakes, which provide over-wintering habitat for many fish species, estuaries, and the connections to the marine environment required by anadromous fishes, are found within this portion of the planning area. Various scenarios could also be studied with this data and the vertical integration technique. AUs corresponding with proposed development, for example, could be locked out, preventing the algorithm from selecting them. Subsequent outputs would then show which species can no longer meet specific goal levels with those AUs removed, or how an optimal reserve design might shift. Extreme seasonal weather fluctuations modify the regions habitats, and the species that use them. Seasonal reserve designs might also be a useful conservation tool.  Planning for the Pacific Northwest Coast The Conservation Area analysis for the Pacific Northwest Coast (PNWC) was the most exhaustive integrative exercise yet attempted in Conservation Area planning. The targets were broken into several groups, terrestrial, estuarine, freshwater aquatics (3 size classes), and near-shore marine. Assessment units were crafted for each group and separate suitability indices were calculated for each. Each target group was analyzed in a stand-alone fashion to see what the ideal automated solution might be for that group. All target groups were then run in a vertically integrated analysis, the solutions decomposed into their constituent layers and compared back to their original stand-alone runs to gauge the sacrifices made by any target group to accommodate integration with the others. Iterative runs also allowed us to weight the groups appropriately (by scaling their suitability and boundary values) so no one target group was dominating the outcomes. The final Conservation Area design met goals for virtually every target, with all targets having an influence in the outcome. The terrestrial group was attributed to small-scale watersheds approximately 2,500 hectares in size. These were chosen because they cover the full extent of the ecoregion, and make ecological sense to many of our partners and reviewers. The aquatic group was represented by three classes of nesting polygonal watersheds, tributary and headwater drainages less than 100 square kilometers (Class 1), small river drainages between 100 - 1000 square kilometers (Class 2), and large river drainages more than 1000 square kilometers (Class 3). These three classes of watershed were all represented by polygons depicting their full contributing area. The Class 3 polygons contain the Class 1 and 2 polygons contributing to them, and the Class 2 polygons encompass the Class 1 polygons which contribute to them. Some watersheds don't drain into others, for example, when a small coastal creek flows directly into the ocean. For the vast majority of watersheds, however, this nesting was a key to the analysis as each polygon was made aware of all the polygons contributing to it, or which it contributed to. The near-shore marine AUs were line segments corresponding to reaches of shore-zone habitat; unique combinations of substrate, wave exposure, and biotic assemblage. Estuaries were represented by polygons. In the US portion of the ecoregion, those polygons were defined by salinity zones and estuarine vegetation. On Vancouver Island they were merely polygonal depictions of the extent of each estuary. Vancouver Island estuaries tend to be quite small, as they often occur at the heads of narrow fjords, and are fed by smaller streams. To give our model the context to discriminate between these estuaries the sum of the shore-zone habitats intersecting each was attributed to the polygons. Each of these planning unit layers had suitability information tailored specifically for the targets within them. Each group was run in a stand-alone analysis, with the "best" output of each (10 runs, 5,000,000 iterations each, boundary modifier 0.1) saved as the benchmark to gauge future solutions during the integration process. All target layers were combined into one analysis using the "vertical integration" technique. We had earlier determined that a boundary modifier of 0.1 was optimal to achieve appropriately sized clumps in our terrestrial solution without many extraneous AUs. However, we wished to ensure that the overlap between layers was maximized in the integrated solution without sweeping lots of extraneous AUs into the solution. Increasing the boundary modifier would unfortunately have that effect. Instead, we held the boundary modifier at 0.1, and increased the boundary values between layers in the boundary relations file. The initial boundaries between layers were set at the overlap of AU polygons in hectares. Iterative runs were performed and the inter-layer boundaries were increased at each run by 20%. This iterative process was repeated until the costs of one of the constituent solution layers began to spike. The run previous to that spike, in this case the fourth iteration, was then used to identify the integrated Conservation Area. As the values of the boundaries between layers increased, the area of overlap between layers also increased, while the costs of the solutions remained fairly flat. The solutions were shifting to allow targets, for which multiple combinations of planning units at similar costs could meet goals, to accommodate integration.  As a comparison, all targets were attributed to a single layer of AUs for a traditional one-layer analysis. Suitability values for these AUs were set at the average of the corresponding terrestrial and aquatic AUs suitability scores. All other weightings and settings were held constant. The outputs for both scenarios were compared for the Olympic sub-section of the PNWC assessment. This subsection was chosen at it had the tightest coincidence between the aquatic and terrestrial sub-sectional boundaries (Figure 4). The goals were generally met well by both analyses (Appendix 1). The combined footprint of the vertically integrated terrestrial and aquatic conservations area within the Olympic Peninsula sub-region was 521,677 hectares. The footprint for the one-layer conservation area was 558,202 hectares, 7% larger (Figure 4). A better comparison of the performance of the different analysis may be the Elwa River. The vast majority of the Elwa River drainage is within Olympic National Park. The uplands surrounding the river are in exquisite condition with large tracts of old growth forest. Unfortunately, one of the largest dams in the PNWC sits low in the watershed, providing hydropower for the region. 25,000 hectares of the Elwa watershed were selected by the one-layer methodology, applying all of the aquatic systems they contained towards the goals. In the vertically integrated approach, none of the Elwa appears in the aquatic solution, but large portions are in the terrestrial portion of the solution. This is a classic example of the blindness of a suitability index crafted for all targets to appropriately assess impacts for an individual target group. Discussion Vertical integration, since its inception 2 years ago, has been used by several planning teams in the United States and Canada. Aquatic planning teams, most specifically, have found it beneficial because it has solved the problem of connectivity. In the one-layer approach, no AU is aware of any other AU it does not touch. The one-layer approach is inherently unable to link many contributing watersheds together to form continuous aquatic conservation areas. Because the vertical integration technique, when used with nesting watersheds, creates relationships between larger size classes and all of their smaller contributing watersheds, it is able to build these connections. There are, however, some considerations a team must be aware of when attempting to utilize this technique. An automated portfolio is a mathematical solution for a conservation area design problem. Planners must realize that any automated output only represents the solution with the smallest value of the objective function. The numeric value of the objective function is largely a dynamic tension between the sums of the suitability scores and sums of the boundary penalties. If either factor is weighted too heavily it will dominate the outcome. Therefore, planners are urged to look at the tabular outputs of their analysis, specifically the component values of the objective function. If, for example, a solution has nearly perfect overlap between selected terrestrial and aquatic AUs, the boundary values in the objective function will probably far exceed the suitability scores. In this case, the team may also notice the aquatic targets are generally far exceeding their goals. Similarly, AU layers with greatest relative costs will have the largest impact on the value of the objective function, and therefore have the greatest influence on the conservation area design. Therefore, when designing your analysis, layers which contain the most robust information, or layers of special conservation interest, may be weighted more heavily to allow them more influence in the outcome. The actual influence a layer has can somewhat be gauged by the cost shift of the other AU layers when compared against their stand-alone runs. In the PNWC analysis, for example, we didnt want the shore-zone segments to have a very large influence on the conservation area design. Costs were scaled down relative to the terrestrial and aquatic AUs. During the integrated runs, the costs of the aquatic and terrestrial components of the vertical solution showed no significant difference compared to the stand-alone runs for those layers. The shore-zone portion of the vertical solution, on the other hand, showed an average 23% increase in costs compared to its stand-alone counterpart, after accounting for the shore-zone cost scale change. The shore-zone component of the vertical solution was being forced into less favorable areas to accommodate integration with the aquatic and terrestrial layers. The geometries of AUs can also greatly influence the outcome. Hexagons, for example, cluster much more easily, and at lower boundary modifier levels, than irregularly shaped AUs like watersheds. In all cases I recommend planners build the terrestrial AU boundaries, experiment with ranges of boundary modifiers and suitability values that produce reasonable outcomes, then create boundary relations and suitability values for other layers based upon the ranges established in the terrestrial analysis. With aquatic analysis, we have learned that using synthetic values for boundary lengths is beneficial. Because any watershed of a type counts toward goals as much as any other watershed of that same type, area need not be a factor in the boundary relations. Basing boundary values for the aquatic layers upon the mean of the terrestrial boundary lengths produces a more robust integrated solution. For example, if the mean of the terrestrial boundaries is 3000, any Class 1 to terrestrial AU, Class1 to Class 2, or Class 1 to Class 3 boundary relation should be set near 3000, Class 2 to Class 3 relations perhaps twice as much. In the stand-alone aquatic solutions, the suitability values can then be scaled up or down until they are appropriately balanced against the boundary values of the objective function. Class 2 and 3 suitability values need not necessarily be the sums of their constituent Class 1s, they can be scaled independently such that the average Class 2 is twice the cost of the average Class1, and the average of the Class 3s three times the cost of the average Class 1. Iterative runs, with careful scrutiny of the objective function constituents, and goal attainment of the solutions, will assist the planner in achieving the appropriate balance. Linear features, like the shore-zone habitats used in the PNWC analysis, can also be used as layer in a vertical analysis. As line features have no true area, boundary relations should be proportional to the length of the segments intersection with other AU layers, and scaled to be appropriately balanced against those other layers. An early criticism, partially based upon fears that vertically integrated solutions would be less efficient, was that if targets are split between multiple AU layers, the algorithm would only receive credit for that portion of the targets in the selected AUs. In other words, if an area is selected only for terrestrial targets, and conservation resources will be applied to those targets, wouldnt the aquatic resources there also benefit, and therefore shouldnt they be counted towards goals as well? As our Elwa example demonstrates, it is not necessarily advantageous to count all targets which occur on the landscape every time an AU is selected. In fact, this is a chief failing of the one-layer methodology; areas are often selected for one group of targets that may be unsuitable for another group. Additionally, the specificity of the vertical outputs is very useful information. The overlap between terrestrial and aquatic solutions is the area where it does make sense to work on both target groups. AUs which appear in only one portion of a solution may have management and conservation strategies applied to them which are specific to those targets. In a world where human and financial resources are tight, tailoring conservation solutions efficiently and appropriately is paramount. Stand alone analysis for terrestrial and aquatic realms are valuable exercises in themselves. They reveal patterns of biodiversity, possible conservation opportunities for targets, and help identify threats to those same resources. If considerations for integration, like ranges of cost and boundary values, are used when building the stand alone analyses, the boundary relations between AU layers are the only additional piece required. All other tables can be cut and pasted together with no additional modification. This is much easier than having to rebuild all data from scratch to fit all targets into a single AU layer. Finally, it should be noted that any automated output is only as good as the information the algorithm was given. Data is a snapshot in time, often a snapshot taken 10 years ago. Peer review of any automated output is critical if we wish the conservation area design to truly meet the needs of the targets over time.  Conservation TargetAmount AvailableGoalProportion of Goal Captured by "One-Layer" Proportion of Goal Captured by "Vertical Integration" Astragalus australis var olympicus95140.000140.000Astragalus microcystis22100.000100.000Carex pluriflora33100.000100.000Cimicifuga elata11100.000100.000Dodecatheon austrofrigidum11100.000100.000Pellaea breweri22100.000100.000Plantago macrocarpa87114.286114.286Saxifraga tischii22100.000100.000Sparganium fluctuans21200.000200.000Synthyris pinnatifida var lanugino1917105.882105.882Accipiter gentilis169111.111144.444Ardea herodias fannini21100.000100.000Brachyramphus marmoratus670338126.627111.243Dicamptodon copei518412.500450.000Euphydryas chalcedona perdiccas1513107.692100.000Falco peregrinus2813146.154146.154Haliaeetus leucocephalus19767179.104150.746Hemphillia burringtoni3110150.000100.000Hemphillia glandulosa glandulosa335300.000240.000Histrionicus histrionicus514800.000800.000Icaricia icarioides blackmorei86116.667133.333Incisalia mossii mossii21200.000200.000Lycaena mariposa charlottensis21100.000200.000Oeneis chryxus valerata108125.000112.500Parnassius smintheus olympianus1313100.000100.000Plebejus acmon spangelatus21200.000200.000Plethodon vandykei209166.667144.444Progne subis31200.000200.000Rana cascadae44100.000100.000Rhyacotriton olympicus7624225.000237.500Speyeria zerene bremnerii54125.000125.000Strix occidentalis caurina232119125.210113.445Oncorhynchus gorbuscha12294136882173.99873.431Oncorhynchus keta pop ?1555327776640.46934.914Oncorhynchus keta pop ?2289526868658.84643.953Oncorhynchus keta pop 4227973068391911.57712.091Oncorhynchus kisutch pop ?69849820954981.26847.260Oncorhynchus kisutch pop ?1953219585966130.269135.326Oncorhynchus kisutch pop 1469883914096521.2662.072Oncorhynchus mykiss pop ?1155963346789146.647157.535Oncorhynchus mykiss pop ?452456135737125.86978.464Oncorhynchus nerka344003440081.72893.534Oncorhynchus nerka6107610799.99599.995Oncorhynchus nerka8407584075100.000100.000Oncorhynchus tshawytscha309270492781151.34056.318Oncorhynchus tshawytscha1042244312673122.810145.621Oncorhynchus tshawytscha486454145936158.147182.357Oncorhynchus tshawytscha1999229996196.82472.409Salvelinus confluentus13522367611160.153168.345North Pacific Coastal Herbaceous Bald And Bluff 233600.000500.000North Pacific Dry And Mesic Alpine Dwarf-shrubland And Meadow227492275718.683648.902North Pacific Hypermaritime Sitka Spruce Forest29579588739124.509133.719North Pacific Maritime Dry-mesic Doug Fir-western Hemlock Forest19596558790170.906135.374North Pacific Maritime Wet-mesic Doug Fir-western Hemlock Forest24184172552148.667136.164North Pacific Montane Riparian Woodland And Shrubland33100.00066.667North Pacific Mountain Hemlock Forest12500325001320.145316.606North Pacific Western Hemlock-silver Fir Forest19680739361312.743285.359Coast Tributaries - Outwash, Low Elevation, Moderate Gradient321163.63654.545Coastal Upland - Glacial Till, Low Elevation, Low To Moderate Gradient4214107.143135.714Coastal Upland - Sandstones, Low Elevation, Moderate Gradient401369.231100.000Olympics - Sandstones, High Elevation, High Gradient124275.000300.000Olympics - Sandstones, Mid Elevation, High Gradient3110190.000150.000Puget Lowlands - Outwash, Low Elevation, Moderate Gradient9560.00040.000Haliaeetus leucocephalus wintering area11100.000100.000 References Andelman, S., I. Ball, F. Davis, and D. Stoms (1999). Sites V 1.0: An Analytical Toolbox for Ecoregional Conservation Portfolios, a Manual Ball, I. R. and H. P. Possingham (2000). Marxan (V1.8.2): Marine Reserve Design Using Spatially Explicit Annealing, a Manual Jorgenson, M. T. and M. Heiner (2003). Ecosystems of North Alaska. Unpublished 1:2.5 million-scale map produced by ABR, Inc., Fairbanks, AK and The Nature Conservancy, Anchorage, AK. Kirkpatrick, S., C. D. Gelatt, and M. P. Vecchi (1983). Optimization by Simulated Annealing. Science, New Series, Vol. 220, No. 4598, 671-680 Otten, R. H. J. M. and L. P. P. P. Van Ginneken (1989). The Annealing Algorithm. Kluwer Acedemic Publishers Possingham, H. P., I. R. Ball and S. Andelman (2000). Mathematic methods for identifying representative reserve networks. In: S. Ferson and M. Burgman (eds) Quantitative methods for conservation biology. Springer-Varlag, New York, pp. 291-305 The Nature Conservancy in Alaska (March, 2004). Alaska-Yukon Arctic Ecoregional Assessment Update # 4: Freshwater Ecosystem Model. The Nature Conservancy, Anchorage, Alaska The Nature Conservancy in Alaska (June, 2004). Alaska-Yukon Arctic Ecoregional Assessment Update # 10: Decision Support Tool. The Nature Conservancy, Anchorage, Alaska The Nature Conservancy in Alaska (August, 2004). Alaska-Yukon Arctic Ecoregional Assessment Update # 11: Application of Ecoregional Data: Teshekpuk Lake. The Nature Conservancy, Anchorage, Alaska PAGE 1 PAGE 9 Figure 1: Both of these selections of AUs have the same area. The right hand grouping has a perimeter more than twice as long as the left grouping. Figure 2: A schematic demonstrating the boundary relations between stacked and horizontally adjacent AUs. Each AU must relate to all other AUs above or below it, and in some cases, from side to side. Appendix 1: Comparison of goal attainment between vertically integrated, and one layer site selection model for the PNWC Figure 4: Comparison of automated Conservation Area designs between vertically integrated and one-layer methodologies for the Olympic Peninsula region of the PNWC. Figure 3: Alaska-Yukon Arctic Assessment: Relative Biodiversity and automated solutions for all targets at 3 goal levels m " $ 5 9 B Y Z _ ` z { g h i " ) * . ˾о}u}u}u}u}u}hqhq5 hq5 hqhqh(C5B*phhvMd5B*phhe[Khe[K5B*phhq5B*phhe[K5B*phhe[Khe[K5B*ph he[K5 hD^5 h(C5 hG5 hvMd5 hb&5hb&hb&5 hO$5 hO$CJhO$.kl}   IJgd\&$a$$a$. 7 8 :    ]^_`jo$%NbǸǗ{vvqjq{ hMohMo hMo6 hO$6hdRh 6QhMo hO$5CJhO$ hO$5h\&h\&5h\&h\&5B*phh\&B*CJaJphh\&h\&B*CJaJphh\&5B*phhO$5B*phhe[K5B*phhe[Khe[K5B*ph hoH5hqhq5(FGHhjnu[\Qs"8|T-/bdeijFGOR he hO$5CJhD^B*phhO$B*phhO$B*phh]h 6QhD^h0 h0h0h_hxYh(EhMohO$hX{h"Nh$=@ !R T ! !.!0!*","l"""#########$$&$,$$$ %*%L%V%%%%%p&r&}&&&&|'}''''''''' (((((!(.(w(޶ުަަަުޢޢޢhaKh2hxYjch_Ujh_UmHnHuh0hehD^jhO$UmHnHuhO$(jhO$UcHdhdhdhۢjhO$U@ "#$%& ######((`+a+----22gd2#$d%d&d'dNOPQw(z((((((((((C)E)))+*,*D*E*N*P*R*U*****D+J+q+s+++++++,,--------A.C.a.c.....////)/-/./5/:/M////////00'030P0R0e0l0h~#h0Kh(E h2h2haheh]hahD^hxYh2h0hO$Ol00011m1w122222222222222222 3 32343333333A4C4i4n4|444444+5-5Y5[555663656q6s66677<7>7u7w7888888999995:]:^:::::::h~#h>rhe hO$6h]hD^hYhahq!2h0KhO$S28 8;;;;<<O?P?ABBBDDFFHHhgd%$a$gd%:;;;;;;_<`<<<<==P?m?|?@@@AA AlAmAB>BABBBCBLBUBBBCCCDTDvDxDDDDDREFFFFGGGGGSH\H]HgHhHHHHHHܾܾܾܾܶ²²¶®ܲ® h3hil_h hil_hyNh>rh$y;h'h Khph>h h3h% hFh%h% h~#h%h%5CJaJh~#h%5CJaJhO$@HHHHI I3IIIIIII&J(JSJbJvJwJJJzKKKKK L8L;Lrh%h'hphD h3hil_h`7VhyNhil_h >OOOOOOOOOOOO'T(TUUwXxXYY.[/[|\}\aacgdE"^"gdYtgdAs$a$gd%,P-PPP|QQQQQQQQQQQQQQQR+RRRSS_SlSrSwSSSSSSJUUUUUUUUUUU8V:V;VBVVVVVVVVVVVWW#W$WLWWWwXxXXXXXXXZZ{]}] ^ ^l^t^h)(hD^hYh(ChAshehq!2 hO$6hoHh=hhO$ hCJOt^^^*_,_____```aaaadde emhwh#x$xdyeyzz{{T|U|i|||||||||| } }(},}-}<}D}F}S}t}ÿ×h$y;hoVlh(Ch%hqf)hh h*h>h hYthEhE5CJaJhhEhYtjvhYtUjhYtUmHnHuhAs hAshAshh 3hD^hO$h)(hq!23cclhmhxhyhkknnttt{{T|U|ij݃ރ !"gdt1gd$a$gdEgdEt}}}}}}}6~C~D~M~l~m~~~~~~~~~ !"07Rn{~59GJL|hiȁUł͂΂q܃݃ރQRļ踴дشh_hgzhUXhph."h!h(ChYth#Nh>h*h/uhnh>hht1h&/&hqf)h$y;hDF%&'(*BChi 78qrƇLJ !WXވ߈EFyzډۉ ?@op݊ފ :;no֋׋ @hO$6CJOJQJh hO$CJh ht1CJhjht1UmHnHuhYtht1 ht1ht1O"#$%&()*>OT $$1$Ifa$$1$Ifgdt1gdYt ܅ޅNF;;;; $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a PH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a.02:BPH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6aBCTVX`hPH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6ahiPH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6aPH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a҆ԆֆކPH==== $$1$Ifa$$1$Ifkd *$$If64  rL*"H t0  #4 6a PH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a #%'/7PH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a78[^aiqPH==== $$1$Ifa$$1$Ifkd!*$$If64  rL*"H t0  #4 6aqrPH==== $$1$Ifa$$1$Ifkd(*$$If64  rL*"H t0  #4 6aƇPH==== $$1$Ifa$$1$Ifkd/*$$If64  rL*"H t0  #4 6aƇLJPH==== $$1$Ifa$$1$Ifkd6*$$If64  rL*"H t0  #4 6a  PH==== $$1$Ifa$$1$Ifkd=*$$If64  rL*"H t0  #4 6a !ADGOWPH==== $$1$Ifa$$1$IfkdD*$$If64  rL*"H t0  #4 6aWXilowPH==== $$1$Ifa$$1$IfkdK*$$If64  rL*"H t0  #4 6aPH==== $$1$Ifa$$1$IfkdR*$$If64  rL*"H t0  #4 6aȈˈΈֈވPH==== $$1$Ifa$$1$IfkdY*$$If64  rL*"H t0  #4 6aވ߈ PH==== $$1$Ifa$$1$Ifkd`*$$If64  rL*"H t0  #4 6a035=EPH==== $$1$Ifa$$1$Ifkdg*$$If64  rL*"H t0  #4 6aEFegiqyPH==== $$1$Ifa$$1$Ifkdn*$$If64  rL*"H t0  #4 6ayzPH==== $$1$Ifa$$1$Ifkdu*$$If64  rL*"H t0  #4 6aƉȉʉ҉ډPH==== $$1$Ifa$$1$Ifkd|*$$If64  rL*"H t0  #4 6aډۉPH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a ),/7?PH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6a?@[]_goPH==== $$1$Ifa$$1$Ifkd*$$If64  rL*"H t0  #4 6aopPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6aPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6aɊˊ͊Պ݊PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a݊ފ PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a &(*2:PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a:;VZ^fnPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6anoPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a‹ȋϋ֋PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a֋׋ PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a #+29@PH==== $$1$Ifa$$1$Ifkd +$$If64  rL*"H t0  #4 6a@A\cjqxPH==== $$1$Ifa$$1$Ifkd +$$If64  rL*"H t0  #4 6a@Axy%&]^WXĎŎ@AGH45Ց֑34ӒԒmnƓɓϓғٓߓ89:;Cl h CJ hdRCJ hcpCJ h$=CJ hO$5CJh_mhO$ hO$CJhhO$6CJOJQJhNxyPH==== $$1$Ifa$$1$Ifkd +$$If64  rL*"H t0  #4 6aό׌ߌPH==== $$1$Ifa$$1$Ifkd +$$If64  rL*"H t0  #4 6a%PH==== $$1$Ifa$$1$Ifkd +$$If64  rL*"H t0  #4 6a%&@GNV]PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a]^qw}PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6aPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6aˍэ׍ߍPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a PH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6a8@GOWPH==== $$1$Ifa$$1$Ifkd+$$If64  rL*"H t0  #4 6aWXqxPH==== $$1$Ifa$$1$Ifkd$+$$If64  rL*"H t0  #4 6aĎPH==== $$1$Ifa$$1$Ifkd++$$If64  rL*"H t0  #4 6aĎŎ܎PH==== $$1$Ifa$$1$Ifkd2+$$If64  rL*"H t0  #4 6a+.08@PH==== $$1$Ifa$$1$Ifkd9+$$If64  rL*"H t0  #4 6a@APH==== $$1$Ifa$$1$Ifkd@+$$If64  rL*"H t0  #4 6aˏҏ؏PH==== $$1$Ifa$$1$IfkdG+$$If64  rL*"H t0  #4 6a*17?GPH==== $$1$Ifa$$1$IfkdN+$$If64  rL*"H t0  #4 6aGHPH==== $$1$Ifa$$1$IfkdU+$$If64  rL*"H t0  #4 6aݐߐPH==== $$1$Ifa$$1$Ifkd\+$$If64  rL*"H t0  #4 6a$,4PH==== $$1$Ifa$$1$Ifkdc+$$If64  rL*"H t0  #4 6a45elrzPH==== $$1$Ifa$$1$Ifkdj+$$If64  rL*"H t0  #4 6ađǑΑՑPH==== $$1$Ifa$$1$Ifkdq +$$If64  rL*"H t0  #4 6aՑ֑ #+3PH==== $$1$Ifa$$1$Ifkdx!+$$If64  rL*"H t0  #4 6a34ruxPH==== $$1$Ifa$$1$Ifkd"+$$If64  rL*"H t0  #4 6aÒ˒ӒPH==== $$1$Ifa$$1$Ifkd#+$$If64  rL*"H t0  #4 6aӒԒ PH==== $$1$Ifa$$1$Ifkd$+$$If64  rL*"H t0  #4 6a[]_fmPH==== $$1$Ifa$$1$Ifkd%+$$If64  rL*"H t0  #4 6amnPH==== $$1$Ifa$$1$Ifkd&+$$If64  rL*"H t0  #4 6aDE”NHFDDDD^kd'+$$If64  PrL*"H t0  #4 6alm”z{ؕߕ uvjŗHr+Zżwshqf)hqf)0JmHnHu hqf)0Jjhqf)0JUhO$h h$y;CJ h>hCJh$y;h$y;CJaJ h$y;CJ hO$CJh h CJ hO$5CJhdRhO$CJhdRhdR6CJhdRhdRCJ hdRCJ h^\2CJ h CJ hcpCJ)”Ôz{ vwjk˜DE ,-.$a$gd_$a$gd gdYt$a$gdYt$a$gd(C$a$gdTh]h&`#$ؙ֙EMКҚٚښ-.68ƲhO$h_ h$tzh_ h_5hYt hYt5h0hD^ hqf)5hShqf) hqf)0Jjhqf)0JUh0JmHnHu/ =!"#$%DdT~l9  l  c HA$bound_length_compR}reoL.<DF}reoL.<JFIFCCl" }!1AQa"q2#BR$3br %&'()*456789:CDEFGHIJSTUVWXYZcdefghijstuvwxyz w!1AQaq"2B #3Rbr $4%&'()*56789:CDEFGHIJSTUVWXYZcdefghijstuvwxyz ?(?_/)go|w|-)(kE"#=7W#Uu?>>jw^U *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O74oC+ *o#O7?&W^/uYg|cCχ<+? sN񏃼YW犼'j7:?|5_k>/hZޗڥ667\,K\k_<+7 YgdE|F,%&G{Yd~ο?.lyGB͙? >_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~_K/Ig:ƣ~< ߴ5>< Uߴ5>< ߴ5>< Uߴ5>< ߴ5>< Uߴ5>< ߴ5>< Uߴ5>< r.Akυ|\G[j.O|7jSAsuYZ}-v Y?Ί`+< 6&x_jS%/W_/< G +^g[~)Ŀ;>EPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEP_X_R0#Y5}_X_R0#Y5')?h* )M_5S ~UR,?#P|EQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQE4?p/ #4!_T>*xHt?_ZLJ񷍥 _ ?tys~;4m;Fw-KWx^%\ZF gY ?ٻeSQ🆮uƞ6-GѮ|cx}:_]T(./k7y@E/Ef?h|߱C@w_5<}"~ǟTWQ]>_o; G>Q_4U:WM;oa/>?ow?6ïx]^l^Þ.OOƉ?^j/xN'?b0| >^O>Ś V^ < G[h(((((((((((((((((((((((((((A_)g )M_h+ Yfd:F\:_bYK??lg_`k_<~((((((((((((((((((((((((((((((((((((((((O'g#A cH4k|j?|y'*O-_#Ѧ|A.t~~ _}ҼG__El:WsrRgA5/*O*/SҴ&yu-c߁5bzƟc'< {~7>0xg <?OV?/~ocR,m~._U, Smbw;~&d- [CUV{p@?X(((+)Y~ߟGvx[j'mOןd^ EI>E迴ߪAs_@4??>-!|)ӼYE>5RҴKf~R' >5-(e~)|:/fEO|YwZ~ #_fkVX}캞e{UԿ_ + )@㳢x#Wn>%$~t4khx&6z}nt%"3 c|[OY?e?4}u/|M h-<]7]c>) JoxK࿇@QEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEWV/Eu_V/Eu@ uac Sय़GWj q_s?ǟ|bK?lh4QEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQExUK,gƽ[#K_rX-O7# f> x]+9xH~"mo#6:^o?ࣚߌ?4k6ZWif5[M3BmOx]kT|?cq/1З\imo/׋:xW'9 t:6\ψ4/g_gNK/z[ý2xSK53Qu:Wѵ{/VuKK}CL4BKK;Qk[Y[{ydh7e?=>&*oo|w~ x#6z_/axWJo</S.(ItmGGw_Pk3 /_?B_*߰g?zǍ? ?8<m}'Jj:>;YGh=#V;oCuo|+5s_g((((((((((((((((((((((((((A_)g )M_h+ Yfd:F\:_bYK??lg_`k_<~(((((((((((((((((((((((((((((((((((((((+߱a,h/?+X?}h>-/ǧ.|#⯇_hu?zt~,|R'5Re)|E/f;OԼYwú~ _#FfhZV_f{{,gqN~Y_(F/^.G?7m -> ~ xY<% (/z.<޳iw-?0]ƾ-N~O/\ga|:J4xo##]|5/ k^wewu?)%x#>~^l`.io$*)wjw?%{:E6ZmcQEQEQEQEQEx ?O_wMXz|NմbLn{kydt?/77ŏ3V_7FMZ]}ľ/aW7mW[-HԦf]:O'3'q/f>*v~o2iR:<^^+mI6zl2^O(m=gNK/z[ý2xn ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( (  (ـ~ο (ـ~ο?N5w3yAWV/Jo֊_?kgX?_)g?((((((((((((((((((((((((((((((((_T? .`ωxF*x\q#:Nk5dӴ}OvO<5ϚGbU>WhOũ~Wx?>*~zog߃+> i._k#wttORtKo4__ ݬ]Wj:^hzfui6^v~}vvai6V6VmiimVG 1h?#N:cxNlt/ovO4'[9~x~5M?W7˿L'9t[e~ ~՚u{K__-Dw!\ibxJ >p'W436b_basj8j^[/OB#<m'񯈾"3h(((((((((((((((((((((((((o+7 YgdE|FSय़GWjo (ـ~ο?N5w3yAW}X_R0#Y5')?h*_( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( (?oO"kwiFa ǎ{ -uxcw7'%cW %"? k.SO}(o&kǟ㟇Ѯ|Mn *r?^7? J"!OxDzO cDZ/ ?ిPM_|ZK|J>.|OPMհOko xM7< >þӴ?co h&/Nlֳ [ [xᯅtoxgJ#I7F.%t#GP%uIt. N[DӲG((((((:kfZuj7z^ꖖivawַ7Kmwis6 ѼnoQ@u?'?᎟[|$įZ鷞-ԷgB .4MKLwx~$t<KC gXڏ_?ǃV3IMGG>'k>#հxឡ8{k~5?4o χ"qo%4~g߶//ƈn" wtoVS]G]I5G5AaAkCQ0uMO:7ocg_2,kí-/Ksh:&]It:\/k665( ( ( ( ( (>Gd~?٫6v2|Aծ¾8#2Xߋcta Iu Wu%]L?SG?i#W#Dx_?Rx{C>xM_R-`H]Y5K kS֍g_ x[Rӿ׎_: _ <xh:?,o?8,in G/K׉5/|@"iC3]Ҵbe+ا+2U2|k~"h&;to'_~?f4}+K?|SxITxKCƯO־+_W'0/$k*K_ĭC=v⋭g:Lj,[zό<__( Kj+_F`' oqxk?n/>27Vgak^ 1 >@,SٟO{n<+KKucyWTt/L;Et#F,m4'I-3K4x,4:8mlla(cHT^(((((((]fi?3 -cj#7xrNе>(Vi7 5+?҂+?_įٯDWo[x7VOΈW^v:rXjFWd`7ih| W436b_basj8j^[/OB#<m'񯈾"3k;>|{?5<]; mwvڞy/~?kfM嵶?Znxմџ; G>Q_ GuYc?(./k7y@E/Ef?h|߱C@w_5<}"~ǟTWQ]>_o; G>Q_WGt_|A|h֓wď D]Kú_|+w>x#ZƩi*y@w'/=@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@|bYK??lg_|bYK??lg_`k_<+7 YgdE|FN5w3yAWV/Jo֊@EPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEPEP?xW^/u_/> ŚΝ O¾Ӯuľ!.l gwk:ΩwkizuK*Oh~;?'xK¿#5D⿁4Y<[w_a G'~Wž*|\?f>m׋:xW߇|Gq:6\7hxO9տUGBS[TԿ;]-h%;}#oC|(׊?u\j<_[6 Vk/XNm̟x+uY . o^'Ɲ&>׼gk<:JkqMsL<]WM/|7 jtBD?߃|5+m#X4/ƿtwj|3GLjae># c5g:WWx[¾+~O>tk>0Qw6H|_ ~ 9"SŞ|+_^|ig~1w4m;~g|G\!׉|=[^i"KtmRN4뫛i&?-RUK*­nKe=WLէR \.'x>79xmX31lfaW+w{P_[ĽKWkot}N?]x Vo5oWR[-(b(((((((K]%|C{.$mw}wsw}wZVAwkz]h_kz}aڥ۵ᚯ\ Gmu5ߺZ' ߅|9ؾ(Zu:o[V xJЀ>+ԿǽKPP~au猾"꺖nIwxGPcx: ^񭥤7<45/ro[ ^WW@Mы|v?I1o?.?y_P_9 7F-G/( &ž\~Q@$k#sE~UÐ?b]r?O_?Ə|'ϮxS? |/U}gA[xrk^hwSlԧmWp$'o1qo^烯~K~K]Z>$9;t/> ^3zKA N/Gů|1񮨺[4?xWǾ׎ _>4xzΝ? ¾#ӭx^!o4wyihΗwujuյ+UԿ_ + )@㳢x#Wn>%$~t4khx&6z}nt%"3 c|[OY?e?4}u/|M h-<]7]c>) JoxK࿇@QEQEQEQEQEQEQEQEQEQEWV/Eu_V/Eu@ uac Sय़GWj q_s?ǟ|bK?lh4QEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEW/7~4~o쏆_?/ c]ȵ? %?@ t|: #þWϿT?x/?1vvDӦӿ&3A /5:ΩyxTmNVbu+k+m˿{ƿĿ y/_ ٿ-N=MF>1Դxn;+(((( ['>!Ǧ&]kĿ~#j//TWl4}/^V':izfj:.X6cwNioi[ij6qk{c{k4v1Koso,>jq$o{mJ.wqk:o6 ω$3⻭g>5oK h>|m h(((((+)~G6D~in[m3k¾2=VF>Om_ɿPl/^߀9ˈuߋ>3\^ OJN Oů1Vn,? W^ׁ| /> ?x;z6 O¾ӭxk^ gi6ikizuE8?k_O3Q> GaK?,t?|9g? |K~YH|_[Yk dW3 b(((((((((((((((((((((((((A_)g )M_h+ Yfd:F\:_bYK??lg_`k_<~((((((((((((((((((((((((((((((((((((((+Eh~=~_#B?WD> kyÏ? o YE?) Ou|jgmWS0[.׆i_C/O^i4o`MoXҬeT/T53Nִ]FWѵ{MSIմ}CL4B; GNk[Y[{yc7V?EX/3|1|CŚǺ|~kV^xŞY6 Wa_g \տhsCҵ^ķ%7&zֿorݽLJ[m94o6^F(((((((7 OO AcYx?ΓqucN֚<ծwO /kQEQEQEQEQEx ?O_wMXz|NմbLn{kydt?/77ŏ3V_7FMZ]}ľ/aW7mW[-HԦf]:V%4;OO> m/S|Ck´էټ3ɤKuiFk~/?a:~"<5bg]ǃ/|K o]OKƑJ<;˙B?߶Oa^,%h?x+/h>/ G|.տ#⯈)迳GtϳtO(((((((((++ Yfd:F+ Yfd:F\:_X)R?#Z+5}uac Sय़GWjo( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( (?8?P/ 2`ώn7m w$ӧ$? 7g-_ŞvǿoxY|_zÍ3M|9Yy < vMW:]~ms4[~յ]Zr?].5]ZտY*l_ſGx~˟_.>?/ Z?oχxOWox W5o>]jt? %H hnuUu{s=IC։ w4-#ZϰhoV_k6A۫6ygh ( ( ( ( 񯋺^)_j3cc?6iqxᾱgTm'A7T^xGO>=H4zsh :^kfhj6iwviwvi$ַ7Esiwm,642@Q@Q@Q@Q@UST4=3QֵFHѴSVյK}?L>KGQ[++Xe+{kxi4fpw_S_xc!Bhz>SIx/oiw:q c_&H"l5=%f%7xs ?hO_V7S/|9ox⇍ wTn~ xo>}/ ]~^~^<DS|5] if/|Sxu9[x?> Xj.⋈y>4e?nVO!w5 zw>6??iy_%iz_/)g q_s?ǟQ@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@Q@V'টg_KX6ϋ>Amcᇉ[7^ ?ꚺ/σ;?/~ x?=sks6ɇ>*$:!h cn[χ<+ &|Hgt$ƥj׵xxX+?[_xK՞LJwheSMψ<% +]Z?w-:l>;oƏ\|g3| ]ߏ4;/xCPuP̶476,Wv7q\X[[}jտ_rǾ :ߎ?glKZ_|7$o%޿7-t/ wqYjzl?*(((((+ƾ.ꚞo|$}xu:uƑ}>/⏎tz M/z?KzPk#~i^2?W5M3C5kZltH5m[T/L仿u![bI#FaOT~-NҼg:KmFH}S6ڶKƭǞ¤="  _a/=m{Ş,𯀼+_xľ_|gŞ15;Þ|9\!/|C\-"TuRN[-U-_Z~'8v4~Z~M+S|Gi >׮?᾿e ~&w?ωz{x,Ws3LB8Kt~ ]◊|Ezn)mxT @7CJ76s$K+q>=J_"]U8|Qxgw'G&0i֡oq'QEQEQEQEQEQEp| iki`@tox?ĩgth7ksk$Eu\Zj>񾓥:$_״ub>_Ma:O;?EkZ:[F<;ۮfk\iww5[z?,~O% xkxUxq&itQEQEQEQEQ |O;_Ul/,5 k?!d6]3,|Wiib֒/4٥zhO?fh?uKoM4vE-6+uK۫;${~Z)>,~|7ҵ|uςnu.o%^ċ)|E|O; x*^i˪k:7~x_ƞ6-"/¾|C=\t_[ eDY:UPE _u誣#VY~οT~_F*տ_@o5EU_W,k?g_?5o%? gU@jK/ ??#oWž a ?wkLl x>?#?J?Z}EPEPEPEPEPEPEPEPEPEPEPEPEPEPEPX)R?#Z+57R,?#W|V/Eu_ uac (ـ~ο?N5w3yAPQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEQEW ?0x|S WƝO%>'?|U oW|O;:> M >&𭅗MCL6OA+?5o*o? gUGFM*>տ੿@`5WEU7誠#V?~_TjT +?0O[ DYEUQSk?hO+?5o*o? gU_o:?k oogϊ Ҿj? ᆣ{^/'}|fO|Ag DAIp(((((((fG\פk2_|_,|Q3Ex׈xþ5E=o >7ਾ)nx/ ^WOob~-σ$o 5G{-D5}^? ?Losυw-/O_'Zs~4 +}wöz'MVxV]f-_W4VQEQEQEQEQ\ax÷># /-O4x^<+OiL^'ִYԴP^4S?5_ﴝOm/cw%|Dfc^W4 |Dρ~8_^o% ~%~#m[NQLjn4[]fžY+ZWѾctX" ( ( ( ( ( ( ( 99 w v埈xw^ uk~ܺO'NlR=Q.,_gۯJG*|asgc- AZHy+^>(/ =L[ou{-.R]0%+ ߴq?jsĿ .[kxG_xgKmuxĶmz|);ѼM _xw( ?1L VlGƏT>>^fM泽i/t}Ÿ3֥:GM'OYY?XiG_)k2@_W׼Q3E/4jxþ"䗺E=((((((((,य़`Gj,य़`Gj0O5S ~UR,?#Wk_<+7 YgdE|F(((((((((((((((((((((((((hg9x#YK+i#𮻬{X5KIw /R~^_xv{5tk{}{KK9on|7Ox^y7E6~?|qY/- _>o_ i'$^K~.п|;i ??Ň>զjWAsBe/E*x?Tv>%if?(]nGK³}N-\3ê>Ҽ'? Rgu)m49'j~~ iN5Y]h+qxK?e[/NxFftL-aֶ_u~I~C~ΊH?io$?!}пG$_>_s<3[/H/ -Dg$? K3 AK%:(|k[Gޣž+/RK[U|A6uMIͧ߀W:.XԭKաFw c#"x ?mu% | cƿ Ju^{MΕ >G-|bӾ"k nٟL}tG-zׂ×ƅj?|DxL