USGS Science⌕ Search
Home/Archive/USGS

USGS · ofr97810

FEMFLOW3D; a finite-element program for the simulation of three-dimensional aquifers; version 1.0

Timothy J. Durbin·Linda D. Bond
View original source ↗Read source document ↗Export metadata ↓

Abstract

This document describes a computer program that simulates three-dimensional ground- water systems using the finite-element method. The program was developed to simulate regional ground-water systems, but it can be applied to small-scale problems as well. This program can be used to simulate both confined and water-table aquifers. The program simulates a linearized three-dimensional free-surface ground-water system with a fixed grid. FEMFLOW3D is applicable to the simulation of various free-surface ground-water systems for which the change in aquifer thickness is small relative to the overall aquifer thickness. The finite-element method provides flexibility in the design of a geometric grid that represents the physical dimensions of an aquifer system. For example, features that can be well represented with a finite-element grid include irregular, random geographic and geologic features; irregular boundaries; and increased detail within localized areas of particular interest within the study area. The structure of the computer program consists of a main program, which serves as a simple driver, and a set of subroutines in which the calculations are performed. The background, mathematical basis, structure, and inputs for each of the subroutines are described in the document where applicable. Each subroutine generally handles (1) a part of the mathematical calculations related to the finite-element method, (2) a specific feature of the hydrologic system, or (3) special features related to the management input or output data. Hydrologic features that can be represented with the program include stream-aquifer interactions, phreatophytic evapotranspiration, highly permeable fault zones, land subsidence, and land-aquifer interactions associated with land-use activities. The program can also represent the primary features associated with complex irrigation systems, such as irrigated agriculture, and can calculate the ground-water recharge that results from these activities. Three boundary conditions, including specified-head boundaries, specified-flux boundaries, and variable-flux boundaries, can be represented with the program. The program also provides a method for identifying aquifer and river-bed parameters that can be used in the calibration of models. This document also includes model validation, source code, and example input and output files. Model validation was performed using four test problems. For each test problem, the results of a model simulation with FEMFLOW3D were compared with either an analytic solution or the results of an independent numerical approach. The source code, written in the ANSI x3.9-1978 FORTRAN standard, and the complete input and output of an example problem are listed in the appendixes.

Explore related subjects

Open-File ReportReport

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Timothy J. Durbin, Linda D. Bond. 1998. FEMFLOW3D; a finite-element program for the simulation of three-dimensional aquifers; version 1.0. https://doi.org/10.3133/ofr97810

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Special Contributing Area Loading Program user’s manual

Information on the Special Contributing Area Loading Program execution and functions are presented in this user’s manual. An appendix presents a potential improvement for the user to consider. The hydrologic routing simulation method to model flow through multiple reservoirs, or sewer system components, is described. The use of Special Contributing Areas is described to run a successful simulation, which includes user input of hydrologic time series of flow components and the necessary formats. Upon completion of a successful Special Contributing Area Loading Program simulation, the program outputs hydrologic time series and a descriptive text file containing the model results for each defined sub-unit, or Special Contributing Area. The output time series contain flows through, and overflows from, the three reservoirs in the series, and the text file contains input and output path locations.

Open-File Report↗
USGS2026

Estimating aftershock risk for entry into earthquake-damaged buildings

We present a simple method to estimate the risk of experiencing strong shaking from aftershocks during entry into earthquake-damaged buildings. We compute wait times until the probability of strong ground shaking from aftershocks reaches a predefined risk threshold; for example, a 0.4 percent probability of experiencing Modified Mercalli Intensity 7 or greater shaking during the planned building entry. We also develop a relation between aftershock probability and the probability of strong shaking, so that users can reference the U.S. Geological Survey aftershock forecast during an ongoing aftershock sequence to determine if the risk threshold has been met. We apply our method to active continental regions (for example, the Western United States), stable continental regions (for example, the Central and Eastern United States), and subduction zones (for example, Cascadia or Alaska).

Open-File Report↗
USGS2026

End-user needs for remote sensing wetlands of the Prairie Pothole Region of North America

The Prairie Pothole Region (PPR) of North America comprises globally important grassland and wetland ecosystems critical for numerous populations of migratory birds. Due to the importance of this region for migratory birds, and particularly waterfowl, and the threats of habitat loss due to intensifying agriculture, there is a mature and diverse system of conservation organizations, agencies, and partnerships that spends hundreds of millions of dollars annually on habitat conservation to support migratory bird populations. Remote sensing can be a powerful tool for observing and evaluating global change at large scales as well as expanding inferences from field studies to the broader landscape with statistical models. However, development and utilization of these tools has lagged behind their demand for several reasons, including concerns over spatial and temporal resolution and accuracy of products; perception of a misalignment with decision-maker needs; technological barriers such as skill sets of conservation professionals, computing resources, data access, and usability. In this report, we summarize the needs of conservation professionals and scientists who use or want to use remote sensing data products to inform science about wetland change and conservation of wetlands in the PPR. We assembled this information through several methods leading up to, during, and following a January 2026 PPR Wetland Remote Sensing Workshop. The workshop included United States and Canadian scientists, conservation professionals, and policy experts. Our goal was to bring together end-users and remote sensing product developers jointly to explore reducing the lag between product development and utilization of products to inform science and conservation. Specifically, we aimed to identify gaps in wetland remote sensing that limit effective monitoring, management, and conservation in the PPR, and to develop a framework that outlines pathways to address these gaps by fostering collaboration, improving communication networks, encouraging discussion, and building on existing and ongoing efforts. This report summarizes our participants’ descriptions of end-user needs and the outcomes of the workshop.

Prairie Pothole region↗