Search NASA⌕ Search

SEARCH · Search NASA

Results for “total correlation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Reproductive and leaf litterfall fluxes in forest ecosystem sites globally (1950-2022)

Forest allocation of net primary productivity (NPP) to reproduction is poorly quantified globally, despite its critical role in forest regeneration and a well-supported trade-off with allocation to growth. Although field measurements of total NPP are rare, our work finds that a proxy for reproductive carbon allocation constructed from leaf (L) and reproductive (R) litterfall fluxes, R/(R+L), is strongly correlated with R/NPP, facilitating analysis across a wide range of sites where biometric estimates of NPP are not available (R² = 0.85; Hanbury-Brown et al., 2022, Ward et al., in prep). To investigate relationships between ecosystem-scale reproductive allocation (RA) and climate, soil fertility, and stand age gradients, we conducted a literature search and synthesized 824 observations of annual average leaf and reproductive litterfall fluxes across forest sites globally. The zip file includes 1) a folder Data/ containing the litterfall data ("GlobalForestRA_data.csv") and metadata ("GlobalForestRA_metadata.doc") files. The data file includes geographic coordinates, long-term mean annual temperature and precipitation (1970-2000, extracted from WorldClim2.1), leaf and reproductive litterfall fluxes, sampling interval and protocols, forest characteristics (dominant leaf morphology, information pertaining to forest age and successional stage, and disturbance history) and soil properties (% sand, %silt, %clay, total phosphorus (P), nitrogen (N), cation exchange capacity (CEC) and pH) extracted from SoilGrids250 and from on-site measurements, where available. The metadata file contains information about each variable reported in the data file, including data sources, processing methods, and all references. The Data folder contains two additional files used to create Figure 1; these are described in greater detail in the README.2) R scripts GloalForestRA_analysis.r and GlobalForestRA_SI.r and a folder /Functions used to produce results, figures, and tables in the manuscript Ward et al. (in press)3) a README file describing how the data and R scripts can be used to reproduce statistical results, figures, and tables found in the manuscript. Ward et al. (in press)This repository can also be found at: https://github.com/r-ward/Global_Analysis_ForestRA.Ward, R.E., Zhang-Zheng, H. Aernethy, K., Adu-Bredu, S., Arroyo, L., Bailey, A. et al. (in press). Forest age rivals climate to explain reproductive allocation patterns in forest ecosystems globally. Ecology Letters. Hanbury-Brown, A.R., Ward, R.E. & Kueppers, L.M. (2022). Forest regeneration within Earth system models: current process representations and ways forward. New Phytol., 235, 20–40.Ward et al. (2025), Forest age rivals climate to explain reproductive allocation patterns in forest ecosystems globally, in prep.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1980

Collisional redistribution of radiation. II - The effects of degeneracy on the equations of motion for the density matrix. III - The equation of motion for the correlation function and the scattered spectrum

The effect of correlations between an absorber atom and perturbers in the binary-collision approximation are applied to degenerate atomic systems. A generalized absorption profile which specifies the final state of the atom after an absorption event is related to the total intensities of Rayleigh scattering and fluorescence from the atom. It is suggested that additional dynamical information to that obtainable from ordinary absorption experiments is required in order to describe redistributed atomic radiation. The scattering of monochromatic radiation by a degenerate atom is computed in a binary-collision approximation; an equation of motion is derived for the correlation function which is valid outside the quantum-regression regime. Solutions are given for the weak-field conditions in terms of generalized absorption and emission profiles that depend on the indices of the atomic multipoles.

Burnett, K.↗
DOE OSTI2015

THE LOCAL GROUP IN THE COSMIC WEB

We explore the characteristics of the cosmic web around Local-Group (LG)-like pairs using a cosmological simulation in the ΛCDM cosmology. We use the Hessian of the gravitational potential to classify regions on scales of ∼2 Mpc as a peak, sheet, filament, or void. The sample of LG counterparts is represented by two samples of halo pairs. The first is a general sample composed of pairs with similar masses and isolation criteria as observed for the LG. The second is a subset with additional observed kinematic constraints such as relative pair velocity and separation. We find that the pairs in the LG sample with all constraints are: (1) preferentially located in filaments and sheets, (2) located in a narrow range of local overdensity 0 < δ < 2, web ellipticity 0.1 < e < 1.0, and prolateness –0.4 < p < 0.4, (3) strongly aligned with the cosmic web. The alignments are such that the pair orbital angular momentum tends to be perpendicular to the smallest tidal eigenvector, e-hat {sub 3}, which lies along the filament direction or the sheet plane. A stronger alignment is present for the vector linking the two halos with the vector e-hat {sub 3}. Additionally, we fail to find a strong correlation between the spin of each halo in the pair with the cosmic web. All of these trends are expected to a great extent from the selection of LG total mass in the general sample. Applied to the observed LG, there is a potential conflict between the alignments of the different satellite planes and the numerical evidence for satellite accretion along filaments; the direction defined by e-hat {sub 3}. This highlights the relevance of achieving a precise characterization for the location of the LG in the cosmic web in the cosmological context provided by ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Novel Method for Preoperative Positioning of Total Ankle Replacement Using 3D Digital Model

Abstract Objective To establish a digital model of the ankle joint through 3D imaging technology and explore the preoperative placement of ankle replacement prostheses. Methods Computed tomography images of intact ankle joints from 54 cases in the outpatient and inpatient departments of our hospital were collected; according to the INBONE® total ankle system surgery process, the surgery model and surgical osteotomy were finished using MIMICS based on computer simulation method. The shortest distance was measured between the center point and the anterior, posterior, medial, and lateral, respectively, to ensure the precise position of the ankle replacement prosthesis by digital simulation surgery. The relationship between the two variables was analyzed by bivariate correlation analysis. Results The dataset of this study included 48 cases of the sub‐data set (26 males and 22 females) and included 27 cases of left ankle and 21 cases of right ankle. The average medial malleolar angle was 18.67°± 2.87°, the average amount of bone resection was 12.13 ± 1.86 cm 3 , the mid‐anterior distance was 1.72 ± 0.19 cm, the mid‐posterior distance was 2.00 ± 0.19 cm, the ratio of mid‐anterior to mid‐posterior was 0.87, the mid‐medial distance was 1.26 ± 0.17 cm, the mid‐lateral distance was 1.19 ± 0.16 cm, and the ratio of mid‐medial to mid‐lateral was 1.06. After osteotomy, the anteroposterior diameter was 3.73 ± 0.32 cm, the transverse diameter was 2.46 ± 0.27 cm, and the ratio of anteroposterior diameter to transverse diameter was 1.53. In the bottom view, the shape after osteotomy is rectangular. The mid‐anterior distance was strongly negatively correlated with age, the mid‐anterior distance and the amount of bone resection, the mid‐medial distance and the amount of bone resection, the mid‐lateral distance and the amount of bone resection, the mid‐lateral distance and the anteroposterior diameter, the anteroposterior diameter and the transverse diameter were all strongly positively correlated. Conclusion The projection point of the lower tibia centerline on the tibial horizontal osteotomy surface is located at a position slightly anterior to the midpoint of the transverse diameter after ankle arthroplasty. The rational positioning of the total ankle replacement is located at both a position slightly anterior to the midpoint of the transverse diameter and midpoint of the anteroposterior diameter, which can be used as a reference method before total ankle arthroplasty surgery.

Wu, Shi‐xun↗
DOE OSTI2020

Constraining level densities through quantitative correlations with cross-section data

The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the accurate calculation of cross sections for the different reaction channels. Many common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account in these models: LD at neutron separation energy coming from average spacings of s- and p-wave resonances ( D 0 and D 1 , respectively) whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only in very specific regions of excitation energy, and for specific spins and parities. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for fitting and the determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD as a starting point as the HFB LD provide a more realistic spin and parity distributions than phenomenological models such as Gilbert-Cameron (GC). We then associate variations predicted by the HFB model with the structure observed in double-differential cross sections at low outgoing neutron energy, a region that is dominated by the LD input. We also use (n, p) on 56 Fe , as an example case where angle-integrated cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of D 0 or D 1 are unknown. Finally, it also predicts neutron-induced inelastic γ cross sections that in some cases can differ significantly from more standard phenomenological LD models such as GC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2022

Decadal and Annual Variations in Meteoric Flux from Ulysses, Wind, and SOFIE Observations

Our solar system is filled with meteoric particles, or cosmic dust, which is either interplanetary or interstellar in origin. Interstellar dust (ISD) enters the heliosphere due to the relative motion of the sun and the interstellar flow. Interplanetary dust (IPD) comes primarily from asteroid collisions or comet sublimation, and comprises the bulk of material entering Earth’s atmosphere. This study examines variations in ISD and the IPD flux at Earth using observations from three different satellite techniques. First are size-resolved in situ meteoroid detections by the Ulysses spacecraft, and second are in situ indirect dust observations by Wind. Third are measurements of meteoric smoke in the mesosphere by the Solar Occultation For Ice Experiment (SOFIE). Wind and Ulysses observations are sorted into the interstellar and interplanetary components. Wind ISD show the anticipated correlation to the 22-yr. solar magnetic cycle, and are consistent with model predictions of ISD. Because Wind does not discriminate particle size, the IPD measurements were interpreted using meteoric mass distributions from Ulysses observations and from different models. Wind observations during 2007-2020 indicate a total meteoric influx at Earth of 22 metric tons per day (t d(exp -1)), in reasonable agreement with long-term averages from SOFIE (25 t d(exp -1)) and Ulysses (32 t d(exp -1)). The SOFIE and Wind influx time series both show an unexpected correlation to the 22-yr. solar cycle. This relationship could be an artifact, or may indicate that IPD responds to changes in the solar magnetic field.

Meteoric influx↗
NASA NTRS1994

Depletions in winter total ozone values over southern England

A study has been made of the recently re-evaluated time series of daily total ozone values for the period 1979 to 1992 for southern England. The series consists of measurements made at two stations, Bracknell and Camborne. The series shows a steady decline in ozone values in the spring months over the period, and this is consistent with data from an earlier decade that has been published but not re-evaluated. Of exceptional note is the monthly mean for January 1992 which was very significantly reduced from the normal value, and was the lowest so far measured for this month. This winter was also noteworthy for a prolonged period during which a blocking anticyclone dominated the region, and the possibility existed that this was related to the ozone anomaly. It was possible to determine whether the origin of the low ozone value lay in ascending stratospheric motions. A linear regression analysis of ozone value deviation against 100hPa temperature deviations was used to reduce ozone values to those expected in the absence of high pressure. The assumption was made that the normal regression relation was not affected by atmospheric anomalies during the winter. This showed that vertical motions in the stratosphere only accounted for part of the ozone anomaly and that the main cause of the ozone deficit lay either in a reduced stratospheric circulation to which the anticyclone may be related or in chemical effects in the reduced stratospheric temperatures above the high pressure area. A study of the ozone time series adjusted to remove variations correlated with meteorological quantities, showed that during the period since 1979, one other winter, that of 1982/3, showed a similar although less well defined deficit in total ozone values.

Lapworth, A.↗
NASA NTRS1984

Aggregating available soil water holding capacity data for crop yield models

The total amount of water available to plants that is held against gravity in a soil is usually estimated as the amount present at -0.03 MPa average water potential minus the amount present at -1.5 MPa water potential. This value, designated available water-holding capacity (AWHC), is a very important soil characteristic that is strongly and positively correlated to the inherent productivity of soils. In various applications, including assessing soil moisture status over large areas, it is necessary to group soil types or series as to their productivity. Current methods to classify AWHC of soils consider only total capacity of soil profiles and thus may group together soils which differ greatly in AWHC as a function of depth in the profile. A general approach for evaluating quantitatively the multidimensional nature of AWHC in soils is described. Data for 902 soil profiles, representing 184 soil series, in Indiana were obtained from the Soil Characterization Laboratory at Purdue University. The AWHC for each of ten 150-mm layers in each soil was established, based on soil texture and parent material. A multivariate clustering procedure was used to classify each soil profile into one of 4, 8, or 12 classes based upon ten-dimensional AWHC values. The optimum number of classes depends on the range of AWHC in the population of oil profiles analyzed and on the sensitivity of a crop to differences in distribution of water within the soil profile.

Seubert, C. E.↗
DOE OSTI2016

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory

Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

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

Reduced quantum dynamics with arbitrary bath spectral densities: Hierarchical equations of motion based on several different bath decomposition schemes

We investigated applications of the hierarchical equation of motion (HEOM) method to perform high order perturbation calculations of reduced quantum dynamics for a harmonic bath with arbitrary spectral densities. Three different schemes are used to decompose the bath spectral density into analytical forms that are suitable to the HEOM treatment: (1) The multiple Lorentzian mode model that can be obtained by numerically fitting the model spectral density. (2) The combined Debye and oscillatory Debye modes model that can be constructed by fitting the corresponding classical bath correlation function. (3) A new method that uses undamped harmonic oscillator modes explicitly in the HEOM formalism. Methods to extract system-bath correlations were investigated for the above bath decomposition schemes. We also show that HEOM in the undamped harmonic oscillator modes can give detailed information on the partial Wigner transform of the total density operator. Theoretical analysis and numerical simulations of the spin-Boson dynamics and the absorption line shape of molecular dimers show that the HEOM formalism for high order perturbations can serve as an important tool in studying the quantum dissipative dynamics in the intermediate coupling regime.

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

Exact isovector pairing in a shell-model framework: Role of proton-neutron correlations in isobaric analog states

Here, we utilize a nuclear shell model Hamiltonian with only two adjustable parameters to generate, for the first time, exact solutions for pairing correlations for light to medium-mass nuclei, including the challenging proton-neutron pairs, while also identifying the primary physics involved. In addition to single-particle energy and Coulomb potential terms, the shell model Hamiltonian consists of an isovector T = 1 pairing interaction and an average proton-neutron isoscalar T = 0 interaction, where the T = 0 term describes the average interaction between nonpaired protons and neutrons. This Hamiltonian is exactly solvable, where, utilizing three to seven single-particle energy levels, we reproduce experimental data for 0 + state energies for isotopes with mass A = 10 through A = 62 exceptionally well including isotopes from He to Ge. Additionally, we isolate effects due to like-particle and proton-neutron pairing, provide estimates for the total and proton-neutron pairing gaps, and reproduce N (neutron) = Z (proton) irregularity. These results provide a further understanding for the key role of proton-neutron pairing correlations in nuclei, which is especially important for waiting-point nuclei on the $\mathcal{rp}$ path of nucleosynthesis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Modelling the sulfate capacity of simulated radioactive waste borosilicate glasses

The capacity of simulated high-level radioactive waste borosilicate glasses to incorporate sulfate has been studied as a function of glass composition. Combined Raman, 57 Fe Mössbauer and literature evidence supports the attribution of coordination numbers and oxidation states of constituent cations for the purposes of modelling, and results confirm the validity of correlating sulfate incorporation in multicomponent borosilicate radioactive waste glasses with different models. A strong compositional dependency is observed and this can be described by an inverse linear relationship between incorporated sulfate (mol% SO 4 2- ) and total cation field strength index of the glass, Σ(z/a 2 ), with a high goodness-of-fit (R 2 ≈ 0.950). Similar relationships are also obtained if theoretical optical basicity, Λ th (R 2 ≈ 0.930) or non-bridging oxygen per tetrahedron ratio, NBO/T (R 2 ≈ 0.919), are used. Results support the application of these models, and in particular Σ(z/a 2 ), as predictive tools to aid the development of new glass compositions with enhanced sulfate capacities.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Perspective: Toward “synthesis by design”: Exploring atomic correlations during inorganic materials synthesis

Synthesis of inorganic extended solids is a critical starting point from which real-world functional materials and their consequent technologies originate. However, unlike the rich mechanistic foundation of organic synthesis, with its underlying rules of assembly (e.g., functional groups and their reactivities), the synthesis of inorganic materials lacks an underpinning of such robust organizing principles. In the latter case, any such rules must account for the diversity of chemical species and bonding motifs inherent to inorganic materials and the potential impact of mass transport on kinetics, among other considerations. Without such assembly rules, there is less understanding, less predictive power, and ultimately less control of properties. Despite such hurdles, developing a mechanistic understanding for synthesis of inorganic extended solids would dramatically impact the range of new material discoveries and resulting new functionalities, warranting a broad call to explore what is possible. Here we discuss our recent approaches toward a mechanistic framework for the synthesis of bulk inorganic extended solids, in which either embryonic atomic correlations or fully developed phases in solutions or melts can be identified and tracked during product selection and crystallization. The approach hinges on the application of high-energy x-rays, with their penetrating power and large Q-range, to explore reaction pathways in situ. We illustrate this process using two examples: directed assembly of Zr clusters in aqueous solution and total phase awareness during crystallization from K–Cu–S melts. These examples provide a glimpse of what we see as a larger vision, in which large scale simulations, data-driven science, and in situ studies of atomic correlations combine to accelerate materials discovery and synthesis, based on the assembly of well-defined, prenucleated atomic correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effect of Pressure on the Electronic, Magnetic, and Optical Properties of the In{sub 0.75}Cr{sub 0.25}P Compound

Using the density functional theory (DFT) approach and full-potential linearized augmented plane waves (FP-LAPW) method, some of the electronic (density of states, band structure, and Fermi surfaces), magnetic (total and local magnetic moments, spin polarization, exchange constants), and optical properties (dielectric function, energy loss function, and refractive index) were investigated. The GGA + U approximation was used for the exchange-correlation term. The effect of hydrostatic pressure up to 20 GP has been studied on these properties. The results show that the In{sub 0.75}Cr{sub 0.25}P compound is a half metal with 100% spin polarization at the Fermi level and the increment of pressure does not change its half-metallic behavior. The half-metal gap achieved for this compound is variable between 1.38 and 1.5 eV. The interaction between Cr and P atoms in this compound is an anti-ferromagnetic type. It was observed that the total magnetic moment for this compound is about 3 μ{sub B}. Local magnetic moments change with increasing pressure, but the total magnetic moment remains unaffected by increasing the pressure. Calculations of the optical properties of the In{sub 0.75}Cr{sub 0.25}P compound showed that the dominant mechanism for the low-energy photons is due to the induced free carriers of the Cr atom, and the optical properties of this compound, in particular the refractive index in the visible region, change with increasing pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Spectroscopy of charmed baryons from lattice QCD

We present the ground and excited state spectra of singly, doubly and triply charmed baryons by using dynamical lattice QCD. A large set of baryonic operators that respect the symmetries of the lattice and are obtained after subduction from their continuum analogues are utilized. Using novel computational techniques correlation functions of these operators are generated and the variational method is exploited to extract excited states. The lattice spectra that we obtain have baryonic states with well-defined total spins up to 7/2 and the low lying states remarkably resemble the expectations of quantum numbers from SU(6) x O(3) symmetry. Various energy splittings between the extracted states, including splittings due to hyperfine as well as spin-orbit coupling, are considered and those are also compared against similar energy splittings at other quark masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1988

Correlations between optical/near-infrared and UV extinction parameters and the prediction of UV extinction from ground-based photometry

An analysis of the variability of absolute UV extinction, A lambda derived through the ratio total-to-selective extinction, R Av/E(B-V), for 31 lines of sight for which UV extinction parameters were derived, is presented. For these data, which sample a wide range of environments and which are characterized by R = 2.5 to 6.0, differences in the shapes of UV extinction curves are largely due to variations in optical/near-UV extinction through changes in R. From this, it is found that UV extinction curves can be reproduced with reasonable accuracy by simply knowing R. Detection of an A2175/E(B-V) = 10 + or - 1 mag implies that the total extinction at 2175 A can be estimated from E(B-V) alone.

Cardelli, Jason A.↗
DOE OSTI2021

Consistent Controls on Trace Metal Micronutrient Speciation in Wetland Soils and Stream Sediments

Trace metal are essential for microbially-mediated biogeochemical processes occurring in anoxic wetland soils and stream bed sediments, but low availability of these elements may inhibit anaerobic element cycling and transformations. Solid-phase speciation is likely a critical control on trace metal availability but has seen limited study in anoxic systems having concentrations similar to geological background levels, where metal limitations may be most prevalent. We have investigated trace metal concentrations and solid-phase speciation in three freshwater subsurface aquatic systems: marsh wetland soils, riparian wetland soils, and the sediments of a streambed. These systems displayed low solid-phase trace metal concentrations, generally at or below geological background levels, which generally followed the trend Zn > Cu ≈ Ni > Co and showed no correlation with major element compositions. All soils and sediments were dominated by quartz but varied in clay mineralogy as well as the organic matter, total sulfur, and total iron contents. X-ray absorption near-edge structure (XANES) spectroscopy shows that sulfur speciation in both wetlands is dominated by organic sulfur. Elemental sulfur and iron sulfides together made up <25% of the sulfur in the wetland soils, but the distribution between inorganic and organic forms was reversed in the stream sediments. Ferrous and ferric iron in clay minerals were common species identified by both XANES and extended X-ray absorption fine structure (EXAFS) spectroscopies at all sites. Iron(III) oxides were substantial components in all but the marsh wetland soils. Quantitative analysis of copper, nickel, and zinc XANES spectra revealed similar metal speciation across all sites. Copper speciation was dominated by sulfides, adsorbed species, and minor amounts of copper bound to organic matter; no metallic copper was detected. Nickel speciation also varied little and was dominated by nickel in clay mineral octahedral sheets and nickel sulfide, with adsorbed species also present. Zinc speciation was slightly more varied, with the marsh wetland soils and stream bed sediments containing adsorbed species, zinc associated with clay mineral structures, and zinc bound to reduced sulfur groups on organic matter, whereas the riparian wetland soils lacked clay-associated zinc but contained zinc sulfide. Trace metals bound to reduced sulfur occurred at every site, with a greater sulfur-bound fraction for copper. The fractional abundance of sulfur-bound species showed no relationship with soil or sediment total sulfur content, which varied by two orders of magnitude. More broadly, the observations in this study suggest that trace metal speciation in freshwater wetland soils and stream sediments is consistently dominated by a small set of recurring components which are distinct for each metal. Furthermore, this may represent a general geochemical phenomenon in anoxic soils and sediments containing trace metals at background concentrations (as low as 3 µg g -1 ) that was not predicted from systems that are contaminated with or naturally-enriched in copper, nickel, or zinc.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Brown watersnakes (Nerodia taxispilota) as bioindicators of mercury contamination in a riverine system

Herein, mercury (Hg) is a contaminant that enters the environment through natural or anthropogenic means. Ecological risk assessments have examined Hg bioaccumulation and effects in many taxa, but little is known about Hg dynamics in reptiles, or their potential use as bioindicator species for monitoring Hg in aquatic systems. Numerous snake species, like North American watersnakes (Nerodia spp.), are piscivorous and are exposed to Hg through their diet. The purpose of this study was to identify factors associated with Hg accumulation in a common watersnake species and compare Hg concentrations of the snakes to those in fish occupying the same habitats. To this end, we sampled brown watersnakes (Nerodia taxispilota) from the Savannah River, a major river system in the southeastern U.S., and compared N. taxispilota Hg accumulation trends to those of bass (Micropterus salmoides), catfish (Ictalurus and Ameiurus spp.), and panfish (Lepomis and Pomoxis spp.) collected from the same reach. Total Hg (THg) in N. taxispilota tail tips ranged from 0.020 to 0.431 mg/kg (wet weight; mean: 0.104 ± 0.008). Snake tail THg was significantly correlated with blood THg, which ranged from 0.003 to 1.140 mg/kg (0.154 ± 0.019). Snake size and site of capture were significantly associated with tail THg. Snake tail THg increased at sites along and downstream of the area of historic Hg pollution, consistent with fish THg. Snake muscle THg was predicted based on tail THg and ranged from 0.095 to 1.160 (0.352 ± 0.022). To gauge Hg biomagnification in N. taxispilota, we compared predicted snake muscle THg concentrations to THg in fish of consumable size. Average biomagnification factors for THg in N. taxispilota were 3.1 (panfish) and 5.4 (catfish), demonstrating N. taxispilota likely biomagnify Hg through their diet. These results reveal N. taxispilota to be an effective bioindicator species for monitoring Hg in aquatic environments.

54 ENVIRONMENTAL SCIENCES↗