Multivariate classification methods in planetary sciences

In this paper the problem of the classification of natural samples is discussed. An updated version of the G-mode multivariate statistical method for the classification of natural samples, applicable to a wide range of research fields, is discussed in this paper. This method allows an automatic classification in terms of homogeneous taxonomic units, without any a priori knowledge of the taxonomic structure of the natural observations; it provides informations on the different levels of classification present in the data set under study (classes and subclasses), on the level of information residing in each variable, on the level of similarity and/or difference among homogeneous classes.An earlier version of this method has been widely applied in planetary sciences, astrophysics and geological sciences. The authors give here a resumé of the most interesting results obtained in these different research fields, from the geochemistry of lunar samples to asteroids taxonomy to remote sensing of planetary surfaces. This method is extremely reliable and versatile, and it is suggested that its use be widespread whenever the problem of the classification of a large set of natural samples occurs.     

Five micron limb-darkening and the structure of the Jovian atmosphere

By observing the transit of various cloud features across the Jovian disk, Terrile and Westphal (1977) have constructed limb-darkening curves for three regions in the 4.6 to 5.1 μm band. Several models currently employed in describing the radiative or dynamical properties of planetary atmospheres are here examined to understand their implications for limb-darkening. The statistical problem of fitting these models to the observed data is reviewed and methods for applying multiple regression analysis are discussed. Analysis of variance techniques are introduced to test the viability of a given physical process as a cause of the observed limb-darkening. The intermediate flux region of the North Equatorial Belt appears to be in only modest departure from radiative equilibrium. The limb-darkening curve for the South Temperate Belt is rich in structure and cannot be satisfactorily ascribed to any single physical mechanism; a combination of several, as yet unidentified, processes is likely involved. The hottest areas of the North and South Equatorial Belts exhibit limb-darkening curves that are typical of atmospheres in convective equilibrium. In this case, we derive a measure of the departure of the lapse rate from the dry adiabatic value (η?1.68), which furnishes strong evidence for a phase transition at unit optical depth in the NEB and SEB. Although the system NH₃H₂S cannot be entirely ruled out, the freezing of an aqueous ammonia solution is shown to be consistent with the parameter fit and solar abundance data, while being in close agreement with Lewis' (1969a) cloud models.     

Serpentinization and the inorganic synthesis of H2 in planetary surfaces

The near-surface inorganic synthesis of molecular hydrogen (H₂) is a fundamental process relevant to the origins and to the sustenance of early life on Earth and potentially other planets. Hydrogen production through the decomposition of water is thought to be a principal reaction that occurs during hydrothermal alteration of olivine, an iron-magnesium silicate abundant near planetary surfaces. We demonstrate that copious amounts of H₂ are produced only when the olivine undergoing alteration (serpentinization) contains 1 to 50 mol% iron over a variety of planetary surface P-T conditions. This suggests that extrasolar Earth-like planets that are hosted by a star with iron contents up to two times the solar value could support life provided they are hydrothermally active and fall within the habitable zone around the star.     

Mars: The evolutionary history of the northern lowlands based on crater counting and geologic mapping

The geologic history of planetary surfaces is most effectively determined by joining geologic mapping and crater counting which provides an iterative, qualitative and quantitative method for defining relative ages and absolute model ages. Based on this approach, we present spatial and temporal details regarding the evolution of the Martian northern plains and surrounding regions.The highland–lowland boundary (HLB) formed during the pre-Noachian and was subsequently modified through various processes. The Nepenthes Mensae unit along the northern margins of the cratered highlands, was formed by HLB scarp-erosion, deposition of sedimentary and volcanic materials, and dissection by surface runoff between 3.81 and 3.65 Ga. Ages for giant polygons in Utopia and Acidalia Planitiae are $～$ 3.75 Ga and likely reflect the age of buried basement rocks. These buried lowland surfaces are comparable in age to those located closer to the HLB, where a much thinner, post-HLB deposit is mapped. The emplacement of the most extensive lowland surfaces ended between 3.75 and 3.4 Ga, based on densities of craters generally $>3\mathrm{km}$ in diameter. Results from the polygonal terrain support the existence of a major lowland depocenter shortly after the pre-Noachian formation of the northern lowlands. In general, northern plains surfaces show gradually younger ages at lower elevations, consistent local to regional unit emplacement and resurfacing between 3.6 and 2.6 Ga. Elevation levels and morphology are not necessarily related, and variations in ages within the mapped units are found, especially in units formed and modified by multiple geological processes. Regardless, most of the youngest units in the northern lowlands are considered to be lavas, polar ice, or thick mantle deposits, arguing against the ocean theory during the Amazonian Period (younger than about 3.15 Ga).All ages measured in the closest vicinity of the steep dichotomy escarpment are also 3.7 Ga or older. The formation ages of volcanic flanks at the HLB (e.g., Alba Mons (3.6–3.4 Ga) and the last fan at Apollinaris Mons, 3.71 Ga) may give additional temporal constraint for the possible existence of any kind of Martian ocean before about 3.7 Ga. It seems to reflect the termination of a large-scale, precipitation-based hydrological cycle and major geologic processes related to such cycling.     

The specific molecular identification of life experiment (SMILE)

We describe a compact, highly integrated instrument concept for the detection and identification of a wide range of molecules associated with extinct/extant life or potential life processes. The Specific Molecular Identification of Life Experiment (SMILE) will be sensitive to the presence of a range of target molecules using both electrical and optical transduction techniques, and incorporates molecular imprinted polymers in addition to traditional biological receptors such as antibodies. A number of versions of the concept are possible depending on available resources e.g. mass, volume, etc. The full concept utilises a novel imaging interferometer where a large number of molecular receptors are deposited on the measurement plane of an imaging interferometer and read out by an imaging detector, enabling multiple targets—biomarkers—within a sample to be measured simultaneously. The optics can also form the basis of an UV-NIR imaging Fourier spectrometer allowing basic mineralogy studies to be conducted using optical properties to assist in the determination of the geological context of the samples. By incorporating micro-fabricated transducer arrays, micro-fluidics and artificial molecular recognition systems, as well as recombinant antibody technology with appropriate integration methods, SMILE forms a compact and robust “Life Marker Chip” which has been proposed for future planetary missions including ESA's ExoMars mission, where the instrument offers the possibility of conducting a direct in situ search for signs of past or present biological activity on Mars. In addition to its role in planetary exploration, derivatives of SMILE have multiple terrestrial applications in fields such as forensic analysis and environmental monitoring.     

Implementation of cartographic symbols for planetary mapping in geographic information systems

The steadily growing international interest in the exploration of planets in our Solar System and many advances in the development of space-sensor technology have led to the launch of a multitude of planetary missions to Mercury, Venus, the Earth's moon, Mars and various Outer-Solar System objects, such as the Jovian and Saturnian satellites. Camera instruments carried along on these missions image surfaces in different wavelength ranges and under different viewing angles, permitting additional data to be derived, such as spectral data or digital terrain models. Such data enable researchers to explore and investigate the development of planetary surfaces by analyzing and interpreting the inventory of surface units and structures. Results of such work are commonly abstracted and represented in thematic, mostly geological and geomorphological, maps. In order to facilitate efficient collaboration among different planetary research disciplines, mapping results need to be prepared, described, managed, archived, and visualized in a uniform way. These tasks have been increasingly carried out by means of computer-based geographic information systems (GIS or GI systems) which have come to be widely employed in the field of planetary research since the last two decades. In this paper we focus on the simplification of mapping processes, putting specific emphasis on a cartographically correct visualization of planetary mapping data using GIS-based environments. We present and discuss the implementation of a set of standardized cartographic symbols for planetary mapping based on the Digital Cartographic Standard for Geologic Map Symbolization as prepared by the United States Geological Survey (USGS) for the Federal Geographic Data Committee (FGDC). Furthermore, we discuss various options to integrate this symbol catalog into generic GI systems, and more specifically into the Environmental Systems Research Institute's (ESRI) ArcGIS environment, and focus on requirements for symbol definitions in the field of planetary mapping. A symbology of this type can be embedded into any modular GIS environment capable in dealing with external stand-alone as well as database-driven management of symbol sets. Using such a uniform GIS-based symbol catalog will give the research community access to map results already cartographically elaborated, enabling them to create digital maps as a secondary data source in subsequent studies.     

Pulsations and planets: The asteroseismology‐extrasolar‐planet connection

The disciplines of asteroseismology and extrasolar planet science overlap methodically in the branch of high‐precision photometric time series observations. Light curves are, amongst others, useful to measure intrinsic stellar variability due to oscillations, as well as to discover and characterize those extrasolar planets that transit in front of their host stars, periodically causing shallow dips in the observed brightness. Both fields ultimately derive fundamental parameters of stellar and planetary objects, allowing to study for example the physics of various classes of pulsating stars, or the variety of planetary systems, in the overall context of stellar and planetary system formation and evolution. Both methods typically also require extensive spectroscopic follow‐up to fully explore the dynamic characteristics of the processes under investigation. In particularly interesting cases, a combination of observed pulsations and signatures of a planet allows to characterize a system's components to a very high degree of completeness by combining complementary information. The planning of the relevant space missions has consequently converged with respect to science cases, where at the outset there was primarily a coincidence in instrumentation and techniques. Whether space‐ or ground‐based, a specific type of stellar pulsations can themselves be used in an innovative way to search for extrasolar planets. Results from this additional method at the interface of stellar pulsation studies and exoplanet hunts in a beyond‐mainstream area are presented (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the properties of O and O2 ions in a hybrid model and in Mars Express IMA/ASPERA-3 data: A case study

We have performed a numerical simulation to analyze the energy spectra of escaping planetary O⁺ and O₂⁺ ions at Mars. The simulated time-energy spectrograms were generated along orbit no. 555 (June 27, 2004) of Mars Express when its Ion Mass Analyzer (IMA)/ASPERA-3 ion instrument detected escaping planetary ions. The simulated time-energy spectrograms are in general agreement with the hypothesis that planetary O⁺ and O₂⁺ ions far from Mars are accelerated by the convective electric field. The HYB-Mars hybrid model simulation also shows that O⁺ ions originating from the ionized hot oxygen corona result in a high-energy (E>1 keV) O⁺ ion population that exists very close to Mars. In addition, the simulation also results in a low-energy (E<0.1 keV) planetary ion population near the pericenter. In the analyzed orbit, IMA did not observe a clear high-energy planetary ion or a clear low-energy planetary ion population near Mars. One possible source for this discrepancy may be the Martian magnetic crustal anomalies because MEX passed over a strong crustal field region near the pericenter, but the hybrid model does not include the magnetic crustal anomalies.     

Global survey of color variations on 433 Eros: Implications for regolith processes and asteroid environments

High-resolution imaging acquired with the Near Earth Asteroid Rendezvous Shoemaker (NEAR Shoemaker) spacecraft is used to elucidate the spectral properties and spatial distribution of color units on Asteroid 433 Eros. Previous workers mapped four distinct types of color units on the surface (bright streaks, dark soils, ponded materials, average regolith). These units exhibit albedo and color boundaries but there is no evidence to indicate they represent distinct rock types. Rather the units are thought to show evidence of complex regolith transport and sorting processes. Here we report the results of a comprehensive study of all viable color MultiSpectral Imager (MSI) data to identify and characterize the distribution and nature of color units across the whole asteroid. Due to a spacecraft upset that resulted in contamination of the MSI optics, color images are affected with a scattered light problem that hampers interpretation of subtle color contrasts, even after a rigorous remediation. To constrain interpretations of the MSI color data we characterize this residual scattered light and demonstrate how complete correction would affect derived color ratios. Results of our comprehensive study are consistent with previous mapping—confirming that bright streaks, average regolith and dark soils fall on a mixing line, consistent with space weathering effects. We find that the ponded deposits do not fall on this putative mixing line. The color and reflectance of the ponded deposits are consistent with some combination of compositional, grain size and maturity variations from the average regolith. Additionally we show that spectral separation of the four units on ratio plots would only increase with full removal of residual scattered light, especially for features that are small in terms of pixels. Global analysis of the Eros color units illustrates complex regolith processes and grain sorting that may hold clues to understanding space weathering processes and the link between asteroids and meteorites.     

Do planetary encounters reset surfaces of near Earth asteroids?

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these ‘space weathering’ processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r^∗, at which planetary encounters can reset the surface. We found that t_sw ∼ 10⁶ yr (at 1 AU) and r^∗ ∼ 5R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 10⁵ yr would require that r^∗ > 20R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 10⁷ yr. We found that t_sw ∝ q², where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.     

The frequency-area distribution of volcanic units on Venus: Implications for planetary resurfacing

The areas of volcanic units on Venus have been measured on the 1:5000000 geological maps published by NASA/USGS. These data were used to obtain a frequency-area distribution. The cumulative frequency-area distribution of 1544 specific occurrence of units cover six orders of magnitude from the largest unit (30 × 10⁶ km²) to the smallest (20 km²). The probability distribution function has been calculated. The medium and large volcanic units correlate well with a power-law (fractal) relation for the dependence of frequency on area with a slope of −1.83. There are fewer small units than the expected values provided by the power-law relation. Our measurements cover 21.02% of the planetary surface, 3.59% of the study area was found to be tessera terrain and is excluded from this study of volcanism. The measurements were restricted to areas where geological maps have been published. The analysis was performed on two independent areas of the planet, with a complete coverage of published maps. In both areas the largest volcanic unit covers a significant portion of the surface (58.75% and 63.64%, respectively). For the total measured volcanic units (excluding tessera), these two largest units (that could correspond to the same unit or not) cover the 61.18% and they are stratigraphically superimposed on older volcanic units which cover 3.37% of the area. The remaining area (35.45%) is occupied by younger volcanic units stratigraphically superimposed on the large volcanic unit(s). These results are based on the independent mapping of a large number of geologists with different ideas about the geodynamical evolution of Venus and different criteria for geological mapping. Despite this fact, the presence of these very large units is incompatible with the equilibrium resurfacing models, because their generation at different ages would destroy the crater randomness. Our frequency-area distribution of the mapped volcanic units supports a catastrophic resurfacing due to the emplacement of the largest unit(s) followed by a decay of volcanism. Our data for the frequency-area distribution of volcanic units provide new support for catastrophic resurfacing models. It is difficult to make our observations compatible with equilibrium, steady-state resurfacing models.     

Obliquity evolution of extrasolar terrestrial planets

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ∼1 AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio (m/Mp) and initial obliquity (γ₀), through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized or larger bodies determine the planetary spin and form satellites. Since the giant impacts would be isotropic, tilted spins (sinγ₀∼1) are more likely to be produced than straight ones (sinγ₀∼0). The ratio m/Mp is dependent on the impact parameters and impactors' mass. However, most of previous studies on tidal evolution of the planet-satellite systems have focused on a particular case of the Earth-Moon systems in which m/Mp?0.0125 and γ₀∼10° or the two-body planar problem in which γ₀=0° and stellar torque is neglected. We numerically integrated the evolution of planetary spin and a satellite orbit with various m/Mp (from 0.0025 to 0.05) and γ₀ (from 0° to 180°), taking into account the stellar torques and precessional motions of the spin and the orbit. We start with the spin axis that almost coincides with the satellite orbit normal, assuming that the spin and the satellite are formed by one dominant impact. With initially straight spins, the evolution is similar to that of the Earth-Moon system. The satellite monotonically recedes from the planet until synchronous state between the spin period and the satellite orbital period is realized. The obliquity gradually increases initially but it starts decreasing down to zero as approaching the synchronous state. However, we have found that the evolution with initially tiled spins is completely different. The satellite's orbit migrates outward with almost constant obliquity until the orbit reaches the critical radius ∼10-20 planetary radii, but then the migration is reversed to inward one. At the reversal, the obliquity starts oscillation with large amplitude. The oscillation gradually ceases and the obliquity is reduced to ∼0° during the inward migration. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. We found that the character change of precession about total angular momentum vector into that about the planetary orbit normal is responsible for the oscillation with large amplitude and the reversal of migration. With the results of numerical integration and analytical arguments, we divided the m/Mp-γ₀ space into the regions of the qualitatively different evolution. The peculiar tidal evolution with initially tiled spins give deep insights into dynamics of extrasolar planet-satellite systems and discussions of surface environments of the planets.     

A Model for the Common Origin of Jupiter Family and Halley Type Comets

A numerical simulation of the Oort cloud is used to explain the observed orbital distributions and numbers of Jupiter-family (JF) and Halley-type (HT) short-period (SP) comets. Comets are given initial orbits with perihelion distances between 5 and 36 au, and evolve under planetary, stellar and Galactic perturbations for 4.5 Gyr. This process leads to the formation of an Oort cloud (which we define as the region of semimajor axes a > 1,000 au), and to a flux of cometary bodies from the Oort cloud returning to the planetary region at the present epoch. The results are consistent with the dynamical characteristics of SP comets and other observed cometary populations: the near-parabolic flux, Centaurs, and high-eccentricity trans-Neptunian objects. To achieve this consistency with observations, the model requires that the number of comets versus initial perihelion distance is concentrated towards the outer planetary region. Moreover, the mean physical lifetime of observable comets in the inner planetary region (q < 2.5 au) at the present epoch should be an increasing function of the comets’ initial perihelion distances. Virtually all observed HT comets and nearly half of observed JF comets come from the Oort cloud, and initially (4.5 Gyr ago) from orbits concentrated near the outer planetary region. Comets that have been in the Oort cloud also return to the Centaur (5 < q < 28 au, a < 1,000 au) and near-Neptune high-eccentricity regions. Such objects with perihelia near Neptune are hard to discover, but Centaurs with characteristics predicted by the model (e.g. large semimajor axes, above 60 au, or high inclinations, above 40°) are increasingly being found by observers. The model provides a unified picture for the origin of JF and HT comets. It predicts that the mean physical lifetime of all comets in the region q < 1.5 au is less than ～200 revolutions.     

A Model of Relation between Fluctuation of Double Differential Total Ionospheric Electron Content and Angular Distance of the Two Satellites Observed by Same-beam VLBI

Time delay and phase fluctuation are produced when the signals of a spacecraft are transmitted through the ionosphere of the earth, which give rise to a great influence on the measurement precision of VLBI (Very Long Baseline Interferometry). Using the 1-year same-beam VLBI data of 2 satellites (Rstar and Vstar) in the Japanese lunar exploration project SELENE, we obtained a model of the relation between the fluctuation of double differential total electron content in the ionosphere and the angular distance of the two satellites. For the 6 baselines, the root mean square r of fluctuation (in units of TECU) and the angular distance of the two satellites θ (in units of ^°) has a relation of r = 0.773θ + 0.562, and for the 4 VLBI stations, the relation is r = 0.554θ + 0.399 from the baselines inversion. The results can serve as a reference for the derivation of differential phase delay and for the occultation observation and study of planetary ionospheres.     

Surface roughness and geological mapping at subhectometer scale from the High Resolution Stereo Camera onboard Mars Express

The quantitative measurement of surface roughness of planetary surfaces at all scales provides insights into geological processes. A characterization of roughness variations at the scale of a few tens of meters is proposed that complements the analysis of local topographic data of the martian surface at kilometer scale, as achieved from the Mars Orbiter Laser Altimeter (MOLA) data, and at the subcentimeter scale using photometric properties derived from multi-angular observations. Relying on a Gabor filtering process, an algorithm developed in the context of image classification for the purpose of texture analysis has been adapted to handle data from the High Resolution Stereo Camera (HRSC). The derivation of roughness within a wavelength range of tens of meters, combined with analyses at even longer wavelengths, gives an original view of the martian surface. The potential of this approach is evaluated for different examples for which the geological processes are identified and the geological units are mapped and characterized in terms of roughness.     

A system of astronomical constants in the relativistic framework

The relation between the units and readings of time and space coordinates of terrestrial and barycentric reference frames is discussed from the viewpoint of general relativity. Attention is paid to the unit of space coordinates since the International Astronomical Union (IAU) regulates only the unit of time in the above two frames. Two definitions of unit of length are examined and their effects on the numerical expression of coordinate transformation, equations of planetary motions, and those for light propagation time are discussed. A clear conflict is found between the IAU (1976) recommendation on the definition of the time-scales in different frames of reference and the statement that all constants in the IAU (1976) new system of astronomical constants are defined in terms of the Internationsl System of units (SI units). One of the above two definitions is proposed to resolve this conflict by the least alteration to current procedures for analysing the recent astrometric observations such as the radar/laser rangings, the range and range-rate, and the very long baseline interferometric observations. Also, an interpretation of numerical values in the IAU (1976) new system of astronomical constants is presented. It is stressed that the definition proposed in this paper requires that a formula slightly different from that in current use be employed in the numerical transformation of readings of coordinates between the terrestrial and barycentric reference frames.     

The Zamama-Thor region of Io: Insights from a synthesis of mapping, topography, and Galileo spacecraft data

We have studied data from the Galileo spacecraft's three remote sensing instruments (Solid-State Imager (SSI), Near-Infrared Mapping Spectrometer (NIMS), and Photopolarimeter-Radiometer (PPR)) covering the Zamama-Thor region of Io's antijovian hemisphere, and produced a geomorphological map of this region. This is the third of three regional maps we are producing from the Galileo spacecraft data. Our goal is to assess the variety of volcanic and tectonic materials and their interrelationships on Io using planetary mapping techniques, supplemented with all available Galileo remote sensing data. Based on the Galileo data analysis and our mapping, we have determined that the most recent geologic activity in the Zamama-Thor region has been dominated by two sites of large-scale volcanic surface changes. The Zamama Eruptive Center is a site of both explosive and effusive eruptions, which emanate from two relatively steep edifices (Zamama Tholi A and B) that appear to be built by both silicate and sulfur volcanism. A ∼100-km long flow field formed sometime after the 1979 Voyager flybys, which appears to be a site of promethean-style compound flows, flow-front SO₂ plumes, and adjacent sulfur flows. Larger, possibly stealthy, plumes have on at least one occasion during the Galileo mission tapped a source that probably includes S and/or Cl to produce a red pyroclastic deposit from the same vent from which silicate lavas were erupted. The Thor Eruptive Center, which may have been active prior to Voyager, became active again during the Galileo mission between May and August 2001. A pillanian-style eruption at Thor included the tallest plume observed to date on Io (at least 500 km high) and new dark lava flows. The plume produced a central dark pyroclastic deposit (probably silicate-rich) and an outlying white diffuse ring that is SO₂-rich. Mapping shows that several of the new dark lava flows around the plume vent have reoccupied sites of earlier flows. Unlike most of the other pillanian eruptions observed during the Galileo mission, the 2001 Thor eruption did not produce a large red ring deposit, indicating a relative lack of S and/or Cl gases interacting with the magma during that eruption. Between these two eruptive centers are two paterae, Thomagata and Reshef. Thomagata Patera is located on a large shield-like mesa and shows no signs of activity. In contrast, Reshef Patera is located on a large, irregular mesa that is apparently undergoing degradation through erosion (perhaps from SO₂-sapping or chemical decomposition of sulfur-rich material) from multiple secondary volcanic centers.     

Fundamental physics and absolute positioning metrology with the MAGIA lunar orbiter

MAGIA is a mission approved by the Italian Space Agency (ASI) for Phase A study. Using a single large-diameter laser retroreflector, a large laser retroreflector array and an atomic clock onboard MAGIA we propose to perform several fundamental physics and absolute positioning metrology experiments: VESPUCCI, an improved test of the gravitational redshift in the Earth?CMoon system predicted by General Relativity; MoonLIGHT-P, a precursor test of a second generation Lunar Laser Ranging (LLR) payload for precision gravity and lunar science measurements under development for NASA, ASI and robotic missions of the proposed International Lunar Network (ILN); Selenocenter (the center of mass of the Moon), the determination of the position of the Moon center of mass with respect to the International Terrestrial Reference Frame/System (ITRF/ITRS); this will be compared to the one from Apollo and Lunokhod retroreflectors on the surface; MapRef, the absolute referencing of MAGIA??s lunar altimetry, gravity and geochemical maps with respect to the ITRF/ITRS. The absolute positioning of MAGIA will be achieved thanks to: (1) the laboratory characterization of the retroreflector performance at INFN-LNF; (2) the precision tracking by the International Laser Ranging Service (ILRS), which gives two fundamental contributions to the ITRF/ITRS, i.e. the metrological definition of the geocenter (the Earth center of mass) and of the scale of length; (3) the radio science and accelerometer payloads; (4) support by the ASI Space Geodesy Center in Matera, Italy. Future ILN geodetic nodes equipped with MoonLIGHT and the Apollo/Lunokhod retroreflectors will become the first realization of the International Moon Reference Frame (IMRF), the lunar analog of the ITRF.     

Interstellar dust in the solar system

Dust is an important component of galactic stucture and the cyclic processing of particulate matter leads to stellar and planetary formation. Though astronomical methods using analysis of dust-penetrating starlight can provide some limited information about the dust, the prospect of its in-situ sampling within the Solar System by spacecraft and its remote sensing by ground-based techniques open up a new field in galactic exploration.     

Composition and thermal inertia of the Mawrth Vallis region of Mars from TES and THEMIS data

Clay mineral-bearing deposits previously discovered on Mars with near infrared (λ=0.3-5 μm) remote sensing data are of major significance for understanding the aqueous history, geological evolution, and past habitability of Mars. In this study, we analyzed the thermal infrared (λ=6-35 μm) surface properties of the most extensive phyllosilicate deposit on Mars: the Mawrth Vallis area. Clay mineral-bearing units, which in visible images appear to be relatively light-toned, layered bedrock, have thermal inertia values ranging from 150 to 460 J m⁻² K⁻¹ s^−1/2. This suggests the deposits are composed of a mixture of rock with sand and dust at 100-meter scales. Dark-toned materials that mantle the clay-bearing surfaces have thermal inertia values ranging from 150 to 800, indicating variable degrees of rockiness or induration of this younger sedimentary or pyroclastic unit. Thermal Emission Spectrometer (TES) spectra of the light-toned rocks were analyzed with a number of techniques, but none of the results shows a large phyllosilicate component as has been detected in the same surfaces with near-infrared data. Instead, TES spectra of light-toned surfaces are best modeled by a combination of plagioclase feldspar, high-silica materials (similar to impure opaline silica or felsic glass), and zeolites. We propose three hypotheses for why the clay minerals are not apparent in thermal infrared data, including effects due to surface roughness, sub-pixel mixing of multiple surface temperatures, and low absolute mineral abundances combined with differences in spatial sampling between instruments. Zeolites modeled in TES spectra could be a previously unrecognized component of the alteration assemblage in the phyllosilicate-bearing rocks of the Mawrth Vallis area. TES spectral index mapping suggests that (Fe/Mg)-clays detected with near infrared data correspond to trioctahedral (Fe²⁺) clay minerals rather than nontronite-like clays. The average mineralogy and geologic context of these complex, interbedded deposits suggests they are either aqueous sedimentary rocks, altered pyroclastic deposits, or a combination of both.