Saturn's Rings: Azimuthal variations,phase curves,and radial profiles in four colors

Four-color photographic photometry of Saturn for the 1977–1979 apparitions has been analyzed to determine the dependence of ring brightness on wavelength, solar phase angle, ring particle orbital phase angle (azimuthal effect), declination of the Earth relative to the ring plane (tilt angle), and radial distance from Saturn. Azimuthal brightness variations up to ±20% relative to the ansae are clearly apparent for the maximum of ring A, but are not detectable for ring B or the outer portion of ring A. The shape of the intensity (I) versus orbital phase angle (θ) curve varies with ring tilt (B) and probably with wavelength, and shows 180° symmetry. As characterized by its slope near the ansae, this curve suggests that the azimuthal effect increases as B decreases from 26 to ≈11°. The phase curves l(α) for the ansae show very little dependence on ring tilt (26° > B > 6°), on wavelength, or on radial distance from Saturn; possibly the curves are somewhat steeper at the smallest tilt angles and for ring A relative to ring B. The radial profile of both rings becomes flatter with decreasing tilt angle and with decreasing wavelength. The latter effect is a natural result of the classical, many-particle-thick ring model.     

On the structure of Saturn's rings and the ‘real’ rotational period for the planet

An analysis of the Lowell Observatory photographic plates of Saturn gave the following results: (1) ring A and B show peculiar brightness distributions around the planet, from which we conclude that both are composed of particles in synchronous rotation. (2) The leading side of the particles in ring A is brighter than the trailing side by about 4%, which may indicate an interaction between such particles and the interplanetary medium. (3) Scans of the rings across the major axis show a small (～0.3″) region of enhanced brightness, from which we derive a value ofT _s =10^h13 _. ^m 8±5 _. ^m 4 for the actual planetary rotational period of Saturn. (4) In order to explain the synchronous rotation, the particles in ring A have to be at least 42 m in diameter.     

Saturn's rings, the Yarkovsky effects, and the Ring of Fire

Saturn's icy ring particles, with their low thermal conductivity, are almost ideal for the operation of the Yarkovsky effects (photon thrust due to temperature gradients across the ring particles). An extremely simple case of the Yarkovsky effects is examined here, in which orbital evolution is computed as though each particle travels around Saturn alone in a circular orbit, so that there are no collisions, shadowing, or irradiance from other particles; nor are resonances, tumbling, or micrometeoroid erosion considered. The orbital evolution for random spin orientations appears to be a competition between two effects: the seasonal Yarkovsky effect, which makes orbits contract, and the Yarkovsky-Schach effect, which makes orbits expand. There are values of the far infrared and visible particle albedos for which (working radially out from the planet) the along-track particle acceleration S is negative, then positive, and then negative again; the region for which S>0 is interpreted as a region where stable rings are possible. Typical timescales for centimeter-sized particles to travel half a Saturn radius are ¹⁰⁷-¹⁰⁸ yr. Collisions, shadowing, and resonances may lengthen the timescales, perhaps considerably. It is speculated here that the C ring may be depleted of particles because of the seasonal Yarkovsky effect, and small particles that are present in the C ring ultimately fall on Saturn, possibly creating a “Ring of Fire” as they enter the planet's atmosphere.     

Interferometric observations of Saturn and its rings at a wavelength of 3 3.71 cm

Interferometric observations of Saturn and its rings made at the Owens Valley Radio Observatory at a wavelength of 3.71 cm ar fit to models of the Saturn brightness structure. The models have allowed us to estimate the brightness temperatures and optical thicknesses of the A, B, and C rings as well as the brightness temperature of the planetary disk. The most accurate results are the ratios of the ring temperatures to the planet temperature of 0.030 ± 0.012, 0.050 ± 0.010, and 0.040 ± 0.014 for the A, B, and C rings, respectively. The best estimates of the ring optical thicknesses are τ_A = 0.2 ± 0.1, τ_B = 0.9 ± 0.2, and τ_C = 0.1 ± 0.1. The actual brightness temperatures, which are affected by the absolute calibration errors, are T_planet = 178 ± 8, T_A = 5.2 ± 2.0, T_B = 9.1 ± 1.8, and T_C = 7.1 ± 2.6°K. The particle single-scattering albedo that would be most consistent with the observations is slightly less than one, but probably greater than 0.95. The observations are consistent with particles which conservatively scatter the thermal emission from Saturn to the Earth and emit no thermal emission of their own. The 3.71-cm optical depths which we have estimated are very close to the visible wavelength optical depths. This similarity indicates that the ring particles must be at least a few centimeters in size, although we feel that the particles may well be much larger than this in view of the closeness of the visible and microwave optical depths. Particles which are nearly conservative scatterers at our wavelength and at least a few centimeters in size must be composed of a material which is either a very good reflector of microwaves or a very poor absorber of them. At this time, water ice seems to be the most likely candidate since it is a very poor absorber of microwaves and has been detected in the rings spectroscopically.     

Infrared radiometry of the rings of Saturn

Broad-band radiometry with a spatial resolution of 5 arc sec is presented of Saturn and its rings. The brightness temperature of the B ring is 96 ± 3°K at 20 μm and 91 ± 3°K at 11 μm. These values constrain the bolometric Bond albedo of the ring particles to be less than 0.6, thus requiring a phase integral of less than unity. From differences in the thermal emission of the ansae, I suggest that the leading side of the particles has higher albedo than the trailing side. A measured drop in temperature of the B ring following eclipse of 2.0 ± 0.5°K is consistent with radii for the ring particles of 2 cm or larger.     

Infrared brightness temperature of Saturn's rings based on the inhomogeneous multilayer assumption

Models of Saturn's rings based on the classical multilayer assumption have been studied in the infrared. Thermal energy balance of Saturn's rings is treated rigorously by solving the infrared radiative transfer equations. It was found that a homogeneous multilayer model is incompatible with the observed infrared brightness variation of the A and B rings, although it can fit that of the C ring. The alternative inhomogeneous multilayer model with dark particles within a bright haze of small icy particles is presented in order to satisfy the available infrared data of the A, B, and C rings. The results based on the inhomogeneous multilayer model may be summarized as follows: The observed infrared brightness data of the three rings are explained in terms of the different optical thickness without having significant differences in the ring-particle properties, such as albedo, spin rate, and sizes. But each ring contains a different amount of bright haze particles and their concentration within the rings depends on whether or not dark particles emit radiation mostly from one hemisphere (slow rotator and/or low thermal inertia). If a dark particle is an isothermal radiator, the possible ranges of A₁ and A₂ for all three rings are given by 0.9 ? A₁ ? 0.95 and 0.0 ? A₂ ? 0.15, where A₁ and A₂ are the bolometric bond albedos of a bright haze and a dark particle, respectively. The possible ranges of the optical thickness ratio X of the dark particle layer to the total ring layer for the rings A, B, and C are given by 0.65 ? X ? 0.75, 0.8 ? X ? 0.9, and 0.8 ? X ? 1.0, respectively. If a dark particle is a slow rotator, we obtain 0.9 ? A₁ ? 0.95 and 0.0 ? A₂ ? 0.4 for all three rings. The ranges of X for the rings A, B, and C are given by 0.35 ? X ? 0.7, 0.65 ? X ? 0.9, and 0.35 ? X ? 1.0, respectively. In this paper, for the first time, a consistent model is presented which is applicable to all three rings from the multilayer point of view.     

Saturn's rings: Infrared brightness variation with solar elevation

The variation in infrared equilibrium brightness temperature of Saturn's A, B, and C rings is modeled as a function of solar elevation B′ with respect to the ring plane. The basic model includes estimates of minimum and maximum interparticle shadowing in a monolayer approximation. Simple laboratory observations of random particle distributions at various illumination angles provide more realistic shadowing functions. Radiation balance calculations yield the physical (kinetic) temperature of particles in equilibrium with radiation from the Sun, Saturn, and neighboring particles. Infrared brightness temperatures as a function of B′ are then computed and compared to the available 20-μm data (Pioneer results are also briefly discussed). The A and B rings are well modeled by an optically thick monolayer, or equivalently, a flat sheet, radiating from one side only. This points to a temperature contrast between the two sides, possibly due to particles with low thermal inertia. Other existing models for the B ring are discussed. The good fit for the monolayer model does not rule out the possibility that the A and B rings are many particles thick. It could well be that a multilayer ring produces an infrared behavior (as a function of tilt angle) similar to that of a monolayer. The C ring brightness increases as B′ decreases. This contrast in behavior can be understood simply in terms of the low C ring optical depth and small amount of interparticle shadowing. High-albedo particles (A?0.5) can fit the C ring infrared data if they radiate mostly from one hemisphere due to slow rotation or low thermal inertia (or both). Alternatively, particles isothermal over their surface (owing to a rapid spin, high inertia, or small size), and significantly darker (A?0.3) than the A and B ring particles, can produce a similar brightness variation with ring inclination. In any case, the C ring particles have significantly hotter physical temperatures than the particles in the A and B rings, whether or not the rings form a monolayer.     

Observation and photometry of an outer ring of saturn

A faint outer ring (E ring), which lies outside the classical rings A, B, C, and F, has been detected out to eight Saturn radii. We first observed it on November 1, 1979, and thereby confirmed the 1966 observation by Feibelman. Our plates were taken with a coronographic design and are specially intended for photometry. They are directly scaled in reflectance by reference to the Saturn disk which is properly attenuated. Photometry of the edge-on ring E lineament shows a strong brightness increase at small phase angles, which is compatible with scattering by particles of several microns in radius. The excess reflectivity in blue compared to the B ring implies a significant contribution of small particles in the scattering process. The E ring shows brightness and radial gradient changes, with condensations, which differ between east and west limbs and are not always the same from night to night. The E ring is probably a flat structure with a condensation centered at a distance of 4 R_s, but without a simple axial symmetry. It is probably shaped by segments or lumps and may have streamerlike structures.     

Two-micron spectrophotometry of Uranus and its rings

Multiaperture K photometry and 2.0- to 2.5-μm spectrophotometry of Uranus and its ring system are presented. The photometric results are used, together with a previously published measurement, to set limits on the geometric albedos of Uranus and the rings at ～2.2 μm: (0.74 ± 0.02) × 10(su?4) ≤ p_K (Uranus) ≤ (1.5 ± 0.3) × 10^?4, and (2.7 ± 0.6) × 10^?2 ≤ p_K (rings) ≤ (3.4 ± 0.1) × 10^?2. Reflectance spectra of Uranus and Uranus plus rings show features in the planet's spectrum which are attributed to gaseous CH₄ absorption, and a 2.20-μm feature in the combined spectrum which may be due to the rings. This feature is tentatively identified with either the 2.26-μm absorption feature of NH₃ frost, or the 2.2-μm OH band exhibited by certain silicate minerals. The results of JHK photometry of Uranus' satellite, Ariel (U1), indicate that the infrared colors of this object are very similar to those of the satellites U2, U3, and U4.     

Brightness of Saturn's rings with decreasing solar elevation

Early ground-based and spacecraft observations suggested that the temperature of Saturn's main rings (A, B and C) varied with the solar elevation angle, B^′. Data from the composite infrared spectrometer (CIRS) on board Cassini, which has been in orbit around Saturn for more than five years, confirm this variation and have been used to derive the temperature of the main rings from a wide variety of geometries while B^′ varied from near −24^° to 0^° (Saturn's equinox).Still, an unresolved issue in fully explaining this variation relates to how the ring particles are organized and whether even a simple mono-layer or multi-layer approximation describes this best. We present a set of temperature data of the main rings of Saturn that cover the ∼23^°—range of B^′ angles obtained with CIRS at low (α∼30^°) and high (α≥120^°) phase angles. We focus on particular regions of each ring with a radial extent on their lit and unlit sides. In this broad range of B^′, the data show that the A, B and C rings’ temperatures vary as much as 29-38, 22-34 and 18-23 K, respectively. Interestingly the unlit sides of the rings show important temperature variations with the decrease of B^′ as well. We introduce a simple analytical model based on the well known Froidevaux monolayer approximation and use the ring particles’ albedo as the only free parameter in order to fit and analyze this data and estimate the ring particle's albedo. The model considers that every particle of the ring behaves as a black body and warms up due to the direct energy coming from the Sun as well as the solar energy reflected from the atmosphere of Saturn and on its neighboring particles. Two types of shadowing functions are used. One analytical that is used in the latter model in the case of the three rings and another, numerical, that is applied in the case of the C ring alone. The model lit side albedo values at low phase are 0.59, 0.50 and 0.35-0.38 for the A, B and C rings, respectively.     

A Study of Saturn's Ring Phase Curves from HST Observations

Solar phase curves between 0.3° and 6.0° and color ratios at wavelengths λ=0.336 μm and λ=0.555 μm for Saturn's rings are presented using recent Hubble Space Telescope observations. We test the hypothesis that the phase reddening of the rings is less due to collective properties of the ring particles than to the individual properties of the ring particles. We use a modified Drossart model, the Hapke model, and the Shkuratov model to model reddening by either intraparticle shadow-hiding on fractal and normal surfaces, multiple scattering, or some combination. The modified Drossart model (including only shadowing) failed to reproduce the data. The Hapke model gives fair fits, except for the color ratios. A detailed study of the opposition effect suggests that coherent backscattering is the principal cause of the opposition surge at very small phase angles. The shape of the phase curve and color ratios of each main ring regions are accurately represented by the Shkuratov model, which includes both a shadow-hiding effect and coherent backscatter enhancement. Our analysis demonstrates that in terms of particle roughness, the C ring particles are comparable to the Moon, but the Cassini division and especially the A and B ring particles are significantly rougher, suggesting lumpy particles such as often seen in models. Another conspicuous difference between ring regions is in the effective size d of regolith grains (d∼λ for the C ring particles, d∼1-10 μm for the other rings).     

Optical reflectance polarimetry of Saturn globe and rings I. Measurements on B ring

The light reflected by the Saturn ring B acquires a polarization which varies with phase angle, wavelength, and position on the rings. This polarization results from the combined effects of two components P₀ and T. The component P₀ keeps its azimuth always parallel to the scattering plane and its amount varies with phase angle and wavelength, as for direct reflection at the surface of solid bodies. The polarization T has its azimuth all around the ring either parallel or perpendicular to the radius vector; unexpectedly its amount varies with time, often significantly in a few days in ways which are not predictable.     

Saturn's rings. I. Optical thickness of rings A,B, D and structure of ring B

A photometric study of high-resolution (～0″.3) plates of Saturn taken at the Lowell Observatory in 1943 and 1945 is presented. N-S scans were taken over both the planet and rings. The excess brightness due to the planet seen through the rings is found by taking the difference between the central meridian (CM) scans and scans displaced by 5″.7. Adopting a value for the albedo of the planet, it is possible to obtain the optical thickness, τ_CM(r). In particular, for the regions of maximum brightness in rings A and B, we find τ_CM(I_{A max}) = 0.38 ± 0.11 and τ_CM(I_{B max}) = 0.61 ± 0.11. Observations by Barnard made in 1890 show evidence of ring D, recently discovered by Guerin (1969). The value for the optical thickness of this ring is τ_D(I_{D max}) = 0.03 ± 0.01. Ring B exhibits a pronounced (7–10%) decrease in brightness from the extremity of the major axis to the CM. After considering several possible explanations, we conclude that the ring particles are nonspherical and are in synchronous rotation around the planet with their long axis toward it. The mean value for the ratio of major to minor axis for the particles at 15″ is $\text{(a/b)}\text{? 1.08}$. Because of the shape and orientation of the particles, the optical thickness at the extremity of the major axis and at the CM are different for any saturnicentric latitude B ≠ 90°. Under these circumstances, only a minimum value for τ at the extremity can be derived.     

Aperture synthesis observations of Saturn and its rings at 2.7-mm wavelength

We present 2.7-mm interferometric observations of Saturn made near opposition in June 1984 and June 1985, when the ring opening angle was 19° and 23°, respectively. By combining the data sets we produce brightness maps of Saturn and its rings with a resolution of 6″. The maps show flux from the ring ansae, and are the first direct evidence of ring flux in the 3-mm wavelength region. Modelfits to the visibility data yield a disk brightness temperature of 156 ± 5°K, a combined A, B, and C ring brightness temperature of 19 ± 3°K, and a combined a ring cusp (region of the rings which block the planet's disk) brightness temperature of 85 ± 5°K. These results imply a normal-to-the-ring optical depth for the combined ABC ringof 0.31 ± 0.04, which is nearly the same value found for wavelenghts from the UV to 6 cm. About 6°K of the ring flux is attributed to scattered planetary emission, leaving an intrinsic thermal component of ∼13°K. These results, together with the ring particle size distributions found by the Voyager radio occultation experiments, are consistent with the idea that the ring particles are composed chiefly of water ice.     

A multilayer model for thermal infrared emission of Saturn’s rings II: Albedo, spins, and vertical mixing of ring particles inferred from Cassini CIRS

Since the Saturn orbit insertion of the Cassini spacecraft in mid-2004, the Cassini composite infrared spectrometer (CIRS) measured temperatures of Saturn’s main rings at various observational geometries. In the present study, we apply our new thermal model (Morishima, R., Salo, H., Ohtsuki, K. [2009]. Icarus 201, 634-654) for fitting to the early phase Cassini data (Spilker, L.J., and 11 colleagues [2006]. Planet. Space Sci. 54, 1167-1176). Our model is based on classical radiative transfer and takes into account the heat transport due to particle motion in the azimuthal and vertical directions. The model assumes a bimodal size distribution consisting of small fast rotators and large slow rotators. We estimated the bolometric Bond albedo, A_V, the fraction of fast rotators in cross section, f_fast, and the thermal inertia, Γ, by the data fitting at every radius from the inner C ring to the outer A ring. The albedo A_V is 0.1-0.4, 0.5-0.7, 0.4, 0.5 for the C ring, the B ring, the Cassini division, and the A ring, respectively. The fraction f_fast depends on the ratio of scale height of fast rotators to that of slow rotators, h_r. When h_r = 1, f_fast is roughly half for the entire rings, except for the A ring, where f_fast increases from 0.5 to 0.9 with increasing saturnocentric radius. When h_r increases from 1 to 3, f_fast decreases by 0.2-0.4 for the B and A rings while no change in f_fast is seen for the optically thin C ring and Cassini division. The large f_fast seen in the outer A ring probably indicates that a large number of small particles detach from large particles in high velocity collisions due to satellite perturbations or self-gravity wakes. The thermal inertia, Γ, is constrained from the efficiency of the vertical heat transport due to particle motion between the lit and unlit faces, and is coupled with the type of vertical motion. We found that in most regions, except for the mid B ring, sinusoidal vertical motion without bouncing is more reasonable than cycloidal motion assuming bouncing at the midplane, because the latter motion gives too large Γ as compared with previous estimations. For the mid B ring, where the optical depth is highest in Saturn’s rings, cycloidal vertical motion is more reasonable than sinusoidal vertical motion which gives too small Γ.     

The tilt effect for Saturn's rings

Multiple-scattering computations are carried out to explain the variation of the observed brightness of the A and B rings of Saturn with declination of the Earth and Sun. These computations are performed by a doubling scheme for a homogeneous plane-parallel scattering medium. We test a range of choices for the phase function, albedo for single scattering, and optical depth of both the rings. Isotropic scattering and several other simple phase functions are ruled out, and we find that the phase function must be moderately peaked in both the forward and backward directions. The tilt effect can be explained by multiple scattering in a homogeneous layer, but, for ring B, this requires a single-scattering albedo in excess of 0.8. The brightest part of ring B must have an optical depth greater than 0.9. We find that the tilt effect for ring A can be reproduced by particles having the same properties as those in ring B with the optical depth for the A ring in the range 0.4 to 0.6.     

Saturn's rings: 3-mm low-inclination observations and derived properties

Recent 3-mm observations of Saturn at low ring inclinations are combined with previous observations of E. E. Epstein, M. A. Janssen, J. N. Cuzzi, W. G. Fogarty, and J. Mottmann (Icarus41, 103–118) to determine a much more precise brightness temperature for Saturn's rings. Allowing for uncertainties in the optical depth and uniformity of the A and B rings and for ambiguities due to the C ring, but assuming the ring brightness to remain approximately constant with inclination, a mean brightness temperature for the A and B rings of 17 ± 4°K was determined. The portion of this brightness attributed to ring particle thermal emission is 11 ± 5°K. The disk temperature of Saturn without the rings would be 156 ± 6°K, relative to B. L. Ulich, J. H. Davis, P. J. Rhodes, and J. M. Hollis' (1980, IEEE Trans. Antennas Propag.AP-28, 367–376) absolutely calibrated disk temperature for Jupiter. Assuming that the ring particles are pure water ice, a simple slab emission model leads to an estimate of typical particle sizes of ≈0.3 m. A multiple-scattering model gives a ring particle effective isotropic single-scattering albedo of 0.85 ± 0.05. This albedo has been compared with theoretical Mie calculations of average albedo for various combinations of particle size distribution and refractive indices. If the maximum particle radius (≈5 m) deduced from Voyager bistatic radar observations (E. A. Marouf, G. L. Tyler, H. A. Zebker, V. R. Eshleman, 1983, Icarus54, 189–211) is correct, our results indicate either (a) a particle distribution between 1 cm and several meters radius of the form r^?s with 3.3 ? s ? 3.6, or (b) a material absorption coefficient between 3 and 10 times lower than that of pure water ice Ih at 85°K, or both. Merely decreasing the density of the ice Ih particles by increasing their porosity will not produce the observed particle albedo. The low ring brightness temperature allows an upper limit on the ring particle silicate content of ≈10% by mass if the rocky material is uniformly distributed; however, there could be considerably more silicate material if it is segregated from the icy material.     

Thermal emission from a multiple scattering model of Saturn's rings

Models of Saturn's B ring have been investigated which include the shadowing mechanism, realistic phase functions for the ring particles, and the effects of multiple scattering and a particle size dispersion. These models are based on the assumption that the rings form a layer many particles thick. A power law relation dn ∝ ?^?s is used for the size dispersion law of the ring particles, where dn is the number of particles with radii between ? and ? + d?. In the calculation of the infrared brightness temperature of the rings, the effect of mutual heating among the ring particles is considered quantitatively for the first time. The parameters of the polydisperse s = 2 model can be chosen to satisfy both optical (λ ? 1.1 μ) and infrared data, but the situation could be much clarified if a good phase curve for the rings were available in the red, if the ring brightness were known accurately for λ > 1 μ, and if it could be established whether the ring particles are rotating synchronously.     

Saturn's rings: Particle size distributions for thin layer models

The small physical thickness of Saturn's rings requires that radio occultation observations be interpreted using scattering models with limited amounts of multiple scatter. A new model in which the possible order of near-forward scatter is strictly limited allows for the small physical thickness, and can be used to relate Voyager 1 observations of 3.6-and 13-cm wavelength microwave scatter from Saturn's rings to the ring particle size distribution function n(a), for particles with radius 0.001 ≤ a ≤ 20 m. This limited-scatter model yields solutions for particle size distribution functions for eight regions in Saturn's rings, which exhibit approximately inverse-cubic power-law behavior, with large-size cutoffs in particle radius ranging from about 5 m in ring C to about 10 m in parts of ring A. The power-law index is about 3.1 in ring C, about 2.8 in the Cassini division, and increases systematically with radial location in ring A from 2.7 at 2.10R_s to slightly more than 3.0 at 2.24R_s. Corresponding mass densities are 32–43 kg/m² in ring C, 188 kg/m² in the Cassini division, and 244–344 kg/m² in ring A, under the assumption that the material density of the particles is 0.9 g/cm³. These values are a factor of 1 to 2 lower than first-order mass loading estimates derived from resonance phenomena. In view of the uncertainties in the measurements and in the linear density wave model, and the strong arguments for icy particles with specific gravity not greater than about 1, we interpret this discrepancy as being indicative of possible differences in the regions studied, or systematic errors in the interpretation of the scattering results, the density wave phenomena, or some combination of the above.     

On obtaining the forward phase functions of Saturn ring features from radio occultation observations

The near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equations whose solution by usual inversion techniques (i.e., least mean squares) is precluded by the near singularity of the forward transformation matrix. The instabilities are reduced by applying a combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix, which yields derived phase functions valid over the range of zero to about 12 mrad. These functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Multiple scattering of the signal is a significant effect, and the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical ring models, (a) the monolayer and (b) the classical discrete random slab, and the fraction of opacity in submeter particles for each model for specific ring features is estimated. Four representative regions of the rings approximately between 1.3 and 1.4R_s, 1.5 and 1.52R_s, 2.0 and 2.02R_s, and 2.08 and 2.16R_s have been studied in detail and single-scatter phase functions produced. Each of these features exhibits effective particle sizes in the range of 3–6 m radius. The approximate fractions of optical thickness due to the submeter particles in each of these regions are 0.58, 0.54, 0.23, and 0.0, respectively, for the many-particle-thick model, and 0.67, 0.67, 0.50, and 0.50, for the monolayer model.