Core-collapse supernovae and the Younger Dryas/terminal Rancholabrean extinctions

Early predictions that some supernovae release large quantities of prompt high energy photons are now corroborated by optical identification of core-collapse supernovae associated with extragalactic GRBS (beamed γ-ray bursts) and XRFS (beamed or un-beamed X-ray flashes). Given the in-galaxy supernova frequency and GRB and XRF recurrence statistics, significant Earth-incident events during the past several million years very likely occurred and nearby events should have affected the Earth and other planetary atmospheres, including terrestrial surface solar UV, the Earth’s climate, and its ecology. The Younger Dryas Stadial (∼12,900 to 11,550 calendar yr BP) began with sharply cooler temperatures in the Earth’s northern hemisphere, regional drought, paleoecological evidence compatible with increased UV, and abrupt increases in cosmogenic ¹⁴C and ¹⁰Be in ice and marine cores and tree rings. In North America, stratigraphic and faunal sequences indicate that a major pulse of mammalian extinctions (at least 23-31 genera) began very close to 12,830 calendar yr BP and was sudden: deposits one century younger are devoid of diverse extinct fauna remains. A 10 s beamed GRB within 2 kpc of the Earth delivers 100 kJ m⁻² fluence to the Earth’s atmosphere, where it causes spallation and catalytic reactions depleting 35-50% O³, and producing excess NO_x species (which favor cooling, drought, and surface fertility), ¹⁴C, and ¹⁰Be. An un-beamed, 10⁵⁰ erg hard photon impulse at ∼250 pc produces similar terrestrial atmospheric effects. A well-characterized massive star supernova, the unusually close Vela event (d = 250 ± 30 pc; total energy of 1-2 × 10⁵¹ erg; age constrained from remnant nebula shock velocities considerations at 13,000-16,000 yr and from the pulsar characteristic age at ∼11,400 yr) may have initiated the Younger Dryas climate change, and caused the extinction of the terminal Rancholabrean fauna.     

Penetration of nearby supernova dust in the inner solar system

We investigate the method by which nearby supernovae – within a few tens of pc of the solar system – can penetrate the solar system and deposit live radioactivities on earth. The radioactive isotopic signatures that could potentially leave an observable geological imprint are in the form of refractory metals; consequently, it is likely they would arrive in the form of supernova-produced dust grains. Such grains can penetrate into the solar system more easily than the bulk supernova plasma, which gets stalled and deflected near the solar system due to the solar wind plasma pressure. We therefore examine the motion of charged grains as they decouple from the supernova plasma and are influenced by the solar magnetic, radiation, and gravitational fields. We characterize the dust trajectories with analytical approximations which display the roles of grain size, initial velocity, and surface voltage. These results are verified with full numerical simulations for wide ranges of dust properties. We find that supernova dust grains traverse the inner solar system nearly undeflected, if the incoming grain velocity – which we take to be that of the incident supernova remnant – is comparable to the solar wind speeds and much larger than the escape velocity at 1 AU. Consequently, the dust penetration to 1 AU has essentially 100% transmission probability and the dust capture onto the earth should have a geometric cross section. Our results cast in a new light the terrestrial deposition of radioisotopes from nearby supernovae in the geological past. For explosions beyond ～10 pc from earth, dust grains can still deliver supernova ejecta to earth, and thus the amount of supernova material deposited is set by the efficiency of dust condensation and survival in supernovae. Turning the problem around, we use observations of live ⁶⁰Fe in both deep-ocean and lunar samples to infer a conservative lower bound iron condensation efficiency of M_dust,Fe/M_tot,Fe ? 4  × 10^?4 for the supernova which apparently produced these species 2–3 Myr ago.     

SN2015bh in NGC 2770: LBV core collapse in merging of binary system components

We present spectroscopy and multicolor photometry for the optical transient PSN J09093496+3307204 in the galaxy NGC2770, which has afterwards been transferred into the supernova phase and got the name SN2015bh. Medium-resolution spectral observations were carried out between February 2015 to May 2016 using the focal reducer SCORPIO at the 6-m Russian telescope BTA. They were followed by by photometric observations at the BTA and six other telescopes with 0.5–1m apertures. Both at the phase of the SN impostor (2015a) and at the supernova phase (2015b), besides Balmer emissions, the strong Fe II emissions are seen in the spectrum; so, these spectra resemble those of Williams Fe II type classical novae. Taking into account circumstellar, interstellar and galactic absorption, A_V = 1.^m 14 ±0.^m 15), we determined maximum absolute magnitudes of the object to be M_V =–15.^m 0 ±0.^m 3 at the 2015a phase and of M^V =–18.^m 14 ±0.^m 30 at the 2015b phase. The light curve at the 2015b phase is similar to those of SN IIL. The supernova progenitor is a luminous blue variable (LBV) star with the powerful Hα emission. We consider several hypotheses of supernovae explosions following optical transients related with LBV. The hypothesis of core collapse of an evolved massive star interrupting the process of its merging with massive companion in a binary system (a failed luminous red nova) was chosen as the preferable one for this event.     

Kinematics of stellar populations with RAVE data

We study the kinematics of the Galactic thin and thick disk populations using stars from the RAVE survey’s second data release together with distance estimates from Breddels et al. (2010). The velocity distribution exhibits the expected moving groups present in the solar neighborhood. We separate thick and thin disk stars by applying the X (stellar-population) criterion of Schuster et al. (1993), which takes into account both kinematic and metallicity information. For 1906 thin disk and 110 thick disk stars classified in this way, we find a vertical velocity dispersion, mean rotational velocity and mean orbital eccentricity of (σ_W, 〈V_Φ〉, 〈e〉)_thin = (18 ± 0.3 km s⁻¹, 223 ± 0.4 km s⁻¹, 0.07 ± 0.07) and (σ_W, 〈V_Φ〉, 〈e〉)_thick = (35 ± 2 km s⁻¹, 163 ± 3 km s⁻¹, 0.31 ± 0.16), respectively. From the radial Jeans equation, we derive a thick disk scale length in the range 1.5-2.2 kpc, whose greatest uncertainty lies in the adopted form of the underlying potential. The shape of the orbital eccentricity distribution indicates that the thick disk stars in our sample most likely formed in situ with minor gas-rich mergers and/or radial migration being the most likely cause for their orbits. We further obtain mean metal abundances of 〈[M/H]〉_thin = +0.03 ± 0.17, and 〈[M/H]〉_thick = −0.51 ± 0.23, in good agreement with previous estimates. We estimate a radial metallicity gradient in the thin disk of −0.07 dex kpc⁻¹, which is larger than predicted by chemical evolution models where the disk grows inside-out from infalling gas. It is, however, consistent with models where significant migration of stars shapes the chemical signature of the disk, implying that radial migration might play at least part of a role in the thick disk’s formation.     

Methane evolution from UV-irradiated spacecraft materials under simulated martian conditions: Implications for the Mars Science Laboratory (MSL) mission

Fifteen organic and three inorganic compounds were tested for methane (CH₄) evolution under simulated martian conditions of 6.9 mbar; UVC (200-280 nm) flux of 4 W m⁻²; 20 °C; simulated optical depth of 0.1; and a Mars gas composition of CO₂ (95.3%), N₂ (2.7%), Ar (1.7%), O₂ (0.13%), and water vapor (0.03%). All three inorganic compounds (i.e., NaCl, CaCO₃, graphite) failed to evolve methane at the minimum detection level 0.5 ppm, or above. In contrast, all organic compounds evolved methane when exposed to UV irradiation under simulated martian conditions. The polycyclic aromatic hydrocarbon, pyrene, released the most methane per unit of time at 0.175 nmol CH₄ g⁻¹ h⁻¹, and a spectral reflectance target material used for the MER rovers and Phoenix lander released the least methane at 0.00065 nmol CH₄ cm⁻² h⁻¹. Methane was also released from UV-killed bacterial endospores of Bacillus subtilis. Although all organic compounds evolved methane when irradiated with UV photons under martian conditions, the concentrations of residual organics, biogenic signature molecules, and dead microbial cells should be relatively low on the exterior surfaces of the MSL rover, and, thus, not significant sources of methane contamination. In contrast, kapton tape was found to evolve methane at the rate of 0.00165 nmol CH₄ cm⁻² h⁻¹ (16.5 nmol m⁻² h⁻¹) under the UV and martian conditions tested. Although the evolution of methane from kapton tape was found to decline over time, the large amount of kapton tape used on the MSL rover (lower bound estimated at 3 m²) is likely to create a significant source of terrestrial methane contamination during the early part of the mission.     

Mid-infrared reflectance spectra and optical constants of six iron oxide/oxyhydroxide phases

We have determined the real and imaginary indices of refraction (n and k) for six iron oxide/oxyhydroxide phases—magnetite, maghemite, goethite, lepidocrocite, akaganéite, and ferrihydrite. A single crystal of magnetite was used to derive bulk n and k values from 100-2000 cm⁻¹ (5-100 μm). Synthetic nanocrystalline samples of maghemite, goethite, lepidocrocite, akaganéite, and ferrihydrite were pressed into compact pellets used to determine bulk n and k values from 100-1200 cm⁻¹ (8.33-100 μm). All values of n and k (the optical constants) were determined from specular reflectance spectra acquired at 2 cm⁻¹ spectral sampling using classical Lorentz-Lorenz dispersion theory. In this paper, we present the optical constants of all six minerals and the oscillator parameters with which they were modeled. Use of these optical constants could aid in radiative transfer models of terrestrial dust as well as Mars, the Moon, and airless bodies in the Solar System.     

Dust extinction and X-ray emission from the starburst galaxy NGC 1482

We present the results based on multiwavelength imaging observations of the prominent dust lane starburst galaxy NGC 1482 aimed to investigate the extinction properties of dust existing in the extreme environment. (B-V) colour-index map derived for the starburst galaxy NGC 1482 confirms two prominent dust lanes running along its optical major axis and are found to extend up to ∼11 kpc. In addition to the main lanes, several filamentary structures of dust originating from the central starburst are also evident. Though, the dust is surrounded by exotic environment, the average extinction curve derived for this target galaxy is compatible with the Galactic curve, with R_V = 3.05, and imply that the dust grains responsible for the optical extinction in the target galaxy are not really different than the canonical grains in the Milky Way. Our estimate of total dust content of NGC 1482 assuming screening effect of dust is ∼2.7 × 10⁵ M_⊙, and provide lower limit due to the fact that our method is not sensitive to the intermix component of dust. Comparison of the observed dust in the galaxy with that supplied by the SNe to the ISM, imply that this supply is not sufficient to account for the observed dust and hence point towards the origin of dust in this galaxy through a merger like event.Our multiband imaging analysis reveals a qualitative physical correspondence between the morphologies of the dust and Hα emission lines as well as diffuse X-ray emission in this galaxy. Spatially resolved spectral analysis of the hot gas along outflows exhibit a gradient in the temperature. Similar gradient was also noticed in the measured values of metallicity, indicating that the gas in the halo is not yet enriched. High resolution, 2-8 keV Chandra image reveals a pair of point sources in the nuclear region with their luminosities equal to 2.27 × 10³⁹ erg s⁻¹ and 9.34 × 10³⁹ erg s⁻¹, and are in excess of the Eddington-limit of 1.5 M_⊙ accreting source. Spectral analysis of these sources exhibit an absorbed-power law with the hydrogen column density higher than that derived from the optical measurements.     

Do planetary encounters reset surfaces of near Earth asteroids?

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these ‘space weathering’ processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r^∗, at which planetary encounters can reset the surface. We found that t_sw ∼ 10⁶ yr (at 1 AU) and r^∗ ∼ 5R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 10⁵ yr would require that r^∗ > 20R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 10⁷ yr. We found that t_sw ∝ q², where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.     

Spatially-resolved high-resolution spectroscopy of Venus 2. Variations of HDO, OCS, and SO2 at the cloud tops

While CO, HCl, and HF, that were considered in the first part of this work, have distinct absorption lines in high-resolution spectra and were detected four decades ago, the lines of HDO, OCS, and SO₂ are either very weak or blended by the telluric lines and have not been observed previously by ground-based infrared spectroscopy at the Venus cloud tops. The H₂O abundance above the Venus clouds is typically below the detection limit of ground-based IR spectroscopy. However, the large D/H ratio on Venus facilitates observations of HDO. Converted to H₂O with D/H ≈ 200, our observations at 2722 cm⁻¹ in the Venus afternoon show a H₂O mixing ratio of ∼1.2 ppm at latitudes between ±40° increasing to ±60° by a factor of 2. The observations in the early morning reveal the H₂O mixing ratio that is almost constant at 2.9 ppm within latitudes of ±75°. The measured H₂O mixing ratios refer to 74 km. The observed increase in H₂O is explained by the lack of photochemical production of sulfuric acid in the night time. The recent observations at the P-branch of OCS at 4094 cm⁻¹ confirm our detection of OCS. Four distributions of OCS along the disk of Venus at various latitudes and local times have been retrieved. Both regular and irregular components are present in the variations of OCS. The observed OCS mixing ratio at 65 km varies from ∼0.3 to 9 ppb with the mean value of ∼3 ppb. The OCS scale height is retrieved from the observed limb darkening and varies from 1 to 4 km with a mean value of half the atmospheric scale height. SO₂ at the cloud tops has been detected for the first time by means of ground-based infrared spectroscopy. The SO₂ lines look irregular in the observed spectra at 2476 cm⁻¹. The SO₂ abundances are retrieved by fitting by synthetic spectra, and two methods have been applied to determine uncertainties and detection limits in this fitting. The retrieved mean SO₂ mixing ratio of 350 ± 50 ppb at 72 km favors a significant increase in SO₂ above the clouds since the period of 1980-1995 that was observed by the SOIR occultations at Venus Express. Scale heights of OCS and SO₂ may be similar, and the SO₂/OCS ratio is ∼500 and may be rather stable at 65-70 km under varying conditions on Venus.     

Average inhomogeneities in Milky Way SNII and the PAMELA anomaly

A model is presented to estimate the fraction of Supernova Type-II events (SNII) occurring inside vs. outside a spiral arm for a given star formation episode. The probability distribution function (PDF) for this fraction is given for use in models similar to those of Shaviv et al. [13] and [11]. The calculated PDF for the SNII fraction, SNII_in/total, defined as the number of SNII inside a spiral arm divided by the total number of SNII from a star formation event, provides a constraint on the magnitude of supernova remnant (SNR) concentrations used in cosmic ray propagation models attempting to explain the PAMELA anomaly. Despite the concentration of star formation within spiral arms, this model predicts the majority of SNII events actually occur in inter-arm regions and calls into question the SNR concentration assumption of Shaviv et al.     

The Unusual Absorption Line Spectrum of Quasar SDSS J2220+0109

We report the Balmer broad absorption lines (BALs) in the quasar SDSS J2220 + 0109 discovered from the SDSS data, and present a detailed analysis of the peculiar absorption line spectrum, including the He I* multiplet at λλ3189, 3889 arising from the metastable 2³s-state helium and the Balmer Hα and Hβ lines from the excited hydrogen H I of n = 2 level, which are rarely seen in quasar spectra, as well as many absorption lines arising from the excited Fe II* of the levels 7 955 cm⁻¹, 13 474 cm⁻¹ and 13 673 cm⁻¹ in the wavelength range 3100∼3300 Å. Ca II H, K absorption line doublets also clearly appear in the SDSS spectrum. All absorption lines show a similar blueshifted velocity structure of Δv ≈ − 1500 ∼ 0 km·s⁻¹ relative to the quasar's systematic redshift determined from the emission lines. Detailed analysis suggests that the Balmer absorption lines should arise from the partially ionized region with a column density of N_HI ≈ 10²¹ cm⁻² for an electron density of ne ∼ 10⁶ cm⁻³; and that the hydrogen n = 2 level may be populated via collisional excitation with Lyα pumping.     

Eclipse reappearances of Io: Time-resolved spectroscopy (1.9-4.2 μm)

We obtained time-resolved, near-infrared spectra of Io during the 60-90 min following its reappearance from eclipse by Jupiter on five occasions in 2004. The purpose was to search for spectral changes, particularly in the well-known SO₂ frost absorption bands, that would indicate surface-atmosphere exchange of gaseous SO₂ induced by temperature changes during eclipse. These observations were a follow-on to eclipse spectroscopy observations in which Bellucci et al. [Bellucci et al., 2004. Icarus 172, 141-148] reported significant changes in the strengths of two strong SO₂ bands in data acquired with the VIMS instrument aboard the Cassini spacecraft. One of the bands (4.07 μm [ν₁ + ν₃]) observed by Bellucci et al. is visible from ground-based observatories and is included in our data. We detected no changes in Io’s spectrum at any of the five observed events during the approximately 60-90 min during which spectra were obtained following Io’s emergence from Jupiter’s shadow. The areas of the three strongest SO₂ bands in the region 3.5-4.15 μm were measured for each spectrum; the variation of the band areas with time does not exceed that which can be explained by the Io’s few degrees of axial rotation during the intervals of observation, and in no case does the change in band strength approach that seen in the Cassini VIMS data. Our data are of sufficient quality and resolution to show the weak 2.198 μm (4549.6 cm⁻¹) 4ν₁ band of SO₂ frost on Io for what we believe is the first time. At one of the events (June 22, 2004), we began the acquisition of spectra ∼6 min before Io reappeared from Jupiter’s shadow, during which time it was detected through its own thermal emission. No SO₂ bands were superimposed on the purely thermal spectrum on this occasion, suggesting that the upper limit to condensed SO₂ in the vertical column above Io’s surface was ∼4 × 10⁻⁵ g cm⁻².     

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.     

The transition from complex crater to peak-ring basin on the Moon: New observations from the Lunar Orbiter Laser Altimeter (LOLA) instrument

Impact craters on planetary bodies transition with increasing size from simple, to complex, to peak-ring basins and finally to multi-ring basins. Important to understanding the relationship between complex craters with central peaks and multi-ring basins is the analysis of protobasins (exhibiting a rim crest and interior ring plus a central peak) and peak-ring basins (exhibiting a rim crest and an interior ring). New data have permitted improved portrayal and classification of these transitional features on the Moon. We used new 128 pixel/degree gridded topographic data from the Lunar Orbiter Laser Altimeter (LOLA) instrument onboard the Lunar Reconnaissance Orbiter, combined with image mosaics, to conduct a survey of craters >50 km in diameter on the Moon and to update the existing catalogs of lunar peak-ring basins and protobasins. Our updated catalog includes 17 peak-ring basins (rim-crest diameters range from 207 km to 582 km, geometric mean = 343 km) and 3 protobasins (137-170 km, geometric mean = 157 km). Several basins inferred to be multi-ring basins in prior studies (Apollo, Moscoviense, Grimaldi, Freundlich-Sharonov, Coulomb-Sarton, and Korolev) are now classified as peak-ring basins due to their similarities with lunar peak-ring basin morphologies and absence of definitive topographic ring structures greater than two in number. We also include in our catalog 23 craters exhibiting small ring-like clusters of peaks (50-205 km, geometric mean = 81 km); one (Humboldt) exhibits a rim-crest diameter and an interior morphology that may be uniquely transitional to the process of forming peak rings. A power-law fit to ring diameters (D_ring) and rim-crest diameters (D_r) of peak-ring basins on the Moon [D_ring = 0.14 ± 0.10(D_r)^1.21±0.13] reveals a trend that is very similar to a power-law fit to peak-ring basin diameters on Mercury [D_ring = 0.25 ± 0.14(D_rim)^1.13±0.10] [Baker, D.M.H. et al. [2011]. Planet. Space Sci., in press]. Plots of ring/rim-crest ratios versus rim-crest diameters for peak-ring basins and protobasins on the Moon also reveal a continuous, nonlinear trend that is similar to trends observed for Mercury and Venus and suggest that protobasins and peak-ring basins are parts of a continuum of basin morphologies. The surface density of peak-ring basins on the Moon (4.5 × 10⁻⁷ per km²) is a factor of two less than Mercury (9.9 × 10⁻⁷ per km²), which may be a function of their widely different mean impact velocities (19.4 km/s and 42.5 km/s, respectively) and differences in peak-ring basin onset diameters. New calculations of the onset diameter for peak-ring basins on the Moon and the terrestrial planets re-affirm previous analyses that the Moon has the largest onset diameter for peak-ring basins in the inner Solar System. Comparisons of the predictions of models for the formation of peak-ring basins with the characteristics of the new basin catalog for the Moon suggest that formation and modification of an interior melt cavity and nonlinear scaling of impact melt volume with crater diameter provide important controls on the development of peak rings. In particular, a power-law model of growth of an interior melt cavity with increasing crater diameter is consistent with power-law fits to the peak-ring basin data for the Moon and Mercury. We suggest that the relationship between the depth of melting and depth of the transient cavity offers a plausible control on the onset diameter and subsequent development of peak-ring basins and also multi-ring basins, which is consistent with both planetary gravitational acceleration and mean impact velocity being important in determining the onset of basin morphological forms on the terrestrial planets.     

Using the youngest asteroid clusters to constrain the space weathering and gardening rate on S-complex asteroids

We have extended our earlier work on space weathering of the youngest S-complex asteroid families to include results from asteroid clusters with ages <10⁶ years and to newly identified asteroid pairs with ages <5 × 10⁵ years. We have identified three S-complex asteroid clusters amongst the set of clusters with ages in the range 10^5-6 years—(1270) Datura, (21509) Lucascavin and (16598) 1992 YC2. The average color of the objects in these clusters agrees with the color predicted by the space weathering model of Willman et al. (Willman, M., Jedicke, R., Nesvorný, D., Moskovitz, N., Ivezi?, Z., Fevig, R. [2008]. Icarus 195, 663-673). SDSS five-filter photometry of the members of the very young asteroid pairs with ages <10⁵ years was used to determine their taxonomic classification. Their types are consistent with the background population near each object. The average color of the S-complex pairs is PC₁ = 0.49 ± 0.03, over 5σ redder than predicted by Willman et al. (Willman, M., Jedicke, R., Nesvorný, D., Moskovitz, N., Ivezi?, Z., Fevig, R. [2008]. Icarus 195, 663-673). This may indicate that the most likely pair formation mechanism is a gentle separation due to YORP spin-up leaving much of the aged and reddened surface undisturbed. If this is the case then our color measurement allows us to set an upper limit of ∼64% on the amount of surface disturbed in the separation process. Using pre-existing color data and our new results for the youngest S-complex asteroid clusters we have extended our space weather model to explicitly include the effects of regolith gardening and fit separate weathering and gardening characteristic time scales of τ_w = 960 ± 160 Myr and τ_g = 2000 ± 290 Myr respectively. The first principal component color for fresh S-complex material is PC₁ = 0.37 ± 0.01 while the maximum amount of local reddening is ΔPC₁ = 0.33 ± 0.06. Our first-ever determination of the gardening time is in stark contrast to our calculated gardening time of τ_g ∼ 270 Myr based on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. A possible resolution for the discrepancy is through a ‘honeycomb’ mechanism in which the surface regolith structure absorbs small impactors without producing significant ejecta. This mechanism could also account for the paucity of small craters on (433) Eros.     

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.     

The aftermath of the July 2009 impact on Jupiter: Ammonia, temperatures and particulates from Gemini thermal infrared spectroscopy

We obtained longitudinally resolved thermal infrared spectra (8-13 μm and 17-25 μm) of Jupiter’s impact debris at the Gemini South Telescope on July 24, 2009; five days after the July 19th collision. These were used to study the mechanisms responsible for the redistribution of thermal energy and material (ammonia and stratospheric particulates) following the impact. Upwelling of (8.5 ± 4.1) × 10¹⁴ g of tropospheric air was sufficient to deposit (6.7 ± 4.1) × 10¹² g of NH₃ over a 6° longitude range above the impact core. The NH₃ was distributed over the 20-80 mbar region with a peak abundance of 1.0 ± 0.6 ppm at 45 mbar. Only a 10th of this abundance was observed over the western ejecta, and it is unlikely that these observations were sensitive to NH₃ entrained in the ballistic plume itself. The pattern of excess thermal energy was markedly different from that of Shoemaker-Levy 9 (SL9), with a localized tropospheric perturbation of 2.0 ± 1.0 K at 200-300 mbar and a broader stratospheric warming of up to 3.5 ± 2.0 K at 10-30 mbar. We find no evidence of residual warmth at p < 1 mbar five days after the impact. The excess thermal energy places lower limits on the total energy of the impact (1.8-15.7 × 10²⁶ ergs), which limits the impactor diameter to 70-510 m (depending on the bulk density chosen for the material).The models of the Gemini spectra required three distinct aerosol features, indicative of the mineralogy of the dark particulate debris, centred at 9.1, 10.0 and 18.5 μm. The retrieved opacities for each of these features were distributed over a larger area (9-10° longitude) and at higher altitudes (above the 10-mbar level) than the stratospheric NH₃, and they are more spatially inhomogeneous. This implies the particulates were either entrained with the rising hot plume or created upon plume re-entry and are subsequently redistributed by stratospheric winds. The three particulate features were consistent with a mixture of amorphous iron and magnesium-rich silicates and silicas in the debris field. A broad 10-μm signature was coincident with peaks expected from material rich in amorphous olivines (but poor in pyroxenes), and similar to silicate features observed during SL9. A narrow 9.1-μm signature was interpreted as a combination of amorphous and crystalline silica. Finally, a broad 18.5-μm emitter was not adequately reproduced by a mixture of simple olivines and pyroxenes and remains to be identified.     

Minimum estimates of the amount and timing of gases released into the martian atmosphere from volcanic eruptions

Volcanism has been a major process during most of the geologic history of Mars. Based on data collected from terrestrial basaltic eruptions, we assume that the volatile content of martian lavas was typically ∼0.5 wt.% water, ∼0.7 wt.% carbon dioxide, ∼0.14 wt.% sulfur dioxide, and contained several other important volatile constituents. From the geologic record of volcanism on Mars we find that during the late Noachian and through the Amazonian volcanic degassing contributed ∼0.8 bar to the martian atmosphere. Because most of the outgassing consisted of greenhouse gases (i.e., CO₂ and SO₂) warmer surface temperatures resulting from volcanic eruptions may have been possible. Our estimates suggest that ∼1.1 × 10²¹ g (∼8 ± 1 m m⁻²) of juvenile water were released by volcanism; slightly more than half the amount contained in the north polar cap and atmosphere. Estimates for released CO₂ (1.6 × 10²¹ g) suggests that a large reservoir of carbon dioxide is adsorbed in the martian regolith or alternatively ∼300 cm cm⁻² of carbonates may have formed, although these materials would not occur readily in the presence of excess SO₂. Up to ∼120 cm cm⁻² (2.2 × 10²⁰ g) of acid rain (H₂SO₄) may have precipitated onto the martian surface as the result of SO₂ degassing. The hydrogen flux resulting from volcanic outgassing may help explain the martian atmospheric D/H ratio. The amount of outgassed nitrogen (∼1.3 mbar) may also be capable of explaining the martian atmospheric ¹⁵N/¹⁴N ratio. Minor gas constituents (HF, HCl, and H₂S) could have formed hydroxyl salts on the surface resulting in the physical weathering of geologic materials. The amount of hydrogen fluoride emitted (1.82 × 10¹⁸ g) could be capable of dissolving a global layer of quartz sand ∼5 mm thick, possibly explaining why this mineral has not been positively identified in spectral observations. The estimates of volcanic outgassing presented here will be useful in understanding how the martian atmosphere evolved over time.     

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.     

Deep Impact, Stardust-NExT and the behavior of Comet 9P/Tempel 1 from 1997 to 2010

We present observational data for Comet 9P/Tempel 1 taken from 1997 through 2010 in an international collaboration in support of the Deep Impact and Stardust-NExT missions. The data were obtained to characterize the nucleus prior to the Deep Impact 2005 encounter, and to enable us to understand the rotation state in order to make a time of arrival adjustment in February 2010 that would allow us to image at least 25% of the nucleus seen by the Deep Impact spacecraft to better than 80 m/pixel, and to image the crater made during the encounter, if possible. In total, ∼500 whole or partial nights were allocated to this project at 14 observatories worldwide, utilizing 25 telescopes. Seventy percent of these nights yielded useful data. The data were used to determine the linear phase coefficient for the comet in the R-band to be 0.045 ± 0.001 mag deg⁻¹ from 1° to 16°. Cometary activity was observed to begin inbound near r ∼ 4.0 AU and the activity ended near r ∼ 4.6 AU as seen from the heliocentric secular light curves, water-sublimation models and from dust dynamical modeling. The light curve exhibits a significant pre- and post-perihelion brightness and activity asymmetry. There was a secular decrease in activity between the 2000 and 2005 perihelion passages of ∼20%. The post-perihelion light curve cannot be easily explained by a simple decrease in solar insolation or observing geometry. CN emission was detected in the comet at 2.43 AU pre-perihelion, and by r = 2.24 AU emission from C₂ and C₃ were evident. In December 2004 the production rate of CN increased from 1.8 × 10²³ mol s⁻¹ to Q_CN = 2.75 × 10²³ mol s⁻¹ in early January 2005 and 9.3 × 10²⁴ mol s⁻¹ on June 6, 2005 at r = 1.53 AU.