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 397 records · Page 22Linked to original sources
USGS2008

Assessment of the Mowry Shale and Niobrara Formation as Continuous Hydrocarbon Systems, Powder River Basin, Montana and Wyoming

A recent U.S. Geological Survey (USGS) oil and gas assessment of the Powder River Basin , Wyoming and Montana, identified the Upper Cretaceous Mowry Shale and Niobrara Formation as the primary hydrocarbon sources for Cretaceous conventional and unconventional reservoirs. Cumulative Mowry-sourced petroleum production is about 1.2 BBO (billion barrels of oil) and 2.2 TCFG (trillion cubic feet of gas) and cumulative Niobrara-sourced oil production is about 520 MMBO (million barrels of oil) and 0.95 TCFG. Burial history modeling indicated that hydrocarbon generation for both formations started at about 0.60 percent Ro at depths of about 8,000 ft. At maximum depths, Ro for the Mowry is about 1.2 to 1.3 percent and about 0.80 percent for the Niobrara. The Mowry and Niobrara continuous reservoirs were assessed using a cell-based methodology that utilized production data. The size of each cell was based on geologic controls and potential drainage areas in analog fields. Current and historical production data were used to determine the estimated ultimate recovery (EUR) distribution for untested cells. Only production data from unconventional fractured shale reservoirs with vertical wells were used. For the Mowry, the minimum, median, and maximum total recovery volumes per cell for untested cells are (1) 0.002, 0.25, and 0.35 MMBO, respectively; and for the Niobrara (2) 0.002, 0.028, and 0.5 MMBO. Sweet spots were identified by lineaments and faults, which are believed to be areas having the greatest petroleum potential; an upper limit of 8,000 ft depth was defined by overpressuring caused by hydrocarbon generation. Mean estimates of technically recoverable undiscovered continuous resource for the Mowry are 198 MMBO, 198 BCF (billion cubic feet of gas), and 11.9 MMBNGL (million barrels of natural gas liquid), and those for the Niobrara are 227 MMBO, 227 BCFG, and 13.6 MMBNGL.

Open-File Report↗
USGS2015

GRIDGEN Version 1.0: a computer program for generating unstructured finite-volume grids

GRIDGEN is a computer program for creating layered quadtree grids for use with numerical models, such as the MODFLOW–USG program for simulation of groundwater flow. The program begins by reading a three-dimensional base grid, which can have variable row and column widths and spatially variable cell top and bottom elevations. From this base grid, GRIDGEN will continuously divide into four any cell intersecting user-provided refinement features (points, lines, and polygons) until the desired level of refinement is reached. GRIDGEN will then smooth, or balance, the grid so that no two adjacent cells, including overlying and underlying cells, differ by more than a user-specified level tolerance. Once these gridding processes are completed, GRIDGEN saves a tree structure file so that the layered quadtree grid can be quickly reconstructed as needed. Once a tree structure file has been created, GRIDGEN can then be used to (1) export the layered quadtree grid as a shapefile, (2) export grid connectivity and cell information as ASCII text files for use with MODFLOW–USG or other numerical models, and (3) intersect the grid with shapefiles of points, lines, or polygons, and save intersection output as ASCII text files and shapefiles. The GRIDGEN program is demonstrated by creating a layered quadtree grid for the Biscayne aquifer in Miami-Dade County, Florida, using hydrologic features to control where refinement is added.

Open-File Report↗
USGS2018

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and Northern Chihuahua, Mexico

Errata **September 28, 2018: The purpose of a USGS Open-file report (OFR) is dissemination of information that must be released immediately to fill a public need or information that is not sufficiently refined to warrant publication in one of the other USGS series. As part of that refinement process, an error was discovered in one of the input data sets of the Rio Grande Transboundary Integrated Hydrologic Model (RGTIHM) that this OFR was based upon. The error involved the assignment of storage properties to “phantom cells.” Phantom cells are required for most variants of MODFLOW that use a structured finite-difference grid when individual stratigraphic layers are represented as separate layers. Using phantom cells is a common practice that allows separate model layers to be maintained without having to combine stratigraphic layers into equivalent model layers or to use an unstructured grid. Typically, phantom cell horizontal hydraulic conductivities and storage properties are set to a small number and vertical hydraulic conductivities are set to a number large enough to allow vertical flow between the vertically adjacent layers. In the RGTIHM, the specific storage properties of the phantom cells for the upper (RGTIHM layers 3 and 4), middle (RGTIHM layers 5 and 6), and lower (RGTIHM layers 7 and 8) members of the Santa Fe Group were inadvertently assigned a value of 1 feet -1 . The revision of these specific storage values to a small number (1.0 x 10 -09 feet -1 ) required additional trial-and-error model calibration and a new sensitivity analysis. After calibration, the overall model fit remained similar to the fit described in the OFR, but the fit for many individual features such as project water available for diversions at the American Canal and Acequia Madre improved due to the reduction in flow coming from lower layers. Overall, there is still an average net depletion of groundwater flow, and the conclusions of the report are not changed. The revised average annual groundwater flow depletion simulated for the period 1953-2014 is -1,480 acre-feet/year for the entire model region, and -3,660 acre-feet/year for the portion of the model in the United States. The final version of the model will be the basis of the USGS Scientific Investigations Report that will supersede this OFR. An updated Model Archive of RGTIHM is available upon request to the USGS California Water Science Center. The corrected version of the model WAS the basis for the USGS Scientific Investigations Report that SUPERSEDED this Open-File Report.** Abstract Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model (BCM), and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model (referred to as One Water). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 8,990 acre-feet per year.

New Mexico, Northern Chihuahua, Texas↗
USGS1976

Results of phytoplankton sampling at National Stream Quality Accounting Network stations in Montana, 1975 water year

Twelve National Stream Quality Accounting Network stations were operated in Montana during the 1975 water year. The network was established for the purpose of acquiring a base of hydrologic data for use by agencies engaged in water-resources planning on a national or regional scale. Among the characteristics analyzed were phytoplankton identification and cell counts. Samples consisted of composites of equal aliquots, collected at the center of each quartile of flow, using modified suspended-sediment samplers and sediment collection techniques. Identification and counting were done using the Sedgwick-Rafter cell method. Cell counts ranged from 21 cells per millilitre at Flathead River at Flathead, British Columbia to 27,000 cells per millilitre at Yellowstone River near Sidney. The class Bacillariophyceae was most abundant in both number and variety at all sampling sites. Anabaena and Aphanizomenon of the phyllum Cyanophyta were found at six stations and two stations respectively. These two genera of blue-green algae often become abundant in enriched waters resulting in nuisance conditions.

Montana↗
USGS2006

Effects of Proposed Additional Ground-Water Withdrawals from the Mississippi River Valley Alluvial Aquifer on Water Levels in Lonoke County, Arkansas

The Grand Prairie Water Users Association, located in Lonoke County, Arkansas, plans to increase ground-water withdrawals from the Mississippi River Valley alluvial aquifer from their current (2005) rate of about 400 gallons per minute to 1,400 gallons per minute (2,016,000 gallons per day). The effect of pumping from a proposed well was simulated using a digital model of ground-water flow. The proposed additional withdrawals were added to an existing pumping cell specified in the model, with increased pumping beginning in 2005, and specified to pump at a total combined rate of 2,016,000 gallons per day for a period of 46 years. In addition, pumping from wells owned by Cabot Water Works, located about 2 miles from the proposed pumping, was added to the model beginning in 2001 and continuing through to the end of 2049. Simulated pumping causes a cone of depression to occur in the alluvial aquifer with a maximum decline in water level of about 8.5 feet in 46 years in the model cell of the proposed well compared to 1998 withdrawals. However, three new dry model cells occur south of the withdrawal well after 46 years. This total water-level decline takes into account the cumulative effect of all wells pumping in the vicinity, although the specified pumping rate from all other model cells in 2005 is less than for actual conditions in 2005. After 46 years with the additional pumping, the water-level altitude in the pumped model cell was about 177.4 feet, which is 41.7 feet higher than 135.7 feet, the altitude corresponding to half of the original saturated thickness of the alluvial aquifer (a metric used to determine if the aquifer should be designated as a Critical Ground-Water Area (Arkansas Natural Resources Commission, 2006)).

Scientific Investigations Report↗
USGS2007

LiDAR-Derived Flood-Inundation Maps for Real-Time Flood-Mapping Applications, Tar River Basin, North Carolina

Flood-inundation maps were created for selected streamgage sites in the North Carolina Tar River basin. Light detection and ranging (LiDAR) data with a vertical accuracy of about 20 centimeters, provided by the Floodplain Mapping Information System of the North Carolina Floodplain Mapping Program, were processed to produce topographic data for the inundation maps. Bare-earth mass point LiDAR data were reprocessed into a digital elevation model with regularly spaced 1.5-meter by 1.5-meter cells. A tool was developed as part of this project to connect flow paths, or streams, that were inappropriately disconnected in the digital elevation model by such features as a bridge or road crossing. The Hydraulic Engineering Center-River Analysis System (HEC-RAS) model, developed by the U.S. Army Corps of Engineers, was used for hydraulic modeling at each of the study sites. Eleven individual hydraulic models were developed for the Tar River basin sites. Seven models were developed for reaches with a single gage, and four models were developed for reaches of the Tar River main stem that receive flow from major gaged tributaries, or reaches in which multiple gages were near one another. Combined, the Tar River hydraulic models included 272 kilometers of streams in the basin, including about 162 kilometers on the Tar River main stem. The hydraulic models were calibrated to the most current stage-discharge relations at 11 long-term streamgages where rating curves were available. Medium- to high-flow discharge measurements were made at some of the sites without rating curves, and high-water marks from Hurricanes Fran and Floyd were available for high-stage calibration. Simulated rating curves matched measured curves over the full range of flows. Differences between measured and simulated water levels for a specified flow were no more than 0.44 meter and typically were less. The calibrated models were used to generate a set of water-surface profiles for each of the 11 modeled reaches at 0.305-meter increments for water levels ranging from bankfull to approximately the highest recorded water level at the downstream-most gage in each modeled reach. Inundated areas were identified by subtracting the water-surface elevation in each 1.5-meter by 1.5-meter grid cell from the land-surface elevation in the cell through an automated routine that was developed to identify all inundated cells hydraulically connected to the cell at the downstream-most gage in the model domain. Inundation maps showing transportation networks and orthoimagery were prepared for display on the Internet. These maps also are linked to the U.S. Geological Survey North Carolina Water Science Center real-time streamflow website. Hence, a user can determine the near real-time stage and water-surface elevation at a U.S. Geological Survey streamgage site in the Tar River basin and link directly to the flood-inundation maps for a depiction of the estimated inundated area at the current water level. Although the flood-inundation maps represent distinct boundaries of inundated areas, some uncertainties are associated with these maps. These are uncertainties in the topographic data for the hydraulic model computational grid and inundation maps, effective friction values (Manning's n), model-validation data, and forecast hydrographs, if used. The Tar River flood-inundation maps were developed by using a steady-flow hydraulic model. This assumption clearly has less of an effect on inundation maps produced for low flows than for high flows when it typically takes more time to inundate areas. A flood in which water levels peak and fall slowly most likely will result in more inundation than a similar flood in which water levels peak and fall quickly. Limitations associated with the steady-flow assumption for hydraulic modeling vary from site to site. The one-dimensional modeling approach used in this study resulted in good agreement between measurements and simulations. T

North Carolina↗
USGS2007

Marine geophysical investigation of selected sites in Bridgeport Harbor, Connecticut, 2006

A marine geophysical investigation was conducted in 2006 to help characterize the bottom and subbottom materials and extent of bedrock in selected areas of Bridgeport Harbor, Connecticut. The data will be used by the U.S. Army Corps of Engineers in the design of confined aquatic disposal (CAD) cells within the harbor to facilitate dredging of the harbor. Three water-based geophysical methods were used to evaluate the geometry and composition of subsurface materials: (1) continuous seismic profiling (CSP) methods provide the depth to water bottom, and when sufficient signal penetration can be achieved, delineate the depth to bedrock and subbottom materials; (2) continuous resistivity profiling (CRP) methods were used to define the electrical properties of the shallow subbottom, and to possibly determine the distribution of conductive materials, such as clay, and resistive materials, such as sand and bedrock; (3) and magnetometer data were used to identify conductive anomalies of anthropogenic sources, such as cables and metallic debris. All data points were located using global positioning systems (GPS), and the GPS data were used for real-time navigation. The results of the CRP, CSP, and magnetometer data are consistent with the conceptual site model of a bedrock channel incised beneath the present day harbor. The channel appears to follow a north-northwest to south-southeast trend and is parallel to the Pequannock River. The seismic record and boring data indicate that under the channel, the depth to bedrock is as much as 42.7 meters (m) below mean low-low water (MLLW) in the dredged part of the harbor. The bedrock channel becomes shallower towards the shore, where bedrock outcrops have been mapped at land surface. CSP and CRP data were able to provide a discontinuous, but reasonable, trace from the channel toward the west under the proposed southwestern CAD cell. The data indicate a high amount of relief on the bedrock surface, as well as along the water bottom. Under the southwestern CAD cell, the sediments are only marginally thick enough for a CAD cell, at about 8 to 15 m in depth. Some of the profiles show small diffractions in the unconsolidated sediments, but no large-scale boulders or boulder fields were identified. No bedrock reflectors were imaged under the southeastern CAD cell, where core logs indicate the rock is as much as 30 m below MLLW. The chirp frequency, tuned transducer, and boomer-plate CSP surveys were adversely affected by a highly reflective water bottom causing strong multiples in the seismic record and very limited depths of penetration. These multiples are attributed to entrapped gas (methane) in the sediments or to very hard bottom conditions. In a limited number of places, the bedrock surface was observed in the CSP record, creating a discontinuous and sporadic image of the bedrock surface. These interpretations generally matched core data at FP-03-10 and FB-06-1. Use of two analog CSP systems, the boomer plate and tuned transducer, did not overcome the reflections off the water bottom and did not improve the depth of penetration. In general, the CRP profiles were used to corroborate the results of the CSP profiles. Relatively resistive zones associated with the locations of seismic reflections were interpreted as bedrock. The shape of the bedrock surface generally was similar in the CRP and CSP profiles. Evaluation of the CRP profiles indicated that the inversions were adversely affected where the depth and (or) ionic concentration of the water column varied. Consequently, the CRP profiles were broken into short intervals that extended just over the area of interest, where the depth to water bottom was fairly constant. Over these short profiles, efforts were made to evaluate the resistivity of the very shallow sediments to determine if there were any large contrasts in the resistivity of the sediments that might indicate differences in the shallow subbottom materials.

Connecticut↗
USGS2015

Groundwater/surface-water interactions in the Bad River Watershed, Wisconsin

A groundwater-flow model was developed for the Bad River Watershed and surrounding area by using the U.S. Geological Survey (USGS) finite-difference code MODFLOW-NWT. The model simulates steady-state groundwater-flow and base flow in streams by using the streamflow routing (SFR) package. The objectives of this study were to: (1) develop an improved understanding of the groundwater-flow system in the Bad River Watershed at the regional scale, including the sources of water to the Bad River Band of Lake Superior Chippewa Reservation (Reservation) and groundwater/surface-water interactions; (2) provide a quantitative platform for evaluating future impacts to the watershed, which can be used as a starting point for more detailed investigations at the local scale; and (3) identify areas where more data are needed. This report describes the construction and calibration of the groundwater-flow model that was subsequently used for analyzing potential locations for the collection of additional field data, including new observations of water-table elevation for refining the conceptualization and corresponding numerical model of the hydrogeologic system. The study area can be conceptually divided into three primary hydrogeologic environments. The first encompasses the southern uplands with relatively low topographic relief, where groundwater-flow is unconfined and occurs primarily in sandy till and glacial outwash overlying Archean-aged crystalline bedrock. The second includes a transitional area of higher topographic relief and shallow depth to bedrock, in the vicinity of ridges formed by steeply dipping, early-Proterozoic aged metasedimentary units of the Marquette Range Supergroup (including the Ironwood Formation), and late-Proterozoic igneous units associated with the Midcontinent Rift System (MRS). Groundwater-flow in this area likely occurs primarily through connected networks of bedrock fractures that are not well characterized, and also in isolated pockets of Quaternary deposits. The third and last hydrogeologic environment includes lowlands along Lake Superior where a deep sandstone aquifer is confined by thick deposits of clay-rich till. Model input was compiled by using both published and unpublished data. Constant flux boundary conditions for the model perimeter were developed from a regional analytic element model described in appendix 1 of this report. Pumping from 26 high-capacity wells within the model area was included. The SFR stream network was developed from the National Hydrography Dataset (NHDPlus Version 2) and hydrography from the Wisconsin Department of Natural Resources (WDNR). Hydraulic conductivity values were determined for each model cell by interpolation from a network of pilot points, within zones representing major hydrogeologic units. Recharge to the groundwater system was estimated on a cell-by-cell basis by using the Soil Water Balance code (SWB), with gridded daily temperature and precipitation data for the period 1980–2011, and GIS coverages of soil and land-surface conditions. Estimated recharge varies considerably, following spatial patterns in the precipitation and soil hydrologic group inputs. The lowest recharge values occur in the Superior lowlands, whereas the highest values occur in the upland areas, especially those underlain by sandy soils, and in the vicinity of bedrock hills. The model was calibrated to groundwater-levels and base flows obtained from the USGS National Water Information System (NWIS) database, and groundwater-levels obtained from the WDNR and Band River Band well-construction databases. Calibration was performed via nonlinear regression by using the parameter-estimation software suite PEST. Groundwater levels and base-flow observations in the calibration dataset were well simulated by the calibrated model, with reasonable values of hydraulic conductivity. The pilot-point parameters that were most constrained by observations during model calibration coincided with the locations containing the most wells (head observations)—especially the population centers of Ashland, Mellen, and other communities along the major highway corridors. Results from the calibrated model illustrate differences in the nature of groundwater-flow within the watershed. In the southern part of the watershed, where bedrock is shallow, groundwater flow paths are relatively short, extending from local recharge areas to adjacent first and second-order streams. In contrast, laterally continuous deposits of clay-rich till covering the Superior Lowlands isolate most smaller streams from the sandstone aquifer, allowing for longer flow paths toward larger streams such as the Bad, Marengo, and White Rivers. Approximately three-quarters of all first-order stream cells were dry in the Superior Lowlands, compared to only half of first-order stream cells in the southern bedrock uplands. The model was used to delineate the groundwatershed for the Bad and Kakagon Rivers. “Groundwatershed” is defined as the area contributing groundwater discharge to one of these streams and their tributaries. The groundwatershed was found to align closely with the surface-watershed, with the most notable exception occurring along the southwestern half of Birch Hill, where surface water drains southwest towards the Potato River, and groundwater flows north and east towards Lake Superior. Similarly, the contributing area of groundwater-flow to the Reservation was delineated. Results indicate the off-Reservation groundwater contributing area to be limited in comparison to the extent of the watershed, extending southward into the highlands underlain by MRS igneous rock units, but not further into the area underlain by the Marquette Range Supergroup. Stable isotope samples were collected from 54 wells within the watershed, to investigate sources of groundwater. Oxygen-18 (δ 18 O) values lower than -13.0 per mil were documented in the sampling, and likely indicate the presence of recharge water from the last glacial period (>9,500 years old) beneath the northern portion of the Reservation, in the vicinity of Odanah, Wisconsin. Finally, a new data-worth analysis of potential new monitoring-well locations was performed by using the model. The relative worth of new measurements was evaluated based on their ability to increase confidence in model predictions of groundwater levels and base flows at 35 locations, under the condition of a proposed open-pit iron mine. Results of the new data-worth analysis, and other inputs and outputs from the Bad River model, are available through an online dynamic web mapping service at ( http://wim.usgs.gov/badriver/ ).

Wisconsin↗
USGS2024

Evidence of nitrate attenuation in intertidal and subtidal groundwater in a subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2015–16

Nitrogen dynamics in intertidal and nearshore subtidal groundwater (subterranean estuary) adjacent to the Seacoast Shores peninsula, Falmouth, Massachusetts, were investigated during 2015–16 by the U.S. Geological Survey. The peninsula is a densely populated residential area with septic systems and cesspools that are substantial sources of nitrogen to groundwater. The study area is in the Eel River, an estuarine saltwater embayment connected to the ocean adjacent to the western shore of the peninsula, that was the subject of an earlier study by Colman and others (2018, https://doi.org/10.3133/sir20185095 ) on nitrogen transport and transformations in groundwater between onshore and offshore locations. The previous study documented the distribution of nitrate concentrations and nitrate attenuation reactions in fresh groundwater beneath the peninsula and the estuary. The current study extended those observations with more detailed sampling and analysis of shallow groundwater from wells near discharge sites beneath the estuary. The current field investigation included sampling of existing wells and installation and sampling of clusters of wells and temporary sampling points in the subterranean estuary, including (1) shallow transects 0.3 to 1.2 meters (m) deep extending from 1 to 13.5 m offshore and (2) deeper wells (from 1.83 to 4.88 m deep) extending from 4.3 to 14.3 m offshore. Measurements of hydraulic-head gradients 2–5 m below the sediment/water interface in the intertidal and nearshore subtidal zones indicated that groundwater flow generally was upwards (towards the estuary) under all tide conditions in October 2016. The magnitude of the gradient was greatest during low tide conditions, indicating that groundwater discharge likely decreased during high tides. Measurements of specific conductance in shallow groundwater in the subterranean estuary in three transects perpendicular to shore were consistent with the existence of saltwater flow cells (infiltration of overlying saline water, mixing with fresh groundwater, and discharge to the overlying saline water) in the intertidal and nearshore subtidal regions. The size of these flow cells was variable in space and time and dependent on the elevation of the tide (spring or neap). At this location in the Eel River subterranean estuary, and offshore to at least 13.5 m, offshore flow of fresh groundwater apparently prevented a deeper saltwater wedge from discharging to the surface. Nitrate concentrations in shallow groundwater (30 to 122 centimeters [cm] depth) were variable in space and time, ranging from not detectable to 600 micromoles per liter (μmol/L) (8.4 milligrams per liter as N), and were highest in June 2016 at depths from 61 to 122 cm below the sediment/water interface and from 4 to 9 m offshore. Nitrate generally was not detectable in saline shallow groundwater at 30-cm depth or at any depth from 30 to 122 cm from 10 to 13.5 m offshore. Dissolved oxygen concentrations were suboxic (less than 16 μmol/L) in 60 percent of the sampled subterranean groundwater beneath the intertidal and subtidal zones. In the remaining sites, the range of dissolved oxygen concentrations was from 18 to 272 μmol/L and the median concentration was 43 μmol/L. Evidence for microbial nitrate reduction (denitrification and possibly anammox) was provided by the distribution of the reaction product nitrogen gas (excess N 2 , or N 2MIC ), as determined from analysis of the dissolved nitrogen gas and argon gas (Ar) concentrations in groundwater samples. Excess nitrogen gas provided evidence for nitrate reduction in shallow groundwater below the subtidal and, to a lesser extent, intertidal zones adjacent to the Seacoast Shores peninsula. These zones, where evidence for nitrate reduction was detected, were in fresh and brackish groundwater near subtidal or intertidal saltwater cells where discharging fresh groundwater mixed with infiltrating saline water. Infiltrating seawater may have supplied organic carbon, one of several potential electron donors that are required for denitrification. Other potential electron donors, such as organic carbon, iron, manganese, hydrogen, methane, ammonium, elemental sulfur, or sulfide phases, may have been supplied by the estuarine sediments. Drainage from surface runoff near the shore also may have supplied organic carbon to fresh groundwater near the intertidal saltwater cell. The highest amounts of nitrate converted to excess nitrogen gas were estimated to be in the range of 230 to 430 μmol/L in nearly fresh groundwater near the subtidal saltwater cell at depths of 61 to 122 cm below the sediment/water interface and from 10 to 13.5 m offshore. Evidence of denitrification within 10 m of the shore was sparse (generally limited to less than 50 μmol/L of N 2 -N) despite the presence of high nitrate concentrations. The spatial distribution of estimated nitrate reduction in the intertidal and nearshore subtidal fresh and brackish groundwater may be related to local variability in the distribution of reactive electron donors in those zones. Variations in the amount of nitrate reduction to nitrogen gas were not clearly related to potential aqueous electron donors such as dissolved organic carbon, nor to potential reaction products such as alkalinity, but may have been controlled by combinations of aqueous and solid-phase reactants. The distribution of relatively shallow fresh groundwater containing nitrate could indicate potential nitrate discharge areas in the lower intertidal zone and uncertain locations farther offshore; however, the data did not extend all the way to the sediment/water interface or to the offshore freshwater limit. This study confirmed substantial loss of nitrate from some of the fresh and brackish groundwater in shallow subestuarine sediments prior to discharge but did not quantify how much nitrate eventually discharged to the estuary.

Massachusetts↗
USGS2012

RIP-ET: A riparian evapotranspiration package for MODFLOW-2005

A new evapotranspiration package for the U.S. Geological Survey's groundwater-flow model, MODFLOW, is documented. The Riparian Evapotranspiration Package (RIP-ET) provides flexibility in simulating riparian and wetland transpiration not provided by the Evapotranspiration (EVT) or Segmented Function Evapotranspiration (ETS1) Packages for MODFLOW 2005. This report describes how the RIP-ET package was conceptualized and provides input instructions, listings and explanations of the source code, and an example. Traditional approaches to modeling evapotranspiration (ET) processes assume a piecewise linear relationship between ET flux and hydraulic head. The RIP-ET replaces this traditional relationship with a segmented, nonlinear dimensionless curve that reflects the eco-physiology of riparian and wetland ecosystems. Evapotranspiration losses from these ecosystems are dependent not only on hydraulic head, but on the plant types present. User-defined plant functional groups (PFGs) are used to elucidate the interaction between plant transpiration and groundwater conditions. Five generalized plant functional groups based on transpiration rates, plant rooting depth, and water tolerance ranges are presented: obligate wetland, shallow-rooted riparian, deep-rooted riparian, transitional riparian and bare ground/open water. Plant functional groups can be further divided into subgroups (PFSGs) based on plant size, density or other characteristics. The RIP-ET allows for partial habitat coverage and mixtures of plant functional subgroups to be present in a single model cell. RIP-ET also distinguishes between plant transpiration and bare-ground evaporation. Habitat areas are designated by polygons; each polygon can contain a mixture of PFSGs and bare ground, and is assigned a surface elevation. This process requires a determination of fractional coverage for each of the plant functional subgroups present in a polygon to account for the mixture of coverage types and resulting transpiration. The fractional cover within a cell has two components: (1) the polygonal fraction of active habitat (excluding area of bare ground, dead trees, or brush) in a cell, and (2) fraction of plant type area or bare ground area in a polygon. RIP-ET determines the transpiration rate for each plant functional group and evaporation from bare ground/open water in a cell, the total ET in the cell, and the total ET rate over the region of simulation.

Techniques and Methods↗
USGS1997

Review of selected features of the natural system model, and suggestions for applications in south Florida

A study was conducted to review selected features of the Natural System Model, version 4.3. The Natural System Model is a regional-scale model that uses recent climatic data and estimates of historic vegetation and topography to simulate pre-canal-drainage hydrologic response in south Florida. Equations used to represent the hydrologic system and the numerical solution of these equations in the model were documented and reviewed. Convergence testing was performed using 1965 input data, and selected other aspects of the model were evaluated. Some conclusions from the evaluation of the Natural System Model include the following observations. Simulations were generally insensitive to the temporal resolution used in the model. However, reduction of the computational cell size from 2-mile by 2-mile to 2/3-mile by 2/3-mile resulted in a decrease in spatial mean ponding depths for October of 0.35 foot for a 3-hour time step. Review of the computer code indicated that there is no limit on the amount of water that can be transferred from the river system to the overland flow system, on the amount of seepage from the river to the ground-water system, on evaporation from the river system, or on evapotranspiration from the overland-flow system. Oscillations of 0.2 foot or less in simulated river stage were identified and attributed to a volume limiting function which is applied in solution of the overland-flow equations. The computation of the resistance coefficient is not consistent with the computation of overland-flow velocity. Ground-water boundary conditions do not always ensure a no-flow condition at the boundary. These inconsistencies had varying degrees of effects on model simulations, and it is likely that simulations longer than 1 year are needed to fully identify effects. However, inconsistencies in model formulations should not be ignored, even if the effects of such errors on model results appear to be small or have not been clearly defined. The Natural System Model can be a very useful tool for estimating pre-drainage hydrologic response in south Florida. The model includes all of the important physical processes needed to simulate a water balance. With a few exceptions, these hydrologic processes are represented in a reasonable manner using empirical, semiempirical, and mechanistic relations. The data sets that have been assembled to represent physical features, and hydrologic and meteorological conditions are quite extensive in their scope. Some suggestions for model application were made. Simulation results from the Natural System Model need to be interpreted on a regional basis, rather than cell by cell. The available evidence suggests that simulated water levels should be interpreted with about a plus or minus 1 foot uncertainty. It is probably not appropriate to use the Natural System Model to estimate pre-drainage discharges (as opposed to hydroperiods and water levels) at a particular location or across a set of adjacent computational cells. All simulated results for computational cells within about 10 miles of the model boundaries have a higher degree of uncertainty than results for the interior of the model domain. It is most appropriate to interpret the Natural System Model simulation results in connection with other available information. Stronger linkages between hydrologic inputs to the Everglades and the ecological response of the system would enhance restoration efforts.

Florida↗
USGS1998

Purification of Piscirickettsia salmonis and partial characterisation of antigens

Piscirickettsia salmonis is the etiological agent of salmonid rickettsial septicemia, an economically significant disease affecting the salmon aquaculture industry. As with other rickettsial pathogens, antigenic analysis of P. salmonis has been limited by the inherent difficulties of purifying an intracellular organism away from host cell material. In this report, we describe the use of diatrizoate meglumine and diatrizoate sodium (DMDS) density gradient centrifugation to purify P. salmonis grown in chinook salmon embryo (CHSE-214) cells. Plaque assay titers and total protein assays confirmed that viable P. salmonis was consistently concentrated in a visible band within the DMDS density gradient at a density of 1.15 to 1.16 g ml -1 . Recovery of purified, viable organisms from DMDS density gradients varied from 0.6 to 3%. Preparations of uninfected CHSE-214 cells, CHSE-214 cells infected with P. salmonis , and gradient-purified P. salmonis were compared using sodium dodecyl sulfate polyacrylamide gel electrophoresis to assess the degree of purification and to identify P. salmonis -specific proteins. Although gradient-purified P. salmonis preparations were not completely free of host cell material, 8 bacterial proteins were identified. Polyclonal rabbit antiserum was used in an immunoblot of proteins from purified P. salmonis to identify 3 major and 5 minor antigens. The major antigens of 56, 30 and 20 kDa were potential candidates for experimental vaccines and development of novel diagnostic assays.

Diseases of Aquatic Organisms↗
USGS2005

A cohabitation challenge to compare the efficacies of vaccines for bacterial kidney disease (BKD) in chinook salmon Oncorhynchus tshawytscha

The relative efficacies of 1 commercial and 5 experimental vaccines for bacterial kidney disease (BKD) were compared through a cohabitation waterborne challenge. Groups of juvenile chinook salmon Oncorhynchus tshawytscha were vaccinated with one of the following: (1) killed Renibacterium salmoninarum ATCC 33209 (Rs 33209) cells; (2) killed Rs 33209 cells which had been heated to 37??C for 48 h, a process that destroys the p57 protein; (3) killed R. salmoninarum MT239 (Rs MT239) cells; (4) heated Rs MT239 cells; (5) a recombinant version of the p57 protein (r-p57) emulsified in Freund's incomplete adjuvant (FIA); (6) the commercial BKD vaccine Renogen; (7) phosphate-buffered saline (PBS) emulsified with an equal volume of FIA; or (8) PBS alone. Following injection, each fish was marked with a subcutaneous fluorescent latex tag denoting its treatment group and the vaccinated fish were combined into sham and disease challenge tanks. Two weeks after these fish were vaccinated, separate groups of fish were injected with either PBS or live R. salmoninarum GL64 and were placed inside coated-wire mesh cylinders (liveboxes) in the sham and disease challenge tanks, respectively. Mortalities in both tanks were recorded for 285 d. Any mortalities among the livebox fish were replaced with an appropriate cohort (infected with R. salmoninarum or healthy) fish. None of the bacterins evaluated in this study induced protective immunity against the R. salmoninarum shed from the infected livebox fish. The percentage survival within the test groups in the R. salmoninarum challenge tank ranged from 59% (heated Rs MT239 bacterin) to 81 % (PBS emulsified with FIA). There were no differences in the percentage survival among the PBS-, PBS/FIA-, r-p57-and Renogen-injected groups. There also were no differences in survival among the bacterin groups, regardless of whether the bacterial cells had been heated or left untreated prior to injection. ?? Inter-Research 2005.

Diseases of Aquatic Organisms↗
USGS2020

Arsenolipids in cultured Picocystis strain ML, and their occurrence in biota and sediment from Mono Lake, California

Primary production in Mono Lake, a hypersaline soda lake rich in dissolved inorganic arsenic, is dominated by Picocystis strain ML. We set out to determine if this photoautotrophic picoplankter could metabolize inorganic arsenic and in doing so form unusual arsenolipids (e.g., arsenic bound to 2- O -methyl ribosides) as reported in other saline ecosystems and by halophilic algae. We cultivated Picocystis strain ML on a seawater-based medium with either low (37 µM) or high (1000 µM) phosphate in the presence of arsenite (400 µM), arsenate (800 µM), or without arsenic additions (ca 0.025 µM). Cultivars formed a variety of organoarsenic compounds, including a phytyl 2- O -methyl arsenosugar, depending upon the cultivation conditions and arsenic exposure. When the cells were grown at low P, the organoarsenicals they produced when exposed to both arsenite and arsenate were primarily arsenolipids (~88%) with only a modest content of water-soluble organoarsenic compounds (e.g., arsenosugars). When grown at high P, sequestration shifted to primarily water-soluble, simple methylated arsenicals such as dimethylarsinate; arsenolipids still constituted ~32% of organoarsenic incorporated into cells exposed to arsenate but < 1% when exposed to arsenite. Curiously, Picocystis strain ML grown at low P and exposed to arsenate sequestered huge amounts of arsenic into the cells accounting for 13.3% of the dry biomass; cells grown at low P and arsenite exposure sequestered much lower amounts, equivalent to 0.35% of dry biomass. Extraction of a resistant phase with trifluoroacetate recovered most of the sequestered arsenic in the form of arsenate. Uptake of arsenate into low P-cultivated cells was confirmed by X-ray fluorescence, while XANES/EXAFS spectra indicated the sequestered arsenic was retained as an inorganic iron precipitate, similar to scorodite, rather than as an As-containing macromolecule. Samples from Mono Lake demonstrated the presence of a wide variety of organoarsenic compounds, including arsenosugar phospholipids, most prevalent in zooplankton ( Artemia ) and phytoplankton samples, with much lower amounts detected in the bottom sediments. These observations suggest a trophic transfer of organoarsenicals from the phytoplankton ( Picocystis ) to the zooplankton ( Artemia ) community, with efficient bacterial mineralization of any lysis-released organoarsenicals back to inorganic oxyanions before they sink to the sediments.

California↗
USGS2005

Toxic Alexandrium blooms in the western Gulf of Maine: The plume advection hypothesis revisited

The plume advection hypothesis links blooms of the toxic dinoflagellate Alexandrium fundyense in the western Gulf of Maine (GOM) to a buoyant plume derived from river outflows. This hypothesis was examined with cruise and moored-instrument observations in 1993 when levels of paralytic shellfish poisoning (PSP) toxins were high, and in 1994 when toxicity was low. A coupled physical-biological model simulated hydrography and A. fundyense distributions. Initial A. fundyense populations were restricted to low-salinity nearshore waters near Casco Bay, but also occurred in higher salinity waters along the plume boundary. This suggests two sources of cells - those from shallow-water cyst populations and those transported to shore from offshore blooms in the eastern segment of the Maine coastal current (EMCC). Observations confirm the role of the plume in A. fundyense transport and growth. Downwelling-favorable winds in 1993 transported the plume and its cells rapidly alongshore, enhancing toxicity and propagating PSP to the south. In 1994, sustained upwelling moved the plume offshore, resulting in low toxicity in intertidal shellfish. A. fundyense blooms were likely nutrient limited, leading to low growth rates and moderate cell abundances. These observations and mechanisms were reproduced by coupled physical-biological model simulations. The plume advection hypothesis provides a viable explanation for outbreaks of PSP in the western GOM, but should be refined to include two sources for cells that populate the plume and two major pathways for transport: one within the low-salinity plume and another where A. fundyense cells originating in the EMCC are transported along the outer boundary of the plume front with the western segment of the Maine coastal current.

Casco Bay, Gulf of Maine↗
USGS2013

Usefulness of hemocytometer as a counting chamber in a computer assisted sperm analyzer (CASA)

Several methods are used to determine sperm cell concentration, such as the haemocytometer, spectrophotometer, electronic cell counter and computer-assisted semen analysers (CASA). The utility of CASA systems has been limited due to the lack of characterization of individual systems and the absence of standardization among laboratories. The aims of this study were to: 1) validate and establish setup conditions for the CASA system utilizing the haemocytometer as a counting chamber, and 2) compare the different methods used for the determination of sperm cell concentration in bull semen. Two ejaculates were collected and the sperm cell concentration was determined using spectrophotometer and haemocytometer. For the Hamilton-Thorn method, the haemocytometer was used as a counting chamber. Sperm concentration was determined three times per ejaculate samples. A difference (P < 0.05) was found between all methods of measuring sperm cell concentration. However, no difference was found between the haemocytometer and CASA system when the haemocytometer was used as a counting chamber (P > 0.05) or between the haemocytometer count and the spectrophotometer. Based on the results of this study, we concluded that the haemocytometer can be used in computerized semen analysis systems as a substitute for the commercially available disposable counting chambers, therefore avoiding disadvantageous high costs and slower procedures.

Animal Reproduction↗
USGS1992

Identification of viral hemorrhagic septicemia virus isolated from Pacific cod Gadus macrocephalus in Prince William Sound Alaska, USA

Ulcerative slun tissues from 2 Pacific cod Gadus rnacrocephalus caught in Prince William Sound, Alaska, USA, were examined for virus by Fish Pathology staff within the F.R.E.D. Division of the Alaska Department of Fish and Game. Six days after inoculation of Epitheliorna papulosum cyprini (EPC) cells at 14"C, diffuse rounding and lifting of cells from the monolayers suggestive of cytopathlc effect became visible in the lower sample dilutions. Ultrastructural examinations of affected EPC cells showed rhabdovirus particles within cytoplasmic vacuoles and on the cell surface membranes. Virus isolates from both cod were subsequently confirmed as viral hemorrhagic septicemia virus (VHSV) by serum neutralizabon and immunoblot assay. This is the first VHSV isolated from Pacific cod, which represents a new host species for the virus. Histologically, cod skin ulcers appeared to be caused by a foreign-body-type inflammatory response to foci of protozoa resembling X cells that also had plasmodial stages. Whether the rhabdovirus was incidental to the slun lesion or played a role in its etiology remains to be determined. The possible relationship between thls virus and the recent occurrences of VHSV in anadromous salmoruds from Washington State, USA, is discussed.

Diseases of Aquatic Organisms↗
USGS2001

Snow cover and snow goose Anser caerulescens caerulescens distribution during spring migration

Arctic geese often use spring migration stopover areas when feeding habitats are partially snow covered. Melting of snow during the stopover period causes spatial and temporal variability in distribution and abundance of feeding habitat. We recorded changes in snow cover and lesser snow goose Anser caerulescens caerulescens distribution on a spring migration stopover area in south-central Alaska during aerial surveys in 1993-1994. Our objectives were to determine whether geese selected among areas with different amounts of snow cover and to assess how temporal changes in snow cover affected goose distribution. We also measured temporal changes in chemical composition of forage species after snow melt. We divided an Arc/Info coverage of the approximately 210 km 2 coastal stopover area into 2-km 2 cells, and measured snow cover and snow goose use of cells. Cells that had 10-49.9% snow cover were selected by snow geese, whereas cells that lacked snow cover were avoided. In both years, snow cover diminished along the coast between mid-April and early May. Flock distribution changed as snow geese abandoned snow-free areas in favour of cells where snow patches were interspersed with bare ground. Snow-free areas may have been less attractive to geese because available forage had been quickly exploited as bare ground was exposed, and because soils became drier making extraction of underground forage more difficult. Fiber content of two forage species increased whereas non-structural carbohydrate concentrations of forage plants appeared to diminish after snow melt, but changes in nutrient concentrations likely occurred too slowly to account for abandonment of snow-free areas by snow geese.

Alaska↗