Tracer applications of noble gas radionuclides in the geosciences

Noble gas radionuclides, including 81Kr (t1/2 = 229,000 years), 85Kr (t1/2 = 10.8 years), and 39Ar (t1/2 = 269 years), possess nearly ideal chemical and physical properties for studies of earth and environmental processes. Recent advances in Atom Trap Trace Analysis (ATTA), a laser-based atom counting method, have enabled routine measurements of the radiokrypton isotopes, as well as the demonstration of the ability to measure 39Ar in environmental samples. Here we provide an overview of the ATTA technique, and a survey of recent progress made in several laboratories worldwide.We review the application of noble gas radionuclides in the geosciences and discuss how ATTA can help advance these fields, specifically: determination of groundwater residence times using 81Kr, 85Kr, and 39Ar; dating old glacial ice using 81Kr; and an 39Ar survey of the main water masses of the oceans, to study circulation pathways and estimate mean residence times. Other scientific questions involving a deeper circulation of fluids in the Earth's crust and mantle are also within the scope of future applications. We conclude that the geoscience community would greatly benefit from an ATTA facility dedicated to this field, with instrumentation for routine measurements, as well as for research on further development of ATTA methods.

Gas-phase Lifetimes of Nucleobase Analogues by Picosecond Pumpionization and Streak Techniques

The picosecond (ps) timescale is relevant for the investigation of many molecular dynamical processes such as fluorescence, nonradiative relaxation, intramolecular vibrational relaxation, molecular rotation and intermolecular energy transfer, to name a few. While investigations of ultrafast (femtosecond) processes of biological molecules, e.g. nucleobases and their analogues in the gas phase are available, there are few investigations on the ps time scale. We have constructed a ps pump-ionization setup and a ps streak camera fluorescence apparatus for the determination of lifetimes of supersonic jet-cooled and isolated molecules and clusters. The ps pump-ionization setup was used to determine the lifetimes of the nucleobase analogue 2-aminopurine (2AP) and of two 2AP˙(H2O)n water cluster isomers with n=1 and 2. Their lifetimes lie between 150 ps and 3 ns and are strongly cluster-size dependent. The ps streak camera setup was used to determine accurate fluorescence lifetimes of the uracil analogue 2-pyridone (2PY), its self-dimer (2PY)2, two isomers of its trimer (2PY)3 and its tetramer (2PY)4, which lie in the 7–12 ns range.

Mass transport around comets and its impact on the seasonal differences in water production rates

Comets are surrounded by a thin expanding atmosphere, and although the nucleus' gravity is small, some molecules and grains, possibly with the inclusion of ices, can get transported around the nucleus through scattering (atoms/molecules) and gravitational pull (grains). Based on the obliquity of the comet, it is also possible that volatile material and icy grains get trapped in regions, which are in shadow until the comet passes its equinox. When the Sun rises above the horizon and the surface starts to heat up, this condensed material starts to desorb and icy grains will sublimate off the surface, possibly increasing the comet's neutral gas production rate on the outbound path. In this paper we investigate the mass transport around the nucleus, and based on a simplified model, we derive the possible contribution to the asymmetry in the seasonal gas production rate that could arise from trapped material released from cold areas once they come into sunlight. We conclude that the total amount of volatiles retained by this effect can only contribute up to a few percent of the asymmetry observed in some comets.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Water-use efficiency and transpiration across European forests during the Anthropocene

The Earth’s carbon and hydrologic cycles are intimately coupled by gas exchange through plant stomata1, 2, 3. However, uncertainties in the magnitude4, 5, 6 and consequences7, 8 of the physiological responses9, 10 of plants to elevated CO2 in natural environments hinders modelling of terrestrial water cycling and carbon storage11. Here we use annually resolved long-term δ13C tree-ring measurements across a European forest network to reconstruct the physiologically driven response of intercellular CO2 (Ci) caused by atmospheric CO2 (Ca) trends. When removing meteorological signals from the δ13C measurements, we find that trees across Europe regulated gas exchange so that for one ppmv atmospheric CO2 increase, Ci increased by ~0.76 ppmv, most consistent with moderate control towards a constant Ci/Ca ratio. This response corresponds to twentieth-century intrinsic water-use efficiency (iWUE) increases of 14 ± 10 and 22 ± 6% at broadleaf and coniferous sites, respectively. An ensemble of process-based global vegetation models shows similar CO2 effects on iWUE trends. Yet, when operating these models with climate drivers reintroduced, despite decreased stomatal opening, 5% increases in European forest transpiration are calculated over the twentieth century. This counterintuitive result arises from lengthened growing seasons, enhanced evaporative demand in a warming climate, and increased leaf area, which together oppose effects of CO2-induced stomatal closure. Our study questions changes to the hydrological cycle, such as reductions in transpiration and air humidity, hypothesized to result from plant responses to anthropogenic emissions.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Kinetic simulation of neutral∕ionized gas and electrically charged dust in the coma of comet 67P∕Churyumov-Gerasimenko

The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near‐nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus. The Rosetta spacecraft is en route to comet 67P/Churyumov‐Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet’s dusty gas environment. In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov‐Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation [1] of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.

Understanding measured water rotational temperatures and column densities in the very innermost coma of Comet 73P/Schwassmann–Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J., 685, 659−677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J., 732) to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann–Wachmann 3 obtained with a very high spatial resolution of ∼30 km using IRCS/Subaru in May 2006 (Bonev, B.P., Mumma, M.J., Kawakita, H., Kobayashi, H., Villanueva, G.L. [2008]. Icarus, 196, 241−248). For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10–20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.

Effect of the Tiger Stripes on the water vapor distribution in Enceladus' exosphere

The jet activity emanating from Enceladus' exosphere south pole region observed by Cassini is a subject of intensive study. The in situ and remote sensing observations performed since 2005 triggered an active modeling campaign. Such modeling is essential for better understanding of the measurements performed by individual instruments as well as to link them for a more complete picture of the volatile and ice grain distribution in Enceladus' exosphere. This paper is focused on the investigation of the effect that diffuse gas sources along the Tiger Stripes have on distribution of the water vapor in Enceladus' exosphere using the updated version of our multiplume model. We have found that accounting for the gas production by Tiger Stripes is critical for interpretation of the Cassini data. According to our calculations, sources along the Tiger Stripes (apart from those originally identified by Spitale and Porco (2007)) must contribute about 23–32% to the total plume source rate, which varies in the range of (6.4–29) ×1027 s−1. The effect of the previously unidentified source suggested in the paper is found to be critical for explaining the Ultraviolet Imaging Spectrograph 2007 and 2010 observations in the whole range of the elapsed times.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Sublimation of water ice mixed with silicates and tholins: Evolution of surface texture and reflectance spectra, with implications for comets

The surfaces of many objects in the Solar System comprise substantial quantities of water ice sometimes mixed with minerals and/or organic molecules. The sublimation of the ice changes the structural and optical properties of these objects. We present laboratory data on the evolution of the structure and the visible and near-infrared spectral reflectance of icy surface analogues of cometary ices, made of water ice, complex organic matter (tholins) and silicates, as they undergo sublimation under low temperature (<-70°C) and pressure (10-⁵mbar) conditions inside the SCITEAS simulation chamber. As the water ice sublimated, we observed in situ the formation of a porous sublimation lag deposit, or sublimation mantle, at the top of the ice. This mantle is a network of filaments made of the non-volatile particles. Organics or phyllosilicates grains, able to interact via stronger inter-particulate forces than olivine grains, can form a foam-like structure having internal cohesiveness, holding olivine grains together. As this mantle builds-up, the band depths of the sub-surface water ice are attenuated until complete extinction under only few millimeters of mantle. Optically thick sublimation mantles are mainly featureless in the near infrared. The absorption bands of the minerals present in the mantle are weak, or even totally absent if minerals are mixed with organics which largely dominate the VIS–NIR reflectance spectrum. During sublimation, ejections of large fragments of mantle, triggered by the gas flow, expose ice particles to the surface. The contrast of brightness between mantled and ice-exposed areas depends on the wavelength range and the dust/ice ratio considered. We describe how the chemical nature of the non-volatiles, the size of their particles, the way they are mixed with the ice and the dust/ice mass ratio influence the texture, activity and spectro-photometric properties of the sublimation mantles. These data provide useful references for interpreting remote-sensing observations of comets and also icy satellites or trans-neptunian objects.