USGS Science⌕ Search

SEARCH · USGS Science

Results for “New Zealand Journal of Geology and Geophysics”

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.

1,689 records · Page 39Linked to original sources
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2025

Spatio-temporal evolution of distributed volcanic fields, case studies—Sierra Chichinautzin and Michoacán-Guanajuato, México

An analysis of 1,375 volcanoes in the Michoacán-Guanajuato (1,148 volcanoes in a 26,200 square-kilometer area) and Sierra Chichinautzin (227 volcanoes in a 3,500 square-kilometer area) volcanic fields in central Mexico identified patterns in the spatial and temporal distribution of past eruptions. A cluster agglomerative hierarchical method and kernel analysis confirmed that the Michoacán-Guanajuato volcanic field comprises four volcanic fields (Valle de Santiago, Uruapan, Apatzingán, and Pátzcuaro volcanic fields) controlled by different fault systems, indicating that it is not a single volcanic field but rather a group of volcanic fields (a “superfield”), each of which has distinct characteristics. In the Sierra Chichinautzin volcanic field, well-constrained isotopic ages were used to build a model of how the spatial distribution of the eruptions has changed over time. Two new 40Ar/39Ar ages from a locally recognized volcanic feature near the town of El Cantil, herein called El Cantil volcano (1,537±17 kilo-annum [ka]) and the volcanic feature at Cerro el Elefante (herein called El Elefante dome) (1,485±92 ka) belong to the oldest volcanic group identified in the Sierra Chichinautzin volcanic field, confirming the timing of the beginning of monogenetic volcanism in the region. Based on the volcanic groups identified in the Sierra Chichinautzin volcanic field, the youngest volcanism (less than 35 ka) is found only in the central-western sector of the field. Principal component analysis determined the directional trends of feeder dikes only for vents <10 ka in the Sierra Chichinautzin volcanic field. Possible magma migration paths through the crust were identified using seismic data from both volcanic fields using an earthquake catalog from 1973 to 2023, which includes 9,016 earthquakes in the Michoacán-Guanajuato volcanic field and 841 in the Sierra Chichinautzin volcanic field. The spatial distribution of the hypocenters does not highlight any trend that could be associated with superficial movement of magma in the Sierra Chichinautzin volcanic field. In the Michoacán-Guanajuato volcanic field, however, eight seismic swarms since 1997 have been detected. These swarms are interpreted to result from ascending magma. Strengthening monitoring systems and reinforcing mitigation measures to address volcanic hazards and risk are important means of preparing for future eruptions in both regions. Analysis such as those herein provide insights as to where an eruption might occur and may help mitigate volcanic hazards.

Michoacán-Guanajuato, Sierra Chichinautzin↗
USGS2025

Map of topographic lineaments interpreted as recent surface ruptures along the Bennett Valley and Southern Maacama Fault Zones, Sonoma County, California

This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.

California↗
USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗
USGS2026

Complex carbonate ore mineralogy in the Mountain Pass carbonatite rare earth element deposit, USA

Economic concentrations of rare earth element (REE) minerals are uncommon in the Earth’s crust, with most occurring in carbonatites. Unlike most igneous rocks composed of silicate minerals, carbonatites are dominated by carbonate minerals, some of which can incorporate significant light REEs (LREEs; La, Ce, Pr, Nd). Technological applications of REEs are numerous and they have been identified as some of the most critical mineral commodities to the global economy. The Mountain Pass carbonatite stock in the Mojave Desert of California is the most economically significant REE deposit in the USA and contains a few to tens of percent (by volume) of the carbonate REE ore mineral bastnäsite. Despite the economic significance of the Mountain Pass deposit, studies of its ore mineralogy are limited. Here we present new carbonate ore mineralogy data for a compositionally diverse suite of carbonatitic rocks from the Mountain Pass stock and related dikes. Whole-rock geochemical data are integrated with mineral-scale textural and chemical data obtained by scanning electron microscopy (SEM), electron probe microanalysis (EPMA), and microRaman spectroscopy. Our results document a complex spectrum of REE-bearing carbonate minerals and intermediate mixed-layer structures. Mineral species include bastnäsite [REE(CO 3 )F], hydroxylbastnäsite [REE(CO 3 )OH], parisite [Ca(REE) 2 (CO 3 ) 3 F 2 ], synchysite [Ca(REE)(CO 3 ) 2 F], röntgenite [Ca 2 (Ce,La) 3 (CO 3 ) 5 F 3 ], and sahamalite [(Mg,Fe 2+ )(REE) 2 (CO 3 ) 4 ]. Carbonate ore mineralogy is heterogeneous within and between samples, including at the intracrystal scale. Complexly zoned crystals exhibit as many as five to six different compositional domains and syntaxial intergrowths, commonly with the more Ca-rich varieties (parisite, synchysite) forming crystal rims that surround relict bastnäsite cores. We attribute the phenocryst variability to changes in the chemistry and temperature of primary carbonatite magmas and evolved/exsolved fluids. Cross-cutting vein textures of calcite, celestine and various REE carbonate minerals, interstitial bastnäsite crystallization, breccia blocks lined by fine-grained bastnäsite, and the presence of hydroxylbastnäsite and partially hydroxylated bastnäsite point to the role of secondary hydrothermal processes in REE mineralization. Fluorcarbonate mineral compositions demonstrate that La and Ce are more structurally abundant in bastnäsite, whereas the more Ca-rich species (parisite, synchysite) contain a greater proportion of REE heavier than Pr (Nd, Sm, Eu, Gd) and Y. Atomic ratios of Pr/(Nd + Pr) are likewise variable, with the highest average value for bastnäsite (0.25) compared to parisite (0.22) and sychysite (0.21). This finding has geometallurgical implications, given that current mining operations are focused on recovery of Nd and Pr for high field strength permanent magnets and the Nd/Pr ratios are a critical factor in ore processing and magnet manufacture.

California↗
USGS2025

Hydrogeologic framework of the Mountain Home area, southern Idaho

In the arid western Snake River Plain around the City of Mountain Home, Idaho, declining groundwater levels concern agricultural, municipal, and other water users who rely on groundwater for sustenance because surface-water resources are limited. The U.S. Geological Survey developed this hydrogeologic framework to provide an updated characterization of groundwater resources in the western Snake River Plain around the City of Mountain Home. The hydrogeologic framework comprises: (1) a conceptual description of hydrogeologic units, (2) a three-dimensional hydrogeologic model and borehole database, (3) a map of groundwater levels and change, and (4) a discussion of groundwater occurrence and movement within the study area. Hydrogeologic units were defined based on existing literature and the borehole database compiled for this study; the five hydrogeologic units are granite, rhyolite, basalt, fine-grained sediments, and coarse-grained sediments. Each unit can bear water, but the main regional aquifer in the study area occurs in the basalt and fine-grained sediment units with depth to water ranging from 150 to 765 feet. A perched groundwater zone near the City of Mountain Home is primarily hosted in basalt and used domestically with most depths to water ranging from 30 to 100 feet. Interflow zones, scoria, and vertical fractures create heterogeneity within the basalt hydrogeologic unit that exerts strong control on groundwater movement, creating horizontal perching conditions and zones of enhanced vertical conductivity that facilitate downward groundwater percolation. In the fine- and coarse-grained sediments and rhyolite units, inferred faults both impede and enhance groundwater movement. The borehole database was constructed by digitizing 540 well-driller reports and was used to build a three-dimensional hydrogeologic framework model which reasonably represents the spatial distribution of hydrogeologic units in the study area. Generally, fine-grained sediments underlie much of the study area, with basalt concentrated in the central and western study area and rhyolite and granite in the uplands to the north. Groundwater levels were measured in 180 wells in March and November 2023; these data were used to develop water-table contour maps and describe groundwater-level change over an irrigation season. Groundwater generally flows south-southwest to the Snake River and groundwater levels declined across most of the study area (from 0.03 to 22.01 feet) between spring and autumn 2023, which is consistent with long-term declines in the Cinder Cone Butte Critical Groundwater Area and Mountain Home Groundwater Management Area. Groundwater levels rose (0.6 to 15.44 feet) over the irrigation season in most wells in the perched groundwater zone near the City of Mountain Home and near the Snake River, indicating the importance of surface-water recharge to groundwater in areas where surface water irrigation occurs. In aggregate, this hydrogeologic framework provides an updated characterization of and new insights into groundwater resources in the study area to help inform water resources management.

Idaho↗
USGS2025

Preliminary field report of landslide hazards following Hurricane Helene

Executive Summary This report reflects our knowledge regarding the widespread landslide activity associated with Hurricane Helene observed during the U.S. Geological Survey’s (USGS) mission assignment to North Carolina in October 2024. The material in this report was originally prepared for the Federal Emergency Management Agency under mission assignment DR-4827-NC. The data and commentary in this report are reflective of a report provided to the Federal Emergency Management Agency (FEMA) on October 18, 2024, as well as information provided in briefings at the Buncombe County Emergency Operations Center. The report has been modified for public dissemination. This assessment was based on systematic visual examination and mapping of landslide locations from aerial and satellite imagery, visual and photographic observations from low-level helicopter overflights and conversations with local landslide experts from the North Carolina Geological Survey and Appalachian Landslide Consultants PLLC, and more than 50 years of combined landslide hazard professional experience of the mission-assigned field team. No systematic field investigations were done by the USGS. While responding to the event, the USGS did not identify any landslides that posed an immediate major threat to recovery personnel in parts of nine counties in North Carolina (Avery, Buncombe, Henderson, McDowell, Mitchell, Polk, Rutherford, Watauga, and Yancey); however, threats from renewed landslide activity may remain heightened in localized areas for months or even years. Known areas of the most abundant landslide occurrence include Bat Cave, Lake Lure, Chimney Rock, Swannanoa, Black Mountain, Fairview, steep areas in Asheville, and the Blue Ridge Parkway. The USGS shared detailed locations of known landslides with the Emergency Operations Centers. The thousands of landslide scars on hillsides and landslide deposits on flatter ground may present some threat to recovery activities. Soil and rocks will continue to erode from newly exposed landslide scars and may pose a threat to people and infrastructure who are immediately nearby. In general, the steeper and taller the landslide scar, the greater the potential threat. This threat is heightened during periods of rainfall and increases with the duration and intensity of rainstorms. Very heavy rainfall, or repeated rainfall events during short periods, could also initiate new landslides on steep slopes. Excavation of landslide deposits, particularly excavation of those deposits directly adjacent to steep slopes, may also pose a threat to nearby people and equipment. An interagency collaborative mapping effort led by the USGS that informed this assessment identified 1,155 landslide locations by the October 2024 briefings, but that number increased to 2,217 in a final reviewed version of the locations published in January 2025. Locations were mapped from satellite imagery, fixed-wing and helicopter surveys, media and social media, and field reports in the 3 weeks following the passage of the remnants of Hurricane Helene. USGS products outlined in this report are publicly available and include geotagged photographs from aerial reconnaissance, hazard models, an interactive view of mapped landslide locations, and landslide safety and education resources.

North Carolina, South Carolina, Tennessee, Virgini↗
USGS1998

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Naval Station Mayport, Florida

Ground-water flow through the surficial aquifer system at Naval Station Mayport near Jacksonville, Florida, was simulated with a two-layer finite-difference model as part of an investigation conducted by the U.S. Geological Survey. The model was calibrated to 229 water-level measurements from 181 wells during three synoptic surveys (July 17, 1995; July 31, 1996; and October 24, 1996). A quantifiable understanding of ground-water flow through the surficial aquifer was needed to evaluate remedial-action alternatives under consideration by the Naval Station Mayport to control the possible movement of contaminants from sites on the station. Multi-well aquifer tests, single-well tests, and slug tests were conducted to estimate the hydraulic properties of the surficial aquifer system, which was divided into three geohydrologic units�an S-zone and an I-zone separated by a marsh-muck confining unit. The recharge rate was estimated to range from 4 to 15 inches per year (95 percent confidence limits), based on a chloride-ratio method. Most of the simulations following model calibration were based on a recharge rate of 8 inches per year to unirrigated pervious areas. The advective displacement of saline pore water during the last 200 years was simulated using a particle-tracking routine, MODPATH, applied to calibrated steady-state and transient models of the Mayport peninsula. The surficial aquifer system at Naval Station Mayport has been modified greatly by natural and anthropogenic forces so that the freshwater flow system is expanding and saltwater is being flushed from the system. A new MODFLOW package (VAR1) was written to simulate the temporal variation of hydraulic properties caused by construction activities at Naval Station Mayport. The transiently simulated saltwater distribution after 200 years of displacement described the chloride distribution in the I-zone (determined from measurements made during 1993 and 1996) better than the steady-state simulation. The advective movement of contaminants from selected sites within the solid waste management units to discharge points was simulated using MODPATH. Most of the particles were discharged to the nearest surface-water feature after traveling less than 1,000 feet in the ground-water system. Most areas within 1,000 feet of a surface-water feature or storm sewer had traveltimes of less than 50 years, based on an effective porosity of 40 percent. Contributing areas, traveltimes, and pathlines were identified for 224 wells at Naval Station Mayport under steady-state and transient conditions by back-tracking a particle from the midpoint of the wetted screen of each well. Traveltimes to contributing areas that ranged between 15 and 50 years, estimated by the steady-state model, differed most from the transient traveltime estimates. Estimates of traveltimes and pathlines based on steady-state model results typically were 10 to 20 years more and about twice as long as corresponding estimates from the transient model. The models differed because the steady-state model simulated 1996 conditions when Naval Station Mayport had more impervious surfaces than at any earlier time. The expansion of the impervious surfaces increased the average distance between contributing areas and observation wells.

Florida↗
USGS2025

Foundational uncertainties in terminal Ediacaran chronostratigraphy revealed by high-precision zircon U-Pb geochronology of the Nama Group, Namibia

The Nama Group of southern Namibia and northwestern South Africa hosts the best-dated mixed carbonate-siliciclastic foreland basin succession of the terminal Ediacaran [ca. 551 million years (Ma) ago to <538 Ma] and is key for resolving the chronology of early metazoan evolution. Numerous silicified volcanic tuff interbeds are present, but differing interpretations regarding the fidelity of their ages lead to different regional stratigraphic correlations, especially for the Urusis Formation of the Schwarzrand Subgroup. An expanded record of the Urusis Formation is found in the Swartpunt area of southern Namibia, which has yielded an important metazoan biota. But the succession in this area is preserved as a series of thrusts at the leading edge of the Gariep orogenic belt and zircon U-Pb data show systematic age repetition. We use regional stratigraphic and structural mapping, integrated with carbonate carbon isotope (δ 13 C carb ) chemostratigraphy and high-precision radioisotope U-Pb zircon geochronology from outcrop and recently acquired drill core to develop a temporally calibrated basin-wide depositional model. This integrated dataset either reflects complex zircon reworking, inheritance, or potential analytical biases (Scenario 1) or the presence of a Gariep-related cryptic décollement within the Spitskop Member that has resulted in stratigraphic repetition (Scenario 2). We investigate the evidence for and against both scenarios and consider their implications for stratigraphic and δ 13 C carb correlations between the Swartpunt area and coeval exposures along the Orange River border with South Africa. Given that these issues are in an area that hosts numerous silicified ash beds and extensive exposure, an inability to confidently discount either scenario highlights a level of compounding uncertainty in zircon U-Pb geochronology that must be considered when attempting to build global chronostratigraphic frameworks. Scenario 1 implies that some of the weighted mean ages and Bayesian eruption ages from the Swartpunt area may be >1 Myr older than the depositional age of their respective ash beds when assuming existing stratigraphic correlations. If this scenario is preferred, then a cautious approach would be to consider all weighted mean zircon U-Pb ages from ash beds to reflect maximum depositional ages. Both scenarios support deposition of the Huns Member >540 Ma in the Swartpunt area if the oldest weighted mean age reported here represents a near-depositional age, which has significant implications for the temporal calibration of important terminal Ediacaran ichnofossil assemblages and future cyclostratigraphic studies. Stratigraphic correlations common to both scenarios allow us to temporally calibrate a basin evolution model for the Nama Group. Temporal trends in initial hafnium isotope (εHf) compositions of zircon grains from ash beds throughout the succession may support progressive crustal thickening associated with underplating of the Damara orogenic belt along the northern periphery of the Kalahari craton from ca. 547 Ma to ca. 538 Ma. The compilation of new and published zircon U-Pb ages may also imply that the locus of carbonate platform development migrated from north to south (present co-ordinates), tracking the migration of foredeep subsidence.

Earth-Science Reviews↗
USGS2000

Shoals and valley plugs in the Hatchie River watershed

Agricultural land use and gully erosion have historically contributed more sediment to the streams of the Hatchie River watershed than those streams can carry. In 1970, the main sedimentation problem in the watershed occurred in the tributary flood plains. This problem motivated channelization projects (U.S. Department of Agriculture, 1970). By the mid-1980's, concern had shifted to sedimentation in the Hatchie River itself where channelized tributaries were understood to contribute much of the sediment. The Soil Conservation Service [Natural Resources Conservation Service (NRCS) since 1996] estimated that 640,000 tons of bedload (sand) accumulates in the Hatchie River each year and identified roughly the eastern two-thirds of the watershed, where loess is thin or absent, as the main source of sand (U.S. Department of Agriculture, 1986a). The U.S. Geological Survey (USGS), in cooperation with the West Tennessee River Basin Authority (WTRBA), conducted a study of sediment accumulation in the Hatchie River and its tributaries. This report identifies the types of tributaries and evaluates sediment, shoal formation, and valley-plug problems. The results presented here may contribute to a better understanding of similar problems in West Tennessee and the rest of the southeastern coastal plain. This information also will help the WTRBA manage sedimentation and erosion problems in the Hatchie River watershed. The source of the Mississippi section of the Hatchie River is in the sand hills southwest of Corinth, Mississippi (fig. 1). This section of the Hatchie River flows northward in an artificial drainage canal, gathering water from tributary streams that also are channelized. The drainage canal ends 2 miles south of the Tennessee State line. The Tennessee section of the Hatchie River winds north and west in a meandering natural channel to the Mississippi River. Although most of the Hatchie River tributaries are also drainage canals, the river's main stem has kept most of its natural character. The Hatchie River flows through a wide valley bottom occupied mostly by riverine wetland. Historically, the valley bottom has supported hardwood forests. Since publication of the first Hatchie River report (U.S. Department of Agriculture, 1970), the channel of the river has become shallower, and flooding has increased (U.S. Department of Agriculture 1986b). These wetter conditions inhibit growth of hardwoods and lead to premature hardwood mortality. The NRCS has predicted that despite efforts to control erosion in the uplands, most of the valley-bottom forest will die. '...swamping may be so prevalent as to change most of the Hatchie River Basin flood plain into a marsh condition, with the only remnants of the present bottomland hardwood timber remaining. (U.S. Department of Agriculture, 1986b). Loss of channel depth has been concentrated in short reaches near tributary mouths. At the mouths of Richland, Porters, Clover, and Muddy Creeks, navigation has become difficult for recreational users (Johnny Carlin, West Tennessee River Basin Authority, oral commun., 1998). As the low-gradient alluvial system of the Hatchie River accumulates sediment, another common outcome has been the formation of valley plugs, areas where 'channels are filled with sediment, and all the additional bedload brought downstream is then spread out over the flood plain until a new channel has been formed' (Happ, 1975). Valley plugs typically form where the slope of a sand-laden tributary decreases downstream, or where the tributary joins its parent stream (Happ and others, 1940; Diehl, 1994, 1997; Smith and Diehl, 2000).

Mississippi, Tennessee↗
USGS1999

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↗
USGS2024

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↗
USGS2024

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island↗
USGS2023

Prediction of the probability of elevated nitrate concentrations at groundwater depths used for drinking-water supply in the Puget Sound basin, Washington, 2004–19

The Puget Sound basin encompasses the 13,700-square-mile area that drains to the Puget Sound and the adjacent marine waters of Washington State. Well more than 4 million people live within the basin, with numbers continuing to increase, who rely on the basin’s natural resources including groundwater. The Puget Sound Partnership was created by a Washington State statute to implement a science-based recovery of the Puget Sound to help address impacts to these resources. As part of the recovery, the partnership developed the Puget Sound Vital Signs as measures of ecosystem health that guide the assessment of progress toward Puget Sound recovery goals. The Puget Sound Partnership Leadership Council adopted a Drinking Water Vital Sign associated with human health and quality of life, recognizing certain indicators as integral to the sustainability of Puget Sound recovery efforts. One such Vital Sign indicator was the vulnerability of groundwater throughout the aquifers of the Puget Sound basin to elevated nitrate concentrations as defined by the probability of exceeding 2 milligrams/liter (mg/L) at a specific location and well depth. The U.S. Geological Survey (USGS) led the effort to characterize groundwater vulnerability. For this study, groundwater vulnerability refers to a probability with which a contaminant applied at or near the land surface can migrate to the aquifer of interest for a given set of land-use practices. Nitrate concentration data were selected for evaluation because elevated nitrate concentrations are typically caused by anthropogenic activities and have been associated with deleterious impacts on human health. To identify groundwater vulnerability to elevated nitrate concentrations, logistic regression was used to relate anthropogenic (human associated) and natural variables to the occurrence of elevated nitrate concentrations in untreated groundwater from large public water supply system wells found within the Washington State Department of Health Sentry database. Variables that were analyzed included well depth, soil hydraulic conductivity, precipitation, population density, fertilizer application amounts, and land-use types. Statistically significant models that predicted the probabilities of groundwater nitrate concentrations greater than 2 mg/L based on the predictor variables were created for the time periods 2000–04, 2005–09, 2010–14, and 2015–19. For all time periods, well depth and a measure of the abundance of urban and agricultural land over or near the well consistently helped explain the vulnerability of the well to elevated nitrate concentrations defined as a probability of exceeding 2 mg/L of nitrate. Precipitation and (or) soil hydraulic conductivity were also important predictor variables in the models. The models for each time period were used to create maps of groundwater vulnerability at 150- and 300-foot depths throughout the Puget Sound basin. As expected, the most vulnerable locations were associated with shallower well depths and increased agriculture and urban land cover. Across all four time periods, groundwater vulnerability throughout the Puget Sound was low, with probabilities of exceeding 2 mg/L concentrations of nitrate at depths at 150 and 300 feet typically less than 50 percent. Results also found a slight decrease in probabilities of elevated nitrate concentrations throughout the basin over time. More specifically, additional statistical tests found that groundwater with probabilities of less than about 60 percent declined from 2000 to 2019 and represented more than 75 percent of the modeled Puget Sound basin aquifer. Wells with greater than 60 percent probability increased over the same time period but represented only about 25 percent of the aquifer. The maps and statistical analysis presented in the study provide valuable and informative evaluation of the vulnerability of groundwater in the Puget Sound basin to elevated nitrate concentrations. The probability maps do not represent measured nitrate concentrations in groundwater, but rather they present the probability that nitrate concentrations exceed 2 mg/L. The models and predictions from this study are a viable indicator for the Puget Sound Partnership’s Healthy Human Population—Drinking Water Vital Sign. The logistic regression modeling approach presented here benefits water managers by allowing them to assess temporal trends in a range of probabilities, explore vulnerability changes as new regional land cover and anthropogenic data are generated, and distinguish vulnerabilities at different depths within the aquifer.

Washington↗
USGS2025

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California↗
USGS2026

Treatability study to evaluate bioremediation of trichloroethene at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, 2020–22

Executive Summary Chlorinated solvents, including trichloroethene (TCE) and other chlorinated volatile organic compounds (cVOCs), are widespread contaminants that can be treated by bioremediation approaches that enhance anaerobic reductive dechlorination. Reductive dechlorination can be enhanced either through the addition of an electron donor (biostimulation) or the addition of a known dechlorinating culture (bioaugmentation) along with an electron donor. Although bioremediation has been applied at many TCE-contaminated groundwater sites, application in source zones at sites where residual dense nonaqueous phase liquid (DNAPL) is present is more limited. In this study, laboratory and field treatability tests were completed to evaluate the potential application of anaerobic bioremediation for a shallow groundwater plume containing TCE in a perched alluvial aquifer at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, which was on the National Priorities List as the New Brighton/Arden Hills Superfund site until 2019. In addition to the presence of residual DNAPL at the site, temporal variability in groundwater flow directions and input of oxygenated recharge were possible complicating factors for the application of enhanced anaerobic biodegradation in the shallow plume. The Site K plume extends beneath the footprint of Building 103, which was demolished in 2006, and soil excavations to a maximum depth of 6 feet (ft) below ground surface in 2014 were known to leave some deeper contaminated soil in place in the TCE source area. Groundwater treatment at the site, formalized as part of the 1997 Record of Decision, has been in operation since 1986 and consists of an extraction trench at the downgradient edge of the plume to collect groundwater, which is then pumped to an on-site air stripper. Groundwater concentrations in the plume have been relatively stable since treatment began, indicating a continued source of TCE in the aquifer. The desire for a destructive remedy that would enhance the removal of cVOCs in the aquifer at Site K and shorten the remediation timeframe led the U.S. Army to request that the U.S. Geological Survey conduct a groundwater treatability study to assess bioremediation. This report describes the U.S. Geological Survey bioremediation treatability study conducted during 2020–22, including pre-design site characterization to assist in formulating the bioremediation approach, laboratory experiments to support the design of the field pilot test, and implementation and 1-year performance monitoring results for the pilot test. Pre-design site characterization included the collection of soil cores for cVOC analysis and lithologic descriptions and the re-installment of three wells to obtain hydrologic measurements and initial groundwater chemistry. Relatively flat head gradients were measured at the site, and substantial decreases in water-level elevations occurred from spring to summer (May–July 2021). Continuous water-level monitoring indicated a rapid response to precipitation. Groundwater flow velocities were consistently less than 0.5 foot per day, and the pilot bioremediation test was therefore designed with short lateral distances (about 5 ft) between injection and individual monitoring points. Soil analyses confirmed that high volatile organic compound contamination was left in place in the source area. The highest concentrations were near or in clay at the base of the perched aquifer. Concentrations of cVOCs measured in the replaced wells were consistent with historical data and had a maximum TCE concentration of 57,700 micrograms per liter (μg/L), indicative of nearby residual DNAPL based on the general rule of observed concentrations exceeding 1 percent of solubility. The primary TCE daughter product detected was 1,2-cis-dichloroethene (cisDCE), which indicated limited reductive dechlorination in the plume. Groundwater in both the source and downgradient areas was relatively reducing during the pre-design characterization, particularly in the source area where methane concentrations greater than 400 μg/L were measured. Initial laboratory tests conducted using native aquifer microorganisms from the three replacement wells showed that anaerobic TCE biodegradation rates were low when biostimulated with the addition of sodium lactate as an electron donor, also known as a carbon donor, and resulted in the production of only cisDCE. Addition of a known dechlorinating culture, WBC-2, however, resulted in rapid biodegradation and production of ethene, verifying complete reductive dechlorination of TCE. Microcosms constructed with aquifer soil collected from the site were used to evaluate other electron donors besides lactate to support reductive dechlorination by WBC-2, including corn syrup as an alternative fast-release compound and whey, soy-based vegetable oil, and 3-D Microemulsion (Regenesis, San Clemente, California) as slow-release compounds. First-order rate constants for total organic chlorine removal in these WBC-2 amended microcosms were greatest with either lactate or vegetable oil as the donor, ranging between 0.061 and 0.047 per day or corresponding half-lives of 11–15 days. Testing of commercial products in other WBC-2-bioaugmented microcosms led to selection for the field pilot test of an emulsified vegetable oil product that also contained some sodium lactate as a fast-release donor. Delaying the addition of WBC-2 relative to the donor in the microcosms resulted in the most rapid overall biodegradation rates. The selected design for the pilot test utilized three separate test plots, each about 30-ft wide and 60-ft long: plots GS1 and GS2 in the source area of the plume and plot GS3 in the downgradient area of the plume near the excavation trench. Each test plot had one injection well, one monitoring well upgradient from the injection point, and 12 surrounding monitoring wells in a grid to capture variable groundwater flow directions. Donor injections, which included a bromide tracer, were completed in October 2021, immediately following baseline sampling, and the WBC-2 culture was injected about 40 days later, between November 30 and December 2, 2021. Performance monitoring conducted until December 2022 included hydrologic measurements and analyses of cVOCs, redox-sensitive constituents, dissolved organic carbon, bromide, volatile fatty acids, compound-specific carbon isotopes, and microbial communities. The biogeochemical data collected during the pilot tests in the three treatment plots showed that enhanced, complete reductive dechlorination of cVOCs in the groundwater was achieved in the GS1 and GS3 plots. In contrast, evidence of distribution of the injected amendments and subsequent biodegradation was limited in GS2, which was in an area of more heterogeneous soil lithology and low water table elevations. The molar composition of volatile organic compounds in the GS1 and GS3 plots was dominated by ethene in wells that were reached by the injected amendments by the end of the monitoring period. In the GS1 and GS3 plots, similar patterns were observed of cVOC concentrations decreasing to near detection levels, or below, at some wells sampled in July and October 2022, whereas ethene became dominant and indicated sustained complete reductive dechlorination. Baseline cVOC concentrations were more than a factor of 10 higher in the groundwater in the GS1 plot than in GS3, but no apparent inhibition of complete dechlorination occurred. As expected from the initial pre-design site data and the laboratory experiments, enhanced dissolution of residual DNAPL coupled to biodegradation was evident in the GS1 plot, where a marked increase in dichloroethene (DCE) above the initial baseline and upgradient TCE and DCE concentrations occurred. DCE concentrations subsequently declined where DNAPL dissolution was evident, concurrent with production of vinyl chloride and then predominantly ethene. Thus, overall biodegradation rates outpaced the DNAPL dissolution and desorption and DCE production in the source area. This success in complete degradation to predominantly ethene was achieved even in areas where the DCE concentrations reached a maximum of about 30,000 μg/L. Compound specific isotope analysis of carbon in TCE, cisDCE, trans-1,2-dichloroethene, and vinyl chloride was conducted to provide another line of evidence of the occurrence and extent of anaerobic biodegradation. Along a flow path in each plot that was affected by the injected amendments, carbon isotopes in the TCE and daughter cVOCs in the groundwater became isotopically heavier, indicating biodegradation. Enhanced biodegradation rates calculated from the field tests in GS1 and GS3 showed half-lives of 36.9–75.3 days for DCE degradation and 9.48–38.5 days for ethene production. Notably, these ethene production rates calculated from the field tests are consistent with the results of WBC-2-bioaugmented microcosms amended with either lactate or vegetable oil, which had half-lives for total organic chlorine removal that ranged from 11 to 15 days. These rates indicated rapid enhanced biodegradation, which is promising for application of a full-scale bioremediation remedy. Ultimately, however, the mass of residual or sorbed TCE in the aquifer that remains accessible for dissolution and biodegradation would likely control the time required for a full-scale bioremediation effort to achieve performance goals for TCE and cisDCE specified in the Record of Decision for Site K. The field pilot tests showed that the relatively low hydraulic head gradients and temporal changes in groundwater flow directions in the shallow aquifer would add complexity to a full-scale bioremediation effort. The radius of influence (ROI) at GS1 and GS3 (16.3 ft and 12.7 ft, respectively) were close to the design ROI of 15 ft. The estimated ROI at GS2 was about four times the design ROI, but may be less reliable at this location owing to groundwater flow direction. In addition, the low temperatures following WBC-2 injection in late November to early December 2021, in combination with the low hydraulic head gradients, were probably major factors in the delay observed before the onset of enhanced biodegradation following injection of the culture. Additional test injections could be beneficial to optimize the timing of donor and culture injections with the variable temperatures and hydraulic head in the shallow aquifer.

Minnesota↗
USGS1973

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida↗