USGS ScienceSearch

SEARCH · USGS Science

Results for “Transactions of the American Fisheries Society”

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,688 records · Page 57Linked to original sources

Sedimentary exhalative (sedex) zinc-lead-silver deposit model

This report draws on previous syntheses and basic research studies of sedimentary exhalative (sedex) deposits to arrive at the defining criteria, both descriptive and genetic, for sedex-type deposits. Studies of the tectonic, sedimentary, and fluid evolution of modern and ancient sedimentary basins have also been used to select defining criteria. The focus here is on the geologic characteristics of sedex deposit-hosting basins that contain greater than 10 million metric tons of zinc and lead. The enormous size of sedex deposits strongly suggests that basin-scale geologic processes are involved in their formation. It follows that mass balance constraints of basinal processes can provide a conceptual underpinning for the evaluation of potential ore-forming mechanisms and the identification of geologic indicators for ore potential in specific sedimentary basins. Empirical data and a genetic understanding of the physicochemical, geologic, and mass balance conditions required for each of these elements are used to establish a hierarchy of quantifiable geologic criteria that can be used in U.S. Geological Survey national assessments. In addition, this report also provides a comprehensive evaluation of environmental considerations associated with the mining of sedex deposits.

Scientific Investigations Report

Techniques for estimating magnitude and frequency of floods on streams in Indiana

Equations are presented for estimating the magnitude and frequency of floods at ungaged sites on unregulated and nonurban streams in Indiana. The equations were developed by multiple-regression, analysis of basin characteristics and peak-flow statistical data from 242 gaged locations in Indiana, Ohio, and Illinois. The State of Indiana was divided into seven areas on the basis of the regression analysis. A set of equations for estimating peak discharges with recurrence intervals of 2, IO, 25, 50, and 100 years was developed for each area. Significant basin characteristics in the equations are drainage area, channel length, channel slope, mean annual precipitation, storage, precipitation intensity, and a runoff coefficient. Standard errors of estimate for the equations range from 24 to 45 percent. Methods are also presented for estimating flood magnitude and frequency at sites on gaged streams. Flood-frequency data based on observed peaks are given for 270 gaged locations. Twenty of these are on regulated streams, and six are on urban streams. Basin characteristics are also included car 245 of the gaged locations on unregulated and nonurban streams. No techniques are given for estimating flood magnitude and frequency at ungaged sites on regulated or urban streams. A rainfall-runoff model was tlsed to synthesize long-term peak data at 11 gaged locations on small streams. Flood-frequency curves developed from the long-term synthetic data were combined with curves based on short-term observed data to provide weighted estimates of flood magnitude and frequency at the rainfall-runoff stations.

Indiana

Base flow and ground water in upper Sweetwater Valley, Tennessee

The upper Sweetwater Valley area has a flow system with complex interaction between surface and ground water. A water budget study indicated that during dry years approximately three-fourths of the annual flow to Sweetwater Creek may be derived from ground-water sources. Hydrograph analysis showed seasonal variation of recharge to the ground-water flow system. Streamflow records were analyzed to estimate the frequency of low flow of Sweetwater Creek at river mile 16.7, and indicated the lowest average flow for 1 day in 20 years to be about 5.1 cubic feet per second. Two periods of base-flow measurements of Sweetwater Creek identified channel reaches with significant gains and losses of streamflow. Base flow measurements also showed interbasin transfer of water among sub-basins of the valley. Major flow surpluses were associated with areas in which the majority of flow originated at a spring. Topographically low areas adjacent to the main stem of Sweetwater Creek generally have surplus flow. Topographically higher areas generally have deficient surface outflow unless significant spring flow occurs in the basin. Ground-water recharge occurs by water draining into sinkholes, faults, and fractures. Ground-water flow is regionally diffused across formation strikes from the topographically low areas unless the water is exposed to highly permeable formations or impervious formations. Ground water infiltrates the highly permeable formations and flows along strike. Ground water encountering impervious formations may discharge at small springs at the contact, or may reroute along the contact if the upgradient rock is sufficiently permeable or has well developed secondary porosity. Groundwater discharges to streams at innumerable seeps and at a few large springs. Areas of ground-water flow up-gradient of large springs are hypothesized as likely areas of significant ground-water reservoirs.

Tennessee

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution

Assessing legacy nitrogen in groundwater using numerical models of the Long Island aquifer system, New York

Nitrogen transported along groundwater flow paths in coastal aquifers can contribute substantially to nitrogen loading into surface water receptors, particularly in hydrologic systems dominated by groundwater discharge. Nitrogen entrained in the aquifer is a function of land use and associated nitrogen sources at the time of groundwater recharge, which may differ considerably from present-day sources. Legacy nitrogen can result in substantial discrepancies between observed present-day nitrogen loading to surface water receptors and loading estimated from present-day sources. Additionally, legacy nitrogen can continue to discharge into surface waters after nitrogen mitigation actions have been undertaken. Here, we use a numerical modeling framework to compare three methods of estimating time-varying historical nitrogen loads to four water bodies (receptors) on eastern Long Island, New York. The methods span a range of data requirements and process complexity, from instantaneous receptor loads calculated from steady-state groundwater contributing areas, to transient loads estimated by explicitly simulating legacy groundwater nitrogen transport over a century with large changes in nitrogen sources and hydrologic conditions. The effects of legacy nitrogen on estimated receptor loads varied temporally and spatially within the study area. Depending on antecedent nitrogen inputs and hydrologic conditions, historical annual nitrogen loads estimated from transient simulations accounting for legacy nitrogen can be quite similar (<10% difference) or substantially different (±100%) from those estimated from simpler instantaneous methods. Continued input of present-day nitrogen sources using methods that account for legacy nitrogen results in asymptotic increases in receptor nitrogen loads over time, indicating that simulated present-day receptor nitrogen loads are not in equilibrium with present-day inputs. For these receptors in disequilibrium, models simulating transient groundwater nitrogen transport could be used to account for legacy nitrogen lag times to help resource managers evaluate the potential effectiveness of proposed nitrogen mitigation actions.

EarthArXiv

Deformation mechanisms in quartz veins and shear zones elucidate the origin of gold mineralization at Pogo, Alaska

Pogo is a quartz vein hosted, ca. 8 Moz gold deposit. Although it has similarities to orogenic and magmatic-hydrothermal deposits, its origin remains enigmatic. Observations from surface exposures, underground workings, and drill core provide new constraints on quartz vein origins with implications for mineralization. Abundant, largely barren metamorphic segregation quartz veins are found throughout interior Alaska and share characteristics with mineralized quartz veins at Pogo. Pogo quartz veins show crystal plastic deformation fabrics consistent with greenschist to amphibolite facies metamorphism and a lack of internal shear. In contrast, gold in Pogo quartz veins occurs in brittle microfaults and fractures with sulphides such as arsenopyrite and pyrrhotite. Major and minor cataclastic shear zones also exist at Pogo. Cataclastic shear zones commonly cut the mineralized quartz veins and porphyroclasts are dominated by quartz with the same plastic deformation fabrics found in the major Pogo quartz veins. The porphyroclasts, and the carbonaceous clay-rich matrix they sit in, also contain gold indicating that the shear zones postdate quartz and earliest gold deposition. These observations suggest that competency contrasts between the quartz veins and their phyllosilicate-rich host rocks controlled transient permeability formed by late, preferential brittle deformation events localizing sulphide ± gold deposition.

Alaska

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California

Characterization and validation of tidally calibrated strains from the Alto Tiberina Near Fault Observatory Strainmeter Array (TABOO-NFO-STAR)

Six horizontal borehole tensor strainmeters (TSM1-6) installed from Fall 2021 to Spring 2022 comprise the Alto Tiberina Near Fault Observatory Strainmeter Array (STAR), providing an unprecedented opportunity to investigate seismic and aseismic deformation from hazardous high- and low-angle normal faults in Italy. Prior to use in tectonic applications, they require in-situ calibration and correction for non-tectonic signals. We tidally calibrate the instruments, characterize the calibration uncertainty, and test the results against environmental and earthquake signals originating from local to teleseismic distances. The STAR sites demonstrably deviate from assumptions common to the standard manufacturer's calibrations, including negative areal coupling at TSM3-6. While the tidally calibrated strains have ~3-56% uncertainty, the calibrated dynamic strains show interstation precision and accuracy to nanostrain levels, and static coseismic offsets in the array footprint are within uncertainty. TSM3 records a complex series of strains that may arise from dynamically triggered near-borehole fracture slip and fluid flow that does not appear to affect its sensitivity to lower strain rate deformation. Future calibration improvement may be afforded with longer stable timeseries, particularly for TSM4. Overall, our analyses demonstrate expanded geodetic capability for detecting deformation in the Alto Tiberina Near Fault Observatory.

Alto Tiberina Near Fault Observatory

Uranium and trace elements in stream sediments as an exploration tool

More than 45 trace elements have been reported in anomalous amounts in the uranium ores. The specific suite of elements associated with any one uranium deposit varies according to deposit type and geologic province. The primary geochemical halo of uranium and associated trace elements in the host rock, together with secondary dispersion halo in soils and alluvium, offers a potential geochemical exploration target. Sediment from streams near low-grade uranium occurrences in arkosic sandstone of the Denver basin, Colorado, and sediment from streams in the igneous and metamorphic terrace near the Midnite mine of eastern Washington have been analyzed for acid-extractable U, Cu,Pb, Zn, Ni, Co, Fe, Mn, Mo, Cd, Cr, Ag, V, and Se, and for total Fe, Mn, Hg, As, and organic C. The -80 mesh fraction of the alluvium downstream from a small uraniferous limonite occurrence in the Denver basin shows anomalous concentrations of U, Pb, As, Ni, Mo, V, Zn, and Se. Anomalous concentrations of U persist farther downstream (as much as 1,000 m) than any of the trace elements. In another stream system southeast of Denver, above-background levels of U, Pb, As, Zn, V, Fe, and Mn occur near a uranium anomaly (defined by other methods) which has no surface mineralization. In the Midnite mine area, sediments of streams draining prospect areas showed anomalous concentrations of uranium (as much as 10 times background), but no accompanying anomalous trace-element concentrations, although analyses of Midnite mine ores show that a broad suite of trace elements are associated with uranium.

Open-File Report

Simulated effects of future water availability and protected species habitat in a perennial wetland, Santa Barbara County, California

This study evaluates the potential water availability in Barka Slough and the effects of changing hydrological conditions on the aquatic habitat of five protected species. Barka Slough is a historically perennial wetland at the downstream western end of the San Antonio Creek Valley watershed (SACVW). A previously published hydrologic model of the SACVW for 1948–2018 was extended to include 2019–2021 and then modified to simulate the future years of 2022–2051. Two models simulating the future years of 2022–2051 were constructed, each with different climate inputs: (1) a repeated historical climate and (2) a 2070-centered Drier Extreme Warming climate (2070 DEW). The model with the 2070 DEW climate had warmer temperatures and an increase in average annual precipitation driven by larger, albeit more infrequent, precipitation events than the model with the historical climate. Simulated groundwater pumpage resulted in cumulative groundwater storage depletion and groundwater-level decline in Barka Slough in both future models. The simulations indicate that Barka Slough may transition from a perennial to an ephemeral wetland. Streamflow, stream disconnection, and depth to groundwater are key habitat metrics for federally listed species in Barka Slough. Future seasonal conditions for each metric are more likely to affect federally listed species’ habitats under 2070 DEW climatic conditions. Future seasonal streamflow volume may negatively impact unarmored threespine stickleback ( Gasterosteus aculeatus williamsoni ) and tidewater goby ( Eucyclogobis newberryi) habitats. Future seasonal stream disconnection may negatively impact the unarmored threespine stickleback habitat. Future groundwater-level decline may negatively impact Gambel’s watercress ( Nasturtium gambelii ) and La Graciosa thistle ( Cirsium scariosum var. loncholepis ) habitats and could influence the ability to use Barka Slough as a restoration or reintroduction site for these species. Results from this study can be used to inform water management decisions to sustain future groundwater availability in the SACVW.

California

Preface to the focus section on intraplate earthquakes

More than a half century after plate tectonics provided an overarching framework to explain earthquakes along active plate boundaries, numerous theories have been proposed to explain where, why, and how often earthquakes occur well away from active plate boundaries, but a paradigm remains elusive. Even the classification of earthquakes away from active plate boundaries as "intraplate" raises issues, with potentially important distinctions between Stable Continental Regions and more actively deforming regions including passive margins and failed rifts. Some of the largest known intraplate earthquakes themselves remain enigmatic, having occurred before the modern instrumental era. Hazard assessments are often data-limited: low fault-slip rates relative to landscape modification rates result in poor discoverability of fault sources, challenging the characterization of source zones and earthquake recurrence; the completeness and homogenization of instrumental earthquake catalogs using uncertain magnitude conversions can lead to uncertainties in earthquake recurrence; and, limited strong-motion observations for large-magnitude events at near-source distances leads to uncertainties in the selection and development of ground-motion models for seismic hazard studies. Data from recent intraplate earthquakes around the world—from the moment magnitude M 7.7 2001 Bhuj, India, earthquake 25 years ago to the 2024 M 4.8 Tewksbury, New Jersey earthquake—have yielded both new insights and new questions. The papers in this special focus discuss many of the long-standing challenges involved with intraplate earthquake investigations and provide a snapshot of the state of the art with current research to advance our understanding.

Seismological Research Letters

Suspended-sediment characteristics in the Housatonic River Basin, western Massachusetts and parts of eastern New York and northwestern Connecticut, 1994-96

Suspended-sediment concentrations, discharges, loads, and yields were determined for eight subbasins in the Housatonic River Basin in western Massachusetts, eastern New York, and northwestern Connecticut from April 1994 through March 1996. Suspended-sediment samples were collected at three continuous-record sediment stations and at four partial-record sediment stations. Suspended-sediment concentrations in samples collected during the period of study ranged from less than 0.5 to 3,400 milligrams per liter, and concurrent streamflows ranged from 0.03 to 126 cubic feet per second per square mile at the seven stations. Median suspended-sediment concentrations in samples collected at each station ranged from 7 to 61 milligrams per liter. Median streamflows during suspended-sediment sampling ranged from 1.86 to 5.88 cubic feet per second per square mile. Instantaneous suspended-sediment yields ranged from less than 0.005 to 185 tons per day per square mile, and medians ranged from 0.03 to 1.12 tons per day per square mile at the seven stations. Total suspended-sediment loads (mass) from April 1994 through March 1996 at the continuous-record sediment stations were 11,603 tons at Housatonic River near Great Barrington, 7,929 tons at Green River, and 54,347 tons at Housatonic River near Ashley Falls. Suspended-sediment load during January 1996 at the Green River station accounted for about 54 percent of the total suspended-sediment load for the Green River during the 2-year study. Suspended-sediment load on January 19 and 20, 1996, at the Green River station accounted for about 50 percent of the January 1996 suspended-sediment load, or about 27 percent of the total suspended-sediment load during the 2-year study. This large suspended-sediment transport was the result of rainfall and snowmelt on January 19 and 20, 1996--the equivalent of a 5- to 6-inch rain storm in the Green River subbasin. Total suspended-sediment loads during the 2-year study at the partial-record sediment stations were 3,052 tons at Williams River, 1,758 tons at Ironworks Brook, and 17,927 tons at Konkapot River. Suspended-sediment yields from April 1994 through March 1996 at the continuous-record sediment stations were 21 (tons/yr)/mi 2 at Housatonic River near Great Barrington, 78 (tons/yr)/mi 2 at Green River, and 58 (tons/yr)/mi 2 at Housatonic River near Ashley Falls. Suspended-sediment yields during the 2- year study at the partial-record sediment stations were 35 (tons/yr)/mi 2 at Williams River, 78 (tons/yr)/mi 2 at Ironworks Brook, and 147 (tons/yr)/mi 2 at Konkapot River. Suspended-sediment yields were estimated for two subbasins in the Housatonic River Basin--Schenob Brook at Sheffield, a partial-record sediment station, and the area adjacent to the Housatonic River between Great Barrington and Ashley Falls. The estimate of suspended-sediment yield at Schenob Brook at Sheffield of 82 (tons/yr)/mi 2 is comparable to the yield determined for the Green River and Ironworks Brook. The estimate of suspended-sediment yield for the area adjacent to the Housatonic River between Great Barrington and Ashley Falls was 395 (tons/yr)/mi 2 . This estimated suspended-sediment yield was 2.7 to 18.8 times greater than that estimated for any of the other subbasins. Several basin and land-use characteristics thought to affect suspended-sediment transport in the subbasins were compared to the suspended-sediment yields. The characteristics that seemed to affect suspended-sediment discharge were dams, which contributes to decreased yields; Hadley, Limerick, Linlithgo, Saco, and Winooski silt loam soils (high erodibility soils), which contributes to increased yields; stratified-drift deposits, which contributes to increased yields; and agricultural and open land, which contributes to increased yields. The effect of stratified-drift deposits on suspended-sediment discharge is thought to be greater when those deposits are of glaciolacustrine (generally clay, silt, and fine sand), rather than glaciofluvial (clay, silt, sand, gravel, and cobbles) origin. The silt loam soils, glaciolacustrine deposits, and agricultural and open land are interrelated, inasmuch as the silt loam soils generally are associated with glaciolacustrine deposits and agricultural activities.

Connecticut, Massachusetts, New York

Spatial differences in soil nutrients along a hydrographic gradient on floodplains in Dongting Lake

The spatial heterogeneity of soil nutrients is crucial for the water bird and whole floodplain wetland ecosystem in large lakes, and it is influenced by the dramatic water level changes and sedimentation progress in West Dongting Lake (WDL). Soil samples were collected at various soil depths along the Yuan River and Li River that feed into WDL. The concentrations of soil total organic carbon (TOC), total nitrogen (TN), total phosphorus (TP), and soil grain size were tested. The stoichiometric ratios of C, N, P, and the mean value of soil grain size ( Mz ) were calculated. The differences of soil TOC, TN, TP and the stoichiometric ratio at different sites and soil depths were compared. Linear regression was used to explore the relationships of Mz and nutrient concentrations, and relationships between TOC, TN, and TP. Redundancy analysis was used to explore the relationship between soil nutrients, heavy metal concentrations, and plant community diversity. The results showed that the distributions of soil TOC, TN, and TP concentrations differed across regions in west Dongting Lake along the Yuan and Li Rivers. Total organic carbon concentration differed at different sedimentation depths. Soil grain size showed negative effect with soil TOC, TN, and TP concentrations in this region. Plant community diversity correlated positively with soil TOC and negatively with Hg. West Dongting Lake was N limited despite the high wet deposition of N. It could potentially be attributed to the insufficient presence of aerobic environments for microbes during intermittent flooding of the floodplain, coupled with feeble mineralization. This study can provide valuable insights for the conservation of water bird habitats and wetland ecosystems.

West Dongting Lake

Status of water-quality conditions in the United States, 2010–20

Degradation of water quality can make water harmful or unusable for humans and ecosystems. Although many studies have assessed the effect of individual constituents or narrow suites of constituents on freshwater systems, no consistent, comprehensive assessment exists over the wide range of water-quality effects on water availability. Using published studies, data, and models completed at regional or national scales in the United States during 2010–20, this chapter moves towards a comprehensive assessment by summarizing how selected anthropogenic and geogenic water-quality constituents affect national-scale water availability for human and ecosystem needs. Several types of human health, agricultural, ecological, and beneficial-use standards or thresholds were used to provide context for categorizing surface-water and groundwater quality. Water availability for human and ecological use is limited by elevated concentrations of geogenic and anthropogenic constituents in surface and groundwater. Elevated concentrations of five geogenic constituents (arsenic, manganese, strontium, radium, and adjusted gross alpha) are common in groundwater and collectively affect the drinking water supply to over 30 million people. Surface water sourced drinking water supplies are impaired in about a third of assessed stream miles, most commonly because of non-mercury metals and salinity. Health-based violations at community water systems may disproportionately affect socially vulnerable communities. Ecological water uses are predominantly limited by nutrients, sediment, temperature, pathogens, salinity, and pesticides. Water availability for human and ecological use is adversely affected by human activities including human contaminant sources (for example, wastewater, agriculture), processes (for example, dredging, groundwater pumping), or permanent landscape modifications (for example, dams, urbanization). Primary contaminant sources vary spatially and include fertilizer and manure, atmospheric deposition, wastewater treatment plants, urban land, and a range of natural sources. Contaminants of emerging concern, contaminants without regulatory thresholds, and mixtures of geogenic and anthropogenic water contaminants also contribute to ecological degradation and human exposure.

conterminous United States

Controls on natural hydrogen generation during serpentinization of mantle rocks

Mantle rocks undergoing serpentinization can generate significant amounts of natural hydrogen, yet the rates and controlling processes remain poorly understood. Here, we constrain the possible hydrogen generation rates in two distinct mantle rock types, the fertile lherzolites of the Western Pyrenees and the depleted harzburgites of Northern California, to relatively low rates of ~0.1 to ~0.5 tonnes H₂ yr⁻¹ km⁻³ of reactive rock. When integrated over the full reactive volumes, this corresponds to total production rates of ~300 to ~600 tonnes H₂ yr⁻¹. By combining three-dimensional geophysical inversion with numerical modelling of fluid-rock processes, we show that hydrogen generation rates are mainly limited by H₂ saturation in the fluid and reaction kinetics. Under these constraints, hydrogen generation in mantle-derived serpentinization systems proceeds slowly, making rapid large-scale replenishment unlikely and suggesting that large, economically relevant accumulations, would require timescales of thousands to tens of thousands of years to develop.

northern California, western Pyrenees

Avak Creek oil occurrence, North Slope, Alaska: Newly discovered oil seep on Native lands, near village of Utqiagvik

An unknown occurrence of oil was detected near Avak Creek on Native lands on the North Slope of Alaska. Determining the source of oil was imperative for allowing stakeholders (Federal, State, and local government agencies and the landowner, an Alaska Native corporation) to make timely and informed decisions and mount a mitigation response, if required. The regional and local geological framework of the Avak Creek site was constructed using seismic surveys, well data, and basin modeling results, to identify local petroleum systems, map structural geometry and faults, define source rock thermal maturity distributions, and infer likely oil-migration pathways. Molecular hydrocarbon fingerprints (biomarkers, diamondoids, compound-specific isotopes) of the oil were compared to those of local and regional oil seeps, exploration well tests, and produced oils. Biomarker acid distributions characterized the history and extent of petroleum biodegradation. Integrating subsurface and geochemical parameters, the oil is interpreted to be a natural seep generated locally, predominantly from the Brookian Lower Cretaceous Hue Shale/gamma-ray zone, rather than an anthropogenic source of pollution. Results highlight sophisticated analytical technologies used to characterize complex, compositionally altered hydrocarbons. Results also advance our understanding of Brookian source rock distribution, subsurface petroleum migration pathways, and Arctic Alaska petroleum systems.

Alaska

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida