Search NASA⌕ Search

SEARCH · Search NASA

Results for “Proxy Signature”

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.

33 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Machine learning and shallow groundwater chemistry to identify geothermal prospects in the Great Basin, USA

This study discovers various geothermal prospects in the Great Basin, USA based on shallow groundwater chemical (geochemical) data. The geochemical data are expected to include hidden (latent) information that is a proxy for geothermal prospectivity. We processed the sparse geochemical data in the Great Basin at 14,341 locations including 18 attributes. Next, a non-negative matrix factorization with customized k-means clustering is applied to the geochemical data matrix that automatically finds three hidden geothermal signatures representing modestly, moderately, and highly confident geothermal prospects. The algorithm also evaluated the probability of occurrence of these types of resources through the studied region. There is a consistency between regional geothermal prospectivity as estimated by our ML methodology and the traditional play fairway analysis conducted over a portion of the study area. We also identify the dominant data attributes associated with each signature. Finally, our ML analyses allow us to reconstruct attributes from sparse into continuous over the study domain. The predicted continuous attributes can be used for future detailed geothermal explorations in the Great Basin.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Differential Measurements of Jet Substructure and Partonic Energy Loss in Au+Au Collisions at √ sNN = 200 GeV

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in p+p and Au+Au collisions at √sNN =200GeV. The substructure observables include SoftDrop groomed momentum fraction (zg), groomed jet radius (Rg), and subjet momentum fraction (zSJ) and opening angle (θSJ). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for p+p collisions and are compared to leading-order Monte Carlo models. The subjet θSJ distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in Au+Au collisions. We compare data from Au+Au collisions to those from p+p which are embedded in minimum-bias Au+Au events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than zg and Rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuumlike. We also report measurements of the differential dijet momentum imbalance (AJ) for jets of varying θSJ. We find no qualitative differences in energy loss signatures for varying angular scales in the range 0.1< θSJ<0.3, leading to the possible interpretation that energy loss in this population of high-momentum dijet pairs, is due to soft medium-induced gluon radiation from a single color charge as it traverses the medium.

Abdallah, M. S.↗
DOE OSTI2025

Behaviour of cadmium isotopes in sulfidic waters and sediments of the Black Sea: Implications for global cadmium cycling and the application of cadmium isotopes as a paleo-oceanographic proxy

The cadmium isotope system has found use as a tracer for biological productivity, redox and organic carbon burial in the oceans. There are, however, very few observational constraints on Cd isotopic behaviour in modern sulfidic marine conditions, limiting our understanding of the modern Cd cycle, and our ability to use Cd isotope measurements of ancient sedimentary deposits as a paleoceanographic proxy. Here we study the behaviour of dissolved Cd in the water column of the Black Sea and its incorporation into sediments on the basin floor. The isotopic composition of dissolved Cd in the upper ∼50 m of the Black Sea water column is controlled by a combination of biological uptake and regeneration along with mixing of river water and Mediterranean seawater. Cadmium declines to <2 % of its peak subsurface concentration within the zone of nitrate-reduction but above the sulfide chemocline. The isotopic composition of dissolved Cd below the sulfide chemocline evolves in a manner that is consistent with experimental Cd sulfide fractionation factors, providing field-based support for these earlier studies. In contrast, the stability of the dissolved Cd isotopic signature above the chemocline, despite a large reduction in dissolved Cd concentrations, is explained by diffusion towards the Cd-sulfide sink below the chemocline. Sediments accumulating in the deep Black Sea below the sulfide chemocline have isotopic compositions that are similar to dissolved Cd in its near-surface aqueous concentration maximum. Sediments accumulating at shallower depths have lighter isotopic compositions that are likely affected by non-sulfidic burial pathways mediated by macrofaunal mixing and diagenesis of shallow burial phases. Our study shows that the isotopic composition of bulk sediments accumulating under sulfidic conditions is similar to dissolved Cd in the upper water column, supporting its use as a paleo-chemical proxy for ancient seawater as long as basin-scale controls on open ocean seawater sources are considered. We show that Cd burial in the deep Black Sea, predominantly as CdS, accounts for ∼0.5–15 % of total annual Cd burial in the global ocean.

Biogeochemical cycles↗
DOE OSTI2021

Dissipation measures in weakly collisional plasmas

ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Declining tropical cyclone frequency under global warming

Assessing the role of anthropogenic warming from temporally inhomogeneous historical data in the presence of large natural variability is difficult and has caused conflicting conclusions on detection and attribution of tropical cyclone (TC) trends. Here, using a reconstructed long-term proxy of annual TC numbers together with high-resolution climate model experiments, we show robust declining trends in the annual number of TCs at global and regional scales during the twentieth century. The Twentieth Century Reanalysis (20CR) dataset is used for reconstruction because, compared with other reanalyses, it assimilates only sea-level pressure fields rather than utilize all available observations in the troposphere, making it less sensitive to temporal inhomogeneities in the observations. It can also capture TC signatures from the pre-satellite era reasonably well. The declining trends found are consistent with the twentieth century weakening of the Hadley and Walker circulations, which make conditions for TC formation less favourable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Characterization of Precipitation-Induced Radon Progeny Deposition Events Using a City-Scale Sensor Network

Networks of radiation detectors provide a platform for real-time radioactive source detection and identification in urban environments. Detection algorithms in these systems must adapt to naturally-occurring changes in background, which requires well-characterized relationships between precipitation events and their corresponding radiological signature. Here, we present a quantitative and qualitative description of rain-induced radon progeny deposition events occurring in Chicago from September 2023 to February 2024. We measure ambient gamma radiation levels, precipitation rate, temperature, pressure, and relative humidity in a network of sensor nodes. For each identified precipitation period, we decompose spectra into static- and radon-associated components as defined by a non-negative matrix factorization (NMF) algorithm. We find a consistent power-law relationship between a precipitation-dependent peak of the radon progeny proxy (RPP) and the peak strength of the radon-associated NMF component for most precipitation events. We conduct a case study of a rainfall period with abnormally high levels of implied radon progeny concentration and describe its temporal and spatial evolution. We hypothesize that this phenomenon is due to the air mass path that intersects a uranium-rich region of Wyoming. Finally, we cluster precipitation events into three distinct categories. One category roughly corresponds to events with deep low-pressure systems and high relative radon concentration, while another is characteristic of light stratiform rain with slightly higher temperatures and intermediate relative radon concentration. The third category appears to contain weak-gradient or lake breeze convection showers with intermittent precipitation and low relative radon concentration. These findings suggest that radiological anomaly detection could be improved by training unique background models corresponding to each category of meteorological event.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A pair of early- and late-forming galaxy cluster samples: A novel way of studying halo assembly bias assisted by a constrained simulation

The halo assembly bias, a phenomenon referring to dependencies of the large-scale bias of a dark matter halo other than its mass, is a fundamental property of the standard cosmological model. First discovered in 2005 from the Millennium Run simulation, it has been proven very difficult to be detected observationally, with only a few convincing claims of detection so far. The main obstacle lies in finding an accurate proxy of the halo formation time. In this study, by utilizing a constrained simulation that can faithfully reproduce the observed structures larger than 2 Mpc in the local universe, for a sample of 634 massive clusters at z ≤ 0.12, we found their counterpart halos in the simulation and used the mass growth history of the matched halos to estimate the formation time of the observed clusters. This allowed us to construct a pair of early- and late-forming clusters, with a similar mass as measured via weak gravitational lensing, and large-scale biases differing at the ≈3σ level, suggestive of the signature of assembly bias, which is further corroborated by the properties of cluster galaxies, including the brightest cluster galaxy and the spatial distribution and number of member galaxies. Our study paves a way to further detect assembly bias based on cluster samples constructed purely on observed quantities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Measuring Black Hole Light Echoes with Very Long Baseline Interferometry

Light passing near a black hole can follow multiple paths from an emission source to an observer due to strong gravitational lensing. Photons following different paths take different amounts of time to reach the observer, which produces an echo signature in the image. The characteristic echo delay is determined primarily by the mass of the black hole, but it is also influenced by the black hole spin and inclination to the observer. In the Kerr geometry, echo images are demagnified, rotated, and sheared copies of the direct image and lie within a restricted region of the image. Echo images have exponentially suppressed flux, and temporal correlations within the flow make it challenging to directly detect light echoes from the total light curve. In this Letter, we propose a novel method to search for light echoes by correlating the total light curve with the interferometric signal at high spatial frequencies, which is a proxy for indirect emission. We explore the viability of our method using numerical general relativistic magnetohydrodynamic simulations of a near-face-on accretion system scaled to M87-like parameters. We demonstrate that our method can be used to directly infer the echo delay period in simulated data. An echo detection would be clear evidence that we have captured photons that have circled the black hole, and a high-fidelity echo measurement would provide an independent measure of fundamental black hole parameters. Our results suggest that detecting echoes may be achievable through interferometric observations with a modest space-based very long baseline interferometry mission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Properties of F Stars with Stable Radial Velocity Timeseries: A Useful Metric for Selecting Low-jitter F Stars

In a companion paper, we have conducted an in-depth analysis of radial velocity jitter of over 600 stars, examining the astrophysical origins including stellar granulation, oscillation, and magnetic activity. In this paper, we highlight a subsample of those stars, specifically the main sequence and “retired” F stars—which we refer to as “MSRF” stars—that show low levels of RV jitter (<10 m s{sup −1}). We describe the observational signatures of these stars that allow them to be identified in radial velocity planet programs, for instance, those performing follow-up of transiting planets discovered by TESS. We introduce a “jitter metric” that combines the two competing effects of RV jitter with age: activity and convection. Using thresholds in the jitter metric, we can select both “complete” and “pure” samples of low jitter F stars. We also provide recipes for identifying these stars using only Gaia colors and magnitudes. Finally, we describe a region in the Gaia color–magnitude diagram where low jitter F stars are most highly concentrated. By fitting a ninth-order polynomial to the Gaia main sequence, we use the height above the main sequence as a proxy for evolution, allowing for a crude selection of low jitter MSRF stars when activity measurements are otherwise unavailable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Tracing the outflow kinematics in Type 2 active galactic nuclei

Here, we used a sample of 577 spectra of active galactic nuclei Type 1.8-2 (z < 0.25) taken from the Sloan Digital Sky Survey to trace the influence of the outflow kinematics on the profiles of different emission lines (Hβ, [O III], Hα, [N II], [S II]). All considered lines were fitted with two Gaussian components: one that fits the core of the line, and another that fits the wings. We provide a procedure for decomposition of the Hα+[N II] wavelength band for spectra where these lines overlap. The influence of the gravitational and non-gravitational kinematics on the line components is investigated by comparing the dispersions of the line components with stellar velocity dispersion. We find that wing components of all the considered emission lines have pure non-gravitational kinematics. The core components are consistent with gravitational kinematics for the Hα, [N II], and [S II] lines, while in the [O III] there is evidence for contribution from non-gravitational kinematics. We adopted the wing components as a proxy for the outflow contribution and investigated the outflow kinematics by analysing the correlations between the widths and shifts of the wing components of different lines. For this purpose, we used the subsets in which wing components are detected in both compared lines, and can be fitted independently. We find strong correlations between wing component shifts, as well as between wing component widths of all considered lines, with the exception of the Hβ wing component width. These correlations indicate that outflow dynamics systemically affect all emission lines in the spectrum. However, it reflects with a different strength in their profiles, which is observed as different widths of the wing components. This is investigated by comparison of the mean widths of the wing components in subsets where wing components are present in all lines. The strongest outflow signature is observed in the [O III] lines, which have the broadest wing components; weaker outflow signatures are found in Hα and [N II], and the weakest is found for [S II]. These results imply that the considered lines arise in different parts of an outflowing region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Integrating Intermediate Traits in Phylogenetic Genotype-to-Phenotype Studies

A major goal of research in evolution and genetics is linking genotype to phenotype. This work could be direct, such as determining the genetic basis of a phenotype by leveraging genetic variation or divergence in a developmental, physiological, or behavioral trait. The work could also involve studying the evolutionary phenomena (e.g., reproductive isolation, adaptation, sexual dimorphism, behavior) that reveal an indirect link between genotype and a trait of interest. When the phenotype diverges across evolutionarily distinct lineages, this genotype-to-phenotype problem can be addressed using phylogenetic genotype-to-phenotype (PhyloG2P) mapping, which uses genetic signatures and convergent phenotypes on a phylogeny to infer the genetic bases of traits. The PhyloG2P approach has proven powerful in revealing key genetic changes associated with diverse traits, including the mammalian transition to marine environments and transitions between major mechanisms of photosynthesis. However, there are several intermediate traits layered in between genotype and the phenotype of interest, including but not limited to transcriptional profiles, chromatin states, protein abundances, structures, modifications, metabolites, and physiological parameters. Each intermediate trait is interesting and informative in its own right, but synthesis across data types has great promise for providing a deep, integrated, and predictive understanding of how genotypes drive phenotypic differences and convergence. We argue that an expanded PhyloG2P framework (the PhyloG2P matrix) that explicitly considers intermediate traits, and imputes those that are prohibitive to obtain, will allow a better mechanistic understanding of any trait of interest. Furthermore, this approach provides a proxy for functional validation and mechanistic understanding in organisms where laboratory manipulation is impractical.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Energetic Electron Precipitation Driven by Electromagnetic Ion Cyclotron Waves from ELFIN’s Low Altitude Perspective

We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data collected by the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at >0.5 MeV which are abrupt (bursty) (lasting ~17 s, or ΔL ~ 0.56) with significant substructure (occasionally down to sub-second timescale). We attribute the bursty nature of the precipitation to the spatial extent and structuredness of the wave field at the equator. Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Case studies employing conjugate ground-based or equatorial observations of the EMIC waves reveal that the energy of moderate and strong precipitation at ELFIN approximately agrees with theoretical expectations for cyclotron resonant interactions in a cold plasma. Using multiple years of ELFIN data uniformly distributed in local time, we assemble a statistical database of ~50 events of strong EMIC wave-driven precipitation. Most reside at L ~ 5 – 7 at dusk, while a smaller subset exists at L ~ 8 – 12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio’s spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of ~ 1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. It should be noted though that this diffusive treatment likely includes effects from nonlinear resonant interactions (especially at high energies) and nonresonant effects from sharp wave packet edges (at low energies). Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven >1 MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to ~ 200-300 keV by much less intense higher frequency EMIC waves at dusk (where such waves are most frequent). At ~100 keV, whistler-mode chorus may be implicated in concurrent precipitation. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two Views of the Radius Gap and the Role of Light Curve Fitting

Recently, several groups have resolved a gap that bifurcates planets between the size of Earth and Neptune into two populations. The location and depth of this feature is an important signature of the physical processes that form and sculpt planets. In particular, planets residing in the radius gap are valuable probes of these processes as they may be undergoing the final stages of envelope loss. Here, we discuss two views of the radius gap by Fulton & Petigura (F18) and Van Eylen et al. (V18). In V18, the gap is wider and more devoid of planets. This is due, in part, to V18's more precise measurements of planet radius R{sub p}. Thanks to Gaia, uncertainties in stellar radii R{sub ⋆} are no longer the limiting uncertainties in determining R{sub p} for the majority of Kepler planets; instead, errors in R{sub p}/R{sub ⋆} dominate. V18's analysis incorporated short-cadence photometry along with constraints on mean stellar density that enabled more accurate determinations of R{sub p}/R{sub ⋆}. In the F18 analysis, less accurate R{sub p}/R{sub ⋆} blurs the boundary the radius gap. The differences in R{sub p}/R{sub ⋆} are largest at high impact parameter (b≳0.8) and often exceed 10%. This motivates excluding high-b planets from demographic studies, but identifying such planets from long-cadence photometry alone is challenging. We show that transit duration can serve as an effective proxy, and we leverage this information to enhance the contrast between the super-Earth and sub-Neptune populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗