How to Model Stream Temperature Using ArcMap

Transcription

1 How to Model Stream Temperature Using ArcMap Take note: Assumption before proceeding: A temperature point file has been attributed with TauDEM variables. There are three processes described in this document. We will: Develop a regression model from predictors and response variables using Excel (you may use other analytical software) Apply the regression equation in ArcMap so that you can make predictions about temperature scenarios for a given stream network Describe some GIS applications for the temperature model 1. Develop a regression model using Excel 1. In Excel, open the.dbf file for the temperature shapefile. Do not modify the.dbf file and save it as an Excel file. Keep the variables that are relevant to modeling stream temperature and delete the other variables if you wish. Re-save the Excel file. (Note: TauDEM reports Downstream Contributing Area (DS_CONT_AR) in meters squared; in ArcMap we convert it to hectares and rename it ds_cont_ha. Below it is named Contributing Area (ha). 1

2 2. Run the regression model. On the Excel toolbar select Tools > Data Analysis > Regression > OK (If Data Analysis is not visible, select Add-Ins > Analysis ToolPak > OK. Data Analysis should then appear). 3. When the dialogue box opens, select the options as shown below: a. Input: i. Input Y Range: highlight the data cells of the response variable. Do not highlight the variable name ii. Input X Range: highlight the data cells of all predictor variables (we highlighted all 6 columns of predictors at once) Do not highlight the variable name b. Output options: i. We chose Output Range and selected a cell in a blank column next to the data columns. You may choose to have the data displayed on a separate worksheet or a separate workbook. Pick the option that you desire c. Residuals: Select any plots that you would like to create d. OK 2

3 The regression output appears as follows: The R Square indicates the proportion of variation explained by the model (i.e., 60%), the coefficients define the regression equation, and the P-values indicate which coefficients are significant 4. The regression equation for this model is: Y predicted temperature = * [Elevation] * [Area_ha] * [Slope] * [Lake Cr Effect] * [Fire Effect] * [Lake Effect] 3

4 2. Apply the Regression Equation in ArcMap 1. Open ArcMap and add the TauDEM stream layer (other layers may be added if desired) 2. Create the new modeled temperature variable: right click on the TauDEM shapefile > Open Attribute Table (Note that the variable names in the Attribute Table are slightly different from the names we used in the Excel spreadsheet) In the Attribute Table, Choose Options > Add Field > Name the new variable PrdSumMWMT (predicted summer MWMT) Type: Double Field Properties Precision: 6 (the number of digits that can be stored in a number field, including decimal digits) Scale: 4 (number of decimal digits) OK 4

5 3. In the Attribute Table, right click on the new variable > Field Calculator. Paste the regression formula into the blank white box. Substitute the Excel variable names for the TauDEM variable names, which are visible in the box labeled Fields. Select OK. You have just assigned every TauDEM stream segment a modeled value for Summer MWMT, based on actual temperature data and geomorphic variables that were derived from TaudDEM and the DEM. 5

6 3. An Application for Modeling Temperature Using ArcMap Take note: There are many reasons to model stream temperature, but here we will describe how to estimate the extent and size (e.g., stream length) of potential aquatic habitat under current conditions. In the Secesh River basin we assume that thermally suitable bull trout habitat is 17.5 C, thus, we wish to identify stream segments in ArcMap that fall within our criteria. 1. In ArcMap, open the TauDEM stream layer 2. Display the stream segments by using different colors for suitable and nonsuitable habitat: a. Right click on the TauDEM shapefile > Properties > Symbology b. In the left box under Show: choose Quantities > Graduate colors c. Fields > Value > select PrdSumMWMT (our modeled stream temperature variable) d. Classification > Classes > 2 > click on Classify 6

7 e. When the Classification dialogue box opens: Classification Method > Manual f. Break Values (red arrow): type 17.5 on the first line and accept the 2 nd line default (Note the distribution of data in the histogram, most of the stream segments in the Secesh River basin are 17.5 C) g. Click OK > OK, again. In this figure, the light blue segments are estimated to be 17.5 C, based on our regression model, and the dark blue stream segments are >17.5 C. 7

8 3. Calculate the stream length of suitable fish habitat (there are several ways to do this, but we show one example below). a. Select the stream segments that are 17.5 C: i. From the toolbar, choose Selection > Select by Attributes 1. Layer: TauDEM 2. Method: Create a new selection 3. Select box: PrdSumMWMT OK All of the segments 17.5 C should be highlighted. b. Right click on the TauDEM shapefile > Open attribute table c. Right click on the variable Length > Statistics The Statistics box displays many useful metrics for Length. The total length of suitable stream segments (based on temperature alone) is 778,513 m. 8

9 Take note: Not all streams 17.5 C are suitable habitat for bull trout. Regardless of temperature, some streams may be too small or too steep to be usable. We filtered our stream layer, again, based on published research that suggests bull trout are typically found in streams where the drainage area is >400 ha and the slope is <15 %. 4. Filter stream network to exclude segments with contributing areas < 400 ha a. Right click on the TauDEM shapefile > Properties > Defintion Query > Query Builder "ds_cont_ha" > 400 AND "Slope" <.15 Click OK > OK The turquoise-colored streams, in the figure below, depict suitable bull trout habitat in the Secesh River basin total stream length is 184,726 m. 9