A submillimetre survey of the kinematics of the Perseus molecular cloud – I. Data

We present submillimetre observations of the   J = 3 → 2  rotational transition of ¹²CO, ¹³CO and C¹⁸O across over 600 arcmin² of the Perseus molecular cloud, undertaken with the Heterodyne Array Receiver Programme (HARP), a new array spectrograph on the James Clerk Maxwell Telescope. The data encompass four regions of the cloud, containing the largest clusters of dust continuum condensations: NGC 1333, IC348, L1448 and L1455. A new procedure to remove striping artefacts from the raw HARP data is introduced. We compare the maps to those of the dust continuum emission mapped with the Submillimetre Common-User Bolometer Array (SCUBA; Hatchell et al.) and the positions of starless and protostellar cores (Hatchell et al.). No straightforward correlation is found between the masses of each region derived from the HARP CO and SCUBA data, underlining the care that must be exercised when comparing masses of the same object derived from different tracers. From the ¹³CO/C¹⁸O line ratio the relative abundance of the two species  ([¹³CO]/[C¹⁸O]∼ 7)  and their opacities (typically τ is 0.02–0.22 and 0.15–1.52 for the C¹⁸O and ¹³CO gas, respectively) are calculated. C¹⁸O is optically thin nearly everywhere, increasing in opacity towards star-forming cores but not beyond  τ₁₈∼ 0.9  . Assuming the ¹²CO gas is optically thick, we compute its excitation temperature, T _ex (around 8–30 K), which has little correlation with estimates of the dust temperature.     

Assimilation of thermal emission spectrometer atmospheric data during the Mars Global Surveyor aerobraking period

The Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft has produced an extensive atmospheric data set, beginning during aerobraking and continuing throughout the extended scientific mapping phase. Temperature profiles for the atmosphere below about 40 km, surface temperatures and total dust and water ice opacities, can be retrieved from infrared spectra in nadir viewing mode. This paper describes assimilation of nadir retrievals from the spacecraft aerobraking period, L_S=190°–260°, northern hemisphere autumn to winter, into a Mars general circulation model. The assimilation scheme is able to combine information from temperature and dust optical depth retrievals, making use of a model forecast containing information from the assimilation of earlier observations, to obtain a global, time-dependent analysis. Given sufficient temperature retrievals, the assimilation procedure indicates errors in the a priori dust distribution assumptions even when lacking dust observations; in this case there are relatively cold regions above the poles compared to a model which assumes a horizontally-uniform dust distribution. One major reason for using assimilation techniques is in order to investigate the transient wave behavior on Mars. Whilst the data from the 2-h spacecraft mapping orbit phase is much more suitable for assimilation, even the longer (45–24 h) period aerobraking orbit data contain useful information about the three-dimensional synoptic-scale martian circulation which the assimilation procedure can reconstruct in a consistent way. Assimilations from the period of the Noachis regional dust storm demonstrate that the combined assimilation of temperature and dust retrievals has a beneficial impact on the atmospheric analysis.     

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.     

Hyperspectral mapping of alteration assemblages within a hydrothermal vug at the Haughton impact structure,Canada

Meteorite impacts on Earth and Mars can generate hydrothermal systems that alter the primary mineralogies of rocks and provide suitable environments for microbial colonization. We investigate a calcite–marcasite‐bearing vug at the ~23 km diameter Haughton impact structure, Devon Island, Nunavut, Canada, using imaging spectroscopy of the outcrop in the field (0.65–1.1 μm) and samples in the laboratory (0.4–2.5 μm), point spectroscopy (0.35–2.5 μm), major element chemistry, and X‐ray diffraction analyses. The mineral assemblages mapped at the outcrop include marcasite; marcasite with minor gypsum and jarosite; fibroferrite and copiapite with minor gypsum and melanterite; gypsum, Fe³⁺ oxides, and jarosite; and calcite, gypsum, clay, microcline, and quartz. Hyperspectral mapping of alteration phases shows spatial patterns that illuminate changes in alteration conditions and formation of specific mineral phases. Marcasite formed from the postimpact hydrothermal system under reducing conditions, while subsequent weathering oxidized the marcasite at low temperatures and water/rock ratios. The acidic fluids resulting from the oxidation collected on flat‐lying portions of the outcrop, precipitating fibroferrite + copiapite. That assemblage then likely dissolved, and the changing chemistry and pH resulting from interaction with the calcite‐rich host rock formed gypsum‐bearing red coatings. These results have implications for understanding water–rock interactions and habitabilities at this site and on Mars.     

Book Review

The aim of this book is to present techniques for the study of motion of solar system objects in highly eccentric orbits. Instead of using the usual anomalies (mean, true, eccentric), the authors define and use a new kind of anomaly, the elliptic anomaly.In this way, it is possible, in a theory using perturbation series expansions, to make the ratio: (accuracy)/(number of needed terms), higher than in the classical techniques. The book consists of six chapters. The first chapter deals with the elliptic anomaly in the two-body problem and the second chapter presents the general technique to construct first-order perturbation theory in elliptic function expansions. The next three chapters deal with applications of the new technique to artificial satellites and asteroids, in highly eccentric orbits. The last chapter describes the basic algorithms of the theory.The tools developed in the book demand the use of computer algebra, which is implemented by means of Mathematica 3.0.The book is well written and the new technique is clearly presented and related to the existing techniques, making it useful to all those who use analytical or semi-analytical methods for the study of highly eccentric motion. Celestial Mechanics at High Eccentricities, Gordon and Breach Publishers, US$95, GBP 59, EUR 79, ISBN 90-5699-212-0     

Hydrodynamic simulations of propagating warps and bending waves in accretion discs

We present the results of a study of propagating warp or bending waves in accretion discs. Three-dimensional hydrodynamic simulations were performed using smoothed particle hydrodynamics (SPH), and the results are compared with calculations based on the linear theory of warped discs.
We examine the response of a gaseous disc to an initially imposed warping disturbance under a variety of physical conditions. We consider primarily the physical regime in which the dimensionless viscosity parameter α < H r , where H r is the disc aspect ratio, so that bending waves are expected to propagate. We also performed calculations for disc models in which α > H r , where the warps are expected to evolve diffusively. Small-amplitude (linear) perturbations are studied in both Keplerian and slightly non-Keplerian discs, and we find that the results of the SPH calculations can be reasonably well fitted by those of the linear theory. The main results of these calculations are: (i) the warp in Keplerian discs when α < H r propagates with little dispersion, and damps at a rate expected from estimates of the code viscosity; (ii) warps evolve diffusively when α > H r ; (iii) the slightly non-Keplerian discs lead to a substantially more dispersive behaviour of the warps, which damp at a similar rate to the Keplerian case, when α < H r .
Initially imposed higher amplitude, non-linear warping disturbances were studied in Keplerian discs. The results indicate that non-linear warps can lead to the formation of shocks, and that the evolution of the warp becomes less wave-like and more diffusive in character.
This work is relevant to the study of the warped accretion discs that may occur around Kerr black holes or in misaligned binary systems, and is mainly concerned with discs in which α < H r . The results indicate that SPH can model the hydrodynamics of warped discs, even when using rather modest numbers of particles.     

Silicic differentiates of abyssal oceanic magmas: Evidence for late-magmatic vapor transport of potassium

Massive, nearly holocrystalline dolerites from DSDP Hole 417D contain from 0.5 to 1.5% of granophyric patches composed mainly of Na-plagioclase and quartz. These patches are compositionally similar to other crystalline silicic rocks from oceanic spreading centers and differ from rarer abyssal silicic glasses. Crystalline varieties withSiO₂60wt.% generally haveNa/K>10, whereas silicic glasses have Na/K in the range 3–6. While crystal fractionation readily accounts for the Na₂O and K₂O contents of abyssal silicic glasses, both the 417D granophyres and other crystalline abyssal silicic rocks have much lower K₂O than that predicted by any reasonable crystal-liquid fractionation model. We propose that high-temperature vapor phase transport is responsible for removal of potassium during late-stage crystallization of these rocks. This allows for the formation of cogenetic silicic glassy and crystalline rocks with greatly different Na/K ratios. These observations and interpretations lead to a more confident assignment of high Na/K silicic rocks of oceanic and ophiolitic environments to a cogenetic origin with basaltic oceanic crust.     

Correction: Guidelines and pitfalls of refraction microtremor surveys

Journal of Seismology -