Jason Kreitler
United States Geological Survey, Boise, ID, 83702
E-mail: [email protected]

R script for calculating climate stability:    ClimateStability.r (MD5: 47588dbeb7967456168c39efd6d3883c)

This code uses a convex hull to approximate climate space for a geographical area. In the main text, the code is used to compare networks of protected areas. The code requires raster data describing climate space based on n variables (3 in the example: DFJTmin, JJATmax, CWD) for historic (hst) and future, and a raster of unique cell IDs for planning units within the study area. In the main text this method was used to calculate climate stability for regional networks as a whole (e.g., iNetLU, vNetLU, iNetReg, vNetReg), and also for the networks within individual landscape units. The output is a .csv table with 5 columns: "planning" the ID for the planning unit, "hullv_hst" the historic/present/current hull volume by planning unit, "hullv_fut" the future value according to a GCM scenario, "hulllv_both" the hull created by the union of future and current climate data, and "climstab" the proportion of future climate space that overlaps current climate space (stability if positive), or a unitless measure of departure (negative values).

Any use of this code should reference main manuscript, Heller et al. 2015. Targeting climate diversity in conservation planning to build resilience to climate change. Ecosphere.

Neither the U.S. Government, the Department of the Interior, nor the USGS, nor any of their employees, contractors, or subcontractors, make any warranty, express or implied, nor assume any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, nor represent that its use would not infringe on privately owned rights.