pinecone/core-2020-05-10-deduplication huggingface dataset:
Probabilistic human health risk assessment of degradation-related chemical mixtures in heterogeneous aquifers: risk statistics, hot spots, and preferential channels

An edited version of this paper was published by AGU. Copyright (2015) American Geophysical Union.The increasing presence of toxic chemicals released in the subsurface has led to a rapid growth of social concerns and the need to develop and employ models that can predict the impact of groundwater contamination on human health risk under uncertainty. Monitored natural attenuation is a common remediation action in many contamination cases. However, natural attenuation can lead to the production of daughter species of distinct toxicity that may pose challenges in pollution management strategies. The actual threat that these contaminants pose to human health depends on the interplay between the complex structure of the geological media and the toxicity of each pollutant byproduct. This work addresses human health risk for chemical mixtures resulting from the sequential degradation of a contaminant (such as a chlorinated solvent) under uncertainty through high-resolution three-dimensional numerical simulations. We systematically investigate the interaction between aquifer heterogeneity, flow connectivity, contaminant injection model, and chemical toxicity in the probabilistic characterization of health risk. We\ud
illustrate how chemical-specific travel times control the regime of the expected risk and its corresponding uncertainties. Results indicate conditions where preferential flow paths can favor the reduction of the overall risk of the chemical mixture. The overall human risk response to aquifer connectivity is shown to be nontrivial for multispecies transport. This nontriviality is a result of the interaction between aquifer heterogeneity and chemical toxicity. To quantify the joint effect of connectivity and toxicity in health risk, we propose a toxicity-based Damk€ohler number. Furthermore, we provide a statistical characterization in terms of low-order moments and the probability density function of the individual and total risks.Peer ReviewedPostprint (published version
Multi-decadal variations in Southern Hemisphere atmospheric ¹⁴C: Evidence against a Southern Ocean sink at the end of the Little Ice Age CO₂ anomaly.

Northern Hemisphere-wide cooling during the Little Ice Age (LIA; CE 1650-1775) is associated with a ~5 ppmv decrease in atmospheric carbon dioxide. Changes in terrestrial and ocean carbon reservoirs have been postulated as possible drivers of this relatively large shift in atmospheric CO₂, potentially providing insights into the mechanisms and sensitivity of the global carbon cycle. Here we report decadally-resolved radiocarbon (¹⁴C) levels in a network of tree rings series spanning CE 1700-1950 located along the northern boundary of, and within, the Southern Ocean. We observe regional dilutions in atmospheric radiocarbon (relative to the Northern Hemisphere) associated with upwelling of ¹⁴CO₂–depleted abyssal waters. We find the inter-hemispheric ¹⁴C offset approaches zero during increasing global atmospheric CO₂ at the end of the LIA, with reduced ventilation in the Southern Ocean and a Northern Hemisphere source of old carbon (most probably originating from deep Arctic peat layers). The coincidence of the atmospheric CO₂ increase and reduction in the inter-hemispheric ¹⁴C offset imply a common climate control. Possible mechanisms of synchronous change in the high latitudes of both hemispheres are discussed
An efficient and general approach for implementing thermodynamic phase equilibria information in geophysical and geodynamic studies

We present a flexible, general, and efficient approach for implementing thermodynamic phase equilibria information (in the form of sets of physical parameters) into geophysical and geodynamic studies. The approach is based on Tensor Rank Decomposition methods, which transform the original multidimensional discrete information into a separated representation that contains significantly fewer terms, thus drastically reducing the amount of information to be stored in memory during a numerical simulation or geophysical inversion. Accordingly, the amount and resolution of the thermodynamic information that can be used in a simulation or inversion increases substantially. In addition, the method is independent of the actual software used to obtain the primary thermodynamic information, and therefore, it can be used in conjunction with any thermodynamic modeling program and/or database. Also, the errors associated with the decomposition procedure are readily controlled by the user, depending on her/his actual needs (e.g., preliminary runs versus full resolution runs). We illustrate the benefits, generality, and applicability of our approach with several examples of practical interest for both geodynamic modeling and geophysical inversion/modeling. Our results demonstrate that the proposed method is a competitive and attractive candidate for implementing thermodynamic constraints into a broad range of geophysical and geodynamic studies. MATLAB implementations of the method and examples are provided as supporting information and can be downloaded from the journal's website.Peer ReviewedPostprint (author's final draft
Erosion rate study at the Allchar deposit (Macedonia) based on radioactive and stable cosmogenic nuclides (26Al, 36Cl, 3He, and 21Ne)

This paper focuses on constraining the erosion rate in the area of the Allchar Sb-As-Tl-Au deposit (Macedonia). It contains the largest known reserves of lorandite (TlAsS2), which is essential for the LORanditeEXperiment (LOREX), aimed at determining the long-term solar neutrino flux. Because the erosion history of the Allchar area is crucial for the success of LOREX, we applied terrestrial in situ cosmogenic nuclides including both radioactive (26Al and 36Cl) and stable (3He and 21Ne) nuclides in quartz, dolomite/calcite, sanidine, and diopside. The obtained results suggest that there is accordance in the values obtained by applying 26Al, 36Cl, and 21Ne for around 85% of the entire sample collection, with resulting erosion rates varying from several tens of m/Ma to ∼165 m/Ma. The samples from four locations (L-8 CD, L1b/R, L1c/R, and L-4/ADR) give erosion rates between 300 and 400 m/Ma. Although these localities reveal remarkably higher values, which may be explained by burial events that occurred in part of Allchar, the erosion rate estimates mostly in the range between 50 and 100 m/Ma. This range further enables us to estimate the vertical erosion rate values for the two main ore bodies Crven Dol and Centralni Deo. We also estimate that the lower and upper limits of average paleo-depths for the ore body Centralni Deo from 4.3 Ma to the present are 250–290 and 750–790 m, respectively, whereas the upper limit of paleo-depth for the ore body Crven Dol over the same geological age is 860 m. The estimated paleo-depth values allow estimating the relative contributions of 205Pb derived from pp-neutrino and fast cosmic-ray muons, respectively, which is an important prerequisite for the LOREX experiment
Northeast Indian stalagmite records Pacific decadal climate change: Implications for moisture transport and drought in India

This is the final version. It is currently under embargo. It was first published by Wiley at http://onlinelibrary.wiley.com/doi/10.1002/2015GL063826/full.Two types of El Niño events are distinguished by sea surface temperature (SST) anomalies
centered in the central or eastern equatorial Pacific. The Central Pacific El Niño events (CP-El Niño) are
more highly correlated with weakening of the central Indian Summer Monsoon and linked to decadal Pacific
climate variability. We present a 50 year, subannually resolved speleothem δ18O record from northeast India
that exhibits a significant correlation with northern Pacific decadal variability and central equatorial Pacific
SSTs. Accordingly, we suggest that δ18O time series in similar northeast Indian speleothems are effective tools
for investigating preinstrumental changes in Pacific climate, including changes in El Niño dynamics. In
contrast to central India, rainfall amounts in northeast India are relatively unaffected by El Niño. However,
back trajectory analysis indicates that during CP-El Niño events moisture transport distance to northeast India
is reduced, suggesting that variations in moisture transport primarily control δ18O in the region.This work was supported through the BanglaPIRE project (NSF OISE-0968354), an award from the Vanderbilt International Office to JLO and SFMB, and awards from the Cave Research Foundation and the Geological Society of America to CGM. SFMB received financial support from the Schweizer National Fond (SNF), Sinergia grant CRSI22 132646/1
Dissolution patterns and mixing dynamics in unstable reactive flow

We study the fundamental problem of mixing and chemical reactions under a Rayleigh-Taylor-type hydrodynamic instability in a miscible two-fluid system. The dense fluid mixture, which is generated at the fluid-fluid interface, leads to the onset of a convective fingering instability and triggers a fast chemical dissolution reaction. Contrary to intuition, the dissolution pattern does not map out the finger geometry. Instead, it displays a dome-like, hierarchical structure that follows the path of the ascending fluid interface and the regions of maximum mixing. These mixing and reaction hot spots coincide with the flow stagnation points, at which the interfacial mixing layer is compressed and deformed. We show that the deformation of the boundary layer around the stagnation points controls the evolution of the global scalar dissipation and reaction rates and shapes the structure of the reacted zones. The persistent compression of the mixing layer explains the independence of the mixing rate from the Rayleigh number when convection dominates.Peer ReviewedPostprint (author's final draft
Two possible source regions for central Greenland last glacial dust

Dust in Greenland ice cores is used to reconstruct the activity of dust-emitting regions and atmospheric circulation. However, the source of dust material to Greenland over the last glacial period is the subject of considerable uncertainty. Here we use new clay mineral and <10 µm Sr–Nd isotopic data from a range of Northern Hemisphere loess deposits in possible source regions alongside existing isotopic data to show that these methods cannot discriminate between two competing hypothetical origins for Greenland dust: an East Asian and/or central European source. In contrast, Hf isotopes (<10 µm fraction) of loess samples show considerable differences between the potential source regions. We attribute this to a first-order clay mineralogy dependence of Hf isotopic signatures in the finest silt/clay fractions, due to absence of zircons. As zircons would also be absent in Greenland dust, this provides a new way to discriminate between hypotheses for Greenland dust sources
Rapid and highly variable warming of lake surface waters around the globe

In this first worldwide synthesis of in situ and satellite-derived lake data, we find that lake summer surface water temperatures rose rapidly (global mean = 0.34°C decade-1) between 1985 and 2009. Our analyses show that surface water warming rates are dependent on combinations of climate and local characteristics, rather than just lake location, leading to the counterintuitive result that regional consistency in lake warming is the exception, rather than the rule. The most rapidly warming lakes are widely geographically distributed, and their warming is associated with interactions among different climatic factors - from seasonally ice-covered lakes in areas where temperature and solar radiation are increasing while cloud cover is diminishing (0.72°C decade-1) to ice-free lakes experiencing increases in air temperature and solar radiation (0.53°C decade-1). The pervasive and rapid warming observed here signals the urgent need to incorporate climate impacts into vulnerability assessments and adaptation efforts for lakes
Lightning as a space-weather hazard: UK thunderstorm activity modulated by the passage of the heliospheric current sheet

Lightning flash rates, RL, are modulated by corotating interaction regions (CIRs) and the polarity of the heliospheric magnetic field (HMF) in near-Earth space. As the HMF polarity reverses at the heliospheric current sheet (HCS), typically within a CIR, these phenomena are likely related. In this study, RL is found to be significantly enhanced at the HCS and at 27 days prior/after. The strength of the enhancement depends on the polarity of the HMF reversal at the HCS. Near-Earth solar and galactic energetic particle fluxes are also ordered by HMF polarity, though the variations qualitatively differ from RL, with the main increase occurring prior to the HCS crossing. Thus, the CIR effect on lightning is either the result of compression/amplification of the HMF (and its subsequent interaction with the terrestrial system) or that energetic particle preconditioning of the Earth system prior to the HMF polarity change is central to solar wind lightning coupling mechanism
Hybrid simulation of Titan's interaction with the supersonic solar wind during Cassini's T96 flyby

By applying a hybrid (kinetic ions and fluid electrons) simulation code, we study the plasma environment of Saturn´s largest moon Titan during Cassini´s T96 flyby on 1 December 2013. The T96 encounter marks the only observed event of the entire Cassini mission where Titan was located in the supersonic solar wind in front of Saturn´s bow shock. Our simulations can quantitatively reproduce the key features of Cassini magnetic field and electron density observations during this encounter. We demonstrate that the large-scale features of Titan´s induced magnetosphere during T96 can be described in terms of a steady state interaction with a high-pressure solar wind flow. About 40 min before the encounter, Cassini observed a rotation of the incident solar wind magnetic field by almost 90Â°. We provide strong evidence that this rotation left a bundle of fossilized magnetic field lines in Titan´s ionosphere that was subsequently detected by the spacecraft.Fil: Feyerabend, Moritz. Georgia Institute Of Techology; Estados UnidosFil: Simon, Sven. Georgia Institute Of Techology; Estados UnidosFil: Neubauer, Fritz M.. Universitat Zu Köln; AlemaniaFil: Motschmann, Uwe. Deutsches Zentrum Fur Luft- Und Raumfahrt; Alemania. Technische Universitat Braunschweig; AlemaniaFil: Bertucci, Cesar. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: Edberg, Niklas J. T.. Instiutet For Rymdfysik; SueciaFil: Hospodarsky, George B.. University Of Iowa; Estados UnidosFil: Kurth, William S.. University Of Iowa; Estados Unido
Attribution of extreme weather and climate events overestimated by unreliable climate simulations

Event attribution aims to estimate the role of an external driver after the occurrence of an extreme weather and climate event by comparing the probability that the event occurs in two counterfactual worlds.\ud
These probabilities are typically computed using ensembles of climate simulations whose simulated probabilities are known to be imperfect. The implications of using imperfect models in this context are largely unknown, limited by the number of observed extreme events in the past to conduct a robust evaluation. Using an idealized framework, this model limitation is studied by generating large number of simulations with variable reliability in simulated probability. The framework illustrates that unreliable climate simulations are prone to overestimate the attributable risk to climate change. Climate model ensembles tend to be overconfident in their representation of the climate variability which leads to systematic increase in the attributable risk to an extreme event. Our results suggest that event attribution approaches comprising of a single climate model would benefit from ensemble calibration in order to account for model inadequacies\ud
similarly as operational forecasting systems.We would like to acknowledge valuable discussions and feedback received from\ud
François Massonnet, Nathalie Schaller,\ud
Chloé Prodhomme, and Fraser Lott. This\ud
work was supported by the EUropean\ud
CLimate and weather Events:\ud
Interpretation and Attribution (EUCLEIA),\ud
funded by the European Union’s Seventh\ud
Framework Programme [FP7/2007-2013]\ud
under grant agreement 607085 and the\ud
ESA Living Planet Fellowship Programme\ud
under the project VERITAS-CCI. We are\ud
further indebted to the s2dverification\ud
(http://cran.r-project.org/web/packages/\ud
SpecsVerification/index.html) and specsverification\ud
(http://cran.r-project.org/\ud
web/packages/s2dverification/index.\ud
html) with which the calculations have\ud
been carried out. The synthetic hindcast\ud
generator has been implemented into\ud
s2dverfiction. No further data were used\ud
in producing this manuscript.Peer ReviewedPostprint (published version
Dependence of the location of the Martian magnetic lobes on the interplanetary magnetic field direction: Observations from Mars Global Surveyor

We use magnetometer data from the Mars Global Surveyor (MGS) spacecraft during portions of the premapping orbits of the mission to study the variability of the Martian-induced magnetotail as a function of the orientation of the interplanetary magnetic field (IMF). The time spent by MGS in the magnetotail lobes during periods with positive solar wind flow-aligned IMF component B IMF suggests that their location as well as the position of the central polarity reversal layer (PRL) are displaced in the direction antiparallel to the IMF cross-flow component B IMF. Analogously, in the cases where B IMF is negative, IMF the lobes are displaced in the direction of B ⟂ . This behavior is compatible with a previously published analytical model of the IMF draping, where for the first time, the displacement of a complementary reversal layer (denoted as IPRL for inverse polarity reversal layer) is deduced from first principles.Fil: Romanelli, Norberto Julio. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: Bertucci, Cesar. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: Gomez, Daniel Osvaldo. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: Mazelle, C.. Université Paul Sabatier; Franci
Effects of Saturn's magnetospheric dynamics on Titan's ionosphere

We use the Cassini Radio and Plasma Wave Science/Langmuir probe measurements of the electron density from the first 110 flybys of Titan to study how Saturn´s magnetosphere influences Titan´s ionosphere. The data is first corrected for biased sampling due to varying solar zenith angle and solar energy flux (solar cycle effects). We then present results showing that the electron density in Titan´s ionosphere, in the altitude range 1600-2400 km, is increased by about a factor of 2.5 when Titan is located on the nightside of Saturn (Saturn local time (SLT) 21-03 h) compared to when on the dayside (SLT 09-15 h). For lower altitudes (1100-1600 km) the main dividing factor for the ionospheric density is the ambient magnetospheric conditions. When Titan is located in the magnetospheric current sheet, the electron density in Titan´s ionosphere is about a factor of 1.4 higher compared to when Titan is located in the magnetospheric lobes. The factor of 1.4 increase in between sheet and lobe flybys is interpreted as an effect of increased particle impact ionization from 200 eV sheet electrons. The factor of 2.5 increase in electron density between flybys on Saturn´s nightside and dayside is suggested to be an effect of the pressure balance between thermal plus magnetic pressure in Titan´s ionosphere against the dynamic pressure and energetic particle pressure in Saturn´s magnetosphere.Fil: Edberg, N. J. T.. University of Iowa; Estados Unidos. Swedish Institute of Space Physics; SueciaFil: Andrews, D. J.. Swedish Institute of Space Physics; SueciaFil: Bertucci, Cesar. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: Gurnett, D. A.. University of Iowa; Estados UnidosFil: Holmberg, M. K. G.. Swedish Institute of Space Physics; SueciaFil: Jackman, C. M.. University Of Southampton; Reino UnidoFil: Kurth, W. S.. University of Iowa; Estados UnidosFil: Menietti, J. D.. University Of Iowa; Estados UnidosFil: Opgenoorth, H. J.. Swedish Institute of Space Physics; SueciaFil: Shebanits, O.. Swedish Institute of Space Physics; SueciaFil: Vigren, E.. Swedish Institute of Space Physics; SueciaFil: Wahlund, J. E.. Swedish Institute of Space Physics; Sueci
GPS as a solar observational instrument: real-time estimation of EUV photons flux rate during strong, medium, and weak solar flares

In this manuscript, the authors show how the Global Navigation Satellite Systems, GNSS (exemplified in the Global Positioning System, GPS), can be efficiently used for a very different purpose from that for which it was designed as an accurate Solar observational tool, already operational from the open global GPS measurements available in real-time, and with some advantages regarding dedicated instruments onboard spacecraft. The very high correlation of the solar extreme ultraviolet (EUV) photon flux rate in the 26–34 mm spectral band, obtained from the solar EUV monitor instrument onboard the SOHO spacecraft during Solar flares, is shown with the GNSS solar flare activity indicator (GSFLAI). The GSFLAI is defined as the gradient of the ionospheric vertical total electron content rate versus the cosine of the Solar zenith angle in the day hemisphere (which filters out nonsolar over ionization), and it is measured from data collected by a global network of dual frequency GPS receivers (giving in this way continuous coverage). GSFLAI for 60 X class flares, 320 M class flares, and 300 C class flares, occurred since 2001, were directly compared with the EUV solar flux rate data to show existing correlations. It was found that the GSFLAI and EUV flux rate present the same linear relationship for all classes of flares, not only the strong and medium intensity ones, X and M class, as in previous works, but also for the weakest C class solar flares, which is a remarkable result.Peer ReviewedPostprint (author's final draft