USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth System Science Data”

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,193 records · Page 35Linked to original sources
USGS2025

Cooling perspectives on the risk of pathogenic viruses from thawing permafrost

Climate change is inducing wide-scale permafrost thaw in the Arctic and subarctic, triggering concerns that long-dormant pathogens could reemerge from the thawing ground and initiate epidemics or pandemics. Viruses, as opposed to bacterial pathogens, garner particular interest because outbreaks cannot be controlled with antibiotics, though the effects can be mitigated by vaccines and newer antiviral drugs. To evaluate the potential hazards posed by viral pathogens emerging from thawing permafrost, we review information from a diverse range of disciplines. This includes efforts to recover infectious virus from human remains, studies on disease occurrence in polar animal populations, investigations into viral persistence and infectivity in permafrost, and assessments of human exposure to the enormous viral diversity present in the environment. Based on currently available knowledge, we conclude that the risk posed by viruses from thawing permafrost is no greater than viruses in other environments such as temperate soils and aquatic systems.

mSystems↗
USGS2025

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS1988

Geology of an area near Brentwood, Williamson County, Tennessee

The lithology and geologic structure of an area near Brentwood, Williamson County, Tennessee, were studied to define the potential aquifers and confining units of the ground-water flow system. Four formations were identified. They are, in descending order, the Bigby-Cannon Limestone, the Hermitage Formation, the Carters Limestone, and the Lebanon Limestone. The Bigby-Cannon Limestone and the Carters Limestone are potential aquifers. The Hermitage Formation and the Lebanon Limestone are confining units. The Bigby-Cannon Limestone and the Hermitage Formation have been affected by recent erosion. Any variation of the Carters Limestone is controlled by pre-Carters erosion of the top of the Lebanon Limestone. The thickness of the Carters Limestone ranges from 65 to 79 feet. Structurally, the area reflects the regional northwest dip. A small scale anti-cline-syncline pair also is evident. This feature is not a result of erosion and also occurs in the T-3 bentonite bed in the Carters Limestone.

Tennessee↗
USGS2026

Production of mineral commodities and geospatial map of the mineral industries and related infrastructure of China

As part of the U.S. Geological Survey’s (USGS) mission to distribute global mineral information and analyze supply chains, this study provides a comprehensive review of the global significance of China’s mineral production and capacity in 2023. Of 77 mineral commodities in the USGS dataset, China produced 74 and was the world’s first-ranked producer for 39 of the 74. Compared to the high share of global mineral production, including up to 98 percent of global gallium production, the country’s share of global mineral reserves was relatively small, ranging from 20 percent (zinc ore) to 52 percent (tungsten ore). China’s imports of metal ores, slag, and ash accounted for 64 percent of global imports of such commodities by value. The country’s exports of base metals and articles of base metal accounted for 17 percent of the global exports. To help nongeographic information system users assess the spatial distribution of mineral mines, processing facilities, and ports for trades in China, this study created a geospatial (also called “georeferenced”) portable document format (GeoPDF) map. In addition, the GeoPDF contains mineral resource tracts (such as antimony, copper, potash, coal, and oil and gas), exploration sites, and energy infrastructure based on the preexisting USGS data.

Open-File Report↗
USGS2024

Interactive effects of salinity and hydrology on radial growth of bald cypress (Taxodium distichum (L.) Rich.) in coastal Louisiana, USA

Tidal freshwater forests are usually located at or above the level of mean high water. Some Louisiana coastal forests are below mean high water, especially bald cypress ( Taxodium distichum (L.) Rich.) forests because flooding has increased due to the combined effects of global sea level rise and local subsidence. In addition, constructed channels from the coast inland act as conduits for saltwater. As a result, saltwater intrusion affects the productivity of Louisiana’s coastal bald cypress forests. To study the long-term effects of hydrology and salinity on the health of these systems, we fitted dendrometer bands on selected trees to record basal area increment as a measure of growth in permanent forest productivity plots established within six bald cypress stands. Three stands were in freshwater sites with low salinity rooting zone groundwater (0.1–1.3 ppt), while the other three had higher salinity rooting zone groundwater (0.2–4.9 ppt). Water level was logged continuously, and salinity was measured monthly to quarterly on the surface and in groundwater wells. Higher groundwater salinity levels were related to decreased bald cypress radial growth, while higher freshwater flooding increased radial growth. With these data, coastal managers can model rates of bald cypress forest change as a function of salinity and flooding.

Louisiana↗
USGS1996

Hydrogeologic investigation and simulation of ground-water flow in the Upper Floridan Aquifer of north-central Florida and southwestern Georgia and delineation of contributing areas for selected city of Tallahassee, Florida, water-supply wells

A 4-year investigation of the Upper Floridan aquifer and ground-water flow system in Leon County, Florida, and surrounding counties of north-central Florida and southwestern Georgia began in 1990. The purpose of the investigation was to describe the ground-water flow system and to delineate the contributing areas to selected City of Tallahassee, Florida, water-supply wells. The investigation was prompted by the detection of low levels of tetrachloroethylene in ground-water samples collected from several of the city's water-supply wells. Hydrologic data and previous studies indicate that; ground-water flow within the Upper Floridan aquifer can be considered steady-state; the Upper Floridan aquifer is a single water-bearing unit; recharge is from precipitation; and that discharge occurs as spring flow, leakage to rivers, leakage to the Gulf of Mexico, and pumpage. Measured transmissivities of the aquifer ranged from 1,300 ft 2 /d (feet squared per day) to 1,300,000 ft2/d. Steady-state ground-water flow in the Upper Floridan aquifer was simulated using a three-dimensional ground- water flow model. Transmissivities ranging from less than 5,000 ft2/d to greater than 11,000,000 ft 2 /d were required to calibrate to observed conditions. Recharge rates used in the model ranged from 18.0 inches per year in areas where the aquifer was unconfined to less than 2 inches per year in broad areas where the aquifer was confined. Contributing areas to five Tallahassee water-supply wells were simulated by particle- tracking techniques. Particles were seeded in model cells containing pumping wells then tracked backwards in time toward recharge areas. The contributing area for each well was simulated twice, once assuming a porosity of 25 percent and once assuming a porosity of 5 percent. A porosity of 25 percent is considered a reasonable average value for the Upper Floridan aquifer; the 5 percent porosity simulated the movement of ground-water through only solution-enhanced bedding plains and fractures. The contributing areas were generally elliptical in shape, reflecting the influence of the sloping potentiometric surface. The contributing areas delineated for a 5 percent porosity were always much larger than those determined using a 25 percent porosity. The lowest average ground-water velocity computed within a contributing area, using a 25 percent porosity, was 1.0 ft/d (foot per day) and the highest velocity was 1.6 ft/d. The lowest average ground-water velocity, determined using a 5 percent porosity, was 2.4 ft/d and the highest was 7.4 ft/d. The contributing areas for each of the five wells was also determined analytically and compared to the model-derived areas. The upgradient width of the simulated contributing areas were larger than the upgradient width of the analytically determined contributing areas for four of the five wells. The model could more accurately delineate contributing areas because of the ability to simulate wells as partially penetrating and by incorporating complex, three-dimensional aquifer characteristics, which the analytical method could not.

Florida, Georgia↗
USGS2025

Fiber-optic distributed temperature sensing of hydrologic processes—Diverse deployments and new applications by the U.S. Geological Survey

Fiber-optic distributed temperature sensing instruments harness the temperature-dependent properties of glass to measure temperature continuously along optical fibers by using precise pulses of laser light. In the mid-2000s, this technology was refined for environmental monitoring purposes such as snowpack-air exchange, groundwater/surface-water exchange, and lake-water stratification. Fiber-optic distributed temperature sensing has revealed unprecedented details about preferential flow processes; however, the method is labor intensive and requires specific training, resulting in limited use by the broader water community. With the ongoing national implementation of the U.S. Geological Survey Next Generation Water Observing System, there has been renewed interest in harnessing the unique spatiotemporal monitoring capabilities of fiber-optic distributed temperature sensing. This fact sheet briefly describes this technology, highlights uses by the U.S. Geological Survey, and discusses current applications and future opportunities.

Fact Sheet↗
USGS2025

International data gaps at the Center for Engineering Strong Motion Data

The Center for Engineering Strong Motion Data (CESMD) is utilized by seismologists, engineers, and disaster management professionals in the US and has historically achieved and distributed waveforms from across the globe for significant earthquakes. The increased access to the waveforms via Web API (Application Programming Interface) offers a unique opportunity to provide the community complete datasets, sampling a variety of tectonic environments and geologic conditions, increasing the number of available ground motion records for use in ground motion models (GMMs) and improving the accuracy of earthquake engineering evaluations. The objective of this study is to programmatically identify gaps in global event data from the past decade and backfill missing data gaps at CESMD. We first compare the CESMD catalog with the Advanced National Seismic System (ANSS) Comprehensive Earthquake Catalog identifying regions and time periods where strong-motion data is limited or inadequate. To backfill datasets at CESMD for significant events, we pinpoint regions and time intervals that lack information, creating a list of events for which we’d like to obtain data. An important facet of this work is identifying the source of data and metadata across earthquake repositories around the world and integrating these data repositories into our current strong-motion data processing workflow. In parallel with these newly processed datasets, we are developing a script to produce data origination citations to include provenance and attribution information to associate with respective datasets at CESMD. We showcase our methodology for identifying and filling data gaps at CESMD using three case studies (the 2018 Anchorage Alaska earthquake sequence, seismicity associated with the 2018 Hawaiian Kilauea volcano eruption, and several earthquakes in Turkey) and then outline our strategy to apply our data gap backfilling methods on an international scale.

Conference Paper↗
USGS2025

Lithium from magma to mine in an early Yellowstone hotspot caldera

Renewable energy technologies rely on the extraction of metals not historically in high demand, such as lithium (Li), for which ore deposit models are incompletely understood. One of the world’s largest Li deposits is hosted in lake sediments of the 16.4 Ma McDermitt caldera, which formed during the early stages of Yellowstone hotspot volcanism in the western United States. Eruptive and posteruptive mobility of Li are major challenges in elucidating deposit formation. Melt inclusions preserved in quartz crystals provide a means to assess pre-eruptive magmatic Li contents. Concentrations of Li determined by ion microprobe for melt inclusions in a McDermitt rhyolite lava are 400−1350 ppm, compared to 20−70 ppm Li in matrix rhyolite glasses. Synthesis with melt inclusion data for eight additional calderas demonstrates a recurrence of Li-rich rhyolitic magmas (200−2000 ppm Li) in the western part of the Yellowstone hotspot track. However, unlike the multicyclic caldera complexes with overlapping fault networks that may have compromised Li retention, the McDermitt caldera remained a closed hydrologic system throughout its evolution. Modeling indicates 100 km 3 of resurgent magma could yield 25−150 Mt Li in a magmatic fluid and supports accumulation of Li-rich magmatic fluid in a closed intracaldera lake, followed by evaporative concentration and sequestration of Li within clay minerals to generate the McDermitt deposit.

Idaho, Nevada, Oregon, Wyoming↗
USGS2024

New developments at the Center for Engineering Strong-Motion Data (CESMD)

The Center for Engineering Strong-Motion Data (CESMD), an internationally utilized joint center of the U.S. Geological Survey (USGS) and the California Geological Survey (CGS), provides a single access point for earthquake strong-motion records and station metadata from the CGS California Strong-Motion Instrumentation Program (CSMIP), the USGS National Strong-Motion Project (NSMP), the USGS Advanced National Seismic System, and other affiliates. The CESMD has been continuously improving its webtools to facilitate the access of strong-motion data and metadata for use in post-earthquake response and for scientific and engineering research applications. The Center provides raw and processed strong-motion data via the Engineering Data Center (EDC) and the Virtual Data Center (VDC) web portals. This paper focuses on the strong-motion products provided by the EDC where more than 48,000 records with peak ground accelerations greater than 0.1% g from over 2400 earthquakes are currently hosted. and on the ongoing efforts to develop data access tools and applications. The new developments and ongoing efforts in the EDC include: 1) enhancements to the CESMD webservices to facilitate access to station metadata, earthquake information, and strong motion records 2) new features to the interactive map interface, improving the visualization and access to earthquake, station, and record information, 3) efforts to develop a new web application tool for data format conversion from a number of data formats, 4) efforts to unify varying waveform data formats into a consistent format, 5) ongoing efforts to compile seismic station site geology, measured or inferred Vs30 values, shear-wave profiles, NEHRP site class, and available structural instrument deployment schematics, and 6) a special studies pages for research topic-specific ground motion datasets that offer uniform processing of records from a variety of sources.

Conference Paper↗
USGS2026

Selected water-quality data from the Cedar River and Cedar Rapids well fields, Cedar Rapids, Iowa, 2017–22

The Cedar River alluvial aquifer is the source of drinking water in Cedar Rapids, Iowa. Production wells are completed in the alluvial aquifer approximately 40 to 80 feet below land surface. The City of Cedar Rapids and the U.S. Geological Survey have studied the groundwater-flow system and water quality of the aquifer in the vicinity of Cedar Rapids since 1992. Results of these studies documented hydrologic conditions, water quality, and geochemistry of the alluvial aquifer and interactions with the Cedar River. Water-quality samples were collected for studies involving well field monitoring, trends, source-water protection, groundwater geochemistry, surface-water–groundwater interaction, and pesticides in groundwater and surface water. Water quality was analyzed for dissolved major ions (boron, bromide, calcium, chloride, fluoride, iron, magnesium, manganese, potassium, silica, sodium, sulfate, and total dissolved solids), dissolved nutrients (ammonia as nitrogen, ammonia plus organic nitrogen as nitrogen, nitrite plus nitrate as nitrogen, nitrite as nitrogen, orthophosphate as phosphorus, and phosphorus), dissolved organic carbon, and selected pesticides. Physical characteristics (alkalinity, dissolved oxygen, pH, specific conductance, and water temperature) were measured on site and recorded for each water sample collected. This report presents the results of routine water-quality data-collection activities from October 2017 through September 2022. Methods of data collection, quality assurance, water-quality analyses, and statistical procedures are presented. Data include the results of water-quality analyses from quarterly sampling from monitoring wells, production wells, two water treatment plants, and the Cedar River at Blairs Ferry Road at Palo, Iowa, streamgage (U.S. Geological Survey station number 05464420), as well as monthly nutrient sampling from the Cedar River and Morgan Creek near Covington, Iowa, streamgage (U.S. Geological Survey station number 05464475).

Iowa↗
USGS2024

Evaluating the impact of uncertainty in ground motion forecasts for post-earthquake impact modeling applications

The US Geological Survey’s (USGS) ShakeMap system provides a rapid characterization of strong ground shaking in areas directly affected by an earthquake. This study focuses on studying the aggregate effects of macroseismic shaking estimates from ShakeMap, expressed in terms of modified Mercalli intensity (MMI), when accounting for the uncertainty in forecasted ground motions. We use a Monte Carlo approach to generate numerous spatially correlated realizations of ground motions by utilizing a combination of circulant embedding and kriging techniques for efficiently handling the correlations. We then assessed the aggregate effects of shaking by looking at bin counts across these realizations. We demonstrate that the aggregate shaking regarding the mean macroseismic intensity estimates (from the ShakeMap output) is a biased representation of the aggregate shaking when shaking uncertainty is included. Incorporating shaking uncertainty can help to improve various downstream earthquake impact applications, such as the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) overall earthquake fatality distribution or estimates of shaking-induced ground failure impacts from consequential earthquakes.

Earthquake Spectra Journal↗
USGS2024

Developing international standards and guidelines for disseminating and validating simulated ground motions

We are leading an effort to develop international standards and guidelines for curating, disseminating, and validating simulated ground-motion data. This effort is organized as a working group within the Consortium of Organizations for Strong Motion Observation Systems (COSMOS). In 2022 we held online workshops on 7 and 8 June and 20 October. The first workshop focused on curating and disseminating simulated ground-motion data. The second workshop focused on validating simulated ground-motions for engineering applications. About 100 people participated in each of the workshops with strong representation from North America, Europe, and western Asia. In the coming year, we intend to form a technical committee to draft the international guidelines and standards while continuing to engage stakeholders from across the globe through online and in-person meetings. The key points from the workshops include: (1) Numerous groups are generating simulated earthquake ground motions and making them openly available; however there is very little coordination among groups to provide consistent interfaces for searching and retrieving data; (2) Standardizing interfaces for metadata and data access should consider agile approaches that can adapt to changing capabilities and user needs while building upon existing efforts; (3) A distributed architecture is preferred by most participants to allow institutions to host and manage their own data while broadcasting their holdings to a combined catalog; (4) Validation of ground-motion simulations applies to the entire workflow for simulating earthquake ground-motions, including the rupture model, seismic velocity model, and seismic wave propagation software; (5) Metrics to evaluate the validation are application dependent; more research is needed to tie structural response characteristics to ground-motion characteristics; and (6) Validation results should provide a clear, transparent, and quantitative assessment of the simulated ground motions.

Conference Paper↗
USGS2024

Geospatial PDF map of the compilation of GIS data for the mineral industries of select countries in the Indo-Pacific region

Introduction In 2024, the U.S. Geological Survey's (USGS) National Minerals Information Center (NMIC) completed the project titled "Compilation of geospatial data for the mineral industries of select countries in the Indo-Pacific." This project aimed to leverage the expertise and capabilities of the NMIC to collect, synthesize, and interpret geospatial data to inform on the extractive resources of select countries in the Indo-Pacific region (area of study) and expand the NMIC's understanding on the impact of mineral industry of these countries in the global economy. The 19 countries of interest in the Indo-Pacific study area include Bangladesh, Bhutan, Brunei, Burma, Fiji, Malaysia, Mongolia, Nauru, New Caledonia, New Zealand, Papua New Guinea, Philippines, Singapore, Solomon Islands, South Korea (Republic of Korea), Sri Lanka, Taiwan, Timor-Leste, and Vietnam. The primary objective of this effort was to create a fully attributed Geographic Information System (GIS) portraying existing mining infrastructure, resources, and production capacities across the Indo-Pacific study area as well as highlight mineral production and processing sites under development and potential areas of future extractive industry operations and development in the region. The compiled GIS geodatabase with supporting documentation including comprehensive metadata was published as a USGS data release titled "Compilation of Geospatial Data (GIS) for the Mineral Industries of Select Countries in the Indo-Pacific." This georeferenced portable document format (GeoPDF) map sheet presents a new geographic information product containing a partial representation of the GIS data. This GeoPDF map provides a visual comparison of the distribution of mineral industry GIS data, which contributes to a deeper understanding of the intersections and complexities of the extractive industries within the select countries in the Indo-Pacific region.

Open-File Report↗
USGS2026

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

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗
USGS2024

SUTRA— A code for simulation of saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of the version 4.0 enhancements—Freeze-thaw capability, saturation and relative-permeability relations, spatially varying properties, and enhanced budget and velocity outputs

Version 4.0 of the Saturated-Unsaturated Transport (SUTRA) software code provides the capability to simulate the freezing and thawing of groundwater during energy transport simulations under saturated and unsaturated conditions. In addition to the types of hydrogeologic processes that SUTRA has been able to simulate in the past, this version can be used to study the effects of the freeze-thaw process on the flow and energy dynamics of hydrogeologic systems. The freeze-thaw simulation capability accounts for the latent heat of fusion and allows thermal property values to vary with changing total-water saturation, liquid-water saturation, and ice saturation. It allows the effective permeability of the porous medium to change as a result of freezing and thawing. This version also provides several user-selectable relations for the dependence of total-water saturation on fluid pressure, the dependence of liquid-water saturation on temperature during freezing and thawing, and the dependence of relative permeability on liquid saturation, as well as three user-selectable formulae for defining the bulk thermal conductivity of a mixture of solid grains, liquid water, ice, and air. For unsaturated simulations without freezing, the selectable total-water saturation relations eliminate the need for the user to program these and their associated relative-permeability functions, as had been required in previous SUTRA versions. Optional nonlinear dependence of fluid density on temperature, which covers the range from supercooled (about −50 degrees Celsius) to superheated (about 400 degrees Celsius), is also provided. Additionally, this version makes it possible to spatially vary parameters that, in previous versions of SUTRA, were required to be spatially uniform: solid-matrix properties, adsorption parameters, and parameters for production of solute mass or energy. Spatial variation is also allowed for the newly included freeze-thaw process parameters. Additional enhancements provide (1) output of water-mass and energy budgets that include values of all component terms in the governing balance equations, and (2) output of Darcy velocities (fluid fluxes), in addition to the velocity output provided by previous SUTRA versions. These enhanced outputs allow fuller interpretation of simulation results, especially for freeze-thaw phenomena. The set of processes simulated by this version of SUTRA are useful for studying a wide range of hydrogeologic system types, conditions, and questions. For cryohydrogeologic simulations, however, this version of the code is limited in that (1) it does not simulate thermomechanical effects of freeze-thaw, (2) pressure changes due to water density change during freezing are neglected, (3) ice saturation cannot exceed the initial porosity of the simulated medium, and (4) cryosuction, the migration of liquid water toward freezing fronts, is neglected. Furthermore, this version does not account for air flow or for water vaporization and sublimation under unsaturated conditions.

Techniques and Methods↗