{"id":22,"date":"2021-12-11T17:48:21","date_gmt":"2021-12-12T00:48:21","guid":{"rendered":"https:\/\/blogs.ubc.ca\/glacierforecast\/?page_id=22"},"modified":"2021-12-15T14:37:21","modified_gmt":"2021-12-15T21:37:21","slug":"22-2","status":"publish","type":"page","link":"https:\/\/blogs.ubc.ca\/glacierforecast\/22-2\/","title":{"rendered":"Methodology"},"content":{"rendered":"<p>Everything in this model was created using ArcGis Pro software.<\/p>\n<h3>Data Collection + Organization<\/h3>\n<p>The first step of making this model is to collect the data layers from their sources. These sources can be found in the\u00a0<a href=\"https:\/\/blogs.ubc.ca\/glacierforecast\/sample-page\/\">References<\/a>\u00a0page. This specific model only shows the results for the southern half of British Columbia, but all the initial data layers represent all of British Columbia. I set my map processing extent to a polygon of southern British Columbia and my projection-based coordinate system is NAD 83 BC Albers Environmental projection. I made sure to transform any data layers into this projection system, this can easily be done by setting the map projection to BC Albers, so any new data will automatically be transformed to this projection system.<\/p>\n<p>To collect the DEM layer, I plotted rough coordinates of British Columbia into the EarthExplorer program, it does not need to be the exact shape of British Columbia shown in Figure 1 below.<\/p>\n<figure id=\"attachment_30\" aria-describedby=\"caption-attachment-30\" style=\"width: 1527px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/blogs.ubc.ca\/glacierforecast\/22-2\/elevation\/\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-30 size-full\" src=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/elevation.png\" alt=\"\" width=\"1527\" height=\"1177\" srcset=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/elevation.png 1527w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/elevation-300x231.png 300w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/elevation-1024x789.png 1024w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/elevation-768x592.png 768w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/elevation-624x481.png 624w\" sizes=\"auto, (max-width: 1527px) 100vw, 1527px\" \/><\/a><figcaption id=\"caption-attachment-30\" class=\"wp-caption-text\">Figure 1: Elevation map of British Columbia. Data source in References<\/figcaption><\/figure>\n<p>Then, using the slope and aspect geoprocessing tools, I extracted the slope and aspect data layers of British Columbia shown in Figures 2 and 3 below.<\/p>\n<p>&nbsp;<\/p>\n<figure id=\"attachment_29\" aria-describedby=\"caption-attachment-29\" style=\"width: 1530px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/blogs.ubc.ca\/glacierforecast\/22-2\/slope\/\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-29 size-full\" src=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/slope.png\" alt=\"\" width=\"1530\" height=\"1181\" srcset=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/slope.png 1530w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/slope-300x232.png 300w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/slope-1024x790.png 1024w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/slope-768x593.png 768w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/slope-624x482.png 624w\" sizes=\"auto, (max-width: 1530px) 100vw, 1530px\" \/><\/a><figcaption id=\"caption-attachment-29\" class=\"wp-caption-text\">Figure 2: Slope map of British Columbia. Data source in References<\/figcaption><\/figure>\n<figure id=\"attachment_23\" aria-describedby=\"caption-attachment-23\" style=\"width: 1535px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/blogs.ubc.ca\/glacierforecast\/22-2\/aspect\/\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-23 size-full\" src=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/aspect.png\" alt=\"\" width=\"1535\" height=\"1182\" srcset=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/aspect.png 1535w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/aspect-300x231.png 300w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/aspect-1024x789.png 1024w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/aspect-768x591.png 768w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/aspect-624x481.png 624w\" sizes=\"auto, (max-width: 1535px) 100vw, 1535px\" \/><\/a><figcaption id=\"caption-attachment-23\" class=\"wp-caption-text\">Figure 3: Aspect map of British Columbia. Data source in References<\/figcaption><\/figure>\n<p>Then, creating a new feature class in the catalog pane, I traced the lines of the BC polygon layer to create the BC coast line layer, which was used to calculate Coast Distance. To do this I used the euclidean distance spatial analyst tool. Coast Distance can be seen in Figure 4 below.<\/p>\n<figure id=\"attachment_26\" aria-describedby=\"caption-attachment-26\" style=\"width: 1529px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/blogs.ubc.ca\/glacierforecast\/22-2\/coast-distan\/\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-26 size-full\" src=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/coast-distan.png\" alt=\"\" width=\"1529\" height=\"1185\" srcset=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/coast-distan.png 1529w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/coast-distan-300x233.png 300w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/coast-distan-1024x794.png 1024w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/coast-distan-768x595.png 768w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/coast-distan-624x484.png 624w\" sizes=\"auto, (max-width: 1529px) 100vw, 1529px\" \/><\/a><figcaption id=\"caption-attachment-26\" class=\"wp-caption-text\">Figure 4: Coast Distance map of British Columbia. Data source in References<\/figcaption><\/figure>\n<p>There is no need to reclassify these 4 data layers (Figures 1, 2, 3 and 4) as they are all going to stay consistent throughout the temporal study of suitability and can be transformed with the Suitability modeller.<\/p>\n<p>Then I needed to access the ClimateBC raster dataset and download the circled data in Figure 5 below.<\/p>\n<figure id=\"attachment_35\" aria-describedby=\"caption-attachment-35\" style=\"width: 1279px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-35 size-full\" src=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/B4F797E7-5278-4B8B-AC00-38E6BCA74F30.jpeg\" alt=\"\" width=\"1279\" height=\"1335\" srcset=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/B4F797E7-5278-4B8B-AC00-38E6BCA74F30.jpeg 1279w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/B4F797E7-5278-4B8B-AC00-38E6BCA74F30-287x300.jpeg 287w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/B4F797E7-5278-4B8B-AC00-38E6BCA74F30-981x1024.jpeg 981w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/B4F797E7-5278-4B8B-AC00-38E6BCA74F30-768x802.jpeg 768w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/B4F797E7-5278-4B8B-AC00-38E6BCA74F30-624x651.jpeg 624w\" sizes=\"auto, (max-width: 1279px) 100vw, 1279px\" \/><figcaption id=\"caption-attachment-35\" class=\"wp-caption-text\">Figure 5: ClimateBC Raster Datasets Available for download<\/figcaption><\/figure>\n<p>Each ZIP file contains 56 different seasonal variables, the two relevant seasonal variables to this study are PAS_wt (winter precipitation as snow mm) and Tave_sm (summer mean temperature C). Then, I needed to delete all the other variables and make sure to label each layer according to the year. These layers don\u2019t need to be reclassified either because the suitability modeller can reclassify and transform the data, standardizing it temporally.<\/p>\n<h3>Analytical Hierarchy Process (AHP) Weighting<\/h3>\n<p>An AHP generator was used to determine the proper weights for each variable and than a sensitivity analysis for one of the normal periods (1961-1990) was conducted against equal weights.<\/p>\n<table>\n<tbody>\n<tr>\n<td width=\"231\"><strong>Factor<\/strong><\/td>\n<td width=\"231\"><strong>Weight<\/strong><\/td>\n<\/tr>\n<tr>\n<td width=\"231\">Summer Mean Temp<\/td>\n<td width=\"231\">31.80<\/td>\n<\/tr>\n<tr>\n<td width=\"231\">Winter Precipitation as Snow<\/td>\n<td width=\"231\">31.80<\/td>\n<\/tr>\n<tr>\n<td width=\"231\">Elevation<\/td>\n<td width=\"231\">16.50<\/td>\n<\/tr>\n<tr>\n<td width=\"231\">Slope<\/td>\n<td width=\"231\">8.50<\/td>\n<\/tr>\n<tr>\n<td width=\"231\">Coast Distance<\/td>\n<td width=\"231\">4.70<\/td>\n<\/tr>\n<tr>\n<td width=\"231\">Aspect<\/td>\n<td width=\"231\">6.70<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<h3>Suitability Modeller<\/h3>\n<p>Five separate suitability models were created for each 30-year normal period, designated from the past to the future, A to E. For the A (1961-1990) model the transformations can be seen below:<\/p>\n<table>\n<tbody>\n<tr>\n<td width=\"208\"><strong>Data Layer<\/strong><\/td>\n<td width=\"208\"><strong>Transformation<\/strong><\/td>\n<td width=\"208\"><strong>Notes<\/strong><\/td>\n<\/tr>\n<tr>\n<td width=\"208\">Elevation<\/td>\n<td width=\"208\">Continuous Transformation: Inverse Mssmall<\/td>\n<td width=\"208\">&lt;0 = 0 suitability<\/p>\n<p>Anything below sea-level has no suitability for Glacier Growth<\/td>\n<\/tr>\n<tr>\n<td width=\"208\">Coast Distance<\/td>\n<td width=\"208\">Continuous Transformation: Mssmall<\/td>\n<td width=\"208\"><\/td>\n<\/tr>\n<tr>\n<td width=\"208\">Slope<\/td>\n<td width=\"208\">Continuous Transformation: Mssmall<\/td>\n<td width=\"208\"><\/td>\n<\/tr>\n<tr>\n<td width=\"208\">Aspect<\/td>\n<td width=\"208\">Continuous Transformation: Inverse Gaussian<\/td>\n<td width=\"208\">This places highest suitability for north facing slopes, and lowest suitability for south facing slopes<\/td>\n<\/tr>\n<tr>\n<td width=\"208\">Summer Mean Temp<\/td>\n<td width=\"208\">Continuous Transformation: Inverse Linear<\/td>\n<td width=\"208\">&gt;200 = 0 suitability<\/p>\n<p>&lt;0 = 100 suitability<\/p>\n<p>If the temp is higher than 20C, max melting can be assumed<\/p>\n<p>If the temp is lower than 0C than there is no melting, therefore max suitability<\/td>\n<\/tr>\n<tr>\n<td width=\"208\">Winter Precipitation as Snow<\/td>\n<td width=\"208\">Continuous Transformation: Linear<\/td>\n<td width=\"208\">&gt;2000mm = 100 Suitability<\/p>\n<p>If the snowfall is higher than 2000 mm, that is significant snowfall and therefore high suitability to glacier growth<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>The reason limits are given to Winter Precip, and Summer Temp layers is to standardize them against the other temporal models. For example: if 1961-1990 min temp is -7C, while 2071-2100 min temp is 3C, if no restrictions were in place than a suitability of 100 would be given to both, even though in one the snow isn\u2019t melting, and in the other melting is a possibility. This therefore standardizes each set of temp and precip datasets against each-other, and makes sure not to skew the suitability models and properly show the progression of glacier growth locations overtime.<\/p>\n<p>Then, I used the suitability modeller five times, each using the same Elevation, Slope, Aspect and Coastal \u00a0Distance layers, while switching the Temp and Precip data layers for the different 30-year normal periods.<\/p>\n<h3>Sensitivity Analysis<\/h3>\n<p>A sensitivity analysis is used to determine uncertainty\/error when using mathematical transformations against data, it is used to determine how the inputs actually affect the results, using two different weights of the same transformation.<\/p>\n<p>My last step was to make a sensitivity analysis. This was done by using the locate tool on the suitability modeller for the AHP weights and Equal weights. Figure 6 shows the flow chart. The reclassify tool in Figure 6 converts all values\/locations to \u201c1\u201d, so that when the two models are added (the plus tool) the locations which overlap are labelled as 2, while the separate values are labelled as 1.<\/p>\n<figure id=\"attachment_25\" aria-describedby=\"caption-attachment-25\" style=\"width: 625px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-25 size-large\" src=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model-1024x279.png\" alt=\"\" width=\"625\" height=\"170\" srcset=\"https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model-1024x279.png 1024w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model-300x82.png 300w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model-768x209.png 768w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model-1536x418.png 1536w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model-624x170.png 624w, https:\/\/blogs.ubc.ca\/glacierforecast\/files\/2021\/12\/Sensitivity-Model.png 1708w\" sizes=\"auto, (max-width: 625px) 100vw, 625px\" \/><figcaption id=\"caption-attachment-25\" class=\"wp-caption-text\">Figure 6: Sensitivity Model Flow<\/figcaption><\/figure>\n<p>&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Everything in this model was created using ArcGis Pro software. Data Collection + Organization The first step of making this model is to collect the data layers from their sources. These sources can be found in the\u00a0References\u00a0page. This specific model only shows the results for the southern half of British Columbia, but all the initial [&hellip;]<\/p>\n","protected":false},"author":89115,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-22","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/pages\/22","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/users\/89115"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/comments?post=22"}],"version-history":[{"count":11,"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/pages\/22\/revisions"}],"predecessor-version":[{"id":77,"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/pages\/22\/revisions\/77"}],"wp:attachment":[{"href":"https:\/\/blogs.ubc.ca\/glacierforecast\/wp-json\/wp\/v2\/media?parent=22"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}