Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Computers and Geology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67Linked to original sources

Clear-Water Contraction Scour at Selected Bridge Sites in the Black Prairie Belt of the Coastal Plain in Alabama, 2006

The U.S. Geological Survey, in cooperation with the Alabama Department of Transportation, made observations of clear-water contraction scour at 25 bridge sites in the Black Prairie Belt of the Coastal Plain of Alabama. These bridge sites consisted of 54 hydraulic structures, of which 37 have measurable scour holes. Observed scour depths ranged from 1.4 to 10.4 feet. Theoretical clear-water contraction-scour depths were computed for each bridge and compared with observed scour. This comparison showed that theoretical scour depths, in general, exceeded the observed scour depths by about 475 percent. Variables determined to be important in developing scour in laboratory studies along with several other hydraulic variables were investigated to understand their influence within the Alabama field data. The strongest explanatory variables for clear-water contraction scour were channel-contraction ratio and velocity index. Envelope curves were developed relating both of these explanatory variables to observed scour. These envelope curves provide useful tools for assessing reasonable ranges of scour depth in the Black Prairie Belt of Alabama.

Scientific Investigations Report↗

Flood-inundation maps for Cayuga Inlet, Sixmile Creek, Cascadilla Creek, and Fall Creek at Ithaca, New York

Digital flood-inundation maps for a 2.9-square-mile area of Ithaca, New York, were created in 2015–18 by the U.S. Geological Survey in cooperation with the City of Ithaca, New York, and the New York State Department of State. The flood-inundation maps depict estimates of the maximum areal extent and depth of flooding corresponding to selected flood frequencies for Cayuga Inlet, Sixmile Creek, Cascadilla Creek, and Fall Creek and selected water-surface elevations of Cayuga Lake. Flood profiles for the stream reaches were computed by combining a one-dimensional step-backwater model for the stream channels and a two-dimensional model for the overbank areas. The resulting hydraulic model was calibrated by using water-surface profiles from five observed storm events. The model was then used to compute 15 water-surface profiles for 5 flood frequencies (50-, 10-, 2-, 1-, and 0.2-percent annual exceedance probabilities, or 2-, 10-, 50-, 100-, and 500-year recurrence intervals) and 3 lake levels (representing average conditions, a 2-year-high condition, and a 100-year-high condition). The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having 0.31‑foot vertical accuracy and 3.3-foot horizontal resolution) to delineate the maximum area flooded at each water level. Flood-inundation maps and geographic information system flood-extent polygons and depth grids are available in the data release associated with this report. These maps can provide emergency management personnel and residents with information that is critical for flood-management planning, flood-response activities, and postflood recovery efforts.

New York↗

Saline water in southeastern New Mexico

Saline waters from formations of several geologic ages are being studied in a seven-county area in southeastern New Mexico and western Texas, where more than 30,000 oil and gas tests have been drilled in the past 40 years. This area of 7,500 sq. miles, which is stratigraphically complex, includes the northern and eastern margins of the Delaware Basin between the Guadalupe and Glass Mountains. Chloride-ion concentrations in water produced from rocks of various ages and depths have been mapped in Lea County, New Mexico, using machine map-plotting techniques and trend analyses. Anomalously low chloride concentrations (1,000–3,000 mg/l) were found along the western margin of the Central Basin platform in the San Andres and Capitan Limestone Formations of Permian age. These low chloride-ion concentrations may be due to preferential circulation of ground water through the more porous and permeable rocks. Data being used in the study were obtained principally from oil companies and from related service companies. The P.B.W.D.S. (Permian Basin Well Data System) scout-record magnetic-tape file was used as a framework in all computer operations. Shallow or non-oil-field water analyses acquired from state, municipal, or federal agencies were added to these data utilizing P.B.W.D.S.-compatible reference numbers and decimal latitude-longitude coordinates. Approximately 20,000 water analyses collected from over 65 sources were coded, recorded on punch cards and stored on magnetic tape for computer operations. Extensive manual and computer error checks for duplication and accuracy were made to eliminate data errors resulting from poorly located or identified samples; non-representative or contaminated samples; mistakes in coding, reproducing or key-punching; laboratory errors; and inconsistent reporting. The original 20,000 analyses considered were reduced to 6,000 representative analyses which are being used in the saline water studies.

New Mexico↗

Practical application of time-lapse camera imagery to develop water-level data for three hydrologic monitoring sites in Wisconsin during water year 2020

Using camera imagery to measure water level (camera-stage) is a well-researched area of study. Previous camera-stage studies have shown promising results when implementing this technology with tight constraints on test conditions. However, there is a need for a more comprehensive evaluation of the extensibility of camera-stage to practical applications. Therefore, the aim of this study was to test a camera-stage method under a wide variety of test conditions to better understand the successes and challenges of using this technology in real-world scenarios. In this study, this approach was tested during Water Year 2020 at three existing U.S. Geological Study (USGS) stream gaging stations in south central Wisconsin that had existing USGS water-level instrumentation. The specific reference objects tested were white pipes and a concrete wall. Since successful application of camera-stage relies on use of suitable images, all captured images in this study were visually inspected to determine suitability for application of camera-stage. Camera-stage measurements were then computed only on images deemed suitable and the results were compared with ground-truth stage values to determine the accuracy. For the purposes of this study, camera-stage values within ±0.10 ft of the actual stage were considered acceptable. One major challenge highlighted was the potential difficulty in obtaining suitable imagery, with the proportion of suitable images varying greatly between the four trials from 38 % to 92 %. The results from applying camera-stage to suitable images were encouraging though, with 79 % to 99 % of evaluated camera-stage values qualifying as acceptable among the four test trials.

Wisconsin↗

Velocity modeling of supercritical pore fluids through porous media under reservoir conditions with applications for petroleum secondary migration and carbon sequestration plumes

Computational methods to characterize secondary migration in porous media traditionally rely on fluid transport equations with assumptions of time invariance, such as flowpath modeling of buoyancy vectors, statistical percolation algorithms, capillary pressure curves, or a form of Darcy’s Law which presumes instantaneous fluid transport. However, in petroleum systems modeling, the timeframe of secondary migration from source to reservoir is important to quantify in relation to other geologic factors such as timing of petroleum generation, fault movement, and seal formation. Additionally, quantifying migration velocities enables an estimation of the distance a plume of geologically sequestered carbon dioxide travels over time, as well as the identification of low-permeability strata appropriate for long-term containment. This study introduces a method to quantify transport velocities of supercritical fluids in low-permeability lithologies for a broad range of rock and fluid properties likely encountered in the sedimentary sequence. A time-dependent form of Darcy’s Law for pressure-driven viscous flow through homogeneous isotropic porous media was used to model flow velocities within a carrier bed. Thermodynamic equations of state were used to determine thermophysical properties of supercritical pore fluids under reservoir pressures ranging from 0–200 MPa (0–29,000 psi) to constrain the momentum equations. Three case studies were examined that (1) estimated fluid flow velocities of methane within the low-permeability Upper Jurassic Haynesville Formation, (2) defined permeability-based flow units to evaluate saline formations for long-term geologic carbon sequestration, and (3) calculated the migration distance of carbon dioxide plumes at the Decatur, Illinois injection and sequestration project.

SEG-AAPG Interpretation↗

Virginia flow-ecology modeling results—An initial assessment of flow reduction effects on aquatic biota

Background The U.S. Geological Survey (USGS), in cooperation with the Virginia Department of Environmental Quality (DEQ), reviewed a previously compiled set of linear regression models to assess their utility in defining the response of the aquatic biological community to streamflow depletion. As part of the 2012 Virginia Healthy Watersheds Initiative (HWI) study conducted by Tetra Tech, Inc., for the U.S. Environmental Protection Agency (EPA) and Virginia DEQ, a database with computed values of 72 hydrologic metrics, or indicators of hydrologic alteration (IHA), 37 fish metrics, and 64 benthic invertebrate metrics was compiled and quality assured. Hydrologic alteration was represented by simulation of streamflow record for a pre-water-withdrawal condition (baseline) without dams or developed land, compared to the simulated recent-flow condition (2008 withdrawal simulation) including dams and altered landscape to calculate a percent alteration of flow. Biological samples representing the existing populations represent a range of alteration in the biological community today. For this study, all 72 IHA metrics, which included more than 7,272 linear regression models, were considered. This extensive dataset provided the opportunity for hypothesis testing and prioritization of flow-ecology relations that have the potential to explain the effect(s) of hydrologic alteration on biological metrics in Virginia streams.

Virginia↗

Time-series model, statistical methods, and software documentation for R–QWTREND—An R package for analyzing trends in stream-water quality

As part of a U.S. Geological Survey water-quality study started in 2018, in cooperation with the International Joint Commission, North Dakota Department of Environmental Quality, and Minnesota Pollution Control Agency, a publicly available software package called R–QWTREND was developed for analyzing trends in stream-water quality. The R–QWTREND package is a collection of functions written in R, an open source language and a general environment for statistical computing and graphics. The package uses a parametric time-series model to express logarithmically transformed concentration in terms of flow-related variability, trend, and serially correlated model errors. Flow-related variability captures natural variability in concentration on the basis of concurrent and antecedent streamflow. The trend identifies systematic changes in concentration in terms of potential step trends, piecewise monotonic trends, or user-specified trends. Maximum likelihood estimation is used to estimate model parameters and determine the best-fit trend model. This report describes the time-series model and statistical methodology behind R–QWTREND and provides formal documentation for installing and using the package.

Open-File Report↗

Occurrence and status of volatile organic compounds in ground water from rural, untreated, self-supplied domestic wells in the United States, 1986-99

Samples of untreated ground water from 1,926 rural, self-supplied domestic wells were analyzed for volatile organic compounds (VOCs) during 1986-99. This information was used to characterize the occurrence and status of VOCs in domestic well water. The samples were either collected as part of the U.S. Geological Survey?s National Water-Quality Assessment (NAWQA) Program occurrence-assessment studies or were compiled by NAWQA from existing ambient ground-water or source-water-quality monitoring programs conducted by local, State, and other Federal agencies. Water samples were collected at the wellhead prior to treatment or storage. In most samples, 55 target VOCs were analyzed, and occurrence and status information generally was computed at an assessment level of 0.2 mg/L (microgram per liter). At least one VOC was detected in 12 percent of samples (232 samples) at an assessment level of 0.2 mg/L. This detection frequency is relatively low compared to the 26 percent detection frequency of at least one VOC in public sup-ply wells sampled by NAWQA, and the difference may be due, in part, to the higher pumping rates, pumping stress factors, and larger contributing areas of public supply wells. Samples with detections of at least one VOC were collected from wells located in 31 of 39 States. Solvents were the most frequently detected VOC group with detections in 4.6 percent of samples (89 samples) at an assessment level of 0.2 mg/L. The geographic distribution of detections of some VOC groups, such as fumigants and oxygenates, relates to the use pattern of com-pounds in that group. With the exception of com-pounds used in organic synthesis, detection frequencies of VOCs by group are proportional to the average half-life of compounds in the group. When the organic synthesis group is excluded from the analysis, a good correlation exists between the detection frequency of VOCs by group and average half-life of compounds in the group. Individually, VOCs were not commonly detected at an assessment level of 0.2 mg/L, with the seven most frequently detected VOCs found in only 1 to 5 percent of samples. Mixtures (two or more compounds) were a common mode of occurrence for VOCs when no assessment level was applied, and mixtures occurred in one-half of all samples that contained at least one VOC. Only 1.4 percent of samples (27 samples) had one or more VOC concentrations that exceeded a federally established drinking-water standard or health criterion. Only 0.1 percent of samples (2 samples) had one or more VOC concentrations that exceeded a taste/odor threshold. Potential point sources of VOCs near domestic wells are numerous. Leaks from under-ground storage tanks and aboveground storage tanks that hold gasoline, diesel fuel, or heating oil have the potential to be major point sources of contaminants to domestic wells. Shock chlorination may be a source of trichloromethane and other trihalomethanes in some domestic wells. Septic systems are believed to be an important source of contaminants to domestic wells, but extensive research on this subject does not exist. VOCs frequently are ingredients in household products such as cleansers and insecticides, and some VOCs have been found in septic systems.

Water-Resources Investigations Report↗

Importance of record length with respect to estimating the 1-percent chance flood

U.S. Geological Survey (USGS) streamflow gages have been established in every State in the Nation, Puerto Rico, and the Trust Territory of the Pacific Islands. From these st reamflow records, estimates of the magnitude and frequency of floods are often developed and used to design transportation and water- conveyance structures to protect lives and property, and to determine flood-insurance rates. Probably the most recognizable flood statistic computed from USGS stream gaging records is the 1- percent (%) chance flood; better known has the 100-year flood. By definition, this is a flood that has a 1% chance of occurring in any given year. The 1% chance flood is a statistical estimate that can be significantly influenced by length of record and extreme flood events captured in that record. Consequently, it is typically recommended that flood statistics be updated on some regular interval such as every 10 years. This paper examines the influence of record length on the 1% chance flood for the Broad River in Georgia and the substantial difference that can occur in the estimate based on record length and the hydrologic conditions under which that record was collected.

Conference Paper↗

The relation between fresh and salty ground water in southern Nassau and southeastern Queens counties, Long Island, New York

An investigation of the geology and ground - water conditions in a part of southwestern Long Island , including the construction of 8 deep and 5 shallow test and observation wells, has been completed by the Geological Survey in cooperation with the Nassau County Department of Public Works and the New York State Water Power and Control Commission. In this area are four main water -bearing units of Late Cretaceous and Pleistocene age. Ground water having a higher than normal salt content is encountered in some wells that penetrate the upper three units. Although the presence of this water constitutes a potential threat, existing fresh - water supplies apparently are in little immediate danger under present conditions. The " salty " ground water has chloride concentrations ranging from about 40 to 18,000 ppm (parts per million), generally more than about 2,000 ppm. It occurs in several bodies, seaward of and beneath the fresh water in the respective aquifers, and is in contact with the fresh water through zones of diffusion. The largest salt- water body, confined under artesian pressure, is in the Jameco gravel and the Magothy (?) formation beneath extreme southwestern Nassau County and adjoining southeastern Queens County . This body is thought to have occupied a substantial part of its present extent since at least 1895. It probably has been and is now encroaching, very slowly, because of increased pumping from the fresh - water bodies and possibly as a result of a postglacial rise of sea level. The rate of encroachment under present conditions probably is less than 100 feet per year. Application of the Ghyben-Herzberg principle to compute depths to the contact between salty and fresh ground water was found to give erroneous results. More accurate depths and a clearer understanding of the relation between heads in adjoining bodies of fresh and salty ground water are obtained by use of a formula adapted from one given by M. King Hubbert.

New York↗

Annual estimates of water and solute export from 42 tributaries to the Yukon River

Annual export of 11 major and trace solutes for the Yukon River is found to be accurately determined based on summing 42 tributary contributions. These findings provide the first published estimates of tributary specific distribution of solutes within the Yukon River basin. First, we show that annual discharge of the Yukon River can be computed by summing calculated annual discharges from 42 tributaries. Annual discharge for the tributaries is calculated from the basin area and average annual precipitation over that area using a previously published regional regression equation. Based on tributary inputs, we estimate an average annual discharge for the Yukon River of 210 km 3  year –1 . This value is within 1% of the average measured annual discharge at the U.S. Geological Survey gaging station near the river terminus at Pilot Station, AK, for water years 2001 through 2005. Next, annual loads for 11 solutes are determined by combining annual discharge with point measurements of solute concentrations in tributary river water. Based on the sum of solutes in tributary water, we find that the Yukon River discharges approximately 33 million metric tons of dissolved solids each year at Pilot Station. Discharged solutes are dominated by cations calcium and magnesium (5.65 × 10 9 and 1.42 × 10 9  g year –1 ) and anions bicarbonate and sulphate (17.3 × 10 9 and 5.40 × 10 9  g year –1 ). These loads compare well with loads calculated independently at the three continuous gaging stations along the Yukon River. These findings show how annual solute yields vary throughout a major subarctic river basin and that accurate estimates of total river export can be determined from calculated tributary contributions.

Hydrological Processes↗

Predicting river travel time from hydraulic characteristics

Predicting the effect of a pollutant spill on downstream water quality is primarily dependent on the water velocity, longitudinal mixing, and chemical/physical reactions. Of these, velocity is the most important and difficult to predict. This paper provides guidance on extrapolating travel-time information from one within bank discharge to another. In many cases, a time series of discharge (such as provided by a U.S. Geological Survey stream gauge) will provide an excellent basis for this extrapolation. Otherwise, the accuracy of a travel time extrapolation based on a resistance equation can be greatly improved by assuming the total flow area is composed of two parts, an active and an inactive area. For 60 reaches of 12 rivers with slopes greater than about 0.0002, travel times could be predicted to within about 10% by computing the active flow area using the Manning equation with n = 0.035 and assuming a constant inactive area for each reach. The predicted travel times were not very sensitive to the assumed values of bed slope or channel width.

Journal of Hydraulic Engineering↗

Landsat data availability from the EROS Data Center and status of future plans

The Department of Interior's EROS Data Center, managed by the U.S. Geological Survey, was established in 1972, in Sioux Falls, South Dakota, to serve as a principal dissemination facility for Landsat and other remotely Sensed data. Through the middle of 1977, the Center has supplied approximately 1.7 million copies of images from the more than 5 million images of the Earth's surface archived at the Center. Landsat accounted for half of these images plus approximately 5,800 computer-compatible tapes of Landsat data were also supplied to users. New methods for processing data products to make them more useful are being developed, and new accession aids for determining data availability are being placed in operation. The Center also provides assistance and training to resource specialists and land managers in the use of Landsat and other remotely sensed data. A Data Analysis Laboratory is operated at the Center to provide both digital and analog multispectral/multitemporal image analysis capabilities in support of the training and assistance programs. In addition to conventionally processed data products, radiometrically enhanced Landsat imagery are now available from the Center in limited quantities. In mid-1978, the Center will convert to an all-digital processing system for Landsat data that will provide improved products for user analysis in production quantities. The Department of Interior and NASA are currently studying concepts that use communication satellites to relay Landsat data between U.S. ground stations, Goddard Space Flight Center and the EROS Data Center which would improve the timeliness of data availability. The Data Center also works closely with the remote sensing programs and Landsat data receiving and processing facilities being developed in foreign countries.

Conference Paper↗

Obtaining valid geologic models from 3-D resistivity inversion of magnetotelluric data at Pahute Mesa, Nevada

We summarize the results of a three-dimensional (3-D) resistivity inversion simulation that we conducted with the intent of characterizing the subsurface 3-D distribution of volcanic composite units of Pahute Mesa, Nevada, without any a priori information on the actual 3-D distribution of the known subsurface geology. The 3-D methodology involved using a 3-D geologic model based on drillhole data and average electrical resistivities of the key hydrostratigraphic units at Pahute Mesa to create a 3-D resistivity forward (“known”) model that depicted the subsurface resistivity structure expected for the input geologic configuration. The calculated magnetotelluric response of the modeled resistivity structure was then assumed to represent observed magnetotelluric data and was used as input into a 3-D resistivity inverse model that was allowed to iteratively estimate in 3-D without any a priori geologic information, in particular, the thickness and resistivity of the volcanic composite units. The resulting 3-D resistivity inversion simulation was compared to the “known” model and the results evaluated. The 3-D inversion was generally able to reproduce the gross resistivity structure of the “known” model, but the simulated conductive volcanic composite unit horizons were often too shallow when compared to the “known” model. Additionally, the chosen computation parameters such as station spacing appear to have resulted in computational artifacts that are difficult to interpret but could potentially be removed with further refinements of the 3-D resistivity inversion modeling technique.

Nevada↗

Evaluation and application of the Purge Analyzer Tool (PAT) to determine in-well flow and purge criteria for sampling monitoring wells at the Stringfellow Superfund site in Jurupa Valley, California, in 2017

The U.S. Geological Survey and U.S. Environmental Protection Agency are developing analytical tools to assess the representativeness of groundwater samples from fractured-rock aquifers. As part of this effort, monitoring wells from the Stringfellow Superfund site in Jurupa Valley in Riverside County, California, approximately 50 miles east of Los Angeles, were field tested to collect information to assist in the evaluation and application of in-well flow as computed by the analytical model called the Purge Analyzer Tool, which computes in-well groundwater travel times for simple piston transport of inflowing groundwater from open intervals of a monitoring well to the pump intake and can provide insight into optimal purging parameters (duration, rate, and pump position) needed for the collection of representative groundwater samples. Field testing of wells included hydraulic, chemistry, and dye tracer analysis to investigate travel times in wells under pumping conditions. The Purge Analyzer Tool was able to replicate dye velocities (travel times) for one of three wells that had appreciable inflow from the aquifer but not the other two wells, which are screened in low-permeability sediments and rock, where flow was dominated by borehole storage. A set of criteria was established to help assess the ability to collect representative groundwater chemistry from monitoring wells; criteria included understanding the height of the static well water column and relative exchange rate between the aquifer and the well.

California↗

Comparison of machine learning approaches used to identify the drivers of Bakken oil well productivity

Geologists and petroleum engineers have struggled to identify the mechanisms that drive productivity in horizontal hydraulically fractured oil wells. The machine learning algorithms of Random Forest (RF), gradient boosting trees (GBT) and extreme gradient boosting (XGBoost) were applied to a dataset containing 7311 horizontal hydraulically fractured wells drilled into the middle member of the Bakken Formation from 2010 through 2017. The initial goal is to use these data‐driven machine learning algorithms to identify the most important explanatory predictors of well productivity within nine subareas and the composite area. Predictor variables representing initial gas production, the initial 180‐day water cut, and vertical depth vary spatially and are identified with geologically favorable areas. Well‐completion predictors include the well lateral length, number of fracture stages, volume of proppant per stage, and the volume of injected fluids per stage. The performance of methods is compared based on a common test sample. The analysis then examines the comparative predictive performance of the three algorithms for 1330 wells that had initiated production after the initial 7311 well sample had been producing. The computations of predictor importance identified the initial 180‐day water cut and the 30‐day initial gas production predictors as having a dominant influence in most subareas and for the composite area. The relative importance of well completion predictor variables, that is, the number of fracture stages per well, volume of injected proppant per stage, volume of injected fluids per stage, and lateral length, varied considerably across the subareas. For the common test or holdout sample, the models calibrated with the XGBoost algorithm had superior predictive power. The predictive power of all the algorithms trained on the data from the original sample suffered some loss when tested with a sample of wells that had started production after the end of that period. Implications of the empirical findings and strategies to mitigate loss of predictive power are discussed in the concluding section.

Statistical Analysis and Data Mining↗

Plan for the design, development, and implementation, and operation of the National Water Information System

The Water Resources Division of the U.S. Geological Survey is developing a National Water Information System (NWIS) that will integrate and replace its existing water data and information systems of the National Water Data Storage and Retrieval System, National Water Data Exchange, National Water-Use Information, and Water Resources Scientific Information Center programs. It will be a distributed data system operated as part of the Division 's Distributed Information System, which is a network of computers linked together through a national telecommunication network known as GEONET. The NWIS is being developed as a series of prototypes that will be integrated as they are completed to allow the development and implementation of the system in a phased manner. It also is being developed in a distributed manner using personnel who work under the coordination of a central NWIS Project Office. Work on the development of the NWIS began in 1983 and it is scheduled for completion in 1990. This document presents an overall plan for the design, development, implementation, and operation of the system. Detailed discussions are presented on each of these phases of the NWIS life cycle. The planning, quality assurance, and configuration management phases of the life cycle also are discussed. The plan is intended to be a working document for use by NWIS management and participants in its design and development and to assist offices of the Division in planning and preparing for installation and operation of the system. (Author 's abstract)

Open-File Report↗

A dynamic water-quality modeling framework for the Neuse River estuary, North Carolina

As a result of fish kills in the Neuse River estuary in 1995, nutrient reduction strategies were developed for point and nonpoint sources in the basin. However, because of the interannual variability in the natural system and the resulting complex hydrologic-nutrient inter- actions, it is difficult to detect through a short-term observational program the effects of management activities on Neuse River estuary water quality and aquatic health. A properly constructed water-quality model can be used to evaluate some of the potential effects of manage- ment actions on estuarine water quality. Such a model can be used to predict estuarine response to present and proposed nutrient strategies under the same set of meteorological and hydrologic conditions, thus removing the vagaries of weather and streamflow from the analysis. A two-dimensional, laterally averaged hydrodynamic and water-quality modeling framework was developed for the Neuse River estuary by using previously collected data. Development of the modeling framework consisted of (1) computational grid development, (2) assembly of data for model boundary conditions and model testing, (3) selection of initial values of model parameters, and (4) limited model testing. The model domain extends from Streets Ferry to Oriental, N.C., includes seven lateral embayments that have continual exchange with the main- stem of the estuary, three point-source discharges, and three tributary streams. Thirty-five computational segments represent the mainstem of the estuary, and the entire framework contains a total of 60 computa- tional segments. Each computational cell is 0.5 meter thick; segment lengths range from 500 meters to 7,125 meters. Data that were used to develop the modeling framework were collected during March through October 1991 and represent the most comprehensive data set available prior to 1997. Most of the data were collected by the North Carolina Division of Water Quality, the University of North Carolina Institute of Marine Sciences, and the U.S. Geological Survey. Limitations in the modeling framework were clearly identified. These limitations formed the basis for a set of suggestions to refine the Neuse River estuary water-quality model.

North Carolina↗