Search USGSSearch

SEARCH · Search USGS

Results for “Energy”

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 199 records · Page 11Linked to original sources

Subsurface energy storage and transport for solar-powered geysers on Triton

The location of active geyser-like eruptions and related features close to the current subsolar latitude on Triton suggests a solar energy source for these phenomena. Solidstate greenhouse calculations have shown that sunlight can generate substantially elevated subsurface temperatures. A variety of models for the storage of solar energy in a sub-greenhouse layer and for the supply of gas and energy to a geyser are examined. "Leaky greenhouse" models with only vertical gas transport are inconsistent with the observed upper limit on geyser radius of ∼ 1.5 kilometers. However, lateral transport of energy by gas flow in a porous N 2 layer with a block size on the order of a meter can supply the required amount of gas to a source region ∼1 kilometer in radius. The decline of gas output to steady state may occur over a period comparable with the inferred active geyser lifetime of five Earth years. The required subsurface permeability may be maintained by thermal fracturing of the residual N2 polar cap. A lower limit on geyser source radius of ∼50 to 100 meters predicted by a theory of negatively buoyant jets is not readily attained.

Science

Sierra Nevada plutonic cycle: Part II, tidal energy and a hypothesis for orogenic-epeirogenic periodicities

The dissipative power of the solid earth tides is the order of 10 19 ergs/sec, or a few percent of terrestrial heat flow. It is proposed that this energy is concentrated along oceanic ridge systems and in the asthenosphere by mechanisms of viscous dissipation involving shear melting. Tidal energy localizes and sustains sources of sea-floor spreading through the melting mechanism, convection and magmatic transfer. Components of this energy enter the continent as magmatic heat either where ridge type sources and continents interact or where lateral motions induce shear zones and viscous dissipation within the continent. Temporal maxima of igneous intrusion into continental crust and related epeirogenic oscillations, spaced at intervals of about 30 m.y., are explained in terms of periodic thermal instabilities in the process of shear melting in the mantle. That is, maxima in rates of magma production in the mantle are relieved by vertical magmatic transfer. This process is coupled with lateral motions of the continent in a way analogous to episodic creep episodes of much shorter period in motions of active fault systems. Calculated periodicities are found to be simultaneously compatible with (l) the Sierra Nevada intrusive epochs of Part I, (2) oscillations in the eustatic curve during the Mesozoic Era, (3) concepts of sea-floor spreading, and (4) the magnitude of tidal power. More profound epeirogenic oscillations, having periods of about 200 m.y., are induced by variations in proportioning of tidal energy dissipation between the solid earth and the epicontinental seas. Thus, the tidal deformations of the earth provide information that leads to a general dynamic theory where magmatism, orogency, epeirogeny, sea-floor spreading and continent migration are systematically interrelated.

California

Introduction to this special section: Geothermal energy

Geothermal energy is a global renewable resource that has the potential to provide a significant portion of baseload energy in many regions. In the United States, it has the potential to provide 8.5% of the electric generation capacity by the middle of the century. In general, geothermal systems require heat, permeability, and water to be viable for energy generation. However, with current technologies, only heat is strictly necessary in a native system. Engineered geothermal systems (EGS) introduce water into the subsurface at elevated pressures and reduced temperatures and enhance permeability through hydraulic and/or shear fracturing. Additionally, although moderate- to high-temperature resources currently dominate geothermal energy production, low-temperature resources have been utilized for direct-use cases. When well balanced and maintained, geothermal resources can produce significant amounts of heat and achieve long-term sustainability on the order of an estimated tens to hundreds of years.

The Leading Edge

Satellite psychrometric formulation of the operational simplified surface energy balance (SSEBop) model for quantifying and mapping evapotranspiration

Remote sensing-based evapotranspiration (ET) can be derived using various methods, from soil moisture accounting to vegetation-index based approaches to simple and complex surface energy balance techniques. Due to the complexity of fully representing and parameterizing ET sub-processes, different models tend to diverge in their estimations. However, most models appear to provide reasonable estimations that can meet user requirements for seasonal water use estimation and drought monitoring. One such model is the Operational Simplified Surface Energy Balance (SSEBop). This study presents a formulation of the SSEBop model using the psychrometric principle for vapor pressure/relative humidity measurements where the “dry-bulb” and “wet-bulb” equivalent readings can be obtained from satellite-based land surface temperature estimates. The difference in temperature between the dry (desired location) and wet limit (reference value) is directly correlated to the soil-vegetation composite moisture status (surface humidity) and thus producing a fractional value (0-1) to scale the reference ET. The reference ET is independently calculated using available weather data through the standardized Penman-Monteith equation. Satellite Psychrometric Approach (SPA) explains the SSEBop model more effectively than the energy balance principle because SSEBop does not solve all terms of the surface energy balance such as sensible and ground-heat fluxes. The SPA explanation demonstrates the psychrometric constant for the air can be readily adapted to a comparable constant for the surface, thus allowing the creation of a “surface” psychrometric constant that is unique to a location and day-of-year. This new surface psychrometric constant simplifies the calculation and explanation of satellite-based ET for several applications in agriculture and hydrology. The SPA formulation of SSEBop was found to be an enhancement of the ET equation formulated in 1977 by pioneering researchers. With only two key parameters, improved model results can be obtained using a one-time calibration for any bias correction. The model can be set up quickly for routine monitoring and assessment of ET at landscape scales and beyond.

Arizona

Composition, shell strength, and metabolizable energy of Mulinia lateralis and Ischadium recurvum as food for wintering surf scoters (Melanitta perspicillata)

Decline in surf scoter (Melanitta perspicillata) waterfowl populations wintering in the Chesapeake Bay has been associated with changes in the availability of benthic bivalves. The Bay has become more eutrophic, causing changes in the benthos available to surf scoters. The subsequent decline in oyster beds (Crassostrea virginica) has reduced the hard substrate needed by the hooked mussel (Ischadium recurvum), one of the primary prey items for surf scoters, causing the surf scoter to switch to a more opportune species, the dwarf surfclam (Mulinia lateralis). The composition (macronutrients, minerals, and amino acids), shell strength (N), and metabolizable energy (kJ) of these prey items were quantified to determine the relative foraging values for wintering scoters. Pooled samples of each prey item were analyzed to determine composition. Shell strength (N) was measured using a shell crack compression test. Total collection digestibility trials were conducted on eight captive surf scoters. For the prey size range commonly consumed by surf scoters (6-12 mm for M. lateralis and 18-24 mm for I. recurvum), I. recurvum contained higher ash, protein, lipid, and energy per individual organism than M. lateralis. I. recurvum required significantly greater force to crack the shell relative to M. lateralis. No difference in metabolized energy was observed for these prey items in wintering surf scoters, despite I. recurvum's higher ash content and harder shell than M. lateralis. Therefore, wintering surf scoters were able to obtain the same amount of energy from each prey item, implying that they can sustain themselves if forced to switch prey.

PLoS ONE

Conservation planning for offsetting the impacts of development: a case study of biodiversity and renewable energy in the Mojave Desert

Balancing society’s competing needs of development and conservation requires careful consideration of tradeoffs. Renewable energy development and biodiversity conservation are often considered beneficial environmental goals. The direct footprint and disturbance of renewable energy, however, can displace species’ habitat and negatively impact populations and natural communities if sited without ecological consideration. Offsets have emerged as a potentially useful tool to mitigate residual impacts after trying to avoid, minimize, or restore affected sites. Yet the problem of efficiently designing a set of offset sites becomes increasingly complex where many species or many sites are involved. Spatial conservation prioritization tools are designed to handle this problem, but have seen little application to offset siting and analysis. To address this need we designed an offset siting support tool for the Desert Renewable Energy Conservation Plan (DRECP) of California, and present a case study of hypothetical impacts from solar development in the Western Mojave subsection. We compare two offset scenarios designed to mitigate a hypothetical 15,331 ha derived from proposed utility-scale solar energy development (USSED) projects. The first scenario prioritizes offsets based precisely on impacted features, while the second scenario offsets impacts to maximize biodiversity conservation gains in the region. The two methods only agree on 28% of their prioritized sites and differ in meeting species-specific offset goals. Differences between the two scenarios highlight the importance of clearly specifying choices and priorities for offset siting and mitigation in general. Similarly, the effects of background climate and land use change may lessen the durability or effectiveness of offsets if not considered. Our offset siting support tool was designed specifically for the DRECP area, but with minor code modification could work well in other offset analyses, and could provide continuing support for a potentially innovative mitigation solution to environmental impacts.

California

Differential energy radiation from two earthquakes in Japan with identical M w : The Kyushu 1996 and Tottori 2000 earthquakes

We examine two closely located earthquakes in Japan that had identical moment magnitudes M w but significantly different energy magnitudes M e . We use teleseismic data from the Global Seismograph Network and strong-motion data from the National Research Institute for Earth Science and Disaster Prevention's K-Net to analyze the 19 October 1996 Kyushu earthquake (M w 6.7, M e 6.6) and the 6 October 2000 Tottori earthquake (M w 6.7, M e 7.4). To obtain regional estimates of radiated energy E S we apply a spectral technique to regional (<200 km) waveforms that are dominated by S and Lg waves. For the thrust-fault Kyushu earthquake, we estimate an average regional attenuation Q(f) 230f 0:65 . For the strike-slip Tottori earthquake, the average regional attenuation is Q(f) 180f 0:6 . These attenuation functions are similar to those derived from studies of both California and Japan earthquakes. The regional estimate of E S for the Kyushu earthquake, 3:8 ?? 10 14 J, is significantly smaller than that for the Tottori earthquake, E S 1:3 ?? 10 15 J. These estimates correspond well with the teleseismic estimates of 3:9 ?? 10 14 J and 1:8 ?? 10 15 J, respectively. The apparent stress (T a = ??E s /M 0 with ?? equal to rigidity) for the Kyushu earthquake is 4 times smaller than the apparent stress for the Tottori earthquake. In terms of the fault maturity model, the significantly greater release of energy by the strike-slip Tottori earthquake can be related to strong deformation in an immature intraplate setting. The relatively lower energy release of the thrust-fault Kyushu earthquake can be related to rupture on mature faults at a subduction environment. The consistence between teleseismic and regional estimates of E S is particularly significant as teleseismic data for computing E S are routinely available for all large earthquakes whereas often there are no near-field data.

Bulletin of the Seismological Society of America

Time allocation by Greater White-fronted Geese: Influence of diet, energy reserves and predation

I determined the amount of time Greater White-fronted Geese ( Anser albifrons frontalis ) allocated to various activities from September to May, 1980-1982 at their primary wintering areas in the Pacific Flyway of North America. The length of time spent on roosts during the day was positively correlated to day length. Geese at roost sites spent the majority of their time sleeping (24-46%), alert (17-40%), walking or swimming (6-24%), and in comfort behaviors (3-25%). The amount of time geese fed each day varied little from early autumn to late spring (4.5-4.9 hr), except during mid-winter when minimum temperatures were below freezing (3.9 hr), and immediately before migration in spring (6.3 hr). The proportion of time devoted to feeding and alert behavior, the two most dominant activities at field sites, varied significantly among seasons and locations. The amount of time geese were actively engaged in foraging each season was more dependent on feeding intensity than the amount of time spent at foraging sites (fields), and varied almost three-fold, from 1.8 hr during late winter to 5.1 hr during late spring. Geese fed in closer proximity to conspecifics, were more frequently disturbed, and spent less time feeding during the hunting season. Exploitation of high energy foods and catabolism of substantial energy reserves probably enabled geese to minimize foraging time during periods of harsh weather and high predation pressure. Seasonal variation in the proportion of time spent feeding corresponded closely to changes in body mass. Greater White-fronted Geese wintering in the Pacific Flyway spent substantially less time feeding than they do in Europe, as geese in California fed primarily on high energy cereal grains, while in Europe they subsist on green vegetation which has relatively less digestible energy than cereal grains.

The Condor

Energy budget and prey requirements of breeding lapland longspurs Calcarius lapponicus near Barrow Alaska, U.S.A.

Breeding Lapland longspurs, Calcarius lapponicus , near Barrow, Alaska, have relatively constant energy demands throughout the summer; the average estimated daily energy budgets (DEBs) were 132 and 118 kJ for the male and female. Thermoregulation accounted consistently for one-quarter to one-third of the total DEB. Flight in the male and incubation in the female were major components of the DEB early in the season, whereas cost of molt was a major component for both sexes near the end of the season. Our estimates of longspur DEB based on a time-activity approach were similar to those based on a cage existence model. Minor differences are explained by increased male territorial activity, by energy savings of the female during incubation, and by contraction of the molt for both sexes within the short summer season. Male and female longspurs were estimated to capture 3000 to 10,000 seeds and insects d -1 (3-20 items min -1 foraging) for self maintenance while in summer residence near Barrow. Each adult was estimated to capture an additional 3000 insects d -1 ( 6-7 insects min -1 foraging) during the peak energy requirements to raise five young. While raising young, the maximum required capture rate of prey per time foraging for each adult occurred during the nestling stage; young are just achieving independence, however, when food supply is at a maximum.

Alaska

Energy map of southwestern Wyoming, Part B: oil and gas, oil shale, uranium, and solar

The U.S. Geological Survey (USGS) has compiled Part B of the Energy Map of Southwestern Wyoming for the Wyoming Landscape Conservation Initiative (WLCI). Part B consists of oil and gas, oil shale, uranium, and solar energy resource information in support of the WLCI. The WLCI represents the USGS partnership with other Department of the Interior Bureaus, State and local agencies, industry, academia, and private landowners, all of whom collaborate to maintain healthy landscapes, sustain wildlife, and preserve recreational and grazing uses while developing energy resources in southwestern Wyoming. This product is the second and final part of the Energy Map of Southwestern Wyoming series (also see USGS Data Series 683, http://pubs.usgs.gov/ds/683/ ), and encompasses all of Carbon, Lincoln, Sublette, Sweetwater, and Uinta Counties, as well as areas in Fremont County that are in the Great Divide and Green River Basins.

Wyoming

Energy Resources Program of the U.S. Geological Survey

Our Nation faces the simultaneous challenges of increasing demand for energy, declining domestic production from existing oil and gas fields, and increasing expectations for environmental protection. The Energy Information Administration (2000) forecasts that worldwide energy consumption will increase 32 percent between 1999 and 2020 because of growth of the world economy. Forecasts indicate that in the same time period, U.S. natural gas consumption will increase 62 percent, petroleum consumption will increase 33 percent, and coal consumption will increase 22 percent. The U.S. Geological Survey provides the objective scientific information our society needs for sound decisions regarding land management, environmental quality, and economic, energy, and strategic policy.

Fact Sheet

U.S. Geological Survey energy and wildlife research annual report for 2019 postcard

This postcard provides details about the U.S. Geological Survey (USGS) Energy and Wildlife Research Annual Report for 2019, which highlights new research on the interactions of energy development with wildlife. Encompassing investigations of conventional and renewable energy development across the United States, from the Arctic Coastal Plain of Alaska to the balmy waters of Florida, the report features progress made by USGS scientists and partners in developing methods to minimize the impacts of energy infrastructure on wildlife. The report is available at https://doi.org/10.3133/cir1458 .

General Information Product

A Regional Approach to Wildlife Monitoring Related to Energy Exploration and Development in Colorado

The U.S. Bureau of Land Management (BLM) is currently developing a National Monitoring Strategy that will guide efforts to create an efficient and effective process for monitoring land health by BLM. To inform the ongoing development of the national strategy, BLM selected two States (Colorado, Alaska) to serve as focal areas on which to base a flexible framework for developing monitoring programs that evaluate wildlife responses to energy development. We developed a three-phase monitoring plan to serve as a template and applied it to the design of a monitoring program for the Colorado focal area (White River and Glenwood Springs Field Offices of the BLM). Phase I is a synthesis and assessment of current conditions that capitalizes on existing but under used data sources. A key component is the use of existing habitat and landscape models to evaluate the cumulative effects of surface disturbance. Phase II is the data collection process that uses information provided in Phase I to refine management objectives and provide a linkage to management decisions. The linkage is established through targeted monitoring, adaptive management, and research. Phase III establishes priorities and strategies for regional and national monitoring, and facilitates coordination among other land management agencies and organizations. The three phases are designed to be flexible and complementary. The monitoring plan guides an iterative process that is performed incrementally, beginning with the highest-priority species and management issues, while building on lessons learned and coordination among administrative levels. The activities associated with each phase can be repeated or updated as new information, data, or tools become available. This allows the development of a monitoring program that expands gradually and allows for rapid implementation. A demonstration application of the three-phase monitoring plan was conducted for a study area encompassing five BLM field offices in Colorado: White River, Glenwood Springs, Kremmling, Grand Junction, and Little Snake. The overall study area was selected to encompass the primary distribution of sagebrush and greater sage-grouse (Centrocercus urophasianus) in Colorado and to provide a larger context for evaluating priority management issues of the White River and Glenwood Springs Field Offices. Within the study area, we selected an additional scale of analysis: the Roan Plateau, an area of increased energy development. We focused our demonstration project on the primary objectives identified by the BLM and other stakeholders: (1) How do we evaluate the cumulative effects of energy development and mitigation activities on species of management interest or concern? (2) How can cumulative effects, including direct and indirect effects, be quantified relative to natural variation? (3) How can we implement compliance and effectiveness monitoring to determine whether lease stipulations are met and restoration/mitigation goals are achieved? We targeted greater sage-grouse because this species was identified as a monitoring priority for the White River Field Office. We also evaluated landscape-level indices to address fragmentation resulting from surface disturbance. The primary purpose of this exercise was to highlight the types of analyses and approaches that could be used to evaluate energy development in an application of the three-phase framework. We provide examples of Phase I products that can be used to guide the refinement of management objectives, development of adaptive management and research frameworks, and planning activities in Phase II. Finally, we report on the BLM's National Sage-Grouse Habitat Conservation Strategy as an example of Phase III activities. The development of an effective monitoring program is a daunting task, and will be difficult to implement by the BLM independently due to the current work load of staff. However, the scientific expertise and experience for developing cr

Open-File Report

Geologic cross section, gas desorption, and other data from four wells drilled for Alaska rural energy project, Wainwright, Alaska, coalbed methane project, 2007-2009

Energy costs in rural Alaskan communities are substantial. Diesel fuel, which must be delivered by barge or plane, is used for local power generation in most off-grid communities. In addition to high costs incurred for the purchase and transport of the fuel, the transport, transfer, and storage of fuel products pose significant difficulties in logistically challenging and environmentally sensitive areas. The Alaska Rural Energy Project (AREP) is a collaborative effort between the United States Geological Survey (USGS) and the Bureau of Land Management Alaska State Office along with State, local, and private partners. The project is designed to identify and evaluate shallow (<3,000 ft) subsurface resources such as coalbed methane (CBM) and geothermal in the vicinity of rural Alaskan communities where these resources have the potential to serve as local-use power alternatives. The AREP, in cooperation with the North Slope Borough, the Arctic Slope Regional Corporation, and the Olgoonik Corporation, drilled and tested a 1,613 ft continuous core hole in Wainwright, Alaska, during the summer of 2007 to determine whether CBM represents a viable source of energy for the community. Although numerous gas-bearing coal beds were encountered, most are contained within the zone of permafrost that underlies the area to a depth of approximately 1,000 ft. Because the effective permeability of permafrost is near zero, the chances of producing gas from these beds are highly unlikely. A 7.5-ft-thick gas-bearing coal bed, informally named the Wainwright coal bed, was encountered in the sub-permafrost at a depth of 1,242 ft. Additional drilling and testing conducted during the summers of 2008 and 2009 indicated that the coal bed extended throughout the area outlined by the drill holes, which presently is limited to the access provided by the existing road system. These tests also confirmed the gas content of the coal reservoir within this area. If producible, the Wainwright coal bed contains sufficient gas to serve as a long-term source of energy for the community.

Open-File Report

Oregon OCS seafloor mapping: Selected lease blocks relevant to renewable energy

In 2014 the U.S. Geological Survey (USGS) and the Bureau of Ocean Energy Management (BOEM) entered into Intra-agency agreement M13PG00037 to map an area of the Oregon Outer Continental Shelf (OCS) off of Coos Bay, Oregon, under consideration for development of a floating wind energy farm. The BOEM requires seafloor mapping and site characterization studies in order to evaluate the impact of seafloor and sub-seafloor conditions on the installation, operation, and structural integrity of proposed renewable energy projects, as well as to assess the potential effects of construction and operations on archaeological resources. The mission of the USGS is to provide geologic, topographic, and hydrologic information that contributes to the wise management of the Nation's natural resources and that promotes the health, safety, and well being of the people. This information consists of maps, databases, and descriptions and analyses of the water, energy, and mineral resources, land surface, underlying geologic structure, and dynamic processes of the earth. For the Oregon OCS study, the USGS acquired multibeam echo sounder and seafloor video data surrounding the proposed development site, which is 95 km2 in area and 15 miles offshore from Coos Bay. The development site had been surveyed by Solmar Hydro Inc. in 2013 under a contract with WindFloat Pacific. The USGS subsequently produced a bathymetry digital elevation model and a backscatter intensity grid that were merged with existing data collected by the contractor. The merged grids were published along with visual observations of benthic geo-habitat from the video data in an associated USGS data release (Cochrane and others, 2015). This report includes the results of analysis of the video data conducted by Oregon State University and the geo-habitat interpretation of the multibeam echo sounder (MBES) data conducted by the USGS. MBES data was published in Cochrane and others (2015). Interpretive data associated with this publication is published in Cochrane (2017). All the data is provided as geographic information system (GIS) files that contain both Esri ArcGIS geotiffs or shapefiles. For those who do not own the full suite of Esri GIS and mapping software, the data can be read using Esri ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed August 29, 2016). Web services, which consist of standard implementations of ArcGIS representational state transfer (REST) Service and Open Geospatial Consortium (OGC) GIS web map service (WMS), also are available for all published GIS data. Web services were created using an ArcGIS service definition file, resulting in data layers that are symbolized as shown on the associated report figures. Both the ArcGIS REST Service and OGC WMS Service include all the individual GIS layers. Data layers are bundled together in a map-area web service; however, each layer can be symbolized and accessed individually after the web service is ingested into a desktop application or web map. Web services enable users to download and view data, as well as to easily add data to their own workflows, using any browser-enabled, standalone or mobile device. Though the surficial substrate is dominated by combinations of mud and sand substrate, a diverse assortment of geomorphologic features are related to geologic processes—one anticlinal ridge where bedrock is exposed, a slump and associated scarps, and pockmarks. Pockmarks are seen in the form of fields of small pockmarks, a lineation of large pockmarks with methanogenic carbonates, and areas of large pockmarks that have merged into larger variously shaped depressions. The slump appears to have originated at the pockmark lineation. Video-supervised numerical analysis of the MBES backscatter intensity data and vector ruggedness derived from the MBES bathymetry data was used to produce a substrate model called a seafloor character raster for the study area. The seafloor character raster consists of three substrate classes: soft-flat areas, hard-flat areas, and hard-rugged areas. A Coastal and Marine Ecological Classification Standard (CMECS) geoform and substrate map was also produced using depth, slope, and benthic position index classes to delineate geoform boundaries. Seven geoforms were identified in this process, including ridges, slump scars, slump deposits, basins, and pockmarks. Statistical analysis of the video data for correlations between substrate, depth, and invertebrate assemblages resulted in the identification of seven biomes: three hard-bottom biomes and four softbottom biomes. A similar analysis of vertebrate observations produces a similar set of biomes. The biome between-group dissimilarity was very high or high. Invertebrates alone represent most of the structure of the whole benthic community into different assemblages. A biotope map was generated using the seafloor character raster and the substrate and depth values of the biomes. Hard substrate biotopes were small in size and were located primarily on the ridge and in pockmarks along the pockmark lineation. The soft-bottom bitopes consisted of large contiguous areas delimited by isobaths.

Open-File Report

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan

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

A finite-element simulation model for saturated-unsaturated, fluid-density-dependent ground-water flow with energy transport or chemically- reactive single-species solute transport

SUTRA (Saturated-Unsaturated Transport) is a computer program which can be used to simulate the movement of fluid and the transport of either energy or dissolved substances in a subsurface environment. The model employs a two-dimensional hybrid finite-element and integrated-finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated by SUTRA: (1) fluid density-dependent saturated or unsaturated groundwater flow, and either (2a) transport of a solute in the groundwater, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay, or, (2b) transport of thermal energy in the groundwater and solid matrix of the aquifer. SUTRA provides, as the primary calculated results, fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA may also be used to simulate simpler subsets of the above process. SUTRA may be employed for areal and cross-sectional models of saturated groundwater flow systems, and for cross-sectional models of unsaturated zone flow. Solute transport simulation using SUTRA may be used to simulate natural or man-induced chemical transport, solute sorption, production and decay. SUTRA may be used for simulation of variable density leachate movement, and for cross-sectional simulation of salt-water intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between fresh water and salt water. SUTRA energy transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. (USGS)

Water-Resources Investigations Report