Photometry of pluto in the last decade and before: evidence for volatile transport?

Photometric observations of Pluto in the BVR filter system were obtained in 1999 and in 1990-1993, and observations in the 0.89-μm methane absorption band were obtained in 2000. Our 1999 observations yield lightcurve amplitudes of 0.30 ± 0.01, 0.26 ± 0.01, and 0.21 ± 0.02 and geometric albedos of 0.44 ± 0.04, 0.52 ± 0.03, and 0.58 ± 0.02 in the B, V, and R filters, respectively. The low-albedo hemisphere of Pluto is slightly redder than the higher albedo hemisphere. A comparison of our results and those from previous epochs shows that the lightcurve of Pluto changes substantially through time. We developed a model that fully accounts for changes in the lightcurve caused by changes in the viewing geometry between the Earth, Pluto, and the Sun. We find that the observed changes in the amplitude of Pluto’s lightcurve can be explained by viewing geometry rather than by volatile transport. We also discovered a measurable decrease since 1992 of ∼0.03 magnitudes in the amplitude of Pluto’s lightcurve, as the model predicts. Pluto’s geometric albedo does not appear to be currently increasing, as our model predicts, although given the uncertainties in both the model and the measurements of geometric albedo, this result is not firm evidence for volatile transport. The maximum of methane-absorption lightcurve occurs near the minimum of the BVR lightcurves. This result suggests that methane is more abundant in the brightest regions of Pluto. Pluto’s phase coefficient exhibits a color dependence, ranging from 0.037 ± 0.01 in the B filter to 0.032 ± 0.01 in the R filter. Pluto’s phase curve is most like those of the bright, recently resurfaced satellites Triton and Europa. Although Pluto shows no strong evidence for volatile transport now (unlike Triton), it is important to continue to observe Pluto as it moves away from perihelion.     

Chemical and physical properties of the variegated Pluto and Charon surfaces

We present new photometric and spectroscopic observations of the Pluto–Charon system carried out at the VLT-ESO (Chile) with two 8-m telescopes equipped with the FORS2, ISAAC and SINFONI instruments. The spectra were obtained in the 0.6–2.45 μm range with a spectral resolution from 300 to 1500. The SINFONI data were obtained using adaptive optics, allowing a complete separation of the two bodies. We derive both objects’ magnitudes in the near infrared and convert them into albedo values. These first near infrared photometric data allow to adjust the different parts of Pluto’s spectrum, provided by the three instruments. We run spectral models in order to give chemical and physical constraints on the surface of Pluto and Charon. We discuss the dilution properties of the methane ice and its implications on Pluto’s surface. The heterogeneities of the pure and diluted methane ice on Pluto’s surface is also investigated. The high signal-to-noise level of the data and our analyses may support the presence of ethane ice on the surface of Pluto, which is one of the main products of the methane irradiation and photolysis. The analyses of the spectra of Charon suggest that the water ice is almost completely in its crystalline form and that the ammonia compound is hydrated on the surface of this satellite.     

The Resonance Region of Plutinos under the Perturbation of Outer Planets

Three resonances, the 3:2 exterior mean motion resonance with Neptune, Kozai resonance and 1:1 super resonance, are known to govern concurrently the stability of the motion of Pluto. We explore numerically the resonance zones in which the three resonance coexist. There might exist plutinos with relatively large masses in these zones. Considering that Pluto's perturbation is important to the long-term evolution of plutinos, the resonance zone is mainly explored in the mirror region of Pluto, which is a mirror image of the region around Pluto with respect to the invariant plane of the solar system. We find six resonance zones in the mirror region. The orbit elements at the centers of the six zones and the corresponding heliocentric distances, longitudes and latitudes on July 1, 2003 have been computed and listed for the reference of observation.     

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.     

A dust cloud around Pluto and Charon

We suggest that Pluto and Charon are immersed in a tenuous dust cloud. The cloud consists of ejecta from Pluto and—especially—Charon, released from their surfaces by impacts of micrometeoroids originating from Edgeworth-Kuiper belt objects. The motion of the ejected grains is dominated by the gravity of Pluto and Charon, which determines a pear-shape of the densest part of the cloud. While the production rates of escaping particles from both sides are comparable, the lifetimes of the Charon particles inside the Hill sphere of Pluto-Charon with respect to the Sun are much longer than of the Pluto ejecta, so that the cloud is composed predominantly of Charon grains. The dust cloud is dense enough to be detected with an in situ dust detector onboard a future space mission to Pluto. The cloud's maximum optical depth of τ≈3×10⁻¹¹ is, however, too low to allow remote sensing observations.     

Satellite searches at Pluto and Mars

We report the results of searches for outer satellites of Pluto and Mars, carried out with the Hale 5-m telescope in 1999 and 2001, respectively. No new satellites were found down to limiting magnitudes of mR=25.0 at Pluto and mV∼22 at Mars, corresponding to diameters of 35 and 1 km, respectively, for an assumed geometric albedo of 0.07. A faint trans-neptunian object, 1999 LB₃₇, was discovered in the Pluto fields; given the depth of our survey, discovery of one such object in the background Kuiper belt is in reasonable agreement with expectations.     

The mass ratio of Charon to Pluto from Hubble Space Telescope astrometry with the fine guidance sensors

The mass ratio of Charon to Pluto is a basic parameter describing the binary system and is necessary for determining the individual masses and densities of these two bodies. Previous measurements of the mass ratio have been made, but the solutions differ significantly (Null et al., 1993; Young et al., 1994; Null and Owen, 1996; Foust et al., 1997; Tholen and Buie, 1997). We present the first observations of Pluto and Charon with a well-calibrated astrometric instrument—the fine guidance sensors on the Hubble Space Telescope. We observed the motion of Pluto and Charon about the system barycenter over 4.4 days (69% of an orbital period) and determined the mass ratio to be 0.122±0.008 which implies a density of 1.8 to 2.1 g cm⁻³ for Pluto and 1.6 to 1.8 g cm⁻³ for Charon. The resulting rock-mass fractions for Pluto and Charon are higher than expected for bodies formed in the outer solar nebula, possibly indicating significant postaccretion loss of volatiles.     

Ejecta exchange and satellite color evolution in the Pluto system, with implications for KBOs and asteroids with satellites

In this Note, I present first-order scaling calculations to examine the efficacy of impacts by Kuiper Belt debris in causing regolith exchange between objects in the Pluto system. It is found that ejecta can escape Nix and Hydra with sufficient velocity to reach one another, as well as Charon, and even Pluto. The degree of ejecta exchanged between Nix and Hydra is sufficient to cover these bodies with much more material than is required for photometrically change. In specific, Nix and Hydra may have exchanged as up to 10s of meters of regolith, and may have covered Charon to depths up to 14 cm with their ejecta. Pluto is likely unaffected by most Nix and Hydra ejecta by virtue of a combination of dynamical shielding from Charon and Pluto's own annual atmospheric frost deposition cycle. As a result of ejecta exchange between Nix, Hydra, and Charon, these bodies are expected to evolve their colors, albedos, and other photometric properties to be self similar. These are testable predictions of this model, as is the prediction that Nix and Hydra will have diameters near 50 km, owing to having a Charon-like albedo induced by ejecta exchange. As I discuss, this ejecta exchange process can also be effective in many KBOs and asteroids with satellites, and may be the reason that very many KBO and asteroid satellite systems have like colors.     

Accurate analytical representation of Pluto modern ephemeris

An accurate development of the latest JPL’s numerical ephemeris of Pluto, DE421, to compact analytical series is done. Rectangular barycentric ICRF coordinates of Pluto from DE421 are approximated by compact Fourier series with a maximum error of 1.3 km over 1900–2050 (the entire time interval covered by the ephemeris). To calculate Pluto positions relative to the Sun, a development of rectangular heliocentric ICRF coordinates of the Solar System barycenter to Poisson series is additionally made. As a result, DE421 Pluto heliocentric positions by the new analytical series are represented to an accuracy of better than 5 km over 1900–2050.     

A search for ethane on Pluto and Triton

We present here a search for solid ethane, C₂H₆, on the surfaces of Pluto and Triton, based on near-infrared spectral observations in the H and K bands (1.4-2.45 μm) using the Very Large Telescope (VLT) and the United Kingdom Infrared Telescope (UKIRT). We model each surface using a radiative transfer model based on Hapke theory (Hapke, B. [1993]. Theory of Reflectance and Emittance Spectroscopy. Cambridge University Press, Cambridge, UK) with three basic models: without ethane, with pure ethane, and with ethane diluted in nitrogen. On Pluto we detect weak features near 2.27, 2.405, 2.457, and 2.461 μm that match the strongest features of pure ethane. An additional feature seen at 2.317 μm is shifted to longer wavelengths than ethane by at least 0.002 μm. The strength of the features seen in the models suggests that pure ethane is limited to no more than a few percent of the surface of Pluto. On Triton, features in the H band could potentially be explained by ethane diluted in N₂, however, the lack of corresponding features in the K band makes this unlikely (also noted by Quirico et al. (Quirico, E., Doute, S., Schmitt, B., de Bergh, C., Cruikshank, D.P., Owen, T.C., Geballe, T.R., Roush, T.L. [1999]. Icarus 139, 159-178)). While Cruikshank et al. (Cruikshank, D.P., Mason, R.E., Dalle Ore, C.M., Bernstein, M.P., Quirico, E., Mastrapa, R.M., Emery, J.P., Owen, T.C. [2006]. Bull. Am. Astron. Soc. 38, 518) find that the 2.406-μm feature on Triton could not be completely due to ¹³CO, our models show that it could not be accounted for entirely by ethane either. The multiple origin of this feature complicates constraints on the contribution of ethane for both bodies.     

Spatial and Compositional Constraints on Non-ice Components and H2O on Pluto's Surface

We present four new near-infrared spectra of Pluto, measured separately from its satellite Charon during four HST/NICMOS observations in 1998, timed to sample four evenly spaced longitudes on Pluto. Being free of contamination by telluric absorptions or by Charon light, the new data are particularly valuable for studies of Pluto's continuum absorption. Previous studies of the major volatile species indicate the existence of at least three distinct terrains on Pluto's surface: N₂-rich, CH₄-rich, and volatile-depleted. The new data provide evidence that each of these three terrains has distinct near-infrared continuum absorption features. CH₄-rich regions appear to show reddish continuum absorption through the near-infrared spectral range. N₂-rich regions have very little continuum absorption. Visually dark, volatile-depleted regions exhibit intermediate continuum albedos with a bluish continuum slope. By analogy with Triton, we expected that careful spectral modeling would reveal strong evidence for the existence of H₂O ice on Pluto's surface, but we found only very weak evidence for its existence in the volatile-depleted regions. These data require H₂O ice to play a much less prominent role on Pluto's surface than it does on Triton's.     

An investigation of Pluto’s troposphere using stellar occultation light curves and an atmospheric radiative-conductive-convective model

We use a radiative-conductive-convective model to assess the height of Pluto’s troposphere, as well as surface pressure and surface radius, from stellar occultation data from the years 1988, 2002, and 2006. The height of the troposphere, if it exists, is less than 1 km for all years analyzed. Pluto has at most a planetary boundary layer and not a troposphere. As in previous analyses of Pluto occultation light curves, we find that the surface pressure is increasing with time, assuming that latitude and longitude variations in Pluto’s atmosphere are negligible. The surface pressure is found to be slightly higher ( μbar in 1988,  μbar in 2002, and 18.5 ± 4.7 μbar in 2006) than in our previous analyses with the troposphere excluded. The surface radius is determined to be . Comparison of the minimum reduced chi-squared values between the best-fit radiative-conductive-convective (i.e., troposphere-included) model and best-fit radiative-conductive (i.e., troposphere-excluded) shows that the troposphere-included model is only a slightly better fit to the data for all 3 years. Uncertainties in the small-scale physical processes of Pluto’s lower atmosphere and consequently the functional form of the model troposphere lend more confidence to the troposphere-excluded results.     

An analysis of Pluto occultation light curves using an atmospheric radiative-conductive model

We use a radiative-conductive model to least-squares fit Pluto stellar occultation light curve data. This model predicts atmospheric temperature based on surface temperature, surface pressure, surface radius, and CH₄ and CO mixing ratios, from which model light curves are to be calculated. The model improves upon previous techniques for deriving Pluto’s atmospheric thermal structure from stellar occultation light curves by calculating temperature (as a function of height) caused by heating and cooling by species in Pluto’s atmosphere, instead of a general assumption that temperature follows a power law with height or some other idealized function. We are able to fit for model surface radius, surface pressure, and CH₄ mixing ratio with one of the 2006 datasets and for surface pressure and CH₄ mixing ratio for other datasets from the years 1988, 2002, 2006, and 2008. It was not possible to fit for CO mixing ratio and surface temperature because the light curves are not sensitive to these parameters. We determine that the model surface radius, under the assumption of a stratosphere only (i.e. no troposphere) model in radiative-conductive balance, is . The CH₄ mixing ratio results are more scattered with time and are in the range of 1.8-9.4 × 10⁻³. The surface pressure results show an increasing trend from 1988 to 2002, although it is not as dramatic as the factor of 2 from previous studies.     

Sublimation’s impact on temporal change of albedo dichotomy on Iapetus

Iapetus, one of the saturnian moons, has an extreme albedo contrast between the leading and trailing hemispheres. The origin of this albedo dichotomy has led to several hypotheses, however it remains controversial. To clarify the origin of the dichotomy, the key approach is to investigate the detailed distribution of the dark material. Recent studies of impact craters and surface temperature from Cassini spacecraft data implied that sublimation of H₂O ice can occur on Iapetus’ surface. This ice sublimation can change the albedo distribution on the moon with time.In this study, we evaluate the effect of ice sublimation and simulate the temporal change of surface albedo. We assume the dark material and the bright ice on the surface to be uniformly mixed with a certain volume fraction, and the initial albedo distribution to incorporate the dark material deposits on the surface. That is, the albedo at the apex is lowest and concentrically increases in a sinusoidal pattern. This situation simulates that dark materials existed around the Iapetus’ orbit billions of years ago, and the synchronously rotating Iapetus swept the material and then deposited it on its surface. The evolution of the surface albedo during 4.0 Gyr is simulated by estimating the surface temperature from the insolation energy on Iapetus including the effect of Saturn’s eccentricity and Iapetus’ obliquity precession, and evaluating the sublimation rate of H₂O ice from the Iapetus’ surface.As a result, we found that the distribution of the surface albedo changed dramatically after 4.0 Gyr of evolution. The sublimation has three important effects on the resultant surface albedo. First, the albedo in the leading hemisphere has significantly decreased to approach the minimum value. Second, the albedo distribution has been elongated along the equator. Third, the edge of the low albedo region has become clear. Considering the effect of ice sublimation, the current albedo distribution can be reconstructed from the sinusoidal albedo distribution, suggesting the apex-antapex cratering asymmetry as a candidate for the origin of the albedo dichotomy. From the model analysis, we obtained an important aspect that the depth of the turn-over layer where the darkening process proceeded for 4 Gyr should be an order of 10 cm, which is consistent with evaluation from the Cassini radar observations.     

Comparing Phoebe’s 2005 opposition surge in four visible light filters

We observed Phoebe for 13 nights over a period of 55 days before, during, and after the 2005 Saturn opposition with the New Mexico State University (NMSU) 1-m telescope at Apache Point Observatory (APO) in Sunspot, NM and characterized the width and magnitude of Phoebe’s opposition surge in BVRI filters. Our observations cover a phase angle range of 4.87° to 0.0509°. We use a Hapke reflectance model incorporating shadow hiding and coherent backscatter to investigate the wavelength dependence of Phoebe’s opposition surge. We find a significant opposition surge magnitude of 55-58% between phase angles of 5° and 0°. We find the strongest opposition surge for phase angles less than 2° in the I-band. The coherent backscatter angular width is on the order of 0.50°. We find Phoebe’s albedo to be spectrally flat within our error limits, with a B-band albedo of 0.0855 ± 0.0031, a V-band albedo of 0.0856 ± 0.0023, an R-band albedo of 0.0843 ± 0.0020, and an I-band albedo of 0.0839 ± 0.0023. We compare Phoebe’s albedo, color, and opposition surge magnitudes and slopes with those of other outer solar system bodies and find similarities to Centaurs, Nereid, Puck, and Comets 19P/Borrelly, 9P/Tempel 1, and 81P/Wild 2. We find that this comparison supports the idea that Phoebe originated in the Kuiper Belt. We also discuss the caveats of using results from a Hapke reflectance model to derive specific surface particle properties.     

Thermal tides on Pluto

Stellar occultations have shown that vertical profiles of density fluctuations in the atmosphere of Pluto typically show wave-like structure with an amplitude of a few percent and vertical wavelengths of a few kilometers. Here we calculate the tidal response of Pluto’s atmosphere to solar-induced sublimation “breathing” from N₂ frost patches. Solutions show global-scale wave-like density structure capable of explaining the observations. The atmospheric response is a combination of eastward and westward migrating tides, together with a zonally symmetric mode. Calculated vertical wavelengths and amplitudes are similar to observations.     

The effect of Nix and Hydra on the putative Pluto–Charon dust cloud

Impact-generated dust clouds around airless bodies have been observed or suggested to be present throughout the solar system, including around the Martian, Galilean and Saturnian satellites. Simulations have assessed Pluto and Charon as sources of a possible dust cloud or torus and found that such a cloud would be dominated by Charon-produced ejecta and would have an optical depth of τ≈10⁻¹¹. These simulations were conducted before the discovery of two additional, small satellites of Pluto, Nix and Hydra. These small moons may yield impact-generated dust in excess of their larger counterparts due to their lower escape velocities, despite their smaller cross sections. In this paper, we extend a previous model of the Pluto–Charon dust cloud to include Nix and Hydra, both as sinks for Pluto- and Charon-generated dust and as sources of impact-generated dust. We find that Nix- and Hydra-generated dust grains outlive Pluto and Charon dust grains significantly and are the dominant contributors of dust in the Pluto–Charon system. Furthermore, we estimate the net geometric optical depth of grains between 0.1 and to be on the order of 10⁻⁷.     

Does pluto have a substantial atmosphere?

The presence of CH₄ ice on Pluto implies that Pluto may have a substantial atmosphere consisting of heavy gases. Without such an atmosphere, sublimation of the CH₄ ice would be so rapid on a cosmogonic time scale that either such an atmosphere would soon develop through the exposure of gases trapped in the CH₄ ice or else the surface CH₄ ice would soon be all sublimated away as other, more stable, ices became exposed. If such stable ices were present from the beginning, the existence of CH₄ frosts would also imply that Pluto's present atmosphere contains a remnant of its primordial atmosphere.     

Martian albedo feature variations with season: Data of 1971 and 1973

This paper is a qualitative and quantitative investigation of recent seasonal and secular albedo feature variations on Mars using the improved red-filter photography of the International Planetary Patrol Program obtained during the 1971 and 1973 apparitons. The boundaries of telescopic albedo features showing significant changes and uncommonly bright areas were measured with an image projection reader and plotted on Mars maps in order to compare changes with Mariner 9 variable feature data obtained in 1972. The Daedalia-Claritas darkening (90–130°W, 25°S) of 1973 was the most significant change observed on Mars since the Laocoontis-Alcyonius darkening of the 1950 decade. The 1973 apparition presented six excellent examples of normal seasonal changes.Photographic photometry of 21 albedo features located in three longitude regions was systematically done on 120 original red-filter negatives, about half representing the Martian southern spring season and the other half the southern summer season. These two sets of photographs differed about 90° in L_s but were nearly equal in viewing geometry. They were identical in image scale, emulsion, color response, and quality. The data, consisting of 3960 density measurements, were analyzed in terms of albedo ratios and their fluctuations. Several types of diagnostic graphs are presented.Seasonal differences of the mean brightness ratios for most of the features were found to be consistent with the contrast changes shown on photographic Mars maps. Seasonal variations of brightnesses with latitude and region were indicated, but not in the strict classical sense of a wave-of-darkening. Most albedo features were more variable during southern spring than they were in summer. Most light albedo features were more active than dark ones during both seasons. There was a general loss of relative contrast in southern summer. No systematic dependence on solar phase angle was detected within the range (38° > i > 16°) spanned by both sets of photographs.