River channels in gis l.jpg
This presentation is the property of its rightful owner.
Sponsored Links
1 / 64

River Channels in GIS PowerPoint PPT Presentation


  • 303 Views
  • Uploaded on
  • Presentation posted in: General

River Channels in GIS. Venkatesh Merwade, Center for Research in Water Resources, University of Texas at Austin. Overview. Fish Habitat Modeling using GIS Standardized 3D representation of river channels River Channel Morphology Model RCMM and Hydraulic Modeling. Instream flow studies.

Download Presentation

River Channels in GIS

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


River channels in gis l.jpg

River Channels in GIS

Venkatesh Merwade, Center for Research in Water Resources, University of Texas at Austin


Overview l.jpg

Overview

  • Fish Habitat Modeling using GIS

  • Standardized 3D representation of river channels

  • River Channel Morphology Model

  • RCMM and Hydraulic Modeling


Instream flow studies l.jpg

Instream flow studies

  • How do we quantify the impact of changing the naturalized flow of a river on species habitat?

  • How do we set the minimum reservoir releases that would satisfy the instream flow requirement?


Objective l.jpg

Objective

  • Objective

    • To model species habitat as a function of flow conditions and help decision making

  • Instream Flow

    • Flow necessary to maintain habitat in natural channel.


Methodology l.jpg

Methodology

  • Species habitat are dependent on channel hydrodynamics – hydrodynamic modeling

  • Criteria to classify species depending on the conditions in the river channel – biological studies

  • Combine hydrodynamics and biological studies to make decisions – ArcGIS


Fish habitat modeling l.jpg

Criterion

Depth & velocity

Species groups

Habitat

Hydrodynamic

Habitat

Model

Model

Descriptions

RMA2

Biological Sampling

GIS

Instream

Flow

DecisionMaking

Fish Habitat Modeling


Data requirement l.jpg

Data Requirement

  • Hydrodynamic Modeling

    • Bathymetry Data (to define the channel bed)

    • Substrate Materials (to find the roughness)

    • Boundary Conditions (for hydrodynamic model)

    • Calibration Data (to check the model)

  • Biological Studies

    • Fish Sampling (for classification of different species)

    • Velocity and depth at sampling points


Study area guadalupe river near seguin tx l.jpg

Study Area (Guadalupe river near Seguin, TX)

1/2 meter Digital Ortho Photography


Depth sounder echo sounder l.jpg

Depth Sounder (Echo Sounder)

The electronic depth sounder operates in a similar way to radar It sends out an electronic pulse which echoes back from the bed. The echo is timed electronically and transposed into a reading of the depth of water.


Acoustic doppler current profiler l.jpg

Acoustic Doppler Current Profiler

Provides full profiles of water current speed and direction in the ocean, rivers, and lakes. Also used for discharge, scour and river bed topography.


Global positioning system gps l.jpg

Global Positioning System (GPS)

Tells you where you are on the earth!


Final setup l.jpg

Final Setup

GPS Antenna

Computer and power setup

Depth Sounder


Final data view l.jpg

Final Data View


2d hydrodynamic model l.jpg

Surface Water Modeling System

(Environmental Modeling Systems, Inc.)

RMA2

(US Army Corps of Engineers)

2D Hydrodynamic Model

  • SMS (Surface Water Modeling System)

    • RMA2 Interface

  • Input Data

    • Bathymetry Data

    • Substrate Materials

    • Boundary Conditions

    • Calibration Data


Sms mesh l.jpg

SMS mesh

Finite element mesh and bathymetric data


Sms results l.jpg

SMS Results


Biological studies tamu l.jpg

Biological Studies (TAMU)

  • Meso Habitat and Micro Habitat

  • Use Vadas & Orth (1998) criterion for Meso Habitats

  • Electrofishing or seining to collect fish samples for Micro Habitat analysis

  • Sample at several flow rates and seasons

  • Measure Velocity and depth at seining points

  • Statistical analysis to get a table for Micro Habitats classification.


Slide18 l.jpg

Deep Pool

Run

Depth [feet]

Medium Pool

Shallow

Pool

Fast Riffle

Slow Riffle

Mesohabitat Criteria: V, D, V/D, FR

(Vadas & Orth, 1998)


Micro habitat table l.jpg

Micro Habitat Table


Hydraulic and biological data l.jpg

Hydraulic and Biological Data

Attribute Table

Bathymetry Points

Habitat Descriptions


Habitat modeling using arcgis l.jpg

Habitat Modeling using ArcGIS


Results l.jpg

Results


Overview24 l.jpg

Overview

  • Fish Habitat Modeling using GIS

  • Standardized 3D representation of river channels

  • River Channel Morphology Model

  • RCMM and Hydraulic Modeling


Channel bathymetry in hydraulic modeling l.jpg

Channel bathymetry in Hydraulic Modeling

Source: RMA2 reference manual, 2002


Channel representation in arc hydro l.jpg

Channel

Channel Representation in Arc Hydro

River channels are represented as a set of cross-sections and profile-lines in Arc Hydro


Gis database for river channels l.jpg

Thalweg

Cross-sections

ProfileLines

3D Network

GIS database for river channels

Measurement points

Surface

Develop generic ways to create all the channel features from measurement points.


Data analysis l.jpg

Centerline/Thalweg

Cross-sections

ProfileLines

Data analysis

Start with points

Extract all the necessary information

Create surface from points

How can we do this…….


Development of geospatial structure for river channels l.jpg

Development of Geospatial Structure for River Channels

Thought Process:

  • Regular FishNet in ArcGIS provides a network of 3D lines, which are not flow oriented

  • If the data are plotted in a flow-oriented system, the regular FishNet becomes flow-oriented.

  • Flow-oriented coordinate system is useful for getting cross-sections and profile-lines.

Regular FishNet


Geospatial structure for river channels methodology l.jpg

Geospatial Structure for River Channels - Methodology

  • Plot the data in a flow-oriented coordinate system (s,n,z).

  • Interpolate the data to create a surface.

  • Create a FishNet from the interpolated surface.

  • Transform the FishNet to (x,y,z).


Measure in arcgis l.jpg

Measure in ArcGIS

A PolylineMZ can store m and z at each vertex along with x and y coordinates.

112.3213

64.0056

0


S n z coordinate system l.jpg

P

s1

Centerlin

e

s2

n1

(s = 0, n = 0)

n2

Bankline

s

P(n1, s1)

Q(n2, s2)

Q

(s,n,z) coordinate system

  • s-coordinate is the flow length along the river channel

  • n-coordinate is the perpendicular distance from the centerline

  • n-coordinate is negative to the LHS and positive to the RHS of the centerline


Defining a thalweg l.jpg

Defining a Thalweg

Input

Output

Steps 3, 4

Step 2

Steps 5,6,7

Step 8

User defines an arbitrary centerline over the measurement points

Thalweg tool creates a surface using the measurement points

Densify the initial centerline to get more points

Normals are drawn at each vertex of the centerline to locate deepest points

All the deepest points replace the vertices of the old centerline

Final result is a 3D polyline defining the thalweg

Old vertices

New vertices


X y z s n z l.jpg

n

+

o

n

-

y

n

s

n

+

x

s

o

(x,y,z)

n

-

s

(s,n,z)

(x,y,z) (s,n,z)


Spatial interpolation l.jpg

Spatial interpolation

Bathymetry Points

  • IDW

  • EIDW

  • Splines

    • Tension

    • Regularized

  • Kriging

    • Ordinary

    • Anisotropic

Interpolated Raster


Spatial interpolation results l.jpg

Spatial Interpolation Results

Anisotropic kriging gave the least RMSE


Fishnet x y z to s n z l.jpg

y

n

x

s

FishNet (x,y,z) to (s,n,z)

FishNet in (s,n,z) is flow-oriented!


Fishnet comparison l.jpg

FishNet comparison

Hydraulic FishNet

Regular FishNet


Profile lines and cross sections in 3d l.jpg

Profile Lines and Cross Sections in 3D

Bird’s eye view!


Slide40 l.jpg

Priority segments are 100s of miles long

Study area is only few miles long

Instream flow studies in Texas

Results from small studies are extrapolated

Are the results valid?? Can we cross-check??


Overview41 l.jpg

Overview

  • Fish Habitat Modeling using GIS

  • Standardized 3D representation of river channels

  • River Channel Morphology Model

  • RCMM and Hydraulic Modeling


Slide42 l.jpg

Goal

  • Based on the knowledge gained from a detailed dataset collected for a reach of river, develop a model for describing the 3D river channel form at regional scale.


Conceptual model l.jpg

Meandering shape

Thalweg location

Cross-section form

C

C

C

C

B

B

B

B

A

A

A

A

Conceptual Model


Channel bathymetry l.jpg

Channel Bathymetry

=

+

  • Channel bathymetry is complex

  • This research is focused on the deterministic component only

Channel Bathymetry

Deterministic Component

Stochastic Component


R iver c hannel m orphology m odel l.jpg

4

River Channel Morphology Model

  • Get the shape (blue line or DOQ)

  • Using the shape, locate the thalweg

  • Using thalweg location, create cross-sections

  • Network of cross-sections and profile lines

1

2

3


Site1 and site2 on brazos river l.jpg

@ 5 miles

@ 30 miles

Site1 and Site2 on Brazos River

The data for Site 1 and Site 2 are available as (x,y,z) points.


Step 1 normalizing the data l.jpg

nL

nR

0

-

+

Z

P(ni, zi)

d

Zd

w = nL + nR

Step 1: Normalizing the data

For any point P(ni,zi), the normalized coordinates are:

nnew = (ni – nL)/w

znew = (Z – zi)/d

For nL = -15, nR = 35, d = 5, Z=10

P (10, 7.5) becomesPnew(0.5, 0.5)


Normalized data l.jpg

Normalized Data

Original cross-section

Modified cross-section

Depth and width going from zero to unity makes life easier without changing the shape of the original cross-section


Shape characterization through radius of curvature l.jpg

Shape characterization through radius of curvature

  • Ifradius of curvature is small, the thalweg is close to the bank and as it increases the thalweg moves towards the center of the channel.

  • If the channel meanders to left, the center of curvature is to the right hand side of the centerline and vice versa.

  • When the center of curvature is to the right, the radius of curvature is considered positive and vice versa

r1

r3

r2


Step 2 locate thalweg using shape l.jpg

Y = 0.076*log(x) + 1.21

0

0.5

1.0

Y = 0.087*log(x) – 0.32

Step 2: locate thalweg using shape


Thalweg and cross section l.jpg

Thalweg and cross-section

  • Cross-section should have an analytical form to relate it to the thalweg location

  • Many probability density functions (pdf) have shapes similar to the cross-section

  • Beta pdf is found feasible

    • its domain is from zero to one

    • it has only two parameters (a,b)


Step 3 cross sections as beta pdfs l.jpg

Step 3: cross-sections as beta pdfs

beta c/s = (beta1 + beta2) * k

a1=5, b1=2, a2=3, b2=3, factor = 0.5

a1=2, b1=2, a2=3, b2=7, factor = 0.6

Create beta cross-sections for different thalweg locations


Cross sections as beta pdf l.jpg

Cross-sections as Beta pdf

Single pdf

Combination of two pdfs

a1=5, b1=2, a2=3, b2=3, factor = 0.5

Simple, only two parameters, 0 < x < 1

A single pdf has a flat tail, which is undesirable.

The condition of unit area under the pdf makes it difficult to maintain z*< 1.

A combination of two beta pdfs offers flexibility to fit any form of cross-sectional shape.


Hydraulic geometry relationships l.jpg

Hydraulic Geometry Relationships

Hydraulic geometry relationships for Brazos River at Richmond.

Hydraulic geometry relationships are developed at USGS gaging stations.

W, d, and v obtained at the gaging stations are then interpolated to get the corresponding values at other locations.

An ideal scenario would be to have gaging stations both upstream and downstream from the point of interest.


Usgs measurements l.jpg

USGS Measurements

http://waterdata.usgs.gov/nwis/measurements


The final framework l.jpg

The final framework

  • Start with a blue line (s), locate the thalweg (t) using the relationship, t = f(s).

  • Using t, describe cross-sections (c) using the relationship, c(a,b) = f(t).

  • The resulting cross-sections have a unit width and unit depth.

  • Rescale the normalized cross-sections using width and depth (hydraulic geometry)


Results57 l.jpg

Results


Lower brazos in texas l.jpg

Lower Brazos in Texas


3d channel representation l.jpg

3D Channel Representation

Cross-sections

Profile-lines

3D Mesh of cross-sections and profile-lines

Set of Volume objects


Overview60 l.jpg

Overview

  • Fish Habitat Modeling using GIS

  • Standardized 3D representation of river channels

  • River Channel Morphology Model

  • RCMM and Hydraulic Modeling


Rcmm and hydraulic modeling l.jpg

RCMM and Hydraulic Modeling

3D Channel Model

  • Blue line to 3D channel using the shape and hydraulic geometry

  • Interaction with external hydraulic models (HEC-RAS) via XML

Blue Line

XML

HEC-RAS

3D Channel

GIS / Hydraulic Model Data Exchange


Hydraulic model attributes l.jpg

Hydraulic Model Attributes

  • Relationships – ReachHasCrossSection

    HydroID of Reach is ReachID of CrossSections


Ftable l.jpg

FTable

  • Linking of 3D channel and hydraulic model can be used to run hydraulic simulations and create FTable in GIS

  • FTable contains useful information on water surface elevations, velocity, volume, residence times

Hydraulic attributes

Reach identifier

Cross-section identifier


Questions l.jpg

Questions


  • Login