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 325 records · Page 18Linked to original sources
USGS2009

A grid-doubling finite-element technique for calculating dynamic three-dimensional spontaneous rupture on an earthquake fault

We present a new finite-element technique for calculating dynamic 3-D spontaneous rupture on an earthquake fault, which can reduce the required computational resources by a factor of six or more, without loss of accuracy. The grid-doubling technique employs small cells in a thin layer surrounding the fault. The remainder of the modelling volume is filled with larger cells, typically two or four times as large as the small cells. In the resulting non-conforming mesh, an interpolation method is used to join the thin layer of smaller cells to the volume of larger cells. Grid-doubling is effective because spontaneous rupture calculations typically require higher spatial resolution on and near the fault than elsewhere in the model volume. The technique can be applied to non-planar faults by morphing, or smoothly distorting, the entire mesh to produce the desired 3-D fault geometry. Using our FaultMod finite-element software, we have tested grid-doubling with both slip-weakening and rate-and-state friction laws, by running the SCEC/USGS 3-D dynamic rupture benchmark problems. We have also applied it to a model of the Hayward fault, Northern California, which uses realistic fault geometry and rock properties. FaultMod implements fault slip using common nodes, which represent motion common to both sides of the fault, and differential nodes, which represent motion of one side of the fault relative to the other side. We describe how to modify the traction-at-split-nodes method to work with common and differential nodes, using an implicit time stepping algorithm.

Geophysical Journal International↗
USGS2006

Comparison of local grid refinement methods for MODFLOW

Many ground water modeling efforts use a finite-difference method to solve the ground water flow equation, and many of these models require a relatively fine-grid discretization to accurately represent the selected process in limited areas of interest. Use of a fine grid over the entire domain can be computationally prohibitive; using a variably spaced grid can lead to cells with a large aspect ratio and refinement in areas where detail is not needed. One solution is to use local-grid refinement (LGR) whereby the grid is only refined in the area of interest. This work reviews some LGR methods and identifies advantages and drawbacks in test cases using MODFLOW-2000. The first test case is two dimensional and heterogeneous; the second is three dimensional and includes interaction with a meandering river. Results include simulations using a uniform fine grid, a variably spaced grid, a traditional method of LGR without feedback, and a new shared node method with feedback. Discrepancies from the solution obtained with the uniform fine grid are investigated. For the models tested, the traditional one-way coupled approaches produced discrepancies in head up to 6.8% and discrepancies in cell-to-cell fluxes up to 7.1%, while the new method has head and cell-to-cell flux discrepancies of 0.089% and 0.14%, respectively. Additional results highlight the accuracy, flexibility, and CPU time trade-off of these methods and demonstrate how the new method can be successfully implemented to model surface water-ground water interactions. Copyright ?? 2006 The Author(s).

Conference Paper↗
USGS2004

Comparison of fluorescence microscopy and solid-phase cytometry methods for counting bacteria in water

Total direct counts of bacterial abundance are central in assessing the biomass and bacteriological quality of water in ecological and industrial applications. Several factors have been identified that contribute to the variability in bacterial abundance counts when using fluorescent microscopy, the most significant of which is retaining an adequate number of cells per filter to ensure an acceptable level of statistical confidence in the resulting data. Previous studies that have assessed the components of total-direct-count methods that contribute to this variance have attempted to maintain a bacterial cell abundance value per filter of approximately 10 6 cells filter -1 . In this study we have established the lower limit for the number of bacterial cells per filter at which the statistical reliability of the abundance estimate is no longer acceptable. Our results indicate that when the numbers of bacterial cells per filter were progressively reduced below 10 5 , the microscopic methods increasingly overestimated the true bacterial abundance (range, 15.0 to 99.3%). The solid-phase cytometer only slightly overestimated the true bacterial abundances and was more consistent over the same range of bacterial abundances per filter (range, 8.9 to 12.5%). The solid-phase cytometer method for conducting total direct counts of bacteria was less biased and performed significantly better than any of the microscope methods. It was also found that microscopic count data from counting 5 fields on three separate filters were statistically equivalent to data from counting 20 fields on a single filter.

Applied and Environmental Microbiology↗
USGS2017

In-vitro replication of Chelonid herpesvirus 5 in organotypic skin cultures from Hawaiian green turtles (Chelonia mydas)

Fibropapillomatosis (FP) is a tumor disease of marine turtles associated with Chelonid herpesvirus 5 (ChHV5) that has historically been refractory to growth in tissue culture. Here, we show for the first time de novo formation of ChHV5-positive intranuclear inclusions in cultured green turtle cells, which is indicative for active lytic replication of the virus. The minimal requirements to achieve lytic replication in cultured cells included 1) either in-vitro culturing of ChHV5-positive tumor biopsies (plugs) or organotypic cultures (rafts) consisting of ChHV5-positive turtle fibroblasts in collagen rafts seeded with turtle keratinocytes and 2) keratinocyte maturation induced by raising raft or biopsy cultures to the air-liquid interface. Virus growth was confirmed by detailed electron microscopic studies revealing intranuclear sun-shaped capsid factories, tubules, various stages of capsid formation, nuclear export by budding into the perinuclear space, tegumentation, and envelopment to complete de novo virus production. Membrane synthesis was also observed as a sign for active viral replication. Interestingly, cytoplasmic particles became associated with keratin filaments, a feature not seen in conventional monolayer cell cultures where most studies of herpesvirus replication have been performed. Our findings draw a rich and realistic picture of ChHV5 replication in cells derived from its natural host and may be crucial not only to better understand ChHV5 circulation but also to eventually complete Koch's postulates for FP. Moreover, the principles described here may serve as model to culture other viruses that are resistant to replication in conventional cell culture.

Hawai'i↗
USGS1991

Nature and time course of acclimation to aluminum in juvenile brook trout ( Salvelinus fontinalis ): II. Gill histology

Gill samples from juvenile brook trout ( Salvelinus fontinalis ) acclimated to low-level aluminum at pH 5.2 showed severe damage by day 4, with necrosis and fusion of secondary lamellae and hyperplasia and hypertrophy of mucous cells. Over the following 20 d, there was a continual process of repair with proliferation and hypertrophy of mucous cells. Qualitative analysis of gill samples plus physiology and mortality data collected in a companion study indicated progressive development (by day 10 onward) of increasing acclimation to Al. Quantitative analysis of gill samples on day 13 showed that mucous cell volume density had tripled and mucous cell area had doubled in Al-exposed fish compared with control fish. A lamellar fusion index showed evidence of fusion in Al-exposed fish by day 4 with recovery to nearly control levels by day 13. Physiological disturbances appear to be directly related to the histological changes observed in the gill epithelium. At the cellular level, changes in either mucous cell production and secretion or changes in mucus chemistry contribute, in part, to acclimation to Al.

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2005

Characteristics of seasonal vegetation cover in the conterminous USA

A data set of the fractional green vegetation cover (FGREEN) for the Conterminous USA was evaluated for regional and seasonal variation. The value of FGREEN was derived monthly for the three most dominant land cover classes per 20 km by 20 km grid cell within the study area. At this grid cell resolution (comprised of 400 1-km pixels), 97 percent of the grid cells included three or fewer land cover classes. FGREEN was found to vary regionally due to local land cover and climate variations. FGREEN was found significantly different between one or more of the land cover classes, for one or more months, in 58 percent of the grid cells included in the study. Monthly FGREEN values for the land cover classes vary sufficiently between the land cover classes to warrant monthly FGREEN data for each of the one to three most dominant land cover classes per grid cell .

Photogrammetric Engineering and Remote Sensing↗
USGS1997

A monoclonal antibody to inclusion body disease of cranes virus enabling specific immunohistochemistry and competitive ELISA

Inclusion body disease of cranes (IBDC) herpesvirus kills some infected cranes and persists in convalescent animals. To enable further study and rapid identification of carrier animals, we developed a monoclonal antibody (MAb) to IBDC virus and used it in immunohistochemistry and a competitive enzyme-linked immunosorbent assay (ELISA). We used conventional techniques to make murine MAbs directed against IBDC virus purified from infected duck embryo cells. Hybridomas reacting in an ELISA with IBDC virus but not uninfected duck embryo cells were characterized by radioimmunoprecipitation, in situ immunohistochemistry, and competitive ELISA with neutralizing and nonneutralizing crane sera. MAb 2C11 immunoprecipitated 59-, 61-, and 110-kD proteins from IBDC virus-infected but not uninfected cells and stained glutaraldehyde-fixed IBDC virus plaques but not surrounding uninfected duck embryo cells in vitro. Antibody 2C11 did not react with duck embryo cells infected with falcon herpesvirus, psittacine herpesvirus, infectious laryngotracheitis, pigeon herpesvirus, or duck plague virus. A competitive ELISA using antibody 2C11 identified most sera that were positive in the neutralization test. This antibody will be useful in further characterizing IBDC virus, its pathogenesis, and its natural history.

Avian Diseases↗
USGS2002

Documentation and analysis of a geographic information system application for combining data layers, using nonpoint-source pollution as an example

A geographical information system application has been developed that allows scientists to combine multiple data layers into a single data layer. This application provides an effective tool for identifying areas where the potential effect of the combination of data layers may be greater than any single data layer. Such a tool is useful in studying an activity that cannot be measured directly. Scientists wanting to identify areas where an activity may have the greatest effect can identify factors that directly or indirectly reflect the effect of the activity being studied. When combined, these factors would identify areas where the potential for the activity to have an effect are greatest. The data layers used to develop the single data layer determine the activity addressed by the application?the application was developed to identify areas where the potential for nonpoint-source pollution to affect areas of Indiana is greater relative to other areas in Indiana. To evaluate the potential in other states or areas, data layers for those states or areas would be used. To address a different activity, even activities not related to water resources, data layers that directly or indirectly reflect the effects of the activity being studied would be used in the application. The application was developed using Environmental System Research Institute's ArcView geographical information system and the ArcView extension, Spatial Analyst. To use the application, the user selects data layers related to the activity being studied, describes the variability of a data element between geographic areas on each data layer, assigns a relative importance factor to each data layer, and focuses the output by identifying a watershed of interest. The application then assigns a rank value to the features within each data layer on the basis of the spatial variation of the data layer. The data layers are converted to grids. The values of the grid cells are the product of the rank values and the importance factors. The cells for each data layer are combined to form a single grid. The combination is summarized to describe the spatial variation for the watershed of interest. The combined data layer does not show the actual potential effect of the selected factors but rather the relative difference in the potential effect among areas when all data layers are considered. An analysis of the application indicates that the selected data layers to be combined should be at the greatest spatial resolution possible; however, all data layers do not have to be at the same spatial resolution. The spatial variation of the data layers should be adequately defined. The size of each grid cell should be small enough to maintain the spatial definition of smaller features within the data layers. The most accurate results are shown to occur when the values for the grid cells representing the individual data layers are summed and the mean of the summed grid-cell values is used to describe the watershed of interest.

Indiana↗
USGS2002

MODFLOW-2000, the U.S. Geological Survey modular ground-water model -- Documentation of the Model-Layer Variable-Direction Horizontal Anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) package

This report documents the model-layer variable-direction horizontal anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) Package of MODFLOW-2000. The LVDA capability allows the principal directions of horizontal anisotropy to be different than the model-grid row and column directions, and for the directions to vary on a cell-by-cell basis within model layers. The HUF Package calculates effective hydraulic properties for model grid cells based on hydraulic properties of hydrogeologic units with thicknesses defined independently of the model layers. These hydraulic properties include, among other characteristics, hydraulic conductivity and a horizontal anisotropy ratio. Using the LVDA capability, horizontal anisotropy direction is defined for model grid cells within which one or more hydrogeologic units may occur. For each grid cell, the HUF Package calculates the effective horizontal hydraulic conductivity along the primary direction of anisotropy using the hydrogeologic-unit hydraulic conductivities, and calculates the effective horizontal hydraulic conductivity along the orthogonal anisotropy direction using the effective primary direction hydraulic conductivities and horizontal anisotropy ratios. The direction assigned to the model layer effective primary hydraulic conductivity is specified using a new data set defined by the LVDA capability, when active, to calculate coefficients needed to solve the ground-water flow equation. Use of the LVDA capability is illustrated in four simulation examples, which also serve to verify hydraulic heads, advective-travel paths, and sensitivities calculated using the LVDA capability. This version of the LVDA capability defines variable-direction horizontal anisotropy using model layers, not the hydrogeologic units defined by the HUF Package. This difference needs to be taken into account when designing model layers and hydrogeologic units to produce simulations that accurately represent a given field problem. This might be a reason, for example, to make model layer boundaries coincide with hydrogeologic-unit boundaries in all or part of a model grid.

Open-File Report↗
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↗
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↗