The Galactic Cosmic Ray Intensity over the Past 106–109 Years as Recorded by Cosmogenic Nuclides in Meteorites and Terrestrial Samples

Concentrations of stable and radioactive nuclides produced by cosmic ray particles in meteorites allow us to track the long term average of the primary flux of galactic cosmic rays (GCR). During the past ∼10 Ma, the average GCR flux remained constant over timescales of hundreds of thousands to millions of years, and, if corrected for known variations in solar modulation, also during the past several years to hundreds of years. Because the cosmic ray concentrations in meteorites represent integral signals, it is difficult to assess the limits of uncertainty of this statement, but they are larger than the often quoted analytical and model uncertainties of some 30%. Time series of concentrations of the radionuclide 10Be in terrestrial samples strengthen the conclusions drawn from meteorite studies, indicating that the GCR intensity on a ∼0.5 million year scale has remained constant within some ±10% during the past ∼10 million years. The very long-lived radioactive nuclide 40K allows to assess the GCR flux over about the past one billion years. The flux over the past few million years has been the same as the longer-term average in the past 0.5–1 billion years within a factor of ∼1.5. However, newer data do not confirm a long-held belief that the flux in the past few million years has been higher by some 30–50% than the very long term average. Neither does our analysis confirm a hypothesis that the iron meteorite data indicate a ∼150 million year periodicity in the cosmic ray flux, possibly related to variations in the long-term terrestrial climate.

Liulin silicon semiconductor spectrometers as cosmic ray monitors at the high mountain observatories Jungfraujoch and Lomnický stit

Currently, most cosmic ray data are obtained by detectors on satellites, aircraft, high-altitude balloons and ground (neutron monitors). In our work, we examined whether Liulin semiconductor spectrometers (simple silicon planar diode detectors with spectrometric properties) located at high mountain observatories could contribute new information to the monitoring of cosmic rays by analyzing data from selected solar events between 2005 and 2013. The decision thresholds and detection limits of these detectors placed at Jungfraujoch (Switzerland; 3475 m a.s.l.; vertical cut-off rigidity 4.5 GV) and Lomnicky stıt (Slovakia; 2633 m a.s.l.; vertical cut-off rigidity 3.84 GV) highmountain observatories were determined. The data showed that only the strongest variations of the cosmic ray flux in this period were detectable. The main limitation in the performance of these detectors is their small sensitive volume and low sensitivity of the PIN photodiode to neutrons.

Cosmic-ray exposure ages of pallasites

We analyzed cosmogenic nuclides in metal and/or silicate (primarily olivine) separated from the main-group pallasites Admire, Ahumada, Albin, Brahin, Brenham, Esquel, Finmarken, Glorieta Mountain, Huckitta, Imilac, Krasnojarsk, Marjalahti, Molong, Seymchan, South Bend, Springwater, and Thiel Mountains and from Eagle Station. The metal separates contained an olivine fraction which although small, <1 wt% in most cases, nonetheless contributes significantly to the budgets of some nuclides (e.g., up to 35% for Ne-21 and Al-26). A correction for olivine is therefore essential and was made using model calculations and/or empirical relations for the production rates of cosmogenic nuclides in iron meteoroids and/or measured elemental concentrations. Cosmic-ray exposure (CRE) ages for the metal phases of the main-group pallasites range from 7 to 180 Ma, but many of the ages cluster around a central peak near 100 Ma. These CRE ages suggest that the parent body of the main-group pallasites underwent a major break-up that produced most of the meteorites analyzed. The CRE age distribution for the pallasites overlaps only a small fraction of the distribution for the IIIAB iron meteorites. Most pallasites and IIIAB irons originated in different collisions, probably on different parent bodies; a few IIIABs and pallasites may have come out of the same collision but a firm conclusion requires further study. CRE ages calculated from noble gas and radionuclide data of the metal fraction are higher on average than the Ne-21 exposure ages obtained for the olivine samples. As the metal and olivine fractions were taken in most cases from different specimens, the depth-dependency of the production rate ratio Be-10/Ne-21 in metal, not accounted for in our calculations, may explain the difference.

Calibration and analysis of the 2007 data

he physics program of the NA61/SHINE (SHINE = SPS Heavy Ion and Neutrino Experiment) experiment at the CERN SPS consists of three subjects. In the first stage of data taking (2007-2009) measurements of hadron production in hadron-nucleus interactions needed for neutrino (T2K) and cosmic-ray (Pierre Auger and KASCADE) experiments will be performed. In the second stage (2009-2010) hadron production in proton-proton and proton-nucleus interactions needed as reference data for a better understanding of nucleus-nucleus reactions will be studied. In the third stage (2009-2013) energy dependence of hadron production properties will be measured in p+p, p+Pb interactions and nucleus-nucleus collisions, with the aim to identify the properties of the onset of deconfinement and find evidence for the critical point of strongly interacting matter. The NA61 experiment was approved at CERN in June 2007. The first pilot run was performed during October 2007. Calibrations of all detector components have been performed successfully and preliminary uncorrected spectra have been obtained. High quality of track reconstruction and particle identification similar to NA49 has been achieved. The data and new detailed simulations confirm that the NA61 detector acceptance and particle identification capabilities cover the phase space required by the T2K experiment. This document reports on the progress made in the calibration and analysis of the 2007 data.

Forcing of stratospheric chemistry and dynamics during the Dalton Minimum

The response of atmospheric chemistry and dynamics to volcanic eruptions and to a decrease in solar activity during the Dalton Minimum is investigated with the fully coupled atmosphere–ocean chemistry general circulation model SOCOL-MPIOM (modeling tools for studies of SOlar Climate Ozone Links-Max Planck Institute Ocean Model) covering the time period 1780 to 1840 AD. We carried out several sensitivity ensemble experiments to separate the effects of (i) reduced solar ultra-violet (UV) irradiance, (ii) reduced solar visible and near infrared irradiance, (iii) enhanced galactic cosmic ray intensity as well as less intensive solar energetic proton events and auroral electron precipitation, and (iv) volcanic aerosols. The introduced changes of UV irradiance and volcanic aerosols significantly influence stratospheric dynamics in the early 19th century, whereas changes in the visible part of the spectrum and energetic particles have smaller effects. A reduction of UV irradiance by 15%, which represents the presently discussed highest estimate of UV irradiance change caused by solar activity changes, causes global ozone decrease below the stratopause reaching as much as 8% in the midlatitudes at 5 hPa and a significant stratospheric cooling of up to 2 °C in the mid-stratosphere and to 6 °C in the lower mesosphere. Changes in energetic particle precipitation lead only to minor changes in the yearly averaged temperature fields in the stratosphere. Volcanic aerosols heat the tropical lower stratosphere, allowing more water vapour to enter the tropical stratosphere, which, via HOx reactions, decreases upper stratospheric and mesospheric ozone by roughly 4%. Conversely, heterogeneous chemistry on aerosols reduces stratospheric NOx, leading to a 12% ozone increase in the tropics, whereas a decrease in ozone of up to 5% is found over Antarctica in boreal winter. The linear superposition of the different contributions is not equivalent to the response obtained in a simulation when all forcing factors are applied during the Dalton Minimum (DM) – this effect is especially well visible for NOx/NOy. Thus, this study also shows the non-linear behaviour of the coupled chemistry-climate system. Finally, we conclude that especially UV and volcanic eruptions dominate the changes in the ozone, temperature and dynamics while the NOx field is dominated by the energetic particle precipitation. Visible radiation changes have only very minor effects on both stratospheric dynamics and chemistry.

Coupled W–Os–Pt isotope systematics in IVB iron meteorites: In situ neutron dosimetry for W isotope chronology

Tungsten isotope compositions of magmatic iron meteorites yield ages of differentiation that are within ±2 Ma of the formation of CAIs, with the exception of IVB irons that plot to systematically less radiogenic compositions yielding erroneously old ages. Secondary neutron capture due to galactic cosmic ray (GCR) irradiation is known to lower the ε182W of iron meteorites, adequate correction of which requires a measure of neutron dosage which has not been available, thus far. The W, Os and Pt isotope systematics of 12 of the 13 known IVB iron meteorites were determined by MC-ICP-MS (W, Os, Pt) and TIMS (Os). On the same dissolutions that yield precise ε182W, stable Os and Pt isotopes were determined as in situ neutron dosimeters for empirical correction of the ubiquitous cosmic-ray induced burn-out of 182W in iron meteorites. The W isotope data reveal a main cluster with ε182W of ∼−3.6, but a much larger range than observed in previous studies including irons (Weaver Mountains and Warburton Range) that show essentially no cosmogenic effect on their ε182W. The IVB data exhibits resolvable negative anomalies in ε189Os (−0.6ε) and complementary ε190Os anomalies (+0.4ε) in Tlacotepec due to neutron capture on 189Os which has approximately the same neutron capture cross section as 182W, and captures neutrons to produce 190Os. The least irradiated IVB iron, Warburton Range, has ε189Os and ε190Os identical to terrestrial values. Similarly, Pt isotopes, which are presented as ε192Pt, ε194Pt and ε196Pt range from +4.4ε to +53ε, +1.54ε to −0.32ε and +0.73ε to −0.20ε, respectively, also identify Tlacotepec and Dumont as the most GCR-damaged samples. In W–Os and W–Pt isotope space, the correlated isotope data back-project toward a 0-epsilon value of ε192Pt, ε189Os and ε190Os from which a pre-GCR irradiation ε182W of −3.42±0.09 (2σ) is derived. This pre-GCR irradiation ε182W is within uncertainty of the currently accepted CAI initial ε182W. The Pt and Os isotope correlations in the IVB irons are in good agreement with a nuclear model for spherical irons undergoing GCR spallation, although this model over-predicts the change of ε182W by ∼2×, indicating a need for better W neutron capture cross section determinations. A nucleosynthetic effect in ε184W in these irons of −0.14±0.08 is confirmed, consistent with the presence of Mo and Ru isotope anomalies in IVB irons. The lack of a non-GCR Os isotope anomaly in these irons requires more complex explanations for the production of W, Ru and Mo anomalies than nebular heterogeneity in the distribution of s-process to r-process nuclides.

Neutron capture on Pt isotopes in iron meteorites and the Hf–W chronology of core formation in planetesimals

The short-lived 182Hf–182W isotope system can provide powerful constraints on the timescales of planetary core formation, but its application to iron meteorites is hampered by neutron capture reactions on W isotopes resulting from exposure to galactic cosmic rays. Here we show that Pt isotopes in magmatic iron meteorites are also affected by capture of (epi)thermal neutrons and that the Pt isotope variations are correlated with variations in 182W/184W. This makes Pt isotopes a sensitive neutron dosimeter for correcting cosmic ray-induced W isotope shifts. The pre-exposure 182W/184W derived from the Pt–W isotope correlations of the IID, IVA and IVB iron meteorites are higher than most previous estimates and are more radiogenic than the initial 182W/184W of Ca–Al-rich inclusions (CAI). The Hf–W model ages for core formation range from +1.6±1.0 million years (Ma; for the IVA irons) to +2.7±1.3 Ma after CAI formation (for the IID irons), indicating that there was a time gap of at least ∼1 Ma between CAI formation and metal segregation in the parent bodies of some iron meteorites. From the Hf–W ages a time limit of <1.5–2 Ma after CAI formation can be inferred for the accretion of the IID, IVA and IVB iron meteorite parent bodies, consistent with earlier conclusions that the accretion of differentiated planetesimals predated that of most chondrite parent bodies.