USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers and Fluids”

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 37 records · Page 2Linked to original sources

Analysis of alternative weir designs for improved passage of select fish at the U.S. Geological Survey streamgaging weir at Blackwells Mills, New Jersey

As the population of New Jersey continues to remain dense, the need for water supply will likely continue to be high, which can lead to water managers needing to make difficult decisions about managing drinking-water supply. Streamgaging weirs like the ones used by the U.S. Geological Survey (USGS) play a critical role in providing accurate and stable streamflow data, but their presence can affect the passage of diadromous fish species such as river herring ( Alosa pseudoharengus [alewife], Alosa aestivalis [blueback herring], and Alosa sapidissima [American shad]). In some situations, weirs existing in rivers and streams are no longer used because they were part of a farm irrigation system or some type of industrial operation. The weir at the USGS streamgage 01402000 Millstone River at Blackwells Mills, New Jersey, was purposefully built as a hydraulic-control structure that provides a precise and stable control for the measurement of stage and computation of continuous streamflow. To satisfy the dual need of maintaining accurate streamflow data and providing improved fish passage for select species of fish during migration season, the USGS proposed the development and evaluation of two alternative weir designs that would meet the criteria established for successful passage of American shad, alewife, and blueback herring during their yearly migration. The designs were also required to maintain adequate control of the upstream pool elevation necessary for the precise computation of streamflow used by State agencies for municipal water-supply purposes for surrounding communities. Two alternative weir design modifications were incorporated at the center of the Blackwells Mills weir and modeled using two-dimensional hydraulic modeling software and three-dimensional computational fluid-dynamics software to simultaneously evaluate conditions for passage of the target fish species and effects to streamflow computations at the streamgage. The models were calibrated to existing conditions around the weir location using surveyed-elevation data and recorded stage, streamflow, and velocity in the Millstone River. The alternative weir designs lowered the weir crest by 1.02 feet (ft) and the resulting simulations showed an effective increase in depth of 0.98 ft at the median streamflow of 251 cubic feet per second (ft 3 /s) and 0.96 ft at the 95-percent exceedance streamflow of 98 ft 3 /s. The alternative weir designs were also found to increase streamflow depth across the shallowest portions of the weir structure at the downstream anti-scour skirt by lowering the skirt about 4 inches, allowing for two or more body depths of water for American shad, alewife, and blueback herring at the median migration streamflow of 251 ft 3 /s. The alternative weir designs also reduced the highest stream velocities across the downstream weir sill and anti-scour skirt from about 9 to 10 feet per second, and the depth-averaged velocity to about 7 to 8 feet per second. The sensitivity of the weir with respect to the computation of streamflow was increased from about 1.8 cubic feet per second per hundredth foot to 1.6 cubic feet per second per hundredth foot for streamflows of about 10–100 cubic feet per second.

New Jersey↗

Documentation of computer program VS2D to solve the equations of fluid flow in variably saturated porous media

This report documents FORTRAN computer code for solving problems involving variably saturated single-phase flow in porous media. The flow equation is written with total hydraulic potential as the dependent variable, which allows straightforward treatment of both saturated and unsaturated conditions. The spatial derivatives in the flow equation are approximated by central differences, and time derivatives are approximated either by a fully implicit backward or by a centered-difference scheme. Nonlinear conductance and storage terms may be linearized using either an explicit method or an implicit Newton-Raphson method. Relative hydraulic conductivity is evaluated at cell boundaries by using either full upstream weighting, the arithmetic mean, or the geometric mean of values from adjacent cells. Nonlinear boundary conditions treated by the code include infiltration, evaporation, and seepage faces. Extraction by plant roots that is caused by atmospheric demand is included as a nonlinear sink term. These nonlinear boundary and sink terms are linearized implicitly. The code has been verified for several one-dimensional linear problems for which analytical solutions exist and against two nonlinear problems that have been simulated with other numerical models. A complete listing of data-entry requirements and data entry and results for three example problems are provided. (USGS)

Water-Resources Investigations Report↗

Induced stresses due to fluid extraction from axisymmetric reservoirs

Earthquakes can be induced by fluid extraction, as well as by fluid injection. Segall (1989) proposed that poroelastic stresses are responsible for inducing earthquakes associated with fluid extraction. Here, I present methods for computing poroelastic stress changes due to fluid extraction for general axisymmetric reservoir geometries. The results of Geertsma (1973) for a thin disk reservoir with uniform pressure drop are recovered as a special case. Predicted surface subsidence agrees very well with measured leveling changes over the deep Lacq gas field in southwestern France. The induced stresses are finite if the reservoir pressure changes are continuous. Computed stress changes are on the order of several bars, suggesting that the preexisting stress states in regions of extraction induced seismicity are very close to frictional instability prior to production. ?? 1992 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

TopoDrive and ParticleFlow--Two computer models for simulation and visualization of ground-water flow and transport of fluid particles in two dimensions

This report serves as a user?s guide for two computer models: TopoDrive and ParticleFlow. These two-dimensional models are designed to simulate two ground-water processes: topography-driven flow and advective transport of fluid particles. To simulate topography-driven flow, the user may specify the shape of the water table, which bounds the top of the vertical flow section. To simulate transport of fluid particles, the model domain is a rectangle with overall flow from left to right. In both cases, the flow is under steady state, and the distribution of hydraulic conductivity may be specified by the user. The models compute hydraulic head, ground-water flow paths, and the movement of fluid particles. An interactive visual interface enables the user to easily and quickly explore model behavior, and thereby better understand ground-water flow processes. In this regard, TopoDrive and ParticleFlow are not intended to be comprehensive modeling tools, but are designed for modeling at the exploratory or conceptual level, for visual demonstration, and for educational purposes.

Open-File Report↗

Ductile creep and compaction: A mechanism for transiently increasing fluid pressure in mostly sealed fault zones

A simple cyclic process is proposed to explain why major strike-slip fault zones, including the San Andreas, are weak. Field and laboratory studies suggest that the fluid within fault zones is often mostly sealed from that in the surrounding country rock. Ductile creep driven by the difference between fluid pressure and lithostatic pressure within a fault zone leads to compaction that increases fluid pressure. The increased fluid pressure allows frictional failure in earthquakes at shear tractions far below those required when fluid pressure is hydrostatic. The frictional slip associated with earthquakes creates porosity in the fault zone. The cycle adjusts so that no net porosity is created (if the fault zone remains constant width). The fluid pressure within the fault zone reaches long-term dynamic equilibrium with the (hydrostatic) pressure in the country rock. One-dimensional models of this process lead to repeatable and predictable earthquake cycles. However, even modest complexity, such as two parallel fault splays with different pressure histories, will lead to complicated earthquake cycles. Two-dimensional calculations allowed computation of stress and fluid pressure as a function of depth but had complicated behavior with the unacceptable feature that numerical nodes failed one at a time rather than in large earthquakes. A possible way to remove this unphysical feature from the models would be to include a failure law in which the coefficient of friction increases at first with frictional slip, stabilizing the fault, and then decreases with further slip, destabilizing it. ?? 1994 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

The frontal method in hydrodynamics simulations

The frontal solution method has proven to be an effective means of solving the matrix equations resulting from the application of the finite element method to a variety of problems. In this study, several versions of the frontal method were compared in efficiency for several hydrodynamics problems. Three basic modifications were shown to be of value: 1. Elimination of equations with boundary conditions beforehand, 2. Modification of the pivoting procedures to allow dynamic management of the equation size, and 3. Storage of the eliminated equations in a vector. These modifications are sufficiently general to be applied to other classes of problems.

Computers and Fluids↗

Analysis of spurious oscillation modes for the shallow water and Navier-Stokes equations

The origin and nature of spurious oscillation modes that appear in mixed finite element methods are examined. In particular, the shallow water equations are considered and a modal analysis for the one-dimensional problem is developed. From the resulting dispersion relations we find that the spurious modes in elevation are associated with zero frequency and large wave number (wavelengths of the order of the nodal spacing) and consequently are zero-velocity modes. The spurious modal behavior is the result of the finite spatial discretization. By means of an artificial compressibility and limiting argument we are able to resolve the similar problem for the Navier-Stokes equations. The relationship of this simpler analysis to alternative consistency arguments is explained. This modal approach provides an explanation of the phenomenon in question and permits us to deduce the cause of the very complex behavior of spurious modes observed in numerical experiments with the shallow water equations and Navier-Stokes equations. Furthermore, this analysis is not limited to finite element formulations, but is also applicable to finite difference formulations.

Computers and Fluids↗

SUTRA: A model for 2D or 3D saturated-unsaturated, variable-density ground-water flow with solute or energy transport

SUTRA (Saturated-Unsaturated Transport) is a computer program that simulates fluid movement and the transport of either energy or dissolved substances in a subsurface environment. This upgraded version of SUTRA adds the capability for three-dimensional simulation to the former code (Voss, 1984), which allowed only two-dimensional simulation. The code employs a two- or three-dimensional finite-element and finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated: 1) fluid density-dependent saturated or unsaturated ground-water flow; and 2) either (a) transport of a solute in the ground water, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay; or (b) transport of thermal energy in the ground water and solid matrix of the aquifer. SUTRA may also be used to simulate simpler subsets of the above processes. A flow-direction-dependent dispersion process for anisotropic media is also provided by the code and is introduced in this report. As the primary calculated result, SUTRA provides fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA flow simulation may be employed for two-dimensional (2D) areal, cross sectional and three-dimensional (3D) modeling of saturated ground-water flow systems, and for cross sectional and 3D modeling of unsaturated zone flow. Solute-transport simulation using SUTRA may be employed to model natural or man-induced chemical-species transport including processes of solute sorption, production, and decay. For example, it may be applied to analyze ground-water contaminant transport problems and aquifer restoration designs. In addition, solute-transport simulation with SUTRA may be used for modeling of variable-density leachate movement, and for cross sectional modeling of saltwater intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between freshwater and saltwater. SUTRA energy-transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal-energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. Mesh construction, which is quite flexible for arbitrary geometries, employs quadrilateral finite elements in 2D Cartesian or radial-cylindrical coordinate systems, and hexahedral finite elements in 3D systems. 3D meshes are currently restricted to be logically rectangular; in other words, they are similar to deformable finite-difference-style grids. Permeabilities may be anisotropic and may vary in both direction and magnitude throughout the system, as may most other aquifer and fluid properties. Boundary conditions, sources and sinks may be time dependent. A number of input data checks are made to verify the input data set. An option is available for storing intermediate results and restarting a simulation at the intermediate time. Output options include fluid velocities, fluid mass and solute mass or energy budgets, and time-varying observations at points in the system. Both the mathematical basis for SUTRA and the program structure are highly general, and are modularized to allow for straightforward addition of new methods or processes to the simulation. The FORTRAN-90 coding stresses clarity and modularity rather than efficiency, providing easy access for later modifications.

Water-Resources Investigations Report↗

A numerical investigation of head waves and leaky modes in fluid- filled boreholes

Although synthetic borehole seismograms can be computed for a wide range of borehole conditions, the physical nature of shear and compressional head waves in fluid-filled boreholes is poorly understood. This paper presents a series of numerical experiments designed to explain the physical mechanisms controlling head-wave propagation in boreholes. These calculations demonstrate the existence of compressional normal modes equivalent to shear normal modes, or pseudo-Rayleigh waves, with sequential cutoff frequencies spaced between the cutoff frequencies for the shear normal modes. Major contributions to head-wave spectra occur in discrete peaks at frequencies just below mode cutoff for both compressional and shear modes. This result is confirmed by calculations with synthetic waveforms at frequencies corresponding to mode cutoff, and by branch-cut integrals designed to yield independent spectra for the compressional mode. For soft formations where shear velocity falls below acoustic velocity in the borehole fluid, leaky compressional normal modes attain properties similar to those of shear normal modes in the case of hard rock. In the limit of vanishing S-wave velocity, this result is formally related to a fluid-fluid waveguide with undamped compressional normal modes. Synthetic waveforms demonstrate that high-amplitude arrivals, traveling at velocities less than the acoustic velocity of the borehole fluid and at frequencies above a few kilohertz, represent the Airy phase of the compressional mode and are not a tube wave. Comparison of synthetic waveforms with waveforms obtained in soft sea sediments indicates that the predicted Airy phase arrivals are present in the experimental data.

Geophysics↗

An experiment in earthquake control at Rangely, Colorado

An experiment in an oil field at Rangely, Colorado, has demonstrated the feasibility of earthquake control. Variations in seismicity were produced by controlled variations in the fluid pressure in a seismically active zone. Precise earthquake locations revealed that the earthquakes clustered about a fault trending through a zone of high pore pressure produced by secondary recovery operations. Laboratory measurements of the frictional properties of the reservoir rocks and an in situ stress measurement made near the earthquake zone were used to predict the fluid pressure required to trigger earthquakes on preexisting fractures. Fluid pressure was controlled by alternately injecting and recovering water from wells that penetrated the seismic zone. Fluid pressure was monitored in observation wells, and a computer model of the reservoir was used to infer the fluid pressure distributions in the vicinity of the injection wells. The results of this experiment confirm the predicted effect of fluid pressure on earthquake activity and indicate that earthquakes can be controlled wherever we can control the fluid pressure in a fault zone.

Colorado↗

A depth-averaged debris-flow model that includes the effects of evolving dilatancy: II. Numerical predictions and experimental tests.

We evaluate a new depth-averaged mathematical model that is designed to simulate all stages of debris-flow motion, from initiation to deposition. A companion paper shows how the model’s five governing equations describe simultaneous evolution of flow thickness, solid volume fraction, basal pore-fluid pressure, and two components of flow momentum. Each equation contains a source term that represents the influence of state-dependent granular dilatancy. Here we recapitulate the equations and analyze their eigenstructure to show that they form a hyperbolic system with desirable stability properties. To solve the equations we use a shock-capturing numerical scheme with adaptive mesh refinement, implemented in an open-source software package we call D-Claw. As tests of D-Claw, we compare model output with results from two sets of large-scale debris-flow experiments. One set focuses on flow initiation from landslides triggered by rising pore-water pressures, and the other focuses on downstream flow dynamics, runout, and deposition. D-Claw performs well in predicting evolution of flow speeds, thicknesses, and basal pore-fluid pressures measured in each type of experiment. Computational results illustrate the critical role of dilatancy in linking coevolution of the solid volume fraction and pore-fluid pressure, which mediates basal Coulomb friction and thereby regulates debris-flow dynamics.

Proceedings of the Royal Society A↗

Integrated model of the shallow and deep hydrothermal systems in the East Mesa area, Imperial Valley, California

Geological, geophysical, thermal, petrophysical and hydrological data available for the East Mesa hydrothermal system that are pertinent to the construction of a computer model of the natural flow of heat and fluid mass within the system are assembled and correlated. A conceptual model of the full system is developed and a subregion selected for quantitative modeling. By invoking the .Boussinesq approximation, valid for describing the natural flow of heat and mass in a liquid hydrothermal system, it is found practical to carry computer simulations far enough in time to ensure that steady-state conditions are obtained. Initial calculations for an axisymmetric model approximating the system demonstrate that the vertical formation permeability of the deep East Mesa system must be very low (kv ~ 0.25 to 0.5 md). Since subsurface temperature and surface heat flow data exhibit major deviations from the axisymmetric approximation, exploratory three-dimensional calculations are performed to assess the effects of various mechanisms which might operate to produce such observed asymmetries. A three-dimensional model evolves from this iterative data synthesis and computer analysis which includes a hot fluid convective source distributed along a leaky fault radiating northward from the center of the hot spot and realistic variations in the reservoir formation properties.

Open-File Report↗

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↗

The hydrogeology of Kilauea volcano

The hydrogeology of Kilauea volcano and adjacent areas has been studied since the turn of this century. However, most studies to date have focused on the relatively shallow, low-salinity parts of the ground-water system, and the deeper hydrothermal system remains poorly understood. The rift zones of adjacent Mauna Loa volcano bound the regional ground-water flow system that includes Kilauea, and the area bounded by the rift zones of Kilauea and the ocean may comprise a partly isolated subsystem. Rates of ground-water recharge vary greatly over the area and discharge is difficult to measure, because streams are ephemeral and most ground-water discharges diffusely at or below sea level. Hydrothermal systems exist at depth in Kilauea's cast and southwest rift zone, as evidenced by thermal springs at the coast and wells in the lower east-rift zone. Available data suggest that dike-impounded, heated ground water occurs at relatively high elevations in the upper east-and southwest-rift zones of Kilauea, and that permeability at depth in the rift zones (probably < 10-15 m2) is much lower than that of unaltered basalt flows closer to the surface (> 10 10 m2). Substantial variations in permeability and the presence of magmatic heat sources influence the structure of the fresh water-salt water interface, so the Ghyben-Herzberg model will often fail to predict its position. Numerical modeling studies have considered only subsets of the hydrothermal system, because no existing computer code solves the coupled fluid-flow, heat- and solute-transport problem over the temperature and salinity range encountered at Kilauea. ?? 1993.

Geothermics↗

Environmental contaminants and biomarker responses in fish from the Columbia River and its tributaries: spatial and temporal trends

Fish were collected from 16 sites on rivers in the Columbia River Basin (CRB) from September 1997 to April 1998 to document temporal and spatial trends in the concentrations of accumulative contaminants and to assess contaminant effects on the fish. Sites were located on the mainstem of the Columbia River and on the Snake, Willamette, Yakima, Salmon, and Flathead Rivers. Common carp (Cyprinus carpio), black bass (Micropterus sp.), and largescale sucker (Catostomus macrocheilus) were the targeted species. Fish were field-examined for external and internal lesions, selected organs were weighed to compute somatic indices, and tissue and fluid samples were preserved for fish health and reproductive biomarker analyses. Composite samples of whole fish, grouped by species and gender, from each site were analyzed for organochlorine and elemental contaminants using instrumental methods and for 2,3,7,8-tetrachloro dibenzo-p-dioxin-like activity (TCDD-EQ) using the H4IIE rat hepatoma cell bioassay. Overall, pesticide concentrations were greatest in fish from lower CRB sites and elemental concentrations were greatest in fish from upper CRB sites. These patterns reflected land uses. Lead (Pb) concentrations in fish from the Columbia River at Northport and Grand Coulee, Washington (WA) exceeded fish and wildlife toxicity thresholds (> 0.4 ??g/g). Selenium (Se) concentrations in fish from the Salmon River at Riggins, Idaho (ID), the Columbia River at Vernita Bridge, WA, and the Yakima River at Granger, WA exceeded toxicity thresholds for piscivorous wildlife (> 0.6 ??g/g). Mercury (Hg) concentrations in fish were elevated throughout the basin but were greatest (> 0.4 ??g/g) in predatory fish from the Salmon River at Riggins, ID, the Yakima River at Granger, WA, and the Columbia River at Warrendale, Oregon (OR). Residues of p,p???-DDE were greatest (> 0.8 ??g/g) in fish from agricultural areas of the Snake, Yakima, and Columbia River basins but were not detected in upper CRB fish. Other organochlorine pesticides did not exceed toxicity thresholds in fish or were detected infrequently. Total polychlorinated biphenyls (PCBs; > 0.11 ??g/g) and TCDD-EQs (> 5 pg/g) exceeded wildlife guidelines in fish from the middle and lower CRB, and ethoxyresorufin O-deethylase (EROD) activity was also elevated at many of the same sites. Temporal trend analysis indicated decreasing or stable concentrations of Pb, Se, Hg, p,p???-DDE, and PCBs at most sites where historical data were available. Altered biomarkers were noted in fish throughout the CRB. Fish from some stations had responded to chronic contaminant exposure as indicated by fish health and reproductive biomarker results. Although most fish from some sites had grossly visible external or internal lesions, histopathological analysis determined these to be inflammatory responses associated with helminth or myxosporidian parasites. Many largescale sucker from the Columbia River at Northport and Grand Coulee, WA had external lesions and enlarged spleens, which were likely associated with infections. Intersex male smallmouth bass (Micropterus dolomieu) were found in the Snake River at Lewiston, ID and the Columbia River at Warrendale, OR. Male bass, carp, and largescale sucker containing low concentrations of vitellogenin were common in the CRB, and comparatively high concentrations (> 0.3 mg/mL) were measured in male fish from the Flathead River at Creston, Montana, the Snake River at Ice Harbor Dam, WA, and the Columbia River at Vernita Bridge, WA and Warrendale, OR. Results from our study and other investigations indicate that continued monitoring in the CRB is warranted to identify consistently degraded sites and those with emerging problems. ?? 2005 Elsevier B.V. All rights reserved.

Columbia River Basin↗

Simulation of fluid, heat transport to estimate desert stream infiltration

In semiarid regions, the contribution of infiltration from intermittent streamflow to ground water recharge may be quantified by comparing simulations of fluid and heat transport beneath stream channels to observed ground temperatures. In addition to quantifying natural recharge, streamflow infiltration estimates provide a means to characterize the physical properties of stream channel sediments and to identify suitable locations for artificial recharge sites. Rates of winter streamflow infiltration along stream channels are estimated based on the cooling effect of infiltrated water on streambed sediments, combined with the simulation of two-dimensional fluid and heat transport using the computer program VS2DH. The cooling effect of ground water is determined by measuring ground temperatures at regular intervals beneath stream channels and nearby channel banks in order to calculate temperature-depth profiles. Additional data inputs included the physical, hydraulic, and thermal properties of unsaturated alluvium, and monthly ground temperatures measurements over an annual cycle. Observed temperatures and simulation results can provide estimates of the minimum threshold for deep infiltration, the variability of infiltration along stream channels, and also the frequency of infiltration events.

Ground Water↗

Biomonitoring of Environmental Status and Trends (BEST) Program: Environmental contaminants and their effects on fish in the Columbia River Basin

This project examined and analyzed 560 fish representing eight species from 16 sites in the Columbia River Basin (CRB) from September 1997 to April 1998. Ten of the 16 sampling locations were historical National Contaminant Biomonitoring Program (NCBP) sites where organochlorine and elemental contaminants in fish had been monitored from 1969 through 1986. Five sites were co-located at U.S. Geological Survey (USGS)-National Stream Quality Accounting Network (NASQAN) stations at which water quality is monitored. The sampling location at Marine Park in Vancouver, Washington did not correspond to either of the established monitoring programs. Eight of the sampling locations were located on the Columbia River; three were on the Snake River; two were on the Willamette River, and one site was on each of the Yakima, Salmon and Flathead Rivers. Common carp ( Cyprinus carpio ), black basses ( Micropterus sp. ), and largescale sucker ( Catostomus macrocheilus ) together accounted for 80% of the fish sampled during the study. Fish were weighed and measured then field-examined for external and internal lesions, and liver, spleen, and gonads were weighed to compute somatic indices. Selected tissues and fluids were obtained and preserved for analysis of fish health and reproductive biomarkers. Composite samples of whole fish from each station were grouped by species and gender and analyzed for persistent organic and inorganic contaminants and for dioxin-like activity using H4IIE rat hepatoma cell bioassay.

Columbia River Basin↗