ArcGIS Role in Maxent Modeling

Transcription

1 ArcGIS Role in Maxent Modeling Christopher Woods Carpinteria, CA

2 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)

3 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

4 About Maxent Runs in a separate application Inputs: CSV file of occurrences ASCII Grid of Environmental Layers Outputs: ASCII Grid of Potential Suitable Habitat

5 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

6 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

7 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

8 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

9 ArcGIS Role in Maxent Modeling

10 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

11 Using ArcGIS to Create Environmental Layers

12 Using ArcGIS to Create Environmental Layers Soil, Geology, Vegetation

13 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 14

15 15

16 16

17 Using ArcGIS to Create Environmental Layers Elevation and Slope

18 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 19

20 20

21 Using ArcGIS to Create Environmental Layers Aspect

22 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

23 Northness and Eastness Variables vs. Aspect in Degrees 1 N E Directionality Score Northness Eastness S W Aspect (Degrees) Eastness = sin ((aspect in degrees * PI)/180) Northness = cos ((aspect in degrees * PI)/180) Zar, J. H Biostatistical Analysis. Prentice Hall, New Jersey. 663 pp.

24 Toolbox to generate the two Aspect Layers 24

25 25

26 26

27 Using ArcGIS to Create Environmental Layers Terrain Ruggedness

28 Terrain Ruggedness Calculated from Elevation Python tool for calculating ruggedness based on paper by Sappington et al Sappington, J.M., K.M. Longshore, and D.B. Thomson Quantifiying Landscape Ruggedness for Animal Habitat Anaysis: A case Study Using Bighorn Sheep in the Mojave Desert. Journal of Wildlife Management. 71(5):

29 Terrain Ruggedness calculated from Elevation 29

30 30

31 Using ArcGIS to Create Environmental Layers Climate

32 Climate Datasets PRISM (Oregon State University) Maximum Temperature Minimum Temperature Annual Precipitation 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 33

34 34

35 35

36 Using ArcGIS to Create Environmental Layers Distance to Water

37 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.

38 Distance to Water Model 38

39 Distance to Water

40 Using ArcGIS to Create Environmental Layers Solar Radiation

41 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); ; 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)}

42 Solar Radiation Index Mean Sun-Earth Distance = 1 astronomical unit (au) Daily Sun-Earth Distance: Equinox = 1 au, Winter Solstice = au, Summer Solstice = 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 degree data set to a degree data set), converted to radians. Declination: Equinox = 0, Winter Solstice = degrees, Summer Solstice = degrees, converted to radians.

43 Solar Radiation Index Equinox: SRI = cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect) Winter Solstice: SRI = (1/ )2* {[sin(latitude)cos(slope) cos(latitude)sin(slope)cos(aspect)] * ( ) + [cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect)] * ( ) Summer Solstice: SRI = (1/ )2* {[sin(latitude)cos(slope) cos(latitude)sin(slope)cos(aspect)] * ( ) + [cos(latitude)cos(slope) + sin(latitude)sin(slope)cos(aspect)] * ( )

44 Using ArcGIS to Create Environmental Layers Conversion to ASCII Raster

45 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

46 Using ArcGIS to Create Environmental Layers Extract Occurrence Data

47 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

48 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

49 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

50 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.

51 Using ArcGIS to Create Environmental Layers Running Maxent

52 52

53 53

54 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 55

56 56

57 57

58 58

59 Using ArcGIS to Create Environmental Layers Pitfalls

60 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

61 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.

62 Acknowledgments Scott Story: Montana Fish Wildlife and Parks Anthony Titolo: BLM Wildlife Habitat Spatial Analysis Lab Lara Juluisson: BLM Wildlife Habitat Spatial Analysis Lab