Origin and thermal evolution of icy satellites

The paper reviews the problem of formation and evolution of the so-called regular satellites of the giant planets, and it consists of two parts: the first describes the possible origin of the satellites, the second studies their evolution, attempting to stress the relations of the present status of the satellites with their evolutionary history.The formation of regular satellite systems around giant planets is probably related to the formation of the central planet. Some characteristics of regular satellite systems are quite similar, and suggest a common origin in a disk present around the central body. This disk can originate through different mechanisms which we will describe, paying attention to the so-called accretion disk model, in which the satellite-forming material is captured. The disk phase links the formation of the primary body with the formation of satellites. The subsequent stages of the disk's evolution can lead first to the formation of intermediate size bodies, and through the collisional evolution of these bodies, to the birth of satellite embryos able to gravitationally capture smaller bodies.Given the scenario in which icy satellites may be formed by homogeneous accretion of planetesimals made of a mixtures of ice and silicates, if no melting occurs during accretion, the satellites have a homogeneous ice-rock composition. For the smaller satellites this homogeneous structure should not be substantially modified; only sporadic local events, such as large impacts, can modify the surface structure of the smaller satellites. For the larger satellites, if some degree of melting appears during accretion, a differentiation of the silicate part occurs, the amount of differentiation and hence the core size depending on the fraction of gravitational potential energy retained during the accumulation process. Melting and differentiation soon after the accretion, for the larger satellites, could also depend on the convective evolution in presence of phase transitions and generate an intermediate rock layer, considerably denser than the underlying, still homogeneous core, and unstable to overturning on a geologic time scale. Moreover the liquid water mantle could be a transient feature because the mantle would freeze over several hundred million years. For these large bodies the stable configuration is expected to be one consisting of a silicate core and a mantle of mixed rock and ice.     

Critical appraisal on the identification of Reference Conditions for the evaluation of ecological quality status along the Emilia-Romagna coast (Italy) using M-AMBI

According to the European Water Framework Directive, the ecological status (ES) of a water body is determined by comparing observation data with undisturbed Reference Conditions (RCs). Defining RCs is crucial when evaluating the ES of a water body as it strongly affects the final outcome of any index application. Identifying RCs by observing real sites is not feasible in many marine environments, such as the Emilia-Romagna coast (Italy, N-Adriatic Sea). We used a statistical approach on a large dataset to derive RCs for the application of the benthic index M-AMBI in this area. We then applied M-AMBI to samples collected along a gradient of presumed environmental disturbance. The results put 14.8% of the Emilia-Romagna samples in “High” ES, 60.2% in “Good”, 23.0% in “Moderate” and 2.0% in “Poor”, showing a spatial gradient of improving quality. These results are in agreement with the extensive ecological knowledge available for this area.     

Experimental multipole-refined and theoretical charge density study of LiGaSi<Subscript>2</Subscript>O<Subscript>6</Subscript> clinopyroxene at ambient conditions

The synthetic LiGaSi₂O₆ clinopyroxene is monoclinic C2/c at room-T. Its experimental electron density, ρ(r), has been derived starting from accurate room-T single-crystal diffraction data. Topological analysis confirms an intermediate ionic-covalent character for Si–O bonding, as found by previous electron-density studies on other silicates such as diopside, coesite and stishovite. The non-bridging Si–O bonds have more covalent character than the bridging ones. The Ga–O bonds have different bonding characters, the Ga–O2 bond being more covalent than the two Ga–O1 bonds. Li–O bonds are classified as pure closed-shell ionic interactions. Similar to spodumene (LiAlSi₂O₆), Li has sixfold coordination, but the bond critical points associated to the two longest bonds are characterized by very low electron density values. Similar to what previously found in spodumene and diopside, O···O interactions were detected from the topological analysis of ρ(r), and indicate a cooperative interaction among the lone pairs of neighbouring oxygen atoms. In particular, this kind of interaction has been obtained for the O1···O1 edge shared between two Ga octahedra. Integration over the atomic basins gives net charges of −1.39(10), 2.82(10), 1.91(10) and 0.82(8) e for O (averaged), Si, Ga and Li atoms, respectively. Periodic Hartree–Fock and DFT calculations confirm the results obtained by multipole refinement of the experimental data. Moreover, the theoretical topological properties of the electron density distribution on the Si₂O₆ group are very similar to those calculated for spodumene. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Askival: An altered feldspathic cumulate sample in Gale crater

Askival is a light-toned, coarsely crystalline float rock, which was identified near the base of Vera Rubin Ridge in Gale crater. We have studied Askival, principally with the ChemCam instrument but also using APXS compositional data and MAHLI images. Askival and an earlier identified sample, Bindi, represent two rare examples of feldspathic cumulate float rocks in Gale crater with >65% relict plagioclase. Bindi appears unaltered whereas Askival shows textural and compositional signatures of silicification, along with alkali remobilization and hydration. Askival likely experienced multiple stages of alteration, occurring first through acidic hydrolysis of metal cations, followed by deposition of silica and possible phyllosilicates at low T and neutral-alkaline pH. Through laser-induced breakdown spectroscopy compositional analyses and normative calculations, we suggest that an assemblage of Fe-Mg silicates including amphibole and pyroxene, Fe phases, and possibly Mg-rich phyllosilicate are present. Thermodynamic modeling of the more pristine Bindi composition predicts that amphibole and feldspar are stable within an upper crustal setting. This is consistent with the presence of amphibole in the parent igneous rocks of Askival and suggests that the paucity of amphiboles in other known Martian samples reflects the lack of representative samples of the Martian crust rather than their absence on Mars.     

Chondritic models of 4 Vesta: Implications for geochemical and geophysical properties

Simple mass‐balance and thermodynamic constraints are used to illustrate the potential geochemical and geophysical diversity of a fully differentiated Vesta‐sized parent body with a eucrite crust (e.g., core size and density, crustal thickness). The results of this analysis are then combined with data from the howardite–eucrite–diogenite (HED) meteorites and the Dawn mission to constrain Vesta's bulk composition. Twelve chondritic compositions are considered, comprising seven carbonaceous, three ordinary, and two enstatite chondrite groups. Our analysis excludes CI and LL compositions as plausible Vesta analogs, as these are predicted to have a negative metal fraction. Second, the MELTS thermodynamic calculator is used to show that the enstatite chondrites, the CV, CK and L‐groups cannot produce Juvinas‐like liquids, and that even for the other groups, depletion in sodium is necessary to produce liquids of appropriate silica content. This conclusion is consistent with the documented volatile‐poor nature of eucrites. Furthermore, carbonaceous chondrites are predicted to have a mantle too rich in olivine to produce typical howardites and to have Fe/Mn ratios generally well in excess of those of the HEDs. On the other hand, an Na‐depleted H‐chondrite bulk composition is capable of producing Juvinas‐like liquids, has a mantle rich enough in pyroxene to produce abundant howardite/diogenite, and has a Fe/Mn ratio compatible with eucrites. In addition, its predicted bulk‐silicate density is within 100 kg m^?3 of solutions constrained by data of the Dawn mission. However, oxidation state and oxygen isotopes are not perfectly reproduced and it is deduced that bulk Vesta may contain approximately 25% of a CM‐like component. Values for the bulk‐silicate composition of Vesta and a preliminary phase diagram are proposed.     

Distribution of iron on Vesta

We have completed a mapping study of 7.6 MeV gamma rays produced by neutron capture by Fe at the surface of the main belt asteroid 4 Vesta as measured by the bismuth germanate scintillator of the Gamma Ray and Neutron Detector (GRaND) on the Dawn spacecraft. The procedures used to determine Fe counting rates are presented, along with a global map, constituting the necessary initial step to quantify Fe abundances. While the final calibration of orbital data to absolute concentrations has not been determined, the range of fully corrected Fe counting rates is compared with that of Fe in howardites. We find that the global distribution of corrected Fe counting rates is generally consistent with mineralogy and composition determined independently by other instruments on the Dawn spacecraft, including measurements of pyroxene absorption bands by the Visible and Infrared Spectrometer and Framing Camera, and an index of diogenitic materials provided by neutron absorption measurements by GRaND. In addition, there is a distinctive low Fe region in the western hemisphere that was not reported by reflectance or optical observations, possibly indicating the presence of a cumulate eucrite component in Vesta's regolith.     

Ganymedean pedestal craters distribution: Implications on thermal and tectonic histories

The distribution analysis of the ganymedean pedestal craters shows a very good correlation between them and the grooved terrains. These craters seem to be the earliest post-grooved impacts. The presence of tectonized pedestal craters is also noted. This type of crater is the only one lying on the grooved terrains and being affected by the grooves. Assuming that a lowering of the terrain's viscosity is the determining factor for the presence of pedestal craters, we conclude that at the time of their formation, the grooved terrains have a lower viscosity than the other terrains. Moreover, using the density ratio between non-pedestal craters and pedestal craters, a relative low viscosity state's duration time is calculated. Based on the density ratio between pedestal craters and tectonized pedestal craters, a grooved terrains formation's time is also calculated. These two times are of the order of 10⁷ yr and 10⁶ yr respectively. This period of low viscosity may be partly due to internal heating.     

U‐Th Zircon Dating by Laser Ablation Single Collector Inductively Coupled Plasma‐Mass Spectrometry (LA‐ICP‐MS)

Zircon crystals in the age range of ca. 10–300 ka can be dated by ²³⁰Th/²³⁸U (U‐Th) disequilibrium methods because of the strong fractionation between Th and U during crystallisation of zircon from melts. Laser ablation inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) analysis of nine commonly used reference zircons (at secular equilibrium) and a synthetic zircon indicates that corrections for abundance sensitivity and dizirconium trioxide molecular ions (Zr₂O₃⁺) are critical for reliable determination of ²³⁰Th abundances in zircon. When corrected for abundance sensitivity and interferences, mean activity ratios of (²³⁰Th)/(²³⁸U) for nine reference zircons analysed on five different days averaged 0.995 ± 0.023 (95% confidence weighted by data‐point uncertainty only, MSWD = 1.6; n = 9), consistent with their U‐Pb ages > 4 Ma that imply equilibrium for all intermediate daughter isotopes (including ²³⁰Th) within the ²³⁸U decay chain. U‐Th zircon ages generated by LA‐ICP‐MS without mitigating (e.g., by high mass resolution) or correcting for abundance sensitivity and molecular interferences on ²³⁰Th are potentially unreliable. To validate the applicability of LA‐ICP‐MS to this dating method, we acquired data from three late Quaternary volcanic units: the 41 ka Campanian Ignimbrite (plutonic clasts), the 161 ka Kos Plateau Tuff (juvenile clasts) and the 12 ka Puy de Dôme trachyte lava (all eruption ages by Ar/Ar, with zircon U‐Th ages being of equal or slightly older). A comparison of the corrected LA‐ICP‐MS results with previously published secondary ion mass spectrometry (SIMS) data for these rocks shows comparable ages with equivalent precision for LA‐ICP‐MS and SIMS, but much shorter analysis durations (~ 2 min vs. ~ 15 min) per spot with LA‐ICP‐MS and much simpler sample preparation. Previously undated zircons from the Yali eruption (Kos‐Nisyros volcanic centre, Greece) were analysed using this method. This yielded a large age spread (~ 45 to > 300 ka), suggesting significant antecryst recycling. The youngest zircon age (~ 45 ± 10 ka) provides a reasonable maximum estimate for the eruption age, in agreement with the previously published age using oxygen isotope stratigraphy (~ 31 ka).     

A study of the properties of a local dust storm with Mars Express OMEGA and PFS data

We present observations of a local dust storm performed by the OMEGA and PFS instruments aboard Mars Express. OMEGA observations are used to retrieve the dust single-scattering albedo in the spectral range 0.4-4.0 μm. The single-scattering albedo shows fairly constant values between 0.6 and 2.6 μm, and a sharp decrease at wavelengths shorter than 0.6 μm, in agreement with previous studies. It presents a small absorption feature due to ferric oxide at 0.9 μm, and a strong absorption feature due to hydrated minerals between 2.7 and 3.6 μm. We use a statistical method, the Independent Component Analysis, to determine that the dust spectral signature is decoupled from the surface albedo, proving that the retrieval of the single-scattering albedo is reliable, and we map the dust optical thickness with a conventional radiative transfer model. The effect of the dust storm on the atmospheric thermal structure is measured using PFS observations. We also simulate the thermal impact of the dust storm using a one-dimensional atmospheric model. A comparison of the retrieved and modeled temperature structures suggests that the dust in the storm should be confined to the 1-2 lowest scale heights of the atmosphere. However, the observed OMEGA reflectance in the CO₂ absorption bands does not support this suggestion.