USGS Science⌕ Search
Home/Archive/USGS

USGS · ofr8750

WATEQ4F - a personal computer Fortran translation of the geochemical model WATEQ2 with revised data base

J.W. Ball·D. Kirk Nordstrom·D.W. Zachmann
View original source ↗Read source document ↗Export metadata ↓

Abstract

A FORTRAN 77 version of the PL/1 computer program for the geochemical model WATEQ2, which computes major and trace element speciation and mineral saturation for natural waters has been developed. The code (WATEQ4F) has been adapted to execute on an IBM PC or compatible microcomputer. Two versions of the code are available, one operating with IBM Professional FORTRAN and an 8087 or 89287 numeric coprocessor, and one which operates without a numeric coprocessor using Microsoft FORTRAN 77. The calculation procedure is identical to WATEQ2, which has been installed on many mainframes and minicomputers. Limited data base revisions include the addition of the following ions: AlHS04(++), BaS04, CaHS04(++), FeHS04(++), NaF, SrC03, and SrHCO3(+). This report provides the reactions and references for the data base revisions, instructions for program operation, and an explanation of the input and output files. Attachments contain sample output from three water analyses used as test cases and the complete FORTRAN source listing. U.S. Geological Survey geochemical simulation program PHREEQE and mass balance program BALANCE also have been adapted to execute on an IBM PC or compatible microcomputer with a numeric coprocessor and the IBM Professional FORTRAN compiler.

Explore related subjects

Open-File ReportReport

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J.W. Ball, D. Kirk Nordstrom, D.W. Zachmann. 1987. WATEQ4F - a personal computer Fortran translation of the geochemical model WATEQ2 with revised data base. https://doi.org/10.3133/ofr8750

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

KEEP EXPLORING

Related USGS reports

USGS2026

Floods of June 20–July 6, 2024, in the Big Sioux River, Rock River, Little Sioux River, Ocheyedan River, and Floyd River Basins, northwestern Iowa

Major flooding occurred on June 20–July 6, 2024, in northwestern Iowa affecting the Big Sioux, Rock, Little Sioux, Ocheyedan, and Floyd River Basins. Heavy rain fell in northwestern Iowa, southwestern Minnesota, and southeastern South Dakota on June 20–22, 2024. Parts of northwestern Iowa recorded 2–6 inches of rainfall and localized amounts exceeding 12 inches. A maximum peak-of-record streamflow of 175,000 cubic feet per second at the U.S. Geological Survey streamgage Big Sioux River at Akron, Iowa (06485500), was recorded on June 22, 2024, and had an annual exceedance probability range of 0.2–0.49 percent. High-water marks were measured at four locations along the Big Sioux River between U.S. Interstate 29 at Sioux City, Iowa, upstream to Iowa Highway 10 north of Hawarden, Iowa, a distance of 75.7 river miles. A maximum peak-of-record streamflow of 157,000 cubic feet per second at the U.S. Geological Survey streamgage Rock River near Rock Valley, Iowa (06483500), was recorded on June 22, 2024, and had an annual exceedance probability of less than 0.2 percent. High-water marks were measured at eight locations along the Rock River between County Road B30 east of Hudson, South Dakota, upstream to Iowa Highway 9 at Rock Rapids, Iowa, a distance of 39.3 river miles. A maximum peak-of-record streamflow of 63,000 cubic feet per second at the U.S. Geological Survey streamgage Little Sioux River at Correctionville, Iowa (06606600), was recorded on June 24, 2024, and had an annual exceedance probability range of 0.2–0.49 percent. High-water marks were measured at 11 locations along the Little Sioux River between Iowa Highway 31 west of Correctionville, Iowa, upstream to U.S. Highway 18 north of Spencer, Iowa, a distance of 134.8 river miles. A maximum streamflow of 24,500 cubic feet per second at the U.S. Geological Survey streamgage Ocheyedan River near Spencer, Iowa (06605000), was recorded on June 22, 2024, and had an annual exceedance probability range of 0.2–0.49 percent. High-water marks were measured at three locations along the Ocheyedan River between County Road M38 west of Spencer, Iowa, upstream to U.S. Highway 18 west of Everly, Iowa, a distance of 12.8 river miles. A maximum streamflow of 41,000 cubic feet per second at the U.S. Geological Survey streamgage Floyd River at Alton, Iowa (06600100), was recorded on June 22, 2024, and had an annual exceedance probability range of 1–1.99 percent. High-water marks were measured at six locations along the Floyd River between Iowa Highway 3 at Le Mars, Iowa, upstream to Iowa Highway 10 at Alton, Iowa, a distance of 27.5 river miles. The high-water marks were used to develop flood profiles for the Big Sioux, Rock, Little Sioux, Ocheyedan, and Floyd Rivers.

Iowa, Minnesota, South Dakota↗
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↗