Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Reduced ecosystem resilience quantifies fine-scale heterogeneity in tropical forest mortality responses to drought

Sensitivity of forest mortality to drought in carbon-dense tropical forests remains fraught with uncertainty, while extreme droughts are predicted to be more frequent and intense. Here, in this study, the potential of temporal autocorrelation of high-frequency variability in Landsat Enhanced Vegetation Index (EVI), an indicator of ecosystem resilience, to predict spatial and temporal variations of forest biomass mortality is evaluated against in situ census observations for 64 site-year combinations in Costa Rican tropical dry forests during the 2015 ENSO drought. Temporal autocorrelation, within the optimal moving window of 24 months, demonstrated robust predictive power for in situ mortality (leave-one- out cross-validation R 2 = 0.54), which allows for estimates of annual biomass mortality patterns at 30 m resolution. Subsequent spatial analysis showed substantial fine-scale heterogeneity of forest mortality patterns, largely driven by drought intensity and ecosystem properties related to plant water use such as forest deciduousness and topography. Highly deciduous forest patches demonstrated much lower mortality sensitivity to drought stress than less deciduous forest patches after elevation was controlled. Our results highlight the potential of high-resolution remote sensing to “fingerprint” forest mortality and the significant role of ecosystem heterogeneity in forest biomass resistance to drought.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Iodine Air Sampling: Installation and Preliminary Results

The F-Area wetlands are a known sink for radioiodine generated at the Savannah River Site (SRS). Wastewater containing iodine-129 ( 129 I) was discharged into the F-Area Seepage Basins for disposal, but the acidity of the wastewater mobilized some of the contaminants from the basin soil through the vadose zone and into the groundwater; ultimately creating a groundwater plume that extends towards Fourmile Branch and its associated wetland areas. Dissolved radioiodine specifically has migrated to the nearby F-Area wetlands, where it persists as a contaminant of concern. Accumulation and resurfacing of 129 I in the wetlands and nearby areas creates a potential secondary source zone where changes in biogeochemical conditions could result in the release of this sequestered 129 I to the air. Monitoring data has indicated that 129 I has a seasonal behavior at several surface water stations, with high concentrations during summertime and low concentrations during wintertime. Still, it is not known which factors control the seasonal behavior of 129 I. Further, it has been reported that 129 I contributed approximately 23% to the 2019 SRS offsite air pathway dose to the representative person for that year (SRNS, 2020), but air measurements of 129 I have never been taken in the F-Area wetlands. Because air emissions of 129 I from the wetlands have never been investigated, it is not known if this could be a potential release pathway that needs to be included in the site conceptual model. To address these unknowns, the Savannah River National Laboratory (SRNL) has an applied research task to investigate whether 129 I in the wetlands is emitted to the air and, if so, whether air concentrations exhibit seasonal variation. This report presents information about the system recently installed in the F-Area wetlands for measuring 129 I, as well as analysis results for samples collected June through August 2023. Quantification of this constituent will provide a better understanding of the contribution of airborne 129 I from the F-Area wetlands to the offsite dose and better inform ongoing and future remedial actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reducing the computational burden of a microgrid energy management system

As renewable technology advances and decreases in cost, microgrids are becoming an appealing means of distributed generation both for isolated communities and integrated with existing electrical grid systems. Due to their small size, however, microgrids may have financial limitations which preclude them from using commercial software to optimize control of their assets. Open-source optimization solvers are a viable alternative, but increase computation time. This work expands on a rolling horizon optimization framework for economic dispatch within an existing residential microgrid located in Hoover, Alabama. The microgrid has an open-source solver requirement and a need for quick solution time on a rolling horizon as opposed to a day-ahead commitment. We present a method of reducing integer variables by relaxation which completes two goals: reduction in computation time for real-time operations, and reduction in daily operational cost for the microgrid. Seasonal data for load and photovoltaic (PV) power was also collected from the microgrid to facilitate simulation testing. Computation time was successfully reduced using multiple variations of the relaxation method, while obtaining solution quality with operational cost similar to or better than the original model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Photonic ring resonator notch filters for astronomical OH line suppression

Photonic ring resonator arrays used as notch filters are a promising novel solution to improve the signal-to-noise ratio of ground-based astronomical observations by suppressing OH emission lines in the near-infrared (NIR) wavelength range (0.9-2.5 μm). We aim to fabricate a series of ring resonators connected by a waveguide, each with its resonance wavelength and full-width-half-maximum (FWHM) matched with one of the OH emission lines. The fundamental structure of ring resonator-based notch filters is a ring resonator evanescently coupled to a waveguide. Due to the optical loss of the ring resonator, the transmission of light through the waveguide will be suppressed at the resonance wavelength of the ring resonator. The degree of suppression can be controlled by varying the coupling strength between the waveguide and the ring resonator, while the resonance wavelength can be controlled by the geometry and size of the ring. However, since a single ring resonator has multiple modes each with a different wavelength separated by its free spectral range (FSR), the additional suppression at other wavelengths may undermine the signal. Therefore, the ideal FSR off the ring resonators needs to be larger than the measurement window. For this application, we aim the FSR to be 1 larger than 30 nm . Our initial measurement results revealed two major challenges for this application. First, the FSR of the ring resonators is limited by the bending loss of ring resonators. In a standard 220 nm silicon-on-insulator material platform, the FSR of ring resonators cannot exceed 30 nm; Second, due to the wafer thickness variation and fabrication inaccuracy, the achievable accuracy in the resonance wavelength of a fabricated single ring resonator is only ~1.5 nm. In this work, we report our novel solutions to the above two challenges. For the first one, we have developed a novel structure based on the Vernier effect to greatly raise the FSR of ring resonators. Serially coupled Vernier ring structures is a typical approach to increase the FSR without encountering bending loss, but it requires restrictive matching of the coupling and optical loss of the rings and any imperfect matching can result in even noisier signals. On the contrary, our structure, as shown in shown in Fig. 1, does not require the restrictive matching. The signals of two output ports in our structure marked in Fig. 1 will be added together to form the final output signal, and Fig. 2 shows the finite-difference time-domain (FDTD) simulation results of the final output signal. As can be seen, the FSR of the entire structure is well-beyond 60 nm with little noise at other wavelengths. For the second one, we have developed an easy-to-implement method to permanently tune the resonance wavelengths of ring resonators within 0.1 nm accuracy. The tuning mechanism is based on the property of hydrogen silsesquioxane (HSQ) resist. As a negative tone electron beam resist, HSQ upon exposure will cross-link and eventually form SiO2. The refractive 2 index of HSQ resist and formed SiO2 are slightly different , enabling the tuning. By having a thin layer of HSQ resist on top of ring resonator structures followed by SiO2 cladding, the resonance wavelengths of ring resonators can be post-tuned Advances in Optical Astronomical Instrumentation 2019, edited by Simon Ellis, Céline d'Orgeville, Downloaded From: https://www.spiedigitallibrary.org/conference-proceedings-of-spie on 23 Sep 2020 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use through electron beam exposure. We have found that, within 5 nm tuning range, the resonance wavelength shift and dosage applied obeys a linear relationship. So a simple calibration beforehand can result in resonance wavelength tuning within 0.1 nm accuracy. The transmission spectra of a 3-ring device before and after tuning is shown in Fig. 3.

Astronomical optics↗
DOE OSTI2021

Advanced Electron Microscopy of Nanophased Synthetic Polymers and Soft Complexes for Energy and Medicine Applications

After decades of developments, electron microscopy has become a powerful and irreplaceable tool in understanding the ionic, electrical, mechanical, chemical, and other functional performances of next-generation polymers and soft complexes. The recent progress in electron microscopy of nanostructured polymers and soft assemblies is important for applications in many different fields, including, but not limited to, mesoporous and nanoporous materials, absorbents, membranes, solid electrolytes, battery electrodes, ion- and electron-transporting materials, organic semiconductors, soft robotics, optoelectronic devices, biomass, soft magnetic materials, and pharmaceutical drug design. For synthetic polymers and soft complexes, there are four main characteristics that differentiate them from their inorganic or biomacromolecular counterparts in electron microscopy studies: (1) lower contrast, (2) abundance of light elements, (3) polydispersity or nanomorphological variations, and (4) large changes induced by electron beams. Since 2011, the Center for Nanophase Materials Sciences (CNMS) at Oak Ridge National Laboratory has been working with numerous facility users on nanostructured polymer composites, block copolymers, polymer brushes, conjugated molecules, organic–inorganic hybrid nanomaterials, organic–inorganic interfaces, organic crystals, and other soft complexes. This review crystalizes some of the essential challenges, successes, failures, and techniques during the process in the past ten years. It also presents some outlooks and future expectations on the basis of these works at the intersection of electron microscopy, soft matter, and artificial intelligence. Machine learning is expected to automate and facilitate image processing and information extraction of polymer and soft hybrid nanostructures in aspects such as dose-controlled imaging and structure analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chemical Control of Spin‐Orbit Coupling and Charge Transfer in Vacancy‐Ordered Ruthenium(IV) Halide Perovskites

Abstract Vacancy‐ordered double perovskites are attracting significant attention due to their chemical diversity and interesting optoelectronic properties. With a view to understanding both the optical and magnetic properties of these compounds, two series of Ru IV halides are presented; A 2 RuCl 6 and A 2 RuBr 6 , where A is K, NH 4 , Rb or Cs. We show that the optical properties and spin‐orbit coupling (SOC) behavior can be tuned through changing the A cation and the halide. Within a series, the energy of the ligand‐to‐metal charge transfer increases as the unit cell expands with the larger A cation, and the band gaps are higher for the respective chlorides than for the bromides. The magnetic moments of the systems are temperature dependent due to a non‐magnetic ground state with J eff =0 caused by SOC. Ru‐ X covalency, and consequently, the delocalization of metal d ‐electrons, result in systematic trends of the SOC constants due to variations in the A cation and the halide anion.

Vishnoi, Pratap↗
DOE OSTI2021

Quantitative assessment of the influence of external magnetic field on clustering of nano-Fe3O4 particles in cementitious paste

In view of active rheology control of cementitious materials, nano-Fe{sub 3}O{sub 4} can be added as responsive particles. Following the concept of magnetorheological fluids, it is assumed that magnetic nanoparticles will form chains or clusters in cementitious paste following magnetic field lines. A quantitative experimental validation of this assumption is presented herein. The clustering of nano-Fe{sub 3}O{sub 4} particles under magnetic fields is studied by mapping iron (Fe) element distribution in cementitious paste using energy dispersive X-ray spectroscopy. By means of image analysis, the Fe-element patterns are quantified by the deviation of Fe-elements in a unit area from the mean value expected in case of a uniform distribution, as expressed by coefficient of variation (COV). The magneto-rheological responses of cementitious pastes are evaluated using small amplitude oscillatory shear technique. Results show that the magneto-rheological effect exhibits a linear relationship with the relative change of COV, providing a quantitative validation of magnetic clustering in cementitious paste.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal Management for FPGA Nodes in HPC Systems

The integration of FPGAs into large-scale computing systems is gaining attention. In these systems, real-time data handling for networking, tasks for scientific computing, and machine learning can be executed with customized datapaths on reconfigurable fabric within heterogeneous compute nodes. At the same time, thermal management, particularly battling the cooling cost and guaranteeing the reliability, is a continuing concern. The introduction of new heterogeneous components into HPC nodes only adds further complexities to thermal modeling and management. The thermal behavior of multi-FPGA systems deployed within large compute clusters is less explored. Here, we first show that the thermal behaviors of different FPGAs of the same generation can vary due to their physical locations in a rack and process variation, even though they are running the same tasks. We present a machine learning–based model to capture the thermal behavior of each individual FPGA in the cluster. We then propose two thermal management strategies guided by our thermal model. First, we mitigate thermal variation and hotspots across the cluster by proactive thermal-aware task placement. Under the tested system and benchmarks, we achieve up to 26.4° C and on average 13.3° C system temperature reduction with no performance penalty. Second, we utilize this thermal model to guide HLS parameter tuning at the task design stage to achieve improved thermal response after deployment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Mechanistic Modeling of Microtopographic Impacts on CO 2 and CH 4 Fluxes in an Alaskan Tundra Ecosystem Using the CLM-Microbe Model

Spatial heterogeneities in soil hydrology have been confirmed as a key control on CO 2 and CH 4 fluxes in the Arctic tundra ecosystem. In this study, we applied a mechanistic ecosystem model, CLM-Microbe, to examine the microtopographic impacts on CO 2 and CH 4 fluxes across seven landscape types in Utqia&#289vik, Alaska: trough, low-centered polygon (LCP) center, LCP transition, LCP rim, high-centered polygon (HCP) center, HCP transition, and HCP rim. We first validated the CLM-Microbe model against static-chamber measured CO 2 and CH 4 fluxes in 2013 for three landscape types: trough, LCP center, and LCP rim. Model application showed that low-elevation and thus wetter landscape types (i.e., trough, transitions, and LCP center) had larger CH 4 emissions rates with greater seasonal variations than high-elevation and drier landscape types (rims and HCP center). Sensitivity analysis indicated that substrate availability for methanogenesis (acetate, CO 2 + H 2 ) is the most important factor determining CH 4 emission, and vegetation physiological properties largely affect the net ecosystem carbon exchange and ecosystem respiration in Arctic tundra ecosystems. Modeled CH 4 emissions for different microtopographic features were upscaled to the eddy covariance (EC) domain with an area-weighted approach before validation against EC-measured CH 4 fluxes. The model underestimated the EC-measured CH 4 flux by 20% and 25% at daily and hourly time steps, suggesting the importance of the time step in reporting CH 4 flux. The strong microtopographic impacts on CO 2 and CH 4 fluxes call for a model-data integration framework for better understanding and predicting carbon flux in the highly heterogeneous Arctic landscape.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Electrical Transport Interplay with Charge Density Waves, Magnetization, and Disorder Tuned by 2D van der Waals Interface Modification via Elemental Intercalation and Substitution in ZrTe 3 , 2H-TaS 2 , and Cr 2 Si 2 Te 6 Crystals

Electrical transport in 2D materials exhibits unique behaviors due to reduced dimensionality, broken symmetries, and quantum confinement. It serves as both a sensitive probe for the emergence of coherent electronic phases and a tool to actively manipulate many-body correlated states. Exploring their interplay and interdependence is crucial but remains underexplored. This review integratively cross-examines the atomic and electronic structures and transport properties of van der Waals-layered crystals ZrTe 3 , 2H-TaS 2 , and Cr 2 Si 2 Te 6 , providing a comprehensive understanding and uncovering new discoveries and insights. A common observation from these crystals is that modifying the atomic and electronic interface structures of 2D van der Waals interfaces using heteroatoms significantly influences the emergence and stability of coherent phases, as well as phase-sensitive transport responses. In ZrTe 3 , substitution and intercalation with Se, Hf, Cu, or Ni at the 2D vdW interface alter phonon–electron coupling, valence states, and the quasi-1D interface Fermi band, affecting the onset of CDW and SC, manifested as resistance upturns and zero-resistance states. We conclude here that these phenomena originate from dopant-induced variations in the lattice spacing of the quasi-1D Te chains of the 2D vdW interface, and propose an unconventional superconducting mechanism driven by valence fluctuations at the van Hove singularity, arising from quasi-1D lattice vibrations. Short-range in-plane electronic heterostructures at the vdW interface of Cr 2 Si 2 Te 6 result in a narrowed band gap. The sharp increase in in-plane resistance is found to be linked to the emergence and development of out-of-plane ferromagnetism. The insertion of 2D magnetic layers such as Mn, Fe, and Co into the vdW gap of 2H-TaS 2 induces anisotropic magnetism and associated transport responses to magnetic transitions. Overall, 2D vdW interface modification offers control over collective electronic behavior, transport properties, and their interplays, advancing fundamental science and nanoelectronic devices.

2D van der Waals interfaces↗
DOE OSTI2025

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY↗
DOE OSTI2023

Advanced monitoring of soil-vegetation co-dynamics reveals the successive controls of snowmelt on soil moisture and on plant seasonal dynamics in a mountainous watershed

Evaluating the interactions between above- and below-ground processes is important to understand and quantify how ecosystems respond differently to atmospheric forcings and/or perturbations and how this depends on their intrinsic characteristics and heterogeneity. Improving such understanding is particularly needed in snow-impacted mountainous systems where the complexity in water and carbon storage and release arises from strong heterogeneity in meteorological forcing and terrain, vegetation and soil characteristics. This study investigates spatial and temporal interactions between terrain, soil moisture, and plant seasonal dynamics at the intra- and inter-annual scale along a 160 m long mountainous, non-forested hillslope-to-floodplain system in the upper East River Watershed in the upper Colorado River Basin. To this end, repeated UAV-based multi-spectral aerial imaging, ground-based soil electrical resistivity imaging, and soil moisture sensors were used to quantify the interactions between above and below-ground compartments. Results reveal significant soil-plant co-dynamics. The spatial variation and dynamics of soil water content and electrical conductivity, driven by topographic and soil intrinsic characteristics, correspond to distinct plant types, with highest plant productivity in convergent areas. Plant productivity in heavy snow years benefited from more water infiltration as well as a shallow groundwater table depth. In comparison, low snowpack years with an early first bare-ground date, which are linked to an early increase in plant greenness, imply a short period of saturated conditions that leads to lower average and maximum greenness values during the growing season. Overall, these results emphasize the strong impact of snowpack dynamics, and terrain and subsurface characteristics on the heterogeneity in plant type and seasonal dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ion-beam assisted sputtering of titanium nitride thin films

Abstract Titanium nitride is a material of interest for many superconducting devices such as nanowire microwave resonators and photon detectors. Thus, controlling the growth of TiN thin films with desirable properties is of high importance. This work aims to explore effects in ion beam-assisted sputtering (IBAS), were an observed increase in nominal critical temperature and upper critical fields are in tandem with previous work on Niobium nitride (NbN). We grow thin films of titanium nitride by both, the conventional method of DC reactive magnetron sputtering and the IBAS method, to compare their superconducting critical temperatures $$T_{c}$$ T c as functions of thickness, sheet resistance, and nitrogen flow rate. We perform electrical and structural characterizations by electric transport and x-ray diffraction measurements. Compared to the conventional method of reactive sputtering, the IBAS technique has demonstrated a 10% increase in nominal critical temperature without noticeable variation in the lattice structure. Additionally, we explore the behavior of superconducting $$T_c$$ T c in ultra-thin films. Trends in films grown at high nitrogen concentrations follow predictions of mean-field theory in disordered films and show suppression of superconducting $$T_c$$ T c due to geometric effects, while nitride films grown at low nitrogen concentrations strongly deviate from the theoretical models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tapping into the physiological responses to mistletoe infection during heat and drought stress

Mistletoes are important co-contributors to tree mortality globally, particularly during droughts. In Australia, mistletoe distributions are expanding in temperate woodlands, while their hosts have experienced unprecedented heat and drought stress in recent years. We investigated whether the excessive water use of mistletoes increased the probability of xylem emboli in a mature woodland during the recent record drought that was compounded by multiple heatwaves. Additionally, we continuously recorded transpiration ($T_{SLA}$) of infected and uninfected branches from two eucalypt species over two summers, monitored stem and leaf water potentials ($\Psi $) and used hydraulic vulnerability curves to estimate percent loss in conductivity (PLC) for each species. Variations in weather (vapor pressure deficit, photosynthetically active radiation, soil water content), host species and % mistletoe foliage explained 78% of hourly $T_{SLA}$. While mistletoe acted as an uncontrollable sink for water in the host even during typical summer days, daily $T_{SLA}$ increased up to 4-fold in infected branches on hot days, highlighting the previously overlooked importance of temperature stress in amplifying water loss in mistletoes. The increased water use of mistletoes resulted in significantly decreased host $\Psi _{\rm{leaf}}$ and $\Psi _{\rm{trunk}}$. It further translated to an estimated increase of up to 11% PLC for infected hosts, confirming greater hydraulic dysfunction of infected trees that place them at higher risk of hydraulic failure. However, uninfected branches of Eucalyptus fibrosa F.Muell. had much tighter controls on water loss than uninfected branches of Eucalyptus moluccana Roxb., which shifted the risk of hydraulic failure towards an increased risk of carbon starvation for E. fibrosa. The contrasting mechanistic responses to heat and drought stress between both co-occurring species demonstrates the complexity of host–parasite interactions and highlights the challenge in predicting species-specific responses to biotic agents in a warmer and drier climate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Designing Recyclable Biomass-Based Polyesters (Final Technical Report)

This project successfully demonstrated the production of sustainable polyesters from biomass-derived monomers, offering properties on par with or superior to fossil-based plastics such as LLDPE and PBAT. Using renewable chemical platforms like furfural and HMF, we synthesized over eight novel monomers and nearly 100 aliphatic-aromatic polyesters. Among these, PPeAT and PDDF emerged as highly promising candidates for flexible packaging applications, exhibiting excellent mechanical performance, tunable thermal properties, and enhanced barrier characteristics. Key advances included the development of ring-opening polymerization strategies using cyclic monomer intermediates to produce high molecular weight materials with excellent control. The resulting polymers achieved biobased content of at least 44 wt.%. We also implemented chemical recycling via catalytic methanolysis, recovering over 90% of monomers from PPeAT at semi-pilot scale. This supports a circular economy approach by enabling material reuse without performance loss. Techno-economic and life cycle assessments confirmed that these polyesters can reduce greenhouse gas emissions by up to 50% compared to petroleum-derived alternatives. Favorable mechanical properties may further reduce material demand, and cost competitiveness was demonstrated under specific market conditions. Feedstock cost and environmental footprint were identified as primary levers for future optimization. In summary, the project met or closely approached all technical targets, including biodegradability, strength, and processing compatibility. While minor issues such as color and melting point variation remain, the results strongly support continued development and scale-up of these sustainable materials.

09 BIOMASS FUELS↗
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↗
DOE OSTI2021

NFD Methodology Report

The NetFlow Dynamics (NFD) model was developed for estimating the availability of a commodity supplied by a national- or regional-scale infrastructure following unexpected disruption of one or more of its components. The large scope of the disruptions of interest produce changes in availability lasting days to weeks. Consequently, the model does not resolve daily variations in system state and does not include the long-term processes that cause infrastructures to evolve as assets are added and removed according to owners ’planning decisions. NFD simulates fluid flow, including petroleum and other incompressible fluids, as well as natural gas and other compressible fluids, through pipeline networks characterized by limits on transmission capacity and storage. It was designed to enable efficient exploration of possible transmission system responses to large-scale disruptions lasting for days or longer. The model formulation reflects constraints on transmission and storage capacity imposed by the physical system assets. Those capacity limits are input parameters and are not derived from more basic system properties such as pipeline diameters and compressor power. A system’s response to a large disruption is controlled by operational decisions as well as damage to physical assets. The NFD model formulation allows users to efficiently consider alternative scenarios about the way remaining capacity might be used so that the analysis result appropriately reflects uncertainties about operator response.

24 POWER TRANSMISSION AND DISTRIBUTION↗