Thursday, January 7, 2010

Voss comments on Sutra and its application in the Shenandoah Valley

Cliff wrote me last night and gave me permission to share his e-mail on the blog. Click the comment below. The USGS wrote cool structural geology solid modeling tools for the Shenandoah Valley project. This underscores that solid modeling is a significant research topic in and of itself. Many of the codes that handle these problems are understandably proprietary, but the result is that they are a bit black-box to the research community.

Wednesday, January 6, 2010

Developers Wanted / Specific Programming Needs

I would like to start a developers group for this program. If you have ideas and/or can make programming contributions, please join the project over at SourceForge. Click on the "comments" below this post to see a write-up of specific programming needs.

Anisotropy and MODFLOW limitations, circa 2003

This article of mine reviews the issue of the effect on flow of structural dip, its related vertical anisotropy, and related limitations in MODFLOW:

Note that it cites this paper on the two choices a MODFLOW modeler has for grid design: the "grid overlay method" and "boundary matching approach."

Jones, N.L., T. J. Budge, A.M. Lemon, and A.K. Zundel. 2002. Generating MODFLOW grids from boundary representation solid models. Ground Water 40, no. 2: 194-200.

Though perhaps not using these names for the grid design, Mary Anderson discussed this in both her modeling courses at the University of Wisconsin, and her textbook.

GeoSolid3D sets up models using the "boundary matching approach."

MODFLOW's LVDA package for generalized anisotropy in a layer

I thought I would pass along some information about a useful, relatively new MODFLOW development, the Layer Variable Direction horizontal Anisotropy (LVDA) package. Modeling the effect of anisotropy on flow has interested me for a long time, so I thought I would share some background on where things were when I wrote a paper on the subject in 2003 (previous post), and what LVDA can do now.

MODFLOW's LVDA package is mainly used for implementing generalized anisotropy in a horizontal plane, i.e. as it would change on a map. Speaking with Ken Kipp and Mary Hill circa 2002, I learned that they (and others) were working on this code for the USGS Death Valley Regional Flow project, where they were interested in simulating the effect of jointing that was changing direction spatially.

At the time of my article (2003), I cited the Kladias and Ruskauf (1997) paper “Implementing spatially variable anisotropy in Modflow” in Ground Water 35, and noted that they, “…outline how to incorporate spatially variable (heterogeneous) anisotropy within the plane of the layer in Modflow. Their modification generalizes the column-to-row anisotropy factor into an array, prescribing a unique value for each cell in the model. This approach does not allow for the orientation of the principal axes to change spatially.” I then mentioned, “The U.S. Geological Survey is currently researching ways to incorporate generalized heterogeneous anisotropy within the plane of the layer in MODFLOW, where both the anisotropic factor and orientation of principal axes can change spatially,” a personal communication at the time from Mary Hill (cited) and Ken Kipp. They were referring to the LVDA package in development at the time.

Kipp's and Hill's (and others) work made it to press as a 2002 OFR. It was not available, at least not printed, at the time of my article, but is now incorporated in MODFLOW 2000

with this reference:
Anderman, Evan R.; Kipp, K. L.; Hill, Mary C.; Valstar, Johan; Neupauer, R. M., 2002. Documentation of the Model-Layer Variable-Direction Horizontal Anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) package. USGS Open-File-Report 02-409.

Generalized anisotropy within a layer is a very exciting development for MODFLOW, however the next two sentences from my 2003 article still apply: “Cross-sectional models constructed within a single layer would thus be able to handle structural dip. However two or three-dimensional, multiple layered models exploiting the boundary-matching method [still will not] handle structural dip until a correction for generalized heterogeneous anisotropy is implemented in the vertical plane.”

Though the LVDA package is a significant advancement, MODFLOW is still not capable of handling dipping anisotropy in a vertical plane when using multiple layers, thus cannot accurately model the effect of dip unless the model is a 2D-single-layer, LVDA, cross sectional model with the layer “flipped up.” As far as I know, Cliff Voss’ Sutra is the only USGS model that can handle generalized anisotropy in 3D.

Welcome to GeoSolid3D

This blog is for open-source research into constructing and using geologic solid models, particularly as they pertain to ground water models, though other applications are certainly open for discussion.

My focus will be on the open-source software code I started at

I will also post "blogs" related to modeling groundwater under conditions of parameter anisotropy, specifically generalized 3D anisotropy, and how these issues relate to solid model construction.