Perspectives and challenges for future telescopes

At the celebration of 400 years of telescopic observations, it is appropriate to consider future ground-based facilities enabled by continuing technological developments, and driven by astronomical questions which are amongst the most fundamental in science and are of enormous interest to the general public.     

The DynaMICCS perspective

The DynaMICCS mission is designed to probe and understand the dynamics of crucial regions of the Sun that determine solar variability, including the previously unexplored inner core, the radiative/convective zone interface layers, the photosphere/chromosphere layers and the low corona. The mission delivers data and knowledge that no other known mission provides for understanding space weather and space climate and for advancing stellar physics (internal dynamics) and fundamental physics (neutrino properties, atomic physics, gravitational moments...). The science objectives are achieved using Doppler and magnetic measurements of the solar surface, helioseismic and coronographic measurements, solar irradiance at different wavelengths and in-situ measurements of plasma/energetic particles/magnetic fields. The DynaMICCS payload uses an original concept studied by Thalès Alenia Space in the framework of the CNES call for formation flying missions: an external occultation of the solar light is obtained by putting an occulter spacecraft 150 m (or more) in front of a second spacecraft. The occulter spacecraft, a LEO platform of the mini sat class, e.g. PROTEUS, type carries the helioseismic and irradiance instruments and the formation flying technologies. The latter spacecraft of the same type carries a visible and infrared coronagraph for a unique observation of the solar corona and instrumentation for the study of the solar wind and imagers. This mission must guarantee long (one 11-year solar cycle) and continuous observations (duty cycle > 94%) of signals that can be very weak (the gravity mode detection supposes the measurement of velocity smaller than 1 mm/s). This assumes no interruption in observation and very stable thermal conditions. The preferred orbit therefore is the L1 orbit, which fits these requirements very well and is also an attractive environment for the spacecraft due to its low radiation and low perturbation (solar pressure) environment. This mission is secured by instrumental R and D activities during the present and coming years. Some prototypes of different instruments are already built (GOLFNG, SDM) and the performances will be checked before launch on the ground or in space through planned missions of CNES and PROBA ESA missions (PICARD, LYRA, maybe ASPIICS).     

Strategies for plane change of Earth orbits using lunar gravity and derived trajectories of family G

The dynamics of the circular restricted three-body Earth-Moon-particle problem predicts the existence of the retrograde periodic orbits around the Lagrangian equilibrium point L1. Such orbits belong to the so-called family G (Broucke, Periodic orbits in the restricted three-body problem with Earth-Moon masses, JPL Technical Report 32–1168, 1968) and starting from them it is possible to define a set of trajectories that form round trip links between the Earth and the Moon. These links occur even with more complex dynamical systems as the complete Sun-Earth-Moon-particle problem. One of the most remarkable properties of these trajectories, observed for the four-body problem, is a meaningful inclination gain when they penetrate into the lunar sphere of influence and accomplish a swing-by with the Moon. This way, when one of these trajectories returns to the proximities of the Earth, it will be in a different orbital plane from its initial Earth orbit. In this work, we present studies that show the possibility of using this property mainly to accomplish transfer maneuvers between two Earth orbits with different altitudes and inclinations, with low cost, taking into account the dynamics of the four-body problem and of the swing-by as well. The results show that it is possible to design a set of nominal transfer trajectories that require ΔV _Total less than conventional methods like Hohmann, bi-elliptic and bi-parabolic transfer with plane change.     

Linking the Chassigny meteorite and the Martian surface rock Backstay: Insights into igneous crustal differentiation processes on Mars

Abstract— In order to use igneous surface lithologies to constrain Martian mantle characteristics, secondary processes that lead to compositional modification of primary mantle melts must be considered. Crystal fractionation of a mantle‐derived magma at the base of the crust followed by separation and ascent of residual liquids to the surface is common in continental hotspot regions on Earth. The possibility that this process also takes place on Mars was investigated by experimentally determining whether a surface rock, specifically the hawaiite Backstay analyzed by the MER Spirit could produce a known cumulate lithology with a deep origin (namely the assemblages of the Chassigny meteorite) if trapped at the base of the Martian crust. Both the major cumulus and melt inclusion mineral assemblages of the Chassigny meteorite were produced experimentally by a liquid of Backstay composition within the pressure range 9.3 to 6.8 kbar with bulk water contents between 1.5 and 2.6 wt%. Experiments at 4.3 and 2.8 kbar did not produce the requisite assemblages. This agreement suggests that just as on Earth, Martian mantle‐derived melts may rise to the surface or remain trapped at the base of the crust, fractionate, and lose their residual liquids. Efficient removal of these residual liquids at depth would yield a deep low‐silica cumulate layer for higher magmatic water content; at lower magmatic water content this cumulate layer would be basaltic with shergottitic affinity.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

Rate of growth of the preserved North American continental crust: Evidence from Hf and O isotopes in Mississippi detrital zircons

Detrital zircons from the Mississippi River have been analyzed for U-Th-Pb, Lu-Hf and O isotopes to constrain the rate of growth of the preserved North American continental crust. One hundred and forty two concordant zircon U/Pb dates on grains mounted in epoxy, obtained by Excimer laser ablation ICP-MS method, resolved six major periods of zircon crystallization: 0-0.25, 0.3-0.6, 0.95-1.25, 1.3-1.5, 1.65-1.95 and 2.5-3.0 Ga. These age ranges match the ages of the recognized lithotectonic units of the North American continent in the hinterland of the Mississippi River. Ninety-six zircons mounted on tape, which show no age zonation and were within 7.5% of concordance, were selected to represent the six U/Pb age time intervals and analyzed for Lu-Hf and O isotope by laser ablation MC-ICP-MS and SHRIMP II, respectively. The δ¹⁸O values of the zircons show a small step increase in the maximum δ¹⁸O values at the Archean-Proterozoic boundary from 7.5‰ in the Archean to 9.5‰, and rarely 13‰, in the Proterozoic and Phanerozoic. However, the average value of δ¹⁸O in zircons changes little with time, showing that the increase in the maximum δ¹⁸O values between 2.5 and 2.0 Ga, which can be attributed to an increase in the sediment content of the source regions of younger granitoids, is largely balanced by an increase in zircons with anomalously low δ¹⁸O, which can be attributed to hydrothermally altered crust in the granitoid source region.εHf_i values for the zircons range from 13.1 to −26.9. Zircons derived from juvenile crust, which we define as having mantle δ¹⁸O (4.5-6.5‰) and lying within error of the Hf depleted mantle growth curve, are rare or absent in the Mississippi basin. The overwhelming majority of zircons crystallized from melted pre-existing continental crust, or mantle-derived magmas that were contaminated by continental crust. The average time difference between primitive crust formation and remelting for each of the recognized lithotectonic time intervals, which is defined as crustal incubation time in this study, is 890 ± 460 Myr. There is also a suggestion that the crustal incubation time increases with decreasing age in the Mississippi basin, which is consistent with the declining role of radioactive heat production in the lower crust with time.The average Hf model age (1.94 Ga), weighted by fraction of zircons in the river load is in reasonable agreement with the Nd model age (1.7 Ga) for the Mississippi River. However, if the zircons are weighted by the area of North America covered by the six recognized periods of zircon crystallization the average model age is 2.35 Ga, which compares favorably with an area weighted Nd model age of 2.36 Ga. Our preferred approach is to use the measured O isotope values to constrain variations in the ¹⁷⁶Lu/¹⁷⁷Hf ratio of the granitic source region from which the zircons crystallized, making the assumption that zircons with mantle-like O isotopic ratios have higher ¹⁷⁶Lu/¹⁷⁷Hf than zircons with higher O isotope values. This method gives an average Hf model age of 2.53 Ga, which is 180 Myr older than the constant ¹⁷⁶Lu/¹⁷⁷Hf calculation.The area weighted zircon Hf model ages show two distinct periods of crust formation for the North American continent, 1.6-2.2 and 2.9-3.4 Ga. At least 50% of the preserved North American continental crust was extracted from the mantle by 2.9 Ga and 90% by 1.6 Ga. Two similar periods of crustal growth are also recognized in Gondwana (Hawkesworth C. J. and Kemp A. I. S. (2006) Using hafnium and oxygen isotopes in zircons to unravel the record of crustal evolution. Chem. Geol.226, 144-162.), suggesting that these may be periods of global continental crustal growth. However, we stress that more data from other continents are required before the hypothesis of episodic global continental growth can be accepted with confidence.     

Differential aerobic and anaerobic oxidation of hydrocarbon gases discharged at mud volcanoes in the Nile deep-sea fan

The present study investigates hydrocarbon oxidation processes at Isis and Amon mud volcanoes (MV’s), in the eastern Nile deep-sea fan. In the water column, molecular and carbon isotopic signatures of light hydrocarbons indicate that gases rapidly dissolve in seawater and are partially oxidized.In the upper sediments, anaerobic oxidation of the light hydrocarbons takes place, as clearly shown by their molecular and isotopic composition. These processes lead to the presence of a distinct Sulfate-Hydrocarbon Interface at 120-145 cm and 20-50 cm below the seafloor, for Isis and Amon MV’s, respectively. In contrast to processes occurring in the water column, a clear preferential oxidation of methane, propane and n-butane over ethane and i-butane is observed in the anoxic sediments. Furthermore, for the first time, fractionation factors have been determined for the anaerobic oxidation of propane and butane, being respectively −4.80‰ and −0.7‰ for δ¹³C, and −43.3‰ for δ²H of propane.     

Examining potential genetic links between Jurassic porphyry Cu–Au ± Mo and epithermal Au ± Ag mineralization in the Toodoggone district of North-Central British Columbia,Canada

The Toodoggone district comprises Upper Triassic to Lower Jurassic Hazelton Group Toodoggone Formation volcanic and sedimentary rocks, which unconformably overlie submarine island-arc volcanic and sedimentary rocks of the Lower Permian Asitka Group and Middle Triassic Takla Group, some of which are intruded by Upper Triassic to Lower Jurassic plutons and dikes of the Black Lake suite. Although plutonism occurred episodically from ca. 218 to 191 Ma, the largest porphyry Cu–Au ± Mo systems formed from ca. 202 to 197 Ma, with minor mineralization occurring from ca. 197 to 194 Ma. Porphyry-style mineralization is hosted by small-volume (<1 km³), single-phase, porphyritic igneous stocks or dikes that have high-K calc-alkaline compositions and are comparable with volcanic-arc granites. The Fin porphyry Cu–Au–Mo deposit is anomalous in that it is 16 m.y. older than any other porphyry Cu–Au ± Mo occurrence in the district and has lower REEs. All porphyry systems are spatially restricted to exposed Asitka and Takla Group basement rocks, and rarely, the lowest member of the Hazelton Group (i.e., the ca. 201 Ma Duncan Member). The basement rocks to intrusions are best exposed in the southern half of the district, where high rates of erosion and uplift have resulted in their preferential exposure. In contrast, low- and high-sulfidation epithermal systems are more numerous in the northern half of the district, where the overlying Hazelton Group rocks dominate exposures. Cogenetic porphyry systems might also exist in the northern areas; however, if they are present, they are likely to be buried deeply beneath Hazelton Group rocks. High-sulfidation epithermal systems formed at ca. 201 to 182 Ma, whereas low-sulfidation systems were active at ca. 192 to 162 Ma. Amongst the studied epithermal systems, the Baker low-sulfidation epithermal deposit displays the strongest demonstrable genetic link with magmatic fluids; fluid inclusion studies demonstrate that its ore fluids were hot (>468°C), saline, and deposited metals at deep crustal depths (>2 km). Sulfur, C, O, and Pb isotope data confirm the involvement of a magmatic fluid, but also suggest that the ore fluid interacted with Asitka and Takla Group country rocks prior to metal deposition. In contrast, in the Shasta, Lawyers, and Griz-Sickle low-sulfidation epithermal systems, there is no clear association with magmatic fluids. Instead, their fluid inclusion data indicate the involvement of low-temperature (175 to 335°C), low-salinity (1 to 11 equiv. wt.% NaCl) fluids that deposited metals at shallow depths (<850 m). Their isotope (i.e., O, H, Pb) data suggest interaction between meteoric and/or metamorphic ore fluids with basement country rocks.     

Composition and origin of the Çaldağ oxide nickel laterite,W. Turkey

The Çalda? nickel laterite deposit located in the Aegean region of W. Turkey contains a reserve of 33 million tons of Ni ore with an average grade of 1.14% Ni. The deposit is developed on an ophiolitic serpentinite body which was obducted onto Triassic dolomites in the Late Cretaceous. The deposit weathering profile is both laterally and vertically variable. A limonite zone, which is the main ore horizon, is located at the base of the profile. A hematite horizon is located above the limonite, which in the south of the deposit is capped by Eocene freshwater limestones and in the north by a siliceous horizon. The deposit is unusual in lacking a significant saprolite zone with little development of Ni-silicates. The boundary between the limonite zone and serpentinite below is sharp with a marked decrease in concentrations of MgO from 13 to 1 wt.% over a distance of 2 mm representing the ‘Mg discontinuity’. Ni concentrations within goethite, the main ore mineral, reach a maximum of ~3 wt.% near the base of the limonite zone. Silica concentrations are high throughout most of the laterite with up to 80 wt.% silica in the upper portion of some profiles. The combination of a serpentinite protolith and a high water table at Çalda?, in association with an aggressive weathering environment in a tropical climate, resulted in the formation of an oxide-dominated deposit. The precipitation of silica may coincide with a change in climate with silica precipitation linked to an increase in seasonality. Additional variations within profile morphology are attributed to transportation during and after laterite development as a result of faulting, pocket type laterite formation and slumping, each of which produces a contrasting set of textural and geochemical features.