Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “The Energy Journal”

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 757 records · Page 42Linked to original sources
USGS2020

A two-stage step-wise framework for fast optimization of well placement in coalbed methane reservoirs

Coalbed methane (CBM) has emerged as a clean energy resource in the global energy mix, especially in countries such as Australia, China, India and the USA. The economical and successful development of CBM requires a thorough evaluation and optimization of well placement prior to field-scale exploitation. This paper presents a two-stage, step-wise optimization framework to obtain the optimal placement of wells for large-scale development of CBM reservoirs. In the first stage, an optimal uniform well pattern is obtained by optimizing well pattern description parameters with the particle swarm optimization (PSO) algorithm. Subsequently, the location and status (active/inactive) of each well are perturbed and optimized within the patterns through the integration of the generalized pattern search (GPS) algorithm and a quality map (QM) representing the production potential. This framework was tested in a synthetic anthracite CBM reservoir in the Qinshui basin (with high gas content and low permeability) and a real field high volatile bituminous reservoir in the Illinois basin (with low gas content and high permeability). The results show that: (i) significant variations in the net present value (NPV) exist with respect to different uniform well patterns (even for cases where the total number of wells are identical), the optima of which can be efficiently determined by the PSO within 100 numerical simulation runs; (ii) the optimization of well perturbations by the GPS results in a more noticeable improvement in NPVs for the synthetic (12.3%) than for the real field model (4.6%); (iii) for the low permeable synthetic model with narrow optimal well spacings (320 m × 200 m), the contribution of the optimization of well perturbation to the NPV increment is heavily dependent on the uniform well placement solution; (iv) for the high permeable real field model with large optimal well spacings (1300 m × 1300 m), the initial uniform well placement has a very minor effect on the subsequent well perturbation solutions in terms of NPV; (v) the proposed framework significantly outperforms the conventional well-by-well concatenation procedure in terms of computational efficiency, robustness and optimal criteria set for production potential.

International Journal of Coal Geology↗
USGS1981

Economics and coal resource appraisal: strippable coal in the Illinois Basin ( USA)

Coal-resource appraisals generally describe the location and general characteristics of coal beds. Estimates are made of the average overburden depth (depth of the coal bed below the surface), bed thickness, and perhaps certain chemical properties of the coal [1]. Although such resource compilations represent an important initial step, neither they nor current estimates provide sufficient information to determine the costs of alternative options for National energy policy. Because coal is expected to provide an increasing part of future overall U.S. energy supply, it is crucial for long term planning that coalresource appraisals convey sufficient information regarding the degree of economic resource scarcity (escalation of unit costs for mining remaining deposits as the best deposits in an area are mined out) expected as coal consumption increases. However, assumptions embodied in most large-scale models of coal supply [13; 20] imply that the coal reserves that exist can be commercially produced as needed. As demonstrated by the embargo of 1973-1974 and more recent difficulties in obtaining crude oil, even a slight temporary commodity shortfall will result in significant economic losses when the commodity supplies a large proportion of the total energy used by an economy. For the United States, crude oil resource estimates did not give warning of a decline in domestic petroleum product. In this paper we argue that coal-resource estimates, as they are now made, will not give warning of future supply difficulties. A method for incorporating an economic dimension into appraisals of strippable coal resources is presented and is applied to a major U.S. coal-producing region, the Illinois part of the Illinois basin. Illinois accounts for nearly 70% of the demonstrated strippable coal reserve base of the Illinois basin [12]. In particular, a long-run incremental cost function (that is unit costs vs. cumulative reserves extracted) is estimated for strippable coal in Illinois. The estimated cost function exhibits an initial range over which costs increase at a constant rate-followed by a range where costs increase very rapidly and the function becomes quite inelastic. This long-run incremental cost function is significant for two reasons. Comparison of the demonstrated reserve base for strippable coal (14.8 billion tons) with cumulative production (since 1920 of 1.1 billion tons) might suggest that future depletion will not be significant when, in fact, the cost function presented here indicates otherwise. Because the only states having more reserves of strippable coal than Illinois are Montana and Wyoming [18], these results have implications for national coal-resource appraisal. Secondly, most models used to study the effects of alternative public policies in coal markets are static in nature [15; 16] and cannot consider the effects of increasing costs. Furthermore, models which are multi-period in nature [11; 13] appear to have little empirical basis for their assumptions regarding the escalation of production costs resulting from depletion. The plan of the paper is as follows. First a description of the analytical approach for constructing the long-run incremental cost-reserve function is discussed. Following this, the descriptions of the Illinois basin and the basic physical data are presented. In the concluding section, the main results are presented and their implications are explored

Southern Economic Journal↗
USGS2024

The post-emergence period for denning polar bears: Phenology and influence on cub survival

Among polar bears ( Ursus maritimus ), only parturient females den for extended periods, emerging from maternal dens in spring after having substantially depleted their energy reserves during a fast that can exceed 8 months. Although den emergence coincides with a period of increasing prey availability, polar bears typically do not depart immediately to hunt, but instead remain at the den for up to a month. This delay suggests that there are likely adaptive advantages to remaining at the den between emergence and departure, but the influence of the timing and duration of this post-emergence period on cub survival has not been evaluated previously. We used temperature and location data from 70 denning bears collared within the Southern Beaufort Sea and Chukchi Sea subpopulations to estimate the phenology of the post-emergence period. We evaluated the influence of various spatial and temporal features on duration of the post-emergence period and evaluated the potential influence of post-emergence duration on litter survival early in the spring following denning. For dens that likely contained viable cubs at emergence ( n = 56), mean den emergence occurred on 16 March (SE = 1.4 days) and mean departure on 24 March (SE = 1.6 days), with dates typically occurring later in the Chukchi Sea relative to Southern Beaufort Sea and on land relative to sea ice. Mean duration of the post-emergence period was 7.9 days (SE = 1.4) for bears that were observed with cubs later in the spring, which was over 4 times longer than duration of those observed without cubs (1.9 days). Litter survival in the spring following denning ( n = 31 dens) increased from 0.5 to 0.9 when duration of the post-emergence period increased by ~4 days and other variables were held at mean values. Our limited sample size and inability to verify cub presence at emergence suggests that future research is merited to improve our understanding of this relationship. Nonetheless, our results highlight the importance of the post-emergence period in contributing to reproductive success and can assist managers in developing conservation and mitigation strategies in denning areas, which will be increasingly important as human activities expand in the Arctic.

Journal of Mammalogy↗
USGS2024

Geese migrating over the Pacific Ocean select altitudes coinciding with offshore wind turbine blades

Renewable energy facilities are a key part of mitigating climate change, but can pose threats to wild birds and bats, most often through collisions with infrastructure. Understanding collision risk and the factors affecting it can help minimize impacts on wild populations. For wind turbines, flight altitude is a major factor influencing collision risk, and altitude-selection analyses can evaluate when and why animals fly at certain altitudes under certain conditions. We used GPS tags to track Pacific Flyway geese (Pacific greater white-fronted goose, tule greater white-fronted goose and lesser snow goose) on transoceanic migrations between Alaska and the Pacific Coast of the contiguous United States, an area where offshore windfarm development is beginning. We evaluated how geographic and environmental covariates affected (1) whether birds were at rest on the water versus in flight (binomial model) and (2) altitude selection when in flight (similar to a step-selection framework). We then used a Monte Carlo simulation to predict the probability of flying at each altitude under various conditions, considering both the fly/rest decision and altitude selection. In both spring and fall, geese showed strong selection for altitudes within the expected rotor-swept zone (20–200 m asl), with 56% of locations expected to be within the rotor-swept zone under mean daylight conditions and 28% at night. This indicates a high possibility that migrating geese may be at risk of collision when passing through windfarms. Although there was some variation across subspecies, geese were most likely to be within the rotor-swept zone with little wind or light tailwinds, low clouds, little to no precipitation and moderate to cool air temperatures. Geese were unlikely to be in the rotor-swept zone at night, when most individuals were at rest on the water. Synthesis and applications . These results could be used to inform windfarm management, including decisions to shut down turbines when collision risk is high. The altitude-selection framework we demonstrate could facilitate further study of other bird species to develop a holistic view of how windfarms in this area could affect the migratory bird community as a whole.

Journal of Applied Ecology↗
USGS2025

Component identification of solid biomass fuels using reflected light microscopy: Interlaboratory study 2

As nations transition toward sustainable energy systems , biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

International Journal of Coal Geology↗
USGS2023

Interlaboratory study: Testing reproducibility of solid biofuels component identification using reflected light microscopy

Considering global market trends and concerns about climate change and sustainability, increased biomass use for energy is expected to continue. As more diverse materials are being utilized to manufacture solid biomass fuels, it is critical to implement quality assessment methods to analyze these fuels thoroughly. One such method is reflected light microscopy (RLM), which has the potential to complement and enhance current standard testing, leading to improving fuel quality assessment and, ultimately, preventing avoidable air pollution. An interlaboratory study (ILS) was conducted to test the reproducibility of biomass fuels component identification using a reflected light microscopy technique. The exercise was conducted on thirty photomicrographs showing biomass and various undesired components (like plastics or mineral matter), which were purposely added (by the ILS organizers) to contaminate wood pellets and charcoal-based grilling fuels. Forty-six participants had various levels of difficulty identifying the marked components, and as a result, the percentage of correct answers ranged from 52.2 to 94.4%. Among the most difficult components to distinguish were petroleum products and inorganic matter. Various reasons led to the misidentification, including insufficient morphological descriptions of the components provided to participants, ambiguities of the nomenclature, limitations of the analytical and exercise method, and insufficient experience of the participants. Overall, the results indicate that RLM has the potential to enhance the quality assessment of biomass fuels. However, they also demonstrate that the petrographic classification used in this exercise requires further refinement before it can be standardized. While a new simplified classification of solid biomass fuels components was created as an outcome of this study, future research is necessary to refine the nomenclature, develop a microscopic morphological description of the components, and verify the accuracy of component identification with a follow-up ILS.

International Journal of Coal Geology↗
USGS2026

VegET evapotranspiration for Africa: Continental-scale simulation, multi-product evaluation, and drought assessment

Study region Continental Africa, encompassing diverse climatic zones—tropical, arid, and temperate—and spanning major transboundary river basins such as the Nile, Niger, Congo, Volta, and Zambezi River Basins. The region exhibits pronounced hydroclimatic gradients and heterogeneous land use systems ranging from rainfed croplands and rangelands to dense tropical forests and irrigated schemes. Study focus Actual evapotranspiration (ETa) is a central component of the terrestrial water balance, governing the redistribution of water and energy between the land surface and the atmosphere. Accurate estimation of ETa at continental scale is critical for hydrological monitoring, water resource management, and climate adaptation, as well as for quantifying water, energy, and carbon fluxes that underpin sustainable development. In this study, we applied the agro-hydrologic VegET v2 model to simulate a new, high-resolution, continental-scale ETa dataset for Africa (2000–2021). The model results were benchmarked against four widely used remote sensing-based products—MODIS16 v6.1, SSEBop v6.1, WaPOR v3, and GLEAM v4.1a—across major climate zones, land use types, and River Basins, providing a comprehensive multi-product evaluation of evapotranspiration dynamics across the continent. New hydrological insights for the region Validation against eddy covariance flux tower observations at eight representative sites confirmed that VegET v2 accurately reproduces the seasonal dynamics of observed ETa, achieving a correlation (r) of 0.8 and an RMSE of 25 mm month⁻¹ —accuracy that is comparable to or higher than accuracies of satellite-based products MODIS16, SSEBop, and GLEAM. This study represents one of the first Africa-wide hydrological simulations of ETa, extending the VegET model beyond basin-scale applications. Intercomparisons reveal that VegET aligns closely with MODIS16, SSEBop, and GLEAM in humid and tropical regions (r = 0.80–0.90; RMSE < 20 mm month⁻¹), while greater discrepancies appear in arid and semi-arid zones, where WaPOR tends to overestimate ETa (RMSE ≥ 28 mm month⁻¹). Despite these differences, VegET effectively captures spatial and temporal ETa variability across rainfed croplands, forests, and savannas, supporting its utility in regional water balance assessments, water accounting, and drought monitoring. A key application of VegET v2 is the Evapotranspiration Deficit Index (ETDI), derived by integrating VegET-based ETa with potential evapotranspiration (PET) to quantify water stress. ETDI successfully captured major drought episodes across Africa, including persistent Sahelian and southern African dry spells, the 2020–2021 winter drought in the Maghreb, and the 2018–2019 austral summer drought in southern Africa, while identifying positive anomalies over central Africa indicative of recurrent wetness. These results underscore VegET’s capability as a hydrologically consistent, operational tool for continental ETa monitoring and drought assessment, offering support for basin-scale water balance studies, food security planning, and climate resilience across Africa’s diverse hydrological environments.

Journal of Hydrology: Regional Studies↗
USGS2019

Using a process-based model of pre-eruptive seismic patterns to forecast evolving eruptive styles at Sinabung Volcano, Indonesia

Most volcanoes worldwide are not monitored in real-time; for those that are, patterns of pre-eruptive earthquakes coupled with conceptual models of magma ascent enable short-term forecasting of eruption onset. Basic event locations, characterization of background seismicity, and recognition of changes in earthquake types and energy release are most important to successful eruption forecasting. During renewed activity at Sinabung volcano, Indonesia, this approach was used by the Center for Volcanology and Geological Hazards Mitigation (CVGHM) and the USGS Volcano Disaster Assistance Program to forecast eruption onset, identify changes in eruptive styles and raise or lower alert levels and extend or contract evacuation zones. After > 400 years of quiescence, Sinabung began erupting in August 2010. The volcano was unmonitored at the onset of these eruptions, which were phreatic, but soon after a monitoring network was installed by CVGHM. Increasing swarms of high-frequency volcano tectonic (VT) earthquakes were used to forecast continuing phreatic eruptions. Volcanic activity decreased in mid-September 2010, while additional intrusions at depth (inferred from continued distal VT swarms) continued through September 2013, when explosive phreatic eruptions recurred. Explosive eruptions were forecast based on increases in real-time seismic amplitude measurement (RSAM) and VT seismicity. Seismicity changed markedly in late November and early December 2013 with the occurrence of deep earthquakes and an overall transition from low-frequency (LF) dominated and irregular (in time and magnitude) earthquakes to more regular LF and hybrid seismicity – a transition that accompanied the continued rise, eventual emergence and growth of a lava dome in the summit crater. This lava dome was first observed on 18 December. In late December 2013 to early January 2014, the eruptive style changed again as additional ascending magma deformed the summit and the dome grew beyond the capacity of the summit crater, resulting in the en masse collapse of the lava dome (2 Mm3) on 11 January and the largest pyroclastic flow to date. The collapse was forecast on the basis of a several order of magnitude increase in RSAM, continued strong distal VT seismicity, an increase in proximal seismicity, and large-scale observed deformation of the summit area. Similarly, a later collapse of a second summit lava dome on 1 February 2014 was forecast on the basis of increased distal seismicity. Here, we demonstrate how a process-based volcano seismicity model was used in combination with real-time data to forecast the time and magnitude of eruptions, as well as changes in eruption style.

Sinabung Volcano↗
USGS2020

Diurnal timing of nonmigratory movement by birds: The importance of foraging spatial scales

Timing of activity can reveal an organism's efforts to optimize foraging either by minimizing energy loss through passive movement or by maximizing energetic gain through foraging. Here, we assess whether signals of either of these strategies are detectable in the timing of activity of daily, local movements by birds. We compare the similarities of timing of movement activity among species using six temporal variables: start of activity relative to sunrise, end of activity relative to sunset, relative speed at midday, number of movement bouts, bout duration, and proportion of active daytime hours. We test for the influence of flight mode and foraging habitat on the timing of movement activity across avian guilds. We used 64570 days of GPS movement data collected between 2002 and 2019 for local (non‐migratory) movements of 991 birds from 49 species, representing 14 orders. Dissimilarity among daily activity patterns was best explained by flight mode. Terrestrial soaring birds began activity later and stopped activity earlier than pelagic soaring or flapping birds. Broad‐scale foraging habitat explained less of the clustering patterns because of divergent timing of active periods of pelagic surface and diving foragers. Among pelagic birds, surface foragers were active throughout the day while diving foragers matched their active hours more closely to daylight hours. Pelagic surface foragers also had the greatest daily foraging distances, which was consistent with their daytime activity patterns. This study demonstrates that flight mode and foraging habitat influence temporal patterns of daily movement activity of birds.

Journal of Avian Biology↗
USGS2020

Evaluation of a roughness length parametrization accounting for wind–wave alignment in a coupled atmosphere–wave model

The importance of wind energy as an alternative energy source has increased over the latest years with more focus on offshore winds. A good estimation of the offshore winds is thus of major importance for this industry. Up to now the effect of the wind–wave (mis)alignment has not yet been taken into account in coupled atmosphere–wave models to study the vertical wind profile and power production estimations of offshore wind farms. In this study the roughness length parametrization of Drennan et al . in 2003, and its extension addressing the wind–wave (mis)alignment proposed by Porchetta et al . in 2019, are investigated in the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) model. This study shows that the yearly mean wind estimation at hub height (100 m) is improved by the roughness length parametrization of Porchetta et al . compared to Drennan. This is mainly due to the increased roughness of the former parametrization compare to the latter, even in aligned wind–wave conditions. This difference in roughness is caused by the dataset used to obtain the constants, deep‐water conditions versus mixed offshore conditions. Moreover, the roughness length parametrization of Porchetta et al . performs better in two of three alignment categories. Furthermore, similar model performances are obtained if we exclude the wind directions from the wind shadow zone of the measurement mast or the wind directions from the recently built Alpha Ventus wind farm, which is in close vicinity of the measurement mast. Investigating different wind conditions shows that the new roughness length parametrization of Porchetta et al . performs best for both offshore and onshore winds. Additionally, we show that the coupled model estimations of the vertical wind are only slightly affected by significant wave height estimations. Similar model performances for different accuracies of significant wave height estimations are presented. One exception is the perpendicular alignment category where the new roughness length of Porchetta et al . outperforms the roughness length of Drennan when investigating the wind estimations related to significant wave heights with a higher accuracy. The roughness length parametrization of Porchetta et al . reduced the power production overestimation of the coupled model from 5.7 to 2.8%. We also show that the standalone atmospheric model including the roughness length of Charnock in 1955 has a degraded performance compared to the coupled model including the roughness length parametrization of Porchetta et al . for yearly average wind profiles.

Quarterly Journal of the Royal Meteorological Soci↗
USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗
USGS1991

High-energy carbonate-sand accumulation, the Quicksands, southwest Florida Keys

High-resolution seismic-reflection profiles of the Quicksands, located along a broad ridge on the platform shelf west of Key West, Florida, indicate a significant deposit of non-oolitic carbonate sand occurs in a belt 47 km long by 28 km wide. The surface of the belt is ornamented by large (5 m), migrating tidal bars, oriented in a north-south direction, on which sand waves, oriented in an east-west direction, are superimposed. Some of the sand waves are awash at low tide. The sand waves are formed by strong reversing tidal currents flowing between the Gulf of Mexico and the Straits of Florida. The waves migrate directly over Pleistocene bedrock to the east, but the deposit thickens to the west and sand waves there overlie non-oolitic Holocene accumulations as thick as 12 m. Westward-dipping accretionary bedding indicates that net migration of the sands is to the west, despite north-south movement of tidal currents. The westward edge of the accumulation has accreted over deeper, muddier deposits. Although tidal currents and resultant bedforms appear identical to those of active ooid deposits in the Bahamas and elsewhere, no oolitically coated grains were found in this study. Thin-section analyses show the principal component (average 48%) of the sands is fragmented plates of species of the green alga Halimeda , followed by particulate coral (average 17%), which increases off the flanks of the main sand body. Short vibracores confirm the presence of cross-bedding.

Florida↗
USGS2011

Acetate availability and its influence on sustainable bioremediation of Uranium-contaminated groundwater

Field biostimulation experiments at the U.S. Department of Energy's Integrated Field Research Challenge (IFRC) site in Rifle, Colorado, have demonstrated that uranium concentrations in groundwater can be decreased to levels below the U.S. Environmental Protection Agency's (EPA) drinking water standard (0.126 μM). During successive summer experiments – referred to as “Winchester” (2007) and “Big Rusty” (2008) - acetate was added to the aquifer to stimulate the activity of indigenous dissimilatory metal-reducing bacteria capable of reductively immobilizing uranium. The two experiments differed in the length of injection (31 vs. 110 days), the maximum concentration of acetate (5 vs. 30 mM), and the extent to which iron reduction (“Winchester”) or sulfate reduction (“Big Rusty”) was the predominant metabolic process. In both cases, rapid removal of U(VI) from groundwater occurred at calcium concentrations (6 mM) and carbonate alkalinities (8 meq/L) where Ca-UO 2 -CO 3 ternary complexes constitute >90% of uranyl species in groundwater. Complete consumption of acetate and increased alkalinity (>30 meq/L) accompanying the onset of sulfate reduction corresponded to temporary increases in U(VI); however, by increasing acetate concentrations in excess of available sulfate (10 mM), low U(VI) concentrations (0.1–0.05 μM) were achieved for extended periods of time (>140 days). Uniform delivery of acetate during “Big Rusty” was impeded due to decreases in injection well permeability, likely resulting from biomass accumulation and carbonate and sulfide mineral precipitation. Such decreases were not observed during the short-duration “Winchester” experiment. Terminal restriction fragment length polymorphism (TRFLP) analysis of 16S rRNA genes demonstrated that Geobacter sp. and Geobacter -like strains dominated the groundwater community profile during iron reduction, with 13 C stable isotope probing (SIP) results confirming these strains were actively utilizing acetate to replicate their genome during the period of optimal U(VI) removal. Gene transcript levels during “Big Rusty” were quantified for Geobacter -specific citrate synthase ( gltA ), with ongoing transcription during sulfate reduction indicating that members of the Geobacteraceae were still active and likely contributing to U(VI) removal. The persistence of reducible Fe(III) in sediments recovered from an area of prolonged (110-day) sulfate reduction is consistent with this conclusion. These results indicate that acetate availability and its ability to sustain the activity of iron- and uranyl-respiring Geobacter strains during sulfate reduction exerts a primary control on optimized U(VI) removal from groundwater at the Rifle IFRC site over extended time scales (>50 days).

Geomicrobiology Journal↗
USGS2005

Anatahan, Northern Mariana Islands: Reconnaissance geological observations during and after the volcanic crisis of spring 1990, and monitoring prior to the May 2003 eruption

Anatahan island is 9.5 km east–west by 3.5 km north–south and truncated by an elongate caldera 5 km east–west by 2.5 km north–south. A steep-walled pit crater ∼1 km across and ∼200 m deep occupies the eastern part of the caldera. The island is the summit region of a mostly submarine stratovolcano. The oldest subaerial rocks (stage 1) are exposed low on the outer flanks and in the caldera walls. These include thick (∼10 m) and thin (2–3 m) lava flows, well-indurated tuffs, and scoria units that make up the bulk of the island. Rock compositions range from basaltic andesite to dacite, and most are plagioclase-phyric. On the steep north and south flanks of the volcano, these rocks are cut by numerous east–west-oriented, few-hundred-m-long lineaments of undetermined origin. Indurated breccias unconformably overlie scarps cut into stage 1 units low on the south flank. Intermediate-age eruptive units (stage 2) include caldera-filling lava flows and pyroclastic deposits and, on the outer flanks, vents and valley-filling lava flows. The youngest pre-2003 volcanic unit on Anatahan (stage 3) is a hydromagmatic surge and fall deposit rich in accretionary lapilli. Prior to 2003, this unit was found over almost the entire island, and in many places original depositional surfaces and outcrops could be found in high-energy environments along the coast, indicating a young (but undetermined) age. During reconnaissance visits in 1990, 1992, 1994, and 2001, geothermal activity (fumaroles as well as pits with boiling, sediment-laden pools) was observed in the southern part of the pit crater. In March and April 1990, increased local seismicity, a large regional earthquake, and reported increased fumarolic activity in the pit crater prompted evacuation of Anatahan village, at the west end of the island. Our first field investigation took place in late April 1990 to assess the level of volcanic unrest, conduct reconnaissance geological observations, collect rock and geothermal water samples, and set up a geophysical monitoring network. Results at this time were inconclusive with respect to determining whether the activity was anomalous. Water in some of the geothermal pits within the pit crater was boiling, and pH values as low as 0.7 were recorded in the field. An electronic distance measurement (EDM) network was installed, and over a ∼1-week period, up to 9 cm of extension occurred across some lines but not others. Seismicity was characterized by intermittent local earthquakes but no sustained swarms or tremor. A brief visit in June 1990 revealed that the shallow lake near the boiling pits was gone, but activity in the pits themselves was similar to that of April 1990. Only minor extension had occurred along a single EDM line since the previous visit, and no earthquakes >M2.5 occurred during the visit. Subsequent 1- to 2-day visits occurred in October 1990, May 1992, May 1994, and June 2001. Activity within the geothermal pits was relatively constant during every visit, although during this 11-year period the level of the water in each pit decreased. In June 2001, a ∼50-m-wide region of mud pots and steaming ground in the central part of the geothermal area had developed. No geologic evidence, however, suggested that an eruption would occur < 2 years afterward. Most of the EDM lines showed slight extension between late 1990 and 1992, followed by very gradual contraction from 1992 to 2001. A more extensive seismic-monitoring system was installed on the Northern Mariana Islands during these visits, and it recorded a small seismic swarm at Anatahan from May to July 1993. The telemetry component of the seismic equipment broke prior to 2001 and had not been repaired by the time of the May 2003 eruption, so no precursory seismic data were recorded to indicate pre-eruption unrest.

Northern Mariana Islands↗
USGS2003

The saturated zone at Yucca Mountain: An overview of the characterization and assessment of the saturated zone as a barrier to potential radionuclide migration

The US Department of Energy is pursuing Yucca Mountain, Nevada, for the development of a geologic repository for the disposal of spent nuclear fuel and high-level radioactive waste, if the repository is able to meet applicable radiation protection standards established by the US Nuclear Regulatory Commission and the US Environmental Protection Agency (EPA). Effective performance of such a repository would rely on a number of natural and engineered barriers to isolate radioactive waste from the accessible environment. Groundwater beneath Yucca Mountain is the primary medium through which most radionuclides might move away from the potential repository. The saturated zone (SZ) system is expected to act as a natural barrier to this possible movement of radionuclides both by delaying their transport and by reducing their concentration before they reach the accessible environment. Information obtained from Yucca Mountain Site Characterization Project activities is used to estimate groundwater flow rates through the site-scale SZ flow and transport model area and to constrain general conceptual models of groundwater flow in the site-scale area. The site-scale conceptual model is a synthesis of what is known about flow and transport processes at the scale required for total system performance assessment of the site. This knowledge builds on and is consistent with knowledge that has accumulated at the regional scale but is more detailed because more data are available at the site-scale level. The mathematical basis of the site-scale model and the associated numerical approaches are designed to assist in quantifying the uncertainty in the permeability of rocks in the geologic framework model and to represent accurately the flow and transport processes included in the site-scale conceptual model. Confidence in the results of the mathematical model was obtained by comparing calculated to observed hydraulic heads, estimated to measured permeabilities, and lateral flow rates calculated by the site-scale model to those calculated by the regional-scale flow model. In addition, it was confirmed that the flow paths leaving the region of the potential repository are consistent with those inferred from gradients of measured head and those independently inferred from water-chemistry data. The general approach of the site-scale SZ flow and transport model analysis is to calculate unit breakthrough curves for radionuclides at the interface between the SZ and the biosphere using the three-dimensional site-scale SZ flow and transport model. Uncertainties are explicitly incorporated into the site-scale SZ flow and transport abstractions through key parameters and conceptual models. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2026

Effects of a low-lipid diet on the gut microbiome and head kidney transcriptome of juvenile Chinook Salmon

Objective Pacific salmon Oncorhynchus spp. reared in production hatcheries are typically fed high-lipid, energy-dense diets to achieve large size and high body condition prior to release. In contrast, juveniles in natural environments tend to consume low-lipid, high-protein diets, and fish reared for research or conservation purposes are sometimes fed diets that are formulated to mimic natural diets and promote wild-like phenotypes. Understanding how these alternative diets affect fish health beyond growth and body condition could ultimately contribute to improving hatchery fish fitness. Methods In this work, we evaluated changes in the fecal microbiome and gene expression of juvenile Chinook Salmon O. tshawytscha on a standard high-lipid hatchery diet versus a low-lipid diet formulated to mimic the nutrition profile of natural-origin fish. To evaluate the time scale at which diet alters the fecal microbiome, we collected longitudinal samples over a 12-week period and switched the diets of a subset of fish twice during the experiment. We used 16S ribosomal RNA gene amplicon sequencing to characterize fecal microbiome differences between fish on the two diets as well as hatchery-reared fish at a production hatchery, hatchery fish that had been captured after release into a stream, and natural-origin, stream-reared fish of similar ages. Additionally, we conducted RNA sequencing on head kidney samples from laboratory-reared fish to evaluate changes in gene expression in this important immune organ. Results We found that the low-lipid diet and the hatchery diet resulted in microbiomes that differed from the microbiome of natural-origin fish and from each other and that diet-driven changes to the microbiome could occur in under 14 d. The low-lipid diet did not result in a microbiome that resembled the microbiome of naturally produced fish. Instead, the low-lipid diet resulted in a microbiome community that was distinct from those of fish reared on the hatchery diet and fish sampled from the wild. The RNA sequencing results indicated differential enrichment of pathways related to immunity, metabolism, and hormone synthesis between fish that were fed the two experimental diets. Conclusions The results suggest that additional environmental factors influence the microbiome more strongly than diet formulation or that the low-lipid diet has a smaller effect on the microbiome than a natural, ­invertebrate-based diet. Given that the gut microbiome and systemic immune function contribute significantly to disease resistance, our findings highlight the importance of understanding how diets fed to fish in captivity may affect fish health beyond growth and body condition metrics.

Journal of Aquatic Animal Health↗
USGS1987

Computer simulations of large asteroid impacts into oceanic and continental sites--Preliminary results on atmospheric, cratering and ejecta dynamics

Computer simulations have been completed that describe passage of a 10-km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics caused by impact of the asteroid into both oceanic and continental sites. The asteroid was modeled as a spherical body moving vertically at 20 km/s with a kinetic energy of 2.6 × 10 30 ergs (6.2 × 10 7 Mt ). Detailed material modeling of the asteroid, ocean, crustal units, sedimentary unit, and mantle included effects of strength and fracturing, generic asteroid and rock properties, porosity, saturation, lithostatic stresses, and geothermal contributions, each selected to simulate impact and geologic conditions that were as realistic as possible. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock wave followed by a highly shock compressed and heated air mass. Rapid expansion of this shocked air created a large low-density region that also expanded away from the impact area. Shock temperatures in air reached ∼20, 000K near the surface of the uplifting crater rim and were as high as ∼2000K at more than 30 km range and 10 km altitude. Calculations to 30 s showed that the shock fronts in the air and in most of the expanding shocked air mass preceded the formation of the crater, ejecta, and rim uplift and did not interact with them. As cratering developed, uplifted rim and target material were ejected into the very low density, shock-heated air immediately above the forming crater, and complex interactions could be expected. Calculations of the impact events showed equally dramatic effects on the oceanic and continental targets through an interval of 120 s. Despite geologic differences in the targets, both cratering events developed comparable dynamic flow fields and by ∼29s had formed similar-sized transient craters ∼39km deep and ∼62km across. Transient-rim uplift of ocean and crust reached a maximum altitude of nearly 40 km at ∼30s and began to decay at velocities of 500 m/s to develop large-tsunami conditions. After ∼30s, strong gravitational rebound drove both craters toward broad flat-floored shapes. At 120 s, transient crater diameters were ∼80km (continental) and ∼105km (oceanic) and transient depths were ∼27km; crater floors consisting of melted and fragmented hot rock were rebounding rapidly upward. By 60 s, the continental crater had ejected ∼2 × 10 14 t, about twice the mass ejected from the oceanic crater. By 120 s, ∼70, 000km 3 (continental) and ∼90, 000km 3 (oceanic) target material were excavated (no mantle) and massive ejecta blankets were formed around the craters. We estimate that in excess of ∼70% of the ejecta would finally lie within ∼3 crater diameters of the impact, and the remaining ejecta (∼10 13 t), including the vaporized asteroid, would be ejected into the atmosphere to altitudes as high as the ionosphere. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to contribute substantial material to the atmosphere.

International Journal of Impact Engineering↗
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗