Advanced Terrain Analysis Concepts

1. Advanced Terrain Analysis Concepts • Specific Catchment Area • The D Surface Flow Model • Topmodel • Terrain based calculations of saturated areas and soil moisture deficit • Generalized flow accumulation functions • Flow algebra • Extending ArcGIS • Scripts and Programming

2. Drainage area can be concentrated or dispersed (specific catchment area) representing concentrated or dispersed flow.

3. Stream line Contour line Upslope contributing area a Specific catchment areaa is the upslope area per unit contour length [m2/m  m]

4. The D Algorithm Tarboton, D. G., (1997), "A New Method for the Determination of Flow Directions and Contributing Areas in Grid Digital Elevation Models," Water Resources Research, 33(2): 309-319.) (http://www.engineering.usu.edu/cee/faculty/dtarb/dinf.pdf)

5. Contributing Area using D Contributing Area using D8

6. TOPMODEL Beven, K., R. Lamb, P. Quinn, R. Romanowicz and J. Freer, (1995), "TOPMODEL," Chapter 18 in Computer Models of Watershed Hydrology, Edited by V. P. Singh, Water Resources Publications, Highlands Ranch, Colorado, p.627-668. “TOPMODEL is not a hydrological modeling package. It is rather a set of conceptual tools that can be used to reproduce the hydrological behaviour of catchments in a distributed or semi-distributed way, in particular the dynamics of surface or subsurface contributing areas.”

7. Hydrological processes within a catchment are complex, involving: • Macropores • Heterogeneity • Fingering flow • Local pockets of saturation The general tendency of water to flow downhill is however subject to macroscale conceptualization

8. Map of saturated areas showing expansion during a single rainstorm. The solid black shows the saturated area at the beginning of the rain; the lightly shaded area is saturated by the end of the storm and is the area over which the water table had risen to the ground surface. [from Dunne and Leopold, 1978] Seasonal variation in pre-storm saturated area [from Dunne and Leopold, 1978]

9. TOPMODEL Key Ideas • Surface saturation and soil moisture deficits based on topography • Slope • Specific Catchment Area • Topographic Convergence • Partial contributing area concept • Saturation from below (Dunne) runoff generation mechanism

10. Topmodel - Assumptions Specific catchment areaa [m2/m  m] (per unit contour length) • The soil profile at each point has a finite capacity to transport water laterally downslope based on slope and transmissivity, qcap=T tan. • The actual lateral discharge is proportional to specific catchment area, q=Ra • Relative wetness at a point and depth to water table is determined by comparing q and qcap • The dynamics of the saturated zone can be approximated by successive steady state representations involving the depth to water table adjusting to accommodate q that results in expansion and contraction of the variable source area contributing to saturation excess runoff tan VSA

11. TOPMODEL and GIS Slope Specific Catchment Area Wetness Index ln(a/tan) from Raster Calculator. Average, l = 6.9

12. Numerical Example • Compute • R=0.0002 m/h • l=6.9 • T=2 m2/hr Given • Ko=10 m/hr • f=5 m-1 • Qb = 0.8 m3/s • a (from GIS) • ne = 0.2 Raster calculator -( [ln(sca/S)] - 6.9)/5+0.46

13. Generalized Flow Accumulation

14. Useful for example to track where sediment or contaminant moves

15. Useful for example to track where a contaminant may come from

16. Useful for a tracking contaminant or compound subject to decay or attenuation

17. Useful for a tracking a contaminant released or partitioned to flow at a fixed threshold concentration

18. Transport limited accumulation Useful for modeling erosion and sediment delivery, the spatial dependence of sediment delivery ratio and contaminant that adheres to sediment

19. Reverse Accumulation Useful for destabilization sensitivity in landslide hazard assessment with Bob Pack

20. Extending ArcGIS via programming • Why Programming • Automation of repetitive tasks • Implementation of functionality not available • Programming functionality • Scripts (AML, VB, Python, Avenue) • Interfaces for application programmers • Model Builder • ArcObjects • COM Integration

21. Three Views of GIS Geodatabase view: Structured data sets that represent geographic information in terms of a generic GIS data model. Geovisualization view: A GIS is a set of intelligent maps and other views that shows features and feature relationships on the earth's surface. "Windows into the database" to support queries, analysis, and editing of the information. Geoprocessing view: Information transformation tools that derives new geographic data sets from existing data sets. adapted from www.esri.com

22. Examples • TauDEM – ArcMap toolbar using Visual Basic/C++ (implementation of functionality not available in ArcGIS) • Visual Basic Programming of simple grid calculations

23. TauDEM Software Functionality • Pit removal (standard flooding approach) • Flow directions and slope • D8 (standard) • D (Tarboton, 1997, WRR 33(2):309) • Flat routing (Garbrecht and Martz, 1997, JOH 193:204) • Drainage area (D8 and D) • Network and watershed delineation • Support area threshold/channel maintenance coefficient (Standard) • Combined area-slope threshold (Montgomery and Dietrich, 1992, Science, 255:826) • Local curvature based (using Peuker and Douglas, 1975, Comput. Graphics Image Proc. 4:375) • Threshold/drainage density selection by stream drop analysis (Tarboton et al., 1991, Hyd. Proc. 5(1):81) • Wetness index and distance to streams • Water Quality Functions

24. TauDEM in ArcGIS Visual Basic ESRI ArcGIS 8.x+9.x Toolbar Visual Basic MapWindow plugin Standalone command line applications C++ COM DLL interface Available from TauDEM C++ library Fortran (legacy) components http://www.engineering.usu.edu/dtarb/ Data Access libraries (grid, shape, image, dbf) ESRI gridio API (Spatial analyst) Vector shape files ASCII text grid Binary direct access grid ESRI binary grid Data formats

25. Implementation Details Spatial Analyst includes a C programming API (Application Programming Interface) that allows you to read and write ESRI grid data sets directly. Excerpt from gioapi.h / * GetWindowCell - Get a cell within the window for a layer, * Client must interpret the type of the output 32 Bit Ptr * to be the type of the layer being read from. * * PutWindowCell - Put a cell within the window for a layer. * Client must ensure that the type of the input 32 Bit Ptr * is the type of the layer being read from. * */ int GetWindowCell(int channel, int rescol, int resrow, CELLTYPE *cell); int PutWindowCell(int channel, int col, int row, CELLTYPE cell);

26. C++ COM Methods used to implement functionality using Microsoft Visual C++ STDMETHODIMP CtkTauDEM::Areadinf(BSTR angfile, BSTR scafile, long x, long y, int doall, BSTR wfile, int usew, int contcheck, long *result) { USES_CONVERSION; //needed to convert from BSTR to Char* or String *result = area( OLE2A(angfile), OLE2A(scafile), x,y,doall, OLE2A(wfile), usew, contcheck); return S_OK; }

27. Visual Basic for the GUI and ArcGIS linkage Private TarDEM As New tkTauDEM … Private Function runareadinf(Optional toadd As Boolean = False) As Boolean Dim i As Long runareadinf = False i = TarDEM.Areadinf(tdfiles.ang, tdfiles.sca, 0, 0, 1, "", 0, 1) If TDerror(i) Then Exit Function If toadd Then AddMap tdfiles.sca, 8 End If runareadinf = True End Function

28. Visual Basic Programming in ArcMAP References ESRI, (1999), ArcObjects Developers Guide: ArcInfo 8, ESRI Press, Redlands, California. Zeiler, M., (2001), Exploring ArcObjects. Vol 1. Applications and Cartography. Vol 2. Geographic Data Management, ESRI, Redlands, CA.

29. AREA 2 3 AREA 1 12 Are there any questions ?