Enceladus: Present internal structure and differentiation by early and long-term radiogenic heating

Pre-Cassini images of Saturn's small icy moon Enceladus provided the first indication that this satellite has undergone extensive resurfacing and tectonism. Data returned by the Cassini spacecraft have proven Enceladus to be one of the most geologically dynamic bodies in the Solar System. Given that the diameter of Enceladus is only about 500 km, this is a surprising discovery and has made Enceladus an object of much interest. Determining Enceladus' interior structure is key to understanding its current activity. Here we use the mean density of Enceladus (as determined by the Cassini mission to Saturn), Cassini observations of endogenic activity on Enceladus, and numerical simulations of Enceladus' thermal evolution to infer that this satellite is most likely a differentiated body with a large rock-metal core of radius about 150 to 170 km surrounded by a liquid water-ice shell. With a silicate mass fraction of 50% or more, long-term radiogenic heating alone might melt most of the ice in a homogeneous Enceladus after about 500 Myr assuming an initial accretion temperature of about 200 K, no subsolidus convection of the ice, and either a surface temperature higher than at present or a porous, insulating surface. Short-lived radioactivity, e.g., the decay of ²⁶Al, would melt all of the ice and differentiate Enceladus within a few million years of accretion assuming formation of Enceladus at a propitious time prior to the decay of ²⁶Al. Long-lived radioactivity facilitates tidal heating as a source of energy for differentiation by warming the ice in Enceladus so that tidal deformation can become effective. This could explain the difference between Enceladus and Mimas. Mimas, with only a small rock fraction, has experienced relatively little long-term radiogenic heating; it has remained cold and stiff and less susceptible to tidal heating despite its proximity to Saturn and larger eccentricity than Enceladus. It is shown that the shape of Enceladus is not that of a body in hydrostatic equilibrium at its present orbital location and rotation rate. The present shape could be an equilibrium shape corresponding to a time when Enceladus was closer to Saturn and spinning more rapidly, or more likely, to a time when Enceladus was spinning more rapidly at its present orbital location. A liquid water layer on Enceladus is a possible source for the plume in the south polar region assuming the survivability of such a layer to the present. These results could place Enceladus in a category similar to the large satellites of Jupiter, with the core having a rock-metal composition similar to Io, and with a deep overlying ice shell similar to Europa and Ganymede. Indeed, the moment of inertia factor of a differentiated Enceladus, C/MR², could be as small as that of Ganymede, about 0.31.     

Triaxial induction tool with electrode sleeve: FD modeling in 3D geometries

We have developed a finite-difference (FD) method to model the response of a multiple-spacing triaxial array induction tool with a multiple-electrode type sleeve. The FD software was verified versus two other independent modeling methods. The model response was checked for the FD grid refinement, while attempting to minimize run time. Many details of the electrode-sleeve geometry were studied using this method that led to the final electrode-sleeve design. The electrode sleeve reduces the response of the transverse couplings to eccentering in a conductive borehole. Quantitative agreement was found between the measurements and calculations. Using the FD software we study the tool response to various effects of 3D geometry.     

Solar Magnetic Carpet I: Simulation of Synthetic Magnetograms

This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve through the processes of emergence, cancellation, coalescence and fragmentation. Model parameters for the emergence of bipoles are based upon the results of observational studies. Using this model, several simulations are considered, where the range of flux with which bipoles may emerge is varied. In all cases the model quickly reaches a steady state where the rates of emergence and cancellation balance. Analysis of the resulting magnetic field shows that we reproduce observed quantities such as the flux distribution, mean field, cancellation rates, photospheric recycle time and a magnetic network. As expected, the simulation matches observations more closely when a larger, and consequently more realistic, range of emerging flux values is allowed (4×10¹⁶ – 10¹⁹ Mx). The model best reproduces the current observed properties of the magnetic carpet when we take the minimum absolute flux for emerging bipoles to be 4×10¹⁶ Mx. In future, this 2D model will be used as an evolving photospheric boundary condition for 3D non-potential modeling.     

Carbon dioxide on the satellites of Saturn: Results from the Cassini VIMS investigation and revisions to the VIMS wavelength scale

Several of the icy satellites of Saturn show the spectroscopic signature of the asymmetric stretching mode of C-O in carbon dioxide (CO₂) at or near the nominal solid-phase laboratory wavelength of 4.2675 μm (2343.3 cm⁻¹), discovered with the Visible-Infrared Mapping Spectrometer (VIMS) on the Cassini spacecraft. We report here on an analysis of the variation in wavelength and width of the CO₂ absorption band in the spectra of Phoebe, Iapetus, Hyperion, and Dione. Comparisons are made to laboratory spectra of pure CO₂, CO₂ clathrates, ternary mixtures of CO₂ with other volatiles, implanted and adsorbed CO₂ in non-volatile materials, and ab initio theoretical calculations of CO₂ * nH₂O. At the wavelength resolution of VIMS, the CO₂ on Phoebe is indistinguishable from pure CO₂ ice (each molecule’s nearby neighbors are also CO₂) or type II clathrate of CO₂ in H₂O. In contrast, the CO₂ band on Iapetus, Hyperion, and Dione is shifted to shorter wavelengths (typically ∼4.255 μm (∼2350.2 cm⁻¹)) and broadened. These wavelengths are characteristic of complexes of CO₂ with different near-neighbor molecules that are encountered in other volatile mixtures such as with H₂O and CH₃OH, and non-volatile host materials like silicates, some clays, and zeolites. We suggest that Phoebe’s CO₂ is native to the body as part of the initial inventory of condensates and now exposed on the surface, while CO₂ on the other three satellites results at least in part from particle or UV irradiation of native H₂O plus a source of C, implantation or accretion from external sources, or redistribution of native CO₂ from the interior.The analysis presented here depends on an accurate VIMS wavelength scale. In preparation for this work, the baseline wavelength calibration for the Cassini VIMS was found to be distorted around 4.3 μm, apparently as a consequence of telluric CO₂ gas absorption in the pre-launch calibration. The effect can be reproduced by convolving a sequence of model detector response profiles with a deep atmospheric CO₂ absorption profile, producing distorted detector profile shapes and shifted central positions. In a laboratory blackbody spectrum used for radiance calibration, close examination of the CO₂ absorption profile shows a similar deviation from that expected from a model. These modeled effects appear to be sufficient to explain the distortion in the existing wavelength calibration now in use. A modification to the wavelength calibration for 13 adjacent bands is provided. The affected channels span about 0.2 μm centered on 4.28 μm. The maximum wavelength change is about 10 nm toward longer wavelength. This adjustment has implications for interpretation of some of the spectral features observed in the affected wavelength interval, such as from CO₂, as discussed in this paper.     

Sulfur chemistry in the Venus mesosphere from SO2 and SO microwave spectra

First measurements of SO₂ and SO in the Venus mesosphere (70-100 km) are reported. This altitude range is distinctly above the ∼60-70 km range to which nadir-sounding IR and UV investigations are sensitive. Since July 2004, use of ground-based sub-mm spectroscopy has yielded multiple discoveries. Abundance of each molecule varies strongly on many timescales over the entire sub-Earth Venus hemisphere. Diurnal behavior is evident, with more SO₂, and less SO, at night than during the day. Non-diurnal variability is also present, with measured SO₂ and SO abundances each changing by up to 2× or more between observations conducted on different dates, but at fixed phase, hence identical sub-Earth Venus local times. Change as large and rapid as a 5σ doubling of SO on a one-week timescale is seen. The sum of SO₂ and SO abundances varies by an order of magnitude or more, indicating at least one additional sulfur reservoir must be present, and that it must function as both a sink and source for these molecules. The ratio SO₂/SO varies by nearly two orders of magnitude, with both diurnal and non-diurnal components. In contrast to the strong time dependence of molecular abundances, their altitude distributions are temporally invariant, with far more SO₂ and SO at 85-100 km than at 70-85 km. The observed increase of SO₂ mixing ratio with altitude requires that the primary SO₂ source be upper mesospheric photochemistry, contrary to atmospheric models which assert upward transport as the only source of above-cloud SO₂. Abundance of upper mesospheric aerosol, with assumption that it is composed primarily of sulfuric acid, is at least sufficient to provide the maximum gas phase (SO + SO₂) sulfur reported in this study. Sulfate aerosol is thus a plausible source of upper mesospheric SO₂.     

The geochemistry of primary and weathered oil shale and coquina across the Julia Creek vanadium deposit (Queensland, Australia)

A significant resource of vanadium and molybdenum exists near Julia Creek, Australia, where the middle Cretaceous organic-rich Toolebuc Formation lies between 0 and 25 m of the surface. We present and discuss a comprehensive geochemical study of the Toolebuc Formation and its enclosing stratigraphy near Julia Creek to understand this ore deposit. V and Mo contents in fresh facies are strongly associated with total organic carbon (TOC) contents, but not with Al or CaCO₃; this suggests that V and Mo were originally concentrated in the organic fraction. However, chemical extractions using H₂O₂ indicate that Mo was originally concentrated in pyrite. The data also suggest that V was mobilised from organic matter during early diagenesis and became associated with clays as little V was extracted by H₂O₂ in the fresh samples. TOC contents in the Toolebuc Formation were removed during weathering, residually enriching trace metals including V and Mo, and as a result, the TOC relationship with V and Mo disintegrates. With weathering, both V and Mo predominantly became associated with iron oxide/hydroxide phases (and possibly other unidentified phases) as these elements in the weathered facies were highly soluble in the sodium citrate–sodium dithionite digestion. Large shale-hosted V and Mo deposits such as Julia Creek offer a potentially viable alternative to the currently mined magnetite-hosted deposits. A thorough understanding of the formation and host mineral phases for V and Mo of these shale deposits, however, is critical to ensure that these valuable metals can be feasibly extracted.     

Erosion rates on different timescales derived from cosmogenic 10Be and river loads: implications for landscape evolution in the Rhenish Massif, Germany

We determined erosion rates on timescales of 10¹–10⁴ years for two catchments in the northeastern Rhenish Massif, in order to unravel the Quaternary landscape evolution in a Variscan mountain range typical of central Europe. Spatially averaged erosion rates derived from in situ produced ¹⁰Be concentrations in stream sediment of the Aabach and M?hne watersheds range from 47 ± 6 to 65 ± 14 mm/ka and integrate over the last 9–13 ka. These erosion rates are similar to local rates of river incision and rock uplift in the Quaternary and to average denudation rates since the Mesozoic derived from fission track data. This suggests that rock uplift is balanced by denudation, i.e., the landscape is in a steady state. Short-term erosion rates were derived from suspended and dissolved river loads subsequent to (1) correcting for atmospheric and anthropogenic inputs, (2) establishing calibration curves that relate the amount of suspended load to discharge, and (3) estimating the amount of bedload. The resulting solid mass fluxes (suspended and bedload) agree with those derived from the sediment volume trapped in three reservoirs. However, resulting geogenic short-term erosion rates range from 9 to 25 mm/ka and are only about one-third of the rates derived from ¹⁰Be. Model simulations in combination with published sediment yield data suggest that this discrepancy is caused by at least three factors: (1) phases with higher precipitation and/or lower evapotranspiration, (2) rare flood events not captured in the short-term records, and (3) prolonged periods of climatic deterioration with increased erosion and sediment transport on hillslopes.     

A roadmap for a dedicated Earth Science Grid platform

Due to its intensive data processing and highly distributed organization, the multidisciplinary Earth Science applications community is uniquely positioned for the uptake and exploitation of Grid technologies. Currently Enabling Grids for E-sciencE, and other large Grid infrastructures are already deployed and capable of operational services. So far however, the adoption and exploitation of Grid technology throughout the Earth Science community has been slower than expected. The Dissemination and Exploitation of GRids in Earth sciencE project, proposed by the European Commission to assist and accelerate this process in a number of different ways, had between its main goals the creation of a roadmap towards Earth Science Grid platform. This paper presents the resulting roadmap.     

Clod Size and Moisture Condition Influence on the Shearing Behavior of Compacted Soil

This paper presents results of the influence of clod size and initial moisture condition on the shearing behavior of a clayey soil with a plasticity index of 22. The clods were divided into different size categories, and then two groups of samples were prepared; samples with large clods (LC = clod size less than 38 mm) and samples with small clods (SC = clod size less than 4.75 mm). Two initial moisture conditions were achieved using two different moisture tempering times of 0 and 14 days. All samples were compacted to the same density and moisture content (95 % of maximum dry density and 2 % dry of optimum moisture content). Triaxial compression tests, including constant water content unsaturated tests and backpressure saturated, isotropically consolidated, undrained compression tests were performed at different confining pressures. At 0-day tempering time, samples prepared using large clods were found to be stiffer and stronger than those prepared using small clods. However, at 14-day tempering time, the strength of both LC and SC samples was similar.     

Implementing quality control on a random number stream to improve a stochastic weather generator

For decades, stochastic modellers have used computerized random number generators to produce random numeric sequences fitting a specified statistical distribution. Unfortunately, none of the random number generators we tested satisfactorily produced the target distribution. The result is generated distributions whose mean even diverges from the mean used to generate them, regardless of the length of run. Non‐uniform distributions from short sequences of random numbers are a major problem in stochastic climate generation, because truly uniform distributions are required to produce the intended climate parameter distributions. In order to ensure generation of a representative climate with the stochastic weather generator CLIGEN within a 30‐year run, we tested the climate output resulting from various random number generators. The resulting distributions of climate parameters showed significant departures from the target distributions in all cases. We traced this failure back to the uniform random number generators themselves. This paper proposes a quality control approach to select only those numbers that conform to the expected distribution being retained for subsequent use. The approach is based on goodness‐of‐fit analysis applied to the random numbers generated. Normally distributed deviates are further tested with confidence interval tests on their means and standard deviations. The positive effect of the new approach on the climate characteristics generated and the subsequent deterministic process‐based hydrology and soil erosion modelling are illustrated for four climatologically diverse sites. Copyright © 2007 John Wiley & Sons, Ltd.