USGS ScienceSearch

SEARCH · USGS Science

Results for “Final Project Report”

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.

822 records · Page 37Linked to original sources

Streamflow, base flow, and ground-water recharge in the Housatonic River basin, western Massachusetts and parts of eastern New York and northwestern Connecticut

Streamflows for selected flow durations from 1 to 99 percent and the August median streamflows were estimated for 11 long-term streamflow-gaging stations in and near the study area. Estimates of streamflow and associated standard errors were determined for selected flow durations from 50 to 99 percent and the August median streamflows for 21 low-flow partial-record stations and for selected flow durations from 1 to 99 percent and the August median streamflows for two partial-record stations and seven short-term discontinued streamflow-gaging stations. Median streamflows per square mile for the 10-, 50-, and 90-percent flow durations and the August median streamflows were 3.90, 1.01, 0.185, and 0.248 cubic feet per second per square mile. Streamflows per square mile at selected flow-duration discharges between 1 and 99 percent at the 41 stations were related to basin characteristics to explain differences in streamflow characteristics. Basin characteristics included basin elevations, extent of stratified-drift deposits, land use, aspect, and underlying bedrock geology types. Most streamflow differences were positively correlated to basin elevation differences, most likely because precipitation increases with elevation, and to stratified-drift deposits, which allow more precipitation to recharge the ground water and to discharge later than do till and bedrock deposits. Mean base flow was computed from continuous records of daily mean discharge at 11 long-term streamflow-gaging stations in and near the study area. Mean annual base flow ranged from 13.4 to 24.5 inches per year. Minimum annual base flow ranged from 45 to 72 percent of mean annual rates at the 11 long-term stations, and the ratio of base flow to streamflow (base-flow index) ranged from 0.55 to 0.80. Base-flow durations between 1 and 99 percent were calculated from streamflow records at the 11 long-term streamflow-gaging stations. Base flow accounted for 45.5 to 85.0 percent of total annual streamflow at the 1- and 99-percent flow durations. Ground-water-recharge rates were computed from continuous records of daily mean discharge at 11 long-term streamflow-gaging stations in and near the study area. Mean annual ground-water-recharge rates ranged from 17.5 to 22.4 inches per year at 10 of the 11 long-term stations. Mean annual ground-water-recharge rates ranged from 2 to 7 inches per year higher than base flow. Minimum annual ground-water-recharge rates ranged from 48 to 72 percent of mean annual ground-water-recharge rates. Mean annual potential ground-water recharge was estimated from monthly climatological data collected at six climatological stations in and near the study area. Mean potential ground-water recharge ranged from about 17.9 to 28.9 inches per year, with a median value of 22.6 inches per year. This median value compares well to that calculated by use of streamflow records at the 11 streamflow-gaging stations (20.0 inches per year). Streamflows per square mile for the 10-, 50-, and 90-percent flow durations at stations in and near the study area were similar to those computed for other unregulated long-term continuous streamflow-gaging stations in central and eastern Massachusetts. Base-flow and ground-water-recharge rates in the study area compared closely to results from other studies in southeastern Massachusetts and Rhode Island, which were based on the same computational methods.

Connecticut, Massachusetts, New York

Water resources of Mellette and Todd counties, South Dakota

Mellette and Todd Counties are located in south-central South Dakota and have a combined area of 2,694 square miles. The White River and its tributaries, which include the Little White River, drain Mellette County and about one-half of Todd County. Tributaries to the Niobrara River, which include the Keya Paha River, drain the other one-half of Todd County. The average discharge of the Little White River is about 56 cubic feet per second as the river enters Todd County and is about 131 cubic feet per second as it discharges to the White River in northern Mellette County. The average discharge of the Keya Paha River just outside Todd County is about 39 cubic feet per second. The average annual runoff for Mellette and Todd Counties ranges from 0.94 to 2.36 inches based on records from nine streamflow-gaging stations in and near the counties. The average annual runoff is 1.62 inches, which compares with the average annual precipitation of about 19 inches. In Todd County, shallow wells completed in the alluvial, Ogallala, Arikaree, and White River aquifers generally can supply water that has low concentrations of dissolved solids, is fresh, and is soft to moderately hard. Ground water from shallow aquifers is limited in Mellette County; therefore, deep wells, often greater than 1,000 feet, are sometimes installed. The Pierre Shale often is used to supply rural domestic and stock wells in Mellette County even though well yields are low and the water has high dissolved solids, is moderately saline, and is very hard. Alluvial aquifers are present in both counties and store an estimated 1.6 million acre-feet of water. The water quality of the alluvial aquifers is dependent on the underlying deposits, and generally the water has low concentrations of dissolved solids, is fresh, and is soft to moderately hard where underlain by the Ogallala and Arikaree Formations; has moderate concentrations of dissolved solids, is slightly saline, and is hard where underlain by the White River Group; and has high concentrations of dissolved solids, is saline, and is very hard where underlain by the Pierre Shale. Also, yields often are lower where the alluvial aquifers are underlain by the Pierre Shale. The Ogallala aquifer is present in only Todd County, and the Arikaree aquifer is present throughout most of Todd County and southwestern and south-central Mellette County. The Ogallala aquifer contains an estimated 17 million acre-feet of water in storage, and the Arikaree aquifer contains an estimated 50 million acre-feet of water in storage. Both aquifers generally are suitable for irrigation, and yields from these aquifers are sometimes greater than 1,000 gallons per minute. Nitrate concentrations in 13 out of 92 water samples collected from the Ogallala aquifer exceeded the Primary Drinking Water Maximum Contaminant Level (MCL) of 10 milligrams per liter. In 11 out of 46 samples collected from the Arikaree aquifer, arsenic concentrations exceeded the MCL of 50 micrograms per liter. The White River aquifer, where present, is usually the shallowest source of ground water in Mellette County. The White River aquifer also is used in northern Todd County where the Ogallala and Arikaree aquifers are not present. The White River aquifer contains an estimated 50 million acre-feet of water in storage. Reported yields from the aquifer range from 1 to 30 gallons per minute, which generally is insufficient to support irrigation in most areas. However, yields are sufficient for livestock-watering and rural-domestic purposes. In both counties, the Pierre Shale is the shallowest bedrock aquifer and is exposed at the land surface throughout most of Mellette County. This aquifer is used primarily in Mellette County. Although the aquifer contains an estimated maximum of 1.5 million acre-feet of water in storage, it is not a viable source of ground water because the aquifer is relatively impermeable, yields are low, and water usually can be obtained from shallower sources, especially in Todd County. Reported yields from the Pierre Shale aquifer range from 1 to 8 gallons per minute. Because few test holes and wells penetrate below the Pierre Shale, little is known about the extent of the deeper bedrock aquifers. All wells completed in the Dakota Sandstone, Inyan Kara, and Minnelusa and Madison aquifers in the counties are used for stock-watering purposes. High concentrations of dissolved solids and hard water are characteristic of the water quality in the bedrock aquifers. Depths to the top of the deeper bedrock aquifers range from 1,270 feet to greater than 2,000 feet below land surface.

South Dakota

Revision of ModelMuse to support the use of PEST software with MODFLOW and SUTRA models

Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.

Techniques and Methods

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota

Per- and polyfluoroalkyl substances in waters associated with oil and gas development in the Denver Basin

Use of per- and polyfluoroalkyl substances (PFAS) in the petroleum industry could be a cause for concern due to the large volumes of produced water (PW) generated during oil and gas extraction, the reuse of these wastes in water-stressed regions, and adverse health outcomes related to PFAS exposures. However, PW PFAS characterization is nearly absent in the literature, and hydraulic fracturing (HF) chemical disclosures often omit the identities of additives as proprietary. Here we evaluate PFAS in PW samples from three petroleum wells in the Denver Basin during their first year of production. Total concentrations of targeted PFAS (Σ 40 PFAS) were < 35 ng/L in PW samples, with short-chain PFAS like perfluorobutanoic acid persisting throughout the sampled duration. Analysis of freshwater inputs for hydraulic fracturing (Σ 40 PFAS ~ 113 ng/L) and mixed fracture fluid (Σ 40 PFAS ~ 69 ng/L) indicated much of the targeted PFAS content was derived from the input water, and not from HF additives, however samples subjected to oxidation indicated the presence of PFAS precursors that would not be detected by targeted analysis. This study highlights that while PFAS content is low in the studied PWs, the potential for redistribution of PFAS in the environment may be a consideration for reuse applications.

Colorado, Nebraska, Wyoming

Surficial geologic map of the Owlshead Mountains 30' x 60' quadrangle, Inyo and San Bernardino Counties, California

The surficial geologic map of the Owlshead Mountains 30' x 60' quadrangle depicts the distribution and characteristics of surficial-deposit materials and neotectonic deformation for an area of approximately 5,000 square kilometers (km 2 ) located in the western Basin and Range Province of eastern California. The map represents a new compilation of the surficial geology that encompasses deposits within the late Pliocene to Quaternary. The map is based primarily on new mapping conducted between 2001 and 2009. Map compilation was supported by field observations distributed across the map area, combined with reference to several published and unpublished mapping sources that mostly emphasized neotectonic deformation. The surficial-deposit units included in the map follow a classification scheme that systematically denotes depositional process, relative age, and any secondary sedimentologic or morphologic characteristics. Identification, correlation, and age estimation of map units are based primarily on the relative degree of development of certain time-dependent characteristics such as surface morphology, including local dissection and surface preservation, surface clast modification, and degree of soil development; these characteristics are implicitly incorporated into unit designations. The map represents a detailed and regionally uniform synthesis of the late Neogene geology for this large area that provides a framework applicable to many interpretative studies, such as regional patterns of deposition and dissection; surface drainage development and evolution; and the distribution, style, and timing of neotectonic deformation.

California

Episodic sediment-discharge events in Cascade Springs, southern Black Hills, South Dakota

Cascade Springs is a group of artesian springs in the southern Black Hills, South Dakota, with collective flow of about 19.6 cubic feet per second. Beginning on February 28, 1992, a large discharge of red suspended sediment was observed from two of the six known discharge points. Similar events during 1906-07 and 1969 were documented by local residents and newspaper accounts. Mineralogic and grain-size analyses were performed to identify probable subsurface sources of the sediment. Geochemical modeling was performed to evaluate the geochemical evolution of water discharged from Cascade Springs. Interpretations of results provide a perspective on the role of artesian springs in the regional geohydrologic framework. X-ray diffraction mineralogic analyses of the clay fraction of the suspended sediment were compared to analyses of clay-fraction samples taken from nine geologic units at and stratigraphically below the spring-discharge points. Ongoing development of a subsurface breccia pipe(s) in the upper Minnelusa Formation and/or Opeche Shale was identified as a likely source of the suspended sediment; thus, exposed breccia pipes in lower Hell Canyon were examined. Upper Minnelusa Formation breccia pipes in lower Hell Canyon occur in clusters similar to the discrete discharge points of Cascade Springs. Grain-size analyses showed that breccia masses lack clay fractions and have coarser distributions than the wall rocks, which indicates that the red, fine-grained fractions have been carried out as suspended sediment. These findings support the hypothesis that many breccia pipes were formed as throats of abandoned artesian springs. Geochemical modeling was used to test whether geochemical evolution of ground water is consistent with this hypothesis. The evolution of water at Cascade Springs could not be suitably simulated using only upgradient water from the Minnelusa aquifer. A suitable model involved dissolution of anhydrite accompanied by dedolomitization in the upper Minnelusa Formation, which is caused by upward leakage of relatively fresh water from the Madison aquifer. The anhydrite dissolution and dedolomitization account for the net removal of minerals that would lead to breccia pipe formation by gravitational collapse. Breccia pipes in the lower Minnelusa Formation are uncommon; however, networks of interconnected breccia layers and breccia dikes are common. These networks, along with vertical fractures and faults, are likely pathways for transmitting upward leakage from the Madison aquifer. It is concluded that suspended sediment discharged at Cascade Springs probably results from episodic collapse brecciation that is caused by subsurface dissolution of anhydrite beds and cements of the upper Minnelusa Formation, accompanied by replacement of dolomite by calcite. It is further concluded that many breccia pipes probably are the throats of artesian springs that have been abandoned and exposed by erosion. The locations of artesian spring-discharge points probably have been shifting outwards from the center of the Black Hills uplift, essentially keeping pace with regional erosion over geologic time. Thus, artesian springflow probably is a factor in controlling water levels in the Madison and Minnelusa aquifers, with hydraulic head declining over geologic time, in response to development of new discharge points. Development of breccia pipes as throats of artesian springs would greatly enhance vertical hydraulic conductivity in the immediate vicinity of spring-discharge points. Horizontal hydraulic conductivity in the Minnelusa Formation also may be enhanced by dissolution processes related to upward leakage from the Madison aquifer. Potential processes could include dissolution resulting from leakage in the vicinity of breccia pipes that are abandoned spring throats, active spring discharge, development of subsurface breccias with no visible surface expression or spring discharge, as well as general areal leakage from the Madison aquifer into the Minnelusa Formation.

South Dakota

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri and Mississippi Rivers on the periphery of Missouri, June 13–22, 2022

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near seven bridges at six highway crossings of the Missouri and Mississippi Rivers on the periphery of Missouri from June 13–22, 2022. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 feet longitudinally and generally extending laterally across the active channel from bank to bank during minor flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a minor flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All bridge sites in this study were surveyed and documented in previous studies. Although partial exposure of substructural support elements was observed at several piers, at most sites the exposure most likely is minimal compared to the overall substructure that remains buried in bed material at these piers. The notable exceptions are piers 12 and 13 at structure L0135 on State Highway 51 at Chester, Illinois, where the bedrock material was fully exposed around the piers. The average difference between the bathymetric surfaces between 2022 and 2018 varied from 0.41 foot higher to 1.86 feet lower. Between 2022 and 2014, the average difference between the bathymetric surfaces varied from 1.02 feet higher to 4.69 feet lower. Only the two sites on the Missouri River and the Caruthersville site were surveyed in 2011; for those sites, the average difference between the bathymetric surfaces varied from 5.83 feet higher to 1.34 feet lower. The most substantial overall net gain of sediment in a reach was between 2011 and 2022 at structure A1700 near Caruthersville, Mo. (site 38). This result was expected because structure A1700 is downstream from the confluences of the Missouri and Ohio Rivers, and therefore subject to the largest streamflows, the largest streamflow fluctuations, and the most substantial sediment flux, as has historically been observed at this site. The presence of riprap blankets, pier size and nose shape, and alignment to flow had a substantial effect on the size of the scour hole observed for a given pier. Piers that were surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers at which no rock or riprap were present. New riprap blankets were surveyed at pier 3 of structure L0098 at Brownville, Nebraska, and at piers 15–18 of structure A1700 near Caruthersville, Mo., that effectively mitigated the scour holes historically observed at these piers. Narrow piers having round or sharp noses that were aligned with flow often had scour holes that were difficult to discern from nearby bed features, whereas piers having wide or blunt noses resulted in larger, deeper scour holes. Several of the structures had piers that were skewed to primary approach flow. Scour holes near these piers consistently displayed greater depth on the side of the pier with impinging flow and deposition on the leeward side of the pier.

Missouri

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida

New Hampshire and Landsat

At its widest point, a mere 80 miles separate the eastern and western borders of New Hampshire. Its northern and southern borders are just 175 miles apart. Even so, few States can boast as much rugged natural beauty per mile as the Nation’s fifth smallest. Nestled within New Hampshire are 93 State parks teeming with moose, Ursus americanus (Pallas, 1780; black bears), coyotes, beavers, river otters, and foxes. The largest section of White Mountain National Forest cuts across north-central New Hampshire, drawing visitors to its lakes, streams, mountain peaks, and hardwood forests. New Hampshire also is home to Lake Winnipesaukee, the State’s largest lake, notable for its floating post offices, the annual “ice-out” contest that sees residents vying to guess the date its surface ice dissipates, and its supporting role in films such as “On Golden Pond” and “What About Bob?” However, the scenic forests of New Hampshire face challenges in the form of invasive species such as Lymantria dispar (Linnaeus, 1758; spongy moth), Adelges piceae (balsam woolly adelgid), and Agrilus planipennis (emerald ash borer). In recent years, New Hampshire’s lakes and streams have seen more cyanobacterial blooms as well. The U.S. Geological Survey Landsat Program offers a consistent, reliable, and historically unmatched source of Earth observations that can aid in the mapping, monitoring, and management of New Hampshire’s land and water resources. Here are a few ways Landsat data have been used in the Granite State.

New Hampshire

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, May 19–26, 2021

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near nine bridges at eight highway crossings of the Missouri River between Kansas City and St. Louis, Missouri, from May 19 to 26, 2021. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,840 feet (ft) longitudinally and generally extending laterally across the active channel from bank to bank during low to moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All the bridge sites in this study were previously surveyed and documented in previous studies. Comparisons between bathymetric surfaces from the previous surveys and those of the current (2021) study do not indicate any consistent correlation between channel-bed elevations and streamflow conditions. The average difference between the bathymetric surfaces varied from 1.59 ft higher to 0.95 ft lower in 2021 than 2017, which corresponds to a gain of 100,200 cubic yards and a loss of 55,800 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 2.74 ft higher to 3.05 ft lower in 2021 than 2013, which corresponds to a gain of 111,500 cubic yards and a loss of 169,200 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 4.52 ft higher to 1.38 ft lower in 2021 than 2011, which corresponds to a gain of 221,100 cubic yards and a loss of 90,300 cubic yards, respectively. The most substantial overall net gain was 221,100 cubic yards between 2011 and 2021 at structures L0550 and A4497 at Jefferson City (site 20). The large net gain likely results from a combination of the mitigation of the scour holes near pier 4 of both bridges and the substantially lower flow in 2021 than in 2011. Alternatively, the most substantial overall net loss was 169,200 cubic yards between 2013 and 2021 at structure A6288 at Hermann (site 21), despite comparable streamflows. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; furthermore, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the side with impinging flow. At structure A6288 at Hermann (site 21), the scour hole near pier 5 was difficult to discern from nearby dunes and ripples, whereas the upstream edge of the footing was visible at pier 4, which likely contributes to the larger scour hole near that pier; the top of the footing may blunt the horseshoe vortex at pier 5, but the exposed front of the footing may exacerbate the vortex at pier 4.

Missouri

Leachate plumes in ground water from Babylon and Islip landfills, Long Island, New York

Landfills operated by the Towns of Babylon and Islip in southwest and central Suffolk County contain urban refuse, incinerated garbage, and scavenger (cesspool) waste; some industrial refuse is deposited at the Babylon site. The Islip landfill was started in 1933, the Babylon landfill in 1947. The landfills are in contact with and discharge leachate into the highly permeable upper glacial aquifer (hydraulic conducitvity 190 to 500 ft/d [feet per day]). The aquifer is 74 feet thick at the Babylon landfill and 170 feet thick at the Islip landfill. The leachate-enriched water occupies the entire thickness of the aquifer beneath both landfills, but hydrologic boundaries retard downward migration of the plumes to deeper aquifers. The Babylon plume is 1,900 feet wide at the landfill and narrows to about 700 feet near its terminus 10,000 feet from the landfill. The Islip plume is 1,400 feet wide at the landfill and narrows to 500 feet near its terminus 5,000 feet from the landfill. Hydrochemical maps and sections show the distribution of the major chemical constituents of the plumes. The most highly leachate-enriched ground water obtained was from the Babylon site; it contained 860 mg/L (milligrams per liter) sodium, 110 mg/L potassium, 565 mg/L calcium, 100 mg/L magnesium, 2,700 mg/L bicarbonate, and 1,300 mg/L chloride. Sulfate was notably absent or in low concentration in most parts of both plumes. Nitrogen in plume water was mostly in the form NH4+, and concentrations as high as 90 mg/L were found; concentrations of nitrogen as N in the plume were less than 10 mg/L. As much as 440 mg/L iron and 190 mg/L manganese were found in the leachate-enriched water. Samples were also tested for arsenic, boron, cadmium, cobalt, chromium, copper, mercury, nickel, lead, selenium, strontium, and zinc. Boron was more or less ubiquitous and was found in concentrations as high as 2 mg/L. Organic carbon was found in concentrations as high as 2,250 mg/L in the most highly leachate-enriched water but attenuated rapidly to less than 20 mg/L. Dissolved-solids concentrations near the landfills were between 400 and 3,000 mg/L at Babylon and between 500 and 1,500 mg/L at Islip. Ground-water temperatures near the landfills exceed those in ambient water by as much as 7°C at Babylon and 16°C at Islip. Heat contributed by the landfills was mostly dissipated within 0.4 mi of the landfill, butat Islip, the warm leachate-enriched water extended 0.5 mi downgradient. The entrance of leachate into the less dense ground water as pulsations after rainfall may explain the presence of high leachate enrichment near the bottom of the aquifer. A comparison of the physical characteristics of leachate-enriched ground water with those of ambient water suggests that the downward movement of leachate results from its greater density. Simulation of the movement and dispersion of the Babylon plume with a mathematical dispersion model indicated the coefficient of longitudinal dispersion to be about 60 ft 2 /d (feet squared per day) and the ground-water velocity to be 1 ft/d. However, the velocity determined from the hydraulic gradient and public-supply wells in the area was 4 ft/d, which would cause a plume four times as long as that predicted by the mathmatical dispersion model. At the Islip site, the plume was one-third the length calculated on the basis of the age of the landfill. The shortness of the plumes has not been explained; it may be a result of the leachate's having been too dilute to form a plume during the early years of the landfills.

New York

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2022

Potentiometric-surface and depth-to-water maps for spring 2022 were created for the Mississippi River Valley alluvial aquifer (MRVA) using groundwater-altitude data from 1,136 wells completed in the MRVA and from the altitude of the top of the water surface in area rivers from 160 streamgages. The potentiometric-surface and depth-to-water maps for 2022 were created to support investigations to characterize the MRVA as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the MRVA for spring 2022 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 340 feet (ft) above the North American Vertical Datum of 1988. The regional direction of groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient was into the depression, and near rivers, where the groundwater gradient can be from aquifer to the river or from the river into the aquifer. There are large depressions in the potentiometric-surface map in the lower one-half of the Cache region and in much of the Grand Prairie and Delta regions. Depth to water in the MRVA, spring 2022, by well ranged from 5.00 ft above land surface to 145.66 ft below land surface.

Mississippi River Valley alluvial aquifer

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2019–2021

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Water availability is an important issue in the Black Mesa area because of the arid climate, past industrial water use, and continued water requirements for municipal use by a growing population. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year, depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents the results of data collected as part of the monitoring program in the Black Mesa area from calendar years 2020–2021 and, additionally, uses streamflow statistics from November and December 2019. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2020, total groundwater withdrawals were estimated to be 2,680 acre-feet (acre-ft), and, in 2021, total withdrawals were estimated to be 2,570 acre-ft. Total withdrawals during 2021 were about 65 percent less than total withdrawals in 2005 because the Peabody Western Coal Company discontinued its use of water to transport coal in a coal slurry pipeline after 2005 and ceased mining operations in 2019. Owing to Navajo Nation and Hopi Reservation access restrictions during the Coronavirus pandemic, water levels were not collected from municipal wells in 2020 or 2021. Water levels measured in 2021 from wells completed in the unconfined areas of the N aquifer within the Black Mesa area showed a decline in 7 of 13 wells when compared with water levels from the prestress period (prior to 1965). The changes in water levels across all 13 wells ranged from +8.4 feet (ft) to −42.4 ft, and the median change was −0.4 ft. Water levels also showed decline in 11 of 12 wells measured in the confined area of the aquifer when compared to the prestress period. The median change for the confined area of the aquifer was −25.9 ft, with changes across all 12 wells ranging from +17.3 ft to −133.7 ft. Spring flow was measured at four springs between 2020 and 2021. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso, Arizona. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s, and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976–2021), Dinnebito Wash near Sand Springs 09401110 (1993–2020), Polacca Wash near Second Mesa 09400568 (1994–2020), and Pasture Canyon Springs 09401265 (2004–2021). Median winter flows (November through February) of each winter were used as an estimate of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Polacca Wash and Pasture Canyon Springs, whereas a decreasing trend was observed at Moenkopi Wash and Dinnebito Wash. In 2020 and 2021, water samples were collected from a total of four springs in the Black Mesa area and analyzed for selected chemical constituents. Results from the four springs were compared with previous analyses from the same springs. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. However, concentrations of chloride from Pasture Canyon Spring show a diminishing trend. Concentrations of dissolved solids, chloride, and sulfate at Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Concentrations of dissolved solids at Burro Spring have varied for the period of record, but there is no statistical trend in the data. However, concentrations of chloride and sulfate from Burro Spring show a trend towards lower concentrations.

Arizona

Determination of per- and polyfluoroalkyl substances in water by direct injection of matrix-modified centrifuge supernatant and liquid chromatography/tandem mass spectrometry with isotope dilution

A direct-injection liquid chromatography/tandem mass spectrometry method was developed to determine 34 per- and polyfluoroalkyl substances (PFAS), including selected branched isomers, in centrifuge supernatant of matrix-modified (amended with approximately 50 percent methanol) water samples. The method has been validated in reagent water, surface water, groundwater, and wastewater effluent. Other water types (for example, drinking water, untreated wastewater, and landfill leachate) have been analyzed by the method but not systematically validated. Recovery of isotope-dilution standards, added to each sample, may be used to assess method performance in nonvalidated matrices on a sample-by-sample basis. Using this method, PFAS concentrations were determined in the range of 2–2,000 nanograms per liter in water samples. This range can be extended by diluting concentrated samples. At circumneutral pH, most compounds are present in the environment in their ionized form, and data are reported as such (for example, perfluorooctanoic acid is referred to as “perfluorooctanoate” [PFOA], perfluorooctane sulfonic acid is referred to as “perfluorooctane sulfonate” [PFOS]). Sample preparation procedures were designed without the use of filtration and with minimum sample handling steps to mitigate procedural losses of target compounds due to sorption to surfaces. Further, isotope-dilution quantification allowed for the correction of bias that may result from procedural losses, matrix-induced signal suppression or enhancement, and other factors. Validation experiments to characterize bias and variability, method detection level, and holding time were done in four distinct water matrices—reagent water, surface water, treated wastewater effluent, and groundwater—at multiple concentration levels. Mean PFAS recoveries met data quality objectives of bias and variability studies in all four validation matrices except for two compounds with low and variable recovery in the reagent water matrix only. Isotope-dilution standards, treated as surrogate compounds, were analyzed in more than 1,500 customer-submitted environmental samples with aggregate recovery of 102.5±6.5 percent (mean±standard deviation). Maximum holding times for all target compounds in the four validation matrices were 28 days for refrigerated samples and 90 days for frozen samples.

Techniques and Methods