ArcGIS Role in Maxent Modeling Christopher Woods Christopher.R.Woods @leidos.com Carpinteria, CA
Modeling Remember that all models are wrong; the practical question is how wrong do they have to be to not be useful Statistician, George E. P. Box Empirical Model-Building and Response Surfaces (1987)
What is Maxent Maxent stands for Maximum Entropy Type of inductive modeling that uses presence only occurrence data along with environmental variables to derive a probability of occurrence or potential suitable habitat. Software can be downloaded from http://www.cs.princeton.edu/~schapire/maxent/
About Maxent Runs in a separate application Inputs: CSV file of occurrences ASCII Grid of Environmental Layers Outputs: ASCII Grid of Potential Suitable Habitat
Why use Maxent? Limited Range or Habitat data exists for species of interest Occurrence records available Locates Potential Suitable Habitat Subject Matter Experts needed to ground truth results to highlight areas needed for more focused study to determine presence or absence
Environmental Layers Used in Maxent PRISM Precipitation Example of environmental raster layers used in Maxent Layers are stacked upon each other within the model PRISM Temperature Ruggedness Aspect Vegetation Distance to Water Soils
Occurrence Data Used in Maxent Occurrence data (point locations) are added to the model Stored in a CSV file containing X,Y positions in the same coordinate system as environmental layers Species Observation
Occurrence Data Used in Maxent Occurrence data (point locations) are added to the model Stored in a CSV file containing X,Y positions in same coordinate system as environmental layers Maxent uses the point occurrence data to look for relationships within the stack of environmental layers Aspect Species Observation Ruggedness Prism Temperature Prism Precipitation
ArcGIS Role in Maxent Modeling
ArcGIS Role Creating Snap Grid for Study Area Create Environmental Layers Align Layers to snap grid using the same projection, pixel size and grid extent Create CSV of Occurrence Data Populate Samples with Data (optional) Convert Environmental Data to ASCII Convert Results from ASCII to Imagine/Tiff/Grid Symbolize Results on Map
Using ArcGIS to Create Environmental Layers
Using ArcGIS to Create Environmental Layers Soil, Geology, Vegetation
LANDFIRE Existing Vegetation, Geology, Soils LANDFIRE Existing Vegetation type, currently raster dataset (30m cell size resolution) Extract using Snap Grid US Geology, currently a vector dataset Convert to Raster and limit extent using the snap grid STATSGO soils, currently a vector SSURGO is better resolution but gaps can exist when using it for a large area (such as a state) Convert to raster and limit the extent using the snap grid
14
15
16
Using ArcGIS to Create Environmental Layers Elevation and Slope
Elevation and Slope Extract Elevation from the National Elevation Dataset or other gridded dataset to Snap Grid Percent Slope is calculated from the Elevation Dataset using the standard surface toolbox within the Spatial Analyst toolset. Specify Snap Grid in Environment
19
20
Using ArcGIS to Create Environmental Layers Aspect
Aspect Aspect is calculated from the Elevation Dataset using the standard surface toolbox within the Spatial Analyst toolset. Specify Snap Grid in Environment Settings Further steps are required to extract North/South Aspect and East/West Aspect
Northness and Eastness Variables vs. Aspect in Degrees 1 N E 0.8 0.6 Directionality Score 0.4 0.2 0-0.2-0.4 Northness Eastness -0.6-0.8-1 S W 0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 Aspect (Degrees) Eastness = sin ((aspect in degrees * PI)/180) Northness = cos ((aspect in degrees * PI)/180) Zar, J. H. 1999. Biostatistical Analysis. Prentice Hall, New Jersey. 663 pp.
Toolbox to generate the two Aspect Layers 24
25
26
Using ArcGIS to Create Environmental Layers Terrain Ruggedness
Terrain Ruggedness Calculated from Elevation Python tool for calculating ruggedness based on paper by Sappington et al. 2007 Sappington, J.M., K.M. Longshore, and D.B. Thomson. 2007. Quantifiying Landscape Ruggedness for Animal Habitat Anaysis: A case Study Using Bighorn Sheep in the Mojave Desert. Journal of Wildlife Management. 71(5): 1419-1426
Terrain Ruggedness calculated from Elevation 29
30
Using ArcGIS to Create Environmental Layers Climate
Climate Datasets PRISM (Oregon State University) Maximum Temperature Minimum Temperature Annual Precipitation 1981-2010 Normals 800m cell size Monthly 4km data Can also use a specific month or bins of months for maximum and minimum temperature Other climate datasets available based on region such as Climate Impacts Group (CIG), National Centers for Environmental Prediction (NCEP), etc.
33
34
35
Using ArcGIS to Create Environmental Layers Distance to Water
Distance to Water Perennial Waters (Water Bodies and Streams) converted to rasters and mosaiced to create a single layer Euclidean Distance is used to measure distance from perennial water throughout the raster Data Sources: State Datasets, National Hydrographic Dataset, etc.
Distance to Water Model 38
Distance to Water
Using ArcGIS to Create Environmental Layers Solar Radiation
Solar Radiation Index The Solar Radiation Index (SRI) was calculated for the Equinox, Winter Solstice and Summer Solstice following the method described by Keating et al in their paper A Simple Solar Radiation Index for Wildlife Habitat Studies. Journal of Wildlife Management 71 (4); 1344 1348; 2007 SRI = (mean sun-earth distance / daily sun-earth distance)2* {[sin(latitude)cos(slope) cos(latitude)sin(slope)cos(aspect)]sin(declination) + [cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect)]cos(declination)cos(hour angle) + cos(declination)sin(slope)sin(aspect)sin(hour angle)}
Solar Radiation Index Mean Sun-Earth Distance = 1 astronomical unit (au) Daily Sun-Earth Distance: Equinox = 1 au, Winter Solstice = 0.983747 au, Summer Solstice = 1.017204 au Hour Angle: All calculations are for the hour surrounding Solar Noon. The Hour Angle during this time = 0. Latitude: $$Ymap command used in Raster Calculator to create a raster with the latitude of each cell center as the cell value. Converted to radians. Slope: Slope raster in degrees converted to radians. Aspect: 180 Aspect raster (converts the 0-360 degree data set to a -180-180 degree data set), converted to radians. Declination: Equinox = 0, Winter Solstice = -23.45 degrees, Summer Solstice = 23.45 degrees, converted to radians.
Solar Radiation Index Equinox: SRI = cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect) Winter Solstice: SRI = (1/0.983747)2* {[sin(latitude)cos(slope) cos(latitude)sin(slope)cos(aspect)] * (- 0.39794863130761038954479576746719) + [cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect)] * (0.91740769935748826377361893702175) Summer Solstice: SRI = (1/1.017204)2* {[sin(latitude)cos(slope) cos(latitude)sin(slope)cos(aspect)] * (0.39794863130761038954479576746719) + [cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect)] * (0.91740769935748826377361893702175)
Using ArcGIS to Create Environmental Layers Conversion to ASCII Raster
Conversion to Raster All Environmental Rasters need to be in ASCII raster format for use in the model. All ASCII rasters need to have the same: Extent Cell size Projection Conversion to ASCII can be done from ArcToolbox Conversion toolbox Helpful to specify the snap grid using extent, cells size, snap raster, etc in the Environment settings just to be safe
Using ArcGIS to Create Environmental Layers Extract Occurrence Data
Occurrence Data Accuracy Timeliness (old data, spatial recording method) Precision (GPS, trapping records, etc.) Relevance road kill data or observations taken from a road (don t necessarily want to have roads modeled as suitable habitat
Extract Occurrence Data to CSV Extract occurrence data in the same projection and extent as the raster environmental layers. Needs to be in a specific structure Name (of what is being modeled) X (Easting) Y (Northing) Additional Columns with can be added for Samples with Data
Samples with Data Used when the analysis cell size being used within Maxent is greater than the native data. Using 90m cell resolution when your elevation or vegetation datasets originally were in 30m resolution The original value is extracted from the native dataset and stored with the occurrence data in a column with the CSV file. This allows for more accuracy and prevents Maxent from extracting the value from a more coarse dataset when the model is run Samples with data are stored in the CSV file Column name should match the environmental raster layer name Must be in the same order that the environmental layers are stacked within the Maxent model
Samples with Data Tools such as Hawth s Tool s Point Intersect Tool or Geospatial Modeling Environment s isectpntrst function Allows for the processing many rasters at once. http://www.spatialecology.com/htools/isect.php http://www.spatialecology.com/gme/isectpntrst.htm
Using ArcGIS to Create Environmental Layers Running Maxent
52
53
Maxent Modeled Output The results of a Maxent model run contain values that range from 0 to 1. Values closer to 0 have less potential suitable habitat while values closer to 1 have higher potential suitability. A threshold is defined that determines what value between 0 to 1 is the cutoff for what is potentially suitable habitat and what is not.
55
56
57
58
Using ArcGIS to Create Environmental Layers Pitfalls
Suggestions for Minimizing Headaches Start first by making a snap raster and use it in the environment settings. Make sure all the rasters have the same extent, projection, alignment Error messages will occur if occurrence data doesn t overlap study area. Maxent will still run but will pop up an error message for each one. Use the ArcToolbox Conversion Tool to convert ASCII grid to Raster. Make sure to choose Float instead of the default Integer Careful of Pyramids. Maxent data can be very pixelated so pyramiding data may give a map a different view point. 60
Maxent Topics Not Able to Cover Due to time constraints and being a GIS conference I focused on the role of ArcGIS within running Maxent modeling. I purposely glossed over the actual running of the model, validation results, replicates, defining thresholds, etc.
Acknowledgments Scott Story: Montana Fish Wildlife and Parks Anthony Titolo: BLM Wildlife Habitat Spatial Analysis Lab Lara Juluisson: BLM Wildlife Habitat Spatial Analysis Lab