Detecting offsets in GPS time series: first results from the Detection of Offsets in GPS Experiment

The accuracy of Global Positioning System (GPS) time series is degraded by the presence of offsets. To assess the effectiveness of methods that detect and remove these offsets, we designed and managed the Detection of Offsets in GPS Experiment. We simulated time series that mimicked realistic GPS data consisting of a velocity component, offsets, white and flicker noises (1/f spectrum noises) composed in an additive model. The data set was made available to the GPS analysis community without revealing the offsets, and several groups conducted blind tests with a range of detection approaches. The results show that, at present, manual methods (where offsets are hand picked) almost always give better results than automated or semi‒automated methods (two automated methods give quite similar velocity bias as the best manual solutions). For instance, the fifth percentile range (5% to 95%) in velocity bias for automated approaches is equal to 4.2 mm/year (most commonly ±0.4 mm/yr from the truth), whereas it is equal to 1.8 mm/yr for the manual solutions (most commonly 0.2 mm/yr from the truth). The magnitude of offsets detectable by manual solutions is smaller than for automated solutions, with the smallest detectable offset for the best manual and automatic solutions equal to 5 mm and 8 mm, respectively. Assuming the simulated time series noise levels are representative of real GPS time series, robust geophysical interpretation of individual site velocities lower than 0.2–0.4 mm/yr is therefore certainly not robust, although a limit of nearer 1 mm/yr would be a more conservative choice. Further work to improve offset detection in GPS coordinates time series is required before we can routinely interpret sub‒mm/yr velocities for single GPS stations.

Detection of nonrandom association of alleles from the distribution of the number of heterozygous loci in a sample.

The distribution of the number of heterozygous loci in two randomly chosen gametes or in a random diploid zygote provides information regarding the nonrandom association of alleles among different genetic loci. Two alternative statistics may be employed for detection of nonrandom association of genes of different loci when observations are made on these distributions: observed variance of the number of heterozygous loci (s2k) and a goodness-of-fit criterion (X2) to contrast the observed distribution with that expected under the hypothesis of random association of genes. It is shown, by simulation, that s2k is statistically more efficient than X2 to detect a given extent of nonrandom association. Asymptotic normality of s2k is justified, and X2 is shown to follow a chi-square (chi 2) distribution with partial loss of degrees of freedom arising because of estimation of parameters from the marginal gene frequency data. Whenever direct evaluations of linkage disequilibrium values are possible, tests based on maximum likelihood estimators of linkage disequilibria require a smaller sample size (number of zygotes or gametes) to detect a given level of nonrandom association in comparison with that required if such tests are conducted on the basis of s2k. Summarization of multilocus genotype (or haplotype) data, into the different number of heterozygous loci classes, thus, amounts to appreciable loss of information.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

Light noble gases in 12 meteorites from the Omani desert, Australia, Mauritania, Canada, and Sweden

We measured the concentrations and isotopic compositions of He, Ne, and Ar in 14 fragments from 12 different meteorites: three carbonaceous chondrites, six L chondrites (three most likely paired), one H chondrite, one R chondrite, and one ungrouped chondrite. The data obtained for the CV3 chondrites Ramlat as Sahmah (RaS) 221 and RaS 251 support the hypothesis of exposure age peaks for CV chondrites at approximately 9 Ma and 27 Ma. The exposure age for Shişr 033 (CR chondrite) of 7.3 Ma is also indicative of a possible CR chondrite exposure age peak. The three L chondrites Jiddat al Harasis (JaH) 091, JaH 230, and JaH 296, which are most likely paired, fall together with Hallingeberg into the L chondrite exposure age peak of approximately 15 Ma. The two L chondrites Shelburne and Lake Torrens fall into the peaks at approximately 40 Ma and 5 Ma, respectively. The ages for Bassikounou (H chondrite) and RaS 201 (R chondrite) are approximately 3.5 Ma and 5.8 Ma, respectively. Six of the studied meteorites show clear evidence for 3He diffusive losses, the deficits range from approximately 17% for one Lake Torrens aliquot to approximately 45% for RaS 211. The three carbonaceous chondrites RaS 221, RaS 251, and Shişr 033 all have excess 4He, either of planetary or solar origin. However, very high 4He/20Ne ratios occur at relatively low 20Ne/22Ne ratios, which is unexpected and needs further study. The measured 40Ar ages fit well into established systematics. They are between 2.5 and 4.5 Ga for the carbonaceous chondrites, older than 3.6 Ga for the L and H chondrites, and about 2.4 Ga for the R chondrite as well as for the ungrouped chondrite. Interestingly, none of our studied L chondrites has been degassed in the 470 Ma break-up event. Using the amount of trapped 36Ar as a proxy for noble gas contamination due to terrestrial weathering we are able to demonstrate that the samples studied here are not or only very slightly affected by terrestrial weathering (at least in terms of their noble gas budget).

Nuclear power plants and childhood leukaemia: lessons from the past and future directions

In the 1980s, leukaemia clusters were discovered around nuclear fuel reprocessing plants in Sellafield and Dounreay in the United Kingdom. This raised public concern about the risk of childhood leukaemia near nuclear power plants (NPPs). Since then, the topic has been well-studied, but methodological limitations make results difficult to interpret. Our review aims to: (1.) summarise current evidence on the relationship between NPPs and risk of childhood leukaemia, with a focus on the Swiss CANUPIS (Childhood cancer and nuclear power plants in Switzerland) study; (2.) discuss the limitations of previous research; and (3.) suggest directions for future research. There are various reasons that previous studies produced inconclusive results. These include: inadequate study designs and limited statistical power due to the low prevalence of exposure (living near a NPP) and outcome (leukaemia); lack of accurate exposure estimates; limited knowledge of the aetiology of childhood leukaemia, particularly of vulnerable time windows and latent periods; use of residential location at time of diagnosis only and lack of data on address histories; and inability to adjust for potential confounders. We conclude that risk of childhood leukaemia around NPPs should continue to be monitored and that study designs should be improved and standardised. Data should be pooled internationally to increase the statistical power. More research needs to be done on other putative risk factors for childhood cancer such as low-dose ionizing radiation, exposure to certain chemicals and exposure to infections. Studies should be designed to allow examining multiple exposures.

Microwave radiometer to retrieve temperature profiles from the surface to the stratopause

TEMPERA (TEMPERature RAdiometer) is a new ground-based radiometer which measures in a frequency range from 51–57 GHz radiation emitted by the atmosphere. With this instrument it is possible to measure temperature profiles from ground to about 50 km. This is the first ground-based instrument with the capability to retrieve temperature profiles simultaneously for the troposphere and stratosphere. The measurement is done with a filterbank in combination with a digital fast Fourier transform spectrometer. A hot load and a noise diode are used as stable calibration sources. The optics consist of an off-axis parabolic mirror to collect the sky radiation. Due to the Zeeman effect on the emission lines used, the maximum height for the temperature retrieval is about 50 km. The effect is apparent in the measured spectra. The performance of TEMPERA is validated by comparison with nearby radiosonde and satellite data from the Microwave Limb Sounder on the Aura satellite. In this paper we present the design and measurement method of the instrument followed by a description of the retrieval method, together with a validation of TEMPERA data over its first year, 2012.

Daily ozone cycle in the stratosphere: global, regional and seasonal behaviour modelled with the Whole Atmosphere Community Climate Model

The Whole Atmosphere Community Climate Model (WACCM) is utilised to study the daily ozone cycle and underlying photochemical and dynamical processes. The analysis is focused on the daily ozone cycle in the middle stratosphere at 5 hPa where satellite-based trend estimates of stratospheric ozone are most biased by diurnal sampling effects and drifting satellite orbits. The simulated ozone cycle shows a minimum after sunrise and a maximum in the late afternoon. Further, a seasonal variation of the daily ozone cycle in the stratosphere was found. Depending on season and latitude, the peak-to-valley difference of the daily ozone cycle varies mostly between 3 and 5% (0.4 ppmv) with respect to the midnight ozone volume mixing ratio. The maximal variation of 15% (0.8 ppmv) is found at the polar circle in summer. The global pattern of the strength of the daily ozone cycle is mainly governed by the solar zenith angle and the sunshine duration. In addition, we find synoptic-scale variations in the strength of the daily ozone cycle. These variations are often anti-correlated to regional temperature anomalies and are due to the temperature dependence of the rate coefficients k2 and k3 of the Chapman cycle reactions. Further, the NOx catalytic cycle counteracts the accumulation of ozone during daytime and leads to an anti-correlation between anomalies in NOx and the strength of the daily ozone cycle. Similarly, ozone recombines with atomic oxygen which leads to an anti-correlation between anomalies in ozone abundance and the strength of the daily ozone cycle. At higher latitudes, an increase of the westerly (easterly) wind cause a decrease (increase) in the sunshine duration of an air parcel leading to a weaker (stronger) daily ozone cycle.

A 13,600-year diatom oxygen isotope record from the South Carpathians (Romania): Reflection of winter conditions and possible links with North Atlantic circulation changes

This study provides a continuous lateglacial and Holocene record of diatom silica oxygen isotope changes (delta O-18(DIAT)) in a subalpine lake sediment sequence obtained from the Retezat Mts (Taul dintre Brazi, 1740 m a.s.l.). This through-flow, shallow, high-altitude lake with a surface area of only 0.4 ha has short water residence time and is predominantly fed by snowmelt and rainwater. Its delta O-18(DIAT) record principally reflects the oxygen isotope composition of the winter and spring precipitation, as diatom blooms occur mainly in the spring and early summer. Hence, changes in delta O-18(DIAT) are interpreted as seasonal scale changes: in the amount of winter precipitation. Low oxygen isotope values (27-28.5 parts per thousand) occurred during the lateglacial until 12,300 cal BP, followed by a sharp increase thereafter. In the Holocene delta O-18(DIAT) values ranged from 29 to 31 parts per thousand until 3200 cal BP, followed by generally lower values during the late Holocene (27-30 parts per thousand). Short-term decreases in the isotopic values were found between 10,140-9570, 9000-8500, 7800-7300, 6300-5800, 5500-5000 and at 8015, 4400, 4000 cal BP. After 3200 cal BP a decreasing trend was visible with the lowest values between 3100-2500 and after 2100 cal BP The general trend in the record suggests that contribution of winter precipitation was generally lower between 11,680 and 3200 cal BP, followed by increased contribution during the last millennia. The late Holocene decrease in delta O-18(DIAT) shows good agreement with the speleothem delta O-18, lake level and testate amoebae records from the Carpathian Mountains that also display gradual delta O-18 decrease and lake level/mire water table level rise after 3200 cal BR Strong positive correlation with North Atlantic circulation and solar activity proxies, such as the Austrian and Hungarian speleothem records, furthermore suggested that short-term increases in the isotopic ratios in the early and mid Holocene are likely connectable to high solar activity phases and high frequency of positive North Atlantic Oscillation indexes that may have resulted in decreased winter precipitation in this region.

Fire in ice: two millennia of boreal forest fire history from the Greenland NEEM ice core

Biomass burning is a major source of greenhouse gases and influences regional to global climate. Pre-industrial fire-history records from black carbon, charcoal and other proxies provide baseline estimates of biomass burning at local to global scales spanning millennia, and are thus useful to examine the role of fire in the carbon cycle and climate system. Here we use the specific biomarker levoglucosan together with black carbon and ammonium concentrations from the North Greenland Eemian (NEEM) ice cores (77.49° N, 51.2° W; 2480 m a.s.l) over the past 2000 years to infer changes in boreal fire activity. Increases in boreal fire activity over the periods 1000–1300 CE and decreases during 700–900 CE coincide with high-latitude NH temperature changes. Levoglucosan concentrations in the NEEM ice cores peak between 1500 and 1700 CE, and most levoglucosan spikes coincide with the most extensive central and northern Asian droughts of the past millennium. Many of these multi-annual droughts are caused by Asian monsoon failures, thus suggesting a connection between low- and high-latitude climate processes. North America is a primary source of biomass burning aerosols due to its relative proximity to the Greenland Ice Cap. During major fire events, however, isotopic analyses of dust, back trajectories and links with levoglucosan peaks and regional drought reconstructions suggest that Siberia is also an important source of pyrogenic aerosols to Greenland.

The diurnal variation in stratospheric ozone from the MACC reanalysis, the ERA-Interim reanalysis, WACCM and Earth observations: characteristics and intercomparison

In this study we compare the diurnal variation in stratospheric ozone derived from free-running simulations of the Whole Atmosphere Community Climate Model (WACCM) and from reanalysis data of the atmospheric service MACC (Monitoring Atmospheric Composition and Climate) which both use a similar stratospheric chemistry module. We find good agreement between WACCM and the MACC reanalysis for the diurnal ozone variation in the high-latitude summer stratosphere based on photochemistry. In addition, we consult the ozone data product of the ERA-Interim reanalysis. The ERA-Interim reanalysis ozone system with its long-term ozone parametrization can not capture these diurnal variations in the upper stratosphere that are due to photochemistry. The good dynamics representations, however, reflects well dynamically induced ozone variations in the lower stratosphere. For the high-latitude winter stratosphere we describe a novel feature of diurnal variation in ozone where changes of up to 46.6% (3.3 ppmv) occur in monthly mean data. For this effect good agreement between the ERA-Interim reanalysis and the MACC reanalysis suggest quite similar diurnal advection processes of ozone. The free-running WACCM model seriously underestimates the role of diurnal advection processes at the polar vortex at the two tested resolutions. The intercomparison of the MACC reanalysis and the ERA-Interim reanalysis demonstrates how global reanalyses can benefit from a chemical representation held by a chemical transport model. The MACC reanalysis provides an unprecedented description of the dynamics and photochemistry of the diurnal variation of stratospheric ozone which is of high interest for ozone trend analysis and research on atmospheric tides. We confirm the diurnal variation in ozone at 5 hPa by observations of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) experiment and selected sites of the Network for Detection of Atmospheric Composition Change (NDACC). The latter give valuable insight even to diurnal variation of ozone in the polar winter stratosphere.

Reducing CT radiation dose with iterative reconstruction algorithms: the influence of scan and reconstruction parameters on image quality and CTDIvol

OBJECTIVES In this phantom CT study, we investigated whether images reconstructed using filtered back projection (FBP) and iterative reconstruction (IR) with reduced tube voltage and current have equivalent quality. We evaluated the effects of different acquisition and reconstruction parameter settings on image quality and radiation doses. Additionally, patient CT studies were evaluated to confirm our phantom results. METHODS Helical and axial 256 multi-slice computed tomography scans of the phantom (Catphan(®)) were performed with varying tube voltages (80-140kV) and currents (30-200mAs). 198 phantom data sets were reconstructed applying FBP and IR with increasing iterations, and soft and sharp kernels. Further, 25 chest and abdomen CT scans, performed with high and low exposure per patient, were reconstructed with IR and FBP. Two independent observers evaluated image quality and radiation doses of both phantom and patient scans. RESULTS In phantom scans, noise reduction was significantly improved using IR with increasing iterations, independent from tissue, scan-mode, tube-voltage, current, and kernel. IR did not affect high-contrast resolution. Low-contrast resolution was also not negatively affected, but improved in scans with doses <5mGy, although object detectability generally decreased with the lowering of exposure. At comparable image quality levels, CTDIvol was reduced by 26-50% using IR. In patients, applying IR vs. FBP resulted in good to excellent image quality, while tube voltage and current settings could be significantly decreased. CONCLUSIONS Our phantom experiments demonstrate that image quality levels of FBP reconstructions can also be achieved at lower tube voltages and tube currents when applying IR. Our findings could be confirmed in patients revealing the potential of IR to significantly reduce CT radiation doses.

Backscattered energetic neutral atoms from the Moon in the Earth's plasma sheet observed by Chandarayaan-1/Sub-keV Atom Reflecting Analyzer instrument

We present the observations of energetic neutral atoms (ENAs) produced at the lunar surface in the Earth's magnetotail. When the Moon was located in the terrestrial plasma sheet, Chandrayaan-1 Energetic Neutrals Analyzer (CENA) detected hydrogen ENAs from the Moon. Analysis of the data from CENA together with the Solar Wind Monitor (SWIM) onboard Chandrayaan-1 reveals the characteristic energy of the observed ENA energy spectrum (the e-folding energy of the distribution function) ∼100 eV and the ENA backscattering ratio (defined as the ratio of upward ENA flux to downward proton flux) <∼0.1. These characteristics are similar to those of the backscattered ENAs in the solar wind, suggesting that CENA detected plasma sheet particles backscattered as ENAs from the lunar surface. The observed ENA backscattering ratio in the plasma sheet exhibits no significant difference in the Southern Hemisphere, where a large and strong magnetized region exists, compared with that in the Northern Hemisphere. This is contrary to the CENA observations in the solar wind, when the backscattering ratio drops by ∼50% in the Southern Hemisphere. Our analysis and test particle simulations suggest that magnetic shielding of the lunar surface in the plasma sheet is less effective than in the solar wind due to the broad velocity distributions of the plasma sheet protons.

Warm Breze from the Starboard Bow: a new Population of neutral Helium in the Heliosphere

We investigate the signals from neutral helium atoms observed in situ from Earth orbit in 2010 by the Interstellar Boundary Explorer (IBEX). The full helium signal observed during the 2010 observation season can be explained as a superposition of pristine neutral interstellar He gas and an additional population of neutral helium that we call the Warm Breeze. The Warm Breeze is approximately 2 times slower and 2.5 times warmer than the primary interstellar He population, and its density in front of the heliosphere is ~7% that of the neutral interstellar helium. The inflow direction of the Warm Breeze differs by ~19° from the inflow direction of interstellar gas. The Warm Breeze seems to be a long-term, perhaps permanent feature of the heliospheric environment. It has not been detected earlier because it is strongly ionized inside the heliosphere. This effect brings it below the threshold of detection via pickup ion and heliospheric backscatter glow observations, as well as by the direct sampling of GAS/Ulysses. We discuss possible sources for the Warm Breeze, including (1) the secondary population of interstellar helium, created via charge exchange and perhaps elastic scattering of neutral interstellar He atoms on interstellar He+ ions in the outer heliosheath, or (2) a gust of interstellar He originating from a hypothetic wave train in the Local Interstellar Cloud. A secondary population is expected from models, but the characteristics of the Warm Breeze do not fully conform to modeling results. If, nevertheless, this is the explanation, IBEX-Lo observations of the Warm Breeze provide key insights into the physical state of plasma in the outer heliosheath. If the second hypothesis is true, the source is likely to be located within a few thousand AU from the Sun, which is the propagation range of possible gusts of interstellar neutral helium with the Warm Breeze characteristics against dissipation via elastic scattering in the Local Cloud. Whatever the nature of the Warm Breeze, its discovery exposes a critical new feature of our heliospheric environment.

First Axion Results from the XENON100 Experiment

We present the first results of searches for axions and axionlike particles with the XENON100 experiment. The axion-electron coupling constant, g Ae , has been probed by exploiting the axioelectric effect in liquid xenon. A profile likelihood analysis of 224.6 live days × 34-kg exposure has shown no evidence for a signal. By rejecting g Ae larger than 7.7×10 −12 (90% C.L.) in the solar axion search, we set the best limit to date on this coupling. In the frame of the DFSZ and KSVZ models, we exclude QCD axions heavier than 0.3 and 80  eV/c 2 , respectively. For axionlike particles, under the assumption that they constitute the whole abundance of dark matter in our galaxy, we constrain g Ae to be lower than 1×10 −12 (90% C.L.) for masses between 5 and 10  keV/c 2 .

OSL-thermochronometry of feldspar from the KTB borehole, Germany

The reconstruction of thermal histories of rocks (thermochronometry) is a fundamental tool both in Earth science and in geological exploration. However, few methods are currently capable of resolving the low-temperature thermal evolution of the upper ∼2 km of the Earth's crust. Here we introduce a new thermochronometer based on the infrared stimulated luminescence (IRSL) from feldspar, and validate the extrapolation of its response to artificial radiation and heat in the laboratory to natural environmental conditions. Specifically, we present a new detailed Na-feldspar IRSL thermochronology from a well-documented thermally-stable crustal environment at the German Continental Deep Drilling Program (KTB). There, the natural luminescence of Na-feldspar extracted from twelve borehole samples (0.1–2.3 km depth, corresponding to 10–70 °C) can be either (i) predicted within uncertainties from the current geothermal gradient, or (ii) inverted into a geothermal palaeogradient of 29 ± 2  °C km−1, integrating natural thermal conditions over the last ∼65 ka. The demonstrated ability to invert a depth–luminescence dataset into a meaningful geothermal palaeogradient opens new venues for reconstructing recent ambient temperatures of the shallow crust (<0.3 Ma, 40–70 °C range), or for studying equally recent and rapid transient cooling in active orogens (<0.3 Ma, >200 °C Ma−1 range). Although Na-feldspar IRSL is prone to field saturation in colder or slower environments, the method's primary relevance appears to be for borehole and tunnel studies, where it may offer remarkably recent (<0.3 Ma) information on the thermal structure and history of hydrothermal fields, nuclear waste repositories and hydrocarbon reservoirs.