Astrodynamical Space Test of Relativity using Optical Devices I (ASTROD I)—a class-M fundamental physics mission proposal for cosmic vision 2015–2025: 2010 Update

ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system; and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals.For this mission, accurate pulse timing with an ultra-stable clock, and a drag-free spacecraft with reliable inertial sensor are required. T2L2 has demonstrated the required accurate pulse timing; rubidium clock on board Galileo has mostly demonstrated the required clock stability; the accelerometer on board GOCE has paved the way for achieving the reliable inertial sensor; the demonstration of LISA Pathfinder will provide an excellent platform for the implementation of the ASTROD I drag-free spacecraft. These European activities comprise the pillars for building up the mission and make the technologies needed ready. A second mission, ASTROD or ASTROD-GW (depending on the results of ASTROD I), is envisaged as a three-spacecraft mission which, in the case of ASTROD, would test General Relativity to one part per billion, enable detection of solar g-modes, measure the solar Lense-Thirring effect to 10 parts per million, and probe gravitational waves at frequencies below the LISA bandwidth, or in the case of ASTROD-GW, would be dedicated to probe gravitational waves at frequencies below the LISA bandwidth to 100?nHz and to detect solar g-mode oscillations. In the third phase (Super-ASTROD), larger orbits could be implemented to map the outer solar system and to probe primordial gravitational-waves at frequencies below the ASTROD bandwidth. This paper on ASTROD I is based on our 2010 proposal submitted for the ESA call for class-M mission proposals, and is a sequel and an update to our previous paper (Appouchaux et al., Exp Astron 23:491?C527, 2009; designated as Paper I) which was based on our last proposal submitted for the 2007 ESA call. In this paper, we present our orbit selection with one Venus swing-by together with orbit simulation. In Paper I, our orbit choice is with two Venus swing-bys. The present choice takes shorter time (about 250?days) to reach the opposite side of the Sun. We also present a preliminary design of the optical bench, and elaborate on the solar physics goals with the radiation monitor payload. We discuss telescope size, trade-offs of drag-free sensitivities, thermal issues and present an outlook.     

Groundwater recharge rates for regional groundwater modelling: a case study using GROWA in the Lower Rhine lignite mining area, Germany

Groundwater recharge rates calculated with the GROWA model have been applied as the recharge boundary condition for the regional groundwater model Rurscholle. This model simulates groundwater dynamics in the Pleistocene aquifers of the Lower Rhine lignite mining area (Germany). GROWA uses an area-differentiated approach to calculate recharge rates depending on runoff-relevant site characteristics, which are represented by a set of baseflow indices. The regional accuracy of the coupled groundwater and GROWA models has been checked using groundwater hydrographs as validation criteria. The results suggest that the current (unadjusted) version of GROWA underestimates the regional groundwater recharge rate by 10–20 mm/yr. The comparative analysis identified areas where recharge calculations could be improved by adjusting the baseflow indices for areas where runoff is dominated by slope, low water-logging and a low degree of sealing. Using the adjusted set of baseflow indices, the mean groundwater recharge rate of the Rurscholle region was modelled as approx. 170 mm/yr. This study highlights the benefit of using a coupled approach and being able to independently calibrate and validate groundwater recharge boundary conditions in regional groundwater models.     

Noble gases and nitrogen in Martian meteorites Dar al Gani 476, Sayh al Uhaymir 005 and Lewis Cliff 88516: EFA and extra neon

Meteorite “finds” from the terrestrial hot deserts have become a major contributor to the inventory of Martian meteorites. In order to understand their nitrogen and noble gas components, we have carried out stepped heating experiments on samples from two Martian meteorites collected from hot deserts. We measured interior and surface bulk samples, glassy and non-glassy portions of Dar al Gani 476 and Sayh al Uhaymir 005. We have also analyzed noble gases released from the Antarctic shergottite Lewis Cliff 88516 by crushing and stepped heating. For the hot desert meteorites significant terrestrial Ar, Kr, Xe contamination is observed, with an elementally fractionated air (EFA) component dominating the low temperature releases. The extremely low Ar/Kr/Xe ratios of EFA may be the result of multiple episodes of trapping/loss during terrestrial alteration involving aqueous fluids. We suggest fractionation processes similar to those in hot deserts to have acted on Mars, with acidic weathering on the latter possibly even more effective in producing elementally fractionated components. Addition from fission xenon is apparent in DaG 476 and SaU 005. The Ar-Kr-Xe patterns for LEW 88516 show trends as typically observed in shergottites - including evidence for a crush-released component similar to that observed in EETA 79001. A trapped Ne component most prominent in the surface sample of DaG 476 may represent air contamination. It is accompanied by little trapped Ar (²⁰Ne/³⁶Ar > 50) and literature data suggest its presence also in some Antarctic finds. Data for LEW 88516 and literature data, on the other hand, suggest the presence of two trapped Ne components of Martian origin characterized by different ²⁰Ne/²²Ne, possibly related to the atmosphere and the interior. Caution is recommended in interpreting nitrogen and noble gas isotopic signatures of Martian meteorites from hot deserts in terms of extraterrestrial sources and processes. Nevertheless our results provide hope that vice-versa, via noble gases and nitrogen in meteorites and other relevant samples from terrestrial deserts, Martian secondary processes can be studied.     

On the asteroid belt's orbital and size distribution

For absolute magnitudes greater than the current completeness limit of H-magnitude ∼15 the main asteroid belt's size distribution is imperfectly known. We have acquired good-quality orbital and absolute H-magnitude determinations for a sample of small main-belt asteroids in order to study the orbital and size distribution beyond H=15, down to sub-kilometer sizes (H>18). Based on six observing nights over a 11-night baseline we have detected, measured photometry for, and linked observations of 1087 asteroids which have one-week time baselines or more. The linkages allow the computation of full heliocentric orbits (as opposed to statistical distances determined by some past surveys). Judged by known asteroids in the field the typical uncertainty in the (a/e/i) orbital elements is less than 0.03 AU/0.03/0.5°. The distances to the objects are sufficiently well known that photometric uncertainties (of 0.3 magnitudes or better) dominate the error budget of their derived H-magnitudes. The detected asteroids range from HR=12-22 and provide a set of objects down to sizes below 1 km in diameter. We find an on-sky surface density of 210 asteroids per square degree in the ecliptic with opposition magnitudes brighter than mR=23, with the cumulative number of asteroids increasing by a factor of 10^0.27/mag from mR=18 down to the mR?23.5 limit of our survey. In terms of absolute H magnitudes, we find that beyond H=15 the belt exhibits a constant power-law slope with the number increasing proportional to ¹⁰0.30H from H?15 to 18, after which incompleteness begins in the survey. Examining only the subset of detections inside 2.5 AU, we find weak evidence for a mildly shallower slope for H=15-19.5. We provide the information necessary such that anyone wishing to model the main asteroid belt can compare a detailed model to our detected sample.     

Injection of Oort Cloud comets: the fundamental role of stellar perturbations

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1% for a < 20 000 AU to about 100% for a ? 100 000 AU.     

Development of a Matrix‐Matched Sphalerite Reference Material (MUL‐ZnS‐1) for Calibration of In Situ Trace Element Measurements by Laser Ablation‐Inductively Coupled Plasma‐Mass Spectrometry

Sphalerite (ZnS) is an abundant ore mineral and an important carrier of elements such as Ge, Ga and In used in high‐technology applications. In situ measurements of trace elements in natural sphalerite samples using LA‐ICP‐MS are hampered by a lack of homogenous matrix‐matched sulfide reference materials available for calibration. The preparation of the MUL‐ZnS1 calibration material containing the trace elements V, Cr, Mn, Co, Ni, Cu, Ga, Ge, As, Se, Mo, Ag, Cd, In, Sn, Sb, Tl and Pb besides Zn, Fe and S is reported. Commercially available ZnS, FeS, CdS products were used as the major components, whereas the trace elements were added by doping with single‐element ICP‐MS standard solutions and natural mineral powders. The resulting powder mixture was pressed to pellets and sintered at 400 °C for 100 h using argon as an inert gas. To confirm the homogeneity of major and trace element distributions within the MUL‐ZnS1 calibration material, measurements were performed using EPMA, solution ICP‐MS, ICP‐OES and LA‐ICP‐MS. The results show that MUL‐ZnS‐1 is an appropriate material for calibrating trace element determination in sphalerite using LA‐ICP‐MS.     

Dynamic instabilities under isotropic drained compression of idealized granular materials

The compressibility behaviour of loose and contracting granular assemblies, normally consolidated and overconsolidated, under isotropic drained compression is investigated in this paper. Short cylindrical samples of water-saturated monodisperse glass beads, initially assembled in loose state by moist-tamping technique, are isotropically compressed in a classical axisymmetric triaxial machine. Very loose glass bead samples experience numerous unexpected events, sometimes cascading, under undetermined triggered effective isotropic stress in loading and in unloading, while the classical compressibility behaviour of granular material is recovered once these events ignored. Each event, resembling the stick–slip instability during shear in triaxial compression, is characterized by a transient dynamic phase I with very fast drop of effective isotropic stress \(\sigma ^{'}\) due to an excess pore pressure development at nearly constant volume and constant axial strain, followed by a quasi-static phase II with gradual increase in axial \(\varepsilon _\mathrm{a}\) (contraction) and volumetric \(\varepsilon _\mathrm{v}\) (compaction) strain, and a full progressive recovery of \(\sigma ^{'}\) to the previous level before event. A short-lived liquefaction with null \(\sigma ^{'}\) measured in the first phase I results in a local collapse state. Collapse events also happen on unsaturated moist and dry states. Rare events even occur during the unloading of subsequent isotropic compression cycles. The effects of triggered isotropic stress are discussed, the instability characteristics measured, the comparison with stick–slip instability made and the hypothesis of micro-structural instability with local collapse of contact networks and rapid micro-structural rearrangement argued.     

Electrical conductivity of NaCl-bearing aqueous fluids to 600 °C and 1 GPa

The electrical conductivity of aqueous fluids containing 0.01, 0.1, and 1 M NaCl was measured in an externally heated diamond cell to 600 °C and 1 GPa. These measurements therefore more than double the pressure range of previous data and extend it to higher NaCl concentrations relevant for crustal and mantle fluids. Electrical conductivity was generally found to increase with pressure and fluid salinity. The conductivity increase observed upon variation of NaCl concentration from 0.1 to 1 M was smaller than from 0.01 to 0.1 M, which reflects the reduced degree of dissociation at high NaCl concentration. Measured conductivities can be reproduced (R ² = 0.96) by a numerical model with log \(\sigma\) = ?1.7060– 93.78/T + 0.8075 log c + 3.0781 log \(\rho\) + log \(\varLambda\) ₀(T, \(\rho\)), where \(\sigma\) is the conductivity in S m^?1, T is temperature in K, c is NaCl concentration in wt%, \(\rho\) is the density of pure water (in g/cm³) at given pressure and temperature, and \(\varLambda\) ₀ (T, \(\rho\)) is the molar conductivity of NaCl in water at infinite dilution (in S cm² mol^?1), \(\varLambda\) ₀ = 1573–1212 \(\rho\) + 537 062/T–208 122 721/T ². This model allows accurate predictions of the conductivity of saline fluids throughout most of the crust and upper mantle; it should not be used at temperatures below 100 °C. In general, the data show that already a very small fraction of NaCl-bearing aqueous fluid in the deep crust is sufficient to enhance bulk conductivities to values that would be expected for a high degree of partial melting. Accordingly, aqueous fluids may be distinguished from hydrous melts by comparing magnetotelluric and seismic data. H₂O–NaCl fluids may enhance electrical conductivities in the deep crust with little disturbance of v _p or v _p/v _s ratios. However, at the high temperatures in the mantle wedge above subduction zones, the conductivity of hydrous basaltic melts and saline aqueous fluids is rather similar, so that distinguishing these two phases from conductivity data alone is difficult. Observed conductivities in forearc regions, where temperatures are too low to allow melting, may be accounted for by not more than 1 wt% of an aqueous fluid with 5 wt% NaCl, if this fluid forms a continuous film or fills interconnected tubes.