{"id":7,"date":"2017-12-04T15:37:33","date_gmt":"2017-12-04T22:37:33","guid":{"rendered":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/?page_id=7"},"modified":"2017-12-05T23:23:25","modified_gmt":"2017-12-06T06:23:25","slug":"methodology","status":"publish","type":"page","link":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/methodology\/","title":{"rendered":"Methodology"},"content":{"rendered":"<p><b>Overview<\/b><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 To answer our three questions we used an MCE (multi-criteria evaluation) to determine the most favourable areas for grizzly movement and settlement. Using different types of land use and land features, both natural and human-made, as our criteria, and using the analytical hierarchical procedure to decide the weightings for use in the weighted sum calculation, this MCE will score the cells in our extent according to how favourable they are for grizzly movement and habitat. We can then highlight areas of interest from the MCE map, such as an extended patch of favourable terrain, which can serve as a potential wildlife corridor (Fig 1).<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-40\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Flowchart.png\" alt=\"\" width=\"279\" height=\"540\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Flowchart.png 279w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Flowchart-155x300.png 155w\" sizes=\"auto, (max-width: 279px) 100vw, 279px\" \/><\/p>\n<p><strong>Fig 1.\u00a0<\/strong>Flow chart of our methodological process to produce a habitat suitability and wildlife corridor analysis of southern B.C. grizzly subpopulations.<\/p>\n<p><b>DATA TRANSFORMATION AND NORMALIZATION<\/b><\/p>\n<p><strong>Combining DEM\u2019s<\/strong><\/p>\n<p><span style=\"font-weight: 400;\">ArcToolbox -&gt; Data Management Tools -&gt; Raster -&gt; Raster Dataset -&gt; Mosaic To New Raster<\/span><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 Our study area, which contains the majority of southern BC, is composed of grid squares <\/span><span style=\"font-weight: 400;\">82M, L, and E and 92O, P, J, I, G, H on the map.<\/span><span style=\"font-weight: 400;\"> To begin working with the data, we needed a base DEM to work with. We started by downloading the DEMs of each grid square and then using the \u2018Mosaic to Raster\u2019 function to produce a combined DEM of the region (Fig 2).<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-56\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/DEM-1.jpg\" alt=\"\" width=\"3300\" height=\"2550\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/DEM-1.jpg 3300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/DEM-1-300x232.jpg 300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/DEM-1-768x593.jpg 768w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/DEM-1-1024x791.jpg 1024w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/DEM-1-624x482.jpg 624w\" sizes=\"auto, (max-width: 3300px) 100vw, 3300px\" \/><\/p>\n<p><strong>Fig 2.\u00a0<\/strong>Map of the combined final DEM.<\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 Then to simplify the acquisition of other layers, we created an extent polygon. This would allow us to accomplish two things: creating polygon\/polyline layers that matched the extent of our DEM minimizes the amount of clipping transformations we have to apply to any subsequent data retrieved, and minimizing the chances of having a layer that exceeded DataBC\u2019s download size limit (Fig 3).<\/span><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 Our intent was to match the size of the combined DEM from the previous section. We did this by drawing a polygon in the Arcmap editor, to roughly the dimensions of our area of interest. The resulting polygon did not quite match the DEM, as it has 2 clipped edges (Fig 3). However, we do not believe that this flaw significantly affected the results of our project, since it only removed portions on the outer margins of our extent.<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-42\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Extent-Polygon.jpg\" alt=\"\" width=\"3300\" height=\"2550\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Extent-Polygon.jpg 3300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Extent-Polygon-300x232.jpg 300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Extent-Polygon-768x593.jpg 768w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Extent-Polygon-1024x791.jpg 1024w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Extent-Polygon-624x482.jpg 624w\" sizes=\"auto, (max-width: 3300px) 100vw, 3300px\" \/><\/p>\n<p><strong>Fig 3.\u00a0<\/strong>Map of the extent polygon utilized for data acquisition and clipping.<\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 To create a layer encompassing human activity and development in BC, we went into the Other Land Cover layer and used the \u2018Select by Attributes\u2019 function to extract the developed land, perennial, and annual agriculture land use classes. This layer was used alongside the \u2018Other Land Use\u2019 layers acquired from DataBC in our MCE (Fig 4; Fig 5).<\/span><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 We then projected all our layers onto the NAD 1983 BC Environmental Albers coordinate system using the \u2018Project\u2019 tool, and at the same time we set the processing environment to match the cell size of the DEM we were using, in order to prevent any issues when converting our polygons to raster. <\/span><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 Until now almost our layers were in polygon or polyline form, so we then converted our polygon layers to raster, which was necessary in order to utilize them within a Multi-Criteria Analysis. <\/span><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 Once we projected our polygon layers, we then converted them to raster using the \u2018Feature Class to Raster\u2019 tool. Afterwards, we used the \u2018Raster Reclassify\u2019 tool to alter the values in the rasters so that they would represent the presence or absence of a feature on the map. We used the tool to set the cells where an objectID was present to a value of 1 and those which did not have an objectID for that layer to a value of 0. The resulting layer would have cells with only values of 1s and 0s for the entire extent of the map (Fig 4)<\/span><span style=\"font-weight: 400;\">\u00a0<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-57\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Land-use-1.jpg\" alt=\"\" width=\"3300\" height=\"2550\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Land-use-1.jpg 3300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Land-use-1-300x232.jpg 300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Land-use-1-768x593.jpg 768w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Land-use-1-1024x791.jpg 1024w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Land-use-1-624x482.jpg 624w\" sizes=\"auto, (max-width: 3300px) 100vw, 3300px\" \/><\/p>\n<p><strong>Fig 4.\u00a0<\/strong>Map of final land-use categories of our B.C. extent, converted to raster format.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-58\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Roadways-Transmission-Lines-Waterbodies-1.jpg\" alt=\"\" width=\"3300\" height=\"2550\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Roadways-Transmission-Lines-Waterbodies-1.jpg 3300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Roadways-Transmission-Lines-Waterbodies-1-300x232.jpg 300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Roadways-Transmission-Lines-Waterbodies-1-768x593.jpg 768w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Roadways-Transmission-Lines-Waterbodies-1-1024x791.jpg 1024w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/Roadways-Transmission-Lines-Waterbodies-1-624x482.jpg 624w\" sizes=\"auto, (max-width: 3300px) 100vw, 3300px\" \/><\/p>\n<p><strong>Fig 5.\u00a0<\/strong>Map of major roadways, transmission lines, and water bodies<\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 We utilized the \u2018Slope Calculator\u2019 tool to give us the slope from the DEM, and then used the \u2018Map Algebra\u2019 tool to score these slopes with values ranging between 1 and 0; with 1 representing the shallowest slope, and 0 representing the steepest slopes.<\/span><\/p>\n<p><strong>AHP Weighting<\/strong><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 To complete our multi-criteria analysis, we used 123AHP to help determine the weightings to be given to different land features in our study area (Fig 6). We based our decisions on how passable a certain feature is for bears to travel through or live in. Our ranking places lakes and urban\/agricultural land as being most unfavorable for bears, with grasslands and forests being the most favourable. <\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-59\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/123AHPtake2.png\" alt=\"\" width=\"871\" height=\"860\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/123AHPtake2.png 871w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/123AHPtake2-300x296.png 300w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/123AHPtake2-768x758.png 768w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/123AHPtake2-624x616.png 624w\" sizes=\"auto, (max-width: 871px) 100vw, 871px\" \/><\/p>\n<p><strong>Fig 6.\u00a0<\/strong>123AHP results for our weighting decisions.<\/p>\n<p><strong>Grizzly Subpopulations<\/strong><\/p>\n<p><span style=\"font-weight: 400;\">\u00a0 \u00a0 \u00a0 \u00a0 To mark the areas of currently occupied habitat, potential wildlife corridors, and possible reintroduction sites we chose zones\/corridors in our maps which scored well in our MCE. To highlight these high-scoring areas we used the editor feature to draw polygons that matched the shape of these areas (Fig 7).<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-61\" src=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/GRIZZLYBEARSUBPOP.png\" alt=\"\" width=\"748\" height=\"852\" srcset=\"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/GRIZZLYBEARSUBPOP.png 748w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/GRIZZLYBEARSUBPOP-263x300.png 263w, https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/files\/2017\/12\/GRIZZLYBEARSUBPOP-624x711.png 624w\" sizes=\"auto, (max-width: 748px) 100vw, 748px\" \/><\/p>\n<p><strong>Fig 7.\u00a0<\/strong>Map of combined subpopulation areas and abundance estimates, overlaid on the DEM.<\/p>\n<p>&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Overview \u00a0 \u00a0 \u00a0 \u00a0 To answer our three questions we used an MCE (multi-criteria evaluation) to determine the most favourable areas for grizzly movement and settlement. Using different types of land use and land features, both natural and human-made, as our criteria, and using the analytical hierarchical procedure to decide the weightings for use [&hellip;]<\/p>\n","protected":false},"author":55497,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-7","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/pages\/7","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/users\/55497"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/comments?post=7"}],"version-history":[{"count":6,"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/pages\/7\/revisions"}],"predecessor-version":[{"id":91,"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/pages\/7\/revisions\/91"}],"wp:attachment":[{"href":"https:\/\/blogs.ubc.ca\/bcgrizzlydynamics\/wp-json\/wp\/v2\/media?parent=7"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}