A laboratory study of the bidirectional reflectance from particulate samples

The Hapke (Hapke, B. [1981]. J. Geophys. Res. 86, 3039-3054) photometric model and its modifications are widely used to characterize telescopic, spacecraft, and laboratory observations of the bidirectional reflectance of particulate surfaces. Following work and methods laid out in a companion paper (Helfenstein, P., Shepard, M.K. [2011]. Icarus, in press), we deconstruct the Hapke model and, separating all empirical and ad hoc parameters (opposition surge, particle phase function, surface roughness), combine them into a single parameter called the surface phase function, F(α). We illustrate how to extract this function from scattering data sets acquired with the Bloomsburg University Goniometer (BUG). We show how this method can be used to rapidly and accurately characterize bidirectional reflectance data sets from laboratory and spacecraft measurements, often giving better fits to the data. We examine samples with strong color contrasts in different wavelengths. This allows us to examine the exact same surface, changing only the albedo to investigate how the amplitude and the detailed shape of the surface phase function might systematically depend on wavelength and albedo. We also examine the changes in scattering behavior that result when samples are compacted and find the surface phase function and single scattering albedo to be significantly changed. We suggest that these observations support the hypothesis that much of the scattering behavior attributed to the single particle phase function is instead cause by the surface micro-structure.     

Physical characteristics of Hayabusa target Asteroid 25143 Itokawa

In March 2001, the Hayabusa spacecraft target, Asteroid 25143 Itokawa, made its final close approach to Earth prior to the spacecraft's launch. We carried out an extensive observing campaign from January to September 2001 to better characterize this near-Earth asteroid. Global physical properties of the surface of Itokawa were characterized by analyzing its photometric properties and behavior. Results included here capitalize on analysis of broadband photometric observations taken with a number of telescopes, instruments, and observers. We employed a Hapke model to estimate the surface roughness, single particle scattering albedo, single particle scattering characteristics, phase integral, and geometric and bond albedo. We find that this asteroid has a higher geometric albedo than average main belt S-class asteroids; this is consistent with results from other observers. The broadband colors of Itokawa further support evidence that this is an atypical S-class asteroid. Broadband colors show spectral characteristics more typically found on large-diameter main-belt asteroids believed to be space-weathered, suggesting the surface of this small diameter, near-Earth asteroid could likewise be space-weathered.     

Gusev photometric variability as seen from orbit by HRSC/Mars-express

Minnaert and two-term phase function Hapke models are used to describe the photometric properties of the martian surface using HRSC (High Resolution Stereo Camera) multi-angular observations acquired along the ongoing Mars-Express mission. These observations can be pieced together to derive integrated phase functions over a wide range of phase angles. The photometric diversity at 675 nm, as seen from orbit, of the martian surface properties across Gusev is depicted with seven units. Three photometric units widespread across the flanks of Apollinaris Patera flank and the floor of Gusev Crater are identified as having high single scattering albedo with rather forward scattering properties, low to intermediate macroscopic roughness and porous or not compacted powdered surface state as indicated by the opposition parameters. Another unit has the highest single scattering albedo, the smoothest surface in terms of macroscopic roughness, associated with an extremely forward scattering behavior. The opposition parameters are consistent with the presence of transparent particles in the surface powder layer. The distribution of this unit appears quite intermittent across the crater and does not seem to indicate any relationship with a given morphological structure. It may correspond to sparse areas where the structure of the surface dust layer is the most preserved. The most pronounced photometric changes are observed in three units associated with the low-albedo features corresponding to dark wind streaks. These units have a low single scattering albedo, are the most backscattering surfaces across Gusev, have a high surface roughness and present variable surface states as shown by the opposition parameters estimates, consistent with the occurrence of large grains organized in more or less packed layers. Clear differences are seen among these units in terms of opposition effect. While one exhibits typical characteristics for the opposition effect, another appears more unusual in terms of lobe width and the third suggests the occurrence of a packed/compressed/narrow size distribution powder surface. The opposition effect thus appears to play a significant role suggesting that the surface state optical properties across Gusev are strongly influenced by the porosity and packing characteristics or grain size distribution of the upper layer of the martian soil. The mapping aspect of the investigation is quite useful to get a better sense of the meaning of the observed photometric variations. Indeed, the Hapke modeling suggests that surface organization (surface roughness, packing state) is more important than the simple physical characterization of the intrinsic optical properties of the constitutive particles. Given the overall spatial patterns derived from the photometric analysis, the variations, at least for the western and central part of Gusev Crater, are likely partly driven by the prevailing wind regimes, considered to be oriented north-northwest/south-southeast and disturbing the very upper surface layer. The present photometric results agree with independent investigations based on thermal inertia, reflectance spectroscopy, in situ photometric and microscopic imaging and support the idea of a thin layer of fine-grained dust, being stripped off in the low albedo units to reveal a dark basaltic substrate comprising coarse-grained materials.     

The roughness of the dark side of Iapetus from the 2004 to 2005 flyby

The roughness of a planetary surface offers clues to its past geologic history. We apply a surface roughness model developed by Buratti and Veverka (Buratti, B.J., Veverka, J. [1985]. Icarus 64, 320-328) to Cassini ISS data from the January 1st, 2005 flyby of Iapetus. This model uses the observed scattering behavior to provide a depth to radius factor q quantifying the size of idealized craters on the surface. Our findings indicate that the surface on the dark side is significantly smoother than the surfaces of other icy low-albedo saturnian satellites. We have found that the average depth to radius on the leading (dark) side is 0.084, corresponding to a Hapke mean slope angle of 6°. As compared to the 13-33° Hapke mean slope angle of other icy satellites (Buratti, B.J., and 10 colleagues [2008]. Icarus 193, 309-322), our results present a clearly different picture for the leading surface of Iapetus, suggesting that the dark deposit contributes to the decrease in macroscopic surface roughness of the leading side. Attempts were made to obtain an average depth to radius value for the trailing (bright) side; however the scans of the bright side from this flyby exhibited large variations in albedo, resulting in results that were physically unrealistic.     

Properties of the Hermean regolith: IV. Photometric parameters of Mercury and the Moon contrasted with Hapke modelling

A comparison of the photometric properties of Mercury and the Moon is performed, based on their integral phase curves and disk-resolved image data of Mercury obtained with the Swedish Vacuum Solar Telescope. Proper absolute calibration of integral V-band magnitude observations reveals that the near-side of the Moon is 10-15% brighter than average Mercury, and 0-5% brighter for the “bolometric” wavelength range 400-1000 nm. As shown, this is supported by recent estimates of their geometric albedos. Hapke photometric parameters of their surfaces are derived from identical approaches, allowing a contrasting study between their surface properties to be performed. Compared to the average near-side Moon, Mercury has a slightly lower single-scattering albedo, an opposition surge with smaller width and of marginally smaller amplitude, and a somewhat smoother surface with similar porosity. The width of the lobes of the single-particle scattering function are smaller for Mercury, and the backward scattering anisotropy is stronger. In terms of the double Henyey-Greenstein b-c parameter plot, the scattering properties of an average particle on Mercury is closer to the properties of lunar maria than highlands, indicating a higher density of internal scatterers than that of lunar particles. The photometric roughness of Mercury is well constrained by the recent study of Mallama et al. (2002, Icarus 155, 253-264) to a value of about 8°, suggesting that the surfaces sampled by the highest phase angle observations (Borealis, Susei, and Sobkou Planitia) are lunar mare-like in their textural properties. However, Mariner 10 disk brightness profiles obtained at intermediate phase angles indicate a surface roughness of about twice this value. The photometric parameters of the Moon are more difficult to constrain due to limited phase angle coverage, but the best Hapke fits are provided by rather small surface roughnesses. Better-calibrated, multiple-wavelength observations of the integral and disk-resolved brightnesses of both bodies, and obtained at higher phase angle values in the case of the Moon, are urgently needed to arrive at a more consistent picture of the contrasting light scattering properties of their surfaces.     

Re-analysis of previous laboratory phase curves: 1. Variations of the opposition effect morphology with the textural properties,and an application to planetary surfaces

Typical variations in the opposition effect morphology of laboratory samples at optical wavelengths are investigated to probe the role of the textural properties of the surface (roughness, porosity and grain size). A previously published dataset of 34 laboratory phase curves is re-analyzed and fit with several morphological models. The retrieved morphological parameters that characterize the opposition surge, amplitude, width and slope (A, HWHM and S respectively) are correlated to the single scattering albedo, the roughness, the porosity and the grain size of the samples. To test the universality of the laboratory samples’ trends, we use previously published phase curves of planetary surfaces, including the Moon, satellites and rings of the giant planets. The morphological parameters of the surge (A and HWHM) for planetary surfaces are found to have a non-monotonic variation with the single scattering albedo, similar to that observed in asteroids (Belskaya, I.N., Shevchenko, V.G. [2000]. Icarus 147, 94–105), which is unexplained so far. The morphological parameters of the surge (A and HWHM) for laboratory samples seem to exhibit the same non-monotonic variation with single scattering albedo. While the non-monotonic variation with albedo was already observed by Nelson et al. (Nelson, R.M., Hapke, B.W., Smythe, W.D., Hale, A.S., Piatek, J.L. [2004]. Planetary regolith microstructure: An unexpected opposition effect result. In: Mackwell, S., Stansbery, E. (Eds.), Proc. Lunar Sci. Conf. 35, p. 1089), we report here the same variation for the angular width.     

Photometry of rough planetary surfaces: The role of multiple scattering

Topographic features affect the scattering properties of planetary surfaces by casting shadows and altering the local incidence and emission angles. Measurements of this phenomenon were obtained on the Cornell goniometer for both high and low albedo surfaces. For the low albedo surface, the decrease in reflected radiation due to topography increases sharply with increasing phase angle, whereas for the high albedo sample the effects are approximately constant between phase angles of 30 and 70°. The observations are in good agreement with a theoretical model in the case of the dark surface. However, for the high albedo surface the model overestimates the effects by about a factor of 2, since it does not include the partial illumination of shadows by multiple scattering. For both high and low albedo surfaces, the effects of topography do not become significant until a phase angle of 30–40°.     

Light Scattering from Regolith: Intensity Versus Particle Size Behavior

Nature of the photometric phase curves of the regolith like surfaces (like those of the asteroids) are believed to be dependent on the single particle characteristics like particle size, shape, composition etc. and physical characteristics of the surface like porosity and roughness. Most of the phase curves have a rapid surge of intensity at small phase angles (typically below 5^°) known as opposition effect, followed by a linear less decreasing trend at larger phase angles. Average intensity of the linear region has been found to be mostly dependent on the average particle size and its composition, in many laboratory observations. Generally, it is difficult to explain the nature of light scattering by an ensemble of irregular shaped inhomogeneous particles with a theoretical model, just by studying the phase curves. In the present work, we have investigated whether the theoretically expected variation of the scattered light intensity (at a given phase angle) with the average particle size of the grains constituting regoliths, for a given material of the particle is in agreement with the experimental results or not? If yes, this can be a simpler but efficient way to study light scattering by regolith like surfaces. For theoretical analysis, Hapke formula has been used with Mie theory for single particle phase function, where we have neglected the influence of porosity and roughness presently. The data are also fitted with an empirical formula. It has been found that this empirical formula may also be used to estimate the unknown average particle size of a real regolith with known composition.     

NEAR Infrared Spectrometer Photometry of Asteroid 433 Eros

We present near-infrared spectrometer (NIS) observations (0.8 to 2.4 μm) of the S-type asteroid 433 Eros obtained by the NEAR Shoemaker spacecraft and report results of our Hapke photometric model analysis of data obtained at phase angles ranging from 1.2° to 111.0° and at spatial resolutions of 1.25×2.5 to 2.75×5.5 km/spectrum. Our Hapke model fits successfully to the NEAR spectroscopic data for systematic color variations that accompany changing viewing and illumination geometry. Model parameters imply a geometric albedo at 0.946 μm of 0.27±0.04, which corresponds to a geometric albedo at 0.550 μm of 0.25±0.05. We find that Eros exhibits phase reddening of up to 10% across the phase angle range of 0-100°. We observe a 10% increase in the 1-μm band depth at high phase angles. In contrast, we observe only a 5% increase in continuum slope from 1.486 to 2.363 μm and essentially no difference in the 2-μm band depth at higher phase angles. These contrasting phase effects imply that there are phase-dependent differences in the parametric measurements of 1- and 2-μm band areas, and in their ratio. The Hapke model fits suggest that Eros exhibits a weaker opposition surge than either 951 Gaspra or 243 Ida (the only other S-type asteroids for which we possess disk-resolved photometric observations). On average, we find that Eros at 0.946 μm has a higher geometric albedo and a higher single-scatter albedo than Gaspra or Ida at 0.56 μm; however, Eros's single-particle phase function asymmetry and average surface macroscopic roughness parameters are intermediate between Gaspra and Ida. Only two of the five Hapke model parameters exhibit a notable wavelength dependence: (1) The single-scatter albedo mimics the spectrum of Eros, and (2) there is a decrease in angular width of the opposition surge with increasing wavelength from 0.8 to 1.7 μm. Such opposition surge behavior is not adequately modeled with our shadow-hiding Hapke model, consistent with coherent backscattering phenomena near zero phase.     

Global photometric properties of Asteroid (4) Vesta observed with Dawn Framing Camera

Dawn spacecraft orbited Vesta for more than one year and collected a huge volume of multispectral, high-resolution data in the visible wavelengths with the Framing Camera. We present a detailed disk-integrated and disk-resolved photometric analysis using the Framing Camera images with the Minnaert model and the Hapke model, and report our results about the global photometric properties of Vesta. The photometric properties of Vesta show weak or no dependence on wavelengths, except for the albedo. At 554 nm, the global average geometric albedo of Vesta is 0.38 ± 0.04, and the Bond albedo range is 0.20 ± 0.02. The bolometric Bond albedo is 0.18 ± 0.01. The phase function of Vesta is similar to those of S-type asteroids. Vesta’s surface shows a single-peaked albedo distribution with a full-width-half-max ∼17% relative to the global average. This width is much smaller than the full range of albedos (from ∼0.55× to >2× global average) in localized bright and dark areas of a few tens of km in sizes, and is probably a consequence of significant regolith mixing on the global scale. Rheasilvia basin is ∼10% brighter than the global average. The phase reddening of Vesta measured from Dawn Framing Camera images is comparable or slightly stronger than that of Eros as measured by the Near Earth Asteroid Rendezvous mission, but weaker than previous measurements based on ground-based observations of Vesta and laboratory measurements of HED meteorites. The photometric behaviors of Vesta are best described by the Hapke model and the Akimov disk-function, when compared with the Minnaert model, Lommel–Seeliger model, and Lommel–Seeliger–Lambertian model. The traditional approach for photometric correction is validated for Vesta for >99% of its surface where reflectance is within ±30% of global average.     

Comparison of USGS and DLR topographic models of Comet Borrelly and photometric applications

Stereo analysis of images obtained during the 2001 flyby of Comet Borrelly by NASA's Deep Space 1 (DS1) probe allows us to quantify the shape and photometric behavior of the nucleus. The shape is complex, with planar facets corresponding to the dark, mottled regions of the surface whereas the bright, smooth regions are convexly curved. The photometric as well as textural differences between these regions can be explained in terms of topography (roughness) at and below the image resolution, without invoking significant variations in single-particle properties; the material on Borrelly's surface could be quite uniform. A statistical comparison of the digital elevation models (DEMs) produced from the three highest-resolution images independently at the USGS and DLR shows that their difference standard deviation is 120 m, consistent with a matching error of 0.20 pixel (similar to reported matching accuracies for many other stereo datasets). The DEMs also show some systematic differences attributable to manual versus automatic matching. Disk-resolved photometric modeling of the nucleus using the DEM shows that bright, smooth terrains on Borrelly are similar in roughness (Hapke roughness θ=20°) to C-type asteroid Mathilde but slightly brighter and more backscattering (single-scattering albedo w=0.056, Henyey-Greenstein phase parameter g=−0.32). The dark, mottled terrain is photometrically consistent with the same particles but with roughnesses as large as 60°. Intrinsically darker material is inconsistent with the phase behavior of these regions. Many local radiance variations are clearly related to topography, and others are consistent with a topographic explanation; one need not invoke albedo variations greater than a few tens of percent to explain the appearance of Borrelly.     

On Hapke photometric model predictions on reflectance of closely packed particulate surfaces

In a recent paper Hapke et al. (Hapke, B., Shepard, M., Nelson, R., Smythe, W., Piatek, J. [2009]. Icarus 199, 210-218) performed bi-directional reflectance measurements on closely-packed particulate surfaces of micrometer-sized particles and compared these with both the Hapke IMSA photometric model, and a numerical radiative transfer algorithm, the MDYZ (Mishchenko, M., Dlugach, J., Yanovitskij, E., Zakharova, N. [1999a]. J. Quant. Spectrosc. Radiat. Trans. 63, 409-432). To account for the effects of close packing, Hapke et al. applied a diffraction truncation scheme to remove the diffraction spike and supplied the renormalized single scattering phase function to the IMSA. They found that the IMSA prediction is a better match with measurement than that of MDYZ. In this work we demonstrate that the diffraction truncation procedure outlined by Hapke et al. contains an error. By following Hapke et al.’s intended truncation scheme, we have found that the IMSA model is not sufficiently anisotropic to describe the reflectance pattern of measurements on surface reflectance of closely packed large spherical particles.     

Photometric analysis of the nucleus of Comet 81P/Wild 2 from Stardust images

The disk-resolved flyby images of the nucleus of Comet 81P/Wild 2 collected by Stardust are used to perform a detailed study of the photometric properties of this cometary nucleus. A disk-integrated phase function from phase angle 11° to about 100° is measured and modeled. A phase slope of 0.0513 ± 0.0002 mag/deg is found, with a V-band absolute magnitude of 16.29 ± 0.02. Hapke’s photometric model yields a single-scattering albedo of 0.034, an asymmetry factor of phase function −0.53, a geometric albedo 0.059, and a V-band absolute magnitude of 16.03 ± 0.07. Disk-resolved photometric modeling from both the Hapke model and the Minnaert model results in 11% model RMS, indicating small photometric variations. The roughness parameter is modeled to be 27 ± 5° from limb-darkening profile. The modeled single-scattering albedo and asymmetry factor of the phase function are 0.038 ± 0.004 and −0.52 ± 0.04, respectively, consistent with those from disk-integrated phase function. The bulk photometric properties of the nucleus of Wild 2 are comparable with those of other cometary nuclei. The photometric variations on the surface of the nucleus of Wild 2 are at a level of or smaller than 15%, much smaller than those on the nucleus of Comet 19P/Borrelly and comparable or smaller than those on the nucleus of Comet 9P/Tempel 1. The similar photometric parameters of the nuclei of Wild 2, Tempel 1, and the non-source areas of fan jets on Borrelly may reflect the typical photometric properties of the weakly active surfaces on cometary nuclei.     

Scattering of light from particulate surfaces: I. A laboratory assessment of multiple-scattering effects

A convenient photometric function for many particulate surfaces is the generalization of the Lommel-Seeliger law derived by Hapke (1963) and Irvine (1966). This generalization accounts for the effects of mutual shadowing among particles, but still assumes that multiple scattering within the surface layer can be neglected—an assumption which is evidently valid for dark surfaces. We describe a series of laboratory measurements which test the range of validity of this basic assumption, and the applicability of the Hapke-Irvine photometric function, for particulate surfaces whose normal reflectances ranges from 0.04 to 1.04. We find that multiple-scattering effects can be neglected, and that the Hapke-Irvine function can be used, for particulate surfaces whose normal reflectance is about 0.3, or less. The function is definitely inapplicable to surfaces whose normal reflectance exceeds 0.4.     

Whole-disk spectrophotometric properties of Mercury: Synthesis of MESSENGER and ground-based observations

Disk-integrated and disk-resolved measurements of Mercury’s surface obtained by both the Mercury Dual Imaging System (MDIS) and the Mercury Atmospheric and Surface Composition Spectrometer (MASCS) onboard the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft were analyzed and compared with previous ground-based observations of Mercury at 11 wavelengths. The spectra show no definitive absorption features and display a red spectral slope (increasing reflectance with increasing wavelength) typical of space-weathered rocky surfaces. The MDIS spectra show evidence of phase reddening, which is not observed in the MASCS spectra. The MDIS spectra are commensurate with ground-based observations to within 10%, whereas the MASCS spectra display greater discrepancies with ground-based observations at near-infrared wavelengths. The derived photometric calibrations provide corrections within 10% for observations taken at phase angles less than ∼100°. The derived photometric properties are indicative of a more compact regolith than that of the lunar surface or of average S-type asteroids. The photometric roughness of the surface is also much smoother than the Moon’s. The calculated geometric albedo (reflectance at zero phase) is higher than lunar values. The lower reflectance of immature units on Mercury compared with immature units on the Moon, in conjunction with the higher geometric albedo, is indicative of more complicated grain structures within Mercury’s regolith.     

Sources of error in the photoclinometric determination of planetary topography: A reappraisal

The technique of photoclinometry has frequently been used to determine planetary topography without proper consideration of possible sources of error. Previous studies of error sources have been limited in extent and have overlooked the importance of factors such as atmospheric scattering and the choice of a surface photometric function. This paper adopts a thorough and more direct approach to error analysis, whereby known topography is compared with photoclinometric profiles derived from synthetic quantised reflectance scans.Instrumental and geometric sources of error are found to exert a minimal influence on profiles in practice, provided that sufficient care is taken in the selection of images and the extraction of scans from those images. Environmental factors — relating to the scattering properties of the surface and, if present, atmosphere — are far more important. It is found that a simple Lommel-Seeliger law is unlikely to be appropriate to the majority of planetary terrains, given its inability to model the effects of multiple scattering or unresolved macroscopic roughness. It is further demonstrated that a Minnaert function or combination of Lommel-Seeliger and Lambert laws may empirically compensate for the first of these phenomena but not the second; in this respect, Hapke's equation is a far superior model of surface optical properties. In the case of an atmosphere, the need to correct for scattering by aerosols or suspended dust becomes more acute as atmospheric opacity increases and as particle scattering becomes more forward-biased. To perform this correction, a model for the combined reflectance of surface and atmosphere must be used when deriving profiles.Two case studies — of a small impact crater on Triton and a dust-mantled basaltic lava flow on Mars - are presented here. Regarding the latter, the implications that errors in photoclinometric flow thickness measurements have for inferred lava rheology are examined. Conservative estimates of errors in yield strength and apparent viscosity easily exceed 100% when one of the simplest photometric models possible — a Lommel-Seeliger law — is used to derive a profile.In the light of these findings, strategies are suggested for improving the results obtained from photoclinometry in the future.     

Statistical biases and errors inherent in photoclinometric surface slope estimation with natural light

Photoclinometry is the most common method used to obtain high-resolution topographic maps of planetary terrain. We derive the likelihood function of photoclinometric surface slope from (1) the probability distribution of the measured photon count of natural sunlight through a Charge-Coupled Device (CCD) including uncertainty due to camera shot noise, camera read noise, small-scale albedo fluctuation and atmospheric haze, and (2) a photometric model relating photocount to surface orientation. We then use classical estimation theory to determine the theoretically exact biases and errors inherent in photoclinometric surface slope and show when they may be approximated by asymptotic expressions for sufficiently high sample size. We show how small-scale albedo variability often dominates biases and errors, which may become an order of magnitude larger than surface slopes when surface reflectance has a weak dependence on surface tilt. We provide bounds on the minimum possible error of any unbiased photoclinometric surface slope estimate, and compute the sample sizes necessary to constrain errors within desired design thresholds.     

Martian phase function: Modeling the visible to near-infrared surface photometric function using HST-WFPC2 data

Images of Mars in the visible to near-infrared acquired from 1996 to 2005 using the Hubble Space Telescope WFPC2 have been used to model the martian surface photometric function at 502, 673, 953, and 1042 nm. These data range in spatial resolution from 12 to 70 km/pixel at the sub-Earth point, and in phase angle coverage from 0.34° to 40.5°. The WFPC2 images have been calibrated to radiance factor or I/F and projected to a cylindrical map for coregistration and comparison to similarly mapped spacecraft data sets of albedo, topography, thermal inertia, composition, and geology. We modeled the observed I/F as a function of phase angle using Minnaert, Lambert, lunar-Lambert, and Hapke photometric functions for numerous regions of interest binned into albedo units defined by Viking and TES albedo maps, and thermal-inertia units defined by TES thermal-inertia maps. Visibly opaque water-ice clouds and data acquired under high dust opacity conditions were excluded from the analysis. Our modeling suggests that under average to low atmospheric dust opacity conditions and over this range of phase angles, the photometric properties of the martian surface at 502, 673, 953, and 1042 nm are best modeled by lunar-Lambert functions with parameters derived for three surface units defined by low, moderate, and high TES bolometric albedos.