USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cells”

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 595 records · Page 33Linked to original sources

Potential interaction of groundwater and surface water including autonomous underwater vehicle reconnaissance at Nolin River Lake, Kentucky, 2016

The U.S. Geological Survey collaborated with the U.S. Army Corps of Engineers, Louisville District, on a synoptic study of water quality at Nolin River Lake during August 2016. The purpose of the study was to develop a better understanding of the potential for interaction between groundwater and surface water at Nolin River Lake, Kentucky. Groundwater can have properties that are measurably different from those in adjacent surface water, and inflows and outflows can be an important component of water quality and quantity. An improved understanding of potential interaction of groundwater and surface water at Nolin River Lake may be used to refine lake-management strategies. This study (1) compiled and interpreted existing information to characterize the hydrogeological setting and implications for potential interaction of groundwater and surface water in the Nolin River Lake watershed; (2) collected transects of onsite water-quality parameters using an autonomous underwater vehicle (AUV) in areas with potential for interaction of groundwater and surface water, including five sites on Nolin River Lake and one site on the Nolin River; and (3) collected discrete water-quality and phytoplankton community data at the same six sites. A review of existing hydrogeologic information did not indicate the presence of karst features adjacent to or beneath Nolin River Lake that would facilitate groundwater interaction with the reservoir. Observations leading to this conclusion include (1) limestone that is adjacent to the shoreline and perhaps beneath the lake, is overlain with siliciclastic rocks and fine-grained sediment that inhibits infiltration and development of karst features that encourage rapid groundwater flow; (2) the geologic deposits surrounding the reservoir are described as having limited or no potential for development of karst features, some exceptions may exist in tributary valleys; (3) very few karst features were mapped within 1 mile of the reservoir or in the area currently occupied by the reservoir; and (4) faults that intersect the reservoir but may not possess hydraulic properties that cause the faults to be conduits for groundwater flow. Groundwater interaction with reservoir tributaries is likely more common in areas of the watershed upstream from Nolin River Lake where karst hydrogeology is prevalent. Results of water-quality surveys using an AUV from August 15 to 19, 2016, did not identify areas of anomalous values that might indicate groundwater inflows through preferential flow zones. Spatial distributions of water-quality parameters were generally uniform within each constant-depth layer. The constant-depth layers were selected to be above, within, and below the thermocline and ranged from the water surface to 25 feet. Surveys near the bottom of the reservoir that might have been more sensitive to groundwater inflows were not done because presurvey data were not available to indicate locations of obstacles that could ensnare the AUV. Water-quality data collected with the AUV did identify water-quality anomalies where stream tributaries were discharging to the reservoir. The discrete water-quality samples indicated uniformity among the five reservoir sites. The riverine site that is immediately upstream from Nolin River Lake, however, had some unique water-quality characteristics relative to sites on the reservoir. The highest concentrations of nitrate plus nitrite as nitrogen (0.145 milligrams per liter [mg/L]), total phosphorous (0.07 mg/L), chlorophyll a (36.1 micrograms per liter), and pheophytin a (10.2 micrograms per liter) were measured at the Nolin River Lake riverine site (site 2NRR20034). The concentrations of nutrients and chlorophyll a at the riverine site did exceed the 25th percentile of median concentrations measured by the U.S. Environmental Protection Agency (EPA) at other lakes and reservoirs in EPA level IV ecoregion 71a. Concentrations of most nutrients and chlorophyll a at the five reservoir sites also exceeded the 25th percentile of median concentrations in EPA level IV ecoregion 72h. The exception was the concentrations of total phosphorus as phosphorus at the reservoir sites that were at or below the 25th percentile of median concentrations measured by EPA (0.03 mg/L). Concentrations of orthophosphate as phosphorus were less than the method detection limit of 0.004 mg/L at all sites. The phytoplankton community in Nolin River Lake was almost exclusively (greater than 90 percent of total phytoplankton abundance) cyanobacteria, also known as blue-green algae. A species of Cylindrospermopsis dominated the cyanobacterial community at the five reservoir sites, while Chroococcus microscopicus was most abundant at the riverine site. Cyanobacterial cell densities ranged from 10,000 to 198,067,460 cells per liter in five areas in the reservoir and from 4,800 to 73,751,253 cells per liter at the riverine site. Multiple potential sources of water to Nolin River Lake include direct precipitation, overland flow, interflow, groundwater, and surface water. Understanding the exact contribution of each of these components to the water budget at Nolin River Lake may help the U.S. Army Corps of Engineers manage the water quality, water quantity, and biological communities in the reservoir. Additional hydrogeologic and water-quality data that builds on the results of this study may refine the inferences of this study; for example, deeper AUV surveys that target the largest fault zones might further the understanding of the potential for groundwater flow through those features. A complete understanding of the reservoir hydrology, however, may require the use of scientific methods intended for water bodies as large as Nolin River Lake, such as aerial infrared photography and imagery; water mass, chemical, and isotopic balance studies; geophysical measurements; and numerical simulations.

Kentucky↗

Scoops3D: software to analyze 3D slope stability throughout a digital landscape

The computer program, Scoops3D, evaluates slope stability throughout a digital landscape represented by a digital elevation model (DEM). The program uses a three-dimensional (3D) method of columns approach to assess the stability of many (typically millions) potential landslides within a user-defined size range. For each potential landslide (or failure), Scoops3D assesses the stability of a rotational, spherical slip surface encompassing many DEM cells using a 3D version of either Bishop’s simplified method or the Ordinary (Fellenius) method of limit-equilibrium analysis. Scoops3D has several options for the user to systematically and efficiently search throughout an entire DEM, thereby incorporating the effects of complex surface topography. In a thorough search, each DEM cell is included in multiple potential failures, and Scoops3D records the lowest stability (factor of safety) for each DEM cell, as well as the size (volume or area) associated with each of these potential landslides. It also determines the least-stable potential failure for the entire DEM. The user has a variety of options for building a 3D domain, including layers or full 3D distributions of strength and pore-water pressures, simplistic earthquake loading, and unsaturated suction conditions. Results from Scoops3D can be readily incorporated into a geographic information system (GIS) or other visualization software. This manual includes information on the theoretical basis for the slope-stability analysis, requirements for constructing and searching a 3D domain, a detailed operational guide (including step-by-step instructions for using the graphical user interface [GUI] software, Scoops3D-i) and input/output file specifications, practical considerations for conducting an analysis, results of verification tests, and multiple examples illustrating the capabilities of Scoops3D. Easy-to-use software installation packages are available for the Windows or Macintosh operating systems; these packages install the compiled Scoops3D program, the GUI (Scoops3D-i), and associated documentation. Several Scoops3D examples, including all input and output files, are available as well. The source code is written in the Fortran 90 language and can be compiled to run on any computer operating system with an appropriate compiler.

Techniques and Methods↗

The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1

The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'. The first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain. All of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR). In addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas. The example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de

Techniques and Methods↗

Documentation of the seawater intrusion (SWI2) package for MODFLOW

The SWI2 Package is the latest release of the Seawater Intrusion (SWI) Package for MODFLOW. The SWI2 Package allows three-dimensional vertically integrated variable-density groundwater flow and seawater intrusion in coastal multiaquifer systems to be simulated using MODFLOW-2005. Vertically integrated variable-density groundwater flow is based on the Dupuit approximation in which an aquifer is vertically discretized into zones of differing densities, separated from each other by defined surfaces representing interfaces or density isosurfaces. The numerical approach used in the SWI2 Package does not account for diffusion and dispersion and should not be used where these processes are important. The resulting differential equations are equivalent in form to the groundwater flow equation for uniform-density flow. The approach implemented in the SWI2 Package allows density effects to be incorporated into MODFLOW-2005 through the addition of pseudo-source terms to the groundwater flow equation without the need to solve a separate advective-dispersive transport equation. Vertical and horizontal movement of defined density surfaces is calculated separately using a combination of fluxes calculated through solution of the groundwater flow equation and a simple tip and toe tracking algorithm. Use of the SWI2 Package in MODFLOW-2005 only requires the addition of a single additional input file and modification of boundary heads to freshwater heads referenced to the top of the aquifer. Fluid density within model layers can be represented using zones of constant density (stratified flow) or continuously varying density (piecewise linear in the vertical direction) in the SWI2 Package. The main advantage of using the SWI2 Package instead of variable-density groundwater flow and dispersive solute transport codes, such as SEAWAT and SUTRA, is that fewer model cells are required for simulations using the SWI2 Package because every aquifer can be represented by a single layer of cells. This reduction in number of required model cells and the elimination of the need to solve the advective-dispersive transport equation results in substantial model run-time savings, which can be large for regional aquifers. The accuracy and use of the SWI2 Package is demonstrated through comparison with existing exact solutions and numerical solutions with SEAWAT. Results for an unconfined aquifer are also presented to demonstrate application of the SWI2 Package to a large-scale regional problem.

Techniques and Methods↗

Documentation of a computer program to simulate lake-aquifer interaction using the MODFLOW ground water flow model and the MOC3D solute-transport model

Heads and flow patterns in surficial aquifers can be strongly influenced by the presence of stationary surface-water bodies (lakes) that are in direct contact, vertically and laterally, with the aquifer. Conversely, lake stages can be significantly affected by the volume of water that seeps through the lakebed that separates the lake from the aquifer. For these reasons, a set of computer subroutines called the Lake Package (LAK3) was developed to represent lake/aquifer interaction in numerical simulations using the U.S. Geological Survey three-dimensional, finite-difference, modular ground-water flow model MODFLOW and the U.S. Geological Survey three-dimensional method-of-characteristics solute-transport model MOC3D. In the Lake Package described in this report, a lake is represented as a volume of space within the model grid which consists of inactive cells extending downward from the upper surface of the grid. Active model grid cells bordering this space, representing the adjacent aquifer, exchange water with the lake at a rate determined by the relative heads and by conductances that are based on grid cell dimensions, hydraulic conductivities of the aquifer material, and user-specified leakance distributions that represent the resistance to flow through the material of the lakebed. Parts of the lake may become "dry" as upper layers of the model are dewatered, with a concomitant reduction in lake surface area, and may subsequently rewet when aquifer heads rise. An empirical approximation has been encoded to simulate the rewetting of a lake that becomes completely dry. The variations of lake stages are determined by independent water budgets computed for each lake in the model grid. This lake budget process makes the package a simulator of the response of lake stage to hydraulic stresses applied to the aquifer. Implementation of a lake water budget requires input of parameters including those representing the rate of lake atmospheric recharge and evaporation, overland runoff, and the rate of any direct withdrawal from, or augmentation of, the lake volume. The lake/aquifer interaction may be simulated in both transient and steady-state flow conditions, and the user may specify that lake stages be computed explicitly, semi-implicitly, or fully-implicitly in transient simulations. The lakes, and all sources of water entering the lakes, may have solute concentrations associated with them for use in solute-transport simulations using MOC3D. The Stream Package of MODFLOW-2000 and MOC3D represents stream connections to lakes, either as inflows or outflows. Because lakes with irregular bathymetry can exist as separate pools of water at lower stages, that coalesce to become a single body of water at higher stages, logic was added to the Lake Package to allow the representation of this process as a user option. If this option is selected, a system of linked pools (sublakes) is identified in each time step and stages are equalized based on current relative sublake surface areas.

Water-Resources Investigations Report↗

Simulation of net infiltration and potential recharge using a distributed-parameter watershed model of the Death Valley region, Nevada and California

This report presents the development and application of the distributed-parameter watershed model, INFILv3, for estimating the temporal and spatial distribution of net infiltration and potential recharge in the Death Valley region, Nevada and California. The estimates of net infiltration quantify the downward drainage of water across the lower boundary of the root zone and are used to indicate potential recharge under variable climate conditions and drainage basin characteristics. Spatial variability in recharge in the Death Valley region likely is high owing to large differences in precipitation, potential evapotranspiration, bedrock permeability, soil thickness, vegetation characteristics, and contributions to recharge along active stream channels. The quantity and spatial distribution of recharge representing the effects of variable climatic conditions and drainage basin characteristics on recharge are needed to reduce uncertainty in modeling ground-water flow. The U.S. Geological Survey, in cooperation with the Department of Energy, developed a regional saturated-zone ground-water flow model of the Death Valley regional ground-water flow system to help evaluate the current hydrogeologic system and the potential effects of natural or human-induced changes. Although previous estimates of recharge have been made for most areas of the Death Valley region, including the area defined by the boundary of the Death Valley regional ground-water flow system, the uncertainty of these estimates is high, and the spatial and temporal variability of the recharge in these basins has not been quantified. To estimate the magnitude and distribution of potential recharge in response to variable climate and spatially varying drainage basin characteristics, the INFILv3 model uses a daily water-balance model of the root zone with a primarily deterministic representation of the processes controlling net infiltration and potential recharge. The daily water balance includes precipitation (as either rain or snow), snow accumulation, sublimation, snowmelt, infiltration into the root zone, evapotranspiration, drainage, water content change throughout the root-zone profile (represented as a 6-layered system), runoff (defined as excess rainfall and snowmelt) and surface water run-on (defined as runoff that is routed downstream), and net infiltration (simulated as drainage from the bottom root-zone layer). Potential evapotranspiration is simulated using an hourly solar radiation model to simulate daily net radiation, and daily evapotranspiration is simulated as an empirical function of root zone water content and potential evapotranspiration. The model uses daily climate records of precipitation and air temperature from a regionally distributed network of 132 climate stations and a spatially distributed representation of drainage basin characteristics defined by topography, geology, soils, and vegetation to simulate daily net infiltration at all locations, including stream channels with intermittent streamflow in response to runoff from rain and snowmelt. The temporal distribution of daily, monthly, and annual net infiltration can be used to evaluate the potential effect of future climatic conditions on potential recharge. The INFILv3 model inputs representing drainage basin characteristics were developed using a geographic information system (GIS) to define a set of spatially distributed input parameters uniquely assigned to each grid cell of the INFILv3 model grid. The model grid, which was defined by a digital elevation model (DEM) of the Death Valley region, consists of 1,252,418 model grid cells with a uniform grid cell dimension of 278.5 meters in the north-south and east-west directions. The elevation values from the DEM were used with monthly regression models developed from the daily climate data to estimate the spatial distribution of daily precipitation and air temperature. The elevation values were also used to simulate atmosp

Water-Resources Investigations Report↗

Conjunctive-use optimization model of the Mississippi River Valley alluvial aquifer of northeastern Arkansas

The Mississippi River Valley alluvial aquifer is a water-bearing assemblage of gravels and sands that underlies about 32,000 square miles of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas. Because of the heavy demands placed on the aquifer, several large cones of depression over 100 feet deep have formed in the potentiometric surface, resulting in lower well yields and degraded water quality in some areas. A ground-water flow model of the alluvial aquifer was previously developed for an area covering 14,104 square miles, extending northeast from the Arkansas River into the northeast corner of Arkansas and parts of southeastern Missouri. The flow model showed that continued ground-water withdrawals at rates commensurate with those of 1997 could not be sustained indefinitely without causing water levels to decline below half the original saturated thickness of the aquifer. To develop estimates of withdrawal rates that could be sustained in compliance with the constraints of critical ground-water area designation, conjunctive-use optimization modeling was applied to the flow model of the alluvial aquifer in northeastern Arkansas. Ground-water withdrawal rates form the basis for estimates of sustainable yield from the alluvial aquifer and from rivers specified within the alluvial aquifer model. A management problem was formulated as one of maximizing the sustainable yield from all ground-water and surface-water withdrawal cells within limits imposed by plausible withdrawal rates, and within specified constraints involving hydraulic head and streamflow. Steady-state flow conditions were selected because the maximized withdrawals are intended to represent sustainable yield of the system (a rate that can be maintained indefinitely). Within the optimization model, 11 rivers are specified. Surface-water diversion rates that occurred in 2000 were subtracted from specified overland flow at the appropriate river cells. Included in these diversions were the planned diversions of 63,339,248 ft3/d for the Bayou Meto project area and 55,078,367 ft3/d for the Grand Prairie project area, which factor in an additional 30 and 40 percent transmission loss, respectively. Streamflow constraints were specified at all 1,165 river cells based on average 7-day minimum flows for 10 years. Sustainable yield for all rivers ranged from 0 (Current, Little Red, and Bayou Meto Rivers) to almost 5 billion cubic feet per day for the Arkansas River. Total sustainable yield from all rivers combined was 12.8 billion cubic feet per day, which represents a substantial source for supplementing ground water to meet the total water demand. Sustainable-yield estimates are affected by the allowable upper limit on withdrawals from wells specified in the optimization model. Ground-water withdrawal rates were allowed to vary as much as 200 percent of the withdrawal rate in 1997. As the overall upper limit on withdrawals is increased, the sustainable yield generally increases. Tests with the optimization model show that without limits on pumping, wells adjacent to sources of water would have optimized withdrawal rates that were orders of magnitude larger than rates corresponding to those of 1997. The sustainable yield from ground water for the entire study area while setting the maximum upper limit as the amount withdrawn in 1997 is 360 million cubic feet per day, which is only about 57 percent of the amount withdrawn in 1997 (635.6 million cubic feet per day). Optimal sustainable yields from within the Bayou Meto irrigation project area and within the Grand Prairie irrigation project area are 18.1 and 9.1 million cubic feet per day, respectively, assuming a maximum allowable withdrawal rate equal to 1997 rates. These values of sustainable yield represent 35 and 30 percent respectively of the amount pumped from these project areas in 1997. Unmet demand (defined as the difference between the optimized withdrawal rate or sustainable yield, and the anticipated demand) was calculated using different demand rates based on multiples of the 1997 withdrawal rate. Assuming that demand is the 1997 withdrawal rate, and that sustainable-yield estimates are those obtained using upper limits of withdrawal rates of 100-, 150-, and 200-percent of 1997 withdrawal rates, then the resulting unmet demand for the entire model area is 275.5, 190.9, and 110 million cubic feet per day, respectively. Whereas, if the demand is specified as 100-, 150-, and 200-percent of the 1997 withdrawal rate, and the sustainable-yield estimates remain the same, then the resulting unmet demand for the entire model area is 275.5, 508.8, and 745.8 million cubic feet per day, respectively. These unmet demands for ground water could be obtained from large sustainable surface-water withdrawals.

Arkansas↗

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida↗

Hydrologic environments and water-quality characteristics at four landfills in Mecklenburg County, North Carolina, 1980-86

A water-quality study was conducted during 1980-86 at four landfills in Mecklenburg County, North Carolina. Each landfill has a three-layered hydrogeologic system typical of the Piedmont, consisting of (1) the regolith; (2) a transition zone; and (3) unweathered, fractured crystalline bedrock. As much as 7.6 inches per year of rainfall enters the ground-water system and has the potential to generate leachate within landfill cells. Ground water and leachate discharge to tributaries within the landfill sites or to streams adjacent to them. Water-quality samples were collected from 53 monitoring wells and 20 surface-water sites. Samples were analyzed for selected physical and biological characteristics, major inorganic ions, nutrients, trace elements, and organic compounds. Selected indicators of water quality, including specific conductance; hardness; and concentrations of chloride, manganese, dissolved solids, total organic carbon, and specific organic compounds were analyzed to determine the effects of each landfill on ground- and surface-water quality. Increases in concentrations of inorganic constituents above background levels were detected in ground water downgradient of the landfills. The increases were generally greatest in samples from wells in close proximity to the older landfill cells. In general, the increases in concentrations in downgradient wells were greater for calcium, magnesium, and chloride than for other major ions. Manganese exhibited the largest relative increase in concentration between upgradient and downgradient wells of any constituent, and manganese concentration data were effective in defining areas with extensive anaerobic biological activity. Differences between upgradient and downgradient concentrations of total organic carbon and specific organic compounds generally were not as apparent. The most frequently identified organic contaminants were the herbicides 2,4-D and 2,4,5-T. Chlorofluoromethanes were identified in three of four ground-water samples analyzed for volatile organic compounds. Landfills affected the water quality of several smaller streams but did not noticeably affect larger ones. Apparent effects on water quality were greatest at the oldest landfill, located on Statesville Road, where waste is in cells that are partly below the water table.

North Carolina↗

Radon-222 concentrations in ground water and soil gas on Indian reservations in Wisconsin

The weighted average radon-222 concentration of indoor air in homes located on Wisconsin Indian Reservations is 5.8 picocuries per liter, which exceeds the U.S. Environmental Protection Agency action limit of 4 picocuries per liter. Ground water is the principle source of drinking water on Wisconsin Indian Reservations and generally accounts for about 5 percent of the total indoor air radon-222 concentrations found in homes. To determine the distribution of radon-222, ground water from 29 private and community Wisconsin Indian Reservation wells and soil gas at a depth of about 3 feet below land surface adjacent to the wells were sampled. Sites with wells were distributed among the 11 Wisconsin Indian Reservations so that each Reservation contained at least 2 sites. The remaining seven sites were divided among the Reservation by acreage held by each tribe. Ground-water samples were collected using a syringe technique after the wells had been purged sufficiently to reach chemical stability. Samples were then sent by overnight mail to the U.S. Geological Survey, National Water Quality Lab for analysis by liquid scintillation counting. Soil gas was collected in Lucas cells, using a small-diameter soil probe and peristaltic pump. The cells were then analyzed by Lucas cell alpha scintillation counting with a portable radon detector. The highest radon-222 concentrations in ground water and soil gas are from sites with wells finished in the crystalline bedrock aquifer. Radon-222 concentrations in 29 ground-water samples collected from the sand and gravel and sedimentary and crystalline bedrock aquifers range from 260 to 22,000 picocuries per liter with a median concentration of 560 picocuries per liter. Only 2 of the 29 ground-water samples have radon-222 concentrations less than the U.S. Environmental Protection Agency proposed standard of 300 picocuries per liter. The highest radon-222 concentrations were found in ground water from wells in Shawano, Menominee, Forest, and Marathon counties. Radon-222 concentrations in soil gas range from 130 to 7,810 picocuries per liter with a median concentration of 560 picocuries per liter. For sites with wells finished in the sand and gravel aquifer, the coefficient of determination (R2) of the regression of concentration of radon-222 in ground water as a function of well depth is 0.003 and the significance level is 0.32, which indicates that there is not a statistically significant relation between radon-222 concentrations in ground water and well depth. The coefficient of determination of the regression of radon-222 in ground water and soil gas is 0.19 and the root mean square error of the regression line is 271 picocuries per liter. Even though the significance level (0.036) indicates a statistical relation, the root mean square error of the regression is so large that the regression equation would not give reliable predictions. Because of an inadequate number of samples, similar statistical analyses could not be performed for sites with wells finished in the crystalline and sedimentary bedrock aquifers.

Wisconsin↗

Description and field analysis of a coupled ground-water/surface-water flow model (MODFLOW/BRANCH) with modifications for structures and wetlands in southern Dade County, Florida

A coupled surface-water model (BRANCH) and ground-water model (MODFLOW) model were tested to simulate the interacting wetlands/surface-water/ ground-water system of southern Dade County. Several options created for the MODFLOW ground- ground-water model were used in representing this field situation. The primary option is the MODBRANCH interfacing software, which allows leakage to be accounted for between the MODFLOW ground-water model and the BRANCH dynamic model for simulation of flow in an interconnected network of open channels. A modification to an existing software routine, which is referred to as BCF2, allows cells in MODFLOW to rewet when dry--a requirement in representing the seasonal wetlands in Dade County. A companion to BCF2 is the modified evapotranspiration routine EVT2. The EVT2 routine changes the cells where evapotranspiration occurs, depending on which cells are wet. The Streamlink package represents direct connections between the canals and wetlands at locations where canals open directly into overland flow. Within the BRANCH model, the capability to represent the numerous hydraulic structures, gated spillways, gated culverts, and pumps was added. The application of these modifications to model surface-water/ground-water interactions in southern Dade County demonstrated the usefulness of the coupled MODFLOW/BRANCH model. Ground-water and surface-water flows are both simulated with dynamic models. Flow exchange between models, intermittent wetting and drying, evapotranspiration, and hydraulic structure operations are all represented appropriately. Comparison was made with a simulation using the RIV1 package instead of MODBRANCH to represent the canals. RIV1 represents the canals by user-defined stages, and computes leakage to the aquifer. Greater accuracy in reproducing measured ground- water heads was achieved with MODBRANCH, which also computes dynamic flow conditions in the canals, unlike RIV1. The surface-water integrated flow and transport two-dimensional model (SWIFT2D) was also applied to the southeastern coastal wetlands for comparison with the wetlands flow approximation made in MODFLOW. MODFLOW simulates the wetlands as a highly conductive upper layer of the aquifer, whereas SWIFT2D solves the hydrodynamic equations. Comparison in this limited test demonstrated no specific advantage for either method of representation. However, much additional testing on a wider variety of geometric and hydraulic situations, such as in areas with greater tidal or other dynamic forcing effects, is needed to make definite conclusions. A submodel of the existing southern Dade County model schematization was used to examine water-delivery alternatives proposed by the U.S. Army Corps of Engineers. For this application, the coupled MODFLOW/BRANCH model was used as a design tool. A new canal and several pumps to be tested to maintain lower water levels in a residential area (while water levels in the Everglades are raised) were added to the model schematization. The pumps were assumed to have infinite supply capacity in the model so that their maximum pumping rates during the simulation could be used to determine pump sizes.

Water-Resources Investigations Report↗

Numerical simulation of ground-water flow in the central part of the western San Joaquin Valley, California

The occurrence of selenium in agricultural drain water in the central part of the western San Joaquin Valley, California, has focused concern on strategies for managing shallow, saline ground water. To assess alternatives to agricultural drains, a three-dimensional, finite-difference numerical model of the regional groundwater flow system was developed. This report documents the mathematical approach used to model the flow system, the data base on which the model is based, and the methods used to calibrate the model. The 550-square-mile study area includes parts of the Panoche Creek alluvial fan and parts of the Little Panoche Creek and Cantua Creek alluvial fans. The model simulates transient flow in the semiconfined and confined zones above and below the Corcoran Clay Member of the Tulare Formation of Pleistocene age. The model incorporates areally distributed ground-water recharge, areally and vertically distributed pumping, regional-collector drains in the Wesdands Water District (operative from 1980 to 1985), on-farm drains in parts of the Panoche, Broadview, and Firebaugh Water Districts, and bare-soil evaporation (which occurs if the water table is within 7 feet of land surface). The model also incorporates texture-based estimates of hydraulic conductivity, where texture is defined as the fraction of coarse-grained deposits present in a given subsurface interval. The numerical model was developed using hydrologic data from 1972 to 1988. Most of the parameters incorporated into the model were evaluated independently of the model, including system geometry, the distribution of texture, the altitudes of the water table and potentiometric surface of the confined zone in 1972 (initial condition), the hydraulic conductivity of coarse-grained deposits derived from the Coast Ranges, the hydraulic conductivity of coarse-grained deposits derived from the Sierra Nevada, specific storage, recharge, pumping, and parameters needed to incorporate drains and bare-soil evaporation. Four parameters were calibration variables: the hydraulic conductivity of fine-grained deposits in the semiconfined zone, the hydraulic conductivity of the Corcoran Clay Member, specific yield, and the transmissivity of the confined zone. The model was calibrated in two phases. In the first phase, a steady-state model of the ground-water flow system in 1984 was used to constrain the relation between the hydraulic conductivity of fine-grained deposits in the semiconfined zone and the hydraulic conductivity of the Corcoran Clay Member, thus reducing the number of independent variables from four to three. In the second phase of calibration, the change in altitude of the water table from 1972 to 1984, the change in altitude of the potentiometric surface of the confined zone from 1972 to 1984, and the number of model cells subject to bare-soil evaporation from 1972 to 1988 were used to evaluate the remaining three variables. The calibrated model reproduces the average change in water-table altitude (1972-84) to within 0.4 foot (average measured change 11.5 feet) and the average change in confined zone head (1972- 84) to within 19 feet (average measured change 120 feet). The simulated time-series record of the total number of model cells subject to bare-soil evaporation (each cell is 1 mile square) is within the range of the measured data. The measured values are at a minimum in October and a maximum in July. The October values ranged from 103 in 1972 to 132 in 1984 (the drains were closed in 1985) to 151 in 1988. The July values ranged from 144 in 1973 to 198 in 1984, to 204 in 1988. The simulated values ranged from 103 in 1972 to 161 in 1984, to 208 in 1988.

Water Supply Paper↗

Effects of nutrient enrichment on Prymnesium parvum population dynamics and toxicity: Results from field experiments, Lake Possum Kingdom, USA

Large fish kills associated with toxic populations of the haptophyte Prymnesium parvum occur worldwide. In the past 5 yr, incidences of P. parvum blooms in inland water bodies of Texas (USA) have increased dramatically, where cell densities in excess of 1 × 10 7 cells l –1 are typically observed. We conducted field experiments (Lake Possum Kingdom) during the fall and early spring of 28 d duration using 24 enclosures of 1.57 m 3 each. The experiments investigated the effect of nutrient enrichment, immigration of P. parvum and addition of barley straw extract on phytoplankton biomass and assemblage structure, P. parvum population density, zooplankton biomass and assemblage structure, bacteria, and toxicity. Nutrient enrichment stimulated P. parvum population growth beyond bloom proportions (>1 × 10 7 cells l –1 ). However, P. parvum did not dominate the assemblage under these conditions, as it does during natural blooms. Instead, euglenophytes and chlorophytes dominated. Toxicity, estimated using fish ( Pimephales promelas ) and cladoceran ( Daphnia magna ) bioassays and which is linked to P. parvum ’s allelopathic and mixotrophic effectiveness, was greatly reduced (eliminated in many cases) under conditions of nutrient enrichment. The suppression of toxicity by nutrient addition suggested that targeted and time-limited nutrient manipulations might be used to mitigate the effects of P. parvum blooms. Immigration of P. parvum into natural assemblages and addition of barley straw extract had no significant effect on plankton dynamics. Inter-Research 2007.

Texas↗

Gill tissue reactions in walleye Stizostedion vitreum vitreum and common carp Cyprinus carpio to glochidia of the freshwater mussel Lampsilis radiata siliquoidea

The glochidia of many freshwater mussels, which are obligate parasites on the gills, fins, and other body parts of specific fishes, attach to a suitable host, become encapsulated, and develop to the free-living juvenile stage. Using light and electron microscopy we compared gill tissue reactions in a suitable host (walleye Stizostedion vitreum vitreum ) and unsuitable host (common carp Cyprinus carpio ) infected with Lampsilis radiata siliquoidea . Encapsulation of glochidia on walleye gills was completed by 6 h post-infection at 20 to 22°C. Capsular formation and compaction were accompanied by a general increase in epithelioid cells. Fibrotic material appeared in capsules at about 48 h and virtually filled capsular cells from about Day 5 to Day 11 post-infection. Liberation of juvenile mussels was accompanied by thinning of the capsule from about Day 11 to Day l7. Although glochidia attached to the gills of common carp, few became encapsulated. By 48 h post-infection, preliminary capsular growth was evident and necrotic cells and cellular debris appeared at the edges of the growth. However, all glochidia were sloughed from carp gills by 60 h. Host specificity of L. radiata siliquoidea apparently depended on a combination of the attachment response of glochidia, differences in the encapsulation process, and tissue reactions in the fish.

Diseases of Aquatic Organisms↗

Characterization of a rhabdovirus isolated from carpione Salmo trutta carpio in Italy

A virus, strain 583, was isolated from carpione Salmo trutta carpio fry exhibiting high mortality. The virus was not neutralized by rabbit antisera against the fish rhabdoviruses viral haemorrhagic septicaemia virus (VHSV), infectious hematopoietic necrosis virus, eel rhabdovirus European X, spring viraemia of carp virus or pike fry rhabdovirus, or against the birnavirus infectious pancreatic necrosis virus. The virus replicated in several fish cell lines incubated at 20 to 25*C and grew optimally in the bluegill fry (BF-2) and fathead minnow (FHM) cell lines. Electron microscopy of infected BF-2 cell cultures revealed the presence of typical rhabdovirus particles, and immunofluorescent staining was observed using various polyclonal and monoclonal antibodies (MAbs) against Egtved virus, the causative agent of viral haemorrhagic septicaemia. The staining by a MAb against the nucleoprotein (N) of VHSV was particularly strong, a MAb against the glycoprotein (G) gave a moderate reaction, whereas a second MAb against the G protein and MAbs against the matrix proteins, M_(1) and M_(2), of VHSV did not react. Fluorescence titres using 3 rabbit antisera against whole Egtved virus varied between negative and moderately positive. Western blotting using polyclonal and monoclonal sera confirmed that both the N and G proteins of the carpione virus shared some epitopes with those of VHSV, but the M_(1) and M_(2) proteins did not. SDS-PAGE showed the structural proteins of the carpione virus produced a pattern typical of members of the Lyssavirus genus of the Rhabdoviridae and the molecular weights were very similar to those of VHSV, except for the M_(2) protein which was somewhat smaller. Infection trials showed the carpione virus induced high mortalities in carpione fry but not in rainbow trout Oncorhynchus mykiss fry. The carpione virus was clearly distinguishable from Egtved virus despite limited serological cross reaction. Since it was also easily distinguishable by immunofluorescence from the other fish rhabdoviruses included in the present study as well as in studies published elsewhere, it is concluded that the virus is a previously undescribed one. It is proposed that the virus be given the preliminary designation 'carpione rhabdovirus'.

Diseases of Aquatic Organisms↗

Isolation, characterization and molecular identification of a novel aquareovirus that infects the endangered fountain darter, Etheostoma fonticola

The fountain darter Etheostoma fonticola (FOD) is a federally endangered fish listed under the US Endangered Species Act. Here, we identified and characterized a novel aquareovirus isolated from wild fountain darters inhabiting the San Marcos River. This virus was propagated in Chinook salmon embryo (CHSE)-214, rainbow trout gonad-2 and fathead minnow cells at 15°C. The epithelioma papulosum cyprini cell line was refractory at all temperatures evaluated. High throughput sequencing technologies facilitated the complete genome sequencing of this virus utilizing ribosomal RNA-depleted RNA extracted from infected CHSE-214 cells. Conventional PCR primer sets were developed for the detection and confirmation of this virus to assist diagnostic screening methods. Phylogenetic analysis suggests this virus belongs to the Aquareovirus A genus. This research provides requisite initial data critical to support hatchery and refugia biosecurity measures for this endangered species.

Diseases of Aquatic Organisms↗

Variations of transcript profiles between sea otters Enhydra lutris from Prince William Sound, Alaska, and clinically normal reference otters

Development of blood leukocyte gene transcript profiles has the potential to expand condition assessments beyond those currently available to evaluate wildlife health, including sea otters Enhydra lutris , both individually and as populations. The 10 genes targeted in our study represent multiple physiological systems that play a role in immuno-modulation, inflammation, cell protection, tumor suppression, cellular stress-response, xenobiotic metabolizing enzymes, and antioxidant enzymes. These genes can be modified by biological, physical, or anthropogenic impacts and consequently provide information on the general type of stressors present in a given environment. We compared gene transcript profiles of sea otters sampled in 2008 among areas within Prince William Sound impacted to varying degrees by the 1989 ‘Exxon Valdez’ oil spill with those of captive and wild reference sea otters. Profiles of sea otters from Prince William Sound showed elevated transcription in genes associated with tumor formation, cell death, organic exposure, inflammation, and viral exposure when compared to the reference sea otter group, indicating possible recent and chronic exposure to organic contaminants. Sea otters from historically designated oiled areas within Prince William Sound 19 yr after the oil spill had higher transcription of genes associated with tumor formation, cell death, heat shock, and inflammation than those from areas designated as less impacted by the spill.

Alaska↗

Microbial aggregates within tissues infect a diversity of corals throughout the Indo-Pacific

Coral reefs are highly diverse ecosystems where symbioses play a pivotal role. Corals contain cell-associated microbial aggregates (CAMA), yet little is known about how widespread they are among coral species or the nature of the symbiotic relationship. Using histology, we found CAMA within 24 species of corals from 6 genera from Hawaii, American Samoa, Palmyra, Johnston Atoll, Guam, and Australia. Prevalence (%) of infection varied among coral genera: Acropora, Porites, and Pocillopora were commonly infected whereas Montipora were not. Acropora from the Western Pacific were significantly more likely to be infected with CAMA than those from the Central Pacific, whereas the reverse was true for Porites. Compared with apparently healthy colonies, tissues from diseased colonies were significantly more likely to have both surface and basal body walls infected. The close association of CAMA with host cells in numerous species of apparently healthy corals and lack of associated cell pathology reveals an intimate agent-host association. Furthermore, CAMA are Gram negative and in some corals may be related to chlamydia or rickettsia. We propose that CAMA in adult corals are facultative secondary symbionts that could play an important ecological role in some dominant coral genera in the Indo-Pacific. CAMA are important in the life histories of other animals, and more work is needed to understand their role in the distribution, evolution, physiology, and immunology of reef corals.

Marine Ecology Progress Series↗