Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “logs”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The effect of composition, temperature and pressure on the oxidation state and coordination environment of copper in silicate melts

Copper is a redox variable element that may occur as Cu 0 , Cu + , and Cu 2+ in the Earth’s crust. The oxidation state will affect its partitioning between coexisting minerals, melts and fluids and hence its behaviour in magmatic processes. Copper bearing silicate glasses were quenched from melts with 23 synthetic compositions (19 CaO-Na 2 O-MgO-Al 2 O 3 -SiO 2 (CNMAS), two “granites” containing K 2 O ± H 2 O, and Fe-bearing “MORB” and “andesite”) equilibrated at oxygen fugacities (fO 2 ), expressed in log units relative to the fayalite-magnetite-quartz (FMQ) buffer, ranging from −0.7 to 14, temperatures from 900 to 1500 °C and pressures from 0 to 2.5 GPa. Cu K-edge X-ray absorption near edge structure (XANES) spectra were recorded from the glasses and a pre-edge feature in the XANES spectra was found to scale with the proportion of Cu + . Cu + /ΣCu (where ΣCu = Cu + + Cu 2+ ) was quantified by fitting the intensity of the pre-edge feature as a function of fO 2 to the thermodynamically expected relationship. Cu + /ΣCu was found to only weakly depend on melt composition, with more basic melts (e.g., basalts rather than granites) preferentially stabilising Cu 2+ . Increasing temperature stabilises Cu + , while increasing pressure had little effect on Cu + /ΣCu in CNMAS melts but preferentially stabilised Cu 2+ in granite melts. Cu + /ΣCu can be predicted in silicate melts by the empirical equation: log(Cu 2+ /Cu + ) = 0.25(ΔFMQ + 8.58 − 25050/T + 940P/T – 0.02P) − 4.73 + 5400/T + 1.99Λ + (280P – 90P2) where T is temperature in K, P is pressure in GPa and is the optical basicity of the composition. The effects of fO 2 , melt composition, temperature and pressure on Cu + /ΣCu indicate that Cu + will be the dominant oxidation state in terrestrial silicate melts (e.g., Cu + /ΣCu = 99% in an andesitic melt at 900 °C, 1 GPa and ΔFMQ = 1). The electron exchange reaction Cu 2+ + Fe 2+ → Cu + + Fe 3+ occurs on cooling and, given the abundance of Fe in natural melts, the oxidation state of Cu in natural glasses is unlikely to correspond to that in the original melt.

Copper↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2021

A comprehensive analysis of transient pressure and rate data from CO 2 storage projects in a depleted pinnacle reef oil field complex, Michigan, USA

Pressure and rate data are commonly recorded as part of a basic monitoring program in CCS projects. This paper discusses the application of multiple analytical techniques to interpret pressure and rate transient data from CO 2 injection and storage operations. The techniques of interest, i.e., injection-falloff analysis, injectivity/productivity index analysis and pressure pulse arrival time analysis, are commonly used in the oil and gas industry to assess reservoir properties, but not well known in the CCS literature (especially the last two). Injection-falloff analysis involves log-log pressure derivative plotting for the falloff data and history-matching of the entire injection-falloff sequence to determine permeability. In the injectivity/productivity index analysis, rate-normalized pressure buildup is plotted against material balance time or ratio of cumulative injection to injection rate to determine the injectivity index (ratio of injection rate to stabilized pressure buildup) which can be related to the permeability-thickness product. The arrival time analysis identifies the arrival of a pressure disturbance (~0.1 psi change from ambient) to determine the hydraulic diffusivity from which permeability can be estimated. The applicability of these techniques is demonstrated via illustrative examples from multiple wells in different pinnacle carbonate reefs undergoing CO 2 -EOR in Northern Michigan. The paper ends with a discussion of the relative merits of each interpretive technique, as well as recommendations that could be useful for other field projects.

42 ENGINEERING↗
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Consequences of product inhibition in the quantification of kinetic parameters

While the potential for product inhibition in catalytic reactions is well known, the impact of neglected inhibition on measured kinetic parameters is often overlooked. The presence of product inhibition, most often caused by the competitive adsorption of products with reactants on catalytic active sites, is difficult to determine a priori for an arbitrary catalytic system. The significance of product inhibition relies on the concentration of the products, their adsorption thermodynamics on catalytically relevant sites, and process parameters such as temperature and pressure. When inhibition is significant, however, apparent activation energies and reaction orders vary from the differential-reactor apparent activation energy by a factor of 1/(1 - δ), where δ is the total inhibition order (e.g., the factor (1 - δ) = 1.6 for a system with product inhibition of -0.6 order). This is illustrated here with the kinetics of NO oxidation over Cu ion clusters (Cu x O y ) in Cu-SSZ-13, for which the product NO 2 inhibits the forward reaction. Furthermore, in the presence of inhibition, when only reactants are fed to a flow reactor or placed in a batch reactor, there is often no practical conversion that is low enough to guarantee differential behavior. Inclusion of products in the feed solves this problem, allowing accurate determination of kinetic parameters such as apparent reaction orders and activation energies. We also demonstrate that evaluation of the necessity of co-feeding products to assure measurement of differential-reactor data in a given catalytic system is straightforward from a plot of the log of the rate (or conversion) versus the log of the space time in a flow reactor or elapsed time in a batch reactor. We encourage inclusion of this test in all kinetic analyses that are reasonably approximated by power law rate expressions.

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

Post-remediation geophysical assessment: Investigating long-term electrical geophysical signatures resulting from bioremediation at a chlorinated solvent contaminated site

There is a growing need to assess long-term impacts of active remediation strategies on treated aquifers. A variety of biogeochemical alterations can result from interactions of the amendment with the aquifer, conceivably leading to a geophysical footprint of the long-term alteration of an aquifer. This concept of post-remediation geophysical assessment was investigated in a shallow, chlorinated solvent-contaminated aquifer six to eight years after amendment delivery. Surface resistivity imaging and cross-borehole resistivity and induced polarization (IP) imaging were performed on a transect that spanned treated and untreated zones of the aquifer. Established relationships between IP parameters and surface electrical conductivity were used to predict vertical profiles of electrolytic conductivity and surface conductivity from the inverted cross-borehole images. Aqueous geochemistry data, along with natural gamma and magnetic susceptibility logs, were used to constrain the interpretation. The electrical conductivity structure determined from surface and borehole imaging was foremost controlled by the electrolytic conductivity of the interconnected pore space, being linearly related to fluid specific conductance. The electrolytic conductivity (and thus the conductivity images alone) did not discriminate between treated and untreated zones of the aquifer. In contrast, inverted phase angles and surface conductivities did discriminate between treated and untreated zones of the aquifer, with the treated zone being up to an order of magnitude more polarizable in places. Supporting aqueous chemistry and borehole logging datasets indicate that this geophysical footprint of the long-term impact of the remediation on the aquifer is most likely associated with the formation of polarizable, dispersed iron sulfide minerals. The study opens the door to the possibility of employing time-lapse electrical geophysical measurements to assist with long-term environmental stewardship of legacy sites undergoing active remediation.

bioremediation↗
DOE OSTI2025

Shock compression of additively manufactured high solids loaded polymer composites with aligned porosity

The effects of aligned porosity on the shock-compression response of additively manufactured (AM) composite structures with collinear and log-cabin geometries are investigated. The composites contain two inorganic and two organic particle populations and a UV-curable polymer binder. X-ray phase contrast imaging (X-PCI) is used as an in-material diagnostic in plate-impact experiments performed at the Dynamic Compression Sector of the Advanced Photon Source at the Argonne National Laboratory. Macroscale shock and particle velocities and microscale pore collapse velocities are measured from X-PCI images obtained from samples impacted in various orientations relative to the AM build directions. The linear shock and particle velocity relation for the structures with collinear build geometry shows directional differences, with the build direction exhibiting the highest slope and the filament print direction having the shallowest slope. Composites with log-cabin geometry, due to their slightly lower porosity and smaller sized pores, show higher shock velocities than the collinear structures at the same particle velocity. Microscale directionality effects are also observed, with the pore collapse velocity dependent on the pore configuration. Specifically, the pore collapse velocity is highest for pores elongated along the impact direction which is the case for the colinear structures impacted in the filament print direction, which is consistent with the directionality trend observed in the data for the shock and particle velocity equation of state. The pore collapse velocity approaches ∼80 % of the shock velocity and is several times the average bulk particle velocity, suggesting a highly hydrodynamic mechanism of pore collapse. Overall, both the macroscale response as measured by the shock and particle velocity relationship, and the microscale response as measured by the pore collapse velocity, are found to be anisotropic and dependent on the AM composite structure.

Wagner, K. B. [Georgia Institute of Technology, At↗
DOE OSTI2020

Theoretical power-law relationship between permeability and formation factor

Permeability k and formation factor F are important petrophysical properties that have broad applications to reservoir engineering. Although well-log analysis has advanced reservoir characterizations, direct measurement of k is still not available. Accordingly, resistivity logs and more particularly formation factor have been frequently used to estimate k. In this study, we present a theoretical framework to relate k to F. We use power-law scaling from continuum percolation theory to develop a theoretic relationship between k and F, which is similar in form to a published power-law equation with empirical coefficients. Our result gives a physical interpretation to this empirical relationship. We also compare power-law scaling laws from percolation theory with data derived from experimental measurements and numerical simulations in the literature. Finally, using the Curve Fitting toolbox of MATLAB, the parameters of the power-law equations are determined, and the significance of the non-universal power laws are discussed.

42 ENGINEERING↗
DOE OSTI2020

Measurement of groomed jet substructure observables in p+p collisions at $\sqrt{s} =$ 200 GeV with STAR

In this letter, measurements of the shared momentum fraction ( z g ) and the groomed jet radius ( R g ), as defined in the SoftDrop algorithm, are reported in p+p collisions at s = 200 GeV collected by the STAR experiment. These substructure observables are differentially measured for jets of varying resolution parameters from R = 0.2 – 0.6 in the transverse momentum range 15 < p T , jet < 60 GeV/ c . These studies show that, in the p T , jet range accessible at s = 200 GeV and with increasing jet resolution parameter and jet transverse momentum, the z g distribution asymptotically converges to the DGLAP splitting kernel for a quark radiating a gluon. The groomed jet radius measurements reflect a momentum-dependent narrowing of the jet structure for jets of a given resolution parameter, i.e., the larger the p T , jet , the narrower the first splitting. For the first time, these fully corrected measurements are compared to Monte Carlo generators with leading order QCD matrix elements and leading log in the parton shower, and to state-of-the-art theoretical calculations at next-to-leading-log accuracy. We observe that PYTHIA 6 with parameters tuned to reproduce RHIC measurements is able to quantitatively describe data, whereas PYTHIA 8 and HERWIG 7, tuned to reproduce LHC data, are unable to provide a simultaneous description of both z g and R g , resulting in opportunities for fine parameter tuning of these models for p+p collisions at RHIC energies. We also find that the theoretical calculations without non-perturbative corrections are able to qualitatively describe the trend in data for jets of large resolution parameters at high p T , jet , but fail at small jet resolution parameters and low jet transverse momenta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Rheological behavior and fiber spinning of polyacrylonitrile (PAN)/Carbon nanotube (CNT) dispersions at high CNT loading

Rheological studies play an important role in polymer processing including fiber spinning. In the current work, rheological behavior has been studied for polyacrylonitrile (PAN)/carbon nanotube (CNT) dispersions in dimethyl formamide (DMF), where CNT loading is as high as 15 wt% of the total solids (polymer + CNT). The presence of CNTs increased the elastic-like and shear thinning behavior of the dispersions. The terminal slope of the log G' versus log G" curve, which was used to characterize the homogeneity of solution or dispersion, is dependent on the CNT concentration but independent of temperature between 0 and 60 °C. PAN/CNT-10 (90 wt% PAN +10 wt% CNT) and PAN/CNT-15 were gel spun at different conditions, with small diameter capillary at high shear rate and large diameter capillary at relatively low shear rate. Additionally, for PAN/CNT-10, fiber spinning with lower molecular weight polymer showed increased continuous jetting time when compared to spinning with higher molecular weight polymer. In addition, the rheological behavior of PAN/cellulose nanocrystal (CNC) dispersion with 20–60 wt% CNC was also compared to that of the PAN/CNT dispersions. PAN/CNC and PAN/CNT share similarities in both rheological behavior and fiber spinning at high filler loading.

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

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

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

Plutonium Speciation and Oxidation State Distributions in the Presence of Citrate

We explored the speciation and kinetics of the Pu(VI)-citrate and Pu(III)-citrate systems (pH m = 2.5–11.0, I = 0.1 M NaCl, T = 23 °C, O 2 (g) < 2 ppm) using ultraviolet–visible–near-infrared (UV–vis–NIR) spectrophotometry, solvent extraction, and PHREEQC modeling. Formation constants were determined for PuO 2 (HcitH)(aq) (log K° 1,1 = 1.09 ± 0.05) and PuO 2 (HcitH)(citH) 3– (log K° 1,2 = -0.20 ± 0.07), and evidence for (PuO 2 ) m (citH –k ) n (OH) x 2m(3+k)n–x was identified under alkaline conditions. Pu(VI) species were found to be less stable in the presence of citrate than in the absence of citrate (t ≤ 168 days); the rate of reduction increased with increasing pH. Further, the direct reduction of Pu(VI) to Pu(IV) was required to fit experimental data in the presence of citrate but did not improve the fit for Pu in the absence of citrate. We also observed increased Pu(III) stability in the presence of citrate (t ≤ 293 days), with higher concentrations of Pu(III) favored at lower pH. Finally, we provide evidence of a radiolysis-driven mechanism for the citrate-mediated reduction of plutonium that involves electron transfer from the oxidative breakdown of citrate. Our work highlights the need to investigate the redox effect of organic ligands on plutonium oxidation states under repository-relevant conditions.

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

Reining in Radium for Nuclear Medicine: Extra-Large Chelator Development for an Extra-Large Ion

Targeted α therapy (TAT) of soft-tissue cancers using the α particle-emitting radionuclide 223 Ra holds great potential because of its favorable nuclear properties, adequate availability, and established clinical use for treating metastatic prostate cancer of the bone. Despite these advantages, the use of 223 Ra has been largely overshadowed by other α emitters due to its challenging chelation chemistry. A key criterion that needs to be met for a radionuclide to be used in TAT is its stable attachment to a targeting vector via a bifunctional chelator. The low charge density of Ra 2+ arising from its large ionic radius weakens its electrostatic binding interactions with chelators, leading to insufficient complex stability in vivo. In this study, we synthesized and evaluated macropa-XL as a novel chelator for 223 Ra. It bears a large 21-crown-7 macrocyclic core and two picolinate pendent groups, which we hypothesized would effectively saturate the large coordination sphere of the Ra 2+ ion. The structural chemistry of macropa-XL was first established with the nonradioactive Ba 2+ ion using X-ray diffraction and X-ray absorption spectroscopy, which revealed the formation of an 11-coordinate complex in a rare anti pendent-arm configuration. Subsequently, the stability constant of the [Ra(macropa-XL)] complex was determined via competitive cation exchange with 223 Ra and 224 Ra radiotracers and compared with that of macropa, the current state-of-the-art chelator for Ra 2+ . A moderate log K ML value of 8.12 was measured for [Ra(macropa-XL)], which is approximately 1.5 log K units lower than the stability constant of [Ra(macropa)]. This relative decrease in Ra 2+ complex stability for macropa-XL versus macropa was further probed using density functional theory calculations. Additionally, macropa-XL was radiolabeled with 223 Ra, and the kinetic stability of the resulting complex was evaluated in human serum. Furthermore, although macropa-XL could effectively bind 223 Ra under mild conditions, the complex appeared to be unstable to transchelation. Collectively, this study sheds additional light on the chelation chemistry of the exotic Ra 2+ ion and contributes to the small, but growing, number of chelator development efforts for 223 Ra-based TAT.

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

A Mesopore-Confined and Graphene Oxide-Localized Ruthenium Catalyst Increases Rates of Mid-Chain Polyolefin Hydrogenolysis

A catalyst architecture with mesoporous silica coating Ru nanoparticles on reduced graphene oxide (mSiO 2 /Ru/rGO) localizes the 2 nm Ru solely at the closed bottoms of 2.9 nm diameter mesopores. mSiO 2 /Ru/rGO catalyzes the rapid, selective hydrogenolysis of polyolefins at wax formation rates (ν wax ) up to 1700 g wax ·g Ru –1 ·h –1 and a turnover frequency (TOF) for C–C bond cleavage of 130 ± 8 min –1 . The ν wax and TOF for mSiO 2 /Ru/rGO are 23× and 16× those of nonporous Ru/rGO, respectively, indicating that faster chain cleavage is also more selective inside mesopores. The methane yield from mSiO 2 /Ru/rGO is ca. 30% of the value obtained from Ru/rGO. Methane decreases, and wax selectivity improves with narrow-pore (2.3 nm) mSiO 2 /Ru/rGO. The reaction rate decreases at lower and higher H 2 pressures from its maximum at 37 bar. Log(TOF)-log(P H2 ) plots reveal rates ∝ [P H2 ] 2.8 or [P H2 ] −2.0 in the lower or higher pressure ranges, respectively. Corresponding plots for Ru/C give rates ∝ [P H2 ] 1 or [P H2 ] −2.4 . Ru is hydrocarbon-covered at low pressure; thus, the higher H 2 order for mSiO 2 /Ru/rGO indicates that mesopores increase the density of cleavable unsaturated moieties on the Ru surface. The lower H 2 order on Ru/C implies a lower density of cleavable groups because saturated segments occupy a larger fraction of the Ru surface. Higher surface occupancy by the hydrocarbon reactant correlates with more methane formation. Therefore, pore confinement, narrower pore diameter in mSiO 2 /Ru/rGO, or higher H 2 pressure leads to lower methane yields. Ru nanoparticles in mSiO 2 /Ru/rGO maintain equivalent activity and selectivity over five recycling tests.

Catalysts↗
DOE OSTI2020

Physical Properties and Gas Hydrate at a Near-Seafloor Thrust Fault, Hikurangi Margin, New Zealand

The Papaku Fault Zone, drilled at International Ocean Discovery Program (IODP) Site U1518, is an active splay fault in the frontal accretionary wedge of the Hikurangi Margin. In logging-while-drilling data, the 33-m-thick fault zone exhibits mixed modes of deformation associated with a trend of downward decreasing density, P-wave velocity, and resistivity. Methane hydrate is observed from ~30 to 585 m below seafloor (mbsf), including within and surrounding the fault zone. Hydrate accumulations are vertically discontinuous and occur throughout the entire logged section at low to moderate saturation in silty and sandy centimeter-thick layers. In this paper, we argue that the hydrate distribution implies that the methane is not sourced from fluid flow along the fault but instead by local diffusion. This, combined with geophysical observations and geochemical measurements from Site U1518, suggests that the fault is not a focused migration pathway for deeply sourced fluids and that the near-seafloor Papaku Fault Zone has little to no active fluid flow.

58 GEOSCIENCES↗
DOE OSTI2023

Linking Stream Chemistry to Subsurface Redox Architecture

As drinking-water scarcity grows worldwide, we need to improve predictions of the quantity and quality of our water resources. An overarching problem for model improvement is that we do not know the geological structure of aquifers in sufficient detail. In this work, we demonstrate that mineral-water reactions imprint structure in the subsurface that impacts the flow and transport of some chemical species. Specifically, pyrite, a ubiquitous mineral, commonly oxidizes and depletes in the upper layers of the weathering profile in most humid watersheds, only remaining at depths of meters. We hypothesize that variations in concentrations (C) of pyrite-derived sulfate released into rivers as a function of discharge (q) reflect the rate-limiting step and depth of pyrite-oxidizing layers. We found that log C - log q behaviors thus differ in small and large watersheds in the Susquehanna River Basin as well as in selected watersheds in the Western United States. Although coal mining changes pyrite oxidation from closed to open system with respect to O 2 , patterns in stream chemistry as a function of discharge are consistent with deep and shallow pyrite oxidation zones in small and large watersheds respectively. Therefore, understanding the subsurface patterns of mineral reactions and how they affect the architecture of aquifers will elucidate patterns of changing river chemistry and our ability to manage water resources in the future under accelerated land use and climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Surface Energy Dynamics and Canopy Structural Properties in Intact and Disturbed Forests in the Southern Amazon

The Brazilian Amazon has been a focus of land development with large swaths of forests converted to agriculture. Forest degradation by selective logging and fires has accompanied the agricultural frontier and has resulted in significant impacts on Amazonian ecosystems. Changes in forest structure resulting from forest disturbances have large impacts on the surface energy balance, including on land surface temperature (LST) and evapotranspiration (ET). This study's objective is to assess the effects of forest disturbances on water fluxes and forest structure in a transitional forest site in the Southern Amazon. We used ET and LST products from MODIS and Landsat 8 and GEDI-derived forest structure data to address our research questions. We found that disturbances induced seasonal water stress, more pronounced in croplands/pastures than in forests (differences up to 20% in the dry season), and more pronounced in second-growth and recently burned areas than in logged and intact forests (differences up to 12% in the dry season). Moreover, ET and LST were negatively related, with more consistent relationships across disturbance classes in the dry season (R 2 : 0.41–0.87) than in the wet season (R 2 : 0.18–0.49). Forest and cropland or pasture classes showed contrasting relationships in the dry season. Finally, we found that forest structure exhibited stronger relationships with ET and LST in the most disturbed forests (R 2 : 0.01–0.43) than in the least disturbed forests (R 2 < 0.05). Our findings help to elucidate degraded forests functioning under a changing climate and to improve estimates of water and energy fluxes in Amazonian degraded forests.

54 ENVIRONMENTAL SCIENCES↗