Pages

Showing posts with label Columbia River. Show all posts
Showing posts with label Columbia River. Show all posts

Monday, February 10, 2014

Vertical Datum Adjustments in HEC-RAS

Written by Daniela Todesco, P.E.  |  WEST Consultants

Copyright © The RAS Solution 2014.  All rights reserved.

If you have built enough HEC-RAS models, you’ve probably discovered that having data in different vertical datums is a bit of a pain. RAS lets you adjust the vertical datum of individual cross sections (under “Geometric Data”, “Cross Section”, “Options”, “Adjust Elevations”) or portions of the model (“Geometric Data”, “Tools”, “Datum Adjustment”). However, I have to confess this: datum adjustments are always a little bit puzzling to me. I need pen and paper and a little numbering conversion to make sure that I’m doing things right when moving from one vertical datum to the next. And I recently had a “eureka” moment when I finally realized that the sentence “NAVD88 is higher than NGVD29 (for most of the locations we work on)” is actually wrong. Let me step back a moment and I’ll explain why.

The two most common vertical datums used in the U.S. are the North American Vertical Datum of 1988 (NAVD88) and the National Geodetic Vertical Datum of 1929 (NGVD29). NOAA’s VERTCON program can be used to compute the conversion factors from one datum to the other.

Tidal datums are, in general, confusing. Vertical elevations can refer to NAVD88, the station datum, Mean High Water, Mean Low Water, Mean Range of Tide, and a variety of other datums. In recent years, NOAA has been updating its reference tables with excellent graphs showing how the different datums relate to each other (like the one below for Astoria, OR - http://tidesandcurrents.noaa.gov/datums.html?id=9439040 -, where all elevations presented are related to the station datum, which is the Columbia River Datum or CRD).

clip_image002

Let’s say that you download data in the Astoria station datum and you want to convert the data in NAVD88. Since the NAVD88 datum is 2.02 ft above the station datum, you need to subtract 2.02 ft from the station datum to obtain elevations in NAVD88. Or, to make it easier to remember, if the data you obtain is in datum X and you want to convert it to datum Y, you need to ADD the vertical distance between the two datums if X is ABOVE Y (or SUBTRACT if X is below Y). In our case, X is the Astoria station datum, Y is NAVD88, X is below Y, and you subtract 2.02 ft to the Astoria datum to obtain elevations in NAVD88 (the distance between the two). Makes sense, right? This also explains why for most locations on the West Coast USA, the elevation in NAVD88 of a certain location is always higher than the elevation in NGVD29 of the same location (that’s because the NGVD29 datum is actually above the NAVD88 datum). For example, in Astoria, you need to subtract 3.4 ft to elevations in NAVD88 to obtain elevations in NGVD29.

Now, to confuse things even more, gage locations along the Columbia River follow an intrinsic datum that actually changes as you move upstream from the river’s mouth. For example, at Astoria, the station datum is BELOW NGVD29, while at Vancouver, the station datum is ABOVE NGVD29 (see graph below, at http://www.thsoa.org/hy05/09_1.pdf).

clip_image004

The USGS sent us a nice graph (see below) that helped me clear things up when reading datum-related station description information. This one is for the Columbia River at Vancouver gage in WA. The short station description reads: “Datum of gage: 1.82 ft above NGVD29”. The long version reads: “Datum of gage is Columbia River Datum (add 1.82 ft to correct to NGVD29)”. Remember: X is the station datum, Y is NGVD 29, X is ABOVE Y, so you need to add the distance to go from X to Y. To add something more to this graph, I would say that in CRD, the 0 of the NGVD29 datum is 1.82 ft below the 0 of the CDR datum, which reflects the previous graph. Hope this helps!

clip_image005

Monday, March 25, 2013

Quasi Two-Dimensional Modeling in HEC-RAS

Written by Chris Goodell, P.E., D. WRE | WEST Consultants
Copyright © RASModel.com. 2013.  All rights reserved.

One of the limitations of HEC-RAS is that it is a one-dimensional model. Simply put, RAS assumes all flow moves along a singular dimension. For a given cross section, all of the flow is assumed to move either downstream, or all of it moves upstream, along the singular dimension (which can be defined as a polyline-does not have to be a straight line). The consequence of this is that there is only one water surface elevation (stage), and one total flow for a given time step at a given cross section. All of the other variables for a given cross section that you see in the profile output table, detailed output table, DSS, etc. are derived from the stage and flow values. This includes the velocity and shear stress distributions over a cross section, which can provide the appearance of a 2-dimensional analysis. But that is all based on a conveyance distribution over geometric segments of the cross section using that single water surface elevation and single total flow.

So why do I bring this up? First, it's always good to know ALL of the limitations of whatever model you're using to predict future outcomes. But I also want to demonstrate the "quasi-2-dimensional" capabilities of HEC-RAS. While planning a hydraulic study in an estuarine environment, you may immediately start thinking about which 2-dimensional model you want to use. But I've seen many great (and creative) applications of HEC-RAS in these 2-dimensional environments that produce very reasonable, if not accurate results.  In short, a quasi-2-d analysis in RAS requires you, the user, to understand up front the likely flow patterns in your study area. This is best accomplished by going out to the field and looking at your site, studying topographic and bathymetric maps, looking at aerial photographs, and simple common sense and experience. Once you've determined your perceived flow paths, all water outside of these flow paths should either be ineffective flow areas, storage areas, or even separate reaches.
  
clip_image002
Here’s an example of an estuarine environment on the Oregon Coast (Yaquina Bay). I haven’t modeled this yet, but if I were, here’s how I would approach my model setup:










1. clip_image004 Draw a stream centerline (blue in the figure) that represents the singular dimension of flow movement-i.e. flow will either move downstream or upstream along in the direction of this line. Cut cross sections at an appropriate spacing, making sure to cover all areas that could get wet during the simulation. Yes, the trib channel south of the main reach is not covered, but I’ll get to that in a second.

2. Define ineffective flow areas. These are areas that you will expect WON’T have flow actively moving along the singular dimension. Be sure to appropriately define expansion and contraction of flow as you draw in the ineffective polygons. All portions of your cross sections that fall within these areas should be set to be ineffective in your RAS model.  


3. clip_image006Areas that could possibly have a different water surface elevation than the nearest cross section should be split out and modeled as an off-line storage area. Connect that Storage Area to the main reach using a Lateral Structure. You’ll have to come up with a stage-storage curve for the storage area, to be able to model it in RAS. This is a very easy and straight-forward exercise in GIS, as long as you have sufficient topographic coverage. Keep in mind, RAS uses the simplified level pool routing method for Storage Areas. Lateral Structures used for this application will not have an actual “structure” associated with it, so the discharge coefficient you use is very subjective. Typically values on the order of 0.5 to 1.5 are used. Calibrate this if you can.


4. clip_image008Alternatively, you can model the tributary as its own reach, connected to the main channel with a junction. This will allow you to model it using the full dynamic St. Venant equations, giving a more physically representative answer in the trib. However, if movement of water through this reach is relatively slow (i.e. typical ebb and flood tides), a storage area will be fine-and easier!  You can get as complex as you want. There are no limitations within RAS to the number of storage areas, ineffective flow areas, lateral structures, and tributary reaches you use. Just keep in mind, the more complex you make it, the more difficult it will be to troubleshoot any instabilities.


The following video is a great example of a quasi-2-d application of HEC-RAS. This very complex model and the video were created by Gary Brunner at the Hydrologic Engineering Center.
HEC-RAS model of the Lower Columbia River Estuary-Courtesy Gary Brunner