Mass loss due to sputtering and thermal processes in meteoroid ablation

Conventional meteoroid theory assumes that the dominant mode of ablation (which we will refer to as thermal ablation) is by evaporation following intense heating during atmospheric flight. Light production results from excitation of ablated meteoroid atoms following collisions with atmospheric constituents. In this paper, we consider the question of whether sputtering may provide an alternative disintegration process of some importance. For meteoroids in the mass range from 10^-3 to and covering a meteor velocity range from 11 to , we numerically modeled both thermal ablation and sputtering ablation during atmospheric flight. We considered three meteoroid models believed to be representative of asteroidal ( mass density), cometary () and porous cometary () meteoroid structures. Atmospheric profiles which considered the molecular compositions at different heights were use in the sputtering calculations. We find that while in many cases (particularly at low velocities and for relatively large meteoroid masses) sputtering contributes only a small amount of mass loss during atmospheric flight, in some cases sputtering is very important. For example, a porous meteoroid at will lose nearly 51% of its mass by sputtering, while a asteroidal meteoroid at will lose nearly 83% of its mass by sputtering. We argue that sputtering may explain the light production observed at very great heights in some Leonid meteors. We discuss methods to observationally test the predictions of these computations. A search for early gradual tails on meteor light curves prior to the commencement of intense thermal ablation possibly represents the most promising approach. The impact of this work will be most dramatic for very small meteoroids such as those observed with large aperture radars. The heights of ablation and decelerations observed using these systems may provide evidence for the importance of sputtering.     

CCD photometry of 23 minor planets

We present CCD photometric observations of 23 main-belt asteroids, of which 8 have never been observed before; thus, the data of these objects are the first in the literature. The majority showed well-detectable light variations, exceeding 0^m1. We have determined synodic periods for 756 Lilliana (936), 1270 Datura (34), 1400 Tirela (1336), 1503 Kuopio (998), 3682 Welther (359), 7505 Furushu (414) and 11436 1969 QR (123), while uncertain period estimates were possible for 469 Argentina (123), 546 Herodias (104) and 1026 Ingrid (53). The shape of the lightcurves of 3682 Welther changed on a short time-scale and showed dimmings that might be attributed to eclipses in a binary system. For the remaining objects, only lower limits of the periods and amplitudes were concluded.     

Near infrared imaging spectroscopy of Venus with the Anglo-Australian Telescope

We have obtained full-disk spatially resolved spectra of the Venus nightside at near-infrared wavelengths during July 2007 using the Anglo-Australian Telescope and Infrared Imager and Spectrograph 2 (IRIS2). The data have been used to map the intensity and rotational temperature of the O₂(a¹Δ_g) airglow band at . The temperatures agree with those obtained in earlier IRIS2 observations and are significantly higher than expected from the Venus International Reference Atmosphere (VIRA) profile. We also report the detection of the corresponding ν=0-1O₂ airglow band at with a similar spatial distribution to the ν=0-0 band. Observations in the thermal window have been used to image surface topography using two different methods of cloud correction. We have also obtained images that can be used to study cloud motion.     

Modeling crater populations on Venus and Titan

We describe a model for crater populations on planets and satellites with dense atmospheres, like those of Venus and Titan. The model takes into account ablation (or mass shedding), pancaking, and fragmentation. Fragmentation is assumed to occur due to the hydrodynamic instabilities promoted by the impactors’ deceleration in the atmosphere. Fragments that survive to hit the ground make craters or groups thereof. Crater sizes are estimated using standard laws in the gravity regime, modified to take into account impactor disruption. We use Monte Carlo methods to pick parameters from appropriate distributions of impactor mass, zenith angle, and velocity. Good fits to the Venus crater populations (including multiple crater fields) can be found with reasonable values of model parameters. An important aspect of the model is that it reproduces the dearth of small craters on Venus: this is due to a cutoff on crater formation we impose, when the expected crater would be smaller than the (dispersed) object that would make it. Hydrodynamic effects alone (ablation, pancaking, fragmentation) due to the passage of impactors through the atmosphere are insufficient to explain the lack of small craters. In our favored model, the observed number of craters (940) is produced by ∼5500 impactors with masses , yielding an age of (1-σ uncertainty) for the venusian surface. This figure does not take into account any uncertainties in crater scaling and impactor population characteristics, which probably increase the uncertainty to a factor of two in age.We apply the model with the same parameter values to Titan to predict crater populations under differing assumptions of impactor populations that reflect present conditions. We assume that the impactors (comets) are made of 50% porous ice. Predicted crater production rates are ≈190 craters . The smallest craters on Titan are predicted to be in diameter, and ≈5 crater fields are expected. If the impactors are composed of solid ice (density ), crater production rates increase by ≈70% and the smallest crater is predicted to be in diameter. We give cratering rates for denser comets and atmospheres 0.1 and 10 times as thick as Titan's current atmosphere. We also explicitly address leading-trailing hemisphere asymmetries that might be seen if Titan's rotation rate were strictly synchronous over astronomical timescales: if that is the case, the ratio of crater production on the leading hemisphere to that on the trailing hemisphere is ≈4:1.     

Tvashtar awakening detected in April 2006 with OSIRIS at the W.M. Keck Observatory

We detected a volcanic outburst in Io's northern hemisphere on 17 April 2006 with the OSIRIS imaging spectrometer at Keck, and confirmed it was still erupting on 2 June 2006. The eruption, which we name 060417A, was located in Tvashtar Paterae, ∼100 km southeast of the February 2000 eruption. The observed temperature was , over a surface area of , providing a total thermal output of .     

A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.     

Saturn eddy momentum fluxes and convection: First estimates from Cassini images

We apply an automated cloud feature tracking algorithm to estimate eddy momentum fluxes in Saturn's southern hemisphere from Cassini Imaging Science Subsystem near-infrared continuum image sequences. Voyager Saturn manually tracked images had suggested no conversion of eddy to mean flow kinetic energy, but this was based on a small sample of <1000 wind vectors. The automated procedure we use for the Cassini data produces an order of magnitude more usable wind vectors with relatively unbiased sampling. Automated tracking is successful in and around the westward jet latitudes on Saturn but not in the vicinity of most eastward jets, where the linearity and non-discrete nature of cloud features produces ambiguous results. For the regions we are able to track, we find peak eddy fluxes and a clear positive correlation between eddy momentum fluxes and meridional shear of the mean zonal wind, implying that eddies supply momentum to eastward jets and remove momentum from westward jets at a rate . The behavior we observe is similar to that seen on Jupiter, though with smaller eddy-mean kinetic energy conversion rates per unit mass of atmosphere (). We also use the appearance and rapid evolution of small bright features at continuum wavelengths, in combination with evidence from weak methane band images where possible, to diagnose the occurrence of moist convective storms on Saturn. Areal expansion rates imply updraft speeds of over the convective anvil cloud area. As on Jupiter, convection preferentially occurs in cyclonic shear regions on Saturn, but unlike Jupiter, convection is also observed in eastward jet regions. With one possible exception, the large eddy fluxes seen in the cyclonic shear latitudes do not seem to be associated with convective events.     

Interaction between eddies and mean flow in Jupiter's atmosphere: Analysis of Cassini imaging data

Beebe et al. [Beebe, R.F., et al., 1980. Geophys. Res. Lett. 17, 1-4] and Ingersoll et al. [Ingersoll, A.P., et al., 1981. J. Geophys. Res. 86, 8733-8743] used images from Voyagers 1 and 2 to analyze the interaction between zonal winds and eddies in Jupiter's atmosphere. They reported a high positive correlation between Jupiter's eddy momentum flux, , and the variation of zonal velocity with latitude, . This correlation implied a surprisingly high rate of conversion of energy from eddies to zonal flow: , a value more than 10% of Jupiter's thermal flux emission. However, Sromovsky et al. [Sromovsky, L.A., et al., 1982. J. Atmos. Sci. 39, 1413-1432] argued that possible biases in the analysis could have caused an artificially high correlation. In addition, significant differences in the derived eddy flux between datasets put into question the robustness of any one result. We return to this long-standing puzzle using images of Jupiter from the Cassini flyby of December 2000. Our method is similar to previous analyses, but utilizes an automatic feature tracker instead of the human eye. The number of velocity vectors used in this analysis is over 200,000, compared to the 14,000 vectors used by Ingersoll et al. We also find a positive correlation between and and derive a global average power per unit mass, , ranging from . Utilizing Ingersoll et al.'s estimate of the mass per unit area involved in the transport, this would imply a rate of energy conversion of . We discuss the implications of this result and employ several tests to demonstrate its robustness.