Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “proxies”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Extrapair offspring of the blue-footed booby show no sign of higher fitness in the first 10 years of life

According to the good genes and genetic compatibility hypotheses, females of socially monogamous species obtain genetic benefits for their offspring by performing extrapair copulations with males of higher quality than their social mates or males with whom they are more genetically compatible. If extrapair offspring do receive genetic benefits in the form of advantageous alleles or more compatible allele combinations, they should outperform their within-pair half-siblings' survival and/or reproductive output. Here, in this study, we followed 52 extrapair and 737 within-pair blue-footed booby, Sula nebouxii, offspring during their first 10 years of life (excluding the embryonic period and the first 10 days after hatching) to assess whether they differed in fledging probability, fledgling body condition, recruitment probability, age at first reproduction, number of breeding events or accumulated breeding success. Extrapair and within-pair offspring did not differ in any of these proxies of fitness. Furthermore, we found that extrapair offspring were equally likely to occur in any hatching position. However, differences in fledgling production over the lifetime could not be ruled out, and because only within-pair production of eggs and fledglings was tallied, the possibility remains that extrapair offspring could produce more extrapair offspring later in life than do within-pair offspring. Furthermore, the possibility of context-dependent genetic benefits occurring only under stressful conditions cannot be discounted because our sample of offspring was obtained in a single exceptionally favourable reproductive season.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Multiphysics modeling of a critical dual-purpose canister in a saturated geological repository

Dual purpose dry cask storage canisters for spent nuclear fuel are designed for storage and transportation, but are not licensed for permanent disposal in a geological repository. If dual purpose canisters were to be used to dispose of spent nuclear fuel in a geological repository, they would be expected to eventually breach and be flooded with groundwater, and it is shown that some fraction of these canisters will achieve criticality. In order to evaluate the consequences of canisters going critical in a repository, an initial capability has been developed for estimating the quasi-static power level of a critical canister using loosely coupled multiphysics simulations. The low power level in a critical canister enables coupling through precomputed physics proxies. This calculated power level is then used to compute the change in the critical canister’s isotopic inventory as a function of time. Three as-loaded canisters are evaluated and two were found to have power levels below 4 kW, with a modest effect on the radiological inventory over time. This effort also shows that although some DPCs will have extremely peaked power shapes, the relatively low power and long-time scales result in relatively homogenized thermohydraulic properties in the water within the DPC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Co-optimizing water-alternating-carbon dioxide injection projects using a machine learning assisted computational framework

In this report, a robust machine-learning-based computational framework that couples multi-layer neural network (MLNN) proxies and a multi-objective particle swarm optimizer (MOPSO) to design water-alternating-carbon dioxide injection (CO 2 -WAG) projects is presented. The proposed optimization protocol considers various objectives, including oil recovery and CO 2 storage volume. Expert MLNN systems are trained and employed as surrogate models of the high-fidelity compositional simulator in the optimization workflow. When multiple objective functions are considered, two approaches are employed to treat the objectives: the weighted sum method and the Pareto-front-based scheme. A field-scale implementation focusing on tertiary recovery in the Morrow B formation at Farnsworth Unit (FWU) is presented. The developed Pareto-optimal solutions indicate the maximal available oil production can be 1.64 × 10 7 barrels and maximal carbon storage can achieve 2.35 × 10 7 tons. Trade-offs factor is defined to divide the constructed Pareto front into 4 sections with the trade-off factors’ value ranges from 0.35 to 49.9. This work also compares the optimum solution found by the aggregative objective function and the solution repository covered by the Pareto front that considers the physical and operational constraints and reduces uncertainties involved by the multi-objective optimization process. Our comparison indicates multiple solutions exist to satisfy the objective criteria of the WAG design, and these results cannot be found using the traditional weighted sum method. The Pareto front solution can provide more options for project designers, but decisions regarding necessary trade-offs must be made using the solution repository to balance the project economics and CO 2 storage amount.

10 SYNTHETIC FUELS↗
DOE OSTI2021

The role of current and emerging technologies in meeting Japan’s mid- to long-term carbon reduction goals

Using Japan as a proxy for a developed nation, we investigated the role of existing and nascent technologies in curbing carbon emissions. We simulated possible pathways to meeting 2030 and 2050 emission targets within the Japanese electricity supply sector using a single-region model in The Integrated MARKAL-EFOM System (TIMES). Critically, our simulations incorporate novel technologies like hydrogen electrolysers, carbon capture, photochemical water splitting, and emerging photovoltaic cells, assess long-term impacts up to the year 2100, and include life-cycle emissions and learning curves for parameters such as investment cost, efficiency, and emission coefficients. Results indicate that a hybrid approach, using nuclear power and hydrogen from renewable energy-based electrolysis, is cost-effective and provides long-term emission reduction along with energy security. Finally, nuclear, wind, solar, and hydrogen from renewables emerge as key emission reduction technologies, while natural gas with carbon capture plays a minor role in achieving emission reduction targets.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Dynamic development of geochemical reaction fronts during hydraulic stimulation of shale

Injection of acidic hydraulic fracture fluid (HFF) into shale formations for unconventional oil/gas production results in chemical reactions in the shale matrix that can alter fluid transport. Here, we report the results of set of experiments designed to evaluate the impact of calcite dissolution as a function of carbonate mineral content on matrix chemical reactivity and pore-space modification concomitant with imbibition. In this study, we tracked acidic HFF transport in four samples of Wolfcamp shale with calcite contents varying from 4% to 59% by monitoring the rate and spatial extent of bromide tracer transport using synchrotron-based X-ray fluorescence microprobe (XFM) imaging. Concurrently, we also carried out XFM imaging of the spatial distribution of Ca in the Wolfcamp shale cores (as a proxy of calcite distribution). Our approach thus yields a direct record of time-resolved selective ion transport resulting from the penetration of acidic HFF and the associated mineral transformations in the shale cores. We show that the variability in calcite content of Wolfcamp shale samples can directly affect the rate and spatial extent of imbibition. Although reaction of the acidic HFF with carbonates in shales enhances calcite dissolution and increases porosity, the spatial extent of calcite dissolution in the shale matrix is limited due to a rapid neutralization of pH. The relative abundance and spatial distribution of calcite control the chemical saturation state of the HFF progressing into the matrix. As a result, calcite has a major impact on the spatial extent and rate of matrix alteration and thus on HFF transport during subsurface reservoir stimulation. Consequently, increased calcite content in the shale matrix inhibits the spatial extent of the pore-volume increase and, by extension, the spatial extent and rate of imbibition. Our results thus show that the overall rates of calcite dissolution approach the rates of acidic HFF transport (i.e., Damköhler number ~1), which could contribute to the efficiency of subsurface reservoir stimulation. A better understanding of HFF-calcite reaction rates is crucial for improving the prediction and optimization of fluid transport across HFF-shale interfaces during hydraulic fracturing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

How much information on precipitation is contained in satellite infrared imagery?

Satellite infrared (IR) brightness temperature providing information on the characteristics of cloud tops is one of the primary input data used for estimating precipitation from satellites. Thus, the natural questions to ask are: Does IR imagery contain more information about precipitation in some regions more than others? What are the factors that lead to such differences? Should we expect the amount of information provided by IR imagery to change with respect to temporal or spatial aggregation? Comprehensive analysis of the accuracy of satellite-based precipitation estimates obtained from any given algorithm may provide qualitative answers to the above questions. However, such an approach only offers conclusions within the limited context of the algorithm used to obtain precipitation estimates; that is, patterns caused by the inherent information content of IR imagery are inextricable from those induced by assumptions embedded in the precipitation estimation algorithm. Here, we address the abovementioned questions from a different perspective using an information-theoretic measure, free of structural assumptions and general to a wide range of relationships, to characterize the average and seasonal dependence between IR imagery provided by satellite sensors in the spectral range (10.7–11.5 μm) and precipitation rates in the contiguous United States (CONUS) across distinct temporal and spatial scales. We analyze a total of more than 1.3 billion pairs of IR and precipitation observations over CONUS, and we observe interesting patterns. First, we show that there is a strong inverse relationship (ρ = - 0.73) between the information content of IR data and the number of no-rain observations; however, its robustness varies regionally and seasonally with less significant correlation during the warm season. Second, we demonstrate that the intuitive relationship of increasing dependence between IR and precipitation as a result of temporal or spatial aggregation exhibits a diminishing returns behavior. For instance, temporal aggregation from 1 to 3 h increases the dependence approximately 7 times as much as temporal aggregation from 12 to 24 h. Furthermore, we quantitatively examine and re-confirm several statements previously reported in the literature such as the strong association between IR and precipitation in convective storms and the low association in orographic rainfall. Finally, we conclude by illustrating the potential of the analysis in diagnosis of operational algorithms for estimating precipitation from IR brightness temperature. The findings of this study pinpoint spatial domains and time scales at which IR is not an adequate proxy for estimating precipitation. Thus, they can potentially guide the development of operational algorithms that utilize satellite Infrared imagery for estimating, downscaling and data fusion of precipitation. Finally, the analysis presented here opens up the possibility of developing new methodologies for diagnosis of satellite-based precipitation estimation algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Disturbance and diversity: Lichen species richness decreases with increasing anthropogenic disturbance

Anthropogenic disturbance is rapidly increasing through habitat degradation, development, and deforestation. Gaps remain in understanding the effects of this disturbance on diverse and ecologically important organisms such as lichens. In North America, studies have focused on epiphytic macrolichens and catastrophic disturbance, largely ignoring microlichens and less severe disturbances. The Intermediate Disturbance Hypothesis predicts these moderate disturbances will lead to higher species richness. Here we examine the effects of anthropogenic disturbance in the form of land management on overall lichen richness, including microlichens, and on the species richness of eight lichen functional groups. The study draws on a comprehensive data set of 872 species, in 208 one-hectare plots throughout the Southern Appalachian Mountains, a global biodiversity hotspot in eastern North America. A habitat quality index based on an established forestry metric was used as a proxy for anthropogenic disturbance and was quantified using a 10-part score including categories such as percent native tree canopy cover and degree of fragmentation. Linear models were used to compare habitat quality scores to overall species richness and to species richness of functional groups. Rather than following the Intermediate Disturbance Hypothesis, all groups uniformly follow a negative linear relationship: as disturbance increases, species richness decreases. This pattern held even for widely variable functional groups such as morphotype. Effective conservation of lichen richness should prioritize the maintenance of existing older, less-disturbed stands within large, contiguously forested areas.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Optimization-based, property-preserving finite element methods for scalar advection equations and their connection to Algebraic Flux Correction

In this paper, we continue our efforts to exploit optimization and control ideas as a common foundation for the development of property-preserving numerical methods. Here we focus on a class of scalar advection equations whose solutions have fixed mass in a given Eulerian region and constant bounds in any Lagrangian volume. Our approach separates discretization of the equations from the preservation of their solution properties by treating the latter as optimization constraints. This relieves the discretization process from having to comply with additional restrictions and makes stability and accuracy the sole considerations in its design. A property-preserving solution is then sought as a state that minimizes the distance to an optimally accurate but not property-preserving target solution computed by the scheme, subject to constraints enforcing discrete proxies of the desired properties. Furthermore, we consider two such formulations in which the optimization variables are given by the nodal solution values and suitably defined nodal fluxes, respectively. A key result of the paper reveals that a standard Algebraic Flux Correction (AFC) scheme is a modified version of the second formulation obtained by shrinking its feasible set to a hypercube. In conclusion, we present numerical studies illustrating the optimization-based formulations and comparing them with AFC

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Implementing a neural network interatomic model with performance portability for emerging exascale architectures

The two main thrusts of computational science are increasingly accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved thread scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution enables simulations up to 20 million atoms on a single CPU node and 4 million atoms with improved performance on a single GPU. Furthermore, we also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~50% and ~5% improvements in performance on a GPU by choosing the right level of parallelism and data layout respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Monitoring and modeling hydrologic conditions in Ukraine for hydropower generation

Study region: The Dnieper and Dniester Rivers of Ukraine. Study focus: The ongoing conflict in Ukraine has caused disruptions to electricity generation, of which hydroelectric sources contribute approximately 9 % to the country’s needs. With the takeover of the Zaporizhzhia nuclear power plant by enemy forces, the loss of the Kakhovka hydroelectric dam, and the future impacts of the conflict on electricity generation unclear, it may be valuable for the Ukrainian government to better understand how it could leverage hydroelectric power sources in the near future. Unfortunately, measurements of river discharge throughout Ukraine ceased data collection in the late 1980’s to early 1990’s. To address this data gap, we developed a protocol that combined satellite-based time-series measurements of river width at seven locations throughout Ukraine from 2013 to 2023 with reanalysis data, climate-model predictions, and hydrologic models to both provide a means of monitoring a proxy for near-real-time discharge and also predict near-term (i.e., 2023–2030) hydrologic patterns for the region. New hydrological insights for the region: We ran new algorithms on 144 WorldView-2 and WorldView-3 satellite images to map rivers and extract width, one of which was validated against river gauge data located along the same river but in a neighboring country. Hydrologic models using two climate scenarios found minimal change in annual discharge at all sites, but magnitude and timing of peak discharge showed a moderate trend. The results suggest that hydropower is underutilized in Ukraine.

13 HYDRO ENERGY↗
DOE OSTI2020

Spatio-temporal impacts of a utility’s efficiency portfolio on the distribution grid

Energy Efficiency has historically focused on delivering savings to offset growth in energy supply. Today’s growing emphasis on decarbonization of the energy supply is driving renewables adoption and increased interest in electrification. As a result, energy efficiency is being assessed not just in its ability to offset load growth, but also for its ability to alleviate location-specific constraints on transmission and distribution infrastructure. This work demonstrates that advanced measurement and verification modeling techniques can be used to estimate the spatio-temporal grid impact of a portfolio of energy efficiency programs. It extends measurement-based methods to an entire Demand Side Management portfolio and uses a single model to predict annual as well as seasonal building energy use with near-zero bias. In addition, new metrics are introduced to assess grid level impacts of energy efficiency. The results show that the efficiency program portfolio delivers savings of over 12% at the territory-wide proxy level, with substation and feeder level savings ranging from 0.4% to 26%, and ~5%-42% respectively. Overall, these savings impacted 1.0%–1.4% of the energy used at these locations in the grid. This work provides a methodology with potential to connect efficiency with distribution planning, carrying implications for non-wires alternatives and targeted delivery of efficiency programs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Impact of operating reserve rules on electricity prices with high penetrations of renewable energy

In competitive wholesale electricity markets, significant effort is devoted to designing markets that set efficient prices for maintaining supply-demand balance. One factor that can impact prices is administratively-set scarcity pricing, which sets prices to a preset level when the market is not able to meet operating reserve or energy requirements. When energy and operating reserves are co-optimized, assumptions surrounding operating reserve requirements and scarcity pricing can impact system-wide price outcomes for both operating reserves and energy. This study uses production cost modeling of an ERCOT-like system to evaluate the impact of operating reserve eligibility, scarcity pricing, and quantity rules on electricity prices, and therefore also on generator revenues. Results reveal economic and operational benefits with allowing open participation in reserve markets, as this enables greater access to the full set of capable resources at lowest cost. Furthermore, both energy and reserve prices are strongly impacted by reserve scarcity pricing events, which reveals that reserve scarcity pricing assumptions can impact price outcomes even for units not providing reserves. This study highlights the importance of operating reserve scarcity pricing rules because of the strong coupling between energy and reserve prices and because these rules serve as proxies for true price responsive demand.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

No detectable redox exchange between sulfur and iron during rapid cooling of basalts

The valence states of iron and sulfur in mantle-derived melts influence and respond to the composition of the mantle as well as the subsequent evolution of those melts in Earth's crust. Characterizing the valence states of iron and sulfur in silicate melts that have quenched to form glasses is therefore critical to understanding the composition of the mantle and the magmatic processes that take place before and during volcanic eruptions. Glasses quenched from silicate melts offer the opportunity to precisely measure sulfur and iron valence states, but whether the transition from melt to glass affects the valence state of sulfur and iron is uncertain. Here, we use the glassy margins of two pillow basalts, one from a mid-ocean ridge and one from a back-arc basin, as a natural experiment to test the effect of quench rate on sulfur and iron valence states. We report micro-X-ray Absorption Near Edge Structure (XANES) measurements of S 6+ /ΣS (S 6+ /[S 6+ +S 2- ]) and Fe 3+ /ΣFe (Fe 3+ /[Fe 2+ +Fe 3+ ]) along transects from the rapidly quenched rims of the pillows down into the slowly cooled crystal-rich interiors. The range of average quench rates estimated for our selected analysis areas range from 1.2 x 10 5 °C/sec at the pillow margin to 1.1 °C/sec in the interior and span quench rates experienced by natural and experimental glasses. On average, S6+/ΣS = 0.09 ± 0.01 and Fe 3+ /ΣFe = 0.147 ± 0.002 in the mid-ocean ridge sample and S 6+ /ΣS = 0.21 ± 0.01and Fe 3+ /ΣFe = 0.153 ± 0.004 in the back-arc sample, where the quoted uncertainty is the 1σ standard deviation of n= 9-16 analyses. We find that the S 6+ /ΣS and Fe 3+ /ΣFe values along transects in each sample show no systematic changes with depth, and by proxy, no systematic changes with quench rate. In each pillow, all measured S 6+ /ΣS and Fe 3+ /ΣFe ratios are within the absolute uncertainty of each technique. We propose that these observations are consistent with either sluggish reaction kinetics between sulfur and iron, or an insignificant temperature effect on sulfur-iron redox exchange equilibria. These observations demonstrate that sulfur and iron valence states in glasses with varying quench rates can be directly compared, and that changes in sulfur and iron valence during quench are below the limits of detection with XANES.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

On the feasibility of using physics-informed machine learning for underground reservoir pressure management

In this work, we evaluate the feasibility of using physics-informed machine learning (PIML) for underground energy-related pressure management. To this end, we develop a PIML framework to manage underground reservoir pressures by training neural networks to determine fluid extraction rates for dedicated extraction wells during fluid injection operations given a range of reservoir conditions (e.g., transmissivity and storativity). We implement an automatically-differentiable analytical physics model of fluid flow in porous media within the PIML framework as a proxy for more complicated models. This allows us to execute a sufficient number of training scenarios to fully evaluate the feasibility of using PIML to support pressure management activities. We quantify the number of physics-model parameters required for automatic differentiation to become more efficient than finite-difference gradient calculations. We use a simple scenario with a single injector, extractor, and critical location for our feasibility analysis. We evaluate the effect of the size of the training dataset (i.e., the number of reservoir condition samples) on the accuracy and efficiency of the PIML framework. For an equivalent number of model evaluations, the larger training dataset took less time to train and produced a neural network that was able to more accurately manage reservoir pressures. We also evaluate the effect of the training dataset batch size (i.e., number of reservoir condition samples used to update the neural network coefficients during training; i.e., how the training dataset is partitioned). While training ran faster with larger batch sizes, they produced neural networks that managed pressures less accurately. We demonstrate the approach on a more complex scenario involving 10 injectors, 10 extractors, and 4 critical locations (a relatively high well density of 20 wells/km2). We provide the number of forward and adjoint model evaluations required in each case as an indication of the feasibility of using PIML for pressure management when more complicated physics models with longer execution times are used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Optical diagnostic design for measuring the radiation front with a mid-leg pumped divertor on DIII-D

The design of an optical diagnostic system to localize the radiation front in a mid-leg pumped divertor configuration on DIII-D is presented. Divertor detachment is a key mechanism for handling power exhaust in tokamaks, and mid-leg pumping offers a promising approach to radiatively dissipate power while maintaining acceptable core performance. To predict the location of the radiation front and design a spectroscopic diagnostic to measure its position, a database of SOLPS-ITER simulations across a range of input powers and gas puffing rates representative of DIII-D operation was employed. These simulations provide self-consistent plasma backgrounds for a Cherab-Raysect synthetic diagnostic framework, which incorporates detailed tokamak geometry and physically accurate ray-tracing. Synthetic simulations of optical sightlines and viewing cone geometry were used to calculate line-integrated emission from the plasma, particularly of the C III 465 nm line, to serve as a proxy for the T e ≈ 7 − 10 eV temperature region associated with the onset and evolution of divertor detachment. Neutral deuterium emission is also evaluated for comparison. The synthetic diagnostic is used to assess the performance of existing DIII-D optical systems, including filterscopes and the Multichord Divertor Spectrometer, and to optimize line-of-sight placement within mechanical and installation constraints. The results provide quantitative guidance for diagnostic implementation on DIII-D and demonstrate the advantage of integrated synthetic diagnostics for divertor design studies and future advanced divertor concepts.

Cherab↗
DOE OSTI2020

Deposition of sulfate aerosols with positive Δ 33 S in the Neoarchean

Anomalous sulfur isotope compositions present in Archean rocks have been intensely scrutinized over the last 20 years because they record key aspects of Earth's atmospheric composition prior to the appearance of free molecular oxygen ca. 2.3 billion years ago. These isotopic compositions can be described as mass anomalous fractionations (MAF) and are produced in the atmosphere as UV light interacts with SO 2 molecules. Most interpretations suggest that atmospheric processes generate a reduced S-phase with a positive ( 33 S-enriched) MAF signature, as measured in pyrites, and an oxidized S-phase with a negative anomaly, as measured in bedded barite deposits. However, recent data for carbonate-associated sulfate (CAS) — a direct proxy for the isotopic composition of sulfur from seawater sulfate — in Neoarchean rocks showed no such negative values, but rather the opposite. To understand if the positive MAF anomalies we measured in Neoarchean CAS reflect secondary processes (diagenetic, metamorphic, handling) instead of original signals of Archean seawater sulfate, we collected additional sample suites with various degrees of preservation and metamorphic alteration across the Campbellrand-Malmani platform in South Africa. Results illustrate that within this comprehensive suite, less-altered samples all contain positive MAF values while secondary processes tend to either remove CAS from the sample and/or decrease the 33 S-enrichment. This positive MAF signal in sulfate is therefore reasonably interpreted as a primary depositional origin, and implies that the assumption that sulfate always carries a negative MAF anomaly throughout the Archean rock record needs to be reconsidered. Our CAS observations suggest that future experiments and calculations should also consider atmospheric and/or sulfur cycling processes that can produce oxidized sulfur with a positive MAF signature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Past fires and post-fire impacts reconstructed from a southwest Australian stalagmite

Speleothem records of past environmental change provide an important opportunity to explore fire frequency and intensity in the past, and the antecedent climatic conditions leading to fire events. Here, fire sensitive geochemical signals in a stalagmite from Yonderup Cave, a shallow cave in Western Australia, are compared to well-documented wildfire events that occurred in recent decades. The results are extended to identify wildfires during the growth interval of the speleothem (1760 CE–2005 CE). Principal component analysis of the stalagmite time-series revealed distinct peaks in a combination of phosphorus and metal (aluminium, zinc, copper and lead) concentrations in response to known fire events, which are interpreted to have come from ash. Varying responses in the geochemical signal in the pre- and post-European colonisation period are likely linked to changes in land management, fire frequency and fire intensities. Of note is what we infer to be a particularly intense fire event concluded to have occurred in 1897 ± 5 CE, which climate sensitive proxies (δ18O and bedrock-derived elements) indicate was preceded by a multi-decadal dry period that began in the late 1860s. The intensity of the deduced fire event is supported by a peak in P that is 6.3 times higher than elsewhere in the speleothem, accompanied by a peak in Zn. It was also inferred that changes to the surface-cave hydrology occurred as a result of this event, indicated by increased input of colloidal organic material onto the stalagmite that was interpreted to be caused by shallow karst bedrock fracturing from the heat-induced deformation that is observed after intense fires. These findings show the potential for speleothems to provide dated records of fire intensity and recurrence intervals. Further development could lead to a better understanding of the climate-fire relationship and the effects of land-management practices on wildfire frequency and intensity.

58 GEOSCIENCES↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗