Measurement of the electron neutrino charged-current interaction rate on water with the T2K ND280 π0 detector

This paper presents a measurement of the charged current interaction rate of the electron neutrino beam component of the beam above 1.5 GeV using the large fiducial mass of the T2K π0 detector. The predominant portion of the νe flux (∼85%) at these energies comes from kaon decays. The measured ratio of the observed beam interaction rate to the predicted rate in the detector with water targets filled is 0.89 ± 0.08 (stat.) ± 0.11 (sys.), and with the water targets emptied is 0.90 ± 0.09 (stat.) ± 0.13 (sys.). The ratio obtained for the interactions on water only from an event subtraction method is 0.87 ± 0.33 (stat.) ± 0.21 (sys.). This is the first measurement of the interaction rate of electron neutrinos on water, which is particularly of interest to experiments with water Cherenkov detectors.

Net accumulation rates derived from ice core stable isotope records of Pío XI glacier, Southern Patagonia Icefield

Pío XI, the largest glacier of the Southern Patagonia Icefield, reached its neoglacial maximum extent in 1994 and is one of the few glaciers in that area which is not retreating. In view of the recent warming it is important to understand glacier responses to climate changes. Due to its remoteness and the harsh conditions in Patagonia, no systematic mass balance studies have been performed. In this study we derived net accumulation rates for the period 2000–2006 from a 50 m (33.2 4 m weq) ice core collected in the accumulation area of Pío XI (2600 m a.s.l., 49°16'40"S, 73°21'14"W). Borehole temperatures indicate near temperate ice, but the average melt percent is only 16 ± 14%. Records of stable isotopes are well preserved and were used for identification of annual layers. Net accumulation rates range from 3.4–7.1 water equivalent (m weq) with an average of 5.8 m weq, comparable to precipitation amounts at the Chilean coast, but not as high as expected for the Icefield. Ice core stable isotope data correlate well with upper air temperatures and may be used as temperature proxy.

Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

TALDICE-1 age scale of the Talos Dome deep ice core, East Antarctica

A new deep ice core drilling program, TALDICE, has been successfully handled by a European team at Talos Dome, in the Ross Sea sector of East Antarctica, down to 1620 m depth. Using stratigraphic markers and a new inverse method, we produce the first official chronology of the ice core, called TALDICE-1. We show that it notably improves an a priori chronology resulting from a one-dimensional ice flow model. It is in agreement with a posteriori controls of the resulting accumulation rate and thinning function along the core. An absolute uncertainty of only 300 yr is obtained over the course of the last deglaciation. This uncertainty remains lower than 600 yr over Marine Isotope Stage 3, back to 50 kyr BP. The phasing of the TALDICE ice core climate record with respect to the central East Antarctic plateau and Greenland records can thus be determined with a precision allowing for a discussion of the mechanisms at work at sub-millennial time scales.

Glacier mass balance reconstruction by sublimation induced enrichment of chemical species on Cerro Tapado (Chilean Andes)

A 36 m long ice core down to bedrock from the Cerro Tapado glacier (5536 m a.s.l, 30°08' S, 69°55' W) was analyzed to reconstruct past climatic conditions for Northern Chile. Because of the marked seasonality in the precipitation (short wet winter and extended dry summer periods) in this region, major snow ablation and related post-depositional processes occur on the glacier surface during summer periods. They include predominantly sublimation and dry deposition. Assuming that, like measured during the field campaign, the enrichment of chloride was always related to sublimation, the chemical record along the ice core may be applied to reconstruct the history of such secondary processes linked to the past climatic conditions over northern Chile. For the time period 1962–1999, a mean annual net accumulation of 316 mm water equivalent (weq) and 327 mm weq loss by sublimation was deduced by this method. This corresponds to an initial total annual accumulation of 539 mm weq. The annual variability of the accumulation and sublimation is related with the Southern Oscillation Index (SOI): higher net-accumulation during El-Niño years and more sublimation during La Niña years. The deepest part of the ice record shows a time discontinuity; with an ice body deposited under different climatic conditions: 290 mm higher precipitation but with reduced seasonal distribution (+470 mm in winter and –180 mm in summer) and –3°C lower mean annual temperature. Unfortunately, its age is unknown. The comparison with regional proxy data however let us conclude that the glacier buildup did most likely occur after the dry mid-Holocene.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

Temperature and precipitation signal in two Alpine ice cores over the period 1961-2001

Water stable isotope ratios and net snow accumulation in ice cores are commonly interpreted as temperature or precipitation proxies. However, only in a few cases has a direct calibration with instrumental data been attempted. In this study we took advantage of the dense network of observations in the European Alpine region to rigorously test the relationship of the annual and seasonal resolved proxy data from two highly resolved ice cores with local temperature and precipitation. We focused on the time period 1961–2001 with the highest amount and quality of meteorological data and the minimal uncertainty in ice core dating (±1 year). The two ice cores were retrieved from the Fiescherhorn glacier (northern Alps, 3900 m a.s.l.), and Grenzgletscher (southern Alps, 4200 m a.s.l.). A parallel core from the Fiescherhorn glacier allowed assessing the reproducibility of the ice core proxy data. Due to the orographic barrier, the two flanks of the Alpine chain are affected by distinct patterns of precipitation. The different location of the two glaciers therefore offers a unique opportunity to test whether such a specific setting is reflected in the proxy data. On a seasonal scale a high fraction of δ18O variability was explained by the seasonal cycle of temperature (~60% for the ice cores, ~70% for the nearby stations of the Global Network of Isotopes in Precipitation – GNIP). When the seasonality is removed, the correlations decrease for all sites, indicating that factors other than temperature such as changing moisture sources and/or precipitation regimes affect the isotopic signal on this timescale. Post-depositional phenomena may additionally modify the ice core data. On an annual scale, the δ18O/temperature relationship was significant at the Fiescherhorn, whereas for Grenzgletscher this was the case only when weighting the temperature with precipitation. In both cases the fraction of interannual temperature variability explained was ~20%, comparable to the values obtained from the GNIP stations data. Consistently with previous studies, we found an altitude effect for the δ18O of −0.17‰/100 m for an extended elevation range combining data of the two ice core sites and four GNIP stations. Significant correlations between net accumulation and precipitation were observed for Grenzgletscher during the entire period of investigation, whereas for Fiescherhorn this was the case only for the less recent period (1961–1977). Local phenomena, probably related to wind, seem to partly disturb the Fiescherhorn accumulation record. Spatial correlation analysis shows the two glaciers to be influenced by different precipitation regimes, with the Grenzgletscher reflecting the characteristic precipitation regime south of the Alps and the Fiescherhorn accumulation showing a pattern more closely linked to northern Alpine stations.

Martian sub-surface ionising radiation: biosignature and geology

The surface of Mars, unshielded by thick atmosphere or global magnetic field, is exposed to high levels of cosmic radiation. This ionising radiation field is deleterious to the survival of dormant cells or spores and the persistence of molecular biomarkers in the subsurface, and so its characterisation is of prime astrobiological interest. Here, we present modelling results of the absorbed radiation dose as a function of depth through the Martian subsurface, suitable for calculation of biomarker persistence. A second major implementation of this dose accumulation rate data is in application of the optically stimulated luminescence technique for dating Martian sediments. We present calculations of the dose-depth profile in the Martian subsurface for various scenarios: variations of surface composition (dry regolith, ice, layered permafrost), solar minimum and maximum conditions, locations of different elevation (Olympus Mons, Hellas basin, datum altitude), and increasing atmospheric thickness over geological history. We also model the changing composition of the subsurface radiation field with depth compared between Martian locations with different shielding material, determine the relative dose contributions from primaries of different energies, and discuss particle deflection by the crustal magnetic fields.

How stable are 20th century calibration models? A high-resolution summer temperature reconstruction for the eastern Swiss Alps back to A.D. 1580 derived from proglacial varved sediments

We found a significant positive correlation between local summer air temperature (May-September) and the annual sediment mass accumulation rate (MAR) in Lake Silvaplana (46°N, 9°E, 1800 m a.s.l.) during the twentieth century (r = 0.69, p < 0.001 for decadal smoothed series). Sediment trap data (2001-2005) confirm this relation with exceptionally high particle yields during the hottest summer of the last 140 years in 2003. On this base we developed a decadal-scale summer temperature reconstruction back to AD 1580. Surprisingly, the comparison of our reconstruction with two other independent regional summer temperature reconstructions (based on tree-rings and documentary data) revealed a significant negative correlation for the pre-1900 data (ie, late ‘Little Ice Age’). This demonstrates that the correlation between MAR and summer temperature is not stable in time and the actualistic principle does not apply in this case. We suggest that different climatic regimes (modern/‘Little Ice Age’) lead to changing state conditions in the catchment and thus to considerably different sediment transport mechanisms. Therefore, we calibrated our MAR data with gridded early instrumental temperature series from AD 1760-1880 (r = -0.48, p < 0.01 for decadal smoothed series) to properly reconstruct the late LIA climatic conditions. We found exceptionally low temperatures between AD 1580 and 1610 (0.75°C below twentieth-century mean) and during the late Maunder Minimum from AD 1680 to 1710 (0.5°C below twentieth-century mean). In general, summer temperatures did not experience major negative departures from the twentieth-century mean during the late ‘Little Ice Age’. This compares well with the two existing independent regional reconstructions suggesting that the LIA in the Alps was mainly a phenomenon of the cold season.

Dosimetric properties of an amorphous silicon EPID for verification of modulated electron radiotherapy

Purpose: To investigate the dosimetric properties of an electronic portal imaging device (EPID) for electron beam detection and to evaluate its potential for quality assurance (QA) of modulated electron radiotherapy (MERT). Methods: A commercially available EPID was used to detect electron beams shaped by a photon multileaf collimator (MLC) at a source-surface distance of 70 cm. The fundamental dosimetric properties such as reproducibility, dose linearity, field size response, energy response, and saturation were investigated for electron beams. A new method to acquire the flood-field for the EPID calibration was tested. For validation purpose, profiles of open fields and various MLC fields (square and irregular) were measured with a diode in water and compared to the EPID measurements. Finally, in order to use the EPID for QA of MERT delivery, a method was developed to reconstruct EPID two-dimensional (2D) dose distributions in a water-equivalent depth of 1.5 cm. Comparisons were performed with film measurement for static and dynamic monoenergy fields as well as for multienergy fields composed by several segments of different electron energies. Results: The advantageous EPID dosimetric properties already known for photons as reproducibility, linearity with dose, and dose rate were found to be identical for electron detection. The flood-field calibration method was proven to be effective and the EPID was capable to accurately reproduce the dose measured in water at 1.0 cm depth for 6 MeV, 1.3 cm for 9 MeV, and 1.5 cm for 12, 15, and 18 MeV. The deviations between the output factors measured with EPID and in water at these depths were within ±1.2% for all the energies with a mean deviation of 0.1%. The average gamma pass rate (criteria: 1.5%, 1.5 mm) for profile comparison between EPID and measurements in water was better than 99% for all the energies considered in this study. When comparing the reconstructed EPID 2D dose distributions at 1.5 cm depth to film measurements, the gamma pass rate (criteria: 2%, 2 mm) was better than 97% for all the tested cases. Conclusions: This study demonstrates the high potential of the EPID for electron dosimetry, and in particular, confirms the possibility to use it as an efficient verification tool for MERT delivery.