Near-infrared studies of the satellites of Saturn and Uranus

New JHK photometry and spectrometry (1.4–2.6 μm) are presented for Enceladus, Hyperion, Phoebe, Umbriel, Titania, and Oberon. From spectral signatures, mainly in the 2-μm region, water ice is verified on Enceladus and identified on Hyperion and the three Uranian satellites. The JHK photometry shows that Phoebe is different from all other satellites and asteroids observed thus far. The new photometry corroborates the earlier conclusion by Cruikshank et al. (1977) Astrophys. J217, 1006–1010] that the Uranian satellites, as a class, have overall surface reflectances different from other water-ice-covered satellites, and the reason for the difference remains unclear. The diameters and the masses of the Uranian satellites are reviewed in light of the probable high albedo representative of ice-covered surfaces and the new dynamical studies by Greenberg, 1975, Greenberg, 1976, Greenberg, 1978.     

Nitrate in the Greenland ice sheet in the years following the 1908 tunguska event

The Tunguska event on 30 June 1908 has been subjected to much speculation within different fields of research. Publication of the results of the 1961 expedition to the Tunguska area (Florensky, 1963) supports that a cometary impact caused the event. Based on this interpretation, calculations of the impactor energy release and explosion height have been reported by Ben-Menahem (1975), and velocity, mass, and density of the impactor by Petrov and Stulov (1975). Park (1978) and Turco et al., 1981, Turco et al., 1982, used these numbers to calculate a production of ca. 30 × 10⁶ tons of NO during atmospheric transit. This paper presents a high-resolution study of nitrate concentration in the Greenland ice sheet in ca. 10 years covering the Tunguska event. No signs of excess nitrate are found in three ice cores from two different sites in Greenland in the years following the Tunguska event. By comparing these results with results for other aerosols generally found in the ice, the lack of excess NO₃^? following the Tunguska event can be interpreted as indicating that the impactor nitrate production calculated by Park (1978) and Turco et al., 1981, Turco et al., 1982 are 1–2 orders of magnitude too high. To explain this it is suggested, from other lines of reasoning, that the impactor density determined by Petrov and Stulov (1975) probably is too low.     

Spatially resolved absolute spectrophotometry of Saturn: 3390 to 8080 Å

Absolute spectrophotometry of four regions on the visible disk of Saturn (north and south polar regions, equatorial band, south “temperate” region) from 3390 to 8080 Å is reported. Spectral resolution is 10 Å in the interval 3390–6055 Å, and 20 Å; aperture size is 1.92 arcsec. The explicit purpose of our observations was to provide ground-based photometric calibration for the Pioneer Saturn Imaging Photopolarimeter (IPP). We also compare our data with earlier spectrophotometric measurements of Saturn (R.L. Younkin and G. Munch, 1963,Mem. Soc. Roy. Sci. Liege7, 123–136; W.M. Irvine and A.P. Lane, 1971,Icarus16, 10–26; T.B. McCord, T.V. Johnson, and J.H. Elias, 1971,Astrophys. J.165, 413–424) and with the M. Podolak and R.EE. Danielson (1977)Icarus30, 479–492) parameterization of “Axel Dust.” The latter reproduces the broad features but not the details of the observed spectral reflectivity (I/F). We find that large depths of clear molecular hydrogen (>14 km-am in the temperate regions) are needed to match the observed upturn in reflectivity shortward of 3800 Å.     

The pole orientation of asteroid 433 Eros determined by photometric astrometry

Refinements to the pole-determination method photometric astrometry (PA) were completed in 1983 (R. C. Taylor and E. F. Tedesco, 1983, Icarus54, 13–22). A goal is to redo the pole analysis for every asteroid whose pole had been determined from earlier versions of PA: Previous PA poles are reviewed in this paper. Asteroid 433 Eros is in that collection and has redone. The result are prograde rotation; a sidereal period of 0.219588 ± 0.000005 day; and a north pole at 22° longitude, +9° latitude. The uncertainty of the pole is 10°. The pole position of Eros determined by C.D. Vesely (1971, In Physical Studies of Minor Planets (T. Gehrels, Ed.), pp. 133–140, NASA SP-267) and Dunlap (1976, Icarus28, 69–78), using earlier versions of photometric astrometry, were within 21 and 7°, respectively, of the present result.     

Neptune's rotation period: A correction and a speculation on the difference between photometric and spectroscopic results

An error in the Hayes and Belton (1977), Icarus32, 383–401) estimate of the rotation period of Neptune is corrected. If Neptune exhibits the same degree of limb darkening as Uranus near 4900 Å, the rotation period is 15.4 ± 3 hr. This value is compatible with a recent spectroscopic determination of Munch and Hippelein (1979) who find a period of 11.2_?1.2^+1.8 hr. However, if, as indirect evidence suggests, the law of darkening on Neptune at these wavelengths is less pronounced than on Uranus, then the above estimates may need to be lengthened by several hours. Recent photometric data are independently analyzed and are found to admit several possible periods, none of which can be confidently assumed to be correct. The period of Neptune most probably falls somewhere in the range 15–20 hr. The Hayes-Belton estimate of the period of Uranus is essentially unaffected by the above-mentioned error and remains at 24 ± 4 hr. All observers agree that the rotation period of Uranus is longer than that of Neptune.     

Interpretation of the 6818.9-Å methane feature observed on Jupiter,Saturn, and Uranus

High-resolution (0.1-Å) spectra of the 6818.9-Å methane feature obtained for Jupiter, Saturn, and Uranus by K. H. Baines, W. V. Schempp, and W. H. Smith ((1983). Icarus56, 534–542) are modeled using a doubling and adding code after J. H. Hansen ((1969). Astrophys. J.155, 565–573). The feature's rotational quantum number is estimated using the relatively homogeneous atmosphere of Saturn, with only J = 0 and J = 1 fitting the observational constraints. The aerosol content within Saturn's northern temperate region is shown to be substantially less than at the equator, indicating a haze only half as optically thick. Models of Jupiter's atmosphere are consistent with the rotational quantum-number assignment. Synthetic line profiles of the 6818.9-Å feature observed on Uranus reveal that a substantial haze exists at or above the methane condensation region with an optical depth eight times greater than previously reported. Seasonal effects are indicated. The methane column abundance is 5 ± 1 km-am. The mixing ratio of methane to hydrogen within the deep unsaturated region of the planet is 0.045 ± 0.025, based on an H₂ column abundance of 240 ± 60 km-am (W. H. Smith, W. Macy, and C. B. Pilcher (1980). Icarus43, 153–160), thus indicating that the methane comprises between one-sixth and one-half of the planet's mass. However, proper reevaluation of H₂ quadrupole features accounting for the haze reported here may significantly reduce the relative methane abundance.     

Infrared spectrum of Io, 2.8–5.2 μm

The reflectance spectrum of Io is presented from 2.8 to 5.2 μm, extending the earlier results of D. P. Cruikshank, T. J. Jones, and C. B. Pilcher (1978, Astrophys. J. 225, L89–L92), and demonstrating the full extent of the broad and deep spectral absorption between 3.5 and 4.8 μm. Laboratory spectra of nitrates and carborates diluted with sulfur do not satisfactorily reproduce the Io spectrum, but new information based on the recently discovered volcanic activity on the satellite lead to consideration of other classes of compounds as reported in a companion paper (F. P. Fanale, R. H. Brown, D. P. Cruikshank, and R. N. Clark, 1979, Nature280, 761–763).     

Variability of carbon monoxide in the mars atmosphere

In January of 1982 we measured a microwave spectrum of CO in the Martian atmosphere utilizing the rotational J = 1 → 2 transition of CO. We have analyzed data and reanalyzed the microwave spectra of R. K. Kakar, J. W. Waters, and W. J. Wilson, (Science196, 1090–1091, 1977, measured in 1975) and J. C. Good and F. P. Schloerb, (Icarus47, 166–172, 1981 measured in 1980) in order to constrain estimates of the temporal variability of CO abundance in the Martian atmosphere. Our values of CO column density from the data of Karar et al., Good and Schloerb, and our own are 1.7 ± 0.9 × 10²⁰, 3.0 ± 1.0 × 10²⁰, and 4.6 ± 2.0 × 10²⁰cm^?2, respectively. The most recent estimate of CO column density from the 1967 infrared spectra of J. Connes, P. Connes, and J.P. Maillard, (Atlas de Spectres Infarouges de Venus, Mars, Jupiter, et Saturne, Editions due Centre National de la Recherche Scientifique, Paris, 1969), is 2.0 ± 0.8 × 10²⁰ cm^?2 (L.D.G. Young and A.T. Young, Icarus30, 75–79, 1977). The large uncertainties given for the microwave measurements are due primarily to uncertainty in the difference between the continuum brightness temperature and atmospheric temperatures of Mars. We have accurately calculated the variation among the observations of the continuum (surface) brightness temperature of Mars, which is primaroly a function of the observed aspect of Mars. A more difficult problem to consider is variability of global atmospheric temperatures among the observations, particularly the effects of global dust storms and the ellipticity of the orbit of Mars. The large bars accompanying our estimates of CO column density from the three sets of microwave measurements are primarily caused by an assumed uncertainty of ±10°K in our atmospheric temperature model due to possible dust in the atmosphere. A qualitative consideration of seasonal variability of global atmospheric temperatures among the measurements suggests that there is not strong evidence for variability of the column abundance of CO on Mars, although variability of 0–100% over a time scale of several years is allowed by the data set. The implication for the variability of Mars O₂ is, crudely, a factor of two less. We found that the altitude distribution of CO in the atmosphere of Mars was not well constrained by any of the spectra, although our spectrum was marginally better fitted by an altitude increasing profile of CO mixing ratios.     

Application of a radiative transfer model to bright icy satellites

A radiative transfer model, derived largely from the work of B.W. Hapke (1981, J. Geophys. Res.86, 3039–3054) and J.D. Goguen (1981, Ph.D. thesis, Cornell University, Ithaca, N.Y.), is fit to Voyager imaging observations of Europa, Mimas, Enceladus, and Rhea. It is possible to place constraints on the single-scattering albedo, the porosity of the optically active upper regolith, the single-particle phase functions, and, in the cases of Europa and Mimas, the mean slope angle of macroscopic surface features. The texture of the surfaces of the Saturnian satellites appears to be similar to the Earth's moon. However, Europa is found to have a distinctly more compact regolith and a more forward-scattering single-particle phase function.     

Interpretation of the infrared polarization of Venus

A contradiction in the sulfuric acid cloud hypothesis of Venus, i.e., nondetection of 4.8 μm polarization by Landau (1975), is examined on the basis of the multiple scattering calculations for the cloud model of Hansen and Hovenier (1974) including an internal heat source. Results show that the polarized thermal component cannot depolarize the scattered sunlight, and therefore a large polarization of about 13% is expected at a phase angle of 110° and wavelength of 4.8 μm, in contrast with Landau's measurements. Our computations are, however, in agreement with the measurements by S. Sato et al. (in “Proceedings, 10th Lunar and Planetary Symposium,” pp. 179–182. Institute of Space and Aeronautical Science, University of Tokyo, July 11–13, 1977).     

The occulation of β Scorpii by Jupiter. IV. Diurnal temperature variations and the methane mixing ratio in the Jovian upper atmosphere

The nightime cooling of the Jovian atmosphere near the occulation level of 10¹⁴cm^?3 is calculated using the models of Strobel (1973) and Strobel and Smith (1973). The amount of cooling is found to depend on χ, the methane mixing ratio; μ the mean molecular weight; and the sunrise temperature. Using the range of sunrise (emersion) temperatures observed by Veverka et al. (1974), the overnight cooling is calculated to be 1.5–5.5°K, if reasonable assumptions are made for χ and μ. The argument may be reversed to show that the agreement in measured sunrise and sunset temperatures obtained by other observers of the β Sco occulation implies that χ cannot be significantly greater than the generally accepted value of 7 ×10^?4.     

An explication of the radiometric method for size and albedo determination

A new radiometric model for disk-integrated photometry of asteroids is presented. With empirical support from photometry of Mercury and the Moon, the model assumes that observed sunward beaming of the infrared emission is due to craters. In contrast to earlier theoretical studies of the lunar emission, the observable flux ratio between a cratered sphere and a smooth sphere is calculated for large ranges in wavelength, temperature, and phase angle. Revised diameters and albedos based on the crater model are given for 84 asteroids. The revised values are in good agreement with Morrison's (1977) radiometric results. It is shown that the systematic discrepancy between radiometric and polarimetric albedos (Zellner and Gradie, 1976) is probably a double-valued function of albedo. Some typical geometric albedos from this paper, Morrison (1977), and Zellner and Gradie (1976), respectively, are: Ceres (0.050 ± 0.005, 0.053 ± 0.004, 0.068), Vesta (0.235 ± 0.032, 0235 ± 0.11, 0.271), mean C type (0.031 ± 0.009, 0.035 ± 0.009, 0.061 ± 0.005), mean S type (0.117 ± 0.030, 0.136 ± 0.032, 0.181 ± 0.23), and mean M type (0.105 ± 0.037, 0.115 ± 0.033, 0.157 ± 0.079). Areas of disagreement between radiometry and polarimetry are underscored, and research to resolve them is suggested.     

A laboratory atlas of the 5ν1 NH3 absorption band at 6475 Å with applications to Jupiter and Saturn

The 5ν₁ absorption band of NH₃ is displayed from 6418 to 6550 Å. The total band intensity has been measured: S_B = 0.66 cm^?1m^?1amagat^?1. Line intensities and self-broadening coefficients have been measured for some of the prominent lines. Our line intensities are in good agreement with those of Rank et al. (1966), but are about 2 times greater than those of Mason (1970). The spectrum displayed was obtained photoelectrically at a pressure of 0.061 atm, and shows many more lines than the spectrum obtained by McBride and Nicholls (1972a) at a pressure of 0.39 atm. Therefore, our new measurements can provide the basis for making a more complete rotational analysis than those of McBride and Nicholls (1972a).Since the total band absorption has previously been measured by others on moderate resolution photoelectric scans of the spectra of Jupiter and Saturn, we can use the band intensity to derive the NH₃ abundance in the atmospheres of these two planets. The NH₃ abundances in a single vertical path obtained by this method are about 10m amagat for Jupiter and 2m amagat for Saturn. These results are in agreement with previous results obtained from higher resolution photographic spectra.     

Absolute spectrophotometry of Neptune: 3390 to 7800 Å

Absolute spectrophotometry of Neptune from 3390 to 7800 Å, with spectral resolution of 10 Å in the interval 3390–6055 and 20 Å in the interval 6055–7800 Å, is reported. The results are compared with filter photometry (Appleby, 1973; Wamsteker, 1973; Savage et al., 1980) and with synthetic spectra computed on the basis of a parameterization proposed by Podolak and Danielson (1977) for aerosol scattering and absorption. A CH₄/H₂ ratio of is derived for the convectively mixed part of Neptune's atmosphere, and constrains optical properties of hypothetical aerosol layers.     

Thermoluminescence of meteorites: Relationships with their K-Ar age and their shock and reheating history

The observations of G. F. Komovsky [Meteoritika21 (1961), 64–69] and A. Liener and J. Geiss (in Thermoluminescence of Geological Materials, Academic Press, New York, 1968), that the thermoluminescence (TL) sensitivity of meteorites correlates with their KAr age, have been confirmed using a suite of 22 ordinary chondrites. In order to interpret this observation, meteorite samples have been exposed to doses of α, β and γ radiation comparable with those experienced over the lifetime of the meteorites and given a dose of protons comparable to the total dose received from cosmic rays. There was no increase in TL sensitivity after these treatments, suggesting that, contrary to the ideas of earlier workers, the TL mechanism does not involve radiation damage. The TL sensitivity of meteorites is therefore time independent. On the other hand, samples of meteorite annealed in a furnace at temperatures between 450 and 1250°C for 1 hr suffered up to an order-of-magnitude decrease in TL sensitivity. Similarly, samples of meteorite artificially shocked to pressures of the order of 400 kbar suffered a comparable decrease in TL sensitivity. It is concluded that the correlation between TL sensitivity and K-Ar age is entirely a result of the low K-Ar age meteorites being shocked or reheated. Data on the thermal and mechanical histories of these meteorites, based on ⁴⁰Ar?³⁹Ar, metallographic, and X-ray diffraction studies, seem to be consistent with this finding.     

Candidates for astrometric or direct-imaging detection of extra-solar planets,and comments on the photometric method

Knowledge of a star's rotation period and ν sin i can be used to select stars that are seen pole-on, and thus are well suited to planetary searches by astrometric or direct-imaging means. A table of such stars is presented. This method is not suitable for discriminating equator-on systems and so cannot be used to select candidates for the photometric method of W. J. Borucki and A. L. Summers (1984, Icarus58, 121–134).     

Interferometry of Saturn and its rings at 1.30-cm wavelength

We present interferometric observations of Saturn and its ring system made at the Hat Creek Radio Astronomy Observatory at a wavelength of 1.30 cm. The data have been analyzed by both model-fitting and aperture synthesis techniques to determine the brightness temperature and optical thickness of the ring system and estimate the amount of planetary limb darkening. We find that the ring optical depth is close to that observed at visible wavelenghts, while the ring brightness temperature is only 7 ± 1°K. These observational constraints require the ring particles to be nearly conservative scatterers at this wavelength. A conservative lower limit to the single-scattering albedo of the particles at 1.30-cm wavelength is 0.95, and if their composition is assumed to be water ice, then this lower limit implies an upper limit of 2.4 m for the radius of a typical ring particle. The aperture synthesis maps show evidence for a small offset in the position of Saturn from that given in the American Ephemeris and Nautical Almanac. The direction and magnitude of this offset are consistent with that found from a similar analysis of 3.71-cm interferometric data which we have previously presented (F.P. Schloerb, D.O. Muhleman, and G.L. Berge, 1979b, Icarus39, 232–250). Limb darkening of the planetary disk has been estimated by solving for the best-fitting disk radius in the models. The best-fitting radius is 0.998 ± 0.004 times the nominal Saturn radius and indicates that the planet is not appreciably limb dark at 1.30 cm. Since our previous 3.71-cm data also indicated that the planet was not strongly limb dark (F.P. Schloerb, D. O. Muhleman, and G.L. Berge, 1979a, Icarus39, 214–230), we feel that the limb darkening is not strongly wavelength dependent between 1.30 and 3.71 cm. The difference between the best-fitting disk radii at 3.71 and 1.30 cm is +0.007 ± 0.007 times the nominal Saturn radius and suggests that the planet is more limb dark at 1.30 cm than at 3.71 cm. Models of the atmosphere which have NH₃ as the principal source of microwave opacity predict that the planet will be less limb dark at 1.30 cm. However, the magnitude of the effect predicted by the NH₃ models is ?0.009 and only marginally different from the observed value.     

Martian CO abundance from the J = 1 → 0 rotational transition: Evidence for temporal variations

Mars was observed in the CO (J = 1 → 0) 2.6-mm wavelength line between 29 March and 1 April, 1980. The data were analyzed using a model atmosphere based on Viking measurements. A least-squares fit of the model to the observed line profile yielded an average CO mixing ratio of (3.2 ± 1.1) × 10^?3. This value is four times larger than that obtained by L. D. Kaplan, J. Connes, and P. Connes, 1969 (Astrophys. J.157, L187-L195) from analysis of an infrared spectrum obtained in 1967 by J. Connes, P. Connes, and J. P. Maillard, 1969 (Atlas of Near Infrared Spectra of Venus, Mars, Jupiter, and Saturn, Centre National de la Recherche Scientifique, Paris). Models of the Martian atmospheric chemistry indicate that this implied temporal variation could easily exist and that it would be due primarily to variations in the abundance of H₂O.     

Chaotic behavior and the origin of the 31 Kirkwood gap

The sudden eccentricity increases discovered by J. Wisdom (Astron J.87, 577–593, 1982) are reproduced in numerical integrations of the planar-elliptic restricted three-body problem, verifying that this phenomenon is real. Maximum Lyapunov characteristic exponents for trajectories near the $\text{3}\text{1}$ commensurability are computed both with the mappings presented in Wisdom (1982) and by numerical integration of the planar-elliptic problem. In all cases the agreement is excellent, indicating that the mappings accurately reflect whether trajectories are chaotic or quasiperiodic. The mappings are used to trace out the chaotic zone near the $\text{3}\text{1}$ commensurability, both in the planar-elliptic problem and to a more limited extent in the three-dimensional elliptic problem. The outer boundary of the chaotic zone coincides with the boundary of the $\text{3}\text{1}$ Kirkwood gap in the actual distribution of asteroids within the errors of the asteroid orbital elements.     

Some consequences of meteoroid impacts on Saturn's rings

High-velocity impacts of interplanetary meteoroids on Saturn's rings are discussed. It is shown that the neutral gas emitted by impact vaporization may be responsible, to a large part, for the observed neutral ring atmosphere. Both the predicted neutral gas injection rate and the gas temperature (or kinetic energy) are compatible with the measurements (see Broadfoot, A. L., B. R. Sandel, D. E. Shemansky, J. B. Holberg, G. R. Smith, D. F. Strobel, J. C. McConnell, S. Kumar, D. M. Hunten, S. K. Atreya, T. M. Dohnahne, H. W. Moos, J. L. Bertaux, J. E. Blamont, R. B. Pomphrey, and S. Linik, Science212, 206–211, 1981). Heavy ejecta particles produce a particulate ring “halo”. The physical properties of this halo are calculated, and it appears to be identical with the tenous particle population discussed by Baum and Kreidl (1982). Erosion of Saturn's ring particles, the resulting mass balance, and regolith formation are estimated. This provides some constraints on surface properties and optical albedo.