Emissivity measurements of analogue materials for the interpretation of data from PFS on Mars Express and MERTIS on Bepi-Colombo

The Planetary Fourier Spectrometer (PFS) onboard the Mars Express Mission provides thermal infrared, hyperspectral images of Mars and in the future the Mercury Radiometer and Thermal Infrared Spectrometer (MERTIS), part of the selected payload for the Bepi-Colombo Mission, will collect analogous data for Mercury. To interpret these remote sensing data it is essential to understand the spectral emittance of planetary analogue materials and a spectral library of emissivity measurements is needed. Here we introduce the emissivity device built at DLR (Berlin). The device is coupled to a Fourier transform infrared spectrometer (Bruker IFS 88), purged with dry air and equipped with a cooled MCT-detector. We discuss theoretical background of our thermal emission measurements and describe our standard experimental procedures, which are being used to create the Berlin emissivity database (BED). This study presents and discuss the 6.3-22 μm thermal emission spectra of fine-grained feldspar separates ranging from <25 to 90-125 μm. We discuss diagnostic potential of the features present in the emission spectra of plagioclase and alkali feldspars and the particle size effects.     

Exponential Gaussian approach for spectral modelling: The EGO algorithm II. Band asymmetry

The present investigation is complementary to a previous paper which introduced the EGO approach to spectral modelling of reflectance measurements acquired in the visible and near-IR range (Pompilio, L., Pedrazzi, G., Sgavetti, M., Cloutis, E.A., Craig, M.A., Roush, T.L. [2009]. Icarus, 201 (2), 781-794). Here, we show the performances of the EGO model in attempting to account for temperature-induced variations in spectra, specifically band asymmetry.Our main goals are: (1) to recognize and model thermal-induced band asymmetry in reflectance spectra; (2) to develop a basic approach for decomposition of remotely acquired spectra from planetary surfaces, where effects due to temperature variations are most prevalent; (3) to reduce the uncertainty related to quantitative estimation of band position and depth when band asymmetry is occurring.In order to accomplish these objectives, we tested the EGO algorithm on a number of measurements acquired on powdered pyroxenes at sample temperature ranging from 80 up to 400 K. The main results arising from this study are: (1) EGO model is able to numerically account for the occurrence of band asymmetry on reflectance spectra; (2) the returned set of EGO parameters can suggest the influence of some additional effect other than the electronic transition responsible for the absorption feature; (3) the returned set of EGO parameters can help in estimating the surface temperature of a planetary body; (4) the occurrence of absorptions which are less affected by temperature variations can be mapped for minerals and thus used for compositional estimates.Further work is still required in order to analyze the behaviour of the EGO algorithm with respect to temperature-induced band asymmetry using powdered pyroxene spanning a range of compositions and grain sizes and more complex band shapes.     

Atmospheric thermal structures of the giant planets

Thermal models of planetary atmospheres can be calculated from assumptions of the energy budget of the atmosphere and from the knowledge of the effective temperature of the studied planet. On the other hand, the retrieval of the thermal atmospheric profiles from infrared measurements by means of the numerical inversion of the radiative transfer equation presents the advantages of not requiring such assumptions. The extent of the atmospheric range which can then be sounded is examined and the vertical resolution of the inferred profiles is discussed. Comparisons of thermal models and retrieved thermal profiles are made for the four giant planets. The retrieved profiles lead to brightness temperature spectra which fit all the available infrared measurements fairly well for Jupiter and Saturn but only part of them for Uranus and Neptune. The values of the planetary effective temperatures calculated from the retrieved profiles show that Jupiter, Saturn, and Neptune have strong internal heating sources while Uranus probably has a very small or null one.     

Refractive index measurements of ammonia and hydrocarbon ices at 632.8 nm

Optical constants in a broad temperature and wavelength range are important input parameters in radiative transfer models used in studies of planetary atmospheres. In the laboratory, the refractive index values of ices at the HeNe laser wavelength (632.8 nm) are often used to monitor the growth rate and thickness of ice films. In this report we present laboratory measurements determining the refractive index at 632.8 nm of ammonia and hydrocarbon ices in the temperature range 80-100 K. Thin ice films are vapor-deposited on a cryogenically cooled mirror located inside a high-vacuum apparatus. The real component of the refractive index of these ice films is determined by a two-angle interferometric technique. Optical modeling calculations of the transmittance and reflectance through the thin ice films assist in the interpretation of the experimental results. We discuss our results and compare them with other measurements available in the literature. The results reported here are relevant to the spectroscopy of icy objects in the solar system; they are needed to perform laboratory characterization of ices, derive optical constants, and model spectra.     

Photometric properties of powdered sulfur

The spectrophotometric (0.39 < λ < 0.7 λm) properties of three particle-size fractions (diameters <10 λm, <150 λm, and 420–850 λm) of sulfur have been investigated in the laboratory. Particle size, temperature, thermal history, and scattering geometry are all shown to influence the spectral reflectance of the normal (S₈) sulfur samples and an “orange-colored” S₈ sample produced by quenching molten sulfur. A scattering law consisting of a linear combination of lunar-like and Lambertian terms adequately describes the data for all particle sizes. Where sulfur is darkest (λ < 0.45 λm), the reflectance decreases with increasing particle size, whereas where sulfur is brightest (λ > 0.45 λm) the reflectance increases with decreasing particle size. In reflected light, the long wavelength edge of the strong ultraviolet absorption retreats smoothly to shorter wavelengths with decreasing temperature at ～1.6Å/°K, a value lower than the 2.2Å/°K value previously reported for transmitted light. Near opposition, sulfur powders are found to follow closely a Minnaert limb darkening law except where the reflectance is low, i.e., in the strong ultraviolet absorption band of the larger particle size fractions. It is clear from our data that quantitative comparisons between disk-integrated observations of Io and laboratory measurements of flat samples of sulfur are not adequate unless temperature effects and changes in scattering geometry are included.     

Interpretation of the lunar microwave brightness temperature spectrum: Feasibility of orbital heat flow mapping

A detailed model of the lunar regolith is analyzed to examine the feasibility of an orbital mapping of heat flow using microwave radiometers. For regolith thermal and electrical properties which are representative of Apollo findings, brightness temperature observations in the bandλ = 5–30 cm would be required for heat flow analysis. Spectral variations shortward of 5 cm are controlled primarily by the temperature dependencies of the thermal conductivity and electrical absorption within the diurnal-varying layer. For wavelengths longer than 30 cm, unwanted emission from high impedance subregolith layers can be significant and size limitations on spacecraft radiometers is a factor. Over the 5- to 30-cm band, lunation-averaged brightness temperature increases of 2–10°K are predicted for heat flow values representative of the Apollo measurements. The magnitude of this increase depends directly on the value of regolith microwave absorption. For absorption values consistent with Apollo laboratory measurements, a spectral increase of 5°K is predicted. This value is considered marginally sufficient for an orbital heat flow measurement. However, important non-heat flow effects must be accounted for. Spectral variations can occur due to surface topography and subsurface scattering. For nadir viewing, surface roughness effects are not expected to be significant and topographic effects are nearly constant with wavelength for λ > 5cm. Volume scattering due to subsurface rock fragments can cause emission darkening of 1–6°K. However, spectral variations will not be large unless the distribution of scatterer sizes is sharply skewed. For the Moon, the most serious spurious effect appears to be emissivity variations due to the near-surface density gradient. A brightness temperature decrease of 10°K is predicted from centimeter to decameter wavelengths. If the transition from porous surface fines to compacted regolith soil occurs rapidly (within the upper 3–5 cm), most of the emissivity decrease will occur in the 5- to 30-cm wavelength band. It is recommended that complementary radar measurements be utilized to augment constraints on regolith emissivity and scattering properties.     

Thermal infrared (8-13 μm) spectra of 29 asteroids: the Cornell Mid-Infrared Asteroid Spectroscopy (MIDAS) Survey

We report the results of the Cornell Mid-IR Asteroid Spectroscopy (MIDAS) survey, a program of ground-based observations designed to characterize the 8-13 μm spectral properties of a statistically significant sample of asteroids from a wide variety of visible to near-IR spectral classes. MIDAS is conducted at Palomar Observatory using the Spectrocam-10 (SC-10) spectrograph on the 200-in Hale telescope. We have measured the mid-infrared spectra of twenty-nine asteroids and have derived temperature estimates from our data that are largely consistent with the predictions of the standard thermal model. We have also generated relative emissivity spectra for the target asteroids. On only one asteroid, 1 Ceres, have we found emissivity features with spectral contrast greater than 5%. Our spectrum of 4 Vesta suggests emissivity variation at the 2-3% level. Published spectra of several of the small number of asteroids observed with ISO (six of which are also included in our survey), which appeared to exhibit much stronger emissivity features, are difficult to reconcile with our measurements. Laboratory work on mineral and meteorite samples has shown that the contrast of mid-IR spectral features is greatly reduced at fine grain sizes. Moreover, the NEAR mission found that 433 Eros is covered by a relatively thick fine-grained regolith. If small bodies in general possess such regoliths, their mid-IR spectral features may be quite subtle. This may explain the evident absence of strong emissivity variation in the majority of the MIDAS spectra.     

Search for sulfates on the surface of Ceres

The formation of hydrated salts is an expected consequence of aqueous alteration of Main Belt objects, particularly for large, volatile‐rich protoplanets like Ceres. Sulfates, present on water‐bearing planetary bodies (e.g., Earth, Mars, and carbonaceous chondrite parent bodies) across the inner solar system, may contribute to Ceres’ UV and IR spectral signature along with phyllosilicates and carbonates. We investigate the presence and stability of hydrated sulfates under Ceres’ cryogenic, low‐pressure environment and the consequent spectral effects, using UV–Vis–IR reflectance spectroscopy. H₂O loss begins instantaneously with vacuum exposure, measured by the attenuation of spectral water absorption bands, and a phase transition from crystalline to amorphous is observed for MgSO₄·6H₂O by X‐ray powder diffraction. Long‐term (>40 h), continuous exposure of MgSO₄·nH₂O (n = 0, 6, 7) to low pressure (10^?3–10^?6 Torr) causes material decomposition and strong UV absorption below 0.5 μm. Our measurements suggest that MgSO₄·6H₂O grains (45–83 μm) dehydrate to 2% of the original 1.9 μm water band area over ~0.3 Ma at 200 K on Ceres and after ~42 Ma for 147 K. These rates, inferred from an Avrami dehydration model, preclude MgSO₄·6H₂O as a component of Ceres’ surface, although anhydrous and minimally hydrated sulfates may be present. A comparison between Ceres emissivity spectra and laboratory reflectance measurements over the infrared range (5–17 μm) suggests sulfates cannot be excluded from Ceres’ mineralogy.     

In situ methods for measuring thermal properties and heat flux on planetary bodies

The thermo-mechanical properties of planetary surface and subsurface layers control to a high extent in which way a body interacts with its environment, in particular how it responds to solar irradiation and how it interacts with a potentially existing atmosphere. Furthermore, if the natural temperature profile over a certain depth can be measured in situ, this gives important information about the heat flux from the interior and thus about the thermal evolution of the body. Therefore, in most of the recent and planned planetary lander missions experiment packages for determining thermo-mechanical properties are part of the payload. Examples are the experiment MUPUS on Rosetta's comet lander Philae, the TECP instrument aboard NASA's Mars polar lander Phoenix, and the mole-type instrument HP³ currently developed for use on upcoming lunar and Mars missions. In this review we describe several methods applied for measuring thermal conductivity and heat flux and discuss the particular difficulties faced when these properties have to be measured in a low pressure and low temperature environment. We point out the abilities and disadvantages of the different instruments and outline the evaluation procedures necessary to extract reliable thermal conductivity and heat flux data from in situ measurements.     

Radiative transfer in the atmosphere of Venus and application to surface emissivity retrieval from VIRTIS/VEX measurements

The radiation measurements of VIRTIS-M-IR (1-5 μm) on Venus Express provide a valuable database for systematic studies of the atmosphere and surface of the Earth’s sister planet. The present paper focuses on the investigation of physical parameters that determine the retrieval accuracy of deep atmosphere and surface features of Venus including compositional conditions, continuum absorption effects, and spectroscopic input data required for radiative transfer simulations. The parameter discussion shall serve as a reference for ongoing and future work on methodical and simulation input data improvements.The high variability of the nightside atmospheric and surface emission window radiances with respect to cloud opacity and surface elevation is modeled and discussed in direct comparison with measurements performed over the northern hemisphere. Venus surface elevation is retrieved using the 1.18 and 1.02 μm emission windows where radiance ratios are well suited to de-cloud the measurement data. In general, the ratio-based VIRTIS topography is in good agreement with the Magellan topography, but differences occur in localized areas. The paper discusses possible origins of such differences including surface emissivity “anomalies”. Surface emissivity variations that may be due to changes in the chemical composition (mineralogy) and surface texture are important indicators of the nature of the surface material.Preliminary radiance retrievals along a number of complete northern orbits reveal a trend towards lower values of highland surface emissivity compared to the surrounding lowlands. Already the Magellan radar experiment suggested compositional variations at moderately high altitudes over the tesserae. They probably indicate a more felsic component giving a hint to older surface forming processes.     

C ring fine structures revealed in the thermal infrared

We analyze data sets obtained with the Composite Infrared Spectrometer (CIRS) onboard the Cassini spacecraft after the Saturn Orbit Insertion (SOI). Using the mid-IR interferometer's FP3 channel (600-1100 cm⁻¹), we derive radial temperature profiles for the C ring with a spatial resolution never achieved before. For the first time, the C ring's plateaus and ringlets can be clearly separated from the optically thinner background and their thermal behavior is studied separately for different viewing geometries. In particular, thermal phase curves derived for the plateaus reveal an interesting surge near 0° phase, not observed in the background. We show that mutual shadowing in the plateaus can explain the existence of the surge but is not sufficient to model the phase curves in detail. By analogy with thermal emission of asteroid surfaces we discuss the possible influence of small scale and large scale roughness of the ring structure itself. Because infrared emissivity cannot be derived without being deconvolved from the ‘structural’ filling factor, we examine temperature and filling factors measurements at opposition where the filling factor is most constrained. The occurrence of higher temperatures in the plateaus than in the background near opposition likely arises from enhanced mutual heating between particles, multiple scattering and surface roughness combined with a higher single-scattering albedo.     

Planetary brightness temperature measurements at 8.6 mm and 3.1 mm wavelengths

New measurements of the Sun, Moon, Mercury, Venus, Mars, Jupiter, and Saturn at 3.1 and 8.6 mm wavelengths are given. The temperatures reported for the planets at 3.1 mm wavelength are higher than previous measurements in this wavelength range and change the interpretation of some planetary spectra. For Mercury, it is found that the mean brightness temperature is independent of wavelength and that a temperature dependent thermal conductivity is not required to match the observations. In the case of Mars, the spectrum is shown to rise in the millimeter region as simple models predict. For Jupiter, the need to recalculate the spectrum with recent models is demonstrated. The flux density scale proposed by Dent (1972) has been revised according to a more accurate determination of the millimeter brightness temperature of Jupiter.     

A coupled microphysical/radiative transfer model of albedo and emissivity of planetary surfaces covered by volatile ices

A new model of albedo and emissivity of planetary surfaces covered by volatile ices in the form of porous slab-like deposits is described. In the model, a radiative transfer model is coupled with a microphysical model in order to link changes in albedo and emissivity to changes in porosity caused by ice metamorphism. Preliminary results for Triton, Pluto, and Io are presented (the martian CO₂ caps will be the subject of a separate publication). The coupled model will aid in the interpretation of ground-based and spacecraft observations and should lead to advances in surface and atmospheric modeling.     

First ENA observations at Mars: Solar-wind ENAs on the nightside

We present measurements with an Energetic Neutral Atom (ENA) imager on board Mars Express when the spacecraft moves into Mars eclipse. Solar wind ions charge exchange with the extended Mars exosphere to produce ENAs that can spread into the eclipse of Mars due to the ions' thermal spread. Our measurements show a lingering signal from the Sun direction for several minutes as the spacecraft moves into the eclipse. However, our ENA imager is also sensitive to UV photons and we compare the measurements to ENA simulations and a simplified model of UV scattering in the exosphere. Simulations and further comparisons with an electron spectrometer sensitive to photoelectrons generated when UV photons interact with the spacecraft suggest that what we are seeing in Mars' eclipse are ENAs from upstream of the bow shock produced in charge exchange with solar wind ions with a non-zero temperature. The measurements are a precursor to a new technique called ENA sounding to measure solar wind and planetary exosphere properties in the future.     

Evolution of density in solar system ICES

Pores present in ices in the solar system do not remain unchanged. In isothermal conditions they shrink while in a thermal gradient they migrate towards the higher temperature and escape so that the ice densifies. This motion has been investigated for pure H₂O- and CO₂-ices in a very simple one-dimensional model assuming uniform thermal conductivity and temperature gradient. The results indicate that the densification of H₂O-ice is so slow that it could be significant only for icy satellites having an internal heat source. On the other hand, CO₂-ice densifies orders of magnitude faster and the effect should be important for the CO₂ component of cometary nuclei. No effect is expected for icy planetary rings.     

The depth of the convective boundary layer on Mars

We are using observations obtained with Mars Express to explore the structure and dynamics of the martian lower atmosphere. We consider a series of radio occultation experiments conducted in May-August 2004, when the season on Mars was midspring of the northern hemisphere. The measurements are widely distributed in latitude and longitude, but the local time remained within a narrow range, 17.0-17.2 h. Most of the atmospheric profiles retrieved from these data contain a distinct, well-mixed convective boundary layer (CBL). We have accurately determined the depth of the CBL and its spatial variations at fixed local time through analysis of these profiles. The CBL extends to a height of 3-10 km above the surface at the season and locations of these measurements. Its depth at fixed local time is clearly correlated with variations in surface elevation on planetary scales, with a weaker dependence on spatial variations in surface temperature. In general, the CBL is deep (8-10 km) where the surface elevation is high, as in Tharsis Montes and Syrtis Major, and shallow (4-6 km) where the surface elevation is low, as in Amazonis and Utopia. This variability results from the combined effects of conditions near the surface and in the atmosphere above the CBL. Convection arises from solar heating of the ground, and the impact of this heat source on thermal structure is largest where the surface pressure and atmospheric density are smallest, at high surface elevations. The vertical extent of the CBL is in turn constrained by the static stability of the overlying atmosphere. These results greatly reduce the long-standing uncertainty concerning the depth of the CBL.     

Northwest Africa 6232: Visible–near infrared reflectance spectra variability of an olivine diogenite

Visible and near-infrared (VNIR) reflectance is an important spectroscopic technique to identify minerals, and their associations, on planetary body surfaces. Howardites, eucrites, and diogenites (HED) are a class of igneous-like meteorites whose genetic connection with asteroid 4 Vesta has since long been inferred and recently confirmed by Dawn mission results. Pyroxene and olivine are the two major mafic minerals present in HED which can be identified with VNIR reflectance measurements. Thus, studying the compositional variability of those phases and their mixtures by means of laboratory spectroscopic measurements on different diogenitic or eucritic samples is one of the prime methods to better understand the evolution of 4 Vesta's crust. Here, we report the VNIR reflectance spectral analysis of a harzburgitic olivine diogenite, Northwest Africa 6232 (probably paired with Northwest Africa 5480), containing variable amounts of olivine as small grains or aggregates. We found that the olivine diogenite spectral parameters (e.g., band position) of powdered samples and polished slabs are in agreement. Moreover, the olivine diogenite band position shifts from synthetic orthopyroxene in accordance with the presence of olivine and chromite. In particular, the presence of a large olivine clast permits us to determine a linear variation of the band position from synthetic orthopyroxene and olivine, but underestimates the presence of olivine in the olivine diogenite spot.     

85-μm fluxes from asteroids: 2 Pallas, 7 Iris, 15 Eunomia,and 45 Eugenia

Broad-bandpass flux measurements with effective wavelengths at 20, 27 and 85 μm are presented for four asteroids—2 Pallas, 7 Iris, 15 Eunomia, and 45 Eugenia. The observations were made during the Far Infrared Sky Survey Experiment (FIRSSE) with a super-fluid helium-cooled telescope flown to exospheric altitudes aboard a sounding rocket. The 20- and 27-μm fluxes are consistent with published 10- and 20-μm photometry. However, the 85-μm fluxes are a factor of 2 to 3 lower than a Planck law extrapolation of the shorter wavelength fluxes. If the thermal radiation from asteroids at these wavelenghts is adequately described by graybody emission, as is expected on the basis of models of asteroid radiation discussed in the literature, then the value of emissivity in the 85-μm bandpass is likely to be substantially lower than the 10- to 20-μm value. The decrease in emissivity between 20 and 85 μm suggested for these asteroids is in qualitative agreement with studies, both observational and theoretical, of emissivity of a lunar-type regolith. Fluxes measured at 20 and/or 27 μm are also reported for an additional 16 asteroids.     

Thermal behavior of horizontally mixed surfaces on Mars

Current methods for deriving thermal inertia from spacecraft observations of planetary brightness temperature generally assume that surface properties are uniform for any given observation or co-located set of observations. As a result of this assumption and the nonlinear relationship between temperature and thermal inertia, sub-pixel horizontal heterogeneity may yield different apparent thermal inertia at different times of day or seasons. We examine the effects of horizontal heterogeneity on Mars by modeling the thermal behavior of various idealized mixed surfaces containing differing proportions of either dust, sand, duricrust, and rock or slope facets at different angles and azimuths. Latitudinal effects on mixed-surface thermal behavior are also investigated. We find large (several 100 J m⁻² K⁻¹ s^−1/2) diurnal and seasonal variations in apparent thermal inertia even for small (∼10%) admixtures of materials with moderately contrasting thermal properties or slope angles. Together with similar results for layered surfaces [Mellon, M.T., Putzig, N.E., 2007. Lunar Planet. Sci. XXXVIII. Abstract 2184], this work shows that the effects of heterogeneity on the thermal behavior of the martian surface are substantial and may be expected to result in large variations in apparent thermal inertia as derived from spacecraft instruments. While our results caution against the over-interpretation of thermal inertia taken from median or average maps or derived from single temperature measurements, they also suggest the possibility of using a suite of apparent thermal inertia values derived from single observations over a range of times of day and seasons to constrain the heterogeneity of the martian surface.     

Thermal observations of Saturn's main rings by Cassini CIRS: Phase, emission and solar elevation dependence

Two and a half years after Saturn orbit insertion (SOI) the Cassini composite infrared spectrometer (CIRS) has acquired an extensive set of thermal measurements (including physical temperature and filling factor) of Saturn's main rings for a number of different viewing geometries, most of which are not available from Earth. Thermal mapping of both the lit and unlit faces of the rings is being performed within a multidimensional observation space that includes solar phase angle, spacecraft elevation and solar elevation. Comprehensive thermal mapping is a key requirement for detailed modeling of ring thermal properties.To first order, the largest temperature changes on the lit face of the rings are driven by variations in phase angle while differences in temperature with changing spacecraft elevation are a secondary effect. Ring temperatures decrease with increasing phase angle suggesting a population of slowly rotating ring particles [Spilker, L.J., Pilorz, S.H., Wallis, B.D., Pearl, J.C., Cuzzi, J.N., Brooks, S.M., Altobelli, N., Edgington, S.G., Showalter, M., Michael Flasar, F., Ferrari, C., Leyrat, C. 2006. Cassini thermal observations of Saturn's main rings: implications for particle rotation and vertical mixing. Planet. Space Sci. 54, 1167-1176, doi: 10.1016/j.pss.2006.05.033]. Both lit A and B rings show that temperature decreases with decreasing rings solar elevation while temperature changes in the C ring and Cassini Division are more muted. Variations in the geometrical filling factor, β, are primarily driven by changes in spacecraft elevation. For the optically thinnest region of the C ring, β variations are found to be nearly exclusively determined by spacecraft elevation. Both a multilayer and a monolayer model provide an excellent fit to the data in this region. In both cases, a ring infrared emissivity >0.9 is required, together with a random and homogeneous distribution of the particles. The interparticle shadowing function required for the monolayer model is very well constrained by our data and matches experimental measurements performed by Froidevaux [1981a. Saturn's rings: infrared brightness variation with solar elevation. Icarus 46, 4-17].