Aquaveo & Water Resources Engineering News

Creating a Confined Aquifer

Does your MODFLOW model contain a confined aquifer? A confined aquifer layer is defined as "an aquifer below the land surface that is saturated with water. Layers of impermeable material are both above and below the aquifer, causing it to be under pressure so that when the aquifer is penetrated by a well, the water will rise above the top of the aquifer."

GSSHA

In MODFLOW, a layer is considered confined when the head in the cell is above the top of the cell. Additionally, any cell located above the water table will be unconfined because the head in those cells will be below the top of the cell.

The layers in your GMS MODFLOW model can be assigned as confined or convertible in any of the flow packages, such as in the LPF package. Other flow packages can be used as well, including the BCF, HUF, and UPW packages. When setting up the MODFLOW model, select the desired flow package in the MODFLOW Packages dialog. The selected flow package must be compatible with the specified MODFLOW version.

In any of the flow package dialogs, under the Layer Type subheading, a layer can be defined as "Confined" or "Convertible". "Convertible" means GMS will automatically assign the layer as confined or unconfined depending on the elevation of the water table in the simulation. Only one layer type can be assigned to each layer. By default, all layers are set to convertible unless specified otherwise.

When a layer is explicitly set to be confined, MODFLOW will use the thickness of the cell, rather than the saturated thickness, to compute a transmissivity value. It will not check for the unconfined condition in the layer.

After defining the layers as confined or convertible, you will have a confined aquifer for your simulation. Try using confined aquifer layers in your GMS models today!

Blog tags: 

Using GSSHA Group Files in WMS

GSSHA Group files in WMS allow multiple GSSHA projects to be saved together in one set of files. It can sometimes be helpful for you to run GSSHA on multiple related projects in order to compare the results.

GSSHA

You can open multiple projects into one WMS session by doing the following:

  1. Open the first project using the GSSHA | Open Project File option.
  2. Open the second, third, ect., project using the same menu.

Once the desired projects are open in WMS, you can save them as a GSSHA Group Project (GGP) file. This keeps everything together in one location.

  1. Save the project as a group by using the GSSHA | Save Group option.
  2. Select the projects to include in the group.
  3. Give the group a file name.

GSSHA can be run on a group file if you do the following:

  1. Select the GSSHA | Run GSSHA Group option.
  2. Select the projects to run.
  3. Give the project a name.
  4. Select whether to suppress screen printing and/or to import the solutions when the GSSHA run finishes.

When GSSHA finishes running on the first project, the Model Wrapper dialog will automatically close. A few moments later, another instance of the Model Wrapper dialog will appear as GSSHA runs on the second project. This one will also automatically close. This continues this way through all of your selected projects. WMS will then import the solutions from the GSSHA runs.

All of the tips here require you to be in the 2D Grid Module—there are additional GSSHA menus in other modules. Additionally, warnings may appear asking if land use and soil type tables should be overwritten. This is up to your discretion. Keep in mind, however, that it is generally best to import projects in an additive manner, with the most simple project being imported first and the most complex project being imported last.

Try out saving and running GSSHA groups in WMS today!

Blog tags: 

How to Relax or Refine a Mesh in SMS

Are you ever tempted to refine your mesh after creation using the mesh editing tools? Manually changing hundreds or thousands of mesh elements using the mesh tools can be very time consuming.

Although possible to use the mesh tools, it is better to recreate a new mesh when wanting to refine or relax the elements. This will maintain the accuracy of your mesh, and reduces possible issues that may arise later on in the project due to inconsistencies that editing manually creates.

To recreate a refined or relaxed mesh:

  1. Create a mesh generator coverage.
    1. If the existing mesh was created from a mesh generator coverage, you can use that coverage or a duplicate of that coverage.
    2. If a mesh generator coverage does not exist, you can create a mesh generator coverage by converting the mesh boundary or nodestrings to a new mesh generator coverage.
  2. Go to your mesh generator coverage and redistribute the vertices so that there are either more or less arc segments depending on what you are trying to accomplish with your model. To relax the mesh, fewer arc segments are desired, and to refine the mesh more arc segments are desired.
  3. Assign the desired meshing type (paving, patch, etc.) to the polygons in the mesh generator coverage and set the bathymetry source to use. Be certain to preview how the mesh will generate with the new spacing for each polygon.
  4. Now generate a new 2D mesh using the edited mesh generator coverage.
  5. Add the new mesh to your model simulation and run the model using the relaxed or refined mesh.

Sometimes you might not want to relax ro refine the entire mesh. To only relax ro refine part of the mesh do the following.

  1. Create a new mesh generator coverage.
  2. Create a polygon in the area of the mesh you wish to relax or refine.
  3. Redistribute the vertices along the arcs in the polygon to so there is a great or fewer number of segments along each arc.
  4. Set the meshing type for the polygon and assign it a bathymetry source. Be certain to preview the mesh for the polygon to be certain you are getting the mesh you want.
  5. Generate a new 2D mesh for the area being refined or relaxed.
  6. Select both 2D meshes and merge them together.
  7. Add the new merged mesh to your model simulation and run the model.

Using these steps when editing a mesh can make a huge difference in your model and save you time. Try out relaxing and refining meshes in SMS today!

Blog tags: 

Assigning Elevations in GMS

GMS offers several options for importing, exporting and manipulating elevation data. With so many options, sometimes it can be confusing when choosing which method to use. If you sometimes struggle with elevation data, you’re in luck, because in this blog post we will be exploring different ways that you can utilize your elevation data to accomplish your goals. Listed below are several ways in which you might be interested in using your elevation data.

  • Assigning elevations in the conceptual model (e.g. to drain nodes)

    Using the Select Points/Nodes Tool , double-click on a point/node such as a drain node. This brings up the Attribute Table dialog. Make sure the Feature Type is set to point/node. Here you can set the type of point/node (such as to a drain or a river) and set the bottom elevation. Things such as river arcs and drain nodes require elevations to run the model in MODFLOW.

  • Interpolating scatter sets/rasters to MODFLOW elevations

    Right-click on the scatter set or raster and select Interpolate To and select MODFLOW layers. This will bring up the Interpolate to MODFLOW Layers dialog. Select the dataset you want to interpolate on the left side of the dialog and the layer you would like to interpolate to on the right-side. With both the dataset and the desired layer selected, click Map. This will add the selection to the Dataset=>MODFLOW data queue. Select Interpolation Options if you want to change the interpolation method. Click OK to exit the dialog and interpolate the scatter set to the layer.

  • Making sure nothing in your conceptual model assigns a polygon elevation that would overwrite the interpolated values

    It is important to note that if you have top and/or bottom elevations assigned as areal properties to a polygon, and you map this coverage to MODFLOW, any scatter points or raster elevations previously interpolated to MODFLOW as the top/bottom elevations that lie within the polygon will be overwritten.

  • Pulling datasets out of MODFLOW (e.g. Layer → 2D Dataset) for manipulation and/or use elsewhere

    Another great feature available in GMS is the ability to pull elevations from a MODFLOW layer to create a 2D dataset. This 2D dataset can then be manipulated and/or used elsewhere for various purposes. This can be done in one of two ways; by using the Layer → 2D Dataset option, or by using the MODFLOW Layers to Scatter option.

  • Layer → 2D Dataset vs. MODFLOW Layers to Scatter (preferred)

    It is possible to create 2D datasets from layer arrays in MODFLOW by going to the Global Options in MODFLOW, and selecting the array you want to create a dataset of (starting heads, top elevations, bottom elevations). In the array dialog box select Layer → 2D Dataset or Grid → 3D Dataset. The dataset that this creates will have ids, i, j, and f values.

    The preferred method for creating datasets from MODFLOW layers is by selecting the layer and then selecting Grid | MODFLOW Layers to 2D Scatter Points. In the MODFLOW Layers → Scatter Points dialog, you can select to create scatter points within a selected coverage and chose the desired coverage. If you use this option, there must be a polygon in the coverage for the points to map to. This dialog also gives many other options that are extremely useful and convenient. When using this method, the 2D scatter set will have x and y coordinates and the f value.

As you can see there are several ways for you to take advantage of the many options available in GMS when working with elevation data. Whether building from a conceptual model, or maybe even building a conceptual model from a MODFLOW simulation, there are many ways to use your elevation data in GMS. Practice using your elevation data in GMS 10.4 today!

Blog tags: 

Converting CAD Data to Feature Objects in WMS

CAD data has been around for a long time and it’s one of the most readily available formats for geological data. WMS can make use of the data in a variety of ways. One of these ways is to convert CAD data into feature objects.

When CAD data for an area is available, DWG and DXF file data can be automatically converted to feature objects in WMS. Lines, points, and polygons in the CAD file can be turned into feature objects on a Map coverage for use in your WMS projects. This is done by doing the following:

  1. Import the CAD data into WMS. It will appear in the Project Explorer.
  2. In the Project Explorer, right-click on the CAD file and select CAD to | Feature Objects.
  3. In the CAD → Feature Objects dialog, select the layers to be converted to feature objects.
  4. Next, in the Clean Options dialog, select options to clean up the feature objects such as removing dangling arcs.
  5. Finally, the Properties dialog will appear letting you designate the type of coverage that will hold the new feature objects.

After you have converted your CAD data to feature objects, there are few items to keep in mind.

  • Typically, CAD data will not designate the stream direction. Make certain to check the direction of arcs generated from CAD data. Use the Reorder Streams command to fix this.
  • Polygons will not be automatically generated from polygons in the CAD data. Any enclosed arcs must be converted to polygons using the Build Polygons command.
  • It may be necessary to use the Clean Options dialog again after conversion.
  • Additional adjustments to the arcs may be necessary using the Map module tools.

Generating feature arcs from CAD data can save you a lot of time and frustration when building your watershed model. Try out converting CAD data into feature objects using WMS today!

Blog tags: 

Moving SRH-2D Material Attributes to Another Project

After setting your SRH-2D material attributes in one SMS project, have you ever wanted to transfer those attributes to another SMS project? Doing this can save you from having to reenter the same material attributes into a separate project. You can even transfer multiple material coverages at once.

There are three project files referred to here: the original project, the project with the saved material coverage attributes, and the third project to which you want to add the saved material attributes. To make sure these instructions are clear, we'll call them Project A, Project B, and Project C, respectively.

To save time on other projects that use the same material attributes, you can use the following steps to export them for later use:

  1. Outside SMS, create a copy (Project B) of your project directory. This prevents you from accidentally damaging the original project (Project A) during the process.
  2. In Project B (the copy), delete any geometric data that isn’t going to be transferred, such as meshes, grids, GIS data, and scatter sets.
  3. Also remove any coverages and material data you do not want transferred.
  4. On the SRH-2D material coverages, remove any polygons and feature objects that you do not want transferred. This is assuming that the project receiving the new material parameters will be assigning the material properties to different polygons.
  5. This is the tricky part. You have to have at least one Map object for a map file to be saved in SMS. Therefore, create a new coverage of any type as a placeholder. Create a single feature point in this placeholder coverage.
  6. Use the Save command to save Project B that only contains the SRH-2D material coverages with your material properties, and one one other coverage with a single feature object.
  7. In another instance of SMS, open Project C.
  8. Use the Open command to open Project B.
  9. SMS should ask if you want to delete this data. If this happens, respond "No". The coverages from Project B (and its material lists and attributes) will be loaded into Project C.
  10. Now to add the material lists from the Project B into the material lists in Project C by right-clicking on the transferred coverage and selecting the Assign Materials to menu command. SMS will bring up the Assign Materials dialog where you can add materials from one coverage to another. You can either add the materials to the list of this coverage or replace the materials list in this coverage.

Now that you know how to transfer material parameters from one project to another, try it out in the SMS today!

Blog tags: 

Incorporating Geology into a MODFLOW Model

Have you created a MODFLOW model and would like to incorporate geological features in between the MODFLOW layers?

The first step is to create the solid to be used. In GMS, solids are representations of stratigraphy used for site characterization and visualization. Solids can be created in any one of three ways:

  1. Convert horizons to a solid.
  2. Convert one or more TINs to a solid.
  3. Manually create a solid by right-clicking in the Project Explorer and selecting the desired type of solid from the New | Solid menu.

Using the first two methods will allow the solid material composition (what the solid is made of) to be automatically interpolated from the horizon or TIN information. Using the third method will allow you to select the material for each solid reated. The third method also requires knowledge of the XYZ coordinates and other attributes, depending on the type of solid being created (cube, cylinder, sphere, or prism).

Once you’ve created your solids, these steps will integrate the solids into the MODFLOW model:

  1. Right-click on the grid and select Classify Material Zones… to bring up the Classify Material Zones dialog.
  2. Select the solids folder you just created.
  3. Select the Classify algorithm you want to use.
  4. Enter the name of the material set being created.

The Centroid algorithm assigns the solid to the cell if it passes over the centroid of the cell. The predominant material algorithm assigns the solid to the cell if it is the predominant material in the cell (the material making up the highest percentage of the cell). This method maintains your grid, while interpolating the materials to the grid, so that you can have multiple materials within a layer.

Try adding solids to your MODFLOW models in GMS today!

Blog tags: 

Working with Rain Gages in GSSHA

Are you needing to add rain gages to your GSSHA model? Using rain gages to define your precipitation in GSSHA is extremely useful due to its ability to spatially model precipitation over a watershed. Combining spatially varying rainfall with the distributed parameters of GSSHA is a great way to create a fairly realistic model for your watershed.

One of the bigger challenges when simulating a storm event is finding reliable data. Although it is important to investigate the most accurate source for your particular watershed, there are websites containing NOAA and GLDAS data for not only the United States, but globally. Having many sources of compiled data all on one site makes Cuahsi’s HydroClient a very useful resource. For those using our software internationally, another possibly useful resource when gathering storm data from the Global Precipitation Climatology Centre. This gives monthly values from 1901-2013, with newer data being added frequently.

Now that you have solid data for your watershed, it is time to define your gages. To use rain gages as your precipitation input:

  1. Create a rain gage coverage.
  2. Create rain gages in their proper locations.
  3. Using the Select tool, double-click on the gage to bring up the Rain Gage Properties dialog.
    1. Set the gage type to GSSHA.
    2. Define the precipitation using either a cumulative or an incremental distribution.
    3. Beneath the Show drop-down, choose GAGES if your data is incremental, and ACCUM if your data is accumulative.
  4. Return to the 2D-Grid module and select GSSHA | Precipitation.
    1. Select Gage as the rainfall event, select Rain Gage, and choose your preferred interpolation method (Inverse distance weighted or Thiessen polygons).

For more information and specifics on working with gages in GSSHA, please take a look at the GSSHA user’s manual. WMS provides a useful and helpful resource when creating a GSSHA model, analyzing and viewing the results. Practice using rain gages as your precipitation source in WMS 11.0 today!

Blog tags: 

New Bridge Scour Tool in SMS 13.0

You are probably aware of the potential destruction caused by bridge scour. Up to 60 percent of bridge failures in the United States from 1950 to 1989 were caused by scour. Bridge scour occurs when the stream bed material around bridge piers is eroded. This can leave the pier unsupported by the stream bed, causing it to collapse.

Recent news has shown bridges of all sizes all around the world collapsing or in danger of collapsing due to scour. Such failures can have a huge negative impact on the economy and also pose a danger to the lives of those who use the bridge regularly. The Federal Highway Administration (FWHA) advised that it is far less expensive to take measures to prevent scour than to replace a bridge that fails due to scour.

Part of the scour mitigation process is evaluating countermeasures that can be taken to prevent such failures. These include onsite surveys and inspections, physical modeling such as that done at the FWHA Hydraulics Laboratory at the Turner-Fairbank Highway Research Center in McLean, Virginia, and computer modeling using tools such as the Surface-water Modeling System (SMS) from Aquaveo.

SMS 13.0 adds powerful new post-processing Bridge Scour tools that can be used to quickly test the scour effect of different pier arrangements, pile sizes and shapes, and vessel impact protection structures such as dolphins and fender rings. The Bridge Scour feature in SMS 13.0 can be used with the FWHA Hydraulic Toolbox to take advantage of its various bridge scour calculators and mapping tools.

One of the most exciting features of Bridge Scour is how it saves time. When creating the contracted section, approach, centerline, bank, abutment toe, and pier arcs, SMS 13.0 automatically assigns the arc type if they are created in order. This translates into saved time and cost, allowing you to be more productive.

Check out the Bridge Scour tool in the SMS 13.0 beta.

Blog tags: 

Troubleshooting MODFLOW

Are you experiencing issues with your MODFLOW simulation? Unable to get your model to converge? Even after properly constructing a model in GMS, you might still find that your model won’t converge or it terminated with an error. Below are some hopefully helpful suggestions on why the model might not be converging and what you can do about it.

To begin, look at what might be causing the convergence issues to occur. The model might have improper aquifer properties which should be reviewed and adjusted if needed.

Another possibility is that there is an unbalanced flow budget. An unbalanced flow budget can manifest itself in two ways. One way is when the inflow is greater than the outflow, then the model can experience sometimes extreme flooding, and the model in turn will not converge. The other way in which there can be an unbalanced flow budget is if the outflow is greater than the inflow. If all the cells in a model are caused to go dry, then the model will not converge. A high outflow may be caused by things such as high conductivity and high pumping rates.

Another possible issue that might cause MODFLOW to have some issues is if you have a specified head for all grid cells in the model. This is because when all cells are Specified Head boundaries, then there is nothing for MODFLOW to compute and the model will terminate with an error.

Some other common issues include: improper initial conditions, improper boundary conditions, wetting and drying issues (as mentioned above) and a highly sensitive model. If the area is known to be highly sensitive, this might cause MODFLOW to not converge due to the speed at which flow can be affected.

Now that we know of some issues which might cause MODFLOW to struggle, we can take a look at three basic troubleshooting steps.

Basic Troubleshooting Steps:

  1. Review the command line output from MODFLOW, and check to see where the issue began to arise. This will enable you to better pinpoint the cause of the error.
  2. Before running MODFLOW, make sure you always run the model checker to see if there are any errors that will prevent convergence. Immediately repair the errors found in the Model Checker. The Model Checker can be found by clicking MODFLOW | Check Simulation…
  3. Look in the MODFLOW output file (*.out) to search for missing values.

Further detail for some specific issues is also given below.

  • If the head is going to drop below the cell, use the MODFLOW-NWT solver.
  • Cells usually go dry due to a low recharge or a high conductivity. Adjust these parameters to better calibrate the model.
  • A transient water table will cause a fluctuation in the heads and the cells may go dry. In this case it is best to use the rewetting option which is available in all flow packages. However when possible MODFLOW-NWT is still the best solution to this problem.
  • Relax the maximum residual or head change criteria. It is best to not increase these values beyond about 1% of the value.

Hopefully with these troubleshooting tips, you can get your MODFLOW simulation up and running in no time. Additional information can be found in the Frequently Asked Questions section of the MODFLOW user manual under the question "My model hasn’t converged. What can I do?" If you are still struggling to get your model to work, consider using Aquaveo's consulting services for expert assistance.

Blog tags: 

Pages