GIS in Water Resources

1. GIS in Water Resources Review for Midterm Exam

2. Data Models • A geographicdata model is a structure for organizing geospatial data so that it can be easily stored and retrieved. Geographic coordinates Tabular attributes

3. Raster and Vector Data Raster data are described by a cell grid, one value per cell Vector Raster Point Line Zone of cells Polygon

4. Workspace Geodatabase Feature Dataset Feature Class Geometric Network Relationship Object Class ArcGIS Geodatabase

5. Geodatabase and Feature Dataset • A geodatabase is a relational database that stores geographic information. • A feature dataset is a collection of feature classes that share the same spatial reference frame.

6. Feature Class • A feature class is a collection of geographic objects in tabular format that have the same behavior and the same attributes. Feature Class = Object class + spatial coordinates

7. Object Class • An object class is a collection of objects in tabular format that have the same behavior and the same attributes. An object class is a table that has a unique identifier (ObjectID) for each record

8. Relationship Relationship between spatial and non-spatial objects Water quality data (non-spatial) Measurement station (spatial)

9. National Hydro Data Programs National Elevation Dataset (NED) National Hydrography Dataset (NHD) NED-Hydrology Watershed Boundary Dataset

10. 1:250,000 Scale Soil Information http://www.ftw.nrcs.usda.gov/stat_data.html

11. National Land Cover Dataset http://landcover.usgs.gov/nationallandcover.html Get the data: http://seamless.usgs.gov/

12. National Water Information System Web access to USGS water resources data in real time http://waterdata.usgs.gov/usa/nwis/

13. Drainage System Hydro Network Flow Time Time Series Hydrography Channel System Arc Hydro Components GIS provides for synthesis of geospatial data with different formats

14. “Burning In” the Streams Synthesis of Raster and Vector data  Take a mapped stream network and a DEM Make a grid of the streams Raise the off-stream DEM cells by an arbitrary elevation increment  Produces "burned in" DEM streams = mapped streams = +

15. AGREE Elevation Grid Modification Methodology

16. Geodesy, Map Projections and Coordinate Systems • Geodesy - the shape of the earth and definition of earth datums • Map Projection - the transformation of a curved earth to a flat map • Coordinate systems - (x,y) coordinate systems for map data

17. Latitude and Longitude in North America Austin: (30°N, 98°W) Logan: (42°N, 112°W) 60 N 30 N 60 W 120 W 90 W 0 N

18. Length on Meridians and Parallels (Lat, Long) = (f, l) Length on a Meridian: AB = ReDf (same for all latitudes) R Dl D R 30 N C B Re Df 0 N Re Length on a Parallel: CD = R Dl = ReDl Cos f (varies with latitude) A

19. Example: What is the length of a 1º increment along • on a meridian and on a parallel at 30N, 90W? • Radius of the earth = 6370 km. • Solution: • A 1º angle has first to be converted to radians • p radians = 180 º, so 1º = p/180 = 3.1416/180 = 0.0175 radians • For the meridian, DL = ReDf = 6370 * 0.0175 = 111 km • For the parallel, DL = ReDl Cos f • = 6370 * 0.0175 * Cos 30 • = 96.5 km • Parallels converge as poles are approached

20. Horizontal Earth Datums • An earth datum is defined by an ellipse and an axis of rotation • NAD27 (North American Datum of 1927) uses the Clarke (1866) ellipsoid on a non geocentric axis of rotation • NAD83 (NAD,1983) uses the GRS80 ellipsoid on a geocentric axis of rotation • WGS84 (World Geodetic System of 1984) uses GRS80, almost the same as NAD83

21. Vertical Earth Datums • A vertical datum defines elevation, z • NGVD29 (National Geodetic Vertical Datum of 1929) • NAVD88 (North American Vertical Datum of 1988) • takes into account a map of gravity anomalies between the ellipsoid and the geoid

22. Coordinate System A planar coordinate system is defined by a pair of orthogonal (x,y) axes drawn through an origin Y X Origin (xo,yo) (fo,lo)

23. Universal Transverse Mercator • Uses the Transverse Mercator projection • Each zone has a Central Meridian(lo), zones are 6° wide, and go from pole to pole • 60 zones cover the earth from East to West • Reference Latitude (fo), is the equator • (Xshift, Yshift) = (xo,yo) = (500000, 0) in the Northern Hemisphere, units are meters

24. UTM Zone 14 -99° -102° -96° 6° Origin Equator -90 ° -120° -60 °

25. Defined for a feature dataset in ArcCatalog Coordinate System Projected Geographic X/Y Domain Z Domain M Domain ArcInfo 8 Reference Frames

26. X/Y Domain (Max X, Max Y) Long integer max value of 231 = 2,147,483,645 (Min X, Min Y) Maximum resolution of a point = Map Units / Precision e.g. map units = meters, precision = 1000, then maximum resolution = 1 meter/1000 = 1 mm on the ground

27. Four Points

28. One degree box and its four lines Geographic Coordinates

29. One Degree Box in USGS Albers Projection

30. USGS Albers Projection

31. Area Calculation in USGS Albers 81.09 km 111.79 km 111.79 km Area = 9130.6 km2 82.26 km 82.26 + 81.09 x 111.79 = 9130.5 km2 2

32. North American Albers Projection

33. Area Calculation in North American Albers 76.64 km 118.17 km 118.17 km Area = 9130.6 km2 77.89 km 77.89 + 76.64 X 118.17 = 9130.4 2 Take home message: Lengths of lines change but area is constant in Albers

34. (x1,y1) y f(x,y) x Two fundamental ways of representing geography are discrete objects and fields. The discrete object view represents the real world as objects with well defined boundaries in empty space. Points Lines Polygons The field view represents the real world as a finite number of variables, each one defined at each possible position. Continuous surface

35. Vector and Raster Representation of Spatial Fields Vector Raster

36. Numerical representation of a spatial surface (field) Grid TIN Contour and flowline

37. Number of columns Cell size Number of rows Grid Datasets • Cellular-based data structure composed of square cells of equal size arranged in rows and columns. • The grid cell size and extent (number of rows and columns), as well as the value at each cell have to be stored as part of the grid definition.

38. Raster Sampling from Michael F. Goodchild. (1997) Rasters, NCGIA Core Curriculum in GIScience, http://www.ncgia.ucsb.edu/giscc/units/u055/u055.html, posted October 23, 1997

39. Extent Spacing The scale triplet Support From: Blöschl, G., (1996), Scale and Scaling in Hydrology, Habilitationsschrift, Weiner Mitteilungen Wasser Abwasser Gewasser, Wien, 346 p.

40. Raster Generalization Largest share rule Central point rule

41. Raster calculation – some subtleties Resampling or interpolation (and reprojection) of inputs to target extent, cell size, and projection within region defined by analysis mask + = Analysis mask Analysis cell size Analysis extent

42. Interpolation Estimate values between known values. A set of spatial analyst functions that predict values for a surface from a limited number of sample points creating a continuous raster. Apparent improvement in resolution may not be justified

43. Topographic Slope • Defined or represented by one of the following • Surface derivative z • Vector with x and y components • Vector with magnitude (slope) and direction (aspect)

44. 67 56 49 52 48 37 58 55 22 Hydrologic Slope - Direction of Steepest Descent 30 30 67 56 49 52 48 37 58 55 22 Slope:

45. 32 64 128 16 1 8 4 2 Eight Direction Pour Point Model Water flows in the direction of steepest descent

46. Filling in the Pits • DEM creation results in artificial pits in the landscape • A pit is a set of one or more cells which has no downstream cells around it • Unless these pits are filled they become sinks and isolate portions of the watershed • Pit filling is first thing done with a DEM

47. 32 64 128 16 1 8 4 2 Flow Direction Grid

48. Stream cell Cell to Cell Grid Network Through the Landscape

49. 1 1 1 1 1 1 1 1 1 1 1 4 3 3 1 4 3 1 1 3 1 1 2 1 1 12 1 2 1 12 1 1 1 2 16 1 1 2 1 16 2 1 3 6 25 3 6 1 2 25 Contributing Area Grid Drainage area threshold > 5 Cells

50. Delineation ofStreams and Watersheds on a DEM