Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface tracking”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

FTK: A Simplicial Spacetime Meshing Framework for Robust and Scalable Feature Tracking

In this work, we present the Feature Tracking Kit (FTK), a framework that simplifies, scales, and delivers various feature-tracking algorithms for scientific data. The key of FTK is our simplicial spacetime meshing scheme that generalizes both regular and unstructured spatial meshes to spacetime while tessellating spacetime mesh elements into simplices. The benefits of using simplicial spacetime meshes include (1) reducing ambiguity cases for feature extraction and tracking, (2) simplifying the handling of degeneracies using symbolic perturbations, and (3) enabling scalable and parallel processing. The use of simplicial spacetime meshing simplifies and improves the implementation of several feature-tracking algorithms for critical points, quantum vortices, and isosurfaces. As a software framework, FTK provides end users with VTK/ParaView filters, Python bindings, a command line interface, and programming interfaces for feature-tracking applications. We demonstrate use cases as well as scalability studies through both synthetic data and scientific applications including tokamak, fluid dynamics, and superconductivity simulations. We also conduct end-to-end performance studies on the Summit supercomputer. FTK is open sourced under the MIT license: https://github.com/hguo/ftk.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Label-free time- and space-resolved exometabolite sampling of growing plant roots through nanoporous interfaces

Abstract Spatial and temporal profiling of metabolites within and between living systems is vital to understanding how chemical signaling shapes the composition and function of these complex systems. Measurement of metabolites is challenging because they are often not amenable to extrinsic tags, are diverse in nature, and are present with a broad range of concentrations. Moreover, direct imaging by chemically informative tools can significantly compromise viability of the system of interest or lack adequate resolution. Here, we present a nano-enabled and label-free imaging technology using a microfluidic sampling network to track production and distribution of chemical information in the microenvironment of a living organism. We describe the integration of a polyester track-etched (PETE) nanofluidic interface to physically confine the biological sample within the model environment, while allowing fluidic access via an underlying microfluidic network. The nanoporous interface enables sampling of the microenvironment above in a time-dependent and spatially-resolved manner. For demonstration, the diffusional flux through the PETE membrane was characterized to understand membrane performance, and exometabolites from a growing plant root were successfully profiled in a space- and time-resolved manner. This method and device provide a frame-by-frame description of the chemical environment that maps to the physical and biological characteristics of the sample.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

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

Anomalous nanoparticle surface diffusion in LCTEM is revealed by deep learning-assisted analysis

Significance Liquid cell transmission electron microscopy (LCTEM) is an emerging technique, which enables nanoscale visualization and tracking of single nanoparticles near interfaces with unprecedented spatial resolution. Here, we studied the diffusion of nanoparticles in LCTEM experiments using techniques powered by deep neural networks and statistical tests. We observed two underlying regimes of diffusive behavior which are governed by the interaction of the electron beam, the nanoparticle, the nearby substrate, and the liquid environment. This understanding forms the foundation to use LCTEM for single-nanoparticle tracking for a broad range of nanoparticles, interfaces, and liquids.

Jamali, Vida↗
DOE OSTI2025

Status of the FCC-ee interaction region design

This paper presents a comprehensive overview of the Machine Detector Interface (MDI) design developed for the FCC-ee Feasibility Study. It highlights novel studies related to the lightweight interaction region, including a mechanical model of the vacuum chambers, integration of the vertex detector, the MDI alignment system, and assessments of machine-induced backgrounds. The small beam pipe radius and thickness, as well as the high power to be dissipated, require state-of-the-art mechanical design. The integration of all mechanical elements and detectors is challenging, necessitating careful studies to allow fulfilling conflicting requirements. The optimisation of the machine detector interface against formidable backgrounds is presented.

Beam-induced backgrounds↗
DOE OSTI2021

GEDet: Detecting Erroneous Nodes with A Few Examples

Detecting nodes with erroneous values in real-world graphs re- mains challenging due to the lack of examples and various error scenarios. We demonstrate GEDet, an error detection engine that can detect erroneous nodes in graphs with a few examples. The GEDet framework tackles error detection as a few-shot node classification problem. We invite the attendees to experience the following unique features. (1) Few-shot detection. Users only need to provide a few examples of erroneous nodes to perform error detection with GEDet. GEDet achieves desirable accuracy with (a) a graph augmentation module, which automatically generates synthetic examples to learn the classifier, and (b) an adversarial detection module, which improves classifiers to better distinguish erroneous nodes from both cleaned nodes and synthetic examples. We show that GEDet significantly improves the state-of-the-art error detection methods. (2) Diverse error scenarios. GEDet profiles data errors with a built-in library of transformation functions from correct values to errors. Users can also easily “plug in” new error types or examples. (3) User-centric detection. GEDet supports (a) an active learning mode to engage users to verify detected results, and adapts the error detection process accordingly; and (b) visual interfaces to interpret and track detected errors.

Guan, Sheng↗
DOE OSTI2024

Electric Vehicle Charging Demand in the Chicago Metropolitan Area through 2030

This report outlines the collaborative efforts between Argonne National Laboratory and Exelon in advancing the Agent-Based Transportation Energy Analysis Model (ATEAM). Aligning with ComEd’s beneficial electrification plan, this study developed eight scenarios to access the temporal and spatial distribution of charging load and demand stemming from the widespread adoption of battery electric vehicles (BEV) adoption, augmented public charging infrastructure deployment, and increased multi-unit dwelling (MUD) charging availability. Enhancements to the ATEAM model encompassed the simulation of multiple days of travel behavior, estimation of total public charging infrastructure needs, user interface refinements, and output tracking at both vehicle and charging station levels. The total electricity consumption for residential and public charging to support over 800,000 BEVs in Chicago in 2029 is projected at approximately 10.2 GWh. Enhanced MUD home charging accessibility (70%) amplifies the home charging load in the study area by 1.5% compared to the baseline scenario (10%). The widespread adoption of BEVs reduces peak charging loads, owing to their inclusion across households with diverse income levels, thus fostering a more dispersed charging activity pattern. However, widespread BEV adoption increases the peak home charging load in areas with lower median household incomes, reflecting a higher BEV concentration in these locales and, subsequently, heightened peak charging demands. In the Widespread BEV adoption scenario, fewer census tracts exhibit elevated peak loads for combined home and public charging, indicating a more even distribution of charging demand across the study area. Predominantly, peak loads for combined charging—both home and public— occur between 2 p.m. and 10 p.m. across all scenarios, encompassing the majority of census tracts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Developing a platform for Fresnel diffractive radiography with 1 μ m spatial resolution at the National Ignition Facility

An x-ray Fresnel diffractive radiography platform was designed for use at the National Ignition Facility. It will enable measurements of micron-scale changes in the density gradients across an interface between isochorically heated warm dense matter materials, the evolution of which is driven primarily through thermal conductivity and mutual diffusion. Here we use 4.75 keV Ti K-shell x-ray emission to heat a 1000 μm diameter plastic cylinder, with a central 30 μm diameter channel filled with liquid D 2 , up to 8 eV. This leads to a cylindrical implosion of the liquid D 2 column, compressing it to ~2.3 g/cm 3 . After pressure equilibration, the location of the D 2 /plastic interface remains steady for several nanoseconds, which enables us to track density gradient changes across the material interface with high precision. For radiography, we use Cu He- α x rays at 8.3 keV. Using a slit aperture of only 1 μm width increases the spatial coherence of the source, giving rise to significant diffraction features in the radiography signal, in addition to the refraction enhancement, which further increases its sensitivity to density scale length changes at the D 2 /plastic interface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Rayleigh-Taylor-instability experiments with elastic-plastic materials

A rotating wheel experimental facility was developed to investigate incompressible Rayleigh-Taylor instability in elastic-plastic materials. A soft solid (mayonnaise) was chosen as the elastic-plastic material for experiments; material properties that include shear modulus and yield strength were fully characterized using a vane spindle type rheometer. Initial perturbations of varying amplitudes and wavelengths were generated on the interface of the soft solid using sinusoidal cutting guides. A backlit imaging technique was used in conjunction with a high-speed camera to track the motion of the interface at various phases of the instability. Results for both two- and three-dimensional perturbations were compared to study the acceleration required for instability and the growth after the interface yielded. Exponential growth rates were observed after instability was reached with trends of increasing growth rates for lower initial amplitudes. It was found that the acceleration required for instability increased when initial amplitude and wavelength decreased. Here, three-dimensional interfaces were found to be more stable. For both cases, a decrease in initial amplitude produced a more stable interface that increased the threshold acceleration required for the instability. Critical amplitude conditions for instability were calculated and compared with various analytical models and other experimental results.

42 ENGINEERING↗
DOE OSTI2021

Pore-scale simulation of drying in porous media using a hybrid lattice Boltzmann: pore network model

In this work, a hybrid method coupling a pseudo-potential lattice Boltzmann model (LBM) and a pore network model (PNM) to simulate drying in porous media is proposed. Based on the watershed method, the porous medium is firstly decomposed into pore regions. According to the liquid–vapor phase distribution at a given time, the pore regions are further divided into four pore types, namely two-phase pores where a liquid–vapor interface exists, buffer pores next to the two-phase pores, single-liquid and single-vapor phase pores. The pseudo-potential LBM is used in the two-phase and buffer pores to simulate liquid drying and track the movement of the interfaces, while the single-phase PNM simulations are conducted in the buffer and single-phase pores to simulate vapor or liquid flow. LBM and PNM are coupled in the buffer pores through exchange of boundary information. The hybrid method is applied to simulate liquid drying in a porous medium. The whole-domain LBM simulation is considered as the reference solution to validate the hybrid method. Liquid saturation variation during the drying process and detailed phase and pressure distributions obtained by the two methods match quite well, demonstrating the accuracy of the hybrid method. For the specific case studied, the hybrid method saves more than 60% computational time compared to the whole-domain LBM simulation. In addition, the speedup of the hybrid method becomes more significant for a larger computational domain. In summary, the hybrid method developed in this work combines the accuracy of LBM and the efficiency of PNM to simulate drying in porous media at pore scale and can lead to significant reduction of computation time, thus allowing the pore-scale consideration of drying in larger porous systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

GANpiler

Focal Area(s): Rather than augment or replace physical models with machine learning, we instead preserve the existing models and augment the underlying code with faster surrogate models created by generative adversarial networks (GAN). To leverage the performance optimization, we also propose a runtime system and user interfaces that allow prediction and tracking of accumulated error as well as dynamic, per-process decision making as to which model (if any) to use for each iteration. Science Challenge: This proposal sits at the nexus of two hard problems. First, climate models based on machine learning will be, by their nature, difficult to trust once their predictions begin diverging from the consensus. Second, compilers and hardware have been making only incremental performance gains for decades. GPGPUs have provided a welcome performance boost for codes that can take advantage of them, but there is no similar technology on the horizon to provide the next performance leap.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗
DOE OSTI2019

Use of In Situ Acoustic Backscatter Systems to Characterize Spent Nuclear Fuel and its Separation in a Thickener - 19288

To aid in the transportation, dewatering and storage of radioactive UK legacy waste sludges at Sellafield Ltd., the application of an Acoustic Backscatter System (ABS) was investigated in this study, as a method for monitoring suspended solids concentration and consolidation of cohesive and aggregated sludges (that form the basis of many wastes at Sellafield). An ABS represents a particularly promising technology, as they can be used with minimal intrusion, to measure both particle size and concentration depending on prior system knowledge, through inversion of the return echo voltage response of pulse-echo signals in the MHz range. To assess the application for continuous, high concentration (∼30 % v/v) thickened wastes, an ABS was utilised to characterise a pseudo-steady state laboratory thickener over hours, using flocculated calcite as a representative test material. Measurements were taken both horizontally through a consolidated bed using a 1 MHz transducer and vertically downwards with a 2 MHz transducer to measure sediment concentration above the settled bed and track bed height via the interface reflection. Characterisation of flocculated aggregate diameters was initially determined to be ∼510 μm, using an in-line Focussed Beam Reflectance Measurement (FBRM) device. The submerged 2 MHz probe indicated that there were only limited changes to the dispersion concentration or size over time. Changes in settled bed density, due to increasing bed height and changes in the underflow flow rate, were qualitatively measured via differences to the acoustic signal attenuation over time. An initial increase in attenuation with bed height was observed as the bed densified due to local compression. Once an equilibrium bed height was reached, the acoustic attenuation remained fairly constant with time, although the bulk underflow density was reduced with the increase in underflow rate. The solids residence time in the bed, determined from a transitive volume-balance model, was found to be ∼3200 s for pseudo-steady state operation at a maintained bed height of ∼0.175 m. Results of the study highlight the potential the ABS as a remote process monitoring tool for both relatively dilute suspensions and concentrated thickened mineral sludges, with potential applications across waste processing sites at Sellafield. (authors)

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

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Electronic Excitation‐Driven β‐Ga 2 O 3 Metastability Transformation and Self‐Organization Mechanism: β→κ/γ/δ Phases

Here, irradiation-driven multiphase self-organization presents emergent opportunities for the customization of nanoscale engineering properties, dynamically tuning strain-field distributions and interfacial electronic structures. Responding to intense electronic excitation-induced energy deposition, the dominant phase transformations, with varying Gibbs free energy $Δ,G^o_f$ are confirmed as β → κ → γ → δ that are located in specific microregions for Gallium (III) oxide (Ga 2 O 3 ), as follows: (i) Surface-localized interstitial accumulation under compressive stress triggers β → δ via semi-coherent interface formation. (ii) Tensile stress within latent tracks drives vacancy-mediated oxygen layer truncation (4/12 periodicity along ⟨0001⟩), stabilizing coherent 4H (ABCB) κ and 3C (ABC) β (ABC) interfaces through strain-compensated octahedral distortion. (iii) Screw dislocation-mediated lattice relaxation induces β → γ via cation disordering (Ga 3+ occupancy at β-interstitial sites), forming metastable spinel γ with mixed occupancy across 16d/8a Wyckoff sites. Irradiation-driven β-Ga 2 O 3 →κ/γ/δ transitions, as mechanistically revealed via inelastic thermal spike (i-TS) calculations and molecular dynamics simulations, induce defect-mediated nonlinear photoresponse, critical for optoelectronic engineering.

Han, Xinqing [Shandong Univ., Jinan (China)]↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗