USGS Science⌕ Search

SEARCH · USGS Science

Results for “Processes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources
USGS2005

Animating the biodynamics of soil thickness using process vector analysis: A dynamic denudation approach to soil formation

This paper expands the dynamic denudation framework of landscape evolution by providing new process insights and details on how soil and its signature morphological feature, the biomantle, form and function in the environment. We examine soils and their biomantles from disparate parts of the world, from the tropics through midlatitudes and hyperarid through perhumid, a range that exhibits varying environments for, and of, life. We then explicate the process pathways that cause soils to thicken and thin, and to even disappear, then reform. We do this by examining thickness relationships, where soil thickness stand biomantle thickness bt are functions of upbuilding u and deepening d minus removal r processes, hence st/bt=f(u+d-r). Upbuilding has two subsets, u1, which includes all exogenous (allochthonous-outside) mineral and/or organic inputs to the soil system, and u2, which includes all endogenous (autochthonous-in situ) processes and productions, including weathering. Exogenous u1 inputs include eolian and slopewash inputs (sedimentations) of mineral and organic materials, mass wasting accumulations and the like. Endogenous u2 processes and productions include the sum of in situ bioturbations, biosynthetic productions, organic accumulations, biovoid productions, weathering and volume increases caused by their sum. Endogenous upbuildings, which dominantly occur in the biomantle, are basically biodynamic bd processes and productions, hence u2=bd. Therefore, if exogenous upbuildings u1 are minimal or zero, then biomantle thickness bt is expressed by bt=f(u2-r) or bt=f(bd-r). Drawing on these relationships, we employ a graphic-conceptual device called process vector analysis in a digital animation (see supplementary materials or cf. https://netfiles.uiuc.edu/jdomier/www/temp/ biomantle.html) that illustrates the main pathways that form both Earth's soil and its unique epidermis, the biomantle. We then discuss the main elements of the animation using still frames that represent thickness turning points. We end by encouraging researchers, instructors and students to view Earth's landforms and soils as integrated dynamic entities that are constantly coevolving elements of a global subaerial-subaqueous continuum. ?? 2004 Elsevier B.V. All rights reserved.

Geomorphology↗
USGS2010

Interactions of solutes and streambed sediment: 2. A dynamic analysis of coupled hydrologic and chemical processes that determine solute transport

Solute transport in streams is determined by the interaction of physical and chemical processes. Data from an injection experiment for chloride and several cations indicate significant influence of solutestreambed processes on transport in a mountain stream. These data are interpreted in terms of transient storage processes for all tracers and sorption processes for the cations. Process parameter values are estimated with simulations based on coupled quasi-two-dimensional transport and first-order mass transfer sorption. Comparative simulations demonstrate the relative roles of the physical and chemical processes in determining solute transport. During the first 24 hours of the experiment, chloride concentrations were attenuated relative to expected plateau levels. Additional attenuation occurred for the sorbing cation strontium. The simulations account for these storage processes. Parameter values determined by calibration compare favorably with estimates from other studies in mountain streams. Without further calibration, the transport of potassium and lithium is adequately simulated using parameters determined in the chloride-strontium simulation and with measured cation distribution coefficients.

Water Resources Research↗
USGS2007

An innovative carbonate coprecipitation process for the removal of zinc and manganese from mining impacted waters

Although mine drainage is usually thought of as acidic, there are many cases where the water is of neutral pH, but still contains metal species that can be harmful to human or aquatic animal health, such as manganese (Mn) and zinc (Zn). Typical treatment of mine drainage waters involves pH adjustment, but this often results in excessive sludge formation and removal of nontoxic species such as magnesium and calcium. Theoretical consideration of the stability of metal carbonate species suggests that the target metals could be removed from solution by coprecipitation with calcium carbonate. The U.S. Geological Survey has developed a limestone-based process for remediation of acid mine drainage that increases calcium carbonate saturation. This treatment could then be coupled with carbonate coprecipitation as an innovative method for removal of toxic metals from circumneutral mine drainage waters. The new process was termed the carbonate coprecipitation (CCP) process. The CCP process was tested at the laboratory scale using a synthetic mine water containing 50 mg/L each of Mn and Zn. Best results showed over 95% removal of both Mn and Zn in less than 2 h of contact in a limestone channel. The process was then tested on a sample of water from the Palmerton zinc superfund site, near Palmerton, Pennsylvania, containing over 300 mg/L Zn and 60 mg/L Mn. Treatment of this water resulted in removal of over 95% of the Zn and 40% of the Mn in the limestone channel configuration. Because of the potential economic advantages of the CCP process, further research is recommended for refinement of the process for the Palmerton water and for application to other mining impacted waters as well. ?? Mary Ann Liebert, Inc.

Environmental Engineering Science↗
USGS2021

Great expectations: Deconstructing the process pathways underlying beaver-related restoration

Beaver-related restoration is a process-based strategy that seeks to address wide-ranging ecological objectives by reestablishing dam building in degraded stream systems. Although the beaver-related restoration has broad appeal, especially in water-limited systems, its effectiveness is not yet well documented. In this article, we present a process-expectation framework that links beaver-related restoration tactics to commonly expected outcomes by identifying the set of process pathways that must occur to achieve those expected outcomes. We explore the contingency implicit within this framework using social and biophysical data from project and research sites. This analysis reveals that outcomes are often predicated on complex process pathways over which humans have limited control. Consequently, expectations often shift through the course of projects, suggesting that a more useful paradigm for evaluating process-based restoration would be to identify relevant processes and to rigorously document how projects do or do not proceed along expected process pathways using both quantitative and qualitative data.

BioScience↗
USGS2026

Bridging ecology and geosciences in riverscapes: Implications for process-based restoration

There has been a growing interest in integrating geological and ecological processes for sustainable river management and restoration. Lotic systems are shaped by diverse physical processes, including geology, geomorphology, hydrology, and interactions with terrestrial processes. However, restoration practices often prioritize specific habitats or river forms without fully considering the underlying physical processes that support biological communities and ecosystem functions. In this paper, we synthesize studies that integrally examine geological and ecological processes across different scales and components of riverscapes, including geohydrological processes, which have been less investigated in riverscape studies. We begin by examining processes at broader spatial scales, including river–watershed and river–riparian interactions, and gradually narrow our focus to the dynamics that occur among habitats within river channels, through which we highlight the significance of conceptualizing rivers as dynamic “networks” rather than linear features. Finally, we identify both scientific and practical challenges that can be addressed to bridge the gap between basic-science implications and their implementation in riverscape restorations.

Ecological Research↗
USGS2023

GeoAI for spatial image processing

The development of digital image processing, as a subset of digital signal processing, depended upon the maturity of photography and image science, introduction of computers, discovery and advancement of digital recording devices, and the capture of digital images. In addition, government and industry applications in the Earth and medical sciences were paramount to the growth of the technology. From the early days when photography was first introduced to science to today, artificial intelligence and deep learning technologies have been intensively used to analyze imagery. Spatial image processing has experienced breakthroughs and evolutions. This chapter presents an overview of the history of image processing, GeoAI-based image processing applications, and the role of GeoAI in advancing image processing methods and research. We also discussed the remaining challenges to using GeoAI for image processing regarding training data annotation, the issues of scale, resolution, and change in space over time.

Book chapter↗
USGS1990

Evaluation of three electronic report processing systems for preparing hydrologic reports of the U.S Geological Survey, Water Resources Division

In 1987, the Water Resources Division of the U.S. Geological Survey undertook three pilot projects to evaluate electronic report processing systems as a means to improve the quality and timeliness of reports pertaining to water resources investigations. The three projects selected for study included the use of the following configuration of software and hardware: Ventura Publisher software on an IBM model AT personal computer, PageMaker software on a Macintosh computer, and FrameMaker software on a Sun Microsystems workstation. The following assessment criteria were to be addressed in the pilot studies: The combined use of text, tables, and graphics; analysis of time; ease of learning; compatibility with the existing minicomputer system; and technical limitations. It was considered essential that the camera-ready copy produced be in a format suitable for publication. Visual improvement alone was not a consideration. This report consolidates and summarizes the findings of the electronic report processing pilot projects. Text and table files originating on the existing minicomputer system were successfully transformed to the electronic report processing systems in American Standard Code for Information Interchange (ASCII) format. Graphics prepared using a proprietary graphics software package were transferred to all the electronic report processing software through the use of Computer Graphic Metafiles. Graphics from other sources were entered into the systems by scanning paper images. Comparative analysis of time needed to process text and tables by the electronic report processing systems and by conventional methods indicated that, although more time is invested in creating the original page composition for an electronically processed report , substantial time is saved in producing subsequent reports because the format can be stored and re-used by electronic means as a template. Because of the more compact page layouts, costs of printing the reports were 15% to 25% less than costs of printing the reports prepared by conventional methods. Because the largest report workload in the offices conducting water resources investigations is preparation of Water-Resources Investigations Reports, Open-File Reports, and annual State Data Reports, the pilot studies only involved these projects. (USGS)

Open-File Report↗
USGS2006

User's guide to the Variably Saturated Flow (VSF) process to MODFLOW

A new process for simulating three-dimensional (3-D) variably saturated flow (VSF) using Richards' equation has been added to the 3-D modular finite-difference ground-water model MODFLOW. Five new packages are presented here as part of the VSF Process--the Richards' Equation Flow (REF1) Package, the Seepage Face (SPF1) Package, the Surface Ponding (PND1) Package, the Surface Evaporation (SEV1) Package, and the Root Zone Evapotranspiration (RZE1) Package. Additionally, a new Adaptive Time-Stepping (ATS1) Package is presented for use by both the Ground-Water Flow (GWF) Process and VSF. The VSF Process allows simulation of flow in unsaturated media above the ground-water zone and facilitates modeling of ground-water/surface-water interactions. Model performance is evaluated by comparison to an analytical solution for one-dimensional (1-D) constant-head infiltration (Dirichlet boundary condition), field experimental data for a 1-D constant-head infiltration, laboratory experimental data for two-dimensional (2-D) constant-flux infiltration (Neumann boundary condition), laboratory experimental data for 2-D transient drainage through a seepage face, and numerical model results (VS2DT) of a 2-D flow-path simulation using realistic surface boundary conditions. A hypothetical 3-D example case also is presented to demonstrate the new capability using periodic boundary conditions (for example, daily precipitation) and varied surface topography over a larger spatial scale (0.133 square kilometer). The new model capabilities retain the modular structure of the MODFLOW code and preserve MODFLOW's existing capabilities as well as compatibility with commercial pre-/post-processors. The overall success of the VSF Process in simulating mixed boundary conditions and variable soil types demonstrates its utility for future hydrologic investigations. This report presents a new flow package implementing the governing equations for variably saturated ground-water flow, four new boundary condition packages unique to unsaturated flow, the Adaptive Time-Stepping Package for use with both the GWF Process and the new VSF Process, detailed descriptions of the input and output files for each package, and six simulation examples verifying model performance.

Techniques and Methods↗
USGS2013

GWM-VI: groundwater management with parallel processing for multiple MODFLOW versions

Groundwater Management–Version Independent (GWM–VI) is a new version of the Groundwater Management Process of MODFLOW. The Groundwater Management Process couples groundwater-flow simulation with a capability to optimize stresses on the simulated aquifer based on an objective function and constraints imposed on stresses and aquifer state. GWM–VI extends prior versions of Groundwater Management in two significant ways—(1) it can be used with any version of MODFLOW that meets certain requirements on input and output, and (2) it is structured to allow parallel processing of the repeated runs of the MODFLOW model that are required to solve the optimization problem. GWM–VI uses the same input structure for files that describe the management problem as that used by prior versions of Groundwater Management. GWM–VI requires only minor changes to the input files used by the MODFLOW model. GWM–VI uses the Joint Universal Parameter IdenTification and Evaluation of Reliability Application Programming Interface (JUPITER-API) to implement both version independence and parallel processing. GWM–VI communicates with the MODFLOW model by manipulating certain input files and interpreting results from the MODFLOW listing file and binary output files. Nearly all capabilities of prior versions of Groundwater Management are available in GWM–VI. GWM–VI has been tested with MODFLOW-2005, MODFLOW-NWT (a Newton formulation for MODFLOW-2005), MF2005-FMP2 (the Farm Process for MODFLOW-2005), SEAWAT, and CFP (Conduit Flow Process for MODFLOW-2005). This report provides sample problems that demonstrate a range of applications of GWM–VI and the directory structure and input information required to use the parallel-processing capability.

Techniques and Methods↗
USGS2012

An assessment of the carbon balance of arctic tundra: comparisons among observations, process models, and atmospheric inversions

Although arctic tundra has been estimated to cover only 8% of the global land surface, the large and potentially labile carbon pools currently stored in tundra soils have the potential for large emissions of carbon (C) under a warming climate. These emissions as radiatively active greenhouse gases in the form of both CO 2 and CH 4 could amplify global warming. Given the potential sensitivity of these ecosystems to climate change and the expectation that the Arctic will experience appreciable warming over the next century, it is important to assess whether responses of C exchange in tundra regions are likely to enhance or mitigate warming. In this study we compared analyses of C exchange of Arctic tundra between 1990–1999 and 2000–2006 among observations, regional and global applications of process-based terrestrial biosphere models, and atmospheric inversion models. Syntheses of the compilation of flux observations and of inversion model results indicate that the annual exchange of CO 2 between arctic tundra and the atmosphere has large uncertainties that cannot be distinguished from neutral balance. The mean estimate from an ensemble of process-based model simulations suggests that arctic tundra acted as a sink for atmospheric CO 2 in recent decades, but based on the uncertainty estimates it cannot be determined with confidence whether these ecosystems represent a weak or a strong sink. Tundra was 0.6 °C warmer in the 2000s compared to the 1990s. The central estimates of the observations, process-based models, and inversion models each identify stronger sinks in the 2000s compared with the 1990s. Similarly, the observations and the applications of regional process-based models suggest that CH 4 emissions from arctic tundra have increased from the 1990s to 2000s. Based on our analyses of the estimates from observations, process-based models, and inversion models, we estimate that arctic tundra was a sink for atmospheric CO 2 of 110 Tg C yr -1 (uncertainty between a sink of 291 Tg C yr -1 and a source of 80 Tg C yr -1 ) and a source of CH 4 to the atmosphere of 19 Tg C yr -1 (uncertainty between sources of 8 and 29 Tg C yr -1 ). The suite of analyses conducted in this study indicate that it is clearly important to reduce uncertainties in the observations, process-based models, and inversions in order to better understand the degree to which Arctic tundra is influencing atmospheric CO 2 and CH 4 concentrations. The reduction of uncertainties can be accomplished through (1) the strategic placement of more CO 2 and CH 4 monitoring stations to reduce uncertainties in inversions, (2) improved observation networks of ground-based measurements of CO 2 and CH 4 exchange to understand exchange in response to disturbance and across gradients of hydrological variability, and (3) the effective transfer of information from enhanced observation networks into process-based models to improve the simulation of CO 2 and CH 4 exchange from arctic tundra to the atmosphere.

Biogeosciences↗
USGS2024

Comparison of sediment and water column nutrient processing rates in agricultural streams of contrasting buffer land use

Watershed nutrient management often focuses on actions that reduce the movement of nitrogen (N) and phosphorus (P) from agricultural lands into streams. One area of management focus is the buffer of land adjacent to streams. Wetlands and forests in this buffer can intercept and retain N and P from the landscape. In addition to directly intercepting agricultural nutrients, natural habitats in the buffer can alter stream geomorphology and influence the in-stream processing and transformation of N and P to less labile and mobile forms. Here, we assess the influence of buffer land cover on in-stream processing of N and P. We measured nutrient dynamics in the water column and sediments of agricultural streams in the Fox River and Duck Creek watersheds (WI, USA) during the growing season. In these streams, water column processing was low, possibly due to a lack of primary producers in the water column. Water column P processing was weakly associated with wetland land cover in the buffer, but buffer land cover had no clear effect on inorganic N processing. On the other hand, sediments were almost always a source of inorganic P and a sink for inorganic N. Sediment P release was higher in streams with more agricultural land cover in the buffer. Sediments in streams with agricultural land cover in the buffer also removed more nitrate, even after accounting for the greater availability of nitrate in those streams. The buffer land cover conditions we quantified occupy a very small portion of the overall watershed (100 m wide, for 1 km upstream of the study site) but nevertheless appear to influence in-stream cycling of N and P. For P management, reducing agricultural land cover in buffers is already a priority due to the ability of wetlands and forests to intercept nutrients, but this study suggests there may be some additional benefit due to changes in in-stream P processing.

Wisconsin↗
USGS2018

Process convolution approaches for modeling interacting trajectories

Gaussian processes are a fundamental statistical tool used in a wide range of applications. In the spatiotemporal setting, several families of covariance functions exist to accommodate a wide variety of dependence structures arising in different applications. These parametric families can be restrictive and are insufficient in some situations. In contrast, process convolutions represent a flexible, interpretable approach to defining the covariance of a Gaussian process and have modest requirements to ensure validity. We introduce a generalization of the process convolution approach that employs multiple convolutions sequentially to form a “process convolution chain”. In our proposed multistage framework, complex dependencies that arise from a combination of different interacting mechanisms are decomposed into a series of interpretable kernel smoothers. We demonstrate an application of process convolution chains to model killer whale movement, in which the paths taken by multiple individuals are not independent but reflect dynamic social interactions within the population. Our proposed model for dependent movement provides inference for the latent dynamic social structure in the study population. Additionally, by leveraging the positive dependence among individual paths, we achieve a reduction in uncertainty for the estimated locations of the killer whales, compared with a model that treats paths as independent.

Environmetrics↗
USGS2024

Processing time and precision of aging structures for Bighead Carp and Silver Carp in the lower Red River catchment in the southern Great Plains

Objective Population demographics of invasive species are commonly evaluated to better develop management actions that are useful for reducing their abundance or controlling the population. Bighead Carp Hypophthalmichthys nobilis and Silver Carp H. molitrix are emblematic invaders in the United States, where they continue to expand their range. There is currently no consensus about which hard structure from these species is best for age estimation. Our study objective was to compare the processing time and precision of five hard structures used for age estimation of both species. Methods We sampled fish in the lower Red River catchment of Oklahoma, Arkansas, and Texas during summer and autumn 2021–2022 and removed the lapillus otolith, left primary pectoral fin ray, postcleithrum, urohyal bone, and anterior-most pterygiophore of the dorsal fin from both Bighead and Silver carp. The structures ( n = 1204) were either embedded in epoxy or thin-sectioned and mounted on slides. Two readers estimated the age of the fish by using each structure and came to a consensus. Processing time was recorded from the onset of laboratory processing to the termination of polishing the cross sections for age estimation. Result Processing of otoliths was comparable to or faster than processing of the other structures and resulted in the highest between-reader agreement. The lowest coefficients of variation in age estimation were represented using lapillus otoliths for Bighead Carp and postcleithra for Silver Carp. Our age bias plots indicated that all other structures underestimated age relative to the lapillus otoliths. Conclusion Our results indicated that using lapillus otoliths for age estimation of these species would have the highest between-reader agreement and would incur no additional laboratory processing time. However, validation is needed to assess whether lapillus otoliths correctly age these species. From a management perspective, use of this structure would facilitate improved population comparisons.

Arkansas, Oklahoma, Texas↗
USGS2016

Processes contributing to resilience of coastal wetlands to sea-level rise

The objectives of this study were to identify processes that contribute to resilience of coastal wetlands subject to rising sea levels and to determine whether the relative contribution of these processes varies across different wetland community types. We assessed the resilience of wetlands to sea-level rise along a transitional gradient from tidal freshwater forested wetland (TFFW) to marsh by measuring processes controlling wetland elevation. We found that, over 5 years of measurement, TFFWs were resilient, although some marginally, and oligohaline marshes exhibited robust resilience to sea-level rise. We identified fundamental differences in how resilience is maintained across wetland community types, which have important implications for management activities that aim to restore or conserve resilient systems. We showed that the relative importance of surface and subsurface processes in controlling wetland surface elevation change differed between TFFWs and oligohaline marshes. The marshes had significantly higher rates of surface accretion than the TFFWs, and in the marshes, surface accretion was the primary contributor to elevation change. In contrast, elevation change in TFFWs was more heavily influenced by subsurface processes, such as root zone expansion or compaction, which played an important role in determining resilience of TFFWs to rising sea level. When root zone contributions were removed statistically from comparisons between relative sea-level rise and surface elevation change, sites that previously had elevation rate deficits showed a surplus. Therefore, assessments of wetland resilience that do not include subsurface processes will likely misjudge vulnerability to sea-level rise.

Ecosystems↗
USGS2024

Measuring and interpreting the surface and shallow subsurface process influences on coastal wetland elevation: A review

A century ago, measuring elevation in tidal wetlands proved difficult, as survey leveling of soft marsh soils relative to a fixed datum was error prone. For 60 years, vertical accretion measures from marker horizons were used as analogs of elevation change. But without a direct measure of elevation, it was not possible to measure the total influence of surface and subsurface processes on elevation. In the 1990s, the surface elevation table (SET) method, which measures the movement of the wetland surface relative to a fixed point beneath the surface (i.e., the SET benchmark base), was combined with the marker horizon method (SET-MH), providing direct, independent, and simultaneous measures of surface accretion and elevation and quantification of surface and shallow subsurface process influences on elevation. SET-MH measures have revealed several fundamental findings about tidal wetland dynamics. First, accretion [ A ] is often a poor analog for elevation change [ E ]. From 50–66% of wetlands experience shallow subsidence ( A > E ), 7–10% shallow expansion ( A < E ), 7% shrink-swell, and for 24–36% A is an analog for E ( A = E ). Second, biological processes within the root zone and physical processes within and below the root zone influence elevation change in addition to surface processes. Third, vegetation plays a key role in wetland vertical dynamics. Plants trap sediment and increase resistance to erosion and compaction. Soil organic matter accumulation can lead to shallow expansion, but reduced plant growth can lead to subsidence, and plant death to soil collapse. Fourth, elevation rates are a better indicator of wetland response to sea-level rise than accretion rates because they incorporate subsurface influences on elevation occurring beneath the marker horizon. Fifth, combining elevation trends with relative sea-level rise (RSLR) trends improves estimates of RSLR at the wetland surface (i.e., RSLR wet ). Lastly, subsurface process influences are fundamental to a wetland’s response to RSLR and plant community dynamics related to wetland transgression, making the SET-MH method an invaluable tool for understanding coastal wetland elevation dynamics.

Estuaries and Coasts↗
USGS2024

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S 2- ) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

California↗
USGS2017

Expanding the role of reactive transport models in critical zone processes

Models test our understanding of processes and can reach beyond the spatial and temporal scales of measurements. Multi-component Reactive Transport Models (RTMs), initially developed more than three decades ago, have been used extensively to explore the interactions of geothermal, hydrologic, geochemical, and geobiological processes in subsurface systems. Driven by extensive data sets now available from intensive measurement efforts, there is a pressing need to couple RTMs with other community models to explore non-linear interactions among the atmosphere, hydrosphere, biosphere, and geosphere. Here we briefly review the history of RTM development, summarize the current state of RTM approaches, and identify new research directions, opportunities, and infrastructure needs to broaden the use of RTMs. In particular, we envision the expanded use of RTMs in advancing process understanding in the Critical Zone, the veneer of the Earth that extends from the top of vegetation to the bottom of groundwater. We argue that, although parsimonious models are essential at larger scales, process-based models offer tools to explore the highly nonlinear coupling that characterizes natural systems. We present seven testable hypotheses that emphasize the unique capabilities of process-based RTMs for (1) elucidating chemical weathering and its physical and biogeochemical drivers; (2) understanding the interactions among roots, micro-organisms, carbon, water, and minerals in the rhizosphere; (3) assessing the effects of heterogeneity across spatial and temporal scales; and (4) integrating the vast quantity of novel data, including “omics” data (genomics, transcriptomics, proteomics, metabolomics), elemental concentration and speciation data, and isotope data into our understanding of complex earth surface systems. With strong support from data-driven sciences, we are now in an exciting era where integration of RTM framework into other community models will facilitate process understanding across disciplines and across scales.

Earth-Science Reviews↗
USGS2023

The dynamic influence of subsurface geological processes on the assembly and diversification of thermophilic microbial communities in continental hydrothermal systems

An accepted paradigm of hydrothermal systems is the process of phase separation, or boiling, of a deep, homogeneous hydrothermal fluid as it ascends through the subsurface resulting in gas rich and gas poor fluids. While phase separation helps to explain first-order patterns in the chemistry and biology of a hot spring’s surficial expression, we know little about the subsurface architecture beneath “phase-separated” pools and the timescales over which phase separation processes occur. Essentially, we have a two-dimensional understanding of a four-dimensional process. By combining geophysical, geochemical, isotopic, and microbiological measurements of two adjacent phase-separated hot springs in Norris Geyser Basin, Yellowstone National Park, we provide a four-dimensional assessment of phase separation processes and their biological manifestation. We uniquely show that Yellowstone’s hydrothermal waters originate from a deep sedimentary aquifer and that both meteoric recharge and shallow reactive transport processes are required to establish the geobiological feedbacks that drive bimodal distributions in the geochemical and microbial composition of hot springs. Specifically, over periods of tens of years, gas-enriched fluids containing volcanic sulfide mix with meteoric waters resulting in microbially-mediated production of sulfuric acid by thermoacidophilic Archaea in the near subsurface. In contrast, over periods of hundreds of years, anoxic residual liquid rises to the surface where it is infused with atmospheric gas fostering Archaea and Bacteria that are largely dependent on oxygen. As such, our results provide formative insight into the causative links between subsurface geological processes, the development of geochemical fluids, and the assembly and diversification of thermophilic microbial communities in hydrothermal systems.

Wyoming↗