Search NASA⌕ Search

SEARCH · Search NASA

Results for “environmental gradients”

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

Saccharomycotina yeasts defy long-standing macroecological patterns

The Saccharomycotina yeasts (“yeasts” hereafter) are a fungal clade of scientific, economic, and medical significance. Yeasts are highly ecologically diverse, found across a broad range of environments in every biome and continent on earth; however, little is known about what rules govern the macroecology of yeast species and their range limits in the wild. Here, we trained machine learning models on 12,816 terrestrial occurrence records and 96 environmental variables to infer global distribution maps at ~1 km2 resolution for 186 yeast species (~15% of described species from 75% of orders) and to test environmental drivers of yeast biogeography and macroecology. We found that predicted yeast diversity hotspots occur in mixed montane forests in temperate climates. Diversity in vegetation type and topography were some of the greatest predictors of yeast species richness, suggesting that microhabitats and environmental clines are key to yeast diversity. We further found that range limits in yeasts are significantly influenced by carbon niche breadth and range overlap with other yeast species, with carbon specialists and species in high-diversity environments exhibiting reduced geographic ranges. Finally, yeasts contravene many long-standing macroecological principles, including the latitudinal diversity gradient, temperature-dependent species richness, and a positive relationship between latitude and range size (Rapoport’s rule). These results unveil how the environment governs the global diversity and distribution of species in the yeast subphylum. These high-resolution models of yeast species distributions will facilitate the prediction of economically relevant and emerging pathogenic species under current and future climate scenarios.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The effects of advanced spectral line shapes on atmospheric carbon dioxide retrievals

Line shape effects such as speed-dependence and collisional (Dicke) narrowing are manifest in high signal-to-noise ratio laboratory spectra for a wide range of molecular species and sample conditions. However, these effects are usually neglected in spectroscopic models used in remote sensing of atmospheric greenhouse gases. Here, in this study, we examine the effects of line profile choice on laboratory spectra and XCO 2 values retrieved from Total Carbon Column Observing Network (TCCON) spectra. To this end, we utilize limiting cases of the Hartmann-Tran profile (HTP) with shape parameters recently measured by cavity ring-down spectroscopy. Significant differences highlighting the importance of including speed-dependence effects in forward models were observed both in scaling retrievals as well as profile retrievals. However, because of uncertainties associated with the atmospheric temperature gradient, it was not possible to extract accurate carbon column profiles even with inclusion of these new line shape parameters and advanced line profiles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Efficient simulation of optical spectra via machine learning and physical decomposition of environmental effects

Simulations of optical spectra can provide key insights to aid experimental interpretation of electronic excitation phenomena. For chromophores in the condensed phase, these spectra, which incorporate the coupling between electronic excitation and molecular and solvent nuclear motions, can be simulated using excitation energies obtained from molecular dynamics simulations of the chromophore and solvent. Here, we present a hybrid scheme that exploits machine learning and physically informed spectral densities to show that as few as 25 ground and excited state energetic gradient calculations can be used to construct models that accurately predict environment-influenced vibronic coupling in optical spectra. We demonstrate our approach for the green fluorescent protein chromophore in water and the cresyl violet chromophore in methanol. We show that our hybrid approach, employing a machine learning model for the high-frequency spectral density and an ab initio parameterized Debye spectral density for the low-frequency, results in a systematic improvement of the optical absorption lineshape, leading to a simple machine learning scheme that can be used for the simulation of spectral densities and optical spectra.

Snider, Andrew [Univ. of California, Merced, CA (U↗
DOE OSTI2022

Impacts of sulfonic acids on fungal manganese oxide production

Microbial manganese (Mn) oxidation plays a critical role in Mn(III/IV) oxide formation in modern day environmental systems. These oxidation processes and resulting biominerals are sensitive to aqueous conditions, such as pH and dissolved constituent concentrations. With Mn and sulfur (S) biogeochemical cycling closely linked in many environmental systems, and dissolved organic sulfur comprising a substantial pool of total sulfur in several environments, the impact of dissolved organic sulfur compounds on Mn redox processes is important to consider. Sulfonic acids, environmentally ubiquitous organosulfur compounds, play substantial roles in S cycling in many natural and contaminated systems. Research to assess the effects of these abundant sulfonic acids on Mn biogeochemical cycling, microbial Mn oxidation processes, and Mn biominerals is needed for understanding and predicting the impact of coupled S and Mn biogeochemical cycles, particularly in environments with dynamic redox gradients or in anthropogenically contaminated systems. Further, with research on microbial and abiotic Mn oxidation processes often using aminosulfonic acids to control pH, understanding the impact of such sulfonic acids on microbial Mn oxidation processes is critical. Several recent studies found that commonly used zwitterionic N substituted aminosulfonic acids, known as Good’s buffers, such as HEPES and MES, can alter abiotic birnessite sheet structures. Here we investigate the impact of two sulfonic acids with broad applications to natural and contaminated sites as well as laboratory settings (HEPES and MES) on fungal Mn oxidation relative to a carbonate buffer and a buffer-free control by three Ascomycete fungi known to oxidize Mn(II): Stagonospora sp. SRC1lsM3a, Paraphaeosphaeria sporulosa AP3s5–JAC2a, and Plectosphaerella cucumerina DS2psM2a2. Structural analyses of the products show that sulfonic acids promote Mn oxidation by P. cucumerina, producing hexagonally symmetric phyllomanganates analogous to hexagonal birnessite or c–disordered H + birnessite [(Ca,Na,K)(Mn 4+ ,Mn 3+ ,$\square$)O 2 nH 2 O], with solid–associated Mn(II) bound to vacancy sites and biomass, while in their absence almost all Mn remains as either aqueous Mn(II) or solid–associated Mn(II) bound to biomass. In contrast, sulfonic acids exert the opposite effect on Mn oxidation by P. sporulosa, with their presence suppressing Mn(II) oxidation to Mn(IV), likely leading to the formation of mycogenic bixbyite (Mn 3+ 2 O 3 ) while the buffer–free control forms a poorly crystalline phyllo- or tectomanganate. Meanwhile, all treatments exert a minimal effect on Mn(II) uptake from solution and Mn oxidation with Stagonospora sp., with all experimental systems and controls forming poorly crystalline, hexagonally symmetric phyllomanganates. The fact that the sulfonic acids here studied exert similar effects on Mn oxidation, but substantially different effects for each fungus, suggests they affect Mn oxidation via mechanistically similar pathways that are likely dependent on interactions with fungal exudates (which vary from species to species) or specific fungal Mn oxidation processes. Interestingly, for all fungi, MES increases Mn(III) in the resulting biominerals, while the carbonate buffer consistently decreases Mn(III). Further, these results clearly demonstrate that sulfonic acids not only alter Mn oxide structures, as has been previously noted in abiotic studies, but can interfere with Mn oxidation reactions themselves, highlighting the incredible sensitivity of both Mn oxide structures and the Mn oxidation process to the aqueous environment.

58 GEOSCIENCES↗
DOE OSTI2024

Jackson, L., Johnson, M.B., Latrach, A., Grimes, D., Martinez, C., and Mclaughlin, J.F., 2024, Multidisciplinary geotechnical data collection, curation, and analysis for conformity with the regulatory framework for geologic carbon storage in Wyoming, USA: Geological Society of America Abstracts with Programs. Vol. 56, No. 5, 2024, doi: 10.1130/abs/2024AM-405024

Title: Multidisciplinary Geotechnical Data Collection, Curation, and Analysis for Conformity with the Regulatory Framework for Geologic Carbon Storage in Wyoming, USA. Text: Construction and operation of wells for geologic sequestration of carbon dioxide necessitate that they are permitted under the Environmental Protection Agency’s Underground Injection Control Class VI requirements. Class VI wells conform to stringent requirements to ensure long-term safety and integrity of the storage site and the protection of Underground Sources of Drinking Water. Entities pursuing Class VI permitting must provide comprehensive geologic site characterization, including regional geologic structure and stratigraphy, aquifer information, reservoir and confining unit geomechanical properties, geochemical analyses, assessment of trapping capacity and mechanisms, and a variety of other of multidisciplinary geotechnical data. The Wyoming Class VI Site Characterization Database Project is focused on developing a geologic site characterization database of geotechnical information, which has been compiled and verified from established, public databases/entities and scientific literature to expedite Class VI permitting in Sweetwater County within the Greater Green River Basin of southern Wyoming. The preliminary suite of compiled data from 14,000 wells includes 8,000 wells with logs and 7,250 wells with formation tops, ~70 wells with core data (e.g., X-Ray diffraction, petrographic, and petrophysical data), ~2,500 water analyses, ~740 seismic events data, and ~520 bottom-hole temperature measurements. Future work on—and stemming from—this project will include new core analyses, calculation and interpolation of subsurface temperature gradients, mechanical earth models, geochemical simulations, storage capacity estimation, stratigraphic column generation and correlation, and construction of subsurface maps. Finally, this work will help to inspire and facilitate subsurface data compilation and curation beyond Sweetwater County, Wyoming.

42 ENGINEERING↗
DOE OSTI2021

Canopy position influences the degree of light suppression of leaf respiration in abundant tree genera in the Amazon Forest

Leaf respiration in the dark (R dark ) and light (R day ) is poorly characterized in diverse tropical ecosystems, and little to no information exists on the degree of light suppression in common tree species within the Amazon basin, and their dependences upon plant functional traits and position within the canopy. We quantified R dark and apparent R day using the Kok method and measured key leaf traits in 26 tree individuals of different species distributed in three different canopy positions: canopy, lower canopy, and understory. To explore the relationships between the leaf traits we used the standardized major axis (SMA). We found that canopy trees had significantly higher rates of R dark and R day than trees in the understory. The difference between R dark and R day (the light suppression of respiration) was greatest in the understory (68 ± 9%, 95% CI) and lower canopy (49 ± 9%, 95% CI) when compared to the canopy (37 ± 10%, 95% CI). We also found that R day was significantly and strongly correlated with R dark (p < 0.001) for all the canopy positions. Also, leaf mass per area (LMA) and leaf Phosphorus concentration (P) had a significant relationship with R dark (p < 0.001; p = 0.003), respectively. In addition, a significant relationship was found for LMA in the canopy and lower canopy positions (p = 0.009; p = 0.048) while P was only significant in the canopy (p = 0.044). Finally, no significant relationship was found between R dark and nitrogen, sugars, and starch. Our results highlight the importance of including representation of the light suppression of leaf respiration in terrestrial biosphere models and also of accounting for vertical gradients within forest canopies and connections with functional traits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The nitrogen gap in soil health concepts and fertility measurements

Soil nitrogen (N) often limits productivity in agroecosystems, prompting fertilizer applications that increase crop yields but can degrade the environment. Nitrogen's dual role in both productivity and environmental quality should center it in soil health frameworks. We use recent evidence to argue that N availability is an emergent property of the integrated soil biogeochemical system and is strongly influenced by plant traits and their interactions with microbes and minerals. Building upon this, we theorize that the sources of plant and microbial N shift across soil health gradients, from inorganic N dependence in ecologically simple systems with poor soil health to a highly networked supply of organic N in healthy soils; ergo, investments in soil health should increase ecological complexity and the pathways by which plants can access N, leading to more resilient nutrient supplies and yields in a variable climate. However, current N assessment methods derive from a historical emphasis on inorganic N pool sizes and are unable to capture the shifting drivers of N availability across soil health gradients. We highlight the need to better understand the plant-microbial-mineral interactions that regulate bioavailable N as a first step to improving our ability to measure it. We conclude it will be necessary to harness agroecosystem complexity, account for plant and microbial drivers, and carefully integrate external N inputs into soils' internal N network to expand the routes by which N from organic pools can be made bioavailable. Finally, by emphasizing N in soil health concepts, we argue that researchers can accelerate advances in N use efficiency and resiliency.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Effects of forest structural and compositional change on forest microclimates across a gradient of disturbance severity

Forest structural diversity and community composition are key in regulating forest microclimates. When disturbance affects structural diversity or composition, forest microclimates may be altered due to changes in soil temperature, soil water content, and light availability. It is unclear however which structural or compositional components, when changed or to what extent, result in microclimatic change. To address this question, we used data from a large scale, manipulative stem-girdling experiment in northern, lower Michigan—the Forest Resilience and Threshold Experiment (FoRTE). FoRTE follows a factorial design with multiple levels of disturbance severity (0, 45, 65, 85%) based on targeted reductions in gross leaf area index via stem-girdling induced mortality. These disturbance severity treatments are applied in two ways: either as top-down (largest trees are killed) or bottom-up (small to medium trees killed) treatments. We examined how multiple components of structural diversity and community composition changed as a product of disturbance severity and type, and then tested for resulting effects on forest microclimates (light availability, soil temperature, and soil water), using a multivariate, Random Forest framework. We found that measures of community composition (species richness, species evenness, and Shannon-Wiener Diversity Index) and stand structure (basal area, standard deviation of DBH, tree size diversity) declined more following disturbance than did measures of canopy cover, heterogeneity, arrangement, or height. However, when changes in each variable from pre- to post-disturbance, measured as log change, were employed in a multivariate, Random Forest regression framework, structural diversity measures of heterogeneity (rugosity, top rugosity), cover (canopy cover), and arrangement (porosity) were the most influential variables, but with differences among bottom-up and top-down treatments We found that the death of large trees from disturbance impacts soil temperature, water, and light environments more substantially and uniformly across disturbance gradients than does the death of smaller trees. Furthermore, our results have implications for both statistical and process-based modeling of forest disturbance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Iron mineral dissolution releases iron and associated organic carbon during permafrost thaw

It has been shown that reactive soil minerals, specifically iron(III) (oxyhydr)oxides, can trap organic carbon in soils overlying intact permafrost, and may limit carbon mobilization and degradation as it is observed in other environments. However, the use of iron(III)-bearing minerals as terminal electron acceptors in permafrost environments, and thus their stability and capacity to prevent carbon mobilization during permafrost thaw, is poorly understood. We have followed the dynamic interactions between iron and carbon using a space-for-time approach across a thaw gradient in Abisko (Sweden), where wetlands are expanding rapidly due to permafrost thaw. We show through bulk (selective extractions, EXAFS) and nanoscale analysis (correlative SEM and nanoSIMS) that organic carbon is bound to reactive Fe primarily in the transition between organic and mineral horizons in palsa underlain by intact permafrost (41.8 ±10.8 mg carbon per g soil, 9.9 to 14.8% of total soil organic carbon). During permafrost thaw, water-logging and O 2 limitation lead to reducing conditions and an increase in abundance of Fe(III)-reducing bacteria which favor mineral dissolution and drive mobilization of both iron and carbon along the thaw gradient. By providing a terminal electron acceptor, this rusty carbon sink is effectively destroyed along the thaw gradient and cannot prevent carbon release with thaw.

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

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Power Plant Testing of Ultrasonic Measurements of Temperature Distributions and Heat Fluxes to Heat Exchange Surfaces

Extreme environmental conditions are common in energy conversion applications, in which the most hardened insertion sensors often do not perform reliably for long. On the other hand, ultrasonic measurements can be acquired noninvasively, with sensitive components kept away from the damaging environment. We have previously developed the ultrasonic method for measuring the spatial distribution of temperatures in solid materials applicable when large thermal gradients are present. In the developed approach, we use the echogenically segmented ultrasound propagation path, structured to contain engineered or naturally occurring echogenic features, to produce a train of echoes in response to an external pulse of ultrasonic excitation. This paper outlines the testing results of this approach applied to the temperature distribution measurements to a waterwall of a 500 MW power-plant utility boiler. The validation results obtained over the prolonged power plant test show that the estimated temperature profile is correctly captured, and the measurement accuracy can be comparable with traditional insertion sensors, such as thermocouples. Overall, the testing has confirmed that the developed approach has matured to become an attractive alternative to conventional sensing in solving challenging problems of long-term temperature measurements in extreme environments. Heat fluxes and thermal stresses in the structure can then be characterized noninvasively using the measured temperature distribution as the basis.

Walton, Kenneth↗
DOE OSTI2021

The genetic basis of the root economics spectrum in a perennial grass

Significance Plants have evolved diverse root form and function across the Earth’s environment, yet only certain types of trait combinations have proved evolutionarily viable. Using a genetic mapping population of the native perennial switchgrass, our study demonstrates multiple genetic linkages among root morphology, growth, and turnover. Switching alleles derived from southern-origin ecotypes to alleles from northern-origin ecotypes increases root turnover but reduces tissue investment in root length construction. The genetic trade-off between construction and turnover likely facilitates the local adaptation of root strategy along the warm to cold climatic gradients of the species range. In practice, our study provides the genetic evidence that increasing switchgrass yield for bioenergy does not directly conflict with enhancing its root-derived carbon sequestration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Explaining Globally Inhomogeneous Future Changes in Monsoons Using Simple Moist Energy Diagnostics

Here, this study examines the annual cycle of monsoon precipitation simulated by models from phase 6 of the Coupled Model Intercomparison Project (CMIP6), then uses moist energy diagnostics to explain globally inhomogeneous projected future changes. Rainy season characteristics are quantified using a consistent method across the globe. Model bias is shown to include rainy season onsets tens of days later than observed in some monsoon regions (India, Australia, and North America) and overly large summer precipitation in others (North America, South America, and southern Africa). Projected next-century changes include rainy season lengthening in the two largest Northern Hemisphere monsoon regions (South Asia and central Sahel) and shortening in the two largest Southern Hemisphere regions (South America and southern Africa). Changes in the North American and Australian monsoons are less coherent across models. To understand these changes, relative moist static energy (MSE) is defined as the difference between local and tropical-mean surface air MSE. Future changes in relative MSE in each region correlate well with onset and demise date changes. Furthermore, Southern Hemisphere regions projected to undergo rainy season shortening are spanned by an increasing equator-to-pole MSE gradient, suggesting their rainfall will be increasingly inhibited by fluxes of dry extratropical air; Northern Hemisphere regions with projected lengthening of rainy seasons undergo little change in equator-to-pole MSE gradient. Thus, although model biases raise questions as to the reliability of some projections, these results suggest that globally inhomogeneous future changes in monsoon timing may be understood through simple measures of surface air MSE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-density CRISPRi screens reveal diverse routes to improved acclimation in cyanobacteria

Cyanobacteria are the oldest form of photosynthetic life on Earth and contribute to primary production in nearly every habitat, from permafrost to hot springs. Despite longstanding interest in the acclimation of these microbes, it remains poorly understood and challenging to rewire. Here, this study uses a high-density, genome-wide CRISPR interference screen to examine the influence of gene-specific transcriptional variation on the growth of Synechococcus sp. PCC 7002 under environmental extremes. Surprisingly, many partial knockdowns enhanced fitness under cold monochromatic conditions. Transcriptional repression of genes for core subunits of the NDH-1 complex, which are important for photosynthesis and carbon uptake, improved growth rates under both red and blue light but at distinct, color-specific optima. Most genes with fitness-improving knockdowns were distinct to each light color, and dual-target transcriptional repression produced nonadditive effects. Findings reveal diverse routes to improved acclimation in cyanobacteria (e.g., attenuation of genes involved in CO 2 uptake, light harvesting, translation, and purine metabolism) and provide an approach for using gradients in sgRNA activity to pinpoint biochemically influential transcriptional changes in cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Climate-Driven Legacies in Simulated Microbial Communities Alter Litter Decomposition Rates

The mechanisms underlying diversity-functioning relationships have been a consistent area of inquiry in biogeochemistry since the 1950s. Though these mechanisms remain unresolved in soil microbiomes, many approaches at varying scales have pointed to the same notion—composition matters. Confronting the methodological challenge arising from the complexity of microbiomes, this study used the model DEMENTpy, a trait-based modeling framework, to explore trait-based drivers of microbiome-dependent litter decomposition. We parameterized DEMENTpy for five sites along a climate gradient in Southern California, United States, and conducted reciprocal transplant simulations analogous to a prior empirical study. The simulations demonstrated climate-dependent legacy effects of microbial communities on plant litter decomposition across the gradient. This result is consistent with the previous empirical study across the same gradient. An analysis of community-level traits further suggests that a 3-way tradeoff among resource acquisition, stress tolerance, and yield strategies influences community assembly. Simulated litter decomposition was predictable with two community traits (indicative of two of the three strategies) plus local environment, regardless of the system state (transient vs. equilibrium). Although more empirical confirmation is still needed, community traits plus local environmental factors (e.g., environment and litter chemistry) may robustly predict litter decomposition across spatial-temporal scales. In conclusion, this study offers a potential trait-based explanation for climate-dependent community effects on litter decomposition with implications for improved understanding of whole-ecosystem functioning across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗