Deep Space 1 photometry of the nucleus of Comet 19P/Borrelly

The NASA-JPL Deep Space 1 Mission (DS1) encountered the short-period Jupiter-family Comet 19P/Borrelly on September 22, 2001, about 8 days after perihelion. DS1's payload contained a remote-sensing package called MICAS (Miniature Integrated Camera Spectrometer) that included a 1024 square CCD and a near IR spectrometer with ∼12 nm resolution. Prior to its closest approach of 2171 km, the remote-sensing package on the spacecraft obtained 25 CCD images of the comet and 45 near-IR spectra (L. Soderblom et al., 2002, Science 296, 1087-1091). These images provided the first close-up view of a comet's nucleus sufficiently unobscured to perform quantitative photometric studies. At closest approach, corresponding to a resolution of 47 meters per pixel, the intensity of the coma was less than 1% of that of the nucleus. An unprecedented range of high solar phase angles (52-89 degrees), viewing geometries that are in general attainable only when a comet is active, enabled the first quantitative and disk resolved modeling of surface photometric physical parameters, including the single particle phase function and macroscopic roughness. The disk-integrated geometric albedo of Borrelly's nucleus is 0.029±0.006, comparable to the dark hemisphere of Iapetus, the lowest albedo C-type asteroids, and the uranian rings. The Bond albedo, 0.009±0.002, is lower than that of any Solar System object measured. Such a low value may enhance the heating of the nucleus and sublimation of volatiles, which in turn causes the albedo to decrease even further. A map of normal reflectance of Borrelly shows variations far greater than those seen on asteroids. The two main terrain types, smooth and mottled, exhibit mean normal reflectances of 0.03 and 0.022. The physical photometric parameters of Borrelly's nucleus are typical of other small dark bodies, particularly asteroids, except preliminary modeling results indicate its regolith may be substantially fluffier. The nucleus exhibits significant variations in macroscopic roughness, with the oldest, darkest terrain being slightly smoother. This result suggests the infilling of low-lying areas with dust and particles that have not been able to leave the comet. The surface of the comet is backscattering, but there are significant variations in the single particle phase function. One region exhibits a flat particle phase function between solar phase angles of 50° and 75° (like cometary dust and unlike planetary surfaces), suggesting that its regolith is controlled by native dust rather than by meteoritic bombardment.     

The nucleus of Comet Borrelly: a study of morphology and surface brightness

Stereo images obtained during the DS1 flyby were analyzed to derive a topographic model for the nucleus of Comet 19P/Borrelly for morphologic and photometric studies. The elongated nucleus has an overall concave shape, resembling a peanut, with the lower end tilted towards the camera. The bimodal character of surface-slopes and curvatures support the idea that the nucleus is a gravitational aggregate, consisting of two fragments in contact. Our photometric modeling suggests that topographic shading effects on Borrelly's surface are very minor (<10%) at the given resolution of the terrain model. Instead, albedo effects are thought to dominate Borrelly's large variations in surface brightness. With 90% of the visible surface having single scattering albedos between 0.008 and 0.024, Borrelly is confirmed to be among the darkest of the known Solar System objects. Photometrically corrected images emphasize that the nucleus has distinct, contiguous terrains covered with either bright or dark, smooth or mottled materials. Also, mapping of the changes in surface brightness with phase angle suggests that terrain roughness at subpixel scale is not uniform over the nucleus. High surface roughness is noted in particular near the transition between the upper and lower end of the nucleus, as well as near the presumed source region of Borrelly's main jets. Borrelly's surface is complex and characterized by distinct types of materials that have different compositional and/or physical properties.     

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.     

Photometric analysis and disk-resolved thermal modeling of Comet 19P/Borrelly from Deep Space 1 data

Images returned from the Deep Space 1 (DS-1) spacecraft during its encounter with Comet 19P/Borrelly are used to study its disk-integrated and disk-resolved photometry and its thermal properties. A disk-integrated phase function was constructed from a combination of DS-1 images and ground-based observations, giving a geometric albedo of 0.072±0.020 and a phase slope of 0.043 mag deg⁻¹. The shape model of Borrelly [Kirk, R.L., Howington-Kraus, E., Soderblom, L.A., Giese, B., Oberst, J., 2004a. Icarus 167, 54-69] and the ephemerides of DS-1 were used to analyze the disk-resolved photometric data with Hapke's theoretical model. It was found that the surface of Borrelly displays large photometric heterogeneities in its photometric parameters. The single-scattering albedo, w, varies by a factor of 2.5 with an average of 0.057±0.009; the asymmetry factor, g, ranges from almost isotropic (−0.1) to strongly backscattering (−0.7) with an average of −0.43±0.07; the roughness parameter, , is less than 35° for most parts of surface but ranges up to 55° in some areas. Its average is 22°±5°. The observed 1-D temperature profile is modeled well by the standard thermal model (STM) for inactive regions and is found to be consistent with a very low thermal inertia. Water sublimation in the source region of the fan jet is observed to decrease the surface temperature from the STM predictions by 20-40 K. The source areas of two collimated jets could not be determined from either photometric model or thermal model. It is evident that the fan jet activity occurring on Borrelly's surface can be correlated to areas of relatively high albedo, weak backscattering, and high roughness.     

Are the circular, dark features on Comet Borrelly's surface albedo variations or pits?

The highest resolution images of Comet 19P/Borrelly show many dark features which, upon casual inspection, appear to be low albedo markings, but which may also be shadows or other photometric variations caused by a depression in the local topography. In order to distinguish between these two possible interpretations we conducted a photometric analysis of three of the most prominent of these features using six of the highest quality images from the September 22, 2001 Deep Space 1 (DS1) flyby. We find that: 1. The radiance in the darkest parts of each feature increases as phase angle decreases, similarly to the radiance behavior of the higher albedo surrounding terrain. The dark features could be either fully illuminated low albedo spots or, alternatively, they could be depressions. No part of any of the three regions was in full shadow. 2. One of the regions has a radiance profile consistent with a rimmed depression, the second, with a simple depression with no rim, and the third with a low albedo spot. 3. The regolith particles are backscattering and carbon black is one of the few candidate regolith materials that might explain this low albedo. We conclude that Borrelly's surface is geologically complex to the limit of resolution of the images with a combination complex topography, pits, troughs, peaks and ridges, and some very dark albedo markings, perhaps a factor of two to three darker than the average 3-4% albedo of the surrounding terrains. Our technique utilizing measured radiance profiles through the dark regions is able to discriminate between rimmed depressions, rimless depressions and simple albedo changes not associated with topography.     

Deep Impact photometry of Comet 9P/Tempel 1

The photometric properties of the nucleus of Comet 9P/Tempel 1 are studied from the disk-resolved color images obtained by Deep Impact (DI). Comet Tempel 1 has typical photometric properties for comets and dark asteroids. The disk-integrated spectrum of the nucleus of Tempel 1 between 309 and 950 nm is linear without any features at the spectral resolution of the filtered images. At V-band, the red slope of the nucleus is 12.5±1% per 100 nm at 63° phase angle, translating to B-V=0.84±0.01, V-R=0.50±0.01, and R-I=0.49±0.02. No phase reddening is confirmed. The phase function of the nucleus of Tempel 1 is constructed from DI images and earlier ground-based observations found from the literature. The phase coefficient is determined to be β=0.046±0.007 mag/deg between 4° and 117° phase angle. Hapke's theoretical scattering model was used to model the photometric properties of this comet. Assuming a single Henyey-Greenstein function for the single-particle phase function, the asymmetry factor of Tempel 1 was fitted to be g=−0.49±0.02, and the corresponding single-scattering albedo (SSA) was modeled to be 0.039±0.005 at 550 nm wavelength. The SSA spectrum shows a similar linear slope to that of the disk-integrated spectrum. The roughness parameter is found to be 16°±8°, and independent of wavelength. The Minnaert k parameter is modeled to be 0.680±0.014. The photometric variations on Tempel 1 are relatively small compared to other comets and asteroids, with a ∼20% full width at half maximum of albedo variation histogram, and ∼3% for color. Roughness variations are evident in one small area, with a roughness parameter about twice the average and appearing to correlate with the complex morphological texture seen in high-resolution images.     

The two faces of Iapetus

Combined photometry and radiometry of Iapetus can be used to investigate the nature of its surface and, in particular, the distribution of albedo that is responsible for the large variations in its visible and infrared brightness as it rotates. We present new 20-μm radiometric observations made in 1971–1973 and discuss these together with the photometric studies by Widorn (in 1949), Mills (in 1971), Noland et al. (in 1972–1973), and Franklin and Cook (in 1972–1974). The linear phase coefficient varies as the satellite rotates from 0.028 to 0.068 mag deg^?1. When corrected for this effect, the photometric variations suggest an albedo distribution characterized by a dark area covering most of the leading hemisphere and a bright trailing hemisphere and bright south polar cap. A combined analysis of the photometry and radiometry yields a radius of 800 to 850 km and mean geometric albedos for the light and dark faces of about 0.35 and 0.07, respectively. The average phase integral of the bright hemisphere is between 1.0 and 1.5. We offer no explanation for the unique photometric properties of this satellite.     

Interpreting photometry of regolith-like surfaces with different topographies: shadowing and multiple scattering

The effects of various types of topography on the shadow-hiding effect and multiple scattering in particulate surfaces are studied. Two bounding cases were examined: (1) the characteristic scale of the topography is much larger than the surface particle size, and (2) the characteristic scale of the topography is comparable to the surface particle size. A Monte Carlo ray-tracing method (i.e., geometric optics approximation) was used to simulate light scattering. The computer modeling shows that rocky topographies generated by randomly distributed stones over a flat surface reveal much steeper phase curves than surface with random topography generated from Gaussian statistics of heights and slopes. This is because rocks may have surface slopes greater than 90°. Consideration of rocky topography is important for interpreting rover observations. We show the roughness parameter in the Hapke model to be slightly underestimated for bright planetary surfaces, as the model neglects multiple scattering on large-scale topographies. The multiple scattering effect also explains the weak spectral dependences of the roughness parameter in Hapke's model found by some authors. Multiple scattering between different parts of a rough surface suppresses the effect of shadowing, thus the effects produced by increases in albedo on the photometric behavior of a surface can be compensated for with the proper decreases in surface roughness. This defines an effective (photometric) roughness for a surface. The interchangeability of albedo and roughness is shown to be possible with fairly high accuracy for large-scale random topography. For planetary surfaces that have a hierarchically arranged large-scale random topography, predictions made with the Hapke model can significantly differ from real values of roughness. Particulate media with surface borders complicated by Gaussian or clumpy random topographies with characteristic scale comparable to the particle size reveal different photometric behaviors in comparison with particulate surfaces that are flat or the scale of their topographies is much larger than the particle size.     

Photometric properties of lunar terrains derived from Hapke's equation

The photometric parameters of B. W. Hapke's (1986, Icarus 67, 264–280) equation are fit to the lunar disk-integrated visual lightcurve and to disk-resolved data of R. W. Shorthill, J. M. Saari, F. E. Baird, and J. R. LeCompte (1969, Photometric Properties of Selected Lunar Features, NASA Contractor Report CR-1429) for dark, average, and bright lunar terrains. The lunar nearside geometric albedo and phase integral computed from the disk-integrated results are consistent with those of earlier investigators. The single scattering albedos of disk-resolved average and bright lunar terrains are systematically larger than that of lunar mare. Average particles in dark terrain regoliths are more backscattering than those in average and bright lunar terrains. The angular width of the opposition surge is greatest in dark terrains and is found to be best explained by modest differences in regolith particle-size distributions which might accompany the normal regolith maturation process. The total amplitude of the opposition surge for dark terrains is larger than for average and bright terrains. This result appears to be a consequence of the fact that in opaque particles, a larger fraction of singly scattered light at zero phase comes from first-surface reflection. The average subcentimeter macroscopic roughness of dark terrains is significantly lower than that of average and bright terrains. The relative magnitude of this difference is consistent with that obtained from radar measurements at decimeter scales.     

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.     

Regolith transport in craters on Eros

Images of Eros from the NEAR Shoemaker spacecraft reveal bright and dark albedo features on steep crater walls unlike markings previously observed on asteroids. These features have been attributed to the downslope movement of space-weathered regolith, exposing less weathered material (Science 292 (2001) 484; Meteor. Planet. Sci. 36 (2001) 1617; Icarus 155 (2002) 145). Here we present observations of the interiors of large craters (>1 km in diameter) to test this hypothesis and constrain the origin of the features. We find that bright regions in these craters correspond to steep slopes, consistent with previous work. The geographic distribution of craters with albedo variations shows no pattern and does not resemble the distribution of ponds, another phenomenon on Eros attributed to regolith movement. Shadows and other indications of topography are not observed at feature boundaries, implying that the transported layer is ?1 m thick. The presence of multiple bright and dark units on long slopes with sharp boundaries between them suggests that mobilized regolith may be halted by frictional or other effects before reaching the foot of the slope. Features on crater walls should darken at the same rate as bright ejecta deposits from crater formation; the lack of observed, morphologically fresh craters with bright interiors or ejecta suggests that the albedo patterns are younger than the most recently formed craters greater than about 100 m in diameter. Smaller or micrometeorite impacts, which would not necessarily leave evident deposits of bright ejecta, remain possible causes of albedo patterns. Although their effectiveness is difficult to assess, electrostatic processes and thermal creep are also candidates.     

Variable features on Mars III: Comparison of Mariner 1969 and Mariner 1971 photography

Mariner 9 (M9) and Mariner 6 and 7 photography of common regions of Mars are compared, with appropriate attention to the photometric properties of the camera systems. The comparison provides a 2.5yr time baseline for study of variable albedo features. We find the development of bright streaks and patches, a phenomenon unobserved through the entire M9 mission; the evolution of dark crater splotches into dark streaks; and a planetwide increase in splotchiness. Yet, a large number of splotches and albedo boundaries remain fixed over the same period. Many of the observations are interpreted in terms of a global fallout and subsequent local redistribution of bright fine particulates raised by global dust storms.     

Time variability of Martian bolometric albedo

A time-sequential set of bolometric albedo maps for Mars has been constructed from Viking Infrared Thermal Mapper data. The maps provide global coverage in longitude for latitudes -60° to +60° at 1° by 1° spatial resolution. Individual maps are constructed under strict geometric constraints for a narrow range of L_s. The set of albedo maps spans a Martian year and includes maps before, during, and after the global dust storms of 1977. Transient brightenings associated with local dust storms or condensate clouds are apparent in some of the maps. During dust-free periods, bolometric albedo maps are generally similar to classical, visual albedo map of Mars. The distribution of bolometric surface albedos is bimodal with typical, clear-sky, Lambert albedos of 0.27 and 0.16 for bright and dark areas, respectively. Atmospheric effects strongly influence apparent surface albedos, especially for dark areas. Neither bright nor dark regions show measurable, long-term variations of bolometric albedos during clear periods.     

Lunar photometric properties at wavelengths 0.5-1.6 μm acquired by SELENE Spectral Profiler and their dependency on local albedo and latitudinal zones

The lunar photometric function, which describes the dependency of the observed radiance on the observation geometry, is used for photometric correction of lunar visible/near-infrared data. A precise photometric correction parameter set is crucial for many applications including mineral identification and reflectance map mosaics. We present, for the first time, spectrally continuous photometric correction parameters for both sides of the Moon for wavelengths in the range 0.5-1.6 μm and solar phase angles between 5° and 85°, derived from Kaguya (SELENE) Spectral Profiler (SP) data. Since the measured radiance also depends on the surface albedo, we developed a statistical method for selecting areas with relatively uniform albedos from a nearly 7000-orbit SP data set. Using the selected data set, we obtained empirical photometric correction parameter sets for three albedo groups (high, medium, and low). We did this because the photometric function depends on the albedo, especially at phase angles below about 20° for which the shadow hiding opposition effect is appreciable. We determined the parameters in 160 bands and discovered a small variation in the opposition effect due to the albedo variation of mafic mineral absorption. The consistency of the photometric correction was checked by comparing observations made at different times of the same area on the lunar surface. Variations in the spectra obtained were lower than 2%, except for the large phase angle data in mare. Lastly, we developed a correction method for low solar elevation data, which is required for high latitude regions. By investigating low solar elevation data, we introduced an additional correction method. We used the new photometric correction to generate a 1° mesh global lunar reflectance map cube in a wavelength range of 0.5-1.6 μm. Surprisingly, these maps reveal that high latitude (?75°) regions in both the north and south have much lower spectral continuum slopes (color ratio r_1547.7nm/r_752.8nm ? 1.8) than the low and medium latitude regions, which implies lower degrees of space weathering.     

Yearly and seasonal variations of low albedo surfaces on Mars in the OMEGA/MEx dataset: Constraints on aerosols properties and dust deposits

The time variations of spectral properties of dark martian surface features are investigated using the OMEGA near-IR dataset. The analyzed period covers two Mars years, spanning from early 2004 to early 2008 (includes the 2007 global dust event). Radiative transfer modeling indicates that the apparent albedo variations of low to mid-latitude dark regions are consistent with those produced by the varying optical depth of atmospheric dust as measured simultaneously from the ground by the Mars Exploration Rovers. We observe only a few significant albedo changes that can be attributed to surface phenomena. They are small-scaled and located at the boundaries between bright and dark regions. We then investigate the variations of the mean particle size of aerosols using the evolution of the observed dark region spectra between 1 and 2.5 μm. Overall, we find that the observed changes in the spectral slope are consistent with a mean particle size of aerosols varying with time between 1 and 2 μm. Observations with different solar zenith angles make it possible to characterize the aerosol layer at different altitudes, revealing a decrease of the particle size of aerosols as altitude increases.     

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.     

Voyager photometry of Europa

Ninety voyager images ranging in phase angle from 3 to 143° and covering the spectral range from 0.34 to 0.58 μm were analyzed to derive the photometric properties of Europa. At small phase angles the disk-integrated phase curve is remarkable in that it shows little or no evidence of an opposition effect (in agreement with earlier Earth-based observations by Millis and Thompson, Icarus26, 408, 1975). The phase integral determined in the Voyager clear filter (centered near 0.47 μm) is 1.09 ± 0.11, in good agreement with previous estimates based on radiometry. The bolometric Bond albedo is 0.62 ± 0.14. The scattering properties of Europa in general, and of the two major terrain types (bright plains and darker mottled terrain) in particular, cannot be represented by a lunar-like photometric law. However, an equation which is a linear superposition of a lunar-like scattering law and a Lambert component provides an adequate simple representation of the scattering properties. The plains are photometrically more homogeneous than the darker mottled terrain. In the Voyager clear filter, the average normal reflectance is 0.71 for the plains on both the leading and trailing hemispheres; for the darker mottled terrain the values are 0.60 on the leading hemisphere, and 0.48 on the trailing one.     

Starspots: The zebra effect

Recent observations of brightness variations on the Sun during the solar cycle have motivated us to re-examine the widely held view that cool, dark starspots, covering a significant fraction of the star, are the centers of magnetic activity on BY Dra stars. We propose that the magnetic regions are better described by a bright facular network, and that the dark areas which give rise to photometric rotational modulation are actually regions where the underlying quiet photosphere is seen. This interpretation is consistent with recent observations of late-type stars that show that bright areas covering much of the star have magnetic fields with strengths of several thousand gauss. It resolves several problems with the current model, including the size, location, and stability of the starspots required to match photometric and Doppler-imaging observations. It also has interesting observational implications for the correlation of photometric rotational modulation and long term brightness variations with other surface activity, and for the positions of magnetically active stars in the H-R diagram.Hubble Fellow.     

Optical properties of dust from Jupiter family comets

To try to define specific physical properties of the dust of Jupiter-family comets (JFCs), we compare the light scattered by them. Amongst the more than 1000 JFCs, less than 200 are numbered, 40 of them being rather bright. In the present work we use data from the latter. In situ observations of three nuclei show low albedo surfaces. The albedo of the dust particles in the coma is low, with generally a red colour. The A(α)fρ product is a measure of cometary activity and secular changes. Images of different regions (jets and fans) give indications on the nucleus rotation and position of the emitting areas, as compared to the position of the rotation axis. Differences in physical properties between the particles in different regions are pointed out by differences in the linear polarization of the scattered light and by spectral variations in brightness and polarization. Jupiter family comets are considered as dust-poor comets. Tails and trails’ studies give an estimation of the size distribution of the particles. However the dust production rates depend on the largest particles (up to centimetre size), which are mainly observed in the trails where large dark compact particles are found. These dark particles are also responsible for the high polarization in the inner most coma of some comets. The meaning, in terms of physical properties, of the linear polarization is discussed through different examples such as 2P/Encke, 9P/Tempel 1 or the fragments of 73P/Schwassmann-Wachmann 3. Cometary outbursts and splitting events show that the properties of the dust ejected from the interior of the nucleus are similar to the ones of more active comets (new or with larger semi-major axis).     

Titan: Preliminary results on surface properties and photometry from VIMS observations of the early flybys

Cassini observations of the surface of Titan offer unprecedented views of its surface through atmospheric windows in the 1-5 μm region. Images obtained in windows for which the haze opacity is low can be used to derive quantitative photometric parameters such as albedo and albedo distribution, and physical properties such as roughness and particle characteristics. Images from the early Titan flybys, particularly T0, Ta, and T5 have been analyzed to create albedo maps in the 2.01 and 2.73 μm windows. We find the average normal reflectance at these two wavelengths to be 0.15±0.02 and 0.035±0.003, respectively. Titan's surface is bifurcated into two albedo regimes, particularly at 2.01 μm. Analysis of these two regimes to understand the physical character of the surface was accomplished with a macroscopic roughness model. We find that the two types of surface have substantially different roughness, with the low-albedo surface exhibiting mean slope angles of ∼18°, and the high-albedo terrain having a much more substantial roughness with a mean slope angle of ∼34°. A single-scattering phase function approximated by a one-term Henyey-Greenstein equation was also fit to each unit. Titan's surface is back-scattering (g∼0.3-0.4), and does not exhibit substantially different backscattering behavior between the two terrains. Our results suggest that two distinct geophysical domains exist on Titan: a bright region cut by deep drainage channels and a relatively smooth surface. The two terrains are covered by a film or a coating of particles perhaps precipitated from the satellite's haze layer and transported by eolian processes. Our results are preliminary: more accurate values for the surface albedo and physical parameters will be derived as more data is gathered by the Cassini spacecraft and as a more complete radiative transfer model is developed from both Cassini orbiter and Huygens Lander measurements.