Trapping and dynamical evolution of interplanetary dust particles in Earth’s quasi-satellite resonance

We used numerical simulations to model the orbital evolution of interplanetary dust particles (IDPs) evolving inward past Earth’s orbit under the influence of radiation pressure, Poynting–Robertson light drag (PR drag), solar wind drag, and gravitational perturbations from the planets. A series of β values (where β is the ratio of the force from radiation pressure to that of central gravity) were used ranging from 0.0025 up to 0.02. Assuming a composition consistent with astronomical silicate and a particle density of 2.5 g cm⁻³ these β values correspond to dust particle diameters ranging from 200 μm down to 25 μm. As the dust particle orbits decay past 1 AU between 4% (for β = 0.02, or 25 μm) and 40% (for β = 0.0025, or 200 μm) of the population became trapped in 1:1 co-orbital resonance with Earth. In addition to traditional horseshoe type co-orbitals, we found about a quarter of the co-orbital IDPs became trapped as so-called quasi-satellites. Quasi-satellite IDPs always remain relatively near to Earth (within 0.1–0.3 AU, or 10–30 Hill radii, R_H) and undergo two close-encounters with Earth each year. While resonant perturbations from Earth halt the decay in semi-major axis of quasi-satellite IDPs their orbital eccentricities continue to decrease under the influence of PR drag and solar wind drag, forcing the IDPs onto more Earth-like orbits. This has dramatic consequences for the relative velocity and distance of closest approach between Earth and the quasi-satellite IDPs. After 10⁴–10⁵ years in the quasi-satellite resonance dust particles are typically less than 10R_H from Earth and consistently coming within about 3R_H. In the late stages of evolution, as the dust particles are escaping the 1:1 resonance, quasi-satellite IDPs can have deep close-encounters with Earth significantly below R_H. Removing the effects of Earth’s gravitational acceleration reveals that encounter velocities (i.e., velocities “at infinity”) between quasi-satellite IDPs and Earth during these close-encounters are just a few hundred meters per second or slower, well below the average values of 2–4 km s⁻¹ for non-resonant Earth-crossing IDPs with similar initial orbits. These low encounter velocities lead to a factor of 10–100 increase in Earth’s gravitationally enhanced impact cross-section (σ_grav) for quasi-satellite IDPs compared to similar non-resonant IDPs. The enhancement in σ_grav between quasi-satellite IDPs and cometary Earth-crossing IDPs is even more pronounced, favoring accretion of quasi-satellite dust particles by a factor of 100–3000 over the cometary IDPs. This suggests that quasi-satellite dust particles may dominate the flux of large (25–200 μm) IDPs entering Earth’s atmosphere. Furthermore, because quasi-satellite trapping is known to be directly correlated with the host planet’s orbital eccentricity the accretion of quasi-satellite dust likely ebbs and flows on 10⁵ year time scales synchronized with Earth’s orbital evolution.     

Absolute properties of the neglected eclipsing B-type binary HD 194495

We present the results of the high-resolution spectroscopic observations of the neglected binary system HD 194495 (B3 IV-V+B4 V). A combined analysis of three different photometric data set (Tycho B_T and V_T photometry, H_p-band data of Hipparcos and V-band data of ASAS3 photometry) and radial velocities indicates that the system has an orbital period of 4.90494 ± 0.00005 days and an inclination of 69 ± 1 degrees. This solution yields masses and radii of M₁ = 7.57 ± 0.08 M_⊙ and R₁ = 5.82 ± 0.03 R_⊙ for the primary and M₂ = 5.46 ± 0.09 M_⊙ and R₂ = 3.14 ± 0.08 R_⊙ for the secondary. Based on the position of the two stars plotted on a theoretical H-R diagram, we find that the age of the system is ?28 Myr, according to stellar evolutionary models. The spectroscopic and photometric results are in agreement with those obtained using theoretical predictions.     

Convective transport of low energy cosmic rays in the interplanetary region

A complete solution has been obtained of the steady-state transport equations, including energy losses, for cosmic-rays in the interplanetary region for conditions in which diffusive transport is negligible and convective effects dominate. The region of validity of the solution will in general be a shell between heliocentric radiiR ₁ andR ₂ (R ₂ may be infinite). The precise range of kinetic energyT and heliocentric radiusr in which the solution is valid is not known but it appears to be applicable in the vicinity of Earth to protons withT≤1 MeV. ForT～0.5 MeV near Earth,R ₁ may be ～0.5 AU andR ₁ will decrease asT, observed near Earth, decreases. The solution is simple in form but quite general; it predicts the differential number densityU (r, T) in terms of that observed at radius a (near Earth, say). Thus it may be quite useful in interpreting and co-ordinating steady-state cosmicray observations atT～1 MeV. The differential and integral intensities, differential anisotropy and differential radial-gradient at (r, T) also are determined. A simple interpretation of the solution is given in terms of energy losses due to adiabatic deceleration of the particles as they are being convected outward from the Sun. This leads to the useful notion of following a particle in (r, T) as it increasesr and decreasesT. Particles convected from the outer corona to Earth decrease their kinetic energy by factor ～500.Following a particle the Compton-Getting factor remains constant. Particles observed at (a, T) in convective transport have come from nearer the Sun; they may be of solar origin but may also be of galactic origin having penetrated tor<R ₁   

Keck and VLT AO observations and models of the uranian rings during the 2007 ring plane crossings

We present observations of the uranian ring system at a wavelength of 2.2 μm, taken between 2003 and 2008 with NIRC2 on the W.M. Keck telescope in Hawaii, and on 15–17 August 2007 with NaCo on the Very Large Telescope (VLT) in Chile. Of particular interest are the data taken around the time of the uranian ring plane crossing with Earth on 16 August 2007, and with the Sun (equinox) on 7 December 2007. We model the data at the different viewing aspects with a Monte Carlo model to determine: (1) the normal optical depth τ₀, the location, and the radial extent of the main rings, and (2) the parameter Aτ₀ (A is the particle geometric albedo), the location, and the radial plus vertical extent of the dusty rings. Our main conclusions are: (i) The brightness of the ? ring is significantly enhanced at small phase and ring inclination angles; we suggest this extreme opposition effect to probably be dominated by a reduction in interparticle shadowing. (ii) A broad sheet of dust particles extends inwards from the λ ring almost to the planet itself. This dust sheet has a vertical extent of ∼140 km, and Aτ₀ = 2.2 × 10⁻⁶. (iii) The dusty rings between ring 4 and the α ring and between the α and β rings are vertically extended with a thickness of ∼300 km. (iv) The ζ ring extends from ∼41,350 km almost all the way inwards to the planet. The main ζ ring, centered at ∼39,500 km from the planet, is characterized by Aτ₀ = 3.7 × 10⁻⁶; this parameter decreases closer to the planet. The ζ ring has a full vertical extent of order 800–900 km, with a pronounced density enhancement in the mid-plane. (v) The η_c ring is optically thin and less than several tens of km in the vertical direction. This ring may be composed of macroscopic material, surrounded by clumps of dust.     

Micrometer-sized ice particles for planetary-science experiments - I. Preparation, critical rolling friction force, and specific surface energy

Coagulation models assume a higher sticking threshold for micrometer-sized ice particles than for micrometer-sized silicate particles. However, in contrast to silicates, laboratory investigations of the collision properties of micrometer-sized ice particles (in particular, of the most abundant H₂O-ice) have not been conducted yet. Thus, we used two different experimental methods to produce micrometer-sized H₂O-ice particles, i.e. by spraying H₂O droplets into liquid nitrogen and by spraying H₂O droplets into a cold nitrogen atmosphere. The mean particle radii of the ice particles produced with these experimental methods are (1.49 ± 0.79) μm and (1.45 ± 0.65) μm. Ice aggregates composed of the micrometer-sized ice particles are highly porous (volume filling factor: ? = 0.11 ± 0.01) or rather compact (volume filling factor: ? = 0.72 ± 0.04), depending on the method of production. Furthermore, the critical rolling friction force of F_Roll,ice = (114.8 ± 23.8) × 10⁻¹⁰ N was measured for micrometer-sized ice particles, which exceeds the critical rolling friction force of micrometer-sized SiO₂ particles . This result implies that the adhesive bonding between micrometer-sized ice particles is stronger than the bonding strength between SiO₂ particles. An estimation of the specific surface energy of micrometer-sized ice particles, derived from the measured critical rolling friction forces and the surface energy of micrometer-sized SiO₂ particles, results in γ_ice = 0.190 J m⁻².     

Radar detection of Asteroid 2002 AA29

Radar echoes from Earth co-orbital Asteroid 2002 AA29 yield a total-power radar cross section of 2.9×10⁻⁵ km² ±25%, a circular polarization ratio of SC/OC=0.26±0.07, and an echo bandwidth of at least 1.5 Hz. Combining these results with the estimate of its visual absolute magnitude, H_V=25.23±0.24, from reported Spacewatch photometry indicates an effective diameter of 25±5 m, a rotation period no longer than 33 min, and an average surface bulk density no larger than 2.0 g cm⁻³; the asteroid is radar dark and optically bright, and its statistically most likely spectral class is S. The H_V estimate from LINEAR photometry (23.58±0.38) is not compatible with either Spacewatch's H_V or our radar results. If a bias this large were generally present in LINEAR's estimates of H_V for asteroids it has discovered or observed, then estimates of the current completeness of the Spaceguard Survey would have to be revised downward.     

Solar activity variations of thermospheric temperatures on Mars and a problem of CO in the lower atmosphere

Long-term MGS drag density observations at 390 km reveal variations of the density with season L_S (by a factor of 2) and solar activity index F_10.7 (by a factor of 3 for F_10.7 = 40-100). According to Forbes et al. (Forbes, J.M., Lemoine, F.G., Bruinsma, S.L., Smith, M.D., Zhang, X. [2008]. Geophys. Res. Lett. 35, L01201, doi:10.1029/2007GL031904), the variation with F_10.7 reflects variations of the exospheric temperature from 192 to 284 K. However, the derived temperature range corresponds to variation of the density at 390 km by a factor of 8, far above the observed factor of 3. The recent thermospheric GCMs agree with the derived temperatures but do not prove their adequacy to the MGS densities at 390 km. A model used by Forbes et al. neglects effects of eddy diffusion, chemistry and escape on species densities above 138 km. We have made a 1D-model of neutral and ion composition at 80-400 km that treats selfconsistently chemistry and transport of species with F_10.7, T_∞, and [CO₂]_80 km as input parameters. Applying this model to the MGS densities at 390 km, we find variation of T_∞ from 240 to 280 K for F_10.7 = 40 and 100, respectively. The results are compared with other observations and models. Temperatures from some observations and the latest models disagree with the MGS densities at low and mean solar activity. Linear fits to the exospheric temperatures are T_∞ = 122 + 2.17F_10.7 for the observations, T_∞ = 131 + 1.46F_10.7 for the latest models, and T_∞ = 233 + 0.54F_10.7 for the MGS densities at 390 km. Maybe the observed MGS densities are overestimated near solar minimum when they are low and difficult to measure. Seasonal variations of Mars’ thermosphere corrected for the varying heliocentric distance are mostly due to the density variations in the lower and middle atmosphere and weakly affect thermospheric temperature. Nonthermal escape processes for H, D, H₂, HD, and He are calculated for the solar minimum and maximum conditions.Another problem considered here refers to Mars global photochemistry in the lower and middle atmosphere. The models gave too low abundances of CO, smaller by an order of magnitude than those observed. Our current work shows that modifications in the boundary conditions proposed by Zahnle et al. (Zahnle, K., Haberle, R.M., Catling, D.C., Kasting, J.F. [2008]. J. Geophys. Res. 113, E11004, doi:10.1029/2008JE003160) are reasonable but do not help to solve the problem.     

A case for the standard model

The fundamental properties of Friedmann Universes, which are attractive because of their simplicity, are linear expansion (except for deceleration), cooling, and evolution. In addition it is assumed that the fundamental constants of physics are constant and that the known laws of physics apply (including GR). An increasing number of observational tests support these premises. In particular the expansion is as linear as can be tested. The present expansion rate (H ₀ = 55 ± 10 km s^-1 Mpc^-1) implies an expansion age of 17.8 ± 3.2 Gyr (forq ₀ = 0) to 11.9 ± 2.1 Gyr (for q₀ = 1/2). This agrees perfectly within the errors, even for a critical Universe, with present age determinations of the oldest objects in the Galaxy which require 13.5 ± 2 Gyr.     

A spatially resolved high spectral resolution study of Neptune’s stratosphere

Using TEXES, the Texas Echelon cross Echelle Spectrograph, mounted on the Gemini North 8-m telescope we have mapped the spatial variation of H₂, CH₄, C₂H₂ and C₂H₆ thermal-infrared emission of Neptune. These high-spectral-resolution, spatially resolved, thermal-infrared observations of Neptune offer a unique glimpse into the state of Neptune’s stratosphere in October 2007, L_S = 275.4° just past Neptune’s southern summer solstice (L_S = 270°). We use observations of the S(1) pure rotational line of molecular hydrogen and a portion of the ν₄ band of methane to retrieve detailed information on Neptune’s stratospheric vertical and meridional thermal structure. We find global-average temperatures of 163.8 ± 0.8, 155.0 ± 0.9, and 123.8 ± 0.8 K at the 7.0 × 10⁻³-, 0.12-, and 2.1-mbar levels with no meridional variations within the errors. We then use the inferred temperatures to model the emission of C₂H₂ and C₂H₆ in order to derive stratospheric volume mixing ratios (hence forth, VMR) as a function of pressure and latitude. There is a subtle meridional variation of the C₂H₂ VMR at the 0.5-mbar level with the peak abundance found at −28° latitude, falling off to the north and south. However, the observations are consistent within error to a meridionally constant C₂H₂ VMR of at 0.5 mbar. We find that the VMR of C₂H₆ at 1-mbar peaks at the equator and falls by a factor of 1.6 at −70° latitude. However, a meridionally constant VMR of at the 1-mbar level for C₂H₆ is also statistically consistent with the retrievals. Temperature predictions from a radiative-seasonal climate model of Neptune that assumes the hydrocarbon abundances inferred in this paper are lower than the measured temperatures by 40 K at 7 × 10⁻³ mbar, 30 K at 0.12 mbar and 25 K at 2.1 mbar. The radiative-seasonal model also predicts meridional temperature variations on the order of 10 K from equator to pole, which are not observed. Assuming higher stratospheric CH₄ abundance at the equator relative to the south pole would bring the meridional trends of the inferred temperatures and radiative-seasonal model into closer agreement.We have also retrieved observations of C₂H₄ emission from Neptune’s stratosphere using TEXES on the NASA Infrared Telescope Facility (IRTF) in June 2003, L_S = 266°. Using the observations from the middle of the planet and an average of the middle three latitude temperature profiles from the 2007 observations (9.5° of L_S later, the seasonal equivalent of 9.5 Earth days within Earth’s seasonal cycle), we infer a C₂H₄ VMR of at 1.5 × 10⁻³ mbar, a value that is 3.25 times that predicted by global-average photochemical models.     

Sublimation’s impact on temporal change of albedo dichotomy on Iapetus

Iapetus, one of the saturnian moons, has an extreme albedo contrast between the leading and trailing hemispheres. The origin of this albedo dichotomy has led to several hypotheses, however it remains controversial. To clarify the origin of the dichotomy, the key approach is to investigate the detailed distribution of the dark material. Recent studies of impact craters and surface temperature from Cassini spacecraft data implied that sublimation of H₂O ice can occur on Iapetus’ surface. This ice sublimation can change the albedo distribution on the moon with time.In this study, we evaluate the effect of ice sublimation and simulate the temporal change of surface albedo. We assume the dark material and the bright ice on the surface to be uniformly mixed with a certain volume fraction, and the initial albedo distribution to incorporate the dark material deposits on the surface. That is, the albedo at the apex is lowest and concentrically increases in a sinusoidal pattern. This situation simulates that dark materials existed around the Iapetus’ orbit billions of years ago, and the synchronously rotating Iapetus swept the material and then deposited it on its surface. The evolution of the surface albedo during 4.0 Gyr is simulated by estimating the surface temperature from the insolation energy on Iapetus including the effect of Saturn’s eccentricity and Iapetus’ obliquity precession, and evaluating the sublimation rate of H₂O ice from the Iapetus’ surface.As a result, we found that the distribution of the surface albedo changed dramatically after 4.0 Gyr of evolution. The sublimation has three important effects on the resultant surface albedo. First, the albedo in the leading hemisphere has significantly decreased to approach the minimum value. Second, the albedo distribution has been elongated along the equator. Third, the edge of the low albedo region has become clear. Considering the effect of ice sublimation, the current albedo distribution can be reconstructed from the sinusoidal albedo distribution, suggesting the apex-antapex cratering asymmetry as a candidate for the origin of the albedo dichotomy. From the model analysis, we obtained an important aspect that the depth of the turn-over layer where the darkening process proceeded for 4 Gyr should be an order of 10 cm, which is consistent with evaluation from the Cassini radar observations.     

Scalar field power-law cosmology with spatial curvature and dark energy-dark matter interaction

We consider a late closed universe of which scale factor is a power function of time using observational data from combined WMAP5+BAO+SN Ia dataset and WMAP5 dataset. The WMAP5 data give power-law exponent, α=1.01 agreeing with the previous study of H(z) data while combined data gives α=0.985. Considering a scalar field dark energy and dust fluid evolving in the power-law universe, we find field potential, field solution and equation of state parameters. Decaying from dark matter into dark energy is allowed in addition to the non-interaction case. Time scale characterizing domination of the kinematic expansion terms over the dust and curvature terms in the scalar field potential are found to be approximately 5.3 to 5.5 Gyr. The interaction affects in slightly lowering the height of scalar potential and slightly shifting potential curves rightwards to later time. Mass potential function of the interacting Lagrangian term is found to be exponentially decay function.     

A motivating exploration on lunar craters and low-energy dynamics in the Earth–Moon system

It is known that most of the craters on the surface of the Moon were created by the collision of minor bodies of the Solar System. Main Belt Asteroids, which can approach the terrestrial planets as a consequence of different types of resonance, are actually the main responsible for this phenomenon. Our aim is to investigate the impact distributions on the lunar surface that low-energy dynamics can provide. As a first approximation, we exploit the hyberbolic invariant manifolds associated with the central invariant manifold around the equilibrium point L ₂ of the Earth–Moon system within the framework of the Circular Restricted Three-Body Problem. Taking transit trajectories at several energy levels, we look for orbits intersecting the surface of the Moon and we attempt to define a relationship between longitude and latitude of arrival and lunar craters density. Then, we add the gravitational effect of the Sun by considering the Bicircular Restricted Four-Body Problem. In the former case, as main outcome, we observe a more relevant bombardment at the apex of the lunar surface, and a percentage of impact which is almost constant and whose value depends on the assumed Earth–Moon distance d_EM. In the latter, it seems that the Earth–Moon and Earth–Moon–Sun relative distances and the initial phase of the Sun θ ₀ play a crucial role on the impact distribution. The leading side focusing becomes more and more evident as d_EM decreases and there seems to exist values of θ ₀ more favorable to produce impacts with the Moon. Moreover, the presence of the Sun makes some trajectories to collide with the Earth. The corresponding quantity floats between 1 and 5 percent. As further exploration, we assume an uniform density of impact on the lunar surface, looking for the regions in the Earth–Moon neighbourhood these colliding trajectories have to come from. It turns out that low-energy ejecta originated from high-energy impacts are also responsible of the phenomenon we are considering.     

Close binary system GO Cyg

In this study, we present long term photometric variations of the close binary system GO Cyg. Modelling of the system shows that the primary is filling Roche lobe and the secondary of the system is almost filling its Roche lobe. The physical parameters of the system are M₁ = 3.0 ± 0.2M_⊙, M₂ = 1.3 ± 0.1M_⊙, R₁ = 2.50 ± 0.12R_⊙, R₂ = 1.75 ± 0.09R_⊙, L₁ = 64 ± 9L_⊙, L₂ = 4.9 ± 0.7L_⊙, and a = 5.5 ± 0.3R_⊙. Our results show that GO Cyg is the most massive system near contact binary (NCB). Analysis of times of the minima shows a sinusoidal variation with a period of 92.3 ± 0.5 yr due to a third body whose mass is less than 2.3M_⊙. Finally a period variation rate of −1.4 × 10⁻⁹ d/yr has been determined using all available light curves.     

Channel morphometry, sediment transport, and implications for tectonic activity and surficial ages of Titan basins

Fluvial features on Titan and drainage basins on Earth are remarkably similar despite differences in gravity and surface composition. We determined network bifurcation (R_b) ratios for five Titan and three terrestrial analog basins. Tectonically-modified Earth basins have R_b values greater than the expected range (3.0-5.0) for dendritic networks; comparisons with R_b values determined for Titan basins, in conjunction with similarities in network patterns, suggest that portions of Titan’s north polar region are modified by tectonic forces. Sufficient elevation data existed to calculate bed slope and potential fluvial sediment transport rates in at least one Titan basin, indicating that 75 mm water ice grains (observed at the Huygens landing site) should be readily entrained given sufficient flow depths of liquid hydrocarbons. Volumetric sediment transport estimates suggest that ∼6700-10,000 Titan years (∼2.0-3.0 × 10⁵ Earth years) are required to erode this basin to its minimum relief (assuming constant 1 m and 1.5 m flows); these lowering rates increase to ∼27,000-41,000 Titan years (∼8.0-12.0 × 10⁵ Earth years) when flows in the north polar region are restricted to summer months.     

Mass-radius curve for extrasolar Earth-like planets and ocean planets

By comparison with the Earth-like planets and the large icy satellites of the Solar System, one can model the internal structure of extrasolar planets. The input parameters are the composition of the star (Fe/Si and Mg/Si), the Mg content of the mantle (Mg# = Mg/[Mg + Fe]), the amount of H₂O and the total mass of the planet. Equation of State (EoS) of the different materials that are likely to be present within such planets have been obtained thanks to recent progress in high-pressure experiments. They are used to compute the planetary radius as a function of the total mass. Based on accretion models and data on planetary differentiation, the internal structure is likely to consist of an iron-rich core, a silicate mantle and an outer silicate crust resulting from magma formation in the mantle. The amount of H₂O and the surface temperature control the possibility for these planets to harbor an ocean. In preparation to the interpretation of the forthcoming data from the CNES led CoRoT (Convection Rotation and Transit) mission and from ground-based observations, this paper investigates the relationship between radius and mass. If H₂O is not an important component (less than 0.1%) of the total mass of the planet, then a relation (R/REarth)=ab(M/MEarth) is calculated with (a,b)=(1,0.306) and (a,b)=(1,0.274) for ¹⁰⁻²MEarth<M<MEarth and MEarth<M<10MEarth, respectively. Calculations for a planet that contains 50% H₂O suggest that the radius would be more than 25% larger than that based on the Earth-like model, with (a,b)=(1.258,0.302) for ¹⁰⁻²MEarth<M<MEarth and (a,b)=(1.262,0.275) for MEarth<M<10MEarth, respectively. For a surface temperature of 300 K, the thickness of the ocean varies from 150 to 50 km for planets 1 to 10 times the Earth's mass, respectively. Application of this algorithm to bodies of the Solar System provides not only a good fit to most terrestrial planets and large icy satellites, but also insights for discussing future observations of exoplanets.     

Secular light curve of 2P/Encke, a comet active at aphelion

We present the secular light curve of Comet 2P/Encke in two phase spaces, the log plot, and the time plot. The main conclusions of this work are: (a) The comet shows activity at perihelion and aphelion, caused by two different active areas: Source 1, close to the south pole, active at perihelion, and Source 2, at the north pole, centered at aphelion. (b) More than 18 physical parameters are measured from the secular light curves, many of them new, and are listed in the individual plots of the comet. Specifically we find for Source 1 the location of the turn on and turn off points of activity, RON=−1.63±0.03 AU, ROFF=+1.49±0.20 AU, TON=−87±5 d, TOFF=+94±15 d, the time lag, LAG(q)=6±1 d, the total active time, TACTIVITY=181±16 d, and the amplitude of the secular light curve, ASEC(1,1)=4.8±0.1 mag. (c) From this information the photometric age and the time-age defined in Ferrín [2005a. Icarus 178, 493-516; 2006. Icarus 185, 523-543], can be calculated, and we find P-AGE = 97 ± 8 comet years and T-AGE = 103 ± 9 comet years (cy). Thus Comet 2P/Encke is an old comet entering the methuselah stage (100 cy < age). (d) The activity at aphelion (Source 2), extends for TACTIVITY=815±30 d and the amplitude of the secular light curve is ASEC(1,Q)=3.0±0.2 mag. (e) From a new phase diagram an absolute magnitude and phase coefficient for the nucleus are determined, and we find RNUC(1,1,0)=15.05±0.14, and β=0.066±0.003. From this data we find a nucleus effective diameter DEFFE=5.12(+2.5;−1.7) km. These values are not much different from previous determinations but exhibit smaller errors. (f) The activity of Source 1 is due to H₂O sublimation because it shows curvature. The activity of Source 2 might also be due to H₂O due to the circumstantial situation that the poles point to the Sun at perihelion and aphelion. (g) We found a photometric anomaly at aphelion, with minimum brightness between +393 and +413 days after perihelion that may be an indication of topography. (h) We have re-reduced the 1858 secular light curve of Kamel [1991. Icarus 93, 226-245]. There are secular changes in 7 physical parameters, and we achieve for the first time, an absolute age calibration. We find that the comet entered the inner Solar System and began sublimating in 1645±40 AD. (i) It is concluded that the secular light curve can place constraints on the pole orientation of the nucleus of some comets, and we measure the ecliptic longitude of the south pole of 2P/Encke equal to 213.2±4.5°, in excellent agreement with other determinations of this parameter, but with smaller error. (j) Using the observed absolute magnitude of 1858 and 2003 and a suitable theoretical model, the extinction date of the comet is determined. We obtain ED=2056±3 AD, implying that the comet's lifetime is 125±12 revolutions about the Sun after entering the inner Solar System.     

Optical properties of Titan and early Earth haze laboratory analogs in the mid-visible

Scattering and absorption of sunlight by aerosols are integral to understanding the radiative balance of any planetary atmosphere covered in a haze, such as Titan and possibly the early Earth. One key optical parameter of an aerosol is its refractive index. We have simulated both Titan and early Earth organic haze aerosols in the laboratory and measured the real and imaginary portion of their refractive index at λ = 532 nm using cavity ringdown aerosol extinction spectroscopy. This novel technique allows analysis on freely-floating particles minutes after formation. For our Titan analog particles, we find a real refractive index of n = 1.35 ± 0.01 and an imaginary refractive index k = 0.023 ± 0.007, and for the early Earth analog particles we find n = 1.81 ± 0.02 and k = 0.055 ± 0.020. The Titan analog refractive index has a smaller real and similar imaginary refractive index compared to most previous laboratory measurements of Titan analog films, including values from Khare et al. (Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callcott, T.A., Williams, M.W. [1984]. Icarus 60, 127-137). These newly measured Titan analog values have implications for spacecraft retrievals of aerosol properties on Titan. The early Earth analog has a significantly higher real and imaginary refractive index than Titan analogs reported in the literature. These differences suggest that, for a given amount of aerosol, the early Earth analog would act as a stronger anti-greenhouse agent than the Titan analog.     

Titan’s atomic hydrogen corona

Based on measurements performed by the Hydrogen Deuterium Absorption Cell (HDAC) aboard the Cassini orbiter, Titan’s atomic hydrogen exosphere is investigated. Data obtained during the T9 encounter are used to infer the distribution of atomic hydrogen throughout Titan’s exosphere, as well as the exospheric temperature.The measurements performed during the flyby are modeled by performing Monte Carlo radiative transfer calculations of solar Lyman-α radiation, which is resonantly scattered on atomic hydrogen in Titan’s exosphere. Two different atomic hydrogen distribution models are applied to determine the best fitting density profile. One model is a static model that uses the Chamberlain formalism to calculate the distribution of atomic hydrogen throughout the exosphere, whereas the second model is a Particle model, which can also be applied to non-Maxwellian velocity distributions.The density distributions provided by both models are able to fit the measurements although both models differ at the exobase: best fitting exobase atomic hydrogen densities of n_H = (1.5 ± 0.5) × 10⁴ cm⁻³ and n_H = (7 ± 1) × 10⁴ cm⁻³ were found using the density distribution provided by both models, respectively. This is based on the fact that during the encounter, HDAC was sensitive to altitudes above about 3000 km, hence well above the exobase at about 1500 km. Above 3000 km, both models produce densities which are comparable, when taking into account the measurement uncertainty.The inferred exobase density using the Chamberlain profile is a factor of about 2.6 lower than the density obtained from Voyager 1 measurements and much lower than the values inferred from current photochemical models. However, when taking into account the higher solar activity during the Voyager flyby, this is consistent with the Voyager measurements. When using the density profile provided by the particle model, the best fitting exobase density is in perfect agreement with the densities inferred by current photochemical models.Furthermore, a best fitting exospheric temperature of atomic hydrogen in the range of T_H = (150-175) ± 25 K was obtained when assuming an isothermal exosphere for the calculations. The required exospheric temperature depends on the density distribution chosen. This result is within the temperature range determined by different instruments aboard Cassini. The inferred temperature is close to the critical temperature for atomic hydrogen, above which it can escape hydrodynamically after it diffused through the heavier background gas.     

Energy balance of the Deep Impact experiment

We present results on the energy balance of the Deep Impact experiment based on analysis of 180 infrared spectra of the ejecta obtained by the Deep Impact spacecraft. We derive an output energy of 16.5 (+9.1/−4.1) GJ. With an input energy of 19.7 GJ, the error bars are large enough so that there may or may not be a balance between the kinetic energy of the impact and that of outflowing materials. Although possible, no other source of energy other than the impactor or the Sun is needed to explain the observations. Most of the energy (85%) goes into the hot plume in the first few seconds, which only represents a very small fraction (<0.01%) of the total ejected mass. The hot plume contains 190 (+263/−71) kg of H₂O, 1.6 ± 0.5 kg of CO₂, 8.2 (+11.3/3.1) kg of CO (assuming a CO/H₂O ratio of 4.3%), 27.9 (+25.0/−8.9) kg of organic material and 255 ± 128 kg of dust, while the ejecta contains ∼10⁷ kg of materials. About 12% of the energy goes into the ejecta (mostly water) and 3% to destroy the impactor. Volatiles species other than H₂O (CO₂, CO or organic molecules) contribute to <7% of the energy balance. In terms of physical processes, 68% of the energy is used to accelerate grains (kinetic energy), 16% to heat them, 6% to sublimate or melt them and 10% (upper limit) to break and compress dust and/or water ice aggregates into small micron size particles. For the hot plume, we derive a dust/H₂O ratio of 1.3 (+1.9/−1.0), a CO₂/H₂O ratio of 0.008 (+0.009/−0.006), an organics/H₂O ratio of 0.15 (+0.29/−0.11) and an organics/dust ratio of 0.11 (+0.30/−0.07). This composition refers to the impact site and is different from that of the bulk nucleus, consistent with the idea of layers of different composition in the nucleus sub-surface. Our results emphasize the importance of laboratory impact experiments to understand the physical processes involved at such a large scale.     

Photometric observations of 107P/Wilson-Harrington

We present lightcurve observations and multiband photometry for 107P/Wilson-Harrington using five small- and medium-sized telescopes. The lightcurve has shown a periodicity of 0.2979 day (7.15 h) and 0.0993 day (2.38 h), which has a commensurability of 3:1. The physical properties of the lightcurve indicate two models: (1) 107P/Wilson-Harrington is a tumbling object with a sidereal rotation period of 0.2979 day and a precession period of 0.0993 day. The shape has a long axis mode (LAM) of L₁:L₂:L₃ = 1.0:1.0:1.6. The direction of the total rotational angular momentum is around λ = 310°, β = −10°, or λ = 132°, β = −17°. The nutation angle is approximately constant at 65°. (2) 107P/Wilson-Harrington is not a tumbler. The sidereal rotation period is 0.2979 day. The shape is nearly spherical but slightly hexagonal with a short axis mode (SAM) of L₁:L₂:L₃ = 1.5:1.5:1.0. The pole orientation is around λ = 330°, β = −27°. In addition, the model includes the possibility of binary hosting. For both models, the sense of rotation is retrograde. Furthermore, multiband photometry indicates that the taxonomy class of 107P/Wilson-Harrington is C-type. No clear rotational color variations are confirmed on the surface.