Search USGSSearch

SEARCH · Search USGS

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

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

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

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

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

This project collected, examined, and analyzed 217 fish representing three species at 10 stations in the U.S. portion of the Yukon River Basin (YRB) from May to October 2002. Four sampling sites were located on the Yukon River; two were located on the Porcupine River, and one site was on each of the Ray, Tanana, Tolavana, and Innoko Rivers. Norther pike ( Esox lucius ), longnose sucker ( Catostomus catostomus ), and burbot ( Lota lota ) were weighed and measured, and examined in the field for external and internal lesions, and liver, spleen, and gonads were weighed to compute somatic indices. Selected tissues and fluids were collected 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 organochlorines and elemental contaminants and for dioxin-like activity using H4IIE rat hepatoma cell bioassay.

Alaska

Biomonitoring of Environmental Status and Trends (BEST) Program: Environmental contaminants, health indicators, and reproductive biomarkers in fish from the Colorado River basin

Seven fish species were collected from 14 sites on rivers in the Colorado River Basin (CDRB) from August to October 2003. Spatial trends in the concentrations of accumulative contaminants were documented and contaminant effects on the fish were assessed. Sites were located on the mainstem of the Colorado River and on the Yampa, Green, Gunnison, San Juan, and Gila Rivers. Common carp (Cyprinus carpio), black bass (Micropterus sp.), and channel catfish (Ictalurus punctatus) were the targeted species. Fish were field-examined for external and internal anomalies, 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 performance-based and instrumental methods. 2,3,7,8-tetrachlorodibenzo-p-dioxin-like activity (TCDD-EQ) was measured using the H4IIE rat hepatoma cell bioassay. Selenium (Se) and mercury (Hg) concentrations were elevated throughout the CDRB, and pesticides concentrations were greatest in fish from agricultural areas in the Lower Colorado River and Gila River. Selenium concentrations exceeded toxicity thresholds for fish (>1.0 ?g/g ww) at all sites except from the Gila River at Hayden, Arizona. Mercury concentrations were elevated (>0.1 ?g/g ww) in fish from the Yampa River at Lay, Colorado; the Green River at Ouray National Wildlife Refuge (NWR), Utah and San Rafael, Utah; the San Juan River at Hogback Diversion, New Mexico; and the Colorado River at Gold Bar Canyon, Utah, Needles, California, and Imperial Dam, Arizona. Concentrations of p,p'-DDE were relatively high in fish from Arlington, Arizona (>1.0 ?g/g ww) and Phoenix, Arizona (>0.5 ?g/g ww). Concentrations of other banned pesticides including toxaphene, total chlordanes, and dieldrin were also greatest at these two sites but did not exceed toxicity thresholds. Current-use or unlisted pesticides such as dacthal, endosulfan, '-HCH, and methoxychlor were also greatest in fish from Gila River. Total polychlorinated biphenyls (PCBs; >0.11 ?g/g ww) and TCDD-EQs (>5 pg/g ww) exceeded wildlife guidelines in fish from the Gila River at Phoenix, Arizona. Hepatic ethoxyresorufin O-deethylase (EROD) activity was also relatively high in carp from the Gila River at Phoenix, Arizona and in bass from the Green River at Ouray NWR, Utah. Altered biomarkers were noted in fish throughout the CDRB. Fish from some stations responded to chronic contaminant exposure as indicated by fish health and reproductive biomarker results. Multiple fish health indicators including altered body and organ weights and high health assessment index scores may be associated with elevated Se concentrations in fish from the Colorado River at Loma, Colorado and Needles, California. Although grossly visible external or internal lesions were found on most fish from some sites, histopathological analysis determined many of these to be inflammatory responses associated with parasites. Edema, exophthalmos, and cataracts were noted in fish from sites with elevated Se concentrations. Reproductive biomarkers including gonad development and maturation, vitellogenin concentrations, and steroid hormone concentrations were anomalous in fish from the Gila River at Hayden and Phoenix, Arizona. In addition, intersex fish were found at seven of 14 sites. The intersex condition was identified in smallmouth bass (M. dolomieu), largemouth bass (M. salmoides), channel catfish, and carp and may indicate exposure to endocrine disrupting compounds. Seven of ten male smallmouth bass from the Yampa River at Lay, Colorado were intersex. Male carp, bass, and channel catfish with low concentrations of vitellogenin were common in the CDRB. Comparatively high vitellogenin concentrations (>0.2 mg/mL) were measured in male fish from the Green River at Ouray NWR, Utah and the Colorado River at Im

Scientific Investigations Report

Biomonitoring of Environmental Status and Trends (BEST) Program: Environmental contaminants, health indicators, and reproductive biomarkers in fish from the Mobile, Apalachicola-Chattahoochee-Flint, Savannah, and Pee Dee River Basins

Largemouth bass ( Micropterus salmoides ) and common carp ( Cyprinus carpio ) were collected from 13 sites in 4 river basins in the southeastern United States to document spatial trends in accumulative contaminants, health indicators, and reproductive biomarkers. Organochlorine residues, 2,3,7,8- tetrachlorodibenzo-p-dioxin-like activity (TCDD-EQ), and elemental contaminants were measured in composite samples of whole fish, grouped by species and gender, from each site. Fish were field-examined for external and internal anomalies, selected organs were weighed to compute somatic indices, and tissue and fluid samples were preserved for fish health and reproductive biomarker analyses. Mercury concentrations in bass samples from all sites exceeded toxicity thresholds for mammals [>0.1 micrograms per gram wet weight (ug/g ww)], fish (>0.2 ug/g ww), and birds (>0.3 ug/g ww) and were greatest (>0.5 ug/g ww) in samples from the Alabama River at Eureka Landing, Alabama; the Mobile River at Bucks, Alabama; the Apalachicola River at Blountstown, Florida; the Savannah River at Sylvania, Georgia; and the Pee Dee River at Bucksport, South Carolina. Selenium concentrations were relatively high (>0.75 ug/g ww) in fish from the Tombigbee River at Lavaca, Alabama; the Mobile River at Bucks; and the Chattahoochee River at Omaha, Georgia compared to those from other sites. Concentrations of 2,2-bis (p-chlorophenyl)- 1,1-dichloroethylene (p,p'-DDE) were high in fish from the Chattahoochee River at Omaha and the Mobile River near Bucks, which was near a 2,2-bis (p-chlorophenyl)-1,1- dichloroethylene (DDT) formulating facility that historically discharged into the lower Mobile River. Toxaphene concentrations in fish from the Flint River near Albany, Georgia (60-100 nanograms per gram (ng/g) ww) may pose a risk to fish. Concentrations of other formerly used (total chlordanes, dieldrin, endrin, aldrin, mirex, and hexachlorobenzene) and currently used (pentachlorobenzene, pentachloroanisole, dacthal, endosulfan, y-HCH, and methoxychlor) organochlorine residues generally were low or did not exceed toxicity thresholds. Total polychlorinated biphenyls concentrations in samples from the Coosa River at Childersburg, Alabama; the Apalachicola River at Omaha; the Apalachicola River at Blountstown; and the Pee Dee River at Bucksport were >480 ng/g ww and may be a risk to piscivorous wildlife. Dioxin-like activity as measured by TCDD-EQ was greatest [>10 picograms per gram (pg/g)] in male fish from the Coosa River at Childersburg and the Mobile River at Bucks. Hepatic ethoxyresorufin O-deethylase activity generally was greatest in carp from the Mobile River Basin [means >10 picomols per minute per milligram of protein (pmol/min/mg)] and in bass from the Tombigbee River at Lavaca and Pee Dee River at Pee Dee, South Carolina (means >65 pmol/min/mg). Altered biomarkers were noted in fish from all basins. The field necropsy and histopathological examination determined that fish from the Mobile River Basin generally were in poorer health than those from the other basins. In bass, health assessment index scores were correlated with mercury and p,p'-DDE concentrations. High health assessment index scores in Mobile River Basin fish were widespread and caused primarily by parasitic infestations, which were most severe in fish from the Tombigbee River at Lavaca and the Alabama River at Eureka Landing. Tumors were present in few fish (n = 5; 0.01%). Ovarian tumors of the same origin (smooth muscle) were present in two older carp from the Chattahoochee River near Omaha, Georgia and may be contaminant related. Reproductive biomarkers including gonadosomatic index, vitellogenin concentrations, and steroid hormone concentrations were anomalous in fish from various sites but were not consistently related to any particular chemical contaminant. Intersex gonads were identified in 47 male bass (42%) representing 12 sites and may indicate exposure to endocrine disrupting comp

Mobile, Apalachicola-Chattahoochee-Flint, Savannah

Addition of simultaneous heat and solute transport and variable fluid viscosity to SEAWAT

SEAWAT is a finite-difference computer code designed to simulate coupled variable-density ground water flow and solute transport. This paper describes a new version of SEAWAT that adds the ability to simultaneously model energy and solute transport. This is necessary for simulating the transport of heat and salinity in coastal aquifers for example. This work extends the equation of state for fluid density to vary as a function of temperature and/or solute concentration. The program has also been modified to represent the effects of variable fluid viscosity as a function of temperature and/or concentration. The viscosity mechanism is verified against an analytical solution, and a test of temperature-dependent viscosity is provided. Finally, the classic Henry-Hilleke problem is solved with the new code. ?? 2006 Elsevier Ltd. All rights reserved.

Computers & Geosciences