Ensemble properties of comets in the Sloan Digital Sky Survey

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u, g, r, i, z) photometry is used to determine the comets’ colors, sizes, surface brightness profiles, and rates of dust production in terms of the Afρ formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(<H) ∝ 10^{(0.49±0.05)H} for H < 18, with evidence of a much shallower fit N(<H) ∝ 10^{(0.19±0.03)H} for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 ± 0.05 in g ? r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.     

The effect of orbital evolution on the Haumea (2003 EL61) collisional family

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family’s age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50–400 m s^?1, approximately 20–45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m s^?1 of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ～50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1σ level, however both models predict more objects at larger relative velocities (>150 m s^?1) than have been identified.     

Constraints on mantle plumes on Venus: Implications for volatile history

The analysis of Venus’ gravity field and topography suggests the presence of a small number of deep mantle plumes (～9). This study predicts the number of plumes formed at the core–mantle boundary, their characteristics, and the production of partial melt from adiabatic decompression. Numerical simulations are performed using a 3D spherical code that includes large viscosity variations and internal heating. This study investigates the effect of several parameters including the core–mantle boundary temperature, the amount of internal heating, and the mantle viscosity. The smallest number of plumes is achieved when no internal heating is present. However, scaling Earth’s radiogenic heating to Venus suggests a value of ～16 TW. Cases with internal heating produce more realistic lid thickness and partial melting, but produce either too many plumes or no plumes if a high mantle temperature precludes the formation of a hot thermal boundary layer. Mantle viscosity must be reduced to at least 10²⁰ Pa s in order to include significant internal heating and still produce hot plumes. In all cases that predict melting, melting occurs throughout the upper mantle. Only cases with high core temperature (>1700 K) produce dry melting. Over time the upper mantle may have lost significant volatiles. Depending on the water content of the lower mantle, deep plumes may contribute to present-day atmospheric water via volcanic outgassing. Assuming 50 ppm water in mantle, 10 plumes with a buoyancy flux of 500 kg/s continuously erupting for 4 myr will outgas an amount of water on the order of that in the lower atmosphere. A higher level of internal heating than achieved to date, as well as relatively low mantle viscosity, may be required to achieve simulations with ～10 plumes and a thinner lid. Alternatively, if the mantle is heating up due to the stagnant lid, the effect is equivalent to having lower rates of internal heating. A temperature increase of 110 K/byr is equivalent to ?13 TW. This value along with the internal heating of 3 TW used in this study may represent the approximate heat budget of Venus’ mantle.     

Free collisions in a microgravity many-particle experiment. I. Dust aggregate sticking at low velocities

Over the past years the processes involved in the growth of planetesimals have extensively been studied in the laboratory. Based on these experiments, a dust-aggregate collision model was developed upon which computer simulations were based to evaluate how big protoplanetary dust aggregates can grow and to analyze which kinds of collisions are relevant in the solar nebula and are worth further studies in the laboratory. The sticking threshold velocity of millimeter-sized dust aggregates is one such critical value that have so far only theoretically been derived, as the relevant velocities could not be reached in the laboratory. We developed a microgravity experiment that allows us for the first time to study free collisions of mm-sized dust aggregates down to velocities of ～0.1 cm s^?1 to assess this part of the protoplanetary dust evolution model. Here, we present the results of 125 free collisions between dust aggregates of 0.5–2 mm diameter. Seven collisions with velocities between 0.2 and 3 cm s^?1 led to sticking, suggesting a transition from perfect sticking to perfect bouncing with a certain sticking probability instead of a sharp velocity threshold. We developed a model to explain the physical processes involved in dust-aggregate sticking, derived dynamical material properties of the dust aggregates from the results of the collisions, and deduced the velocity below which dust aggregates always stick. For millimeter-sized porous dust aggregates this velocity is 8 × 10^?5 m s^?1.     

Updated results of a search for main-belt comets using the Canada–France–Hawaii Telescope Legacy Survey

The results of a search for main-belt comets using Canada–France–Hawaii Telescope Legacy Survey data are updated. The remaining observations in the Very Wide segment of data, taken in the g’ or r’ filters, are visually inspected for cometary activity. The number of main-belt objects in the original and new data sets are 11,438 and 13,802, respectively, giving a total number of 25,240. This is the largest, and least biased, search for main-belt comets to date. One object is observed to show cometary activity, and a new upper limit for strongly active main-belt comets is derived to be 40 ± 18.     

Reassessing the formation of the inner Oort cloud in an embedded star cluster

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser, R., Duncan, M.J., Levison, H.F. [2006]. Icarus 184, 59–82) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster – either a Hernquist or Plummer potential – and chose other parameters based on the latest observations of embedded clusters from the literature. These other parameters result in the stars being on radial orbits and the cluster collapses. Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a ‘concentration radius’ for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters. The increasing trend implies that small star clusters form more compact inner Oort clouds than large clusters. We are unable to constrain the number of stars that resided in the cluster since most clusters yield inner Oort clouds that could be compatible with the current structure of the outer Solar System. The typical formation efficiency of the inner Oort cloud is 1.5%, significantly lower than previous estimates. We attribute this to the more violent dynamics that the Sun experiences as it rushes through the centre of the cluster during the latter’s initial phase of violent relaxation.     

Thermal evolution of Pluto and implications for surface tectonics and a subsurface ocean

Determining whether or not Pluto possesses, or once possessed, a subsurface ocean is crucial to understanding its astrobiological potential. In this study we use a 3D convection model to investigate Pluto’s thermal and spin evolution, and the present-day observational consequences of different evolutionary pathways. We test the sensitivity of our model results to different initial temperature profiles, initial spin periods, silicate potassium concentrations and ice reference viscosities. The ice reference viscosity is the primary factor controlling whether or not an ocean develops and whether that ocean survives to the present day. In most of our models present-day Pluto consists of a convective ice shell without an ocean. However if the reference viscosity is higher than 5 × 10¹⁵ Pa s, the shell will be conductive and an ocean should be present. For the nominal potassium concentration the present-day ocean and conductive shell thickness are both about 165 km; in conductive cases an ocean will be present unless the potassium content of the silicate mantle is less than 10% of its nominal value. If Pluto never developed an ocean, predominantly extensional surface tectonics should result, and a fossil rotational bulge will be present. For the cases which possess, or once possessed, an ocean, no fossil bulge should exist. A present-day ocean implies that compressional surface stresses should dominate, perhaps with minor recent extension. An ocean that formed and then re-froze should result in a roughly equal balance between (older) compressional and (younger) extensional features. These predictions may be tested by the New Horizons mission.     

The mass,orbit, and tidal evolution of the Quaoar–Weywot system

Here we present new adaptive optics observations of the Quaoar–Weywot system. With these new observations we determine an improved system orbit. Due to a 0.39 day alias that exists in available observations, four possible orbital solutions are available with periods of ～11.6, ～12.0, ～12.4, and ～12.8 days. From the possible orbital solutions, system masses of 1.3–1.5 ± 0.1 × 10²¹ kg are found. These observations provide an updated density for Quaoar of 2.7–5.0 g cm^?3. In all cases, Weywot’s orbit is eccentric, with possible values ～0.13–0.16. We present a reanalysis of the tidal orbital evolution of the Quaoar–Weywot system. We have found that Weywot has probably evolved to a state of synchronous rotation, and has likely preserved its initial inclination over the age of the Solar System. We find that for plausible values of the effective tidal dissipation factor tides produce a very slow evolution of Weywot’s eccentricity and semi-major axis. Accordingly, it appears that Weywot’s eccentricity likely did not tidally evolve to its current value from an initially circular orbit. Rather, it seems that some other mechanism has raised its eccentricity post-formation, or Weywot formed with a non-negligible eccentricity.     

Orbital period investigations of two short-period early-type overcontact binaries BH Cen and V701 Sco in two extremely young galactic clusters IC 2944 and NGC 6383

Both V701 Sco and BH Cen are two early-type short-period overcontact systems (P = 0.^d762 and P = 0.^d792, respectively). V701 Sco is a member of the young galactic cluster NGC 6383, while BH Cen is a component of a younger galactic cluster IC 2944 where star formation is in process. They provide good opportunity to understand the formation and evolution of binary stars. In the present paper, orbital period changes of the two binaries are investigated. It is discovered that the orbital period of BH Cen shows a long-term increase with a rate of dP/dt = +1.70(±0.39) × 10⁻⁷ days/year while it undergoes a cyclic oscillation with a period of 44.6 years and an amplitude of A₃ = 0.^d0216. For V701 Sco, its O-C curve reveals a periodic change with a period of 41.2 years and amplitude of A₃ = 0.^d0158. The mass ratio of BH Cen is 0.84, but V701 Sco contains twin B1-1.5V type stars with a mass ratio of unit. The continuous period increase of BH Cen is caused by the mass transfer from the less massive component to the more massive one at a rate of dM₂/dt = 3.5 × 10⁻⁶ days/year.The cyclic period changes of both systems can be plausibly explained as the results of light-travel time effects suggesting that they are triple systems. The astrophysical parameters of the unseen tertiary components in the two systems have been determined. We think that the invisible tertiary components in both binaries played an important role in the formations and evolutions of the overcontact configurations by bringing angular momentum out from the central systems. For BH Cen, this process created the initial short period and will support its evolution into an overcontact configuration via a Case A mass transfer within the life time of the extremely young cluster IC 2944. For V701 Sco, two identical zero-age main-sequence components in an overcontact configuration suggest that it may have been formed by fission, possibly by the fission of the third body. The fact that no long-term continuous period variations were found for V701 Sco may suggest that an overcontact binary with the mass ratio of unity can be in an equilibrium revealing that the original configuration of the binary was overcontact as is its present state. It has been reported that faint stars in the two extremely young clusters are relatively scare. From the present study, it is shown that faint stars in young clusters are usually formed as companions of OB stars (including binaries). It is very difficult to detect them because of their low luminosity when compared with the more luminous OB stars.     

Star forming activities in W75N and DR21 – Are the two regions in collision?

We have observed the massive star formation region W75N in ¹²CO J = 3 ? 2 with KOSMA. The profile of ¹²CO J = 3 ? 2 indicated that besides the 9 km s^?1 component, there is another component of ?3 km s^?1, which is associated with another star formation region, DR21N, located to the north of DR21. We derived the physical and dynamical parameters of the core and high velocity gas associated with the two components separately. Star forming activities were investigated, including outflows and infall analysis. The two regions overlap in space and are not connected in velocity. We found that the cloud–cloud collision scenario may not apply for the DR21/W75N case.     

Pressure dependent trace gas trapping in amorphous water ice at 77 K: Implications for determining conditions of comet formation

The nature of cometary volatile materials is subject to debate. Theoretical models of cometary nuclei and laboratory studies suggest that these objects could be made of amorphous water ice in addition to other volatile molecules and refractory grains. This water ice structure has the ability to encapsulate the gases of surrounding environment, reflecting the physical and chemical conditions during their deposition. Therefore, the knowledge of the chemical composition of volatile molecules trapped in amorphous water ice provides a tool for probing the formation environment of cometary ice grains. Experimental studies of gas trapping efficiency in amorphous water ice have been previously conducted mostly under kinetic conditions, where dynamic pumping and temperature gradients prevented rigorous calibrations. In this work, we investigated the trapping efficiencies of Ar, CO, CH₄, Kr and N₂ by depositing water vapor as ice in the presence of trace gases in a volume submerged in liquid nitrogen at 77 K. The gas trapping efficiencies were determined simply by monitoring the pressure difference of the trace gases before and after the deposition of a known amount of water molecules as amorphous ice.Our results show that the trapped gas to water molecule ratio in amorphous ice is controlled primarily by the partial pressure of the gas during water ice deposition, and is independent of the ice deposition rate as well as the gas to water ratio in the vapor phase. The trapping efficiencies of gases decrease in the order of Kr > CH₄ > CO > Ar > N₂ in accordance with previous studies. Assuming that the water ice structure of comets is at least partially amorphous water ice at the time of their formation, these results suggest that the total pressure and composition of the surrounding environment of amorphous ice formation are significant controlling factors of trace gas concentrations in cometary ice. This further indicates that the evolution of the solar nebula and timing of cometary ice condensation can also be important parameters in linking the volatile contents of comets and their formation process.     

Thermal and transport properties of mantle rock at high pressure: Applications to super-Earths

In the present study, the temperature- and pressure-dependent transport and thermal properties, i.e., viscosity, phonon thermal conductivity, thermal expansivity and heat capacities, as well as electronic and radiative thermal conductivities, have been derived for the mantles of super-Earths. These properties are necessary to understand the interior dynamics and the thermal evolution of those planets. We assume that the mantles consist of MgSiO₃ perovskite (pv), but we discuss the effects of the post-perovskite transition, and we elaborate on an addition of periclase MgO and incorporated Fe. However, MgO is found to only significantly influence the phonon thermal conductivity – the viscosities, heat capacities and thermal expansivities of pv and MgO remain comparable. We use the Keane theory of solids, which takes into account the behavior of solid matter at the infinite pressure limit, adopt the Keane equations of state, and adjust for pv and MgO by comparison with experimental high-pressure and high-temperature data. We find the theory of the infinite pressure limit of Keane to be in excellent agreement with recent ab initio studies and experiments. To calculate the melting curve, we further use the Lindemann–Stacey scaling law and fit it to available experimental data. The best data fitting melting temperature for pv reaches 5700 K at 135 GPa and increases to 20,000 K at 1.1 TPa, corresponding to the core-mantle boundary of a 10 Earth mass super-Earth (10M_Earth). We find the pv adiabatic temperature (with a potential temperature of 1700 K) to reach 2570 K at 135 GPa and 5000 K at 1.1 TPa. To calculate the pressure-and temperature-dependent viscosity, we use the semi-empirical homologous temperature scaling to relate enthalpy change, and hence viscosity, to the melting temperature. We find that the resulting activation volume of pv decreases from 2.8 cm³/mol at 25 GPa to 1.4 cm³/mol at 1.1 TPa-resulting in a viscosity increase by ～15 orders of magnitude through the adiabatic mantle of a 10M_Earth planet. Furthermore, the thermal expansivity (of pv and MgO) decreases by a factor of eight, and the total thermal conductivity (phonon, radiative and electronic) of an Earth-like pv/MgO composite increases by a factor of seven through an adiabatic mantle of a 10M_Earth super-Earth. At higher temperatures, i.e., for super-adiabatic temperature profiles, the electronic and radiative thermal conductivities strongly increase and dominate the conductive heat transport. All findings indicate an increase of heat transfer solely by conduction in the lower mantles of super-Earths. Thus our results disagree with Earth-biased full-mantle convection assumptions made by previous models for super-Earths, and additionally raise questions about the differentiation of massive rocky exoplanets and their ability to generate magnetic fields or sustain plate tectonics.     

Radiative heating of carbonaceous near-Earth objects as a cause of thermal metamorphism for CK chondrites

Metamorphic CK carbonaceous chondrites display matrix textures that are best explained by a transient thermal event with temperatures in the 550–950 K range and durations in the order of days to years, longer than what is commonly admitted for shock events but shorter than what is required for nuclide decay. We propose that radiative heating of small carbonaceous meteoroids with perihelia close to the Sun could account for the petrological features observed in CK chondrites. Numerical thermal modeling, using favorable known NEOs orbital parameters (perihelion distances between 0.07 and 0.15 AU) and physical properties of CV and CK chondrites (albedo in the range 0.01–0.1, 25% porosity, thermal diffusivity of 0.5–1.5 W m^?1 K^?1), shows that radiative heating can heat carbonaceous meteoroids in the meter size range to core temperatures up to 1050 K, consistent with the metamorphic temperatures estimated for CK chondrites. Sizes of known CV and CK chondrites indicate that all these objects were small meteoroids (radii from a few cm to 2.5 m) prior to their atmospheric entry. Simulations of dynamic orbits for NEO objects suggest that there are numerous such bodies with suitable orbits and properties, even if they are only a small percentage of all NEOs. Radiative heating would be a secondary process (superimposed on parent-body processes) affecting meteoroids formed by the disruption of an initially homogeneous CV3-type parent body. Different petrologic types can be accounted for depending on the sizes and heliocentric distances of the objects in such a swarm.     

A whole-moon thermal history model of Europa: Impact of hydrothermal circulation and salt transport

A whole-moon numerical model of Europa is developed to simulate its thermal history. The thermal evolution covers three phases: (i) an initial, roughly 0.5 Gyr-long period of radiogenic heating and differentiation, (ii) a long period from 0.5 Gyr to 4 Gyr with continuing radiogenic heating but no tidal dissipative heating (TDH), and (iii) a final period covering the last 0.5 Gyr until the present, during which TDH is active. Hydrothermal plumes develop after the initial period of heating and differentiation and transport heat and salt from Europa’s silicate mantle to its ice shell. We find that, even without TDH, vigorous hydrothermal convection in the rocky mantle can sustain flow in an ocean layer throughout Europa’s history. When TDH becomes active, the ice shell melts quickly to a thickness of about 20 km, leaving an ocean 80 km or more deep. Parameterized convection in the ice shell is non-uniform spatially, changes over time, and is tied to the deeper ocean–mantle dynamics. We also find that the dynamics are affected by salt concentrations. An initially non-uniform salt distribution retards plume penetration, but is homogenized over time by turbulent diffusion and time-dependent flow driven by initial thermal gradients. After homogenization, the uniformly distributed salt concentrations are no longer a major factor in controlling plume transport. Salt transport leads to the formation of a heterogeneous brine layer and salt inclusions at the bottom of the ice shell; the presence of salt in the ice shell could strongly influence convection in that layer.     

Erosive flood events on the surface of Mars: Application to Mangala and Athabasca Valles

We address key factors involved in determining water flow conditions in outflow channels on Mars, including the temperature of the sub-surface water being released and the environmental conditions of low temperature, low atmospheric pressure, and low acceleration due to gravity. We suggest how some of the assumptions made in previous work may be improved. Our model considers the thermodynamic effects of simultaneous evaporation and freezing of water, and fluid dynamical processes including changes in flow rheology caused by assimilation of cold rock and ice eroded at the channel bed, and ice crystal growth due to water freezing. We model how far initially turbulent water could flow in a channel before it erodes and entrains enough material to become laminar, and subsequently ceases to erode the bed. An ice raft will begin to form on the flood while transition occurs between turbulent and laminar flow. Estimates are given for water transit times, ～17–19 h, initial water depths, 50–62 m, and average flow speeds, 5–12 m s^?1, in the Mangala and Athabasca Valles. We show that these two outflow channels, and by implication others like them, could plausibly have been formed in single water release events. Resulting mean erosion rates are approximately 0.7 mm s^?1, a factor of three greater than previous estimates based on combinations of estimates of flood duration and required water volumes. This is explained by the consideration of the effects of eroded ice and the physics of thermal erosion in the present study.     

Low-velocity collisions between centimeter-sized snowballs: Porosity dependence of coefficient of restitution for ice aggregates analogues in the Solar System

Understanding the collisional behavior of ice dust aggregates at low velocity is a key to determining the formation process of small icy bodies such as icy planetesimals, comets and icy satellites, and this collisional behavior is also closely related to the energy dissipation mechanism in Saturn’s rings. We performed head-on collision experiments in air by means of free-falling centimeter-sized sintered snowballs with porosities from 44% to 80% at impact velocities from 0.44 m s^?1 to 4.12 m s^?1 at ?10 °C. In cases of porosity larger than 70%, impact sticking was the dominant collision outcome, while bouncing was dominant at lower porosity. Coefficients of restitution of snow in this velocity range were found to depend strongly on the porosity rather than the impact velocity and to decrease with the increase of the porosity. We successfully measured the compaction volume of snowballs after the impact, and it enabled us to estimate the dynamic compressive strength of snow with the assumption of the energy conservation between kinetic energy and work for deformation, which was found to be consistent with the upper limit of static compressive strength. The velocity dependence of coefficients of restitution of snow was analyzed using a Johnson’s model, and a diagram for collision outcomes among equal-sized sintered snowballs was successfully drawn as a function of porosity and impact velocity.     

A Study of Potential Earth-like Planets in the GJ 876 Planetary System

We perform numerical simulations to investigate potential Earth-like planets in the GJ 876 planetary system. We show that the secular resonances ν₁ and ν₂ (resulting respectively from the inner and outer giant planets) can excite the eccentricities of Earth-like planets with orbits 0.21 AU ≤ a < 0.50 AU and cause them to be ejected out of the system in a short time. However, in the dynamical sense, Earth-like planets potentially exist in the region 0.50 AU ≤ a ≤ 1.00 AU, in stable low-eccentricity orbits which may last up to 10⁵ yr.     

Chandra X-ray observations of two unusual BAL quasars

We report sensitive Chandra X-ray non-detections of two unusual, luminous Iron Low-Ionization Broad Absorption Line Quasars (FeLoBALs). The observations do detect a non-BAL, wide-binary companion quasar to one of the FeLoBAL quasars. We combine X-ray-derived column density lower limits (assuming solar metallicity) with column densities measured from ultraviolet spectra and CLOUDY photoionization simulations to explore whether constant-density slabs at broad-line region densities can match the physical parameters of these two BAL outflows, and find that they cannot. In the “overlapping-trough” object SDSS J0300+0048, we measure the column density of the X-ray absorbing gas to be N_H ? 1.8 × 10²⁴ cm^?2. From the presence of Fe ii UV78 absorption but lack of Fe ii UV195/UV196 absorption, we infer the density in that part of the absorbing region to be n_e ? 10⁶ cm^?3. We do find that a slab of gas at that density might be able to explain this object’s absorption. In the Fe iii-dominant object SDSS J2215–0045, the X-ray absorbing column density of N_H ? 3.4 × 10²⁴ cm^?2 is consistent with the Fe iii-derived N_H ? 2 × 10²² cm^?2 provided the ionization parameter is log U > 1.0 for both the n_e = 10¹¹ cm^?3 and n_e = 10¹² cm^?3 scenarios considered (such densities are required to produce Fe iii absorption without Fe iiabsorption). However, the velocity width of the absorption rules out its being concentrated in a single slab at these densities. Instead, this object’s spectrum can be explained by a low density, high ionization and high temperature disk wind that encounters and ablates higher density, lower ionization Fe iii-emitting clumps.     

The O2 nightglow in the martian atmosphere by SPICAM onboard of Mars-Express

We present observations of the O₂(a¹Δ_g) nightglow at 1.27 μm on Mars using the SPICAM IR spectrometer onboard of the Mars Express orbiter. In contrast to the O₂(a¹Δ_g) dayglow that results from the ozone photodissociation, the O₂(a¹Δ_g) nightglow is a product of the recombination of O atoms formed by CO₂ photolysis on the dayside at altitudes higher than 80 km and transported downward above the winter pole by the Hadley circulation. The first detections of the O₂(a¹Δ_g) nightglow in 2010 indicate that it is about two order of magnitude less intense than the dayglow (Bertaux, J.-L., Gondet, B., Bibring, J.-P., Montmessin, F., Lefèvre, F. [2010]. Bull. Am. Astron. Soc. 42, 1040; Clancy et al. [2010]. Bull. Am. Astron. Soc. 42, 1041). SPICAM IR sounds the martian atmosphere in the near-IR range (1–1.7 μm) with the spectral resolution of 3.5 cm^?1 in nadir, limb and solar occultation modes. In 2010 the vertical profiles of the O₂(a¹Δ_g) nightside emission have been obtained near the South Pole at latitudes of 82–83°S for two sequences of observations: L_s = 111–120° and L_s = 152–165°. The altitude of the emission maximum varied from 45 km on L_s = 111–120° to 38–49 km on L_s = 152–165°. Averaged vertically integrated intensity of the emission at these latitudes has shown an increase from 0.22 to 0.35 MR. Those values of total vertical emission rate are consistent with the OMEGA observations on Mars-Express in 2010. The estimated density of oxygen atoms at altitudes from 50 to 65 km varies from 1.5 × 10¹¹ to 2.5 × 10¹¹ cm^?3. Comparison with the LMD general circulation model with photochemistry (Lefèvre, F., Lebonnois, S., Montmessin, F., Forget, F. [2004]. J. Geophys. Res. 109, E07004; Lefèvre et al. [2008]. Nature 454, 971–975) shows that the model reproduces fairly well the O₂(a¹Δ_g) emission layer observed by SPICAM when the large field of view (>20 km on the limb) of the instrument is taken into account.     

Main Belt Comet P/2008 R1 Garradd: Duration of activity

In this work we investigated changes of the water emission from a model comet of the size and orbital elements of Comet P/2008 R1 (Garradd). We performed simulations for model cometary nuclei of different compositions and two different orientations in space. Our simulations indicate, that the emission of water decreases from one orbital period to another one, but in some cases slowly. When the rotation axis of the nucleus lies in the orbital plane the seasonal maximum of water production during the first two orbital periods can be as high as about 10²⁶ mol s^?1, but decreases by two orders of magnitude during only 50 orbital periods. The highest rate of water production after many orbital periods is expected when the rotation axis is perpendicular to the orbital plane – the seasonal maximum of water production can be about 5 × 10²⁵ mol s^?1 during the first two orbital periods after activation of the comet and no more than 0.8 × 10²⁵ mol s^?1 500 orbital periods later. The upper estimate for the production of water derived from observations of P/2008 R1 (Garradd) by Jewitt et al. (Jewitt, D., Yang, B., Haghighipour, N. [2009]. Astron. J. 137, 4313–4321) is 5 × 10²⁵.