Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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 1,513 records · Page 84Linked to original sources
USGS2020

Groundwater quality in relation to drinking water health standards and hydrogeologic and geochemical characteristics for 47 domestic wells in Potter County, Pennsylvania, 2017

As part of a regional effort to characterize groundwater in rural areas of Pennsylvania, water samples from 47 domestic wells in Potter County were collected from May through September 2017. The sampled wells had depths ranging from 33 to 600 feet in sandstone, shale, or siltstone aquifers. Groundwater samples were analyzed for physicochemical properties that could be evaluated in relation to drinking-water health standards, geology, land use, and other environmental factors. Laboratory analyses included concentrations of major ions, nutrients, bacteria, trace elements, volatile organic compounds (VOCs), ethylene and propylene glycol, alcohols, gross-alpha/beta-particle activity, uranium, radon-222, and dissolved gases. A subset of samples was analyzed for radium isotopes (radium-226 and -228) and for the isotopic composition of methane. Results of this 2017 study show that groundwater quality generally met most drinking-water standards that apply to public water supplies. However, a percentage of samples exceeded maximum contaminant levels (MCLs) for total coliform bacteria (69.6 percent), Escherichia coli (30.4 percent), arsenic, and barium; and secondary maximum contaminant levels (SMCLs) for field pH, manganese, sodium, iron, total dissolved solids, aluminum, and chloride. All of the analyzed VOCs were below limits of detection and associated drinking water criteria. Radon-222 activities exceeded the proposed drinking-water standard of 300 picocuries per liter in 80.9 percent of the samples. The field pH of the groundwater ranged from 4.6 to 9.0. Generally, the lower pH samples had greater potential for elevated concentrations of dissolved metals, including beryllium, copper, lead, nickel, and zinc, whereas the higher pH samples had greater potential for elevated concentrations of total dissolved solids, sodium, fluoride, boron, and uranium. Near-neutral samples (pH 6.5 to 7.5) had greater hardness and alkalinity concentrations than other samples with pH values outside this range. Calcium/bicarbonate waters were the predominant hydrochemical type for the sampled aquifers, with mixed water types for many samples, including variable contributions from calcium, magnesium, and sodium combined with bicarbonate, sulfate, chloride, and nitrate. Water from 45 wells had concentrations of methane greater than the 0.0002 milligrams per liter (mg/L) detection limit. One sample had the maximum value of 11 mg/L, which exceeds the Pennsylvania action level of 7 mg/L. Additionally, three other samples had concentrations of methane greater than 4 mg/L. Outgassing of such levels of methane from the water to air within a confined space can result in a potential hazard. The elevated concentrations of methane generally were associated with suboxic groundwater (dissolved oxygen less than 0.5 mg/L) that had near-neutral to alkaline pH with relatively elevated concentrations of iron, manganese, ammonia, lithium, fluoride, and boron. Other constituents, including barium, sodium, chloride, and bromide, commonly were elevated, but not limited to, those well-water samples with elevated methane. Low levels of ethane (as much as 1.2 mg/L) were present in eight samples with the highest methane concentrations. Five samples were analyzed for methane isotopes. The isotopic and hydrocarbon compositions in these five samples suggest the methane may be of microbial origin or a mixture of thermogenic and microbial gas, but differed from the compositions reported for mud-gas logging samples collected during drilling of gas wells. The concentrations of sodium (median 8.2 mg/L), chloride (median 7.64 mg/L), and bromide (median 0.02 mg/L) for the 47 groundwater samples collected for this study ranged widely and were positively correlated with one another and with specific conductance and associated measures of ionic strength. Sixty percent of the Potter County well-water samples had chloride concentrations less than 10 mg/L. Samples with higher chloride concentrations had variable bromide concentrations and corresponding chloride/bromide ratios that are consistent with sources such as road-deicing salt and septic effluent (low bromide) or brine (high bromide). Brines are naturally present in deeper parts of the regional groundwater system and, in some cases, may be mobilized by gas drilling. It is also possible that valley wells were drilled close to or into the brine-freshwater interface, so brine signatures do not necessarily indicate contamination due to drilling. The chloride, bromide, and other constituents in road-deicing salt or brine solutions tend to be diluted by mixing with fresh groundwater in shallow aquifers used for water supply. Although 1 of 8 groundwater samples with the highest methane concentrations (greater than 0.2 mg/L) had concentrations of chloride and bromide with corresponding chloride/bromide ratios that indicated mixing with road-deicing salt, the other 7 of 8 samples with elevated methane had concentrations of chloride and bromide with corresponding chloride/bromide ratios that indicated mixing with a small amount of brine (0.02 percent or less) similar in composition to those reported for gas and oil well brines in Pennsylvania. In several eastern Pennsylvania counties where gas drilling is absent, groundwater with comparable chloride/bromide ratios and chloride concentrations have been reported. Approximately 50 percent of Potter County well-water samples, including two samples with the fourth (72.9 mg/L) and fifth (47.0 mg/L) highest chloride concentrations, have chloride/bromide ratios that indicate predominantly anthropogenic sources of chloride, such as road-deicing salt or septic effluent.

Pennsylvania↗
USGS2020

Procedure for calculating estimated ultimate recoveries of wells in the Wolfcamp shale of the Midland Basin, Permian Basin Province, Texas

In 2016, the U.S. Geological Survey published an assessment of technically recoverable continuous oil and gas resources of the Wolfcamp shale in the Midland Basin, Permian Basin Province, Texas. Estimated ultimate recoveries (EURs) were calculated with production data from IHS Markit TM using DeclinePlus software in the Harmony interface. These EURs were a major component of the quantitative resource assessment. For five of the six assessment units in the study, an industry operator in the Midland Basin provided information that was used to differentiate the Wolfcamp horizontal well landing zones. The IHS Markit TM production database does not distinguish between the Wolfcamp A, B, C, and D well landing zones. These different units of the Wolfcamp have different production patterns that are important for calculation of EURs. The calculated mean EURs for each assessment unit ranged from 99,000 barrels of oil in the Wolfcamp C to 142,000 barrels of oil in the Wolfcamp A.

Texas↗
USGS2020

Steps taken for calculating estimated ultimate recoveries of wells in the Eagle Ford Group and associated Cenomanian–Turonian strata, U.S. Gulf Coast, Texas, 2018

In 2018, the U.S. Geological Survey published an assessment of technically recoverable continuous oil and gas resources of the Eagle Ford Group and associated Cenomanian–Turonian strata in the U.S. Gulf Coast of Texas. Estimated ultimate recoveries (EURs) were calculated with production data from IHS Markit TM using DeclinePlus software in the Harmony interface. These EURs were a major component of the aforementioned quantitative resource assessment fact sheet. The calculated mean EURs for each oil assessment unit (AU) ranged from 113,000 barrels of oil in the Cenomanian–Turonian Mudstone Continuous Oil AU to 223,000 barrels of oil in the Submarine Plateau-Karnes Trough Continuous Oil AU. The calculated mean EURs for each gas AU ranged from 2.261 billion cubic feet of gas in the Submarine Plateau-Karnes Trough Continuous Gas AU to 3.116 billion cubic feet of gas in the Eagle Ford Marl Continuous Gas AU.

Texas↗
USGS2021

Simulation of water-table response to sea-level rise and change in recharge, Sandy Hook unit, Gateway National Recreation Area, New Jersey

The Sandy Hook Unit, Gateway National Recreation Area (hereafter Sandy Hook) in New Jersey is a 10-kilometer-long spit visited by thousands of people each year who take advantage of the historical and natural resources and recreational opportunities. The historical and natural resources are threatened by global climate change, including sea-level rise (SLR), changes in precipitation and groundwater recharge, and changes in the frequency and severity of coastal storms. Fresh groundwater resources are important to the ecosystems of Sandy Hook. The Bayside Holly Forest, one of only two known old-growth American holly ( Ilex opaca ) maritime forests, is particularly vulnerable to global climate change because of the proximity of the water table to land surface in low-lying areas and the potential for saltwater intrusion and inundation. The shallow groundwater-flow system on Sandy Hook is dominated by recharge from precipitation, fresh groundwater discharge to evapotranspiration (ET), discharge to surface seeps, and submarine groundwater discharge (groundwater discharging directly to the ocean). A three-dimensional groundwater-flow model that simulates the shallow groundwater-flow system and interaction with surrounding saltwater boundaries was constructed to simulate multi-density groundwater flow, treating the freshwater/saltwater transition zone as a sharp interface that represents the half-seawater surface. Groundwater-flow simulations completed for this study include a Baseline scenario, three SLR scenarios (0.2, 0.4, and 0.6 meter [m]), two Recharge scenarios—a 10-percent Increased Recharge scenario and a 10-percent Decreased Recharge scenario—and a scenario with 0.6 m of SLR and 10-percent increase in recharge. The Recharge scenarios indicate the system is not sensitive to a 10-percent increase or decrease in recharge from the Baseline scenario. In the SLR scenarios, SLR causes the water table to rise, resulting in increased fresh groundwater discharge to ET and seeps, and reduced submarine discharge compared to the Baseline scenario. The increased discharge to ET and seeps causes the magnitude of water-table rise to be less than that of SLR, which in turn causes the thickness of the freshwater lens to thin, reducing the depth to the half-seawater surface. Water-table rise associated with SLR diminishes the thickness of the unsaturated zone; comparing the Baseline and the 0.6-m SLR scenarios, the area where the simulated water table is above land surface increases by 50.6 hectares, from about 0.9 to 7.4 percent of the land area of Sandy Hook. Areas where the simulated water table is above land surface are likely to be emergent wetlands and contain freshwater if they are tens of meters or more from the shoreline. The steady-state simulations indicate that the percentage of land where the half-seawater surface is less than 9 m below the water table increases from about 2.5 percent (20 hectares) to about 9 percent (74 hectares) with 0.6 m of SLR. In low-lying areas close to the Sandy Hook Bay shoreline, the half-seawater surface is simulated to be as much as 20 m closer to the water table with SLR of 0.6 m. Transient salinization, if any, of shallow groundwater from increased frequency or severity of storm-driven inundation is not included in the analysis. Natural resources on Sandy Hook, particularly the Bayside Holly Forest, may be adversely affected by the rising water table associated with SLR. Freshwater emergent wetlands may increase in area at the expense of other ecosystem assemblages occurring in or on the edges of low-lying enclosed depressions. Cultural resources close to the water table, such as existing basements of structures, may be adversely affected.

New Jersey↗
USGS2021

Characterization of groundwater quality and discharge with emphasis on selenium in an irrigated agricultural drainage near Delta, Colorado, 2017–19

Selenium is a water-quality constituent of concern for aquatic ecosystems in the lower Gunnison River Basin. Selenium is derived from bedrock of the Mancos Shale and is mobilized and transported to groundwater and surface water by application of irrigation water. Although it is recognized that groundwater contributes an appreciable amount of selenium to surface water, few studies have addressed interactions between the two. The U.S. Geological Survey in cooperation with the Colorado Water Conservation Board conducted a study during 2017–19 to characterize the quality and quantity of groundwater discharging to an agricultural drainage near Delta, Colorado, locally known as Sunflower Drain. Water quality in the study area is characterized by high dissolved solids with elevated concentrations of selenium and nitrate resulting from dissolution of soluble salts in the Mancos Shale. Selenium concentrations have decreased by 50 percent since the early 2000s, possibly in response to irrigation system improvements. Stable water isotopes indicate streamflow is dominated by canal water during the irrigation season (April to October) and, during the nonirrigation season (November to March), is dominated by groundwater that has undergone some degree of evaporation. Pesticide and pharmaceutical compounds were infrequently detected, and results indicate they were derived from sources outside the study area such that they do not appear to be useful as tracers of groundwater sources. Stable isotopes of nitrate indicate that nitrate originates from the Mancos Shale, and the isotopic composition is enriched by denitrification in the groundwater system. Using a mass-balance approach, estimated groundwater discharge rates to Sunflower Drain ranged from 0.15 to 0.27 cubic feet per second per mile with one losing reach identified. Selenium, sulfate, and nitrate concentrations in groundwater estimated by mass-balance calculations were similar to concentrations measured in the Poly 17 observation well, located in a largely irrigated area in east tributary. One tributary reach had higher concentrations of selenium, sulfate, and nitrate likely reflecting localized inputs of more concentrated groundwater, similar to the concentrations in the Poly 7 observation well, which is downgradient from a residential area in the west tributary. Three pilot studies were conducted, including fiber optic distributed temperature sensing to detect groundwater discharge zones in the stream channel, a passive seismic technique to estimate depth to bedrock, and use of radon-222 as a geochemical tracer of groundwater discharge. All three techniques show promise as additional approaches for investigating groundwater discharge surface-water systems in irrigated drainage areas on Mancos Shale. The factors that affect groundwater movement mainly include when and where irrigation water is transported and applied, and the distribution of bedrock of the Mancos Shale and overlying alluvial deposits. The average groundwater recharge rate for the study area was estimated at 8.1 inches per year, based on mass balance calculations from synoptic survey data. Along the western tributary of Sunflower Drain, there was evidence that spills from the East Canal may recharge the groundwater aquifer adjacent to the stream channel. Groundwater movement to the stream channel may be controlled by the topography of the alluvial/bedrock interface or focused along human-made features, such as tile drains and ditches constructed around irrigated fields. On larger scales, the top of bedrock was also important, creating a topographic constriction that caused a zone of groundwater discharge. The groundwater system is complex, and further study could better define the system, possibly through application of a groundwater flow model and more extensive studies using some of the exploratory methods evaluated in this study.

Colorado↗
USGS2021

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗
USGS2021

Precipitation-runoff processes in the Merced River Basin, Central California, with prospects for streamflow predictability, water years 1952–2013

The U.S. Geological Survey, in cooperation with the California Department of Water Resources (DWR), has constructed a new spatially detailed Precipitation-Runoff Modeling System (PRMS) model for the Merced River Basin, California, which is a tributary of the San Joaquin River in California. Operated through an Object User Interface (OUI) with Ensemble Streamflow Prediction (ESP) and daily climate distribution preprocessing functionality, the model is calibrated primarily to simulate (and eventually, forecast) year-to-year variations of inflows to Lake McClure during the critical April–July snowmelt season. The model is intended to become part of a suite of methods used by DWR for estimating daily streamflow from the Merced River Basin, especially during the snowmelt season. This study describes the results of the application of an analysis tool that simulates responses to climate and land-use variations at a higher spatial resolution than previously available to DWR. A geographic information system was used to delineate the model domain, that is, areas draining to a single outlet at U.S. Geological Survey streamflow-gaging station 11270900, Merced River below Merced Falls Dam, near Snell, CA (also known as California Data Exchange Center station MRC), and subdrainage areas, including four draining to internal gages used as calibration targets. Using this delineation, three contiguous subbasins were recognized and, along with the model domain and nested calibration targets, are the simulation units evaluated in this report. An auto-calibration tool, LUCA (Let Us CAlibrate), was used for each calibration node, from headwaters to basin outlet, and then parameters were manually adjusted to complete the calibration. The main objective was to match April–July snowmelt seasonal discharge values of simulated streamflow to observed (measured or reconstructed) discharge values. Calibration or validation periods used site-specific streamflows—mostly from October 1, 1988, through September 30, 2013—but differed according to the period-of-record available for the measurements collected at internal gages or reconstructed flows for the single outlet. The accuracy of the Merced PRMS streamflow simulations varied seasonally, as compared to observed values. Based on statistical results, the Merced PRMS model satisfactorily simulated snowmelt seasonal streamflows. April–July calibrations for all areas had small negative bias (not greater than 7 percent) and low relative error (less than 8 percent). Less satisfactory performance for other seasons was attributed to several factors: (1) high uncertainty in low or zero flows in summer and fall, (2) lack of accounting for basin withdrawals and anthropogenic water use, (3) unavailability and (or) inaccuracy of observed (measured) meteorological input data, and (4) uncertainty in reconstructed streamflow data. With some additional refinement, the Merced PRMS model may be used for forecasting seasonal and longer-term streamflow variations; evaluating forecasted and past climate and land cover changes; providing water-resource managers with a consistent and documented method for estimating streamflow at ungaged sites within the basin; and aiding environmental studies, hydraulic design, water management, and water-quality projects in the Merced River Basin.

California↗
USGS2024

Hydrogeologic framework and extent of saltwater intrusion in Kings, Queens, and Nassau Counties, Long Island, New York

In 2016, the U.S. Geological Survey began a multiyear cooperative study with the New York State Department of Environmental Conservation to evaluate the sustainability of Long Island’s sole-source aquifer system through hydrogeologic mapping, compilation of groundwater chloride concentrations, and groundwater flow modeling. In the initial phase of the islandwide study, the hydrogeologic framework and extent of saltwater intrusion in aquifers in Kings, Queens, and Nassau Counties on western Long Island, N.Y., were investigated. The aquifer system underlying western Long Island has been under stress from pumping of public, irrigation (golf course), and industrial supply wells. Saltwater intrusion has occurred from surrounding embayments (East River, Long Island Sound, Jamaica Bay, and the Atlantic Ocean) due to pumping. Eighteen boreholes were drilled and cores taken during 2019–21 to collect hydrogeologic, geochemical, and geophysical data to delineate the complex subsurface hydrogeology and extent of saltwater intrusion within the study area. Evaluation of the new cores, reexamination of legacy core descriptions, and analysis of borehole geophysical logs was used to refine the previously published hydrogeologic framework of Pleistocene and Cretaceous unconsolidated sediments in the area, including delineation of a previously undefined hydrogeologic unit between the Magothy aquifer and the Raritan confining unit, herein named the “upper Raritan aquifer.” The upper Raritan aquifer was first recognized in southeastern Nassau County from an analysis of about 50 closely spaced boreholes with high-resolution core descriptions and gamma-ray (gamma) logs. Further analysis of borehole logs across the study area indicated that the upper Raritan aquifer was also present in Kings and Queens Counties. Nuclear magnetic resonance (NMR) logging was used for the first time on Long Island to provide estimates of the hydraulic properties of the major aquifer and confining units. Unlike other geophysical logs that record responses to the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR logs record responses to the presence of hydrogen protons in the formation fluid to determine water fraction and pore-size distribution. NMR log analysis provided estimates of the clay-bound, capillary-bound, and mobile water fractions and hydraulic conductivity of aquifers and confining units penetrated by five wells in Nassau County. Pumpage for public-supply and industrial wells on Long Island began in the 1870s with small, localized suppliers of populated areas in Kings and Queens Counties. By 1904–16, pumpage for public water supply in Kings County averaged 21 million gallons per day, and averaged 37 million gallons per day in Queens County, mostly from the upper glacial aquifer. Saltwater intrusion was reported as early as the beginning of the 20th century and included the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifer systems. By 1936, pumping in central Kings County created a major cone of depression in the water table extending to the south shore of much of Kings County and into southwestern Queens County. Saltwater intrusion has caused the shutdown of public-supply wells in Kings, Queens, and Nassau Counties. A large saltwater intrusion wedge in the Lloyd aquifer was indicated in southern Queens County in the early part of the 20th century, and the saltwater interface may have been onshore predevelopment. Most of Kings and Queens Counties are intruded with saltwater in both the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifers systems. Saltwater increased during the 20th century and continues to increase to the present (2023) in the Lloyd-North Shore aquifer system in Great Neck and Manhasset Neck in northern Nassau County. A major wedge of saltwater intrusion in the upper glacial-Jameco-Magothy aquifer in southwestern Nassau County appears to be increasing.

New York↗
USGS2024

Water-quality monitoring strategy for Mount Hope Bay and the Taunton River Estuary, southeastern Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection (MassDEP), began a study in 2018 to develop a water-quality monitoring strategy (WQMS) for Mount Hope Bay and the Taunton River Estuary in southeastern Massachusetts. MassDEP is interested in water-quality data in Mount Hope Bay and the Taunton River Estuary to characterize current water-quality conditions, assess nutrient-related effects, and capture conditions before and after planned upgrades to wastewater treatment facilities. The WQMS provides an overview of the environmental setting of Mount Hope Bay and the Taunton River Estuary and of dissolved oxygen and nutrient-related water-quality issues, reviews historical and existing monitoring data, and provides recommendations for future monitoring designed to meet five MassDEP management objectives: (1) support MassDEP’s review of coastal and marine dissolved oxygen criteria, (2) assess conditions within MassDEP waterbody assessment units in Mount Hope Bay and the Taunton River Estuary with respect to selected criteria in the Massachusetts Surface Water Quality Standards, (3) assess conditions along the freshwater/saltwater interface (salt wedge) of the Taunton River Estuary and delineate the boundary between the freshwater and saltwater waterbody assessment units, (4) estimate data requirements needed to determine nutrient loads flowing into Mount Hope Bay and the Taunton River Estuary, and (5) evaluate data requirements needed to support hydrodynamic and water-quality models for Mount Hope Bay and the Taunton River Estuary. This WQMS may be used by MassDEP to develop a statewide approach for monitoring estuaries in Massachusetts.

Massachusetts↗
USGS2024

Assessment of the interconnection between Tampa Bay and the Floridan aquifer system: Historical groundwater data compilation and analysis, 1976–2022

The U.S. Geological Survey used existing data collected after the last major navigational channel modification in the mid-1980s to investigate groundwater levels and chloride concentrations in wells in the Floridan aquifer system and other aquifers beneath and near Tampa Bay. Tampa Bay is located on the west-central coast of Florida and provides access for commercial shipping. In 2021, the U.S. Army Corps of Engineers began to investigate alternatives to improve the efficiency of the deep-draft navigation channels within Tampa Bay, reduce costs, improve safety, and analyze the implications of modifying navigational channels. The Floridan aquifer system underlies Tampa Bay and is the primary source of public water supply in the region. Modifications to the channels have the potential to expose the Floridan aquifer system to the saltwater in Tampa Bay, with the potential to increase salinity in wells in the region. Other factors affecting the groundwater levels and the location of the freshwater/saltwater interface were also examined, including changes in sea level, groundwater extraction, and variations in climate. Groundwater levels and well-construction reports were used to identify whether different aquifer units are well-connected. Twelve wells had available data before and after the last major channel modifications, which took place in the 1980s, with six datasets of chloride concentration available in areas along the northern and eastern coastline of Tampa Bay, which is nearest to historical dredging activities. Of these six, Kendall’s τ and p -values indicated increasing trends in chloride concentration for three datasets (TR 11-2, TR 10-2, 51), no trend in chloride concentration for two datasets (TR 9-3, 50), and a decreasing trend in chloride concentration for one dataset (TR 9-1). The upward trends in chloride concentration observed for TR 10-2 and 51 are likely the result of changes in local groundwater withdrawals. Well TR 11-2 had a gradual increasing trend in chloride concentration, fresh groundwater throughout the period of record, and a 3- to 4-foot increase in hydraulic head during the period of record, possibly caused by the construction and control of the Tampa Bypass Canal, resulting in changes to the regional potentiometric surface.

Florida↗
USGS2026

Applying a metapopulation model to compare native and invasive fish passage scenarios in the Upper Mississippi River

The U.S. Geological Survey assessed the ability of metapopulation models to evaluate population-level effects of fish passage projects in the Upper Mississippi River for six native species: Ictalurus furcatus (Valenciennes in Cuvier and Valenciennes, 1840; blue catfish), Ictiobus cyprinellus (Valenciennes in Cuvier and Valenciennes, 1844; bigmouth buffalo), Acipenser fulvescens (Rafinesque, 1817; lake sturgeon), Polyodon spathula (Walbaum, 1792; American paddlefish [paddlefish]), Ictiobus bubalus (Rafinesque, 1818; smallmouth buffalo), Scaphirhynchus platorynchus (Rafinesque, 1818; shovelnose sturgeon), and one invasive species: Hypophthalmichthys molitrix (Valenciennes in Cuvier and Valenciennes, 1844; silver carp). The study had four objectives: to (1) compile parameter values and identify gaps, (2) parameterize species-specific models, (3) apply the models to multiple scenarios, and (4) develop a point-and-click interface. The literature review revealed a paucity of data, for all species except silver carp and paddlefish; however, length–weight relations, female length–maturity curves, growth curves, spawning and recruitment estimates, and movement data among pools were identified. These data allowed the parameterization of native species metapopulation models, in addition to the use of an existing silver carp model. The models were applied to scenarios representing increased passage, decreased passage, combined passage and deterrents effects, and increased native passage with decreased invasive passage. Increased passage scenarios (depicting potential passage projects) produced greater movement and larger upriver populations. Decreased passage scenarios (depicting potential deterrent projects) yielded reduced movement. Similar population dynamic patterns were observed across species for each scenario because of a lack of species-specific empirical movement data. Simulations for longer lived, slow growing species required longer response times, including time to reach equilibria. A point-and-click model was created to allow others to use the model. Data availability for movement limited the ability of the model to be used to understand fish-movement biology and thereby inform management strategies. More fish-movement data could improve model projections of these scenarios.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗
USGS2006

Determination of the δ 34 S of sulfate in water; RSIL lab code 1951

The purpose of the Reston Stable Isotope Laboratory (RSIL) lab code 1951 is to determine the δ( 34 S/ 32 S), abbreviated as δ 34 S, of dissolved sulfate. Dissolved sulfate is collected in the field and precipitated with BaCl 2 at pH 3 to 4 as BaSO 4 in the laboratory. However, the dissolved organic sulfur (DOS) is oxidized to SO 2 , and the carbonate is acidified to CO 2 . Both are degassed from the water sample before the sulfate is precipitated. The precipitated BaSO 4 is filtered and dried before introduction into an elemental analyzer (EA) Carlo Erba NC 2500. The EA is used to convert sulfur in a BaSO 4 solid sample into SO 2 gas, and the EA is connected to a continuous flow isotope-ratio mass spectrometer (CF-IRMS), which determines the differences in the isotope-amount ratios of stable sulfur isotopes ( 34 S/ 32 S) of the product SO 2 gas. The combustion is quantitative; no isotopic fractionation is involved. Samples are placed in a tin capsule and loaded into the Costech Zero Blank Autosampler of the EA. Under computer control, samples are dropped into a heated tube reaction tube that combines the oxidation and reduction reactions. The combustion takes place in a helium atmosphere containing an excess of oxygen gas at the oxidation zone at the top of the reaction tube. Combustion products are transported by a helium carrier through the reduction zone at the bottom of the reaction tube to remove excess oxygen and through a separate drying tube to remove any water. The gas-phase products, mainly CO 2 , N 2 , and SO 2 , are separated by a gas chromatograph. The gas is then introduced into the isotope-ratio mass spectrometer (IRMS) through a Finnigan MAT (now Thermo Scientific) ConFlo II interface, which also is used to inject SO 2 reference gas and helium for sample dilution. The IRMS is a Thermo Scientific Delta V Plus CF-IRMS. It has a universal triple collector with two wide cups and a narrow cup in the middle. It is capable of measuring mass/charge ( m/z ) 64 and 66 simultaneously. The ion beams from SO 2 are as follows: m/z 64 = SO 2 = 32 S 16 O 16 O; m/z 66 = SO 2 = 34 S 16 O 16 O primarily.

Techniques and Methods↗
USGS2006

Determination of the δ 15 N of total nitrogen in solids; RSIL lab code 2893

The purpose of the Reston Stable Isotope Laboratory (RSIL) lab code 2893 is to determine the δ( 15 N/ 14 N), abbreviated as δ 15 N , of total nitrogen in solid samples. A Carlo Erba NC 2500 elemental analyzer (EA) is used to convert total nitrogen in a solid sample into N 2 gas. The EA is connected to a continuous flow isotope-ratio mass spectrometer (CF-IRMS), which determines relative difference in the isotope-amount ratios of stable nitrogen isotopes ( 15 N/ 14 N)of the product N 2 gas. The combustion is quantitative; no isotopic fractionation is involved. Samples are placed in a tin capsule and loaded into the Costech Zero Blank Autosampler of the EA. Under computer control, samples are dropped into a heated reaction tube that contains an oxidant, where the combustion takes place in a helium atmosphere containing an excess of oxygen gas. Combustion products are transported by a helium carrier through a reduction tube to remove excess oxygen and convert all nitrous oxides into N 2 and through a drying tube to remove water. The gas-phase products, mainly CO 2 and N 2 , are separated by a gas chromatograph. The gas is then introduced into the isotope-ratio mass spectrometer (IRMS) through a Finnigan MAT (now Thermo Scientific) ConFlo II interface, which also is used to inject N 2 reference gas and helium for sample dilution. The IRMS is a Thermo Scientific Delta V Plus CF-IRMS. It has a universal triple collector, two wide cups with a narrow cup in the middle, capable of measuring mass/charge ( m/z ) 28, 29, 30, simultaneously. The ion beams from N 2 are as follows: m/z 28 = N 2 = 14 N 14 N; m/z 29 = N 2 = 14 N 15 N primarily; m/z 30 = NO = 14 N 16 O primarily, which is a sign of contamination or incomplete reduction.

Techniques and Methods↗
USGS2006

Determination of the δ 34 S of Total Sulfur in Solids: RSIL Lab Code 1800

The purpose of Reston Stable Isotope Laboratory Lab (RSIL) Code 1800 is to determine the δ( 34 S/ 32 S), abbreviated as δ 34 S, of total sulfur in a solid sample. A Carlo Erba NC 2500 elemental analyzer (EA) is used to convert total sulfur in a solid sample into SO 2 gas. The EA is connected to a continuous flow isotope-ratio mass spectrometer (CF-IRMS), which determines the relative difference in stable sulfur isotope-amount ratio ( 34 S/ 32 S) of the product SO 2 gas. The combustion is quantitative; no isotopic fractionation is involved. Samples are placed in tin capsules and loaded into a Costech Zero-Blank Autosampler on the EA. Under computer control, samples are dropped into a heated tube reaction tube that combines both the oxidation and the reduction reactions. The combustion takes place in a He atmosphere that contains an excess of oxygen gas at the oxidation zone at the top of the reaction tube. Combustion products are transported by a He carrier through the reduction zone at the bottom of the reaction tube to remove excess oxygen and through a separate drying tube to remove any water. The gas-phase products, mainly CO 2 , N 2 , and SO 2 , are separated by a gas chromatograph (GC). The gas is then introduced into the isotope-ratio mass spectrometer (IRMS) through a Thermo-Finnigan ConFlo II interface, which also is used to inject SO 2 reference gas and He for sample dilution. The IRMS is a Thermo-Finnigan DeltaPlus CF-IRMS. It has a universal triple collector with two wide cups and a narrow cup in the middle. It is capable of measuring mass/charge ( m/z ) 64 and 66 simultaneously. The ion beams from SO 2 are as follows: m/z 64 = SO 2 = 32 S 16 O 16 O; and m/z 66 = SO 2 = 34 S 16 O 16 O primarily.

Techniques and Methods↗
USGS2006

Determination of the δ 34 S of low-concentration sulfate in water; RSIL lab code 1949

The purpose of the Reston Stable Isotope Laboratory (RSIL) lab code 1949 is to determine the δ( 34 S/ 32 S), abbreviated as δ 34 S, of dissolved sulfate having a concentration less than 20 milligrams per liter. Dissolved sulfate is collected on an anion-exchange resin in the field, eluted in the laboratory with 3 M KCl, and precipitated with BaCl 2 at pH 3 to 4 as BaSO 4 . The precipitated BaSO 4 is filtered and dried before introduction into an elemental analyzer (EA) Carlo Erba NC 2500. The EA is used to convert sulfur in a BaSO 4 solid sample into SO 2 gas, and the EA is connected to a continuous flow isotope-ratio mass spectrometer (CF-IRMS), which determines differences in the isotope-amount ratios of stable sulfur isotopes ( 34 S/ 32 S) of the product SO 2 gas. The combustion is quantitative; no isotopic fractionation is involved. Samples are placed in a tin capsule and loaded into the Costech Zero Blank Autosampler of the EA. Under computer control, samples are dropped into a heated reaction tube that combines the oxidation and reduction reactions. The combustion takes place in a helium atmosphere containing an excess of oxygen gas at the oxidation zone at the top of the reaction tube. Combustion products are transported by a helium carrier through the reduction zone at the bottom of the reaction tube to remove excess oxygen and through a separate drying tube to remove any water. The gas-phase products, mainly CO 2 , N 2 , and SO 2 , are separated by a gas chromatograph. The gas is then introduced into the isotope-ratio mass spectrometer (IRMS) through a Finnigan MAT (now Thermo Scientific) ConFlo II interface, which is also used to inject SO 2 reference gas and helium for sample dilution. The IRMS is a Thermo Scientific Delta V Plus CF-IRMS. It has a universal triple collector with two wide cups and a narrow cup in the middle. It is capable of measuring mass/charge ( m/z ) 64 and 66 simultaneously. The ion beams from SO 2 are as follows: m/z 64 = SO 2 = 32 S 16 O 16 O; m/z 66 = SO 2 = 34 S 16 O 16 O primarily.

Techniques and Methods↗
USGS2013

SLAMMER: Seismic LAndslide Movement Modeled using Earthquake Records

This program is designed to facilitate conducting sliding-block analysis (also called permanent-deformation analysis) of slopes in order to estimate slope behavior during earthquakes. The program allows selection from among more than 2,100 strong-motion records from 28 earthquakes and allows users to add their own records to the collection. Any number of earthquake records can be selected using a search interface that selects records based on desired properties. Sliding-block analyses, using any combination of rigid-block (Newmark), decoupled, and fully coupled methods, are then conducted on the selected group of records, and results are compiled in both graphical and tabular form. Simplified methods for conducting each type of analysis are also included.

Techniques and Methods↗
USGS2021

A numerical model for the cooling of a lava sill with heat pipe effects

Understanding the cooling process of volcanic intrusions into wet sediments is a difficult but important problem, given the presence of extremely large temperature gradients and potentially complex water-magma interactions. This report presents a numerical model to study such interactions, including the effect of heat pipes on the cooling of volcanic intrusions. Udell (1985) has shown that heat pipes may develop in heated saturated granular media under laboratory conditions. In previous work, Baker and others (2015) calculated temperatures in the vicinity of a volcanic sill that intruded into wet sediment, showing an unexpected temperature profile in which peak temperatures remained near constant over a region extending a meter above and below the sill. This is challenging to explain with conduction or convection heating methods but is predicted if the heat transfer is performed primarily by a heat pipe. We have numerically modeled the cooling of a lava sill under similar circumstances, using the experimental findings of Udell (1985) to estimate the characteristics of the heat pipe. We have constructed a model using Microsoft C#.NET, complete with an intuitive graphical user interface. The model is available from the U.S. Geological Survey and is capable of being run on Microsoft Windows 7 and higher with modest hardware. We find that the resulting overall temperature profile has some key similarities to the profile inferred by Baker and others (2015). Future models including more detailed convective heat transfer physics will be necessary to fully reproduce the effects of boiling in sediments.

Techniques and Methods↗
USGS2015

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

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

Techniques and Methods↗