Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Processes”

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 37 records · Page 2

Autotrophic Ecosystems on the Early Earth

Ophiolite sequences, sections of lower oceanic crust and upper mantle that have been thrust onto continental craton, are located in northern and central California and provide easily accessible areas that serve as good analogs for similar, more extensive areas of the early Earth. We have begun investigating and characterizing these sites in order to understand better the processes that may be responsible for the water chemistry, mineralogy and biology of similar environments on the early Earth. The geophysical and geochemical processes in these terranes provide niches for unique communities of extremeophiles and likely provide a good analog to the location that first gave rise to life on Earth. The ophiolites found in northern and central California include the Trinity, Josephine, Coast Range and Point Sal, all of which are approximately 160 million years old. Fluids from serpentinizing springs are generally alkaline with high pH and H2 contents, indicating that the mafic rock compositions control the fluid composition through water-rock reactions during relatively low-grade hydrothermal processes. There are significant amounts of primary mineralogy remaining in the rocks, meaning that substantial alteration processes are still occurring in these terranes. The general reaction for serpentinization of olivine is given by one of the authors. olivine + H2O = serpentine + brucite + magnetite + H2. We have analyzed the mineralogical composition of several rock samples collected from the Coast Range Ophiolite near Clear Lake, CA by electron microprobe. The remnant primary mineralogy is fairly urnform in composition, with an olivine composition of Fo(sub 90), and with pyroxene compositions of En(sub 90) for orthopyroxene and En(sub 49)Wo(sub 48)Fs(sub 03) for the clinopyroxene. Other primary phases observed include chromites and other spinels. Examination of petrographic thin sections reveals that serpentinization reactions have occurred in these locations. The serpentine resulting from aqueous alteration of olivine resides in veins that are see to cross cut the primary mineral grains. There are several generations of alteration products, comprised mostly of serpentines that are magnesium rich, with magnetite, brucite and carbonates observed as accessory minerals. The formation of carbonates can be taken to indicate the presence of CO2 in the altering fluids. We collected samples from a spring in the Coast Range Ophiolite in order to determine whether the geochemical environment serves as a habitat for chemotrophic microorganisms. DNA was extracted from the sediment samples and the 16s rRNA gene was PCR amplified using universal Archaeal primers. Denaturing gradient gel electrophoresis (DGGE) was used to determine the community of Archaea thriving in these samples. Our results indicate that there were 8 different genera of Archaea from a single sample. A sequence was obtained from one of these eight. The sequence is of an organism similar to Halorubrum tibetense, and alkalophilic Archaeon. This result suggests that these environments are likely hosts for communities of organisms that are adapted for the unique chemistry provided by the alkaline spring.

Schulte, M.↗

The Origin of Carbon-Bearing Volatiles in a Continental Hydrothermal System in the Great Basin: Water Chemistry and Isotope Characterizations

Hydrothermal systems on Earth are active centers in the crust where organic molecules can be synthesized biotically or abiotically under a wide range of physical and chemical conditions [1-3]. Not only are volatile species (CO, CO2, H2, and hydrocarbons) a reflection of deep-seated hydrothermal alteration processes, but they also form an important component of biological systems. Studying carbon-bearing fluids from hydrothermal systems is of specific importance to understanding (bio-)geochemical processes within these systems. With recent detection of methane in the martian atmosphere [4-7] and the possibility of its hydrothermal origin [8, 9], understanding the formation mechanisms of methane may provide constraints on the history of the martian aqueous environments and climate.

Fu, Qi↗

The Inner Solar System's Habitability Through Time

Earth, Mars, and Venus, irradiated by an evolving Sun, have had fascinating but diverging histories of habitability. Although only Earth’s surface is considered to be habitable today, all three planets might have simultaneously been habitable early in their histories to microbial life. We consider how physical processes that have operated similarly or differently on these planets determined the creation and evolution of their atmospheres and surfaces over time. These include the geophysical and geochemical processes that determined the style of their interior dynamics and the presence or absence of a magnetic field; the surface-atmosphere exchange processes that acted as a source or sink for atmospheric mass and composition; the Sun-planet interactions that controlled escape of gases to space; and the atmospheric processes that interacted with these to determine climate and habitability. The divergent evolutions of the three planets provide an invaluable context for thinking about the search for life outside the Solar System.

Anthony D. Del Genio↗

Geochemical Changes in Response to CO 2 Injection in a CO 2 -EOR Complex in Northern Michigan (Volume II.1)

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the geochemical monitoring program under MRCSP to use stable and radiogenic isotope geochemistry in concert with analysis of general geochemical parameters of fluids and gases and analysis of core samples to determine geochemical processes occurring in the reef structure because of CO 2 injection. Specifically, brine and gas samples were collected and analyzed to determine changes occurring between reefs prior to and following CO 2 injection. The analytical results for general geochemical parameters were modeled with chemical equilibrium models to determine if the injection of CO 2 resulted in the mineral dissolution or precipitation. Finally, core samples were collected and analyzed to determine if there was evidence of dissolution features or mineral precipitation.

01 COAL, LIGNITE, AND PEAT↗

Temporal evolution of Pu and Cs sediment contamination in a seasonally stratified pond

There remains a lack of knowledge regarding ecosystem transfer, transport processes, and mechanisms, which influence the long-term mobility of Pu-239 and Cs-137 in natural environments. Monitoring the distribution and migration of trace radioisotopes as ecosystem tracers has the potential to provide insight into the underlying mechanisms of geochemical cycles. This study investigated the distribution of anthropogenic radionuclides Pu-239 and Cs-137 along with total organic carbon, iron, and trace element in contaminated sediments of Pond B at the Savannah River Site (SRS). Pond B received reactor cooling water from 1961 to 1964, and trace amounts of Pu-239 and Cs-137 during operations. Our study collected sediment cores to determine concentrations of Pu-239, Cs-137, and major and minor elements in solid phase, pore water and an electrochemical method was used on wet cores to determine dissolved elemental concentrations. More than 50 years after deposition, Pu-239 and Cs-137 in sediments are primarily located in the upper 5 cm in area where deposition of particulate-bound contaminants was prevalent and located between 5 and 10 cm in areas of high sedimentation, showing a limited migration of Pu-239 and Cs-137. A Factor analysis demonstrated different sediment facies across the pond resulting in a range of geochemical processes controlling accumulation of Pu and Cs. Highest concentrations appear to be controlled by particulate input from the influent canal, dominated by clay, silt, and sand minerals bearing Fe. Elevated Pu-239 in the sediments were observed in areas with high organic matter and higher deposition rate relative to the Pond B system near the outlet indicating strong association of Pu with OM and particulates. Therefore, organic matter cycling likely plays a role in Pu redistribution between sediment and overlying pond water, and deposition in organic rich sediments accumulating near the outlet. Though Pu appears to have been distributed throughout the pond, Cs-137 concentrations remained the highest near the influent canal.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Abundances and Distributions of Volatile Elements (Na, S, Cl, K) Within Lunar Volcanic Glass Spheres

This research developed and applied high-precision analytical methods (electron microprobe) to determine the major- and minor-element abundances (Si, Ti, A], Cr, Fe, Mn, Mg, Ca, Na, K) in three suites of Apollo 14 volcanic green glasses. Since the precision of these data was a factor of 5-10 better than previously published results, the data provided the strongest constraints yet achieved on the geochemical processes that operated at depths less than 400 kilometers inside the Moon. The results of this work were submitted to the 29th Lunar and Planetary Science Conference, formed the basis for a Master's Thesis, and were published in the top-tier geochemical journal. Additional high-precision data on the trace abundances of sodium (Na), sulfur (S), and potassium (K) in the Apollo 17 volcanic orange glasses (74220) were also acquired during this grant-period. These data show diffusion profiles within individual 200-micron diameter, glass spherules that were apparently generated during the fire-fountain eruption of this magma in the presence of a transient vapor phase. Although these results have not yet been published, publication is planned in 2001 after additional numerical modeling has been completed. This work has also catalyzed on-going collaborative work with Professor Timothy Grove (MIT) to analyze and understand the origin of the Apollo 15 volcanic green glasses. That work will be completed for publication in 2001.

Delano, John W.↗

Mineralogy of dark clasts in primitive versus differentiated meteorites

The presence of dark lithic clasts within meteorites can provide information concerning asteroidal regolith processes, the extent of interactions between asteroids, and the relationship between meteorite types, micrometeorites, and interplanetary dust particles. Accordingly, we have been seeking and characterizing dark clasts found within carbonaceous chondrites, unequilibrated ordinary chondrites, howardites, and eucrites. We find that unequilibrated chondrites in this study contain fine-grained, anhydrous unequilibrated inclusions, while the howardites often contain inclusions from geochemically processed, hydrous asteroids (type 1 and 2 carbonaceous chondrites). Eucrites and howardities contain unusual clasts, not easily classified.

Zolensky, M. E.↗

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗

Compilation and Evaluation of Data from Groundwater Impounded within the U12n and U12t Tunnels, Rainier and Aqueduct Mesas, Nevada National Security Site

The Department of Defense (DOD) Defense Threat Reduction Agency (DTRA) and the U.S. Department of Energy (DOE) Office of Environmental Management Nevada Program (EM NV) conducted routine physical and geochemical monitoring of water impounded behind impermeable structures placed within the U12n and U12t tunnels on the Nevada National Security Site (NNSS) (formerly the Nevada Test Site [NTS]). Desert Research Institute (DRI) conducted a semiquantitative evaluation of pressure data to determine the extent to which the tunnels have been inundated and a qualitative evaluation of monitored geochemical parameters to determine if trends exist that were indicative of hydrologic and/or geochemical processes operating behind the tunnel plugs. Data were also evaluated to determine if monitored radiological parameters exceeded primary standards associated with the Environmental Protection Agency’s (EPA) Safe Drinking Water Act (SDWA) (40 C.F.R. § 141.66).

54 ENVIRONMENTAL SCIENCES↗

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗

Quantifying the Impact of Electric Fields on the Local Structure and Migration of Potassium Ions at the Orthoclase (001) Surface

Developing an understanding of the response of mineral/water interfaces to applied electric fields is central to detecting and interpreting signatures of interfacial processes in the subsurface. Here we focus density functional theory calculations on understanding K + cation binding and diffusion across the (001) surface of orthoclase feldspar under various applied electric fields with and without surface hydration. The calculations reveal how water ligands labilize surface K + cations for diffusion while also increasing their sensitivity to electric field effects on the binding energy at different surface sites. The calculations also show how the direction and strength of the electric field systematically affect surface cation mobility, sorption, and hydration behavior. Specifically, electric fields directed toward the surface increase the potassium binding energies, facilitate K + diffusion across multiple binding sites, and increase the propensity of surface K + ions to lose their hydration shells. Furthermore, the findings help fill a major knowledge gap into the impact of electric fields on mineral/water interface structure and dynamics more generally, featuring in this case a commonly found type of feldspar involved in a multitude of atmospheric and geochemical processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗

Development of Cr, Se, U, Sb, and Te isotopes as indicators of redox reactions, contaminant fate, and contaminant transport in aqueous systems: A review

Cr, Se, U, Sb, and Te are toxic, redox-active elements that are more mobile and environmentally problematic in their oxidized forms, and less mobile and bioavailable in their reduced forms. This chapter reviews the development of Cr, Se, U, Sb, and Te isotope ratio measurements as new indicators of redox reactions and contaminant migration. Reliable analytical methods exist, but are still evolving. Understanding of isotopic fractionation induced by various (bio)geochemical processes has been explored in dozens of publications, yet is far from complete: Reduction reactions, the major driver of isotopic variation, have been relatively well studied. However, the magnitude of fractionation is variable and the systematics of that variation are still being explored. Isotopic fractionation induced by oxidation reactions is not well understood. Non-redox reactions, which involve smaller changes in bonding of these elements, tend to induce less isotopic fractionation, but can nonetheless cause significant isotopic shifts. Field applications of Cr, Se, U isotope ratios have demonstrated that they are useful as indicators of reduction in natural systems. A few studies suggest they are also useful as indicators of oxidation and contaminant sources. The physical and chemical complexity of groundwater systems hinders accurate quantitative interpretation of Cr, Se, U isotope data using simple models. Numerical models have been developed that capture the behavior of complex, coupled systems and enable the most effective extraction of information from field data sets.

58 GEOSCIENCES↗

Natural attenuation of uranium in a fluvial Wetland: Importance of hydrology and speciation

A nuclear fuel fabrication facility released 43,500 kg of uranium into a riparian wetland located on the Savannah River Site between 1955 and 1988. Studies were undertaken to evaluate hydrological and geochemical processes influencing uranium accumulation in the wetland. Gamma-radiation-mapping surveys were conducted by systematically walking over the contaminated wetland with backpacks equipped with global positioning systems and NaI gamma detectors. Based on maps compiled from >700,000 gamma spectra and eight sediment uranium depth profiles, it was determined that 94% of the released uranium remained in the wetland. The uranium in the wetland is concentrated in five multi-hectare areas along the stream, accounting for ~11% of the land area adjacent to the stream. While land type (upland or wetland) and topography provided a reasonable first approximation of where much of the uranium was deposited, hydrological watershed modeling revealed that the stream velocity was especially slow through many of the hot spots. Here, using autoradiography combined with SEM/EDX measurements of contaminated sediments, surprisingly few hot particles were detected. Instead, uranium was evenly distributed throughout the sampled sediment, suggesting that dissolved uranium had bound to sediment particles that became suspended and later deposited in low energy (low flow velocity) portions of the stream. EXAFS suggested that U atoms were present as individual ions in disordered complexes within the sediment. Furthermore, linear combination analyses suggested that the predominant component of the U(VI) was adsorbed to sediment minerals (~70%) and a minor component (~30%) was associated with organic matter phases. Furthermore, these studies show that wetlands can be extraordinarily effective at binding and retaining uranium, thereby providing a natural barrier to the transport of uranium out of a watershed. However, significant anthropogenic or climatic changes to wetlands, such as those associated with flooding, forest fires, or land use, may disrupt the complex hydrological and biogeochemical balance necessary to maintain long-term immobilization of uranium.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗