Proton flux and radiation dose from galactic cosmic rays in the lunar regolith and implications for organic synthesis at the poles of the Moon and Mercury

Galactic cosmic rays are a potential energy source to stimulate organic synthesis from simple ices. The recent detection of organic molecules at the polar regions of the Moon by LCROSS (Colaprete, A. et al. [2010]. Science 330, 463–468, http://dx.doi.org/10.1126/science.1186986), and possibly at the poles of Mercury (Paige, D.A. et al. [2013]. Science 339, 300–303, http://dx.doi.org/10.1126/science.1231106), introduces the question of whether the organics were delivered by impact or formed in situ. Laboratory experiments show that high energy particles can cause organic production from simple ices. We use a Monte Carlo particle scattering code (MCNPX) to model and report the flux of GCR protons at the surface of the Moon and report radiation dose rates and absorbed doses at the Moon’s surface and with depth as a result of GCR protons and secondary particles, and apply scaling factors to account for contributions to dose from heavier ions. We compare our results with dose rate measurements by the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) experiment on Lunar Reconnaissance Orbiter (Schwadron, N.A. et al. [2012]. J. Geophys. Res. 117, E00H13, http://dx.doi.org/10.1029/2011JE003978) and find them in good agreement, indicating that MCNPX can be confidently applied to studies of radiation dose at and within the surface of the Moon. We use our dose rate calculations to conclude that organic synthesis is plausible well within the age of the lunar polar cold traps, and that organics detected at the poles of the Moon may have been produced in situ. Our dose rate calculations also indicate that galactic cosmic rays can induce organic synthesis within the estimated age of the dark deposits at the pole of Mercury that may contain organics.     

Oxygen and hydrogen isotope fractionation in serpentine-water and talc-water systems from 250 to 450 °C, 50 MPa

Oxygen and hydrogen isotope fractionation factors in the talc-water and serpentine-water systems have been determined by laboratory experiment from 250 to 450 °C at 50 MPa using the partial exchange technique. Talc was synthesized from brucite + quartz, resulting in nearly 100% exchange during reaction at 350 and 450 °C. For serpentine, D-H exchange was much more rapid than ¹⁸O-¹⁶O exchange when natural chrysotile fibers were employed in the initial charge. In experiments with lizardite as the starting charge, recrystallization to chrysotile enhanced the rate of ¹⁸O-¹⁶O exchange with the coexisting aqueous phase. Oxygen isotope fractionation factors in both the talc-water and serpentine-water systems decrease with increasing temperature and can be described from 250 to 450 °C by the relationships: 1000 ln  = 11.70 × 10⁶/T² − 25.49 × 10³/T + 12.48 and 1000 ln  = 3.49 × 10⁶/T² − 9.48 where T is temperature in Kelvin. Over the same temperature interval at 50 MPa, talc-water D-H fractionation is only weakly dependent on temperature, similar to brucite and chlorite, and can be described by the equation: 1000 ln = 10.88 × 10⁶/T² − 41.52 × 10³/T + 5.61 where T is temperature in Kelvin. Our D-H serpentine-water fractionation factors calibrated by experiment decrease with temperature and form a consistent trend with fractionation factors derived from lower temperature field calibrations. By regression of these data, we have refined and extended the D-H fractionation curve from 25 to 450 °C, 50 MPa as follows: 1000 ln  = 3.436 × 10⁶/T² − 34.736 × 10³/T + 21.67 where T is temperature in Kelvin. These new data should improve the application of D-H and ¹⁸O-¹⁶O isotopes to constrain the temperature and origin of hydrothermal fluids responsible for serpentine formation in a variety of geologic settings.     

N and O isotope effects during nitrate assimilation by unicellular prokaryotic and eukaryotic plankton cultures

In order to provide biological systematics from which to interpret nitrogen (N) and oxygen (O) isotope ratios of nitrate (¹⁵N/¹⁴N, ¹⁸O/¹⁶O, respectively) in the environment, we previously investigated the isotopic fractionation of nitrate during its assimilation by mono-cultures of eukaryotic algae (Granger et al., 2004). In this study, we extended our analysis to investigate nitrate assimilation by strains of prokaryotic plankton. We measured the N and O isotope effects, ¹⁵ε and ¹⁸ε, during nitrate consumption by cultures of prokaryotic strains and by additional eukaryotic phytoplankton strains (where ε is the ratio of reaction rate constants of the light vs. heavy isotopologues, ^lightk and ^heavyk; ε = ^lightk/^heavyk − 1 × 1000, expressed in per mil). The observed ¹⁵ε ranged from 5‰ to 8‰ among eukaryotes, whereas it did not exceed 5‰ for three cyanobacterial strains, and was as low as 0.4‰ for a heterotrophic α-protoeobacterium. Eukaryotic phytoplankton fractionated the N and O isotopes of nitrate to the same extent (i.e., ¹⁸ε ∼ ¹⁵ε). The ¹⁸ε:¹⁵ε among the cyanobacteria was also ∼1, whereas the heterotrophic α-proteobacterial strain, which showed the lowest ¹⁵ε, between 0.4‰ and 1‰, had a distinct ¹⁸ε:¹⁵ε of ∼2, unlike any plankton strain observed previously. Equivalent N vs. O isotope discrimination is thought to occur during internal nitrate reduction by nitrate reductase, such that the cellular efflux of the fractionated nitrate into the medium drives the typically observed ¹⁸ε:¹⁵ε of ∼1. We hypothesize that the higher in the ¹⁸ε:¹⁵ε of the α-proteobacterium may result from isotope discrimination by nitrate transport, which is evident only at low amplitude of ε. These observations warrant investigating whether heterotrophic bacterial assimilation of nitrate decreases the community isotope effects at the surface ocean.     

1) I. A. H. S. / A. I. H. S.

Abstract

A simple method is used to study the response of runoff in the Sahel to climate change. The statistical characteristics of rainfall are calculated over the western part of the Sahel for the period 1961–1990, using the BADOPLU network. Daily rainfall is simulated using a Markov process with Weibull distribution for rainfall depths. Runoff is modelled using a conceptual SCS model and the curve numbers are calculated for West Africa. Climate change is provided by simulations using the Arpège GCM (Scenario A1B), and a perturbation method is used on the parameters which describe the rainfall. Changes in rainfall are assumed to occur through increases in frequency, not intensity. Using Arpège, runoff is mainly found to increase, in depth and in number of events, by the end of the 21st century. Changes in evaporation and land use are not included in the analysis. The impact of this 21st century potential climate change (rainfall) on the runoff is found to be of the same magnitude as the impact of changes in land use.     

Sedimentological constraints on hydrometeor-enhanced particle deposition: 1992 Eruptions of Crater Peak, Alaska

Water is a dominant component of volcanic clouds and has fundamental control on very fine particle deposition. Particle size characteristics of distal tephra-fall (100s km from source volcano) have a higher proportion of very fine particles compared to predictions based on single particle settling rates. In this study, sedimentological analyses of fallout from for the 18 August and 16–17 September 1992 eruptions of Crater Peak, Alaska, are combined with satellite observations, and cloud trajectory and microphysics modeling to investigate meteorological influences on particle sedimentation. Total grain size distributions of tephra fallout were reconstructed for both Crater Peak eruptions and indicate a predominance of fine particles < 125 μm. Polymodal analysis of the deposits has identified a particle subpopulation with mode ~ 15–18 μm involved in particle aggregation. Accounting for the magmatic water source only, calculated ice water content of the 3.7 hour old September 1992 Spurr cloud was ~ 4.5 × 10^− 2 g m^− 3 (based on an estimated cloud thickness of ~ 1000 m from trajectory modeling). Hydrometeor formation on particles in the volcanic cloud and subsequent sublimation may induce a cloud base instability that leads to rapid bulk (en masse) sedimentation of very fine particles through a mammatus-like mechanism.     

The role of prior model calibration on predictions with ensemble Kalman filter

This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.     

A regional model for studies of atmosphere-ice-ocean interaction in the western Arctic

Summary Asa step in the development of a fully coupled regional model of the atmosphere-ice-ocean system, atmospheric and sea ice models have been adapted to a western Arctic domain centered on the Bering Strait. Lateral boundary conditions derived from operational analyses drive the models through simulations on grids having horizontal resolutions of 21 km and 7 km. Sensitivities to the presence of sea ice are large after only 48 hours, by which time the surface temperatures in the Bering and Chukchi Seas are 10–15°C higher without sea ice than with sea ice. The temperatures, in turn, modify the fields of sea level pressure, surface wind and precipitation. By influencing the surface wind stress through the static static stability, the surface state feeds back to the surface momentum exchange, ice/ocean transport, and the rate of formation of new ice. The results also show a resolution-dependence of the surface winds, precipitation rates and new ice formation rates, particularly in areas in which the coastal configuration and topography are spatially complex. The experiments will be augmented by the implementation of an ocean model on the same grids.With 12 Figures