Name: integrating remote sensing with species distribution models; mapping tamarisk invasions using the software for assisted habitat modeling (sahm)
Uploaded: 2016-10-11T15:00:00
Duration: 12 min 26 s
Description: Watch this Scientific Journal Video about Integrating Remote Sensing with Species Distribution Models; Mapping Tamarisk Invasions Using the Software for Assisted Habitat Modeling SAHM at JoVE.com

The authors would like to thank the U.S. Geological Survey, Natural Resource Ecology Laboratory at Colorado State University, Colorado State Forest Service and the Tamarisk Coalition for logistical support, data, use of facilities and expertise. Additionally, we thank Shelly Simmons, Lane Carter, John Moore, and Chandra Reed for their contributions to this work. Thomas J. Stohlgren was partially supported by the Bioenergy Alliance Network of the Rockies (BANR), USDA UV-B Monitoring and Research Program and USDA CSREES/NRI 2008-35615-04666. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

1. Field Data Collection
<ol>
	<li>Derive field data for tamarisk from a vector polygon dataset collected by the Tamarisk Coalition in a state-wide survey in 2005 and 200617. 
	NOTE: The data were obtained through an intensive ground survey, where field technicians mapped all tamarisk stands along the Arkansas River using Global Positioning Systems (GPS) and aerial photographs.</li>
	<li>Within these polygons, generate 499 random points (i.e., presences) to train the models. Place two sets of constraints on the random points: (1) each point is required to be &#62; 30 m from any polygon edge to ensure it is within the tamarisk stand at a Landsat 5 TM resolution; and (2) each point is required to be &#8805; 60 m from any neighboring point, to ensure a robust distribution across the study site and minimize spatial autocorrelation.</li>
	<li>Compile field data in a MS Excel file in three columns with the headings labeled &#34;response&#34;, &#34;X&#34;, and &#34;Y&#34;, where the response values are (1) for presence, UTM easting for X, and UTM northing for Y. Save this file in .csv format for use in SAHM (FieldData module).</li>
	<li>Generate an additional 100 random points within the Tamarisk Coalition&#39;s tamarisk polygons and another 100 random points outside the polygons as independent test data to evaluate model results. Save this file in .csv format for use in SAHM (FieldData module).</li>
</ol>
2. Predictor Variables
<ol>
	<li>Download Landsat 5 surface reflectance L4-5 TM images (Path 32, Row 34) from the U.S. Geological Survey&#39;s Global Visualization Viewer/EarthExplorer (http://earthexplorer.usgs.gov/). The scenes include the extent sampled by the Colorado Tamarisk Coalition in 2005 and 200617. To determine the months used for the models, select scenes that are largely cloud-free (i.e., &#60; 10% cloud cover) and representative of months thought to have distinguishable tamarisk phenology; these were October 2004, April 2005, May 2005, June 2005, July 2005, September 2005, April 2006, May 2006, June 2006, July 2006, August 2006, September 2006, and November 2006.</li>
	<li>Download the Remote Sensing Indices Derivation Tool from 
	https://github.com/rander38/Remote-Sensing-Indices-Derivation-Tool.</li>
	<li>Run the Python script in the tool, either GDAL or Arcpy version; GDAL is recommended.</li>
	<li>Select the appropriate satellite sensor, desired indices, and set the input image file and the output folder where files will be stored (Figure 1). We exported the individual bands and used the NDVI, SAVI, and Tasseled Cap brightness, greenness, and wetness indices derived from each of the Landsat TM scenes. Note that any indices can be modified or added by altering the Sensors_Formulas_RSIDT.ini file.</li>
</ol>
<img alt="Figure 1" src="/files/ftp_upload/54578/54578fig1.jpg" /> 
Figure 1. Remote Sensing Indices Derivation Tool GUI.
<ol start="5">
	<li>Run the tool and verify the output files visually in ArcMap v. 10.0 (ESRI, Redlands, CA) or another geographic information system software.</li>
</ol>
3. Software for Assisted Habitat Modeling (SAHM) (Figure 2)
<img alt="Figure 1" src="/files/ftp_upload/54578/54578fig2.jpg" /> 
Figure 2. The entire SAHM workflow includes input data, preprocessing, preliminary model analysis and decision, correlative models, and output routines.
<ol>
	<li>For running SAHM, first download files (including VisTrails) from the U.S. Geological Survey website at https://my.usgs.gov/catalog/RAM/SAHM. Refer to the User&#8217;s Guide at the same website for detailed instructions on downloading and installing SAHM. Note that the website also has a SAHM tutorial and data for additional instructions.</li>
	<li>To develop the tamarisk species distribution models, use the SAHM_tutorial_2.0.vt file that comes with the package download (in the SAHM examples folder). Within the History view, select the Independent locations workflow. The other workflow examples can be selected and depend on study objectives; descriptions are provided with each. Select Pipeline.</li>
	<li>Set the output folder, by going to Packages and then SAHM &#62; Change Session Folder. Throughout the workflow development process, detailed descriptions of each step and options may be found by selecting the Documentation tab located on the right side of the SAHM viewer screen. All modules listed in the following methods can be found on the left side of the SAHM window under the SAHM tab.</li>
	<li>Next, direct SAHM to the field data that will be used to train species distribution models.
	<ol>
		<li>Click on the TemplateLayer module. Browse to the raster that will be used as a mask and to define the projection, cell size and extent of the analyses.</li>
		<li>Click on the FieldData module in the left side of the workflow. Browse to .csv (i.e., training.csv) file of field data (presence points or presence and absence points) within the field data module.</li>
		<li>Click on the PredictorListFile module and browse to .csv file list (e.g., file containing the full path to all predictors to use in the model- refer to user guide).</li>
	</ol>
	</li>
	<li>Next, conduct preprocessing steps.
	<ol>
		<li>Click on the FieldDataQuery module and fill in response column with column heading for response (i.e., the column names in the FieldData.csv), X and Y columns.</li>
		<li>Click on the MDSBuilder module. Set the backgroundPointField to 10,000. 
		NOTE: If you are using presence and absence data for the species being modeled, you do not need to change the backgroundPointField; you would include these locations with a response (0) in the Field Data.csv. It is optional to set a backgroundProbSurf if you want to constrain background point selection within an area by pointing to a raster surface with values ranging from 0 to 100 (these values represent the probability that a randomly generated point will be retained should it fall within a particular cell). For this study, use a backgroundProbSurf with values of 100 within a 5,000 m buffer of the Arkansas River and 0 for areas outside this buffer (based on the total area sampled by the tamarisk coalition).</li>
	</ol>
	</li>
	<li>Next, specify the species distribution modeling algorithms to be used.
	<ol>
		<li>Note that the BoostedRegressionTree, GLM, MARS, and RandomForest modules are already set up in the Independent locations workflow. Add the MAXENT module to the workflow to test all five models. Connect it to the CovariateCorrelationAndSelection module. 
		NOTE: Begin with default settings for all models; these can be modified based on study objectives (see model documentation for more details).</li>
		<li>Add a ModelOutputViewer module and connect it to the MAXENT module; change column to 5 and row to 1. The ModelOutputViewer produces a spreadsheet that can be used for model result comparisons.</li>
		<li>Click on the OutputName module and type a subfolder name.</li>
	</ol>
	</li>
	<li>Next, add a module that creates an ensemble of model outputs. This module produces two output maps; one with the average continuous probability of all the included outputs and a second with the count of the number of models with a positive binary probability.
	<ol>
		<li>Add an EnsembleBuilder module to the workflow. It is optional to set a Threshold Metric; for this study, choose AUC and leave the Threshold Value at the default value of 0.75. This will ensure that only models with an AUC value greater than or equal to 0.75 are included in the ensemble map outputs. Connect the BoostedRegressionTree, GLM, MARS, RandomForest, and MAXENT modules to the EnsembleBuilder.</li>
	</ol>
	</li>
	<li>Next, direct the models to the independent test data.
	<ol>
		<li>Click on the other FieldData module (right side of the workflow) and browse to the .csv file that contains model validation data. These are the 200 presence and absence points generated in step 1.4 of the PROTOCOL.</li>
		<li>Click on the FieldDataQuery module and make sure the response, x, and y columns match the columns in the field data .csv.</li>
		<li>Add an ApplyModel module and connect it to the MAXENT module. Add a ModelOutputViewer module and connect it to this ApplyModel module; change the column to 5 and the row to 1. In the menu, select Packages -&#62; SAHM -&#62; Change Processing Mode. Since you want to run more than one model, select single models sequentially (n - 1 cores each). This will expedite the execution time of the models by taking advantage of multiple computer cores.</li>
	</ol>
	</li>
	<li>Next, execute the species distribution models.
	<ol>
		<li>Save the file .vt and then click on Execute.</li>
		<li>When the CovariateCorrelationAndSelection widget (Figure 3) appears, de-select one of each correlated pair of variables where the correlation coefficient is |r| &#8805; 0.7 (based on % Deviance Explained from a univariate generalized additive model seen in on the left side of this widget (Figure 3) and informed ecological decisions; in this study prioritize the selection of at least one covariate for each month to capture phenological variation in tamarisk). The Number of Plots in the viewing window can be changed by typing in a number (the default is 8) and clicking on Update.</li>
		<li>After finalizing covariate selection, select OK at the bottom of the CovariateCorrelationAndSelection widget (Figure 3). For this study, the following 9 variables were retained: July_30_2006_Brightness, June_09_2005_SAVI, Sept_16_2006_SAVI, May_24_2005_B4, Oct_28_2004_NDVI, April_22_2005_Brightness, April_09_2006_SAVI, Aug_31_2006_B4, and Nov_19_2006_SAVI. The species distribution models will execute after selecting OK.</li>
	</ol>
	</li>
</ol>
<img alt="Figure 3" src="/files/ftp_upload/54578/54578fig3.jpg" /> 
Figure 3. Covariate correlation and selection SAHM interface.
<ol start="10">
	<li>Results output. 
	NOTE: After the models have completed, a VisTrails Spreadsheet appears for model comparison (Figure 4)
	<ol>
		<li>Compare AUC plots, text outputs, response curves, calibration, confusion matrix, residuals, and response curves across models. 
		NOTE: Output from SAHM models include a CovariateCorrelationOutputMDS, confusion matrix, residual plot, calibration plot, model evaluation plot, variable importance plot, a response curves folder, an expanded output folder, bin map, mess map, MoD map, prob map, residual map, and an output.txt; see the SAHM User Guide for more description.</li>
	</ol>
	</li>
</ol>
<img alt="Figure 4" src="/files/ftp_upload/54578/54578fig4.jpg" /> 
Figure 4. VisTrails spreadsheets can be used to evaluate model output. This is the AUC model comparison for the training data; from left to right the models are BRT, GLM, MARS, RF, and Maxent, respectively.

<ol>
	<li>Evangelista, P. H., Stohlgren, T. J., Morisette, J. T., Kumar, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mapping+invasive+tamarisk+(Tamarix):+a+comparison+of+single-scene+and+time-series+analyses+of+remotely+sensed+data.">Mapping invasive tamarisk (Tamarix): a comparison of single-scene and time-series analyses of remotely sensed data.</a> Remote Sensing. 1, 519-533 (2009).</li><li>DiTomaso, J. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Impact,+biology,+and+ecology+of+saltcedar+(Tamarix+spp.)+in+the+southwestern+United+States.">Impact, biology, and ecology of saltcedar (Tamarix spp.) in the southwestern United States.</a> Weed Technology. 12, 326-336 (1998).</li><li>Evangelista, P., Kumar, S., Stohlgren, T., Crall, A., Newman, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Modeling+above-ground+biomass+of+Tamarisk+ramosissima+in+the+Arkansas+River+Basin+of+Southeastern+Colorado,+USA.">Modeling above-ground biomass of Tamarisk ramosissima in the Arkansas River Basin of Southeastern Colorado, USA.</a> Western North American Naturalist. 67 (4), 503-509 (2007).</li><li>Hirano, A., Madden, M., Welch, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Hyperspectral+image+data+for+mapping+wetland+vegetation.">Hyperspectral image data for mapping wetland vegetation.</a> Wetlands. 23 (2), 436-448 (2003).</li><li>Ge, S., Carruthers, R., Gong, P., Herrera, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Texture+analysis+for+mapping+Tamarix+pariviflora+using+aerial+photographs+along+Cache+Creek,+California.">Texture analysis for mapping Tamarix pariviflora using aerial photographs along Cache Creek, California.</a> Environmental Monitoring and Assessment. 114, 65-83 (2006).</li><li>Hamada, Y., Stow, D. A., Coulter, L. L., Jafolla, J. C., Hendricks, L. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Detecting+tamarisk+species+(Tamarisk+spp.)+in+riparian+habitats+of+Southern+California+using+high+spatial+resolution+hyperspectral+imagery.">Detecting tamarisk species (Tamarisk spp.) in riparian habitats of Southern California using high spatial resolution hyperspectral imagery.</a> Remote Sensing of Environment. 109, 237-248 (2007).</li><li>York, P., Evangelista, P., Kumar, S., Graham, J., Flather, C., Stohlgren, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+Habitat+Overlap+Analysis+derived+from+Maxent+for+Tamarisk+and+the+Southwestern+Willow+Flycatcher.">A Habitat Overlap Analysis derived from Maxent for Tamarisk and the Southwestern Willow Flycatcher.</a> Frontiers of Earth Science. 5 (2), 120-129 (2011).</li><li>Myneni, R. B., Ramakrishna, R., Nemani, R., Running, S. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Estimation+of+global+leaf+area+index+and+absorbed+par+using+radiative+transfer+models.">Estimation of global leaf area index and absorbed par using radiative transfer models.</a> Geoscience and Remote Sensing. 35 (6), 1380-1393 (1997).</li><li>Todd, S. W., Hoffer, R. M., Milchunas, D. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Biomass+estimation+on+grazed+and+ungrazed+rangelands+using+spectral+indices.">Biomass estimation on grazed and ungrazed rangelands using spectral indices.</a> International Journal of Remote Sensing. 19 (3), 427-438 (1998).</li><li>Rouse, J. W., Haas, R. H., Schell, J. A., Deering, D. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Monitoring+vegetation+systems+in+the+Great+Plains+with+ERTS.">Monitoring vegetation systems in the Great Plains with ERTS.</a> Proceedings of the Third Earth Resources Technology Satellite-1 Symposium. , (1974).</li><li>Huete, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+Soil-Adjusted+Vegetation+Index+(SAVI).">A Soil-Adjusted Vegetation Index (SAVI).</a> Remote Sensing of Environment. 25, 295-309 (1988).</li><li>Kauth, R. J., Thomas, G. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+tasselled+cap+-+a+graphic+description+of+the+spectral-temporal+development+of+agricultural+crops+as+seen+in+Landsat.">The tasselled cap - a graphic description of the spectral-temporal development of agricultural crops as seen in Landsat.</a> Proceedings of the Symposium on Machine Processing of Remotely Sensed Data. , 41-51 (1976).</li><li>Cohen, W. B., Spies, T. A., Fiorella, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Estimating+the+age+and+structure+of+forests+in+a+multi-ownership+landscape+of+western+Oregon,+USA.">Estimating the age and structure of forests in a multi-ownership landscape of western Oregon, USA.</a> International Journal of Remote Sensing. 16, 721-746 (1995).</li><li>Jin, S., Sader, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+time+series+tasselled+cap+wetness+and+the+normalized+difference+moisture+index+in+detecting+forest+disturbances.">Comparison of time series tasselled cap wetness and the normalized difference moisture index in detecting forest disturbances.</a> Remote Sensing of Environment. 94 (3), 364-372 (2005).</li><li>Crist, E. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+TM+Tasseled+Cap+equivalent+transformation+for+reflectance+factor+data.">A TM Tasseled Cap equivalent transformation for reflectance factor data.</a> Remote Sensing of Environment. 17 (3), 301-306 (1985).</li><li>Lindauer, I. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+comparison+of+the+plant+communities+of+the+South+Platte+and+Arkansas+River+drainages+in+eastern+Colorado.">A comparison of the plant communities of the South Platte and Arkansas River drainages in eastern Colorado.</a> The Southwestern Naturalist. 28 (3), 249-259 (1983).</li><li>Guisan, A., Zimmerman, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Predictive+habitat+distribution+models+in+ecology.">Predictive habitat distribution models in ecology.</a> Ecological Modeling. 135, 147-186 (2000).</li><li>Friedman, J. H., Hastie, T., Tibshirani, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Additive+logistic+regression:+a+statistical+view+of+boosting.">Additive logistic regression: a statistical view of boosting.</a> Annals of Statistics. 28 (2), 337-407 (2000).</li><li>Breiman, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Random+forests.">Random forests.</a> Machine Learning. 45 (1), 5-32 (2001).</li><li>Friedman, J. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Multivariate+adaptive+regression+splines.">Multivariate adaptive regression splines.</a> Annals of Statistics. 19 (1), 1-141 (1991).</li><li>McCullagh, P., Nelder, J. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Generalized Linear Models, 2nd ed. , (1989).</li><li>Phillips, S. J., Anderson, R. P., Schapire, R. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Maximum+entropy+modeling+of+species+geographic+distributions.">Maximum entropy modeling of species geographic distributions.</a> Ecological Modelling. 190 (3-4), 231-259 (2006).</li><li>Araujo, M. B., New, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ensemble+forecasting+of+species+distributions.">Ensemble forecasting of species distributions.</a> Trends in Ecology and Evolution. 22, 42-47 (2007).</li><li>Elith, J., Graham, C. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Do+they?+How+do+they?+Why+do+they+differ?+On+finding+reasons+for+differing+performances+of+species+distribution+models.">Do they? How do they? Why do they differ? On finding reasons for differing performances of species distribution models.</a> Ecography. 32, 66-77 (2009).</li><li>Freire, J., Silva, C., Callahan, S., Santos, E., Schedegger, C., Moreau, L., Foster, I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Managing+rapidly-evolving+scientific+workflows.">Managing rapidly-evolving scientific workflows.</a> International Provenance and Annotation Workshop (IPAW). , 10-18 (2006).</li><li>Morisette, J. T., Jarnevich, C. S., Holcombe, T. R., Talbert, C. B., Ignizio, D., Talbert, M. K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=VisTrails+SAHM:+visualization+and+workflow+management+for+species+habitat+modeling.">VisTrails SAHM: visualization and workflow management for species habitat modeling.</a> Ecography. 36 (2), 129-135 (2013).</li><li>Fielding, A. H., Bell, J. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+review+of+methods+for+the+assessment+of+prediction+errors+in+conservation+presence/absence+models.">A review of methods for the assessment of prediction errors in conservation presence/absence models.</a> Environmental Conservation. 24, 38-49 (1997).</li><li>Allouche, O., Tsoar, A., Kadmon, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Assessing+the+accuracy+of+species+distribution+models:+prevalence,+kappa+and+the+true+skill+statistic+(TSS).">Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS).</a> Journal of Applied Ecology. 43 (6), 1223-1232 (2006).</li><li>Stohlgren, T. J., Ma, P., Kumar, S., Rocca, M., Morisette, J., Jarnevich, C. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ensemble+habitat+mapping+of+invasive+plant+species.">Ensemble habitat mapping of invasive plant species.</a> Risk Analysis. 30, 224-235 (2010).</li><li>Dormann, C. F., Purschke, O., Marquez, J. R. G., Lautenbach, S., Schrader, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Components+of+uncertainty+in+species+distribution+analysis:+A+case+study+of+the+great+grey+shrike.">Components of uncertainty in species distribution analysis: A case study of the great grey shrike.</a> Ecology. 89, 3371-3386 (2008).</li><li>Marmion, M., Parviainen, M., Luoto, M., Heikkinen, R. K., Thuiller, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Evaluation+of+consensus+methods+in+predictive+species+distribution+modelling.">Evaluation of consensus methods in predictive species distribution modelling.</a> Diversity and Distributions. 15, 59-69 (2009).</li></ol>

The authors have no competing financial interests or conflict of interest.

Our results demonstrate fitting BRT, RF, MARS, GLM, and Maxent with presence points for tamarisk and a time-series of remotely sensed Landsat satellite imagery data can distinguish tamarisk on the landscape and is an effective alternative to traditional single-scene classification methods. It is clear from our results that June is a particularly important time for detecting tamarisk within our study area; this agrees with Evangelista et al. 20091 which indicated June Wetness was the most important predictor for tamarisk occurrence in this area based on a Maxent model fit with a time series of Landsat imagery.

The other spectral indices and bands that were included in the BRT, RF, MARS, and Maxent models may further distinguish tamarisk from soil substrate, other deciduous trees including cottonwood (Populus spp.) and willow (Salix spp.), or irrigated agriculture that is common in the lower Arkansas River basin. Other GIS layers, such as topography, soil types, or climate data could also be considered as covariates and included in these models, but we recommend keeping these to a minimum if the goal is to detect current species distribution on the landscape rather than predict potential occurrence or suitable habitat.
The models tested for our research provided strong analytical capacity and multiple options for evaluating results. Having all of these correlative models within a single framework, such as SAHM, allows the formalization and tractable recording of the modeling process. Pre- and post-processing of response and predictor variables are standardized in SAHM, allowing better and efficient model comparisons, while workflows record each step of the analyses facilitating modification, iteration and replication.
Ensemble mapping aims to combine the strengths of several correlative models, while minimizing the weakness of any one model30. We believe this was the case in our study; however, we caution that models that underperform (i.e., under-predict or over-predict) can weaken the overall results. The limited use of ensemble mapping in the literature has had favorable results, but most of these approaches have attempted to &#34;predict&#34; species occurrence rather than &#34;detect.&#34; Furthermore, ensemble mapping allows for a visual assessment of uncertainty among the different modeling methods, identifying levels of model agreement. Most often it is the choice of modeling method (e.g., GLM versus BRT) that has greatest quantifiable impact on model results rather than other decisions in the modeling process such as location data uncertainty31. Although we believe that our best tamarisk map is where all five models are in agreement, further testing and using various methods of ensemble mapping is recommended (e.g., weighted by AUC)32, and best validated through independent field observations. In summary, these methods can easily be tailored to model the distribution of other species using environmental variables derived for a given study region in SAHM.

Riparian and wetland ecosystems throughout the southwestern United States are being threatened by the invasion of tamarisk (Tamarix spp.), a non-native woody shrub introduced from Eurasia in the 1800s1. Tamarisk has many physiological mechanisms that allow the genus to exploit water resources, out-compete native species, and alter ecosystem processes1-2. Mapping tamarisk distributions for assessing environmental impacts and formulating effective control strategies are high priorities for resource managers. Although ground surveys remain regularly used, they are impractical for extremely large areas due to the associated costs of labor, time, and logistics.
Satellite remote sensing has played an important, but limited, role in the detection and mapping of tamarisk infestations. Conventional classification analyses and remote sensing software have had marginal success3-5. Several recent studies have explored non-traditional approaches to detect invasive plants using remote sensing data1,6. Tamarisk, like many invasive plants, exhibits phenological variation throughout the growing season that differs from native riparian species&#39; phenology. In some areas, for example, tamarisk leaf-out is before some native riparian plants, and tamarisk retains its foliage longer than other native species. By using spectral bands and spectral indices derived from a time-series of satellite data throughout the growing season, we can distinguish tamarisk from native plants based on these phenological differences1,6. Building on the work of Evangelista et al. 20091, in this study we incorporated individual bands 1-7 from a time-series of Landsat 5 Thematic Mapper (TM) satellite imagery and derived normalized difference vegetation index (NDVI), soil-adjusted vegetation index (SAVI), and tasseled cap transformations from these bands. Normalized difference vegetation index (NDVI) is one of the most commonly used spectral indices for estimating vegetation biomass, canopy cover, and leaf area indices8-9, and is a non-linear transformation of the ratio between the visible (red) and near-infrared bands10. Soil-adjusted vegetation index (SAVI) is a modified NDVI used to minimize the effects of soil background on vegetation indices11. Tasseled cap transformations are weighted composites of the six Landsat bands into three orthogonal bands that measure soil brightness (tasseled cap, band 1), vegetation greenness (tasseled cap, band 2), and soil/vegetation wetness (tasseled cap, band 3) and are often used to distinguish vegetation composition, age class, and structure12-14. We used the coefficients reported in Crist (1985)15 for all tasseled cap transformations.
In this study, we test five species distribution models with a time-series of spectral bands and vegetation indices derived from Landsat 5 TM to map tamarisk along the lower Arkansas River in southeastern Colorado, USA. The Arkansas River, spanning 2,364 km (1,469 mi), is the second largest tributary in the Missouri-Mississippi system. Its watershed covers 435,123 km2 (168,002 mi2) with headwaters in the Colorado Rocky Mountains. From its origin at 2,965 m, the Arkansas drops considerably in elevation, leveling out near Pueblo, CO, and meandering through agricultural lands and short-grass prairie. The river is subject to seasonal flooding and is relied on for municipal and agricultural water use in Rocky Ford, La Junta, and Lamar, before continuing into Kansas, Oklahoma, and Arkansas where it flows into the Mississippi River. Tamarisk was first observed on the Arkansas River by R. Niedrach in 1913 near the present-day town of Lamar16. Today, it has been estimated that tamarisk covers more than 100 km2 between Pueblo and the Kansas state line, with an additional 60 km2 along the tributaries of the Arkansas River17. The study area includes irrigation ditches, wetlands, agricultural land, and the confluences of several tributaries; all with varying degrees of tamarisk infestation. Ranching and agriculture are the primary land-uses adjacent to the riparian corridors consisting largely of alfalfa, hay, corn, and winter wheat.
Species distribution models rely on geo-referenced occurrences (i.e., latitude, longitude) to identify relationships between a species&#39; occurrence and its environment18. The environmental data can include multiple remote sensing and other spatial layers. The five species distribution models we tested include boosted regression trees (BRT)19, random forests (RF)20, multivariate adaptive regression splines (MARS)21, a generalized linear model (GLM)22, and Maxent23. These five model algorithms are among the most commonly employed for species distribution modeling, and a number of studies have demonstrated their efficacy24-25. We used the Software for Assisted Habitat Modeling (SAHM) v. 2.0 modules to execute the five models, which are contained in VisTrails v.2.2.226 visualization and processing software. There are several advantages to using SAHM for comparative modeling. In addition to the formalization and tractable recording of modeling processes, SAHM allows users to work with multiple species distribution model algorithms that, individually, have disparate interfaces, software and file formatting27. SAHM produces consistent threshold-independent and threshold-dependent evaluation metrics to evaluate model performance. One of these is Area Under the Receiver Operating Characteristic Curve (AUC), a threshold independent metric that evaluates ability of a model to discriminate presence from background28. An AUC value of 0.5 or less indicates model predictions are not better or worse than random; values between 0.5 and 0.70 indicate poor performance; and values increasing from 0.70 to 1.0 indicate progressively higher performance. Another metric is percent correctly classified (PCC), a threshold dependent metric that weighs sensitivity and specificity based on a user-defined threshold metric; sensitivity measures the percentage of observed presences classified as suitable and specificity measures the percentage of background locations classified as unsuitable. Yet another metric is True Skill Statistic (TSS = sensitivity + specificity - 1), which places more weight on model sensitivity than specificity, with values ranging between -1 and 1 where values &#62; 0 indicate better model performance than chance29.
To map tamarisk using model output, we constructed binary classifications using the threshold that equalizes sensitivity and specificity to define the presence or absence of tamarisk. These individual model derived maps were then summed to create an ensemble map30. Ensemble maps combine the predictions of individual species distribution models to produce a classified map that ranks the collective agreement of the models tested. For example, an ensemble cell value of one indicates that only one model classified that cell as suitable habitat, whereas a value of five indicates that all five models classified the cell as suitable habitat. One advantage to this approach is that ensemble maps yield a lower mean error than any individual model. It also allows users to visually compare the performance of each model tested. Our overall goal was to provide a detailed description of these methods that can be tailored to model the current distribution of species on the landscape.

<table border="0" cellpadding="0" cellspacing="0" width="949">
	<tbody>
		<tr height="20">
			<td height="20" style="height:20px;width:259px;">Earth Explorer</td>
			<td style="width:65px;">USGS</td>
			<td style="width:533px;">http://earthexplorer.usgs.gov</td>
			<td style="width:93px;">Open Access: Yes</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Remote Sensing Indices Derivation Tool</td>
			<td>github</td>
			<td>https://github.com/rander38/Remote-Sensing-Indices-Derivation-Tool</td>
			<td>Open Access: Yes</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Software for Assisted Habitat Modeling</td>
			<td>USGS</td>
			<td>https://my.usgs.gov/catalog/RAM/SAHM</td>
			<td>Open Access: Yes</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">ArcGIS v.10.3&#160;</td>
			<td>Esri</td>
			<td>https://www.arcgis.com/features/</td>
			<td>Open Access: No</td>
		</tr>
	</tbody>
</table>

integrating remote sensing with species distribution models; mapping tamarisk invasions using the software for assisted habitat modeling (sahm)

Early detection of invasive plant species is vital for the management of natural resources and protection of ecosystem processes. The use of satellite remote sensing for mapping the distribution of invasive plants is becoming more common, however conventional imaging software and classification methods have been shown to be unreliable. In this study, we test and evaluate the use of five species distribution model techniques fit with satellite remote sensing data to map invasive tamarisk (Tamarix spp.) along the Arkansas River in Southeastern Colorado. The models tested included boosted regression trees (BRT), Random Forest (RF), multivariate adaptive regression splines (MARS), generalized linear model (GLM), and Maxent. These analyses were conducted using a newly developed software package called the Software for Assisted Habitat Modeling (SAHM). All models were trained with 499 presence points, 10,000 pseudo-absence points, and predictor variables acquired from the Landsat 5 Thematic Mapper (TM) sensor over an eight-month period to distinguish tamarisk from native riparian vegetation using detection of phenological differences. From the Landsat scenes, we used individual bands and calculated Normalized Difference Vegetation Index (NDVI), Soil-Adjusted Vegetation Index (SAVI), and tasseled capped transformations. All five models identified current tamarisk distribution on the landscape successfully based on threshold independent and threshold dependent evaluation metrics with independent location data. To account for model specific differences, we produced an ensemble of all five models with map output highlighting areas of agreement and areas of uncertainty. Our results demonstrate the usefulness of species distribution models in analyzing remotely sensed data and the utility of ensemble mapping, and showcase the capability of SAHM in pre-processing and executing multiple complex models.

Statistical evaluations of BRT, RF, MARS, GLM, and Maxent based on the independent test dataset indicated all five models performed relatively well in detecting tamarisk; there was little difference between threshold independent and threshold dependent evaluation metrics among models. AUC values were &#62; 0.88, percent correctly classified values were&#160;&#62; 77%, sensitivities and specificities were&#160;&#62; 0.77, and TSS were&#160;&#62; 0.54 (Table 1). An ensemble of binary model outputs revealed much model agreement in areas along the Arkansas River (Figure 5). The MESS (multivariate environmental similarity surface) map outputs for each model indicated the available environment of the study area was well sampled (Figure 6), further increasing our confidence in the ensemble approach.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Model</td>
			<td>AUC</td>
			<td>PCC</td>
			<td>Sensitivity</td>
			<td>Specificity</td>
			<td>TSS</td>
		</tr>
		<tr>
			<td>BRT</td>
			<td>0.91</td>
			<td>85</td>
			<td>0.85</td>
			<td>0.85</td>
			<td>0.70</td>
		</tr>
		<tr>
			<td>RF</td>
			<td>0.92</td>
			<td>85</td>
			<td>0.85</td>
			<td>0.85</td>
			<td>0.70</td>
		</tr>
		<tr>
			<td>MARS</td>
			<td>0.90</td>
			<td>82</td>
			<td>0.82</td>
			<td>0.82</td>
			<td>0.64</td>
		</tr>
		<tr>
			<td>GLM</td>
			<td>0.88</td>
			<td>77</td>
			<td>0.77</td>
			<td>0.77</td>
			<td>0.54</td>
		</tr>
		<tr>
			<td>Maxent</td>
			<td>0.92</td>
			<td>84</td>
			<td>0.83</td>
			<td>0.84</td>
			<td>0.67</td>
		</tr>
	</tbody>
</table>
Table 1. Threshold independent (AUC) and threshold dependent (PCC, Sensitivity, Specificity, and TSS) evaluation metrics for BRT, RF, MARS, GLM, and Maxent models fit to an independent test dataset of tamarisk presence and absence.
<img alt="Figure 5" src="/files/ftp_upload/54578/54578fig5.jpg" /> 
Figure 5. Ensemble results combining BRT, GLM, MARS, RF, and Maxent binary output maps in ArcGIS. Areas are colored by number of models in agreement, from 0 (no color) to 5 (red). Note the colored area in the northwest corner of the prediction; this line is an artifact of Landsat imagery; therefore model results should be taken with caution in this region. <a href="https://www.jove.com/files/ftp_upload/54578/54578fig5large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 6" src="/files/ftp_upload/54578/54578fig6.jpg" /> 
Figure 6. Multivariate environmental similarity surface (MESS) output. <a href="https://www.jove.com/files/ftp_upload/54578/54578fig6large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Of the nine predictors used, June 30 2006 Brightness was the most important variable for all five models (Table 2). This was the only variable retained by GLM based on stepwise Akaike information criterion (AIC; this is the default for GLM model selection in SAHM), however it is important to note this model also included a squared term of this variable. RF and Maxent retain all variables by default.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Predictor</td>
			<td>BRT</td>
			<td>RF</td>
			<td>MARS</td>
			<td>GLM</td>
			<td>Maxent</td>
		</tr>
		<tr>
			<td>July 30 2006 Brightness</td>
			<td>41.60</td>
			<td>34.11</td>
			<td>76.78</td>
			<td>100</td>
			<td>67.27</td>
		</tr>
		<tr>
			<td>August 31 2006 Band 4</td>
			<td>6.35</td>
			<td>5.87</td>
			<td>5.16</td>
			<td>0</td>
			<td>2.82</td>
		</tr>
		<tr>
			<td>June 09 2005 SAVI</td>
			<td>13.67</td>
			<td>14.09</td>
			<td>9.14</td>
			<td>0</td>
			<td>9.75</td>
		</tr>
		<tr>
			<td>April 22 2005 Brightness</td>
			<td>6.29</td>
			<td>6.30</td>
			<td>0</td>
			<td>0</td>
			<td>0.43</td>
		</tr>
		<tr>
			<td>October 28 2004 NDVI</td>
			<td>5.66</td>
			<td>8.25</td>
			<td>0</td>
			<td>0</td>
			<td>2.94</td>
		</tr>
	</tbody>
</table>
Table 2. Relative importance of predictors in each model.

Watch this Scientific Journal Video about Integrating Remote Sensing with Species Distribution Models; Mapping Tamarisk Invasions Using the Software for Assisted Habitat Modeling SAHM at JoVE.com

Integrating Remote Sensing with Species Distribution Models; Mapping Tamarisk Invasions Using the Software for Assisted Habitat Modeling SAHM

We demonstrate the utility of remotely sensed data and the newly developed Software for Assisted Habitat Modeling (SAHM) in predicting invasive species occurrence on the landscape. An ensemble of predictive models produced highly accurate maps of tamarisk (Tamarix spp.) invasion in Southeastern Colorado, USA when assessed with subsequent field validations.

Integrating Remote Sensing with Species Distribution Models; Mapping Tamarisk Invasions Using the Software for Assisted Habitat Modeling (SAHM)

Early detection of invasive plant species is vital for the management of natural resources and protection of ecosystem processes. The use of satellite ...

Research

JoVE Journal

Environment

An Optimized Enrichment Technique for the Isolation of Arthrobacter Bacteriophage Species from Soil Sample Isolates

An Optimized Enrichment Technique for the Isolation of Arthrobacter Bacteriophage Species from Soil Sample Isolates

Bacteriophage isolation from environmental samples has been performed for decades using principles set forth by pioneers in microbiology. The isolation of ...

Protocols for Robust Herbicide Resistance Testing in Different Weed Species

Robust protocols to test putative herbicide resistant weed populations at whole plant level are essential to confirm the resistance status. The presented ...

Remote Sensing Evaluation of Two-spotted Spider Mite Damage on Greenhouse Cotton

The objective of this study was to evaluate a ground-based multispectral optical sensor as a remote sensing tool to assess foliar damage caused by the ...

Experimental Methods of Dust Charging and Mobilization on Surfaces with Exposure to Ultraviolet Radiation or Plasmas

Electrostatic dust transport has been hypothesized to explain a number of observations of unusual planetary phenomena. Here, it is demonstrated using ...

Methods for Image-based Surveys of Benthic Macroinvertebrates and Their Habitat Exemplified by the Drop Camera Survey for the Atlantic Sea Scallop

Underwater imaging has long been used in the field of marine ecology but decreasing costs of high-resolution cameras and data storage have made the ...

Use of Principal Components for Scaling Up Topographic Models to Map Soil Redistribution and Soil Organic Carbon

Landscape topography is a critical factor affecting soil formation and plays an important role in determining soil properties on the earth surface, as it ...

Combining Eye-tracking Data with an Analysis of Video Content from Free-viewing a Video of a Walk in an Urban Park Environment

As individuals increasingly live in cities, methods to study their everyday movements and the data that can be collected becomes important and valuable. ...

Composition and Distribution Analysis of Bioaerosols Under Different Environmental Conditions

Variable microorganisms in particulate matter (PM) under different environmental conditions may have significant impacts on human health. In this study, ...

Measuring the Shape and Size of Activated Sludge Particles Immobilized in Agar with an Open Source Software Pipeline

Experimental bioreactors, such as those treating wastewater, contain particles whose size and shape are important parameters. For example, the size and ...

Identifying Per- and Polyfluorinated Chemical Species with a Combined Targeted and Non-Targeted-Screening High-Resolution Mass Spectrometry Workflow

Historical and emerging per- and polyfluoroalkyl substances (PFASs) have garnered significant interest from the public and government agencies from the ...