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 901 records · Page 50Linked to original sources
USGS1996

Use of a ground-water flow model with particle tracking to evaluate ground-water vulnerability, Clark County, Washington

A ground-water flow model was used in conjunction with particle tracking to evaluate ground-water vulnerability in Clark County, Washington. Using the particle-tracking program, particles were placed in every cell of the flow model (about 60,000 particles) and tracked backwards in time and space upgradient along flow paths to their recharge points. A new computer program was developed that interfaces the results from a particle-tracking program with a geographic information system (GIS). The GIS was used to display and analyze the particle-tracking results. Ground-water vulnerability was evaluated by selecting parts of the ground-water flow system and combining the results with ancillary information stored in the GIS to determine recharge areas, characteristics of recharge areas, downgradient impact of land use at recharge areas, and age of ground water. Maps of the recharge areas for each hydrogeologic unit illustrate the presence of local, intermediate, or regional ground-water flow systems and emphasize the three-dimensional nature of the ground-water flow system in Clark County. Maps of the recharge points for each hydrogeologic unit were overlaid with maps depicting aquifer sensitivity as determined by DRASTIC (a measure of the pollution potential of ground water, based on the intrinsic characteristics of the near-surface unsaturated and saturated zones) and recharge from on-site waste-disposal systems. A large number of recharge areas were identified, particularly in southern Clark County, that have a high aquifer sensitivity, coincide with areas of recharge from on-site waste-disposal systems, or both. Using the GIS, the characteristics of the recharge areas were related to the downgradient parts of the ground-water system that will eventually receive flow that has recharged through these areas. The aquifer sensitivity, as indicated by DRASTIC, of the recharge areas for downgradient parts of the flow system was mapped for each hydrogeologic unit. A number of public-supply wells in Clark County may be receiving a component of water that recharged in areas that are more conducive to contaminant entry. The aquifer sensitivity maps illustrate a critical deficiency in the DRASTIC methodology: the failure to account for the dynamics of the ground-water flow system. DRASTIC indices calculated for a particular location thus do not necessarily reflect the conditions of the ground-water resources at the recharge areas to that particular location. Each hydrogeologic unit was also mapped to highlight those areas that will eventually receive flow from recharge areas with on-site waste-disposal systems. Most public-supply wells in southern Clark County may eventually receive a component of water that was recharged from on-site waste-disposal systems.Traveltimes from particle tracking were used to estimate the minimum and maximum age of ground water within each model-grid cell. Chlorofluorocarbon (CFC)-age dating of ground water from 51 wells was used to calibrate effective porosity values used for the particle- tracking program by comparison of ground-water ages determined through the use of the CFC-age dating with those calculated by the particle- tracking program. There was a 76 percent agreement in predicting the presence of modern water in the 51 wells as determined using CFCs and calculated by the particle-tracking program. Maps showing the age of ground water were prepared for all the hydrogeologic units. Areas with the youngest ground-water ages are expected to be at greatest risk for contamination from anthropogenic sources. Comparison of these maps with maps of public- supply wells in Clark County indicates that most of these wells may withdraw ground water that is, in part, less than 100 years old, and in many instances less than 10 years old. Results of the analysis showed that a single particle-tracking analysis simulating advective transport can be used to evaluate ground-water vulnerability for any part of a ground-wate

Open-File Report↗
USGS1993

Mineralogy, mineral chemistry, and paragenesis of gold, silver, and base-metal ores of the North Amethyst vein system, San Juan Mountains, Mineral County, Colorado

Gold-rich adularia-sericite-type mineralization occurs near the southern margin of the San Luis caldera, at the intersection of the Equity fault and the northern extension of the Amethyst fault system. Mineralized rock is confined primarily to steeply dipping structures in silicified rhyolite and dacite. Intense sericitic alteration occurs at higher levels in the vein system, and wall rock adjacent to some veins is bleached. The ores are multiply brecciated, and vein filling locally shows sedimentary textures. Textural, mineralogical, and chemical criteria indicate that there are at least two partially coextensive associations of mineral assemblages separated by periods of brecciation and sedimentation. An older gold-bearing association consists of two fine-grained ore stages, both of which contain electrum, uytenbogaardtite, tetrahedrite, silver sulfosalts, silver sulfides and base-metal sulfides, and a manganese-rich stage containing the assemblages (1) manganese silicate + manganese carbonate minerals + quartz and (2) magnetite + hematite + pyrite + quartz. A younger crosscutting association contains calcite, adularia, fluorite, and quartz, plus the assemblages (1) coarse-grained basemetal sulfides and (2) hematite + chlorite + quartz. Quartz, manganese-rich calcite, and trace pyrite line late-stage vugs. Mineralogic, lead-isotopic, and fluid-inclusion characteristics of the younger association are similar to those of ores of the southern and central parts of the Creede mining district. In contrast, the gold and manganese-silicate assemblages of the older association are rare to absent in the southern and central parts of the district. The local and early occurrence of the manganese and gold assemblages may indicate that they formed in a small hydrothermal cell that predated the extensive hydrothermal system from which ores of the central and southern parts of the Creede district are proposed to have been deposited (Bethke, 1988). If similar early-stage cells were present in the southern and central parts of the district, they may have been replaced or overprinted by later assemblages, and they may remain to be discovered. In the latter case, mineral assemblages that formed at early stages in the paragenesis hold the most promise for gold exploration.

Colorado↗
USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗
USGS2017

Beryllium

Beryllium is a mineral commodity that is used in a variety of industries to make products that are essential for the smooth functioning of a modern society. Two minerals, bertrandite (which is supplied domestically) and beryl (which is currently supplied solely by imports), are necessary to ensure a stable supply of high-purity beryllium metal, alloys, and metal-matrix composites and beryllium oxide ceramics. Although bertrandite is the source mineral for more than 90 percent of the beryllium produced globally, industrial beryl is critical for the production of the very high purity beryllium metal needed for some strategic applications. The current sole domestic source of beryllium is bertrandite ore from the Spor Mountain deposit in Utah; beryl is imported mainly from Brazil, China, Madagascar, Mozambique, and Portugal. High-purity beryllium metal is classified as a strategic and critical material by the Strategic Materials Protection Board of the U.S. Department of Defense because it is used in products that are vital to national security. Beryllium is maintained in the U.S. stockpile of strategic materials in the form of hot-pressed beryllium metal powder. Because of its unique chemical properties, beryllium is indispensable for many important industrial products used in the aerospace, computer, defense, medical, nuclear, and telecommunications industries. For example, high-performance alloys of beryllium are used in many specialized, high-technology electronics applications, as they are energy efficient and can be used to fabricate miniaturized components. Beryllium-copper alloys are used as contacts and connectors, switches, relays, and shielding for everything from cell phones to thermostats, and beryllium-nickel alloys excel in producing wear-resistant and shape-retaining high-temperature springs. Beryllium metal composites, which combine the fabrication ability of aluminum with the thermal conductivity and highly elastic modulus of beryllium, are ideal for producing aircraft and satellite structural components that have a high stiffness-to-weight ratio and low surface vibration. Beryllium oxide ceramics are used in a wide range of applications, including missile guidance systems, radar applications, and cell phone transmitters, and they are critical to medical technologies, such as magnetic resonance imaging (MRI) machines, medical lasers, and portable defibrillators. The United States is expected to remain self-sufficient with respect to most of its beryllium requirements, based on information available at the time this chapter was prepared (2013). The United States is one of only three countries that currently process beryllium ores and concentrate them into beryllium products, and these three countries supply most of the rest of the world with these products. Exploration for new deposits in the United States is limited because domestic beryllium production is dominated by a single producer that effectively controls the domestic beryllium market, which is relatively small and specialized, and the market cannot readily accommodate new competition on the raw material supply side.

Professional Paper↗
USGS2017

Gallium

Gallium is a soft, silvery metallic element with an atomic number of 31 and the chemical symbol Ga. Gallium is used in a wide variety of products that have microelectronic components containing either gallium arsenide (GaAs) or gallium nitride (GaN). GaAs is able to change electricity directly into laser light and is used in the manufacture of optoelectronic devices (laser diodes, light-emitting diodes [LEDs], photo detectors, and solar cells), which are important for aerospace and telecommunications applications and industrial and medical equipment. GaAs is also used in the production of highly specialized integrated circuits, semiconductors, and transistors; these are necessary for defense applications and high-performance computers. For example, cell phones with advanced personal computer-like functionality (smartphones) use GaAs-rich semiconductor components. GaN is used principally in the manufacture of LEDs and laser diodes, power electronics, and radio-frequency electronics. Because GaN power transistors operate at higher voltages and with a higher power density than GaAs devices, the uses for advanced GaN-based products are expected to increase in the future. Gallium technologies also have large power-handling capabilities and are used for cable television transmission, commercial wireless infrastructure, power electronics, and satellites. Gallium is also used for such familiar applications as screen backlighting for computer notebooks, flat-screen televisions, and desktop computer monitors. Gallium is dispersed in small amounts in many minerals and rocks where it substitutes for elements of similar size and charge, such as aluminum and zinc. For example, gallium is found in small amounts (about 50 parts per million) in such aluminum-bearing minerals as diaspore-boehmite and gibbsite, which form bauxite deposits, and in the zinc-sulfide mineral sphalerite, which is found in many mineral deposits. At the present time, gallium metal is derived mainly as a byproduct of the processing of bauxite ore for aluminum; lesser amounts of gallium metal are produced from the processing of sphalerite ore from three types of deposits (sediment-hosted, Mississippi Valley-type, and volcanogenic massive sulfide) for zinc. The United States is expected to meet its current and expected future needs for gallium through imports of primary, recycled, and refined gallium, as well as through domestic production of recycled and refined gallium. The U.S. Geological Survey estimates that world resources of gallium in bauxite exceed 1 billion kilograms, and a considerable quantity of gallium could be present in world zinc reserves.

Professional Paper↗
USGS2017

Niobium and tantalum

Niobium and tantalum are transition metals that are almost always found together in nature because they have very similar physical and chemical properties. Their properties of hardness, conductivity, and resistance to corrosion largely determine their primary uses today. The leading use of niobium (about 75 percent) is in the production of high-strength steel alloys used in pipelines, transportation infrastructure, and structural applications. Electronic capacitors are the leading use of tantalum for high-end applications, including cell phones, computer hard drives, and such implantable medical devices as pacemakers. Niobium and tantalum are considered critical and strategic metals based on the potential risks to their supply (because current production is restricted to only a few countries) and the significant effects that a restriction in supply would have on the defense, energy, high-tech industrial, and medical sectors. The average abundance of niobium and tantalum in bulk continental crust is relatively low—8.0 parts per million (ppm) niobium and 0.7 ppm tantalum. Their chemical characteristics, such as small ionic size and high electronic field strength, significantly reduce the potential for these elements to substitute for more common elements in rock-forming minerals and make niobium and tantalum essentially immobile in most aqueous solutions. Niobium and tantalum do not occur naturally as pure metals but are concentrated in a variety of relatively rare oxide and hydroxide minerals, as well as in a few rare silicate minerals. Niobium is primarily derived from the complex oxide minerals of the pyrochlore group ((Na,Ca,Ce) 2 (Nb,Ti,Ta) 2 (O,OH,F) 7 ), which are found in some alkaline granite-syenite complexes (that is, igneous rocks containing sodium- or potassium-rich minerals and little or no quartz) and carbonatites (that is, igneous rocks that are more than 50 percent composed of primary carbonate minerals, by volume). Tantalum is derived mostly from the mineral tantalite ((Fe,Mn)(Ta,Nb) 2 O 6 ), which is found as an accessory mineral in rare-metal granites and pegmatites that are also enriched in lithium and cesium (termed lithium-cesium-tantalum (LCT)-type pegmatites). Brazil and Canada are the leading nations that produce niobium mineral concentrates, but Brazil is by far the leading producer, accounting for about 90 percent of production, which comes mostly from weathered material derived from carbonatites. Brazil and Canada also have the largest identified niobium resources; additional resources, although they are less well reported, occur in Angola, Australia, China, Greenland, Malawi, Russia, and South Africa. Australia and Brazil have been the leading producers of tantalum mineral concentrates, although recently Ethiopia and Mozambique have also been significant suppliers of tantalum. Artisanal mining of columbite-tantalite (also called coltan) is practiced in many countries, particularly Burundi, the Democratic Republic of the Congo (Congo [Kinshasa]), Nigeria, Rwanda, and Uganda. Brazil has about 40 percent of the identified tantalum resources; other countries and regions with identified tantalum resources include, in decreasing order of resources, Australia, Asia, Russia and the Middle East, Africa, North America, and Europe. Identified niobium and tantalum resources in the United States are small, low grade, and difficult to recover and process, and are thus not commercially recoverable at current prices. Consequently, the United States meets its current and expected future needs for niobium and tantalum through imports of primary mineral concentrates and alloys and through recovery from foreign and domestic alloy scrap that contain the metals. Environmentally, the main issues related to niobium and tantalum mining are land disruptions, the volume of waste materials and their disposal, and the radioactivity of some tailings and waste materials that contain thorium and uranium. Because of the relative biological inertness of niobium and tantalum, human and ecological health concerns are generally minimal under most natural conditions. Demand for both niobium and tantalum is expected to increase as the world economy continues to recover from the downturn that began in 2008. Increased demand for niobium is linked to increased consumption of microalloyed steel, which is used in the manufacture of cars, buildings, ships, and refinery equipment. Demand for these steels will likely increase with continued economic development in such countries as Brazil, China, and India. In addition, increased global demand for cars, cell phones, computers, superconducting magnets, and other high-tech devices will likely spur increased demand for both niobium and tantalum. The estimated global reserves and resources of niobium and tantalum are large and appear more than sufficient to meet global demand for the foreseeable future, possibly the next 500 years. The sale of “conflict coltan” attributed to rebel forces waging a civil war in Congo (Kinshasa) has been of recent concern and has highlighted the need for a transparent and traceable global supply chain that can exclude illegal columbite-tantalite from the conventional market while discerning legitimate artisanal mine production in central Africa.

Professional Paper↗
USGS2009

Terrestrial Ecosystems - Land Surface Forms of the Conterminous United States

As part of an effort to map terrestrial ecosystems, the U.S. Geological Survey has generated land surface form classes to be used in creating maps depicting standardized, terrestrial ecosystem models for the conterminous United States, using an ecosystems classification developed by NatureServe . A biophysical stratification approach, developed for South America and now being implemented globally, was used to model the ecosystem distributions. Since land surface forms strongly influence the differentiation and distribution of terrestrial ecosystems, they are one of the key input layers in this biophysical stratification. After extensive investigation into various land surface form mapping methodologies, the decision was made to use the methodology developed by the Missouri Resource Assessment Partnership (MoRAP). MoRAP made modifications to Hammond's land surface form classification, which allowed the use of 30-meter source data and a 1-km2 window for analyzing the data cell and its surrounding cells (neighborhood analysis). While Hammond's methodology was based on three topographic variables, slope, local relief, and profile type, MoRAP's methodology uses only slope and local relief. Using the MoRAP method, slope is classified as gently sloping when more than 50 percent of the area in a 1-km2 neighborhood has slope less than 8 percent, otherwise the area is considered moderately sloping. Local relief, which is the difference between the maximum and minimum elevation in a neighborhood, is classified into five groups: 0-15 m, 16-30 m, 31-90 m, 91-150 m, and >150 m. The land surface form classes are derived by combining slope and local relief to create eight landform classes: flat plains (gently sloping and local relief =< 15 m), smooth plains (gently sloping and 15 m < local relief =< 30 m), irregular plains (gently sloping and 30 m < local relief =< 90 m), escarpments (gently sloping and local relief > 90 m), low hills (not gently sloping and local relief =< 30 m), hills (not gently sloping and 30 m < local relief =< 90 m), breaks/foothills (not gently sloping and 90 m < local relief =< 150 m), and low mountains (not gently sloping and local relief > 150 m). However, in the USGS application of the MoRAP methodology, an additional local relief group was used (> 400 m) to capture additional local topographic variation. As a result, low mountains were redefined as not gently sloping and 151 m < local relief < 400 m, and a new land surface form class, high mountains/deep canyons, was identified as not gently sloping and local relief > 400 m. The final application of the MoRAP methodology was implemented using the USGS 30-meter National Elevation Dataset and an existing USGS slope dataset that had been derived by calculating the slope from the NED in Universal Transverse Mercator (UTM) coordinates in each UTM zone, and then combining all of the zones into a national dataset. This map shows a smoothed image of the nine land surface form classes based on MoRAP's methodology. Additional information about this map and any data developed for the ecosystems modeling of the conterminous United States is available online at http://rmgsc.cr.usgs.gov/ecosystems/.

Scientific Investigations Map↗
USGS2009

Terrestrial Ecosystems - Topographic Moisture Potential of the Conterminous United States

As part of an effort to map terrestrial ecosystems, the U.S. Geological Survey has generated topographic moisture potential classes to be used in creating maps depicting standardized, terrestrial ecosystem models for the conterminous United States, using an ecosystems classification developed by NatureServe. A biophysical stratification approach, developed for South America and now being implemented globally, was used to model the ecosystem distributions. Substrate moisture regimes strongly influence the differentiation and distribution of terrestrial ecosystems, and therefore topographic moisture potential is one of the key input layers in this biophysical stratification. The method used to produce these topographic moisture potential classes was based on the derivation of ground moisture potential using a combination of computed topographic characteristics (CTI, slope, and aspect) and mapped National Wetland Inventory (NWI) boundaries. This method does not use climate or soil attributes to calculate relative topographic moisture potential since these characteristics are incorporated into the ecosystem model though other input layers. All of the topographic data used for this assessment were derived from the USGS 30-meter National Elevation Dataset (NED ) including the National Compound Topographic Index (CTI). The CTI index is a topographically derived measure of slope for a raster cell and the contributing area from upstream raster cells, and thus expresses potential for water flow to a point. In other words CTI data are 'a quantification of the position of a site in the local landscape', where the lowest values indicate ridges and the highest values indicate stream channels, lakes and ponds. These CTI values were compared to independent estimates of water accumulation by obtaining geospatial data from a number of sample locations representing two types of NWI boundaries: freshwater emergent wetlands and freshwater forested/shrub wetlands. Where these shorelines (the interface between the NWI wetlands and adjacent land) occurred, the CTI values were extracted and a histogram of their statistical distributions was calculated. Based on an evaluation of these histograms, CTI thresholds were developed to separate periodically saturated or flooded land, mesic uplands (moderately moist), and uplands. After the range of CTI values for these three different substrate moisture regimes was determined, the CTI values were grouped into three initial topographic moisture potential classes. As a final step in the generation of this national data layer, the uplands classification was subdivided into either very dry uplands or dry uplands. Very dry uplands were defined as uplands with relatively steep, south-facing slopes, and identification of this class was based on the slope and aspect datasets derived from the NED. The remaining uplands that did not meet these additional criteria were simply re-classified as dry uplands. The final National Topographic Moisture Potential dataset for the conterminous United States contains four classes: periodically saturated or flooded land (CTI = 18.5), mesic uplands (12 =< CTI < 18.5), dry uplands (CTI < 12), and very dry uplands (CTI < 12, Slope > 24 degrees and 91 degrees =< Aspect =< 314 degrees). This map shows a smoothed and generalized image of the four topographic moisture potential classes. Additional information about this map and any of the data developed for the ecosystems modeling of the conterminous United States is available online at http://rmgsc.cr.usgs.gov/ecosystems/.

Scientific Investigations Map↗
USGS2012

California State Waters Map Series — Hueneme Canyon and vicinity, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within the 3-nautical-mile limit of California's State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Hueneme Canyon and vicinity map area lies within the eastern Santa Barbara Channel region of the Southern California Bight. The area is part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation - at least 90&deg; - since the early Miocene has been proposed for the Western Transverse Ranges, and the region is presently undergoing north-south shortening. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area, which is offshore of the Oxnard plain and west of and along the trend of the south flank of the Santa Monica Mountains, lies at the east end of the Santa Barbara littoral cell, characterized by west-to-east littoral transport of sediment derived mainly from coastal watersheds. The Hueneme Canyon and vicinity map area in California's State Waters is characterized by two major physiographic features: (1) the nearshore continental shelf, and (2) the Hueneme and Mugu Submarine Canyon system, which, in the map area, includes Hueneme Canyon and parts of three smaller, unnamed headless canyons incised into the shelf southeast of Hueneme Canyon. The shelf is underlain by tens of meters of interbedded upper Quaternary shelf, estuarine, and fluvial deposits that formed as sea level fluctuated in the last several hundred thousand years. Hueneme Canyon extends about 15 km offshore from its canyon head near the dredged navigation channel of the Port of Hueneme. The canyon is relatively deep (about 150 m at the California's State Waters limit) and steep (canyon flanks as steep as 25&deg; to 30&deg;). Historically, Hueneme Canyon functioned as the eastern termination of the Santa Barbara littoral cell by trapping all eastward littoral drift, not only feeding the large Hueneme submarine fan but acting as the major conduit of sediment to the deep Santa Monica Basin; however, recent dredging programs needed to maintain Channel Islands Harbor and the Port of Hueneme have moved the nearshore sediment trapped by jetties and breakwaters to an area southeast of the Hueneme Canyon head. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft sediment and isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Hueneme Canyon and vicinity map area are related directly to the geomorphology and sedimentary processes that are the result of its Quaternary geologic history. The two basic megahabitats in the map area are Shelf (continental shelf) and Flank (continental slope). The flat seafloor of the continental shelf in the Hueneme Canyon and vicinity map area is dynamic, as indicated by mobile sand sheets and coarser grained scour depressions. The active Hueneme Canyon provides considerable relief to the continental shelf in the map area, and its irregular morphology of eroded walls, landslide scarps, and deposits and gullies provide promising habitat for groundfish, crabs, shrimp, and other marine benthic organisms. Most invertebrates observed in the map area during camera ground-truth field operations are found on the edge of Hueneme Canyon, which may be an important area of recruitment and retention to other invertebrates and fishes. The smaller, more subtle, nonactive headless canyons located primarily on the continental slope also offer relief that provides habitat for groundfish and other organisms.

California↗
USGS2015

California State Waters Map Series: Offshore of Refugio Beach, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California&rsquo;s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Refugio Beach map lies within the western Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation&mdash;at least 90&deg;&mdash;since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 0.5 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The coastal zone of this map area lies at the steep flank of the Santa Ynez Mountains. The crest of the range, which lies about 8 km from the shoreline (north of the map area), has a maximum elevation of about 780 m. This area is largely open space, partly used for livestock grazing, with no significant towns or population centers. Highway 101 crosses the map area, adjacent to and within a few hundred meters of the shoreline. The most significant developments are the recreational state beaches at El Capitan Beach and Refugio Beach. The beaches are subject to erosion each winter during storm-wave attack, and then they undergo gradual recovery or accretion during the more gentle wave climate of the late spring, summer, and fall months. The Offshore of Refugio Beach map area lies in the west-central part of the Santa Barbara littoral cell, which is characterized by west-to-east transport of sediment from Point Arguello on the northwest to Hueneme and Mugu Canyons on the southeast. Longshore drift rates have been reported to range from about 160,000 to 800,000 tons/yr, averaging 400,000 tons/yr. Sediment supply to the western and central part of the littoral cell, including the map area, is mainly from relatively small coastal watersheds. Within the map area, these coastal watersheds include (from east to west) Ca&ntilde;ada del Capitan, Tajiquas Creek, Arroyo Hondo, Ca&ntilde;ada del Molino, and several unnamed canyons and creeks. The Santa Ynez and Santa Maria Rivers, the mouths of which are 80 to 120 km northwest of the map area, are not considered significant sediment sources because Point Conception and Point Arguello provide obstacles to downcoast sediment transport and also because much of their sediment load is trapped in dams. The Ventura and Santa Clara Rivers, the mouths of which are about 70 km to the southeast of the map area, are not sediment sources for the map area. The offshore part of the map area consists of relatively flat and shallow continental shelf, which dips gently seaward (about 0.8&deg; to 1.0&deg;) so that water depths at the shelf break, roughly coincident with the California&rsquo;s State Waters limit, are about 80 to 100 m. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary marine and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. In the map area, the shelf break is at depths of about 90 m and lies about 5.6 to 6.4 km offshore. Beyond the shelf break, the slope is steep (as much as about 7&deg;) and unstable. Several submarine landslides, including the large (130 km 2 ) Goleta landslide complex, have been documented offshore of Goleta, a few kilometers east of the map area. This compound slump complex may have been initiated more than 200,000 year ago, but it also includes three recent failures that may have been generated 8,000 to 10,000 years ago. A local, 10-m-high tsunami may have been generated from these failure events. Small folds related to local faulting are superimposed on the homocline that makes up the south flank of the Santa Ynez Mountains. One of these superimposed anticlines hosts the Molino gas field, which was discovered in 1962 and subsequently developed through directional drilling from onshore wells. The Oligocene Sespe and Vaqueros Formations are the reservoirs in the Molino gas field, and the map area includes numerous seafloor hydrocarbon seeps. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft, unconsolidated sediment interspersed with isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Refugio Beach map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats, primarily are composed of soft sediment interrupted by a few carbonate mounds. This homogeneous seafloor of sediment and low-relief bedrock provides promising habitat for groundfish, crabs, shrimp, and other marine benthic organisms.

California↗
USGS2004

Simulation of ground-water flow in the Cedar River alluvial aquifer flow system, Cedar Rapids, Iowa

The Cedar River alluvial aquifer is the primary source of municipal water in the Cedar Rapids, Iowa, area. Since 1992, the U.S. Geological Survey, in cooperation with the City of Cedar Rapids, has investigated the hydrogeology and water quality of the Cedar River alluvial aquifer. This report describes a detailed analysis of the ground-water flow system in the alluvial aquifer, particularly near well field areas. The ground-water flow system in the Cedar Rapids area consists of two main components, the unconsolidated Quaternary deposits and the underlying carbonate bedrock that has a variable fracture density. Quaternary deposits consist of eolian sand, loess, alluvium, and glacial till. Devonian and Silurian bedrock aquifers overlie the Maquoketa Shale (Formation) of Ordovician age, a regional confining unit. Ground-water and surface-water data were collected during the study to better define the hydrogeology of the Cedar River alluvial aquifer and Devonian and Silurian aquifers. Stream stage and discharge, ground-water levels, and estimates of aquifer hydraulic properties were used to develop a conceptual ground-water flow model and to construct and calibrate a model of the flow system. This model was used to quantify the movement of water between the various components of the aluvial aquifer flow system and provide an improved understanding of the hydrology of the alluvial aquifer. Ground-water flow was simulated for the Cedar River alluvial aquifer and the Devonian and Silurian aquifers using the three-dimensional finite-difference ground-water flow model MODFLOW. The model was discretized into 223 rows and 354 columns of cells. Areal cell sizes range from about 50 feet on a side near the Cedar River and the Cedar Rapids municipal wells to 1,500 feet on a side near the model boundaries and farthest away from the Cedar Rapids municipal well fields. The model is separated into five layers to account for the various hydrogeologic units in the model area. Model results indicate that the primary sources of inflow to the modeled area are infiltration from the Cedar River (53.0 percent) and regional flow in the glacial and bedrock materials (34.1 percent). The primary sources of outflow from the modeled area are discharge to the Cedar River (45.4 percent) and pumpage (44.8 percent). Current steady-state pumping rates have increased the flow of water from the Cedar River to the alluvial aquifer by 43.8 cubic feet per second. Steady-state and transient hypothetical pumpage scenarios were used to show the relation between changes in pumpage and changes in infiltration of water from the Cedar River. Results indicate that more than 99 percent of the water discharging from municipal wells infiltrates from the Cedar River, that the time required for induced river recharge to equilibrate with municipal pumpage may be 150 days or more, and that ground-water availability in the Cedar Rapids area will not be significantly affected by doubling current pumpage as long as there is sufficient flow in the Cedar River to provide recharge.

Iowa↗
USGS2004

Hydrogeology and Hydrologic Landscape Regions of Nevada

In 1999, the U.S. Environmental Protection Agency initiated a rule to protect ground water in areas other than source-water protection areas. These other sensitive ground water areas (OSGWAs) are aquifers that are not currently but could eventually be used as a source of drinking water. The OSGWA program specifically addresses existing wells that are used for underground injection of motor vehicle waste. If the injection well is in a ground-water protection area or an OSGWA, well owners must either close the well or apply for a permit. The Nevada Division of Environmental Protection will evaluate site-specific information and determine if the aquifer associated with a permit application is susceptible to contamination. A basic part of evaluating OSGWAs is characterizing the hydrogeology of aquifer systems including the lithology, hydrologic properties, soil permeability, and faulting, which partly control the susceptibility of ground water to contamination. Detailed studies that evaluate ground-water susceptibility are not practical in a largely unpopulated State like Nevada. However, existing and new information could be extrapolated to other areas of the State if there is an objective framework to transfer the information. The concept of hydrologic landscape regions, which identify areas with similar hydrologic characteristics, provides this framework. This report describes the hydrogeology and hydrologic landscape regions of Nevada. Consolidated rocks that form mountain ranges and unconsolidated sediments that fill the basins between the ranges are grouped into hydrogeologic units having similar lithology and assumed to have similar hydrologic properties. Consolidated rocks and unconsolidated sediments are the two major hydrogeologic units and comprise 51 and 49 percent of the State, respectively. Consolidated rocks are subdivided into 8 hydrogeologic units. In approximate order of decreasing horizontal hydraulic conductivity, consolidated-rock hydrogeologic units consist of: (1) carbonate rocks, Quaternary to Tertiary age; (2) basaltic, (3) rhyolitic, and (4) andesitic volcanic flows; (5) volcanic breccias, tuffs, and volcanic rocks older than Tertiary age; (6) intrusive and metamorphic rocks; (7) consolidated and semi-consolidated tuffaceous rocks and sediments; and (8) clastic rocks consisting of sandstone and siltstone. Unconsolidated sediments are subdivided into four hydrogeologic units on the basis of flow regime, topographic slope, and mapped stream channels. The four units are (1) alluvial slopes, (2) valley floors, (3) fluvial deposits, and (4) playas. Soil permeability was grouped into five descriptive categories ranging from very high to very low, which generally correspond to mapped geomorphic features such as playas and alluvial slopes. In general, soil permeability is low to moderate in northern, northeastern, and eastern Nevada and high to very high in western, southwestern, and southern Nevada. Within a particular basin, soil permeability decreases downslope from the bedrock contact. The type of parent rock, climate, and streamflow velocities are factors that likely cause these spatial patterns. Faults in unconsolidated sediments usually are barriers to ground-water flow. In consolidated rocks, permeability and ground-water flow is reduced in directions normal to the fault zone and increased in directions parallel to the fault zone. With time, mineral precipitation may seal fractures in consolidated rocks, reducing the permeability. However, continued movement along the fault may form new fractures, resulting in a fault alternating from a zone of preferred flow to a flow barrier during geologic time. The effect of faults on ground-water flow at a particular location is difficult to determine without a site- specific investigation. Hydrologic landscape regions were delineated by overlaying a grid of 100-foot (30-meter) cells over the State, estimating the value of five variables for each cell, an

Nevada↗
USGS2004

Fractionation and characterization of organic matter in wastewater from a swine waste-retention basin

Organic matter in wastewater sampled from a swine waste-retention basin in Iowa was fractionated into 14 fractions on the basis of size (particulate, colloid, and dissolved); volatility; polarity (hydrophobic, transphilic, hydrophilic); acid, base, neutral characteristics; and precipitate or flocculates (floc) formation upon acidification. The compound-class composition of each of these fractions was determined by infrared and 13C-NMR spectral analyses. Volatile acids were the largest fraction with acetic acid being the major component of this fraction. The second most abundant fraction was fine particulate organic matter that consisted of bacterial cells that were subfractionated into extractable lipids consisting of straight chain fatty acids, peptidoglycans components of bacterial cell walls, and protein globulin components of cellular plasma. The large lipid content of the particulate fraction indicates that non-polar contaminants, such as certain pharmaceuticals added to swine feed, likely associate with the particulate fraction through partitioning interactions. Hydrocinnamic acid is a major component of the hydrophobic acid fraction, and its presence is an indication of anaerobic degradation of lignin originally present in swine feed. This is the first study to combine particulate organic matter with dissolved organic matter fractionation into a total organic matter fractionation and characterization.

Scientific Investigations Report↗
USGS2004

Environmental factors and chemical and microbiological water-quality constituents related to the presence of enteric viruses in ground water from small public water supplies in southeastern Michigan

A study of small public ground-water-supply wells that produce water from discontinuous sand and gravel aquifers was done from July 1999 through July 2001 in southeastern Michigan. Samples were collected to determine the occurrence of viral pathogens and microbiological indicators of fecal contamination (“indicators”), determine whether indicators are adequate predictors of the presence of enteric viruses, and determine the factors that affect the presence of enteric viruses. Small systems are those that serve less than 3,300 people. Samples were analyzed for specific enteric viruses by reverse transcriptase-polymerase chain reaction (RT-PCR), for culturable viruses by cell culture, and for the indicators total coliforms, Escherichia coli (E. coli) , enterococci, and F-specific and somatic coliphage. Ancillary environmental and water-quality data were collected or compiled. A total of 169 regular samples and 32 replicate pairs were collected from 38 wells. Replicate pairs were samples collected at the same well on the same date. One well was sampled 6 times, 30 wells were sampled five times, 6 wells were sampled twice, and 1 well was sampled once. By use of RT-PCR, enterovirus was found in four wells (10.5 percent) and hepatitis A virus (HAV) in five wells (13.2 percent). In two of these wells, investigators found both enterovirus and HAV, but on different sampling dates. Culturable viruses were found one time in two wells (5.9 percent), and neither of these wells was positive for viruses by use of RT-PCR on any sampling date. If results for all viruses are combined, 9 of the 38 small public-supply wells were positive for enteric viruses (23.7 percent) by either cell culture or RT-PCR. One or more indicators were found in 18 of 38 wells. Total coliforms, E. coli , enterococci, and F-specific and somatic coliphage were found in 34.2, 10.5, 15.8, 5.9, and 5.9 percent, respectively, of the wells tested. In only 3 out of 18 wells were samples positive for an indicator on more than one date at the same well. The co-occurrence of enteric viruses and any indicator was 55.6 percent; five out of the nine virus-positive wells were also found to be positive for an indicator. Two wells with detections of viruses had a detection of total coliforms, one well had a detection of E. coli , one of enterococci, and one of F-specific coliphage. On a per sample basis, of 11 samples that were positive for enteric viruses, indicator bacteria co-occurred in only 2 samples, and coliphage were not present in any. More virus-positive samples were found at sites served by septic systems than those served by sewerlines. Sampling condition (ground water or a mixture of tank and ground water), distance to septic system, type of and distance to nearest surface-water body, well characteristics, or land use were not related to the presence of viruses or indicators. Among continuous water-quality variables, statistically significant relations were found between total coliforms and dissolved organic carbon and between total coliforms and iron. There was a statistically significant relation between chloride concentrations >20 mg/L and detections of total coliforms. Presence of nitrate and nitrite was related to the presence of other indicators ( E. coli , enterococci, and F-specific and somatic coliphage) or enteric viruses, but not to total coliforms. The data indicated that chloride-to-bromide (C1:Br) ratios may be useful as a screening tool for total coliforms and enteric viruses but not for E. coli , enterococci, and F-specific and somatic coliphage. This study provides evidence for fecal contamination of ground water from small public-supply wells, at least on an intermittent basis. Collecting data on multiple lines of evidence would be needed to reliably predict the presence of enteric viruses and protect public health. Future data collection toward this end could include repeat sampling several times a year for different indicators, measuring dissolved-organic carbon, nitrate plus nitrite, and (or) chloride concentrations, or determining C1:Br ratios. The presence of a site served by a septic system is an indication that the well may be more vulnerable to contamination than a site served by a sewerline.

Michigan↗
USGS2008

Documentation of Computer Program INFIL3.0 - A Distributed-Parameter Watershed Model to Estimate Net Infiltration Below the Root Zone

This report documents the computer program INFIL3.0, which is a grid-based, distributed-parameter, deterministic water-balance watershed model that calculates the temporal and spatial distribution of daily net infiltration of water across the lower boundary of the root zone. The bottom of the root zone is the estimated maximum depth below ground surface affected by evapotranspiration. In many field applications, net infiltration below the bottom of the root zone can be assumed to equal net recharge to an underlying water-table aquifer. The daily water balance simulated by INFIL3.0 includes precipitation as either rain or snow; snowfall accumulation, sublimation, and snowmelt; infiltration into the root zone; evapotranspiration from the root zone; drainage and water-content redistribution within the root-zone profile; surface-water runoff from, and run-on to, adjacent grid cells; and net infiltration across the bottom of the root zone. The water-balance model uses daily climate records of precipitation and air temperature 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 model does not simulate streamflow originating as ground-water discharge. Drainage-basin characteristics are represented in the model by a set of spatially distributed input variables uniquely assigned to each grid cell of a model grid. The report provides a description of the conceptual model of net infiltration on which the INFIL3.0 computer code is based and a detailed discussion of the methods by which INFIL3.0 simulates the net-infiltration process. The report also includes instructions for preparing input files necessary for an INFIL3.0 simulation, a description of the output files that are created as part of an INFIL3.0 simulation, and a sample problem that illustrates application of the code to a field setting. Brief descriptions of the main program routine and of each of the modules and subroutines of the INFIL3.0 code, as well as definitions of the variables used in each subroutine, are provided in an appendix.

Scientific Investigations Report↗
USGS2008

Davis Pond freshwater prediversion biomonitoring study: freshwater fisheries and eagles

In January 2001, the construction of the Davis Pond freshwater diversion structure was completed by the U.S. Army Corps of Engineers. The diversion of freshwater from the Mississippi River is intended to mitigate saltwater intrusion from the Gulf of Mexico and to lessen the concomitant loss of wetland areas. In addition to the freshwater inflow, Barataria Bay basin would receive nutrients, increased flows of sediments, and water-borne and sediment-bound compounds. The purpose of this biomonitoring study was, therefore, to serve as a baseline for prediversion concentrations of selected contaminants in bald eagle ( Haliaeetus leucocephalus ) nestlings (hereafter referred to as eaglets), representative freshwater fish, and bivalves. Samples were collected from January through June 2001. Two similarly designed postdiversion studies, as described in the biological monitoring program, are planned. Active bald eagle nests targeted for sampling eaglet blood (n = 6) were generally located southwest and south of the diversion structure. The designated sites for aquatic animal sampling were at Lake Salvador, at Lake Cataouatche, at Bayou Couba, and along the Mississippi River. Aquatic animals representative of eagle prey were collected. Fish were from three different trophic levels and have varying feeding strategies and life histories. These included herbivorous striped mullet ( Mugil cephalus ), omnivorous blue catfish ( Ictalurus furcatus ), and carnivorous largemouth bass ( Micropterus salmoides ). Three individuals per species were collected at each of the four sampling sites. Freshwater Atlantic rangia clams ( Rangia cuneata ) were collected at the downstream marsh sites, and zebra mussels ( Dreissena spp.) were collected on the Mississippi River. The U.S. Geological Survey (USGS) Biomonitoring of Environmental Status and Trends (BEST) protocols served as guides for fish sampling and health assessments. Fish are useful for monitoring aquatic ecosystems because they accumulate pesticides and other contaminants. Biomarker data on individual fish, generated at the USGS National Wetlands Research Center (Lafayette, La.), included percent white blood cells in whole blood, spleen weight to body weight ratio, liver weight to body weight ratio, condition factor, splenic macrophage aggregates, and liver microsomal 7-ethoxyresorufin-o-deethylase (EROD) activity. Fish age was estimated by comparing total lengths with values from the same species in the Southeast United States as determined from the literature. Contaminant analyses were coordinated by the U.S. Fish and Wildlife Service (USFWS) Analytical Control Facility (Laurel, Md.), where residues of organochlorine (OC) pesticides, total polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons (AHs), and trace elements were determined. The organic contaminant data were generated at the Mississippi State University Chemical Lab (Mississippi State, Miss.), and the inorganic contaminant data were generated by the Texas A&M University Geochemical and Environmental Research Group (College Station, Tex.). Statistical tests were performed to assess relationships among contaminants, fish age, fish species, and collection sites. Trends in interspecific differences among fish in concentrations of contaminants were noted. Striped mullet (hereafter mullet) frequently displayed the highest chemical concentrations. Levels of contaminants were generally higher in samples obtained from the Mississippi River than in those collected from the diversion area and were higher in mussels and clams (hereafter bivalves) than in fish. Because the Mississippi River sampling site for mullet and largemouth bass was downriver of the structure and south of New Orleans and the catfish site was upriver, the downriver data may not be directly reflective of the results from the receiving waters at the Davis Pond structure. Compared to the Caernarvon freshwater prediversion study in 1990 that assessed possible influx of contaminants with the freshwater diversion, contaminant levels in fishes and bivalves in this study were generally lower, yet three nontoxic inorganic elements in Davis Pond fish samples exhibited ranges of concentrations that were more than two times higher than did those from Caernarvon. Levels in bivalves were different between diversions but about equal in the numbers of trace elements showing high levels per location. Contaminant values were compared to those listed in various literature and agency sources, both regional and national, including the National Contaminant Biomonitoring Program (NCBP), in which the 85th percentile and above represents what is considered to be an elevated contaminant concentration and cause for concern. Generally, bivalves were at the high end of their ranges for both organic and inorganic contaminants. In this study, OCs were detectable in 67 percent of fish from the Mississippi River site, ranging from 0.15 to 1.09 μg/g wet weight (ww) or fresh weight (fw), and in 11 percent of the fish from the marsh sites, ranging from 0.06 to 0.612 μg/g ww. Bivalves from the Mississippi River had OC levels of 0.096 μg/g ww, whereas none were detectable in bivalves at the marsh sites. In this study, p,p ’-dichlorodiphenyldichloroethylene ( p,p’ -DDE) (a biodegradation product of DDT [dichlorodiphenyl trichloroethane]) and total PCBs were the most frequently detected OCs and were primarily from the Mississippi River. For total OC content, using adjusted least squares means, some significant interactions were noted between fish species and sites. PAHs were detected in aquatic animals at all sites (range of 0.017–17.534 μg/g ww), as were AHs (range of 0.423–4.549 μg/g ww); the highest levels of PAHs and AHs were found in bivalves from the Mississippi River. When analysis of variance (α = 0.05) was performed with data from aquatic animals, there were only two significant relationships between PAHs, AHs, and OCs between species, site, and age or the interaction among these variables. There was an interaction between fish species and n-decane (an AH) in that mullet and largemouth bass had significantly higher levels than did catfish ( P = 0.0175). When general linear means were used to investigate associations of inorganic contaminants among fish species, site, and age or any interactions among these variables, no significant results were noted for arsenic, cadmium, lead, beryllium, boron, molybdenum, or nickel. The range of mercury in fish in this study was 0.04–0.14 μg/g ww (0.14– 0.48 μg/g dry weight [dw]), with the most elevated levels detected in predatory largemouth bass at the sampling point farthest downstream from the structure and within the marsh area. Mercury was positively correlated with fish age ( P = 0.0152), where levels were estimated to increase 0.0253 parts per million (ppm) dw per year. In the Mississippi River, catfish showed significantly higher levels of mercury than did mullet or largemouth bass ( P = 0.00167). Among fish species, mullet displayed the highest levels in fish of aluminum, barium, manganese, and iron, all considered to have low toxicity in hydrologic systems. An interaction between fish and site was seen with aluminum ( P = 0.0031), where concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with barium ( P = 0.0009), chromium ( P = <0.0001), manganese ( P = 0.0004), strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059). For iron ( P = 0.0.0001), mullet and largemouth bass at both the Mississippi River and Lake Salvador showed higher levels than did catfish, and these two species showed higher levels at two of the four sites. An interaction between fish and site was also seen with chromium ( P = <0.0001) in that concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059), metals for which deleterious effects have been demonstrated in other ecosystems. The NCBP program lists the 85th percentile for zinc at 34.2 μg/g fw (117.9 μg/g dw). In the Davis Pond prediversion biomonitoring study (hereafter the current study), one fish (MUL31RIVER, fish ID 8) showed values higher than that (125.4 μg/g dw or 37.54 μg/g ww), and the Mississippi River bivalve sample (MUSSRIVER) had a value of 140 μg/g dw (41.2 μg/g ww). In the current study, approximately 86 percent of the fish had measurable selenium levels, yet none reached the 85th percentile. The 85th percentile for selenium from the NCBP was 0.73 μg/g ww. Significantly higher levels of selenium were seen in mullet than in largemouth bass and catfish ( P = 0.0023). The NCBP 85th percentile for lead is 0.22 μg/g ww (0.76 μg/g dw). In the current study, the range of concentrations of lead was as much as 18.3 ppm dw (MUL31RIVER, fish ID 8), with the three most elevated values (range of 3.46–5.31 μg/g ww) coming from mullet from the Mississippi River. Biomarker data are measurable and directly reflect the condition of the animal, and measuring more than one biomarker in an individual increases confidence in health assessments. In the current study, biomarkers included macrophage aggregates (MAs), liver (hepatosomatic index [HSI]) and spleen (splenosomatic index [SSI]) weight to body weight ratios, percent white blood cells (WBCs) in whole blood, and condition factor. Few significant differences were noted with any of the biomarkers between sites, and there were no relationships between species and sites. For improved use of biomarker assessments, an increase in fish sample size would be useful for postdiversion sampling, as would comparisons of fish of the same sex and reproductive condition. During the current study, success for eagle nests in the diversion area and reference sites was similar as determined by numbers of nestlings fledged. When temperatures were below average during winter 2000, nests in both regions similarly failed. At each nest, the primary evidence of food items was small mammals. Eaglets (n = 6) generally appeared healthy, and whole blood concentrations of organic contaminants exceeded detection limits with three incidences of p,p’ -DDE (0.002–0.006 μg/L ww) and one incidence of oxychlordane (0.002 μg/L ww). The levels of p,p’ -DDE were well below those that have been inversely correlated with productivity and success rates of nesting bald eagles on a regional scale. The low values found in the whole blood samples for OC pesticides and PCBs were even lower when corrected for plasma volume. Aluminum values were 3.66 and 5.75 μg/L in two samples, zinc ranged from 5.21 to 6.77 μg/L ww in six samples, and silicon ranged from 1.7 to 4.6 μg/L in four samples. Selenium was detectable in each bird with the range at 0.332–0.566 μg/L ww, and strontium ranged from 0.0581 to 0.0975 μg/L ww. Mercury was detectable in blood samples from each bird and ranged from 0.0254 to 0.0845 μg/L ww, whereas lead was detectable in four samples and ranged from 0.0042 to 0.0136 μg/L ww. Although no detectable levels of total PCBs were found (also correlated with decreased reproductive productivity), 70 percent of the aquatic animals from the Mississippi River contained total PCBs (range 0.13–0.79 μg/L), whereas only about 7 percent of the aquatic animals sampled from the marsh area contained PCBs. Suggestions for postdiversion sampling include lowering the analytical detection limit for some metals, sampling aquatic animals over the course of a single season, obtaining a higher sample number of mature fish of one species (for example, blue catfish) within a range of total lengths for biomarker analyses, obtaining otoliths for estimating fish ages, assessing dioxins in eaglet blood, examining triazines in water, and obtaining all Mississippi River fish samples as close to the Davis Pond structure intake as possible. Because contaminants found in blood of eaglets reflect their prey species and because of the contaminant levels found in fish in the current study, eaglets may not be consuming primarily these species; therefore, obtaining juvenile nutria ( Myocastor coypus ) or turtle species for contaminant analyses might be considered, as well as collecting greater blood volume and using plasma to measure OCs and PCBs. Data obtained postdiversion will be compared with prediversion data to monitor changes.

Louisiana↗
USGS2009

The Mississippi Embayment Regional Aquifer Study (MERAS): Documentation of a groundwater-flow model constructed to assess water availability in the Mississippi embayment

The Mississippi Embayment Regional Aquifer Study (MERAS) was conducted with support from the Groundwater Resources Program of the U.S. Geological Survey Office of Groundwater. This report documents the construction and calibration of a finite-difference groundwater model for use as a tool to quantify groundwater availability within the Mississippi embayment. To approximate the differential equation, the MERAS model was constructed with the U.S. Geological Survey's modular three-dimensional finite-difference code, MODFLOW-2005; the preconditioned conjugate gradient solver within MODFLOW-2005 was used for the numerical solution technique. The model area boundary is approximately 78,000 square miles and includes eight States with approximately 6,900 miles of simulated streams, 70,000 well locations, and 10 primary hydrogeologic units. The finite-difference grid consists of 414 rows, 397 columns, and 13 layers. Each model cell is 1 square mile with varying thickness by cell and by layer. The simulation period extends from January 1, 1870, to April 1, 2007, for a total of 137 years and 69 stress periods. The first stress period is simulated as steady state to represent predevelopment conditions. Areal recharge is applied throughout the MERAS model area using the MODFLOW-2005 Recharge Package. Irrigation, municipal, and industrial wells are simulated using the Multi-Node Well Package. There are 43 streams simulated by the MERAS model. Each stream or river in the model area was simulated using the Streamflow-Routing Package. The perimeter of the model area and the base of the flow system are represented as no-flow boundaries. The downgradient limit of each model layer is a no-flow boundary, which approximates the extent of water with less than 10,000 milligrams per liter of dissolved solids. The MERAS model was calibrated by making manual changes to parameter values and examining residuals for hydraulic heads and streamflow. Additional calibration was achieved through alternate use of UCODE-2005 and PEST. Simulated heads were compared to 55,786 hydraulic-head measurements from 3,245 wells in the MERAS model area. Values of root mean square error between simulated and observed hydraulic heads of all observations ranged from 8.33 feet in 1919 to 47.65 feet in 1951, though only six root mean square error values are greater than 40 feet for the entire simulation period. Simulated streamflow generally is lower than measured streamflow for streams with streamflow less than 1,000 cubic feet per second, and greater than measured streamflow for streams with streamflow more than 1,000 cubic feet per second. Simulated streamflow is underpredicted for 18 observations and overpredicted for 10 observations in the model. These differences in streamflow illustrate the large uncertainty in model inputs such as predevelopment recharge, overland flow, pumpage (from stream and aquifer), precipitation, and observation weights. The groundwater-flow budget indicates changes in flow into (inflows) and out of (outflows) the model area during the pregroundwater-irrigation period (pre-1870) to 2007. Total flow (sum of inflows or outflows) through the model ranged from about 600 million gallons per day prior to development to 18,197 million gallons per day near the end of the simulation. The pumpage from wells represents the largest outflow components with a net rate of 18,197 million gallons per day near the end of the model simulation in 2006. Groundwater outflows are offset primarily by inflow from aquifer storage and recharge.

Alabama, Arkansas, Kentucky, Louisiana, Mississipp↗
USGS2011

Numerical simulation of the groundwater-flow system in the Chambers-Clover Creek Watershed and Vicinity, Pierce County, Washington

A groundwater-flow model was developed to contribute to an improved understanding of water resources in the Chambers-Clover Creek Watershed. The model covers an area of about 491 square miles in western Pierce County, Washington, and is bounded to the northeast by the Puyallup River valley, to the southwest by the Nisqually River valley, and extends northwest to Puget Sound, and southeast to Tanwax Creek. The Puyallup and Nisqually Rivers occupy large, relatively flat alluvial valleys that are separated by a broad, poorly drained, upland area that covers most of the model area. Chambers and Clover Creeks drain much of the central uplands and flow westward to Puget Sound. The model area is underlain by a northwest-thickening sequence of unconsolidated glacial (till and outwash) and interglacial (fluvial and lacustrine) deposits. Ten unconsolidated hydrogeologic units in the model area form the basis of the groundwater-flow model. Groundwater flow in the Chambers-Clover Creek Watershed and vicinity was simulated using the groundwater-flow model, MODFLOW-2000. The finite-difference model grid comprises 146 rows, 132 columns, and 11 layers. Each model cell has a horizontal dimension of 1,000 by 1,000 feet, and the model contains a total of 123,602 active cells. The thickness of model layers varies throughout the model area and ranges from 1.5 feet in the A3 aquifer unit to 1,567 feet in the G undifferentiated unit. Groundwater flow was simulated for both steady-state and transient conditions. Steady-state conditions were simulated using average recharge, discharge, and water levels for the 24-month period September 2006-August 2008. Transient conditions were simulated for the period September 2006-August 2008 using 24 monthly stress periods. Resource managers and local stakeholders intend to use the model to evaluate a range of water resource issues under both steady-state and transient conditions. Initial conditions for the transient model were developed from a 3-year "lead-in" period that used recorded precipitation and river levels, and temporal extrapolations of other boundary conditions. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. The model as calibrated to steady-state conditions has a standard deviation for heads and flows of 28.42 feet and 2.12 cubic feet per second, respectively; the model as calibrated to transient conditions has a standard deviation for heads and flows of 23.01 feet and 2.67 cubic feet per second, respectively. Simulated steady-state inflow to the model area from precipitation and secondary recharge was 477,266 acre-feet per year (acre-ft/yr) (79 percent of total simulated inflow), and simulated inflow from stream and lake leakage was 129,778 acre-ft/yr (21 percent of total simulated inflow). Simulated outflow from the model primarily was through discharge to streams, lakes, springs, seeps, and Puget Sound (559,192 acre-ft/yr; 92 percent of total simulated outflow), and withdrawals from wells (47,863 acre-ft/yr; 8 percent of total simulated outflow). Six scenarios were formulated and simulated using the calibrated model to provide representative examples of how the model can be used to evaluate the effects on groundwater levels and stream baseflows of potential changes in groundwater withdrawals, in consumptive use, and in recharge.

Washington↗