View source for Suitable Locations for the Reintroduction of the Black Tailed Prairie Dog

=='''Disclaimer'''==

This tutorial was created for the [http://www.carleton.ca/calendars/ugrad/current/courses/GEOM/4008.html GEOM 4008] course at [http://www.carleton.ca/ Carleton University]. The main goal of this tutorial is to demonstrate how to use Free and Open Source Software for Geospatial (FOSS4G) tools in obtaining a suitability map. Although a specific example is given here, similar procedures can be used for many other suitability studies. 


=='''Introduction'''==

The black-tailed prairie dog (''Cynomys'' ''ludovicianus'') is a large, burrow-dwelling rodent, closely related to the squirrel. This species, having an average length of 40 cm and weighing between 1.0 and 1.5 kg, is very social, usually traveling and living in colonies ([[#References|Parks Canada, 2009]]). Whereas the distribution of their habitat once spanned the Great Plains from Mexico all the way to Canada, black-tailed prairie dogs (BTPDs) have now been classified as a species of "Special Concern" by the Committee on the Status of Endangered Wildlife in Canada ([[#References|COSEWIC, 2010]]). 


The habitat of this prairie dog consists of areas with open grasslands with sparse vegetation. This includes short, mixed and tall grassland types. As well, areas that are dry and flat (or gently sloping) are usually characteristics of BTPD habitats. It is also important that the habitats have easy access to a water body and that they be located away from urbanized areas.


In the late 1800's, the amount of North American prairie dogs exceeded 5 billion individuals, spanning an area of approximately 283 million ha. ([[#References|Hof et al., 2002]]). Now however, estimates suggest that the range of the black-tailed prairie dog has decreased to only 2% of their former distribution ([[#References|Calgary Zoo, 2010]]). As well, it is estimated that there has been a 98-99% reduction in these species across the United States ([[#References|Hof et al., 2002]]). The primary causes of this are thought to be due to habitat loss from urbanization and agriculture expansion, habitat fragmentation, hunting activities as well as susceptibility to disease ([[#References|COSEWIC, 2010]]). 


This reduction in population of the BTPDs is cause for concern because this species is considered a keystone species. Many other species are dependant on them for food, burrows, and for their effects on plant communities ([[#References|Hof et al., 2002]]). It is for these reasons that the goal of this project is to outline the most suitable locations for the reintroduction of the black-tailed prairie dog in Canada. 



=='''Goal'''==

The objective of this study is to outline suitable locations for the re-introduction of the black-tailed prairie dog (BTPD) in Canada using Quantum GIS (qGIS) and Geographic Resources Analysis Support System (GRASS) software. 


Given the habitat requirements for the BTPDs, a suitable site for their reintroduction will be located in an area (1) with slight to no slope; (2) of little precipitation; (3) of short, mixed or tall grassland; (4) in close proximity to a water body and (5) away from urbanized areas.&nbsp; Data layers pertaining to the above mentioned characteristics will be obtained, manipulated and then combined to depict the areas that would be most suitable for further reintroduction efforts. The main functions that will be used throughout this tutorial are creating [[#Rivers|buffers]], [[#Slope|rescaling]] and [[#Urbanized Areas|reclassifying]] data values and [[#Suitability|weighting and combining]] the layers. 



=='''Data'''==
{| style="border-spacing:0;"
| style="border-top:0.035cm solid #808080;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| <center>'''Layer'''</center>
| style="border-top:0.035cm solid #808080;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| <center>'''Type'''</center>
| style="border-top:0.035cm solid #808080;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| <center>'''Source'''</center>
| style="border:0.035cm solid #808080;padding:0.049cm;"| <center>'''Projection'''</center>

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Provincial Boundaries of Canada
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| polygon
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Canadian Wildlife Service, Environment Canada - Prairie and Northern Region
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| NAD83 Geographic

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Digital Elevation Model (DEM)
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| grid
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Natural Resources Canada

Shuttle Radar Topography Mission (SRTM)
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| NAD83 UTM 12N

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Precipitation

&nbsp;
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| point
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Agriculture and AgriFood Canda National Agriclimate Information Service (NAIS)

&nbsp;
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| North American Lambert Conformal Conic

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Vegetation

(Grass/Non-grass)
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| grid
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Environment Canada

PFRA

U.S. Geological Survey

Montana State Library - Natural Resource Information System
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| NAD83 UTM 12N

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Rivers
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| polygon
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Canadian Wildlife Service, Environment Canada - Prairie and Northern Region
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| NAD83 UTM 12N

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Cities
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| polygon
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Canadian Wildlife Service, Environment Canada - Prairie and Northern Region
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| NAD83 UTM 12N

|-
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Grasslands National Park proposed boundary
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| polygon
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:none;padding:0.049cm;"| Parks Canada - Grasslands National Park
| style="border-top:none;border-bottom:0.035cm solid #808080;border-left:0.035cm solid #808080;border-right:0.035cm solid #808080;padding:0.049cm;"| NAD83 UTM 12N

|}
&nbsp;

=='''Software'''==

The software used throughout this tutorial is Quantum GIS (qGIS). This free and open source software can be downloaded [http://www.qgis.org/ here] and can run on Windows, Mac and Linux operating systems. Make sure that, along with the qGIS program, the necessary GRASS plug-in as well as the GRASS software itself are also downloaded (GRASS is included by default in many but not all qGIS distributions for Windows; on other operating systems GRASS often needs to be installed separately).

=='''Methods'''==

==='''Creating a new mapset'''===

First, load the DEM layer into qGIS. This is the base layer that contains the projection information and represents the region for analysis that will be used throughout this tutorial. Next, to create a new mapset, the directory for the database needs to be set. The GRASS location - where the collection of maps for a particular project will be stored - will then need to be created. Next, the projection will need to be specified. For this project, the Universal Transverse Mercator (UTM) NAD 83 zone 12N will be used. The EPSG code for this projection is 26912. The default GRASS Region should be set to the current qGIS extent (that of the DEM). Finally, enter the name of the new mapset (the default name is PERMANENT). 


Once this newly created mapset has been opened, the data analysis can begin. 



==='''Manipulating the Data Layers:'''===

It should be noted that for each command used, the general code is displayed followed by the example code used specifically for this tutorial. In addition, since all of the commands described here are geared for the specific example given within this tutorial, further information and additional conditions for each function can be found in the [http://grass.osgeo.org/gdp/html_grass64/index.html GRASS manual] or in the book by  [[#References| Neteler and Mitasova]].


===='''Slope'''====


A suitable site for the reintroduction of the BTPD is in an area of flat terrain (i.e. an area having very little to no slope). To address this, the DEM layer is used. This raster layer is loaded into GRASS using the r.in.gdal command. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''r.in.gdal '''<nowiki>in=<layer name> out=<layer name></nowiki>''&nbsp;&nbsp;</center>

<center>'''''r.in.gdal''''' ''in=dem_original out=DEM''&nbsp;&nbsp;
|}
</center>


This command allows many different formats to be imported including, but not limited to, ASCII Grids, ERMapper Compressed Wavelets (.ecw) TIFF, JPEG, and Arc/Info Binary Grid Format (.adf). 


Once this DEM layer has been added to the GRASS layer, it will need to be converted to obtain a slope value. The r.slope.aspect command is used to accomplish this. This command uses raster elevation values to create a raster map layer of slope, aspect, curvatures and/or partial derivatives. For this example, we are creating a layer that depicts slope values as a percentage. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''r.slope.aspect''''' ''elevation<nowiki>=<</nowiki>layer name> slope<nowiki>=<</nowiki>layer name> format<nowiki>=<</nowiki>type>''&nbsp;&nbsp;</center>

<center>'''''r.slope.aspect''''' ''elevation=DEM slope=Slope_Percent format=percent''
|}
</center>


Using the r.info command, the minimum and maximum values, along with other basic information about a user-specified raster map layer is obtained. Because the values of slope range from a minimum of 0% to a maximum of approximately 89%, the slope layer needs to be rescaled to range between values of 0 and 10. To do this, the r.rescale command is used. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''r.rescale '''<nowiki>input=<layer name> from=<min,max> output=<layer name> to=<min,max></nowiki>''&nbsp;&nbsp;</center>

<center>'''''r.rescale '''input=Slope_Percent from=0,89 output=Slope_Rescaled to=0,10''
|}
</center>


==== '''Precipitation '''====

Next, the precipitation restriction is addressed. Because the data is given in a point layer vector format, an interpolation method is needed. First however, the projection for this data layer needs to be changed to match that of the DEM. To do this, the utility programs that come with the GDAL/OGR libraries are used. As part of the built-in toolkit, the ogr2ogr command converts data between file formats. For more information on these libraries, please [http://www.gdal.org/ogr/index.html click here]. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;&nbsp;&nbsp;'''''ogr2ogr''''' ''<nowiki><input layer> -t_srs ‚"<projection code>" <output layer></nowiki>''</center>

<center>&nbsp;&nbsp;'''''ogr2ogr''''' ''precipitation_original -t_srs "EPSG:26912" precipitation.shp''&nbsp;&nbsp;
|}
</center>


Now that the precipitation layer has the correct projection, this vector layer is loaded into GRASS using the v.in.ogr command. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''v.in.ogr '''<nowiki>in=<layer name> out=<layer name></nowiki>''&nbsp;&nbsp;</center>

<center>'''''v.in.ogr''''' ''in=precipitation.shp out=Precip''
|}
</center>


To create an interpolated raster surface from the vector points, the v.surf.idw command is used. This command requires the user to specify the column within which the values are to be obtained and the number of nearest data points to use in the calculation. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''v.surf.idw '''<nowiki>input=<layer name> =point layer=1 column=<column name> npoints=<integer> output=<layer name></nowiki>''&nbsp;&nbsp;</center>

<center>'''''v.surf.idw '''input=Precip =point layer=1 column=P npoints=12 output=IDW_Precip''
|}</center>


While this tutorial uses an Inverse Distance Weighting (IDW) interpolation method with 12 sample points, altering the number of points used or using a different interpolation method (such as kriging or splines) can also be used. 


Finally, the IDW_Precipitation layer needs to be rescaled so that its values range from 0 to 10. The same procedure conducted for the slope layer is used. 



==== '''Vegetation'''====

Habitats for the BTPD tend to be located in areas composed of short, mixed or tall grasslands. As the vegetation data for this tutorial is given as presence/absence data (with grasslands having a value of 1 and non-grass areas having a value of 0), this data needs to be reclassified. Once the vegetation layer is added to GRASS using the r.in.gdal command, a new text file needs to be created. This text file contains the reclassification rules. For the purposes of this tutorial, the reclassification rules are as follows:


<center>0 = 10</center>

<center>1 = 0</center>


It should, however, be noted that had the data been continuous instead of binary, the reclassification rules could have been:


<center>0.0 thru 0.5 = 10</center>

<center>0.6 thru 1.0 = 0</center>


Next the r.reclass command is used. This command creates a new raster layer whose category values are based on the reclassification rules outlined in a specified file. 



<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;&nbsp;&nbsp;'''''r.reclass '''<nowiki>input=<layer name> rules=<rules file> output=<layer name></nowiki>''</center>

<center>&nbsp;&nbsp;'''''r.reclass '''input=Vegetation rules=Vegetation_rules.txt output=Vegetation_Reclass''&nbsp;&nbsp;

|}</center>


===='''Rivers'''====

BTPD habitats should be located in close proximity to a water body.&nbsp; Therefore, this vector layer depicting rivers needs to first be loaded into GRASS, converted into a raster, have a buffer created around it and then be rescaled to match the scale of the other layers. 


Because the river layer is a vector layer, the v.in.ogr command is used. Then, the v.to.rast command is used to rasterize the vector map. This function requires a column name to be specified. This column must contain integer values, as these values will become the attributes for the raster map. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''v.to.rast '''<nowiki>int=<layer name> out=<layer name> col=<column name></nowiki>''&nbsp;&nbsp;</center>

<center>'''''v.to.rast '''int=Rivers out=Rivers_Raster col=Column_1''
|}</center>


Next, the buffers around the rivers is created using the r.buffer command. This function creates a new raster map layer showing the distance/proximity/ buffer zones around all cell values that are non-NULL. 
The distances are calculated and drawn around the initial category value based on the specified distance values. For the purposes of this tutorial, the buffer distances are designated as 2 km, 4 km, 6 km, 8 km and 10 km away from the rivers. It should be noted that the units can be specified as either meters, feet, miles, kilometers or nautmiles. The actual command for this function is as follows:


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''r.buffer '''<nowiki>input=<layer name> output=<layer name> distances =<buffer distance 1, buffer distance 2, ... , buffer distance n> units=<units></nowiki>''&nbsp;&nbsp;</center>

<center>'''''r.buffer '''input=Rivers_Raster output=Rivers_Buffer distances =2,4,6,8,10 units=kilometers''
|}
</center>


Lastly, the Rivers_Buffer layer is rescaled to have its values range from 0 to 10 using the same method used for the slope and precipitation layers. 



===='''Urbanized Areas'''====

The last layer to be addressed is the cities layer. This layer requires the most processing. First, this vector layer needs to be loaded into GRASS using v.in.ogr. Next, it is converted into a raster layer using the same command used to convert the Rivers layer. When using the v.to.rast command, be sure to specify which column the new raster layer attributes should be based on. Once the raster layer is created, buffer zones are to be created. This time, the buffer zones are to be as follows: 0 km to 10 km, 10 km to 20 km, 20 km to 50 km and 50 km to 800 km.


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;'''''r.buffer '''<nowiki>input=<layer name> output=<layer name> distances =<buffer distance 1, buffer distance 2, ... , buffer distance n> units=<units></nowiki>''&nbsp;&nbsp;</center>

<center>'''''r.buffer '''input=Cities_Raster output=Cities_Buffer distances =10,30,50,800 units=kilometers''
|}
</center>


Since the optimal location for the reintroduction of the BTPDs would be away from urbanized areas, this City_Buffer layer needs to be reclassified so that values would be lower at locations further away from cities. Using r.info Cities_Buffer, the minimum and maximum values for this layer are obtained (they range from 1 to 5). This information is then used to create the reclassification rules seen below:


<center>5 = 0</center>

<center>4 = 3</center>

<center>3 = 6</center>

<center>2 = 8 </center>

<center>1 = 10</center>


The r.reclass command is then used. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;&nbsp;'''''r.reclass '''<nowiki>input=<layer name> rules=<rules file> output=<layer name></nowiki>''&nbsp;&nbsp;</center>

<center>&nbsp;&nbsp;'''''r.reclass '''input=Cities_Buffer rules=Cities_rules.txt output=Cities_Reclass''&nbsp;&nbsp;
|}
</center>



===='''Suitability'''====

The last step in this tutorial involves weighting and combining all of the above layers to create a single raster layer that displays the most suitable locations for the reintroduction of the BTPD. First, weights are assigned to each layer. These weights are as follows:


<center>'''Vegetation = 0.35'''</center>

<center>'''Slope = 0.25'''</center>

<center>'''Cities = 0.15'''</center>

<center>'''Precipitation = 0.15'''</center>

<center>'''Rivers = 0.1'''</center>


Next the r.mapcalc command is used. This function allows simple mathematical equations to be performed on up to 6 existing raster maps. 


<center>
{|style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; border-color: #000; padding: 0"
|-
|style="border-style: solid; border-width: 0 1px 1px 0"|
&nbsp;&nbsp;&nbsp;&nbsp;'''''r.mapcalculator''''' ''<nowiki>amap=<layer name> bmap=<layer name> cmap=<layer name> dmap=<layer name> emap=<layer name> fmap=<layer name> formula=<formula> outfile=<layer name></nowiki>''&nbsp;&nbsp;</center>

<center>&nbsp;&nbsp;'''''r.mapcalculator''''' ''amap=Vegetation_Reclass bmap=Slope_Rescaled cmap=Cities_Reclass dmap=Precip_Rescaled emap=Rivers_Rescaled formula=(0.35*A)+(0.25*B)+(0.15*C)+(0.15*D)+(0.1*E) outfile=Suitability_weights''&nbsp;&nbsp;
|}
</center>



The resulting raster map shows the optimal locations for the reintroduction of the black-tailed prairie dogs. The lower the value, the more suitable the location. While the results of this particular example are depicted in Figure 2, the results will vary greatly depending on the layers used and on the weighting scheme implemented. 

<center>
[[File:Figure_2.jpg|border|''Figure 2: Resulting suitability map for the reintroduction of the black-tailed prairie dog. Values range from high suitability (green) to low suitability (pink).''  ]]
</center>
<center>''Figure 2: Resulting suitability map for the reintroduction of the black-tailed prairie dog. Values range from high suitability (yellow) to low suitability (pink).''
</center>



=='''References'''==

<nowiki>Calgary Zoo (2010) Black-tailed Prairie Dog [Internet] Available from: </nowiki>[http://www.calgaryzoo.org/index.php?option=com_content&task=view&id=506&Itemid=604 http://www.calgaryzoo.org/index.php?option=com_content&task=view&id=506&Itemid=604] Accessed 2010 Nov. 26


<nowiki>COSEWIC (2010) Wildlife Species Search [Internet] Available from: </nowiki>[http://www.cosewic.gc.ca/eng/sct1/searchdetail_e.cfm http://www.cosewic.gc.ca/eng/sct1/searchdetail_e.cfm] Accessed 2010 Oct. 26


<nowiki>GRASS (2010) Welcome to GRASS GIS [Internet] Available from: </nowiki>[http://grass.osgeo.org/ http://grass.osgeo.org/] Accessed 2010 Nov. 27


''Hof, J., Bevers, M., Uresk, D. W., Schenbeck, G. L. ''(2002) Optimizing habitat location for black-tailed prairie dogs in southwestern South Dakota; ‚Ecological Modelling, Vol. 147, pp 11-21.


Neteler, M. and Mitasova, H. (2002). Open Source GIS: A GRASS GIS Approach. Kluwer Academic Publishers (Norwell, MA). 


<nowiki>Parks Canada (2009) Grasslands National Park of Canada‚ The Black-tailed Prairie Dog [Internet] Available from: </nowiki>[http://www.pc.gc.ca/eng/pn-np/sk/grasslands/edu/edu1/a.aspx http://www.pc.gc.ca/eng/pn-np/sk/grasslands/edu/edu1/a.aspx] Accessed 2010 Oct. 26

[[Category:Tutorials]]