USGS Science⌕ Search

SEARCH · USGS Science

Results for “Renewable 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.

129 records · Page 2Linked to original sources
USGS2025

AAPG Energy and Minerals Division Tight Oil and Gas Committee Activities and Commodity Report for 2021-2022: Tuscaloosa Marine Shale, Gulf Coast basin, Louisiana and Mississippi

The Upper Cretaceous Tuscaloosa marine shale (TMS) potential production area encompasses 20.4 million acres across central Louisiana (LA), southern Mississippi (MS), southwestern Alabama (AL), and a small southwestern section of the Florida panhandle (Hackley et al., 2018). It remains a minor and largely undeveloped unconventional shale oil play with production from the TMS confined along the east-west LA-MS State boundary (Fig. 1). The mean undiscovered, technically recoverable resources in the TMS are estimated at 1.5 billion barrels of oil and 4.6 trillion cubic feet of gas (Hackley et al., 2018). Geochemical analyses of source rock solvent extracts and oil samples indicate that, in the play area, the TMS is the primary source of shale oil produced from the TMS as well as of accumulated oils in the underlying conventional reservoirs of the lower Tuscaloosa (Hackley et al., 2020). The API gravity of TMS oils ranges from approximately 34 to 46 degrees (Hackley et al., 2020; Croke et al., 2020).

AAPG Wiki↗
USGS2025

Evaluating the potential of co-located non-petroleum energy resources in the Gulf Coast using play fairwayaAnalysis

Geological resources critical to the energy transition, such as sedimentary geothermal, carbon storage potential, and lithium in brines, commonly struggle for economic feasibility as stand-alone developments but can have greater viability when the potential for more than one of these resources exist within the same reservoir or stacked in different stratigraphic intervals. There may also be instances where development of one resource inhibits development of others and decisions need to be made on how to best prioritize the use of those resources in the subsurface. Subsurface data sets were analyzed and integrated to evaluate the distribution of non-petroleum energy and related resources in the onshore and nearshore U.S. Gulf Coast. Temperature, pressure, brine composition (lithium content), and reservoir quality data for thirty-four depositional units have been compiled and visualized to high-grade areas where multiple resource opportunities likely coexist. For sedimentary geothermal, possible resource areas are defined as low potential (<90°C), moderate potential (90–150°C), and high potential (>150°C). For CO 2 storage, high potential areas exhibit supercritical CO 2 conditions less than 80% of the fracture gradient. Areas with pore pressure between 80% and 100% of the fracture gradient are considered to exhibit moderate potential and areas where the fracture gradient is equal to or greater than 100% are assigned low potential. Lithium resources in brines were defined by lithium concentrations as low potential (<100ppm), moderate potential (100-200ppm), or high potential (>200ppm). Reservoir quality affects the viability of all three of these resources and is evaluated using depositional environment maps of each unit. The resulting play fairway maps can be used for regional scale screening evaluations of these resources and to identify areas of interest where more detailed, prospect-scale studies can be undertaken.

Gulf Coast↗
USGS2025

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal↗
USGS2025

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2026

Neutron scattering reveals fractionation of natural gas mixtures in unconventional petroleum reservoir pores: Perspectives on energy resource recovery and storage

In unconventional petroleum reservoirs hydrocarbon fluids are hosted by both mineral and organic matter pores. These pores can have diameters that range from microns to less than a single nanometer and, for unconventional reservoirs, there is evidence that small pores ( <20 nm diameter) may constitute a large proportion of the available space. Understanding subsurface volumes and how fluids behave in them can be helpful for predicting hydrocarbon production and storage in the subsurface. One area with knowledge gaps regarding hydrocarbon behavior in small pores is the possibility for mixtures to fractionate (i.e., unmix) based on pore size or pore type. Mixture fractionation as a function of pore size could impact recovery of hydrocarbons, drive compositional shifts during production, and limit fluid storage within candidate reservoirs. To investigate natural gas fractionation in small geologic pores, we applied total neutron scattering to probe methane-ethane mixtures at reservoir pressures (up to ≈30 MPa) and temperature (60°C) within a sample from the Upper Cretaceous Niobrara Formation. Neutron scattering data reveal only minor fractionation occurs between methane and ethane in 20-nm diameter sample mesopores. Increased fractionation is observed for sample micropores, with up to 72% (±1% at 1-sigma) methane found in 2 nm diameter pores following injection of a 50%-50% methane-ethane mixture. These data provide rarely available direct experimental observations of hydrocarbon mixture behavior under nanoconfinement in a sample from an important unconventional petroleum reservoir. Our results are discussed in the context of evaluating hydrocarbon resources in unconventional reservoir meso- and micropores, reconciling observed gas composition changes during production, and more broadly, understanding subsurface pore volumes within an energy storage framework.

Fuel↗
USGS2025

The 1912 Ms 7.2 earthquake in the Denali region of central Alaska

The 2002 M w 7.9 earthquake in central Alaska ruptured 340 km on three faults—Susitna Glacier thrust fault, Denali fault, Totschunda fault—crossing both the Richardson Highway and the Alaska Pipeline. Its occurrence prompted renewed interest in historical large earthquakes that possibly originated on the Denali fault. One of these earthquakes was a M s 7.2 event on July 7, 1912, which we revisit with two approaches: (1) probabilistic relocation of the epicenter using globally recorded arrival times, and (2) compilation and reassessment of shaking intensity reports to estimate a macroseismic epicenter. Our preferred instrumental epicenter is west of the Parks Highway and in agreement with the maximum‐reported shaking, which was from the Parker–Browne expedition of Denali. We also relocated a M s 6.4 aftershock, whose epicenter is 11 km from the mainshock. Candidate faults for the 1912 earthquake include the Denali fault, the McLeod Creek thrust fault, and the Kantishna Hills thrust fault. Future analysis of active faults, paleoseismic results, 1912 instrumental data, and 1912 felt reports may help in interpreting the fault and mechanism of the 1912 earthquake.

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

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗
USGS2025

Earthquake stress drop, source spectra from laboratory experiments, and the apparent unsteadiness of frictional slip

Over a range of co-seismic total energy release, high sample rate slip, stress, and slip velocity measurements were made during frictional failure to determine seismic source properties. The near-field measurements were used to estimate stress drops using spectral approaches that are analogous to the standard seismological moment-corner determined Δσ and acceleration amplitude Δσ arms stress drops. Rapid frictional sliding produces slip velocity amplitude spectra that vary with frequency -1 so long as the energy dissipated by friction approaches the total energy release. Under these conditions Δσ and Δσ arms are similar in size. Since these are nearly direct measurements of fault slip velocity and because they produce flat acceleration spectra between the corner and an apparent fmax, friction in these experiments seems to produce dynamic acceleration and deceleration at all resolved frequencies. The physical origin of this unsteady motion is not known exactly but likely reflects rapid dynamic variations in fault strength or applied stress. On the other hand, when the total energy stored prior to failure is artificially made somewhat larger relative to that which can be dissipated by friction the velocity spectra vary with frequency -2 immediately below the corner. At these conditions Δσ and Δσ arms diverge. The excess stored energy is preferentially partitioned into radiated energy upon failure. Collectively the experiments imply that to produce unsteady, white noise accelerations that are observed for natural earthquakes, requires that at least 95% of the energy released does not escape the source region to be radiated to the far-field.

Bulletin of the Seismological Society of America↗
USGS2024

Trends and environmental impacts of virtual water trade

Virtual water describes water embedded in the production of goods and offers meaningful insights about the complex interplay between water, trade and sustainability. In this Review, we examine the trends, major players, traded products and key drivers of virtual water trade (VWT). Roughly 20% of water used in global food production is traded virtually rather than domestically consumed. As such, agriculture dominates VWT, with livestock products, wheat, maize, soybean, oil palm, coffee and cocoa contributing over 70% of total VWT. These products are also driving VWT growth, the volume of which has increased 2.9 times from 1986 to 2022. However, the countries leading VWT contributions (with China, the United States, the Netherlands, Germany and India accounting for 34% of the global VWT in 2022) have remained relatively stable over time, albeit with China becoming an increasingly important importer. VWT can mitigate the effects of water scarcity and food insecurity, although there are concerns about the disconnect between consumers and the environmental impacts of their choices, and unsustainable resource exploitation. Indeed, approximately 16% of unsustainable water use and 11% of global groundwater depletion are virtually traded. Future VWT analyses must consider factors such as water renewability, water quality, climate change impacts and socioeconomic implications.

Nature Reviews Earth & Environment↗
USGS2025

Reduction of red bed sedimentary rocks in connection with energy metal ore formation: A case study from the Sinbad seep, Mesa County, Colorado

The Paradox Basin’s Sinbad seep is a modern analog for ancient bleaching of red bed sediments by introduced alkaline, reducing brines. This bleaching, involving reductive alteration of former red beds, is essential ground preparation that enables the altered rocks to trap Cu, U, and V from later oxidized fluids, forming ore deposits. Study of Sinbad thus offers insights into these metallic mineralization processes in the Paradox and other sedimentary basins. The Sinbad seep occurs where shallow groundwater interacts with organic-rich petroleum source rocks, flows up a fault, and discharges into Salt Creek. Ratios of Na/Cl, Na/Br, and Cl/Br, plus high sulfate concentrations, indicate that the salinity of the seep water originated from dissolution of halite and gypsum during topographically driven flow of meteoric water across a diapir of the Pennsylvanian Paradox Formation, rather than from deep basinal brines. In contact with the organic-rich shales of the Paradox, these SO 4 -rich waters become reduced through bacterial reduction, producing H 2 S. The waters react with red beds of the Permian Cutler Formation, causing pervasive bleaching. The bleaching at Sinbad is characterized by iron-conservative reduction of ferric iron in diagenetic hematite and detrital ilmenite and magnetite to Fe sulfides. The ferrous iron is retained in these sulfides likely as sorbed Fe. Associated alteration includes precipitation of quartz and feldspar overgrowths, partial dissolution of detrital quartz and feldspar caused by pressure solution, formation of clay and authigenic rutile, and precipitation of carbonate and gypsum. The Fe sulfides rapidly degenerate to a mixture of jarosite and iron oxides on weathering at the surface. Among the major-oxide elements and most trace elements, there is no statistical difference between unbleached and bleached samples. The exceptions are uranium and sulfur, which are somewhat greater in bleached samples. Fission track radiography illuminates that uranium is preferentially concentrated in iron oxide cements, iron sorbed into detrital clasts, and finely disseminated iron oxide in illite cement in bleached samples. A leaching experiment suggests that uranium may be more easily available for mobilization from rocks that have undergone bleaching alteration, making them potential U sources. In addition, bleached rocks form effective traps for U, V, and some Cu mineralization.

Colorado↗
USGS2000

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida↗
USGS2025

Estimated ultimate recovery (EUR) Prediction for Eagle Ford Shale using integrated datasets and artificial neural networks

The estimated ultimate recovery (EUR) is an important parameter for forecasting oil and gas production and informing decisions regarding field development strategies. In this study, we combined site-specific geologic, completion, and operational parameters with the predictive capabilities of machine learning (ML) models to predict EURs of the wells for the Eagle Ford Marl Continuous Oil Assessment Unit. We developed an extensive dataset of wells that have produced from the lower and upper Eagle Ford Shale intervals and reduced the model complexity using principal component analysis. We tested the ML models and estimated the sensitivities of ML-predicted EURs to changes in the values of different input variables. The results of applying the optimized ML model to the Eagle Ford suggest that the approach developed in this study could be promising. The ML estimates of the EURs fit the DCA-based values with an R 2 ~ 0.9 and a mean absolute error of ~36 × 10 3 bbl. In the lower Eagle Ford Shale, the EUR estimates were found to be most sensitive to changes in porosity, net thickness of the interval, clay volume, and the API gravity of the oil; and that in the upper Eagle Ford Shale they were most sensitive to changes in the total organic carbon and water saturation, which suggests that it could be important to consider these parameters in assessing these intervals or close analogs.

Louisiana, Mississippi, Texas↗
USGS2026

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado↗
USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2025

Spatiotemporal variations in strain release and seismic rupture in multifault systems: An example from Panamint Valley, southeastern California

Geometrically complex, multifault ruptures have been observed in recent, damaging earthquakes in southeastern California, sparking renewed efforts to identify physical conditions that promote or inhibit fault discontinuity-spanning coseismic ruptures. The likelihood of ruptures propagating across fault discontinuities is thought to be partly controlled by fault geometries, rupture direction, and the history of strain release. However, these parameters vary in space and time over multiple earthquake cycles, making it difficult to forecast the likelihood that an earthquake on one fault will trigger rupture on a nearby fault. Here we use tectono-geomorphic mapping of a geometrically complex fault zone in Panamint Valley, southeastern California, to assess spatiotemporal variations of paleo-rupture patterns and geometries of fault discontinuities over multiple earthquake cycles. First, we identify ten generations of late Pleistocene to Holocene alluvium using geomorphic parameters and luminescence dating to constrain ages of alluvium and bracket late Holocene earthquake timing. Then, we quantify slip kinematics using high-resolution structure from motion digital surface models. We find the Panamint Valley transtensional relay (PVTR) hosted four late Holocene earthquakes, bracketed to ~5.8–3.4 ka, ~3.8–2.2 ka, ~2.4–0.6 ka, and ~0.64–0.16 ka, with ~0.6–1.1 m of slip per event, correlative to Mw ≈ 6.7–6.9 earthquakes. Additionally, we find similarities in earthquake timing on the Ash Hill, PVTR, and Panamint Valley faults and similarities in the slip magnitude and slip kinematics between the Ash Hill and PVTR faults, implying that the PVTR may co-rupture with nearby faults. Paleo-rupture patterns indicate that seismogenic strain transfer may occur through the PVTR, along different combinations of fault segments and jump distances, over multiple earthquake cycles. These data highlight the utility of tectono-geomorphic mapping in evaluating paleo-rupture patterns and suggest that the PVTR may act to propagate and/or arrest rupture between the Ash Hill and Panamint Valley faults.

California↗
USGS1976

An environmental assessment of impacts of coal development on the water resources of the Yampa River basin, Colorado and Wyoming; Phase I work plan

Coal resources of the western United States are being developed at ever-increasing rates, causing concerns of the effects of mining and associated activities on the environment. The Yampa River basin in northwestern Colorado and south-central Wyoming is undergoing economic development of its coal, oil and gas, and uranium resources. The Yampa River basin assessment is a 21-year program of the U.S. Geological Survey. It is designed primarily to assess the availability and quality of the basin's water resources. The basin assessment also will evaluate potential environmental and selected socioeconomic impacts of energy-resource development plans proposed by mining and power companies. This report serves as a project work plan for the basin assessment's first-phase work activities. Economic development in the Yampa River basin will result in higher levels of residuals (noneconomic byproducts) being discharged into the environment. Residuals-management encompasses the technological and economic interrelationships that should be considered in order to comply with Federal and State environmental-control regulations regarding air quality, water quality, and disposal of solid wastes. The basin assessment, the first phase that is described by this report, is designed to provide Federal, State, and local decision-makers with the basic environmental information needed to formulate and to evaluate policies for the development of energy resources. The techniques applied and conclusions reached in the Yampa River basin assessment should aid similar studies of other energy-rich basins in the western United States.

Colorado, Wyoming↗