Search

The Rosetta spacecraft successfully studied the 67P/Churyumov-Gerasimenko comet from 2014 to 2016, when the mission ended with an attempt to land the satellite on the comet nucleus. Despite the presence of the scientific camera OSIRIS on board the spacecraft, several researchers used images of the navigational camera (NavCam) for supporting their analysis of the target. The radiometric calibration of the NavCam allowed the construction of an extended catalogue of NavCam images for which the average radiance and reflectance of the comet nucleus was calculated. Here, that data set is used to fit Hapke’s model for the reflectance to deliver a set of parameters that describe the photometric properties of the comet nucleus. The results allowed for the calculation of the geometric and bond albedo of the target, as well as information on the single scattering albedo, average roughness angle, asymmetric factor and shadow-hiding opposition effect as a function of the phase angle.

This work is the second in a series of manuscripts devoted to the investigation of the spectrophotometric properties of Saturn’s rings from Cassini-VIMS (Visible and Infrared Mapping Spectrometer) observations. The dataset used for this analysis is represented by ten radial spectrograms of the rings which have been derived in Filacchione et al. (2014) by radial mosaics produced by VIMS. Spectrograms report the measured radiance factor I/F of the main rings of Saturn as a function of both radial distance (from 73500 to 141375 km) and wavelength (0.35–5.1 µm) for different observation geometries (phase angle ranging in the 2°–132° interval). We take advantage of a Monte Carlo ray-tracing routine (Ciarniello et al., 2014) to characterize the photometric behavior of the rings at each wavelength and derive the spectral Bond albedo of ring particles. This quantity is used to infer the composition of the regolith covering ring particles by applying Hapke’s theory. Four different regions, characterized by different optical depths, and respectively located in the C ring, inner B ring, mid B ring and A ring, have been investigated. Results from photometric modeling indicate that, in the VIS-NIR spectral range, B ring particles are intrinsically brighter than A and C ring particles, with the latter having the lowest albedo, while the single particle phase function of the ring’s particles is compatible with an Europa-like or Callisto-like formulation, depending on the investigated region. Spectral modeling of the inferred Bond albedo indicates that ring spectrum can be reproduced by water ice grains with inclusion of organic materials (tholin) as a UV absorber intimately mixed with variable amounts of other compounds in pure form (carbon, silicates) or embedded in water ice grains (nanophase hydrated iron oxides, carbon, silicates, crystalline hematite, metallic iron, troilite). The abundance of tholin decreases with radial distance from C ring (0.2–0.6%) to A ring (0.06%) for the selected regions. Its distribution is compatible with an intrinsic origin and is possibly related to the different plasma environment of the different ring regions. The identification of the other absorber(s) and its absolute volumetric abundance is uncertain, depending on the adopted grain size and mixing modality (intraparticle or intimate). However, assuming a common composition of the other absorber in the ring plane, we find that its abundance anti-correlates with the optical depth of the investigated regions, being maximum in the thinnest C ring and minimum in the thickest mid B ring. In the case of the C ring, an additional population of low-albedo grains is required to match the positive spectral slope of the continuum in the 0.55–2.2 µm interval, represented by an intraparticle mixture of water ice and a spectrum similar to troilite or metallic iron. The distribution of the darkening compounds is interpreted as the result of a contamination by exogenous material, which is more effective in the less dense regions of the rings because of their lower surface mass density of pure water ice.

The errors in measurements of the energy of the reflected solar radiation and the thermal radiation emitted from Earth entering space in all directions from the top of the atmosphere are analyzed. The potentials of measurements onboard the Lunar Observatory (LO) and spacecraft (SC) on geostationary and solar-synchronous orbits and at the L1 Lagrange point of the Sun–Earth system (SEL1) are compared. To take into account the radiation from the edge of the Earth’s disk, which cannot be observed from a geostationary orbiter at all phase angles, theoretical models should be constructed. The intensity distribution of radiation reflected to space through all other directions in dependence on the surface type, the radiation incidence angle, and the observation angle is also modeled. The actual inaccuracy of these methods is approximately 1%. At the SEL1 point, a SC may simultaneously observe the whole surface of the planet only at a phase angle α = 0° and only in the periods when the Moon does not eclipse the Earth’s disk. To take into account the radiation reflected to space at all other phase angles, it is necessary to construct theoretical models. This fact, as well as the lack of observations during lunar eclipses of the Earth, is connected with large errors in the required quantities. The system of two optical robotic telescopes mounted on the LO will sequentially survey the Earth’s surface almost in all ranges of phase angles. The annual averages of the Bond albedo and the emitted thermal radiation of the Earth as a planet determined at the LO and the corresponding deviations of the annual average global energy budget of the planet from the equilibrium state will be almost an order of magnitude more accurate than those determined from the data of any orbiter.

M. J. Hooton et al., [Context] The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin–orbit angle Ψ – a notoriously difficult parameter to measure – from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of Ψ for some systems., [Aims] We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b, which is characteristically similar to well-studied contemporaries such as KELT-9 b and WASP-33 b., [Methods] We obtained optical CHaracterising ExOPlanet Satellite (CHEOPS) transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 μm full-phase curve to model the asymmetric transits, occultations, and phase-dependent flux modulation. For the latter, we employed a novel physics-driven approach to jointly fit the phase modulation by generating a single 2D temperature map and integrating it over the two bandpasses as a function of phase to account for the differing planet–star flux contrasts. The reflected light component was modelled using the general ab initio solution for a semi-infinite atmosphere., [Results] When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of Ψ = 72.1−2.4+2.5 deg. From the occultations and phase variations, we derived dayside and nightside temperatures of 3062−68+66 K and 1720 ± 330 K, respectively.Our retrieval suggests that the dayside emission spectrum closely follows that of a blackbody. As the CHEOPS occultation is too deep to be attributed to blackbody flux alone, we could separately derive geometric albedo Ag = 0.171−0.068+0.066 and spherical albedo As = 0.266−0.100+0.097 from the CHEOPS data, and Bond albedoAB = 0.057−0.101+0.083 from the Spitzer phase curve.Although small, the Ag and As indicate that MASCARA-1 b is more reflective than most other ultra-hot Jupiters, where H− absorption is expected to dominate., [Conclusions] Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of Ψ. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally resolved James Webb Space Telescope phase curve data., Y.A. acknowledge the support of the Swiss National Fund under grant 200020_172746. S.H. gratefully acknowledges CNES funding through the grant 837319. D.K. acknowledges partial financial support from the Center for Space and Habitability (CSH), the PlanetS National Center of Competence in Research (NCCR), and the Swiss National Science Foundation and the Swiss-based MERAC Foundation. A.C.C. and T.G.W. acknowledge support from STFC consolidated grant number ST/M001296/1. P.M. acknowledges support from STFC research grant number ST/M001040/1. This work was also partially supported by a grant from the Simons Foundation (PI Queloz, grant number 327127). B.-O.D. acknowledges support from the Swiss National Science Foundation (PP00P2-190080). S.S. has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 833925, project STAREX). A.Br. was supported by the SNSA. This work was supported by FCT – Fundação para a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020 – Programa Operacional Competitividade e Internacionalizacão by these grants: UID/FIS/04434/2019, UIDB/04434/2020, UIDP/04434/2020, PTDC/FIS-AST/32113/2017 and POCI-01-0145-FEDER- 032113, PTDC/FIS-AST/28953/2017 and POCI-01-0145-FEDER-028953, PTDC/FIS-AST/28987/2017 and POCI-01-0145-FEDER-028987, O.D.S.D. is supported in the form of work contract (DL 57/2016/CP1364/CT0004) funded by national funds through FCT. We acknowledge support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grants ESP2016-80435-C2-1-R, ESP2016-80435-C2-2-R, PGC2018-098153-B-C33, PGC2018-098153-B-C31, ESP2017-87676-C5-1-R, MDM-2017-0737 Unidad de Excelencia Maria de Maeztu-Centro de Astrobiologí a (INTA-CSIC), as well as the support of the Generalitat de Catalunya/CERCA programme. The MOC activities have been supported by the ESA contract No. 4000124370. S.C.C.B. acknowledges support from FCT through FCT contracts nr. IF/01312/2014/CP1215/CT0004. X.B., S.C., D.G., M.F., and J.L. acknowledge their roles as ESA-appointed CHEOPS science team members. This project was supported by the CNES. The Belgian participation to CHEOPS has been supported by the Belgian Federal Science Policy Office (BELSPO) in the framework of the PRODEX Program, and by the University of Liège through an ARC grant for Concerted Research Actions financed by the Wallonia-Brussels Federation. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (project FOUR ACES. grant agreement no. 724427). D.E. acknowledges financial support from the Swiss National Science Foundation for project 200021_200726. C.M.P. and M.F. gratefully acknowledge the support of the Swedish National Space Agency (DNR 65/19, 174/18). D.G. gratefully acknowledges financial support from the CRT foundation under Grant No. 2018.2323 ‘Gaseousor rocky? Unveiling the nature of small worlds’. M.G. is an F.R.S.-FNRS Senior Research Associate. KGI is the ESA CHEOPS Project Scientist and is responsible for the ESA CHEOPS Guest Observers Programme. She does not participate in, or contribute to, the definition of the Guaranteed Time Programme of the CHEOPS mission through which observations described in this paper have been taken, nor to any aspect of target selection for the programme. This work was granted access to the HPC resources of MesoPSL financed by the Region Ile de France and the project Equip@Meso (reference ANR-10-EQPX-29-01) of the programme Investissements d’Avenir supervised by the Agence Nationale pour la Recherche. Acknowledges support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grant PGC2018-098153-B-C33, as well as the support of the Generalitat de Catalunya/CERCA programme. S.G.S. acknowledges support from FCT through FCT contract nr. CEECIND/00826/2018 and POPH/FSE (EC). This project has been supported by the Hungarian National Research, Development and Innovation Office (NKFIH) grants GINOP-2.3.2-15-2016-00003, K-119517, K-125015, and the City of Szombathely under Agreement No. 67.177-21/2016., , , ,

Context: We report the discovery of TOI-519 b (TIC 218795833), a transiting substellar object (R = 1.07 RJup) orbiting a faint M dwarf (V = 17.35) on a 1.26 d orbit. Brown dwarfs and massive planets orbiting M dwarfs on short-period orbits are rare, but more have already been discovered than expected from planet formation models. TOI-519 is a valuable addition into this group of unlikely systems, and adds towards our understanding of the boundaries of planet formation. Aims: We set out to determine the nature of the Transiting Exoplanet Survey Satellite (TESS ) object of interest TOI-519 b. Methods: Our analysis uses a SPOC-pipeline TESS light curve from Sector 7, multicolour transit photometry observed with MuSCAT2 and MuSCAT, and transit photometry observed with the LCOGT telescopes. We estimate the radius of the transiting object using multicolour transit modelling, and set upper limits for its mass, effective temperature, and Bond albedo using a phase curve model that includes Doppler boosting, ellipsoidal variations, thermal emission, and reflected light components. Results: TOI-519 b is a substellar object with a radius posterior median of 1.07 RJup and 5th and 95th percentiles of 0.66 and 1.20 RJup, respectively, where most of the uncertainty comes from the uncertainty in the stellar radius. The phase curve analysis sets an upper effective temperature limit of 1800 K, an upper Bond albedo limit of 0.49, and a companion mass upper limit of 14 MJup. The companion radius estimate combined with the Teff and mass limits suggests that the companion is more likely a planet than a brown dwarf, but a brown-dwarf scenario is more likely a priori given the lack of known massive planets in 1 day orbits around M dwarfs with Teff < 3800 K, and the existence of some (but few) brown dwarfs., Accepted to A&A

We are compiling a service, which will produce almost real-time estimates of the Earth Bond albedo in visual and near-infrared wavelengths using the multi-wavelength images from the NOAA Deep Space Climate Observatory’s (DSCOVR) Earth Polychromatic Imaging Camera (EPIC). The service will be based on the known land and sea surface types of the Earth and the angular distribution models (ADMs) from the Clouds and the Earth’s Radiant Energy System (CERES) project (some details in Ihalainen, M.Sc. thesis, 2019, University of Helsinki, Finland). Instead of applying the CERES ADM albedo estimates as such, we will only use the ADM shapes for different land and sea surface types, but scale them using the corresponding observations from the EPIC camera. As the EPIC camera will produce imaging of the Earth disk from the L1 point, in about every two hours and in 10 narrow band filters between 317–780 nm, we can follow the temporal changes in the Bond albedo both in short times scales and over longer periods. We will employ a custom-made classifier for cloud coverage over the Earth disk using the 10 bands in the EPIC images to improve the Bond albedo estimation with temporally suitable cloud coverage ADMs over the land or sea surface pixels. The service will be run automatically in a dedicated server, and will produce the Bond albedo estimates via a public web interface. Figure 1. An example image from EPIC camera, taken at Jan. 1st, 2019. The three first images, counting from the top left corner, are gray-scale images at wavelength bands 317 nm, 388 nm, and 780 nm. The fourth image is false-color RGB composite from those channels. Acknowledgements: Research supported by the Academy of Finland (project 298137).

On 11th April 2015 Cassini's Composite Infrared Spectrometer (CIRS) made a series of observations of Tethys’ daytime anti-Saturn hemisphere over a nine-hour time period. During this time the sub-spacecraft position was remarkably stable (0.3° S to 3.9° S; 153.2° W to 221.8° W), and so these observations provide unprecedented coverage of diurnal temperature variations on Tethys’ anti-Saturn hemisphere. In 2012 a thermal anomaly was discovered at low latitudes on Tethys’ leading hemisphere; it appears cooler during the day and warmer at night than its surroundings (Howett et al., 2012) and is spatially correlated with a decrease in the IR3/UV3 visible color ratio (Schenk et al., 2011). The cause of this anomaly is believed to be surface alteration by high-energy electrons, which preferentially bombard low-latitudes of Tethys’ leading hemisphere (Schenk et al., 2011; Howett et al., 2012; Paranicas et al. 2014; Schaible et al., 2017). The thermal anomaly was quickly dubbed “Pac-Man” due to its resemblance to the 1980s video game icon. We use these daytime 2015 CIRS data, along with two sets of nighttime CIRS observations of Tethys (from 27 June 2007 and 17 August 2015) to make maps of bolometric Bond albedo and thermal inertia variations across the anti-Saturn hemisphere of Tethys (including the edge of its Pac-Man region). These maps confirm the presence of the Pac-Man thermal anomaly and show that while Tethys’ bolometric Bond albedo varies negligibly outside and inside the anomaly (0.69 ± 0.02 inside, compared to 0.71 ± 0.04 outside) the thermal inertia varies dramatically (29 ± 10 J m−2 K−1 s−1/2 inside, compared to 9 ± 4 J m−2 K−1 s−1/2 outside). These thermal inertias are in keeping with previously published values: 25 ± 3 J m−2 K−1 s−1/2 inside, and 5 ± 1 J m−2 K−1 s−1/2 outside the anomaly (Howett et al., 2012). A detailed analysis shows that on smaller spatial-scales the bolometric Bond albedo does vary: increasing to a peak value at 180° W. For longitudes between ∼100° W and ∼160° W the thermal inertia increases from northern to southern latitudes, while the reverse is true for bolometric Bond albedo. The thermal inertia on Tethys generally increases towards the center of its leading hemisphere but also displays other notable small-scale variations. These thermal inertia and bolometric Bond albedo variations are perhaps due to differences in competing surface modification by E ring grains and high-energy electrons which both bombard Tethys’ leading hemisphere (but in different ways). A comparison between the observed temperatures and our best thermal model fits shows notable discrepancies in the morning warming curve, which may provide evidence of regional variations in surface roughness effects, perhaps again due to variations in surface alteration mechanisms.

International audience; Asteroids (24) Themis and (65) Cybele have an absorption feature at 3.1 μm reported to be directly linked to surface water ice. We searched for water vapor escaping from these asteroids with the Herschel Space Observatory Heterodyne Instrument for the Far Infrared. While no H 2 O line emission was detected, we obtain sensitive 3σ water production rate upper limits of Q(H 2 O)

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~\mu$m and $0.30 \pm 0.02 \%$ at $4.5~\mu$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~\mu$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~\mu$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters., Comment: 14 pages, 8 figures. Accepted for publication in AJ

We have made high precision polarimetric observations of the polluted white dwarf G29-38 with the HIgh Precision Polarimetric Instrument 2. The observations were made at two different observatories -- using the 8.1-m Gemini North Telescope and the 3.9-m Anglo AustralianTelescope -- and are consistent with each other. After allowing for a small amount of interstellar polarization, the intrinsic linear polarization of the system is found to be 275.3 +/- 31.9 parts-per-million at a position angle of 90.8 +/- 3.8 degrees in the SDSS g' band. We compare the observed polarization with the predictions of circumstellar disc models. The measured polarization is small in the context of the models we develop which only allows us to place limits on disc inclination and Bond albedo for optically thin disc geometries. In this case either the inclination is near face-on or the albedo is small -- likely in the range 0.05 to 0.15 -- which is in line with other debris disc measurements. A preliminary search for the effects of G29-38's pulsations in the polarization signal produced inconsistent results. This may be caused by beating effects, indicate a clumpy dust distribution, or be a consequence of measurement systematics., 15 pages, 6 figures, 4 tables. Accepted to MNRAS

On January 1st 2019, the New Horizons spacecraft flew by the classical Kuiper belt object (486958) Arrokoth (provisionally designated 2014 MU69), possibly the most primitive object ever explored by a spacecraft. The I/F of Arrokoth is analyzed and fit with a photometric function that is a linear combination of the Lommel-Seeliger (lunar) and Lambert photometric functions. Arrokoth has a geometric albedo of p_V = 0.21_(-0.04)^(+0.05) at a wavelength of 550 nm and ~0.24 at 610 nm. Arrokoth's geometric albedo is greater than the median but consistent with a distribution of cold classical Kuiper belt objects whose geometric albedos were determined by fitting a thermal model to radiometric observations. Thus, Arrokoth's geometric albedo adds to the orbital and spectral evidence that it is a cold classical Kuiper belt object. Maps of the normal reflectance and hemispherical albedo of Arrokoth are presented. The normal reflectance of Arrokoth's surface varies with location, ranging from ~0.10-0.40 at 610 nm with an approximately Gaussian distribution. Both Arrokoth's extrema dark and extrema bright surfaces are correlated to topographic depressions. Arrokoth has a bilobate shape and the two lobes have similar normal reflectance distributions: both are approximately Gaussian, peak at ~0.25 at 610 nm, and range from ~0.10-0.40, which is consistent with co-formation and co-evolution of the two lobes. The hemispherical albedo of Arrokoth varies substantially with both incidence angle and location, the average hemispherical albedo at 610 nm is 0.063 +/- 0.015. The Bond albedo of Arrokoth at 610 nm is 0.062 +/- 0.015., In press in Icarus

The exploration of the Pluto-Charon system by the New Horizons spacecraft represents the first opportunity to understand the distribution of albedo and other photometric properties of the surfaces of objects in the Solar System's "Third Zone" of distant ice-rich bodies. Images of the entire illuminated surface of Pluto and Charon obtained by the Long Range Reconnaissance Imager (LORRI) camera provide a global map of Pluto that reveals surface albedo variegations larger than any other Solar System world except for Saturn's moon Iapetus. Normal reflectances on Pluto range from 0.08-1.0, and the low-albedo areas of Pluto are darker than any region of Charon. Charon exhibits a much blander surface with normal reflectances ranging from 0.20-0.73. Pluto's albedo features are well-correlated with geologic features, although some exogenous low-albedo dust may be responsible for features seen to the west of the area informally named Tombaugh Regio. The albedo patterns of both Pluto and Charon are latitudinally organized, with the exception of Tombaugh Regio, with darker regions concentrated at the Pluto's equator and Charon's northern pole The phase curve of Pluto is similar to that of Triton, the large moon of Neptune believed to be a captured Kuiper Belt Object (KBO), while Charon's is similar to that of the Moon. Preliminary Bond albedos are 0.25+/-0.03 for Charon and 0.72+/-0.07 for Pluto. Maps of an approximation to the Bond albedo for both Pluto and Charon are presented for the first time. Our work shows a connection between very high albedo (near unity) and planetary activity, a result that suggests the KBO Eris may be currently active., Comment: 32 pages, 3 tables, 8 figures Accepted for Publication in Icarus, 2016, v02

During accretion, the young rocky planets are so hot that they become endowed with a magma ocean. From that moment, the mantle convective thermal flux control the cooling of the planet and an atmosphere is created by outgassing. This atmosphere will then play a key role during this cooling phase. Studying this cooling phase in details is a necessary step to explain the great diversity of the observed telluric planets and especially to assess the presence of surface liquid water. We used here a radiative-convective 1D atmospheric model (H2O, CO2) to study the impact of the Bond albedo on the evolution of magma ocean planets. We derived from this model the thermal emission spectrum and the spectral reflectance of these planets, from which we calculated their Bond albedos. Compared to Marcq et al. (2017), the model now includes a new module to compute the Rayleigh scattering, and state of the art CO2 and H2O gaseous opacities data in the visible and infrared spectral ranges. We show that the Bond albedo of these planets depends on their surface temperature and results from a competition between Rayleigh scattering from the gases and Mie scattering from the clouds. The colder the surface temperature is, the thicker the clouds are, and the higher the Bond albedo is. We also evidence that the relative abundances of CO2 and H2O in the atmosphere strongly impact the Bond albedo. The Bond albedo is higher for atmospheres dominated by the CO2, better Rayleigh scatterer than H2O. Finally, we provide the community with an empirical formula for the Bond albedo that could be useful for future studies of magma ocean planets., 16 pages, 10 figures, 5 tables. Accepted for publication in the Journal Icarus the 27 August 2018

133P/Elst-Pizarro is the firstly recognized main-belt comet, but we still know little about its nucleus. Firstly we use mid-infrared data of Spitzer-MIPS, Spitzer-IRS and WISE to estimate its effective diameter $D_{\rm eff}=3.9^{+0.4}_{-0.3}$ km, geometric albedo $p_{\rm v}=0.074\pm0.013$ and mean Bond albedo $A_{\rm eff,B}=0.024\pm0.004$. The albedo is used to compute 133P's temperature distribution, which shows significant seasonal variation, especially polar regions, ranging from $\sim40$ to $\sim200$ K. Based on current activity observations, the maximum water gas production rate is estimated to be $\sim1.4\times10^{23}\rm~s^{-1}$, being far weaker than $\sim10^{26}\rm~s^{-1}$ of JFC 67P at similar helio-centric distance $\sim2.7$ AU, indicating a thick dust mantle on the surface to lower down the gas production rate. The diameter of the sublimation area may be $100$ Myr, indicating that 133P may be more likely to be a relatively old planetesimals or a member of an old family than a recently formed fragment of some young family., Comment: 11 pages, 9 figures, accepted for publication in The Astronomical Journal

We report a comprehensive analysis of the global spectrophotometric properties of Ceres using Dawn Framing Camera images collected from April to June 2015 during the RC3 and Survey mission phases. The single-scattering albedo of Ceres at 555 nm is 0.14$\pm$0.04, the geometric albedo is 0.096$\pm$0.006, and the Bond albedo is 0.037$\pm$0.002. The asymmetry factors calculated from the best-fit two-term Henyey-Greenstein (HG) single-particle phase function (SPPF) show a wavelength dependence, suggesting that the phase reddening of Ceres is dominated by single-particle scattering rather than multiple scattering or small-scale surface roughness. The Hapke roughness parameter of Ceres is derived to be 20$^\circ\pm$6$^\circ$ with no wavelength dependence. The phase function of Ceres shows appreciably strong scattering around 90$^\circ$ phase angle that cannot be fitted with a single-term HG SPPF, suggesting possible stronger forward scattering than other asteroids previously analyzed with spacecraft data. We speculate that such a scattering characteristic of Ceres might be related to its unique surface composition. We grouped the reflectance data into a 1$^\circ$ latitude-longitude grid and fitted each grid independently to study the spatial variations of photometric properties. The albedo and color maps are consistent with previous studies. The SPPF over the surface of Ceres shows stronger backscattering associated with lower albedo and vice versa, consistent with the general trend among asteroids. The Hapke roughness parameter does not vary much across the surface of Ceres, except for the ancient Vendimia Planitia region that has a slightly higher roughness. Based on the wavelength dependence of the SPPF of Ceres, we hypothesize that its regolith grains either contain a considerable fraction of $\lessapprox\mu$m-sized particles, or are strongly affected by internal scatterers of this size., Comment: 43 pages, 3 tables, 17 figures, accepted by Icarus

Most known terrestrial planets orbit small stars with radii less than 60% that of the Sun. Theoretical models predict that these planets are more vulnerable to atmospheric loss than their counterparts orbiting Sun-like stars. To determine whether a thick atmosphere has survived on a small planet, one approach is to search for signatures of atmospheric heat redistribution in its thermal phase curve. Previous phase curve observations of the super-Earth 55 Cancri e (1.9 Earth radii) showed that its peak brightness is offset from the substellar point $-$ possibly indicative of atmospheric circulation. Here we report a phase curve measurement for the smaller, cooler planet LHS 3844b, a 1.3 Earth radius world in an 11-hour orbit around a small, nearby star. The observed phase variation is symmetric and has a large amplitude, implying a dayside brightness temperature of $1040\pm40$ kelvin and a nightside temperature consistent with zero kelvin (at one standard deviation). Thick atmospheres with surface pressures above 10 bar are ruled out by the data (at three standard deviations), and less-massive atmospheres are unstable to erosion by stellar wind. The data are well fitted by a bare rock model with a low Bond albedo (lower than 0.2 at two standard deviations). These results support theoretical predictions that hot terrestrial planets orbiting small stars may not retain substantial atmospheres., Comment: Published in Nature on Aug 19, 2019. Co-authors Koll, Morley and Hu contributed equally to this manuscript

The potential habitability of known exoplanets is often categorized by a nominal equilibrium temperature assuming a Bond albedo of either 0.3, similar to Earth, or 0. As an indicator of habitability, this leaves much to be desired, because albedo on other planets can be very different, and because surface temperature exceeds equilibrium temperature due to the atmospheric greenhouse effect. We use an ensemble of 3-dimensional general circulation model simulations to show that for a range of habitable planets, much of the variability of Bond albedo, equilibrium temperature, and even surface temperature can be predicted with useful accuracy from incident stellar flux and stellar temperature, two known external parameters for every confirmed exoplanet., 45 pages including 8 figures, 5 tables; submitted to The Astrophysical Journal

NASA's OSIRIS-REx spacecraft arrived at its sampling target, asteroid (101955) Bennu, in December 2018 and started a series of global observation campaigns. Here we investigate the global photometric properties of Bennu as observed by the OSIRIS-REx Visible and InfraRed Spectrometer (OVIRS) over the time period December 9, 2018, to September 26, 2019. In this study we used observations obtained over wavelengths ranging from 0.4 to 3.7 μm, with a solar phase angle range of 5.3° to 132.6°. Our aim is to characterize the global average disk-resolved photometric properties of Bennu with multiple models. The best-fit model is a McEwen model with an exponential phase function and an exponential polynomial partition function. We use this model to correct the OVIRS spectra of Bennu to a standard reference viewing and illumination geometry at visible to infrared wavelengths for the purposes of global spectral mapping. We derive a bolometric Bond albedo map in which Bennu's surface values range from 0.021 to 0.027. We find a phase reddening effect, and our model is effective at removing this phase reddening. Our average model albedo shows a blueish spectrum with a > 10% absorption feature centered at 2.74 μm. Of all comparisons with previously visited asteroids and comets, only 28P/Neujmin, 2P/Encke, and (162173) Ryugu are darker than Bennu. We find that Bennu is a few percent brighter than Ryugu in the wavelengths respectively observed by the OSIRIS-REx and Hayabusa2 missions (from 0.48 to 0.86 μm). We also compare our spectroscopic photometry of Bennu with the OSIRIS-REx imaging photometry and with ground-based predictions.

The observations of the Cassini spacecraft gave new perspectives to our understanding of the rings of Saturn and their thermal studies due to their large variety of observation geometries. However, we still know little about the ring particles' physical and, particularly, thermal properties, which is the reason why a number of observations are not easy to explain. Nevertheless, the observed temperature profiles of the rings can be reproduced by considering the rings' general structural properties like local optical depth, τ, and, for the optically-thick rings, wake properties like the ratios H/W or wake thickness to wake width, and S/W or wake separation to wake width; as well as their main energy input sources (e. g., direct and reflected solar energy and Saturn's infrared energy). In this work, we assume that the main rings (C, B, Cassini Division and A rings) are monolayers and study the temperature variations with solar elevation angle, B′, of ten of their regions with a simple thermal balance equation that depends on the rings' shadowing behaviour due to the illumination geometry (B′) and the arrangement of particles across the rings (τ). The shadowing is a key feature, which we reproduce with a numerical function obtained via ray tracing from arrays of spherical particles (C(B′, τ)). Using the former model, we fit Cassini Composite Infrared Spectrometer (CIRS) temperature data with the Bond albedo, AV, and the so-called rotation parameter, f (which represents the rotation status of the ring particles) as only parameters. These parameters are two of the most important physical properties that define the rings' temperatures. With our method, we are able to reproduce the CIRS data temperature trends of the selected regions reasonably well. For evaluation purposes we compare the derived Bond albedo values with those obtained with the Cassini Imaging Science Subsystem (ISS) at wavelengths in the range 451 ≲ λ ≲ 650 nm. The comparison suggests a dominance of slow rotators. Additionally, we evaluate the differences obtained when using either a numerical or an analytical shadowing function (in particular, that by ( Altobelli et al., 2008 )) in the thermal model. The differences between these two approaches reflect some details of the structure of the ring regions we study.

Here, the author presents an analysis of the evolution of global temperature and Bond albedo anomalies in the period 2001-2013. Based on the radiative balance at the top of the atmosphere, an equation is presented to estimate the temperature anomaly following a particular albedo anomaly, and vice versa. The estimated values are compared with the observed values.

The Earth's albedo is a fundamental climate parameter for understanding the radiation budget of the atmosphere. It has been traditionally measured not only from space platforms but also from the ground for 16 years from Big Bear Solar Observatory by observing the Moon. The photometric ratio of the dark (earthshine) to the bright (moonshine) sides of the Moon is used to determine nightly anomalies in the terrestrial albedo, with the aim of quantifying sustained monthly, annual, and/or decadal changes. We find two modest decadal scale cycles in the albedo, but with no significant net change over the 16 years of accumulated data. Within the evolution of the two cycles, we find periods of sustained annual increases, followed by comparable sustained decreases in albedo. The evolution of the earthshine albedo is in remarkable agreement with that from the Clouds and the Earth's Radiant Energy System instruments, although each method measures different slices of the Earth's Bond albedo.

In 2011 a thermally anomalous region was discovered on Mimas, Saturn's innermost major icy satellite (Howett et al., 2011). The anomalous region is a lens-like shape located at low latitudes on Mimas' leading hemisphere. It manifests as a region with warmer nighttime temperatures, and cooler daytime ones than its surroundings. The thermally anomalous region is spatially correlated with a darkening in Mimas' IR/UV surface color (Schenk et al., 2011) and the region preferentially bombarded by high-energy electrons (Paranicas et al., 2012, Paranicas et al., 2014; Nordheim et al., 2017). We use data from Cassini's Composite Infrared Spectrometer (CIRS) to map Mimas' surface temperatures and its thermophysical properties. This provides a dramatic improvement on the work in Howett et al. (2011), where the values were determined at only two regions on Mimas (one inside, and another outside of the anomalous region). We use all spatially-resolved scans made by CIRS' focal plane 3 (FP3, 600 to 1100 cm−1) of Mimas' surface, which are largely daytime observations but do include one nighttime one. The resulting temperature maps confirm the presence and location of Mimas' previously discovered thermally anomalous region. No other thermally anomalous regions were discovered, although we note that the surface coverage is incomplete on Mimas' leading and anti-Saturn hemisphere. The thermal inertia map confirms that the anomalous region has a notably higher thermal inertia than its surroundings: 98 ± 42 J m−2 K−1 s-1/2 inside of the anomaly, compared to 34 ± 32 J m−2 K−1 s-1/2 outside. The albedo inside and outside of the anomalous region agrees within their uncertainty: 0.45 ± 0.08 inside compared to 0.41 ± 0.07 outside the anomaly. Interestingly the albedo appears brighter inside the anomaly region, which may not be surprising given this region does appear brighter at some UV wavelengths (0.338 μm, see Schenk et al., 2011). However, this result should be treated with caution because, as previously stated, statistically the albedo of these two regions is the same when their uncertainties are considered. These thermal inertia and albedo values determined here are consistent with those found by Howett et al. (2011), who determined the thermal inertia inside the anomaly to be 66 ± 23 J m−2 K−1 s-1/2 and

In this paper we quantify and attribute by inspection the constituent elements of the power intensity radiant flux transmission for the atmosphere of the Earth, as recorded in the following two published sources; Oklahoma Climatological Survey and Kiehl and Trenberth. The purpose of our analysis is to establish the common elements of the approach used in the formulation of these works, and to conduct an assessment of the two approaches by establishing a common format for their comparison. By applying the standard analysis of a geometric infinite series feed-back loop to an equipartition (half up and half down) diabatic distribution used for the atmospheric radiant flux to all elements of the climate model; our analysis establishes the relative roles of radiant and mass-motion carried energy fluxes that are implicitly used by the authors in their respective analyses. Having established the key controls on energy flux within each model, we then conduct for the canonical model a series of “what-if” scenarios to establish the limits of temperature rise that can be achieved for specific variations in the controls used to calculate the global average temperature. Our analysis establishes that, for the current insolation and Bond albedo, the maximum temperature that can be achieved for a thermally radiant opaque atmosphere is a rise to 29°C. This global average temperature is achieved by a total blocking of the surface-to-space atmospheric window. In order to raise the global average atmospheric temperature to the expected value of 36°C for a putative Cretaceous hothouse world, it is therefore necessary to reduce the planetary Bond albedo. The lack of continental icecaps, and the presence of flooded continental shelves with epeiric seas in a global eustatic high stand sea level, is invoked as an explanation to support the modelling concept of a reduced global Bond albedo during the Cretaceous period. The geological evidence for this supposition is mentioned with reference to published sources.

The radiant energy budget and internal heat are fundamental properties of giant planets, but precise determination of these properties remains a challenge. Here, we report measurements of Jupiter’s radiant energy budget and internal heat based on Cassini multi-instrument observations. Our findings reveal that Jupiter’s Bond albedo and internal heat, 0.503 ± 0.012 and 7.485 ± 0.160 W m−2 respectively, are significantly larger than 0.343 ± 0.032 and 5.444 ± 0.425 Wm−2, the previous best estimates. The new results help constrain and improve the current evolutionary theories and models for Jupiter. Furthermore, the significant wavelength dependency of Jupiter’s albedo implies that the radiant energy budgets and internal heat of the other giant planets in our solar system should be re-examined. Finally, the data sets of Jupiter’s characteristics of reflective solar spectral irradiance provide an observational basis for the models of giant exoplanets.

This study develops an equilibrium radiative model based on a quasi-adiabatic atmosphere that quantifies the average surface flux of a planetary body as a function of absorbed solar radiation P and the bulk emissivity of the atmosphere with respect to surface radiation &nbsp, &epsilon, The surface flux is then given by &nbsp, P/(1&minus, ), and I define the scaling term 1/(1&minus, ) as the greenhouse factor. The model is applied to all of the rocky planets in the solar system to determine their greenhouse factors, and accounts for the diversity of planetary surface fluxes. The model is modified to allow for a top of atmosphere non-equilibrium state, which when compared with a recent observation-based model of the Earth energy budget, predicts the Earth&rsquo, s radiative fluxes to within the uncertainty ranges of that model. The model developed in this study is able to quantify the changes in Earth&rsquo, s surface flux caused by changes in bond albedo and atmospheric bulk emissivity by using the surface temperature, ocean heat content, incoming solar radiation and outgoing longwave radiation records. The model indicates an increase in absorbed solar radiation over the time period from 1979&ndash, 2015 of the order of 3 W/m2, which was caused by a decrease in planetary bond albedo. The time-series albedo generated by the model is in agreement with Clouds and Earth&rsquo, s Radiant Energy System (CERES) derived albedo over the period from 2000&ndash, 2015. The model also indicates a slight decrease in atmospheric bulk emissivity over the same period. Since atmospheric bulk emissivity is a function of the sum of all of the greenhouse gas species, a simultaneous decrease in atmospheric water vapor may offset the effect of the well-documented increase in the non-condensing greenhouse gases over the period, and result in an overall net decrease in bulk emissivity. Atmospheric water vapor datasets partially support the conclusion, with the International Satellite Cloud Climatology Project (ISCCP) data supporting a decrease. The NASA Water Vapor Project (NVAP-M) data supports a decrease in atmospheric water content over the period 1998&ndash, 2008, but not over the longer period of 1988&ndash, 2008. The model indicates that the decrease in planetary albedo was the driver for the increased surface flux over the stated period.

We analyzed the photometric properties of Asteroid (21) Lutetia based on images captured by Rosetta during its flyby. We utilized the images recorded in the F17 filter (λ = 631.6 nm) of the Wide Angle Camera (WAC) and in the F82 & F22 filters (λ = 649.2 nm) of the Narrow Angle Camera (NAC) of the OSIRIS imaging system onboard the spacecraft. We present the results of Hapke and Minnaert modeling using disk-integrated and disk-resolved data derived from the surface of the asteroid. At 631.6 nm and 649.2 nm, the geometric albedo of Lutetia is 0.194 ± 0.002. The Bond albedo is 0.076 ± 0.002 at 649.2 nm, and 0.079 ± 0.002 at 631.6 nm. The roughness parameter is 28 ° ± 1 ° , the opposition surge parameters B0 and h are 1.79 ± 0.08 and 0.041 ± 0.003, respectively, and the asymmetry factor of the phase function is −0.28 ± 0.01. The single-scattering albedo is 0.226 ± 0.002 at 631.6 and 649.2 nm. The modeled Hapke parameters of Asteroid Lutetia are close to those of typical S-type asteroids. The Minnaert k parameter of Lutetia at opposition ( 0.526 ± 0.002 ) is comparable with other asteroids and comets. Albedo ratio images indicate no significant variation across the surface of Lutetia, apart from the so called NPCC region on Lutetia where a pronounced variation is seen at large phase angle. The small width of the albedo distribution of the surface (∼7% at half maximum) and the similarity between phase ratio maps derived from the measurements and from the modeling suggests that the light scattering property over the whole visible and illuminated surface of the asteroid is widely uniform. The comparison between the reflectance measurement of Lutetia and the available laboratory samples suggests that the regolith on Lutetia is concrete with possible grain size distribution of150 μm or larger.

We used ground-based photometric phase curve data of the OSIRIS-REx target Asteroid (101955) Bennu and low phase angle data from Asteroid (253) Mathilde as a proxy to fit Bennu data with Minnaert, Lommel-Seeliger, (RObotic Lunar Orbiter) ROLO, Hapke, and McEwen photometric models, which capture the global light scattering properties of the surface and subsequently allow us to calculate the geometric albedo, phase integral, spherical Bond albedo, and the average surface normal albedo for Bennu. We find that Bennu has low reflectance and geometric albedo values, such that multiple scattering is expected to be insignificant. Our photometric models relate the reflectance from Bennu’s surface to viewing geometry as functions of the incidence, emission, and phase angles. Radiance Factor functions (RADFs) are used to model the disk-resolved brightness of Bennu. The Minnaert, Lommel-Seeliger, ROLO, and Hapke photometric models work equally well in fitting the best ground-based photometric phase curve data of Bennu. The McEwen model works reasonably well at phase angles from 20° to 70°. Our calculated geometric albedo values of 0.047 - 0.014 + 0.012 , 0.047 - 0.014 + 0.005 , and 0.048 - 0.022 + 0.012 for the Minnaert, the Lommel-Seeliger, and the ROLO models respectively are consistent with the geometric albedo of 0.045 ± 0.015 computed by Emery et al. (Emery, J.P. et al. [2014]. Icarus 234, 17–35) and Hergenrother et al. (Hergenrother, C.W. et al. [2014]. http://arxiv.org/abs/1409.4704 >). Also, our spherical Bond albedo values of 0.016 - 0.004 + 0.005 , 0.015 - 0.001 + 0.003 , and 0.015 - 0.005 + 0.007 for the Minnaert model, Lommel-Seeliger, and ROLO models respectively are consistent with the value of 0.017 ± 0.002 presented by Emery et al. (Emery, J.P. et al. [2014]. Icarus 234, 17–35). On the other hand, the semi-physical models such as the Hapke model, where several assumptions and approximations were necessary, and the McEwen model are not supported by the global disk-integrated data, indicating that disk-resolved measurements will be necessary to constrain these models, as expected.

We consider fifty transiting short-period giant planets for which eclipse depths have been measured at multiple infrared wavelengths. The aggregate dayside emission spectrum of these planets exhibits no molecular features, nor is brightness temperature greater in the near-infrared. We combine brightness temperatures at various infrared wavelengths to estimate the dayside effective temperature of each planet. We find that dayside temperatures are proportional to irradiation temperatures, indicating modest Bond albedo and no internal energy sources, plus weak evidence that dayside temperatures of the hottest planets are disproportionately high. We place joint constraints on Bond albedo, $A_{B}$, and day-to-night transport efficiency, $\varepsilon$, for six planets by combining thermal eclipse and phase variation measurements (HD 149026b, HD 189733b, HD 209458b, WASP-12b, WASP-18b, and WASP-43b). We confirm that planets with high irradiation temperatures have low heat transport efficiency, and that WASP-43b has inexplicably poor transport; these results are statistically significant even if the precision of single-eclipse measurements has been overstated by a factor of three. Lastly, we attempt to break the $A_{B}$-$\varepsilon$ degeneracy for nine planets with both thermal and optical eclipse observations, but no thermal phase measurements. We find a systematic offset between Bond albedos inferred from thermal phase variations ($A_{B} \approx 0.35$) and geometric albedos extracted from visible light measurements ($A_{g} \approx 0.1$). These observations can be reconciled if most hot Jupiters have clouds that reflect 30-50 per cent of incident near-infrared radiation, as well as optical absorbers in the cloud particles or above the cloud deck., Comment: 12 pages, 7 figures, 4 tables; submitted to MNRAS, post-peer review

Thermal phase curves of short-period planets on circular orbits provide joint constraints on the fraction of incoming energy that is reflected (Bond albedo) and the fraction of absorbed energy radiated by the night hemisphere (heat recirculation efficiency). Many empirical studies of hot Jupiters have implicitly assumed that the dayside is the hottest hemisphere and the nightside is the coldest hemisphere. For a given eclipse depth and phase amplitude, an orbital lag between a planet's peak brightness and its eclipse---a phase offset---implies that planet's nightside emits greater flux. To quantify how phase offsets impact the energy budgets of short-period planets, we compile all infrared observations of the nine planets with multi-band eclipse depths and phase curves. Accounting for phase offsets shifts planets to lower Bond albedo and greater day--night heat transport, usually by $\lesssim 1\sigma$. For WASP-12b, the published phase variations have been analyzed in two different ways, and the inferred energy budget depends sensitively on which analysis one adopts.\ Our fiducial scenario supports a Bond albedo of $0.27^{+0.12}_{-0.13}$, significantly higher than the published optical geometric albedo, and a recirculation efficiency of $0.03^{+0.07}_{-0.02}$, following the trend of larger day--night temperature contrast with greater stellar irradiation. If instead we adopt the alternative analysis, then WASP-12b has a Bond albedo consistent with zero and a much higher recirculation efficiency. To definitively determine the energy budget of WASP-12b, new observational analyses will be necessary., Comment: 11 pages, 6 figures, 2 tables; Accepted in ApJ. Minor corrections to text/references. (Improved treatment of stellar brightness temperatures and systematic uncertainties; better discussion of WASP-12b; updated Figures 2, 3, and 5, and fit parameters in Tables 1 and 2)

One of the primary questions when characterizing Earth-sized and super-Earth-sized exoplanets is whether they have a substantial atmosphere like Earth and Venus or a bare-rock surface like Mercury. Phase curves of the planets in thermal emission provide clues to this question, because a substantial atmosphere would transport heat more efficiently than a bare-rock surface. Analyzing phase curve photometric data around secondary eclipse has previously been used to study energy transport in the atmospheres of hot Jupiters. Here we use phase curve, Spitzer time-series photometry to study the thermal emission properties of the super-Earth exoplanet 55 Cancri e. We utilize a semi-analytical framework to fit a physical model to the infrared photometric data at 4.5 micron. The model uses parameters of planetary properties including Bond albedo, heat redistribution efficiency (i.e., ratio between radiative timescale and advective timescale of the atmosphere), and atmospheric greenhouse factor. The phase curve of 55 Cancri e is dominated by thermal emission with an eastward-shifted hot spot. We determine the heat redistribution efficiency to be ~1.47, which implies that the advective timescale is on the same order as the radiative timescale. This requirement cannot be met by the bare-rock planet scenario because heat transport by currents of molten lava would be too slow. The phase curve thus favors the scenario with a substantial atmosphere. Our constraints on the heat redistribution efficiency translate to an atmospheric pressure of ~1.4 bar. The Spitzer 4.5-micron band is thus a window into the deep atmosphere of the planet 55 Cancri e., Comment: Accepted for publication on AJ

Maps of the variation in bolometric Bond albedo and thermal inertia across Rhea and Dione have been produced using various day and nighttime observations, taken by Cassini's Composite Infrared Spectrometer (CIRS). The albedo maps show the same trend that has been previously observed on these satellites from reflected sunlight and thermal observations: a higher albedo on their leading hemispheres. The thermal inertia maps show two previously unknown anomalous high thermal inertia regions: low latitudes on Dione's leading hemisphere and the bright ejecta blanket of Rhea's Inktomi crater. The thermal inertia on Dione increases modestly from a background value of 8 J m(exp -2) K(exp -1) s(exp -1/2) to 11 J m(exp -)2 K(exp -1) s(exp -1/2) in a region preferentially bombarded by high-energy electrons. We believe that this region on Dione is probably analogous to the thermally anomalous regions recently discovered at equivalent locations on Mimas and Tethys, dubbed Pac-Man features. The smaller magnitude of Dione's thermal anomaly, compared to that of Mimas and Tethys, provides additional evidence that surface alteration by high-energy electrons produces these anomalies, as the high-energy electron flux decreases with increasing distance from Saturn. However, unlike on Mimas and Tethys, the thermally anomalous region on Dione does not display a spatially correlated decrease in the IR/UV (0.930 micrometers/0.338 micrometers) color ratio, implying that the minimum electron energy threshold of the IR/UV ice surface darkening mechanism is not met on Dione. The average of the mapped bolometric Bond albedos on the leading and trailing hemispheres of Dione are 0.49 +/- 0.11 and 0.44 +/- 0.13 respectively. On Rhea the thermal inertia increases from a background value of 10 J m(exp -2) K (exp -1) s(exp -1/2) to 19 J m(exp -2) K(exp -1) s(exp -1/2) on the ejecta blanket of Inktomi crater, probably due to a mixture of high and low thermal inertia material in this region, such as a fine- and large-grain ice mixture. This is the first time a thermal inertia anomaly has been associated with an impact crater on an icy saturnian satellite. On Rhea the average of the mapped bolometric Bond albedos on the leading and trailing hemispheres are 0.59 +/- 0.11 and 0.56 +/- 0.13 respectively.

We present two numerical models designed to reproduce the temperatures of the illuminated Saturn rings as measured by the CASSINI-CIRS instrument. Our models are constrained by all available temperature measurements performed on the illuminated rings since SOI. Both models reproduce well the variations of temperature under any illumination and observation geometry. One model is derived from a purely numerical data mining approach, relying on the implementation of a Neural Network that treats the data set globally. This model is used as a test of coverage completeness of the observational parameter space, driving our ability to characterize the rings thermal response. The second (analytical) model is derived using simple physical considerations, by treating the rings as a surface rather than as a collection of individual particles, combined with an empirical anisotropy function to describe the temperature resulting from the Sun’s and Saturn’s heating. The thermal response of this ring-surface is parameterized by its Bond albedo and emissivity, thermal relaxation time and a set of geometrical parameters quantifying the anisotropy of the temperature measurements depending on azimuth and elevation of the observer with respect to the ring plane, as well as on the solar elevation. Both models provide formulae to predict the ring temperature, that will ease the benchmarking of future physical models against data. The physical model is applied to fit the temperature of tens of different radial slices, allowing us to constrain the combined emissivity and albedo, thermal relaxation time and anisotropy parameters of the ring slabs with the highest radial resolution achieved so far with CIRS. Using for the first time all observation geometries available for illuminated rings, we are confident that our values are as unbiased as possible against observation geometry. The thermal relaxation time appears to be short, a few tens of minutes, and independent of the radial distance across the whole ring. A study of the temperature anisotropy suggests inter-particle shadowing is important in the B ring and in the outer A ring regions.

The airborne spectral observations of the upward and downward irradiances are revisited to investigate the dependence of the near-ground albedo as a function of wavelength in the entire solar spectrum for different surfaces (sand, water, snow) and under different conditions (clear or cloudy sky). The radiative upward and downward fluxes were determined by a diffraction spectrometer flown on a research aircraft that was performing multiple flight paths near the ground. The results obtained show that the near-ground albedo does not generally increase with increasing wavelengths for all kinds of surfaces as is widely believed today. Particularly, in the case of water surfaces it was found that the albedo in the ultraviolet region is more or less independent of the wavelength on a long-term basis. Interestingly, in the visible and near-infrared spectra the water albedo obeys an almost constant power-law relationship with wavelength. In the case of sand surfaces it was found that the sand albedo is a quadratic function of wavelength, which becomes more accurate if the ultraviolet wavelengths are neglected. Finally, it was found that the spectral dependence of snow albedo behaves similarly to that of water, i.e. both decrease from the ultraviolet to the near-infrared wavelengths by 20–50%, despite the fact that their values differ by one order of magnitude (water albedo being lower). In addition, the snow albedo vs. ultraviolet wavelength is almost constant, while in the visible near-infrared spectrum the best simulation is achieved by a second-order polynomial, as in the case of sand, but with opposite slopes.

We have studied the photometric properties of the nucleus of a hyperactive comet, 103P/Hartley 2, at visible wavelengths using the DIXI flyby images with both disk-integrated and disk-resolved analyses. The disk-integrated phase function of the nucleus has a linear slope of 0.046 ± 0.002 mag/deg and an absolute magnitude of 18.4 ± 0.1 at V -band. The nucleus displays an overall linear, featureless spectrum between 400 nm and 850 nm. The linear spectral slope is 7.6 ± 3.6% per 100 nm, corresponding to broadband solar-illuminated color indices B – V of 0.75 ± 0.05 and V – R of 0.43 ± 0.04. Disk-resolved photometric analysis with a Hapke model returns a best-fit single-scattering albedo of 0.036 ± 0.006, an asymmetry factor of the Henyey–Greenstein single-particle phase function of −0.46 ± 0.06, and a photometric roughness of 15 ± 10°. The model yields a geometric albedo of 0.045 ± 0.009 and a Bond albedo of 0.012 ± 0.002. The overall photometric variations of the nucleus are small, with an equivalent albedo variation of 15% FWHM, and a color variation of 12% FWHM. Some areas near the terminator visible in the inbound images show an albedo of more than twice the global average value, and a much bluer color than the average nucleus. The overall photometric properties and variations of the nucleus of Hartley 2 are similar to those of the nuclei of Comets Wild 2 and Tempel 1 as studied from previous spacecraft flyby missions at similar resolutions.

The transiting planet CoRoT-1b is thought to belong to the pM-class of planets, in which the thermal emission dominates in the optical wavelengths. We present a detection of its secondary eclipse in the CoRoT white channel data, whose response function goes from ~400 to ~1000 nm. We used two different filtering approaches, and several methods to evaluate the significance of a detection of the secondary eclipse. We detect a secondary eclipse centered within 20 min at the expected times for a circular orbit, with a depth of 0.016+/-0.006%. The center of the eclipse is translated in a 1-sigma upper limit to the planet's eccentricity of ecosomega, 6 pages, to appear in A&A, submitted 18 march 2009, accepted 7 July 2009

Optical phase curves have become one of the common probes of exoplanetary atmospheres, but the information they encode has not been fully elucidated. Building on a diverse body of work, we upgrade the Flexible Modeling System (FMS) to include scattering in the two-stream, dual-band approximation and generate plausible, three-dimensional structures of irradiated atmospheres to study the radiative effects of aerosols or condensates. In the optical, we treat the scattering of starlight using a generalisation of Beer's law that allows for a finite Bond albedo to be prescribed. In the infrared, we implement the two-stream solutions and include scattering via an infrared scattering parameter. We present a suite of four-parameter general circulation models for Kepler-7b and demonstrate that its climatology is expected to be robust to variations in optical and infrared scattering. The westward and eastward shifts of the optical and infrared phase curves, respectively, are shown to be robust outcomes of the simulations. Assuming micron-sized particles and a simplified treatment of local brightness, we further show that the peak offset of the optical phase curve is sensitive to the composition of the aerosols or condensates. However, to within the measurement uncertainties, we cannot distinguish between aerosols made of silicates (enstatite or forsterite), iron, corundum or titanium oxide, based on a comparison to the measured peak offset ($41^\circ \pm 12^\circ$) of the optical phase curve of Kepler-7b. Measuring high-precision optical phase curves will provide important constraints on the atmospheres of cloudy exoplanets and reduce degeneracies in interpreting their infrared spectra., Comment: 11 pages, 8 figures, 1 table. Accepted by MNRAS

New Cassini observations show that Saturn’s moon Tethys, like Mimas, has a region of anomalously high thermal inertia at low latitudes centered on its leading hemisphere. Derivation of the thermophysical properties of the surface across three representative regions indicate that the bolometric Bond albedo across Tethys’ leading hemisphere remains constant, within error, whilst the thermal inertia increases from 5 ± 1 J s −1 m −1 K −1 outside of the anomalous region to 25 ± 3 J s −1 m −1 K −1 inside. The thermally anomalous region is spatially correlated with a decrease in the IR/UV surface coloration. The discovery greatly strengthens the hypothesis that high-energy electrons, which preferentially bombard the leading hemispheres on both satellites, produce dramatic alterations in surface texture.

Io’s sublimation atmosphere is inextricably linked to the SO 2 surface frost temperature distribution which is poorly constrained by observations. We constrain Io’s surface thermal distribution by a parametric study of its thermophysical properties in an attempt to better model the morphology of Io’s sublimation atmosphere. Io’s surface thermal distribution is represented by three thermal units: sulfur dioxide (SO 2 ) frosts/ices, non-frosts (probably sulfur allotropes and/or pyroclastic dusts), and hot spots. The hot spots included in our thermal model are static high temperature surfaces with areas and temperatures based on Keck infrared observations. Elsewhere, over frosts and non-frosts, our thermal model solves the one-dimensional heat conduction equation in depth into Io’s surface and includes the effects of eclipse by Jupiter, radiation from Jupiter, and latent heat of sublimation and condensation. The best fit parameters for the SO 2 frost and non-frost units are found by using a least-squares method and fitting to observations of the Hubble Space Telescope’s Space Telescope Imaging Spectrograph (HST STIS) mid- to near-UV reflectance spectra and Galileo PPR brightness temperature. The thermophysical parameters are the frost Bond albedo, α F , and thermal inertia, Γ F , as well as the non-frost surface Bond albedo, α NF , and thermal inertia, Γ NF . The best fit parameters are found to be α F ≈ 0.55 ± 0.02 and Γ F ≈ 200 ± 50 J m −2 K −1 s −1/2 for the SO 2 frost surface and α NF ≈ 0.49 ± 0.02 and Γ NF ≈ 20 ± 10 J m −2 K −1 s −1/2 for the non-frost surface. These surface thermophysical parameters are then used as boundary conditions in global atmospheric simulations of Io’s sublimation-driven atmosphere using the direct simulation Monte Carlo (DSMC) method. These simulations are unsteady, three-dimensional, parallelized across 360 processors, and include the following physical effects: inhomogeneous surface frosts, plasma heating, and a temperature-dependent residence time on the non-frost surface. The DSMC simulations show that the sub-jovian hemisphere is significantly affected by the daily solar eclipse. The simulated SO 2 surface frost temperature is found to drop only ∼5 K during eclipse due to the high thermal inertia of SO 2 surface frosts but the SO 2 gas column density falls by a factor of 20 compared to the pre-eclipse column due to the exponential dependence of the SO 2 vapor pressure on the SO 2 surface frost temperature. Supersonic winds exist prior to eclipse but become subsonic during eclipse because the collapse of the atmosphere significantly decreases the day-to-night pressure gradient that drives the winds. Prior to eclipse, the supersonic winds condense on and near the cold nightside and form a highly non-equilibrium oblique shock near the dawn terminator. In eclipse, no shock exists since the gas is subsonic and the shock only reestablishes itself an hour or more after egress from eclipse. Furthermore, the excess gas that condenses on the non-frost surface during eclipse leads to an enhancement of the atmosphere near dawn. The dawn atmospheric enhancement drives winds that oppose those that are driven away from the peak pressure region above the warmest area of the SO 2 frost surface. These opposing winds meet and are collisional enough to form stagnation point flow. The simulations are compared to Lyman-α observations in an attempt to explain the asymmetry between the dayside atmospheres of the anti-jovian and sub-jovian hemispheres. Lyman-α observations indicate that the anti-jovian hemisphere has higher column densities than the sub-jovian hemisphere and also has a larger latitudinal extent. A composite “average dayside atmosphere” is formed from a collisionless simulation of Io’s atmosphere throughout an entire orbit. This composite “average dayside” atmosphere without the effect of global winds indicates that the sub-jovian hemisphere has lower average column densities than the anti-jovian hemisphere (with the strongest effect at the sub-jovian point) due primarily to the diurnally averaged effect of eclipse. This is in qualitative agreement with the sub-jovian/anti-jovian asymmetry in the Lyman-α observations which were alternatively explained by the bias of volcanic centers on the anti-jovian hemisphere. Lastly, the column densities in the simulated average dayside atmosphere agree with those inferred from Lyman-α observations despite the thermophysical parameters being constrained by mid- to near UV observations which show much higher instantaneous SO 2 gas column densities. This may resolve the apparent discrepancy between the lower “average dayside” column densities observed in the Lyman-α and the higher instantaneous column densities observed in the mid- to near UV.

The single scattering albedo omega(sub O lambda) in atmospheric radiative transfer is the ratio of the scattering coefficient to the extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, thus the single scattering albedo, are functions of wavelength lambda and droplet size r. This note shows that for water droplets at weakly absorbing wavelengths, the ratio omega(sub O lambda)(r)/omega(sub O lambda)(r (sub O)) of two single scattering albedo spectra is a linear function of omega(sub O lambda)(r). The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo spectrum omega(sub O lambda)(r) via one known spectrum omega(sub O lambda)(r (sub O)). We provide a simple physical explanation of the discovered relationship. Similar linear relationships were found for the single scattering albedo spectra of non-spherical ice crystals.

Titan shows its surface through many methane windows in the 1–5 μ m region. Windows at shorter wavelengths also exist, polluted by scattering off of atmospheric haze that reduces the surface contrast. At visible wavelengths, the surface of Titan has been observed by Voyager I, the Hubble Space Telescope, and ground-based telescopes. We present here global surface mapping of Titan using the visible wavelength channels from Cassini's Visual and Infrared Mapping Spectrometer (VIMS). We show global maps in each of the VIMS-V channels extending from 0.35 to 1.05 μ m . We find methane windows at 0.637, 0.681, 0.754, 0.827, 0.937, and 1.046 μ m and apply an RGB color scheme to the 0.754, 0.827 and 0.937 μ m windows to search for surface albedo variations. Our results show that Titan appears gray at visible wavelengths; hence scattering albedo is a good approximation of the Bond albedo. Maps of this genre have already been made and published using the infrared channels of VIMS. Ours are the first global maps of Titan shortward of 0.938 μ m . We compare the older IR maps to the new VIMS-V maps to constrain surface composition. For instance Tui Regio and Hotei Regio, referred to as 5 ‐ μ m bright spots in previous papers, do not distinguish themselves at all visible wavelengths. The distinction between the dune areas and the bright albedo spots, however, such as the difference between Xanadu and Senkyo, is easily discernible. We employ an empirically derived algorithm to remove haze layers from Titan, revealing a better look at the surface contrast.

Motivated by the work of Guillot, we present a semi-analytical formalism for calculating the temperature–pressure profiles in hot Jovian atmospheres which includes the effects of clouds/hazes and collision-induced absorption. Using the dual-band approximation, we assume that stellar irradiation and thermal emission from the hot Jupiter occur at distinct wavelengths (‘shortwave’ versus ‘longwave’). For a purely absorbing cloud/haze, we demonstrate its dual effect of cooling and warming the upper and lower atmosphere, respectively, which modifies, in a non-trivial manner, the condition for whether a temperature inversion is present in the upper atmosphere. The warming effect becomes more pronounced as the cloud/haze deck resides at greater depths. If it sits below the shortwave photosphere, the warming effect becomes either more subdued or ceases altogether. If shortwave scattering is present, its dual effect is to warm and cool the upper and lower atmospheres, respectively, thus counteracting the effects of enhanced longwave absorption by the cloud/haze. We make a tentative comparison of a four-parameter model to the temperature–pressure data points inferred from the observations of HD 189733b and estimate that its Bond albedo is approximately 10 per cent. Besides their utility in developing physical intuition, our semi-analytical models are a guide for the parameter space exploration of hot Jovian atmospheres via three-dimensional simulations of atmospheric circulation.

The energy balance at the surface of an airless planetary body is strongly influenced by the bolometric Bond albedo and the surface thermal inertia. Both of these values may be calculated through the application of a thermal model to measured surface temperatures. The accuracy of either, though, increases if the value of the other is better constrained. In this study, we used the improved global bolometric Bond albedo map of Iapetus derived from Cassini VIMS and ISS and Voyager ISS data in conjunction with Cassini CIRS temperature data to reevaluate surface thermal inertia across Iapetus. Results showed the thermal inertia of the dark terrain varies between 11 and 14.8 J m −2 K −1 s −1/2 while the light material varies between 15 and 25 J m −2 K −1 s −1/2 . Using an approximation to the thermal properties of the dark overburden derived from our thermal inertia results, we can implement our thermal model to provide estimates on the dark material thickness, which was found to lie between 7 cm and 16 cm. In order to develop an accurate global thermal model, a weighted function that approximates the surface thermal inertia across Iapetus was developed and verified via our measurements. The global bolometric Bond albedo map, surface thermal inertia map, and the thermal model are then used to synthesize global temperature maps that may be used to study the stability of volatiles.

We present a new rough-surface thermophysical model (Advanced Thermophysical Model or ATPM) that describes the observed directional thermal emission from any atmosphereless planetary surface. It explicitly incorporates partial shadowing, scattering of sunlight, self-heating and thermal–infrared beaming (re-radiation of absorbed sunlight back towards the Sun as a result of surface roughness). The model is verified by accurately reproducing ground-based directional thermal emission measurements of the lunar surface using surface properties that are consistent with the findings of the Apollo missions and roughness characterized by an rms slope of ∼32°. By considering the wide range of potential asteroid surface properties, the model implies a beaming effect that cannot be described by a simple parameter or function. It is highly dependent on the illumination and viewing angles as well as surface thermal properties and is predominantly caused by macroscopic rather than microscopic roughness. Roughness alters the effective Bond albedo and thermal inertia of the surface as well as moving the mean emission away from the surface normal. For accurate determination of surface properties from thermal–infrared observations of unresolved bodies or resolved surface elements, roughness must be explicitly modelled, preferably aided with thermal measurements at different emission angles and wavelengths.

A 2-month series of quasi-simultaneous imaging photometric observations of the Moon and the Sun has been performed at Maidanak Observatory (Uzbekistan). New absolute values of lunar albedo have been obtained. Maps of lunar apparent albedo and equigonal albedo at phase angles 1.7–73° at wavelength 603 nm are presented. The standard deviation of our data from a best-fitted phase curve is 2%. The average ratio of the Clementine albedo to ours is 1.41. While the ratio of ROLO albedo to ours is 0.87, our data are in agreement with independent measurements of absolute albedo by Saiki et al. (Saiki, K., Saito, K., Okuno, H., Suzuki, A., Yamanoi, Y., Hirata N., Nakamura, R. [2008]. Earth Planets Space 60, 417–424) at a phase angle near 7°. A phase ratio imaging near opposition (1.6°/2.7°) shows almost the same ratio for maria and highlands, though bright craters (e.g., Tycho, Copernicus, Aristarchus) clearly reveal smaller slopes of phase function. This is an unexpected result, as the craters are bright and one could anticipate a manifestation of the coherent backscattering effect resulting in the opposition spike increasing at so small phase angles.

The authors describe a new radiative transfer model of the Venus atmosphere (RTM) that includes optical properties from nine gases and four cloud modes between 0.1 and 260 μm. A multiple-stream discrete ordinate flux solver is used to calculate solar and atmospheric infrared fluxes with a prescribed temperature profiles and calculate radiative–convective equilibrium temperatures using the model. Components of the RTM are validated using observations from Pioneer Venus and Venus Express. A visible bond albedo of 0.74 and subsolar surface visible flux of 50 W m−2 [4.0% of the top-of-atmosphere (TOA) insolation] are calculated for a suitable temperature and composition profile derived from the Venus International Reference Atmosphere. Solar fluxes are simulated over a range of latitudes and good agreement is found with results from the Pioneer Venus probes and Venera landers. TOA infrared fluxes are compared with Venus Express observations and found to compare well at all observed wavelengths. The RTM is used to calculate radiative heating rates and these calculated heating rates are compared with those prescribed in a modern Venus GCM. Modifications are suggested to improve the prescribed thermal forcing used in recent GCMs. Using a small family of numerical and physical configurations, little sensitivity to vertical resolution is found in the model. For suitable global mean solar forcing a surface temperature of 750 K at radiative–convective equilibrium is calculated, in good agreement with observations and other recent modeling efforts.