Oral Histories in Meteoritics and Planetary Science – XX: Dale Cruikshank

Abstract– In this interview, Dale Cruikshank ( Fig. 1 ) explains how as an undergraduate at Iowa State University he was a summer student at Yerkes Observatory where he assisted Gerard Kuiper in work on his Photographic Lunar Atlas. Upon completing his degree, Dale went to graduate school at the University of Arizona with Kuiper where he worked on the IR spectroscopy of the lunar surface. After an eventful 1968 trip to Moscow via Prague, during which the Soviets invaded Czechoslovakia, Dale assumed a postdoc position with Vasili Moroz at the Sternberg Astronomical Institute and more observational IR astronomy. Upon returning to the United States and after a year at Arizona, Dale assumed a position at the University of Hawai’i that he held for 17 years. During this period Dale worked with others on thermal infrared determinations of the albedos of small bodies beyond the asteroid Main Belt, leading to the recognition that low‐albedo material is prevalent in the outer solar system that made the first report of complex organic solids on a planetary body (Saturn’s satellite Iapetus). After moving to Ames Research Center, where he works currently, he continued this work and became involved in many outer solar system missions. Dale has served the community through his involvement in developing national policies for science‐driven planetary exploration, being chair of the DPS 1990–1991 and secretary/treasurer for 1982–1985. He served as president of Commission 16 (Physics of Planets) of the IAU (2001–2003). He received the Kuiper prize in 2006.

Figure 1 Open in figure viewer PowerPoint Dale P. Cruikshank.     

Oral Histories in Meteoritics and Planetary Science—XXIII: Dieter Stöffler

In this interview, Dieter Stöffler (Fig. 1) describes how his interest in meteorites and impact craters dates from his Ph.D. studies at the University of Tübingen when it was learned that the Ries crater was formed by impact. A paper by Dieter's advisor, Wolf von Engelhardt, also triggered an interest in meteorites. After graduation, Dieter helped to establish a laboratory for high pressure mineralogy and he examined rocks from the Ries crater, which led to the concept of progressive shock metamorphism. The group also worked on newly returned Apollo samples and guided astronauts over the crater. A year at the NASA Ames Research Center taught Dieter about experimental impact research with a light‐gas gun. After a few more years at Tübingen, Dieter obtained a professorship at the University of Münster where he created the Institute of Planetology, got involved in planning space missions including comet sample return, and continued high pressure mineralogy in collaboration with colleagues in Freiburg. Through several decades of research, Dieter and colleagues have documented the effects of shock on all the major rock‐forming minerals and devised widely accepted schemes for the classification of shocked rocks. After the unification of Germany, Dieter became Director of the Natural History Museum in Berlin, during which he made much progress rebuilding the laboratories and the collections. Dieter also helped to create a museum and research center in the Ries crater. He received the Barringer Award of the Meteoritical Society in 1994 and several prestigious awards in Germany.     

Thermal modeling of shock melts in Martian meteorites: Implications for preserving Martian atmospheric signatures and crystallization of high‐pressure minerals from shock melts

The distribution of shock melts in four shergottites, having both vein and pocket geometry, has been defined and the conductive cooling time over the range 2500 °C to 900 °C calculated. Isolated 1 mm² pockets cool in 1.17 s and cooling times increase with pocket area. An isolated vein 1 × 7 mm in Northwest Africa (NWA) 4797 cools to 900 °C in 4.5 s. Interference between thermal haloes of closely spaced shock melts decreases the thermal gradient, extending cooling times by a factor of 1.4 to 100. This is long enough to allow differential diffusion of Ar and Xe from the melt. Small pockets (1 mm²) lose 2.2% Ar and 5.2% Xe during cooling, resulting in a small change in the Ar/Xe ratio of the dissolved gas over that originally trapped. With longer cooling times there is significant fractionation of Xe from Ar and the Ar/Xe ratio increases rapidly. The largest pockets show less variation of Ar/Xe and likely preserve the original trapped gas composition. Considering all of the model calculations, even the smallest isolated pockets have cooling times greater than the duration of the pressure pulse, i.e., >0.01 s. The crystallization products of these shock melts will be unrelated to the peak shock pressure experienced by the meteorite.     

Numerical simulations of granular dynamics: I. Hard-sphere discrete element method and tests

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular “atmosphere” that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.     

A numerical model of cohesion in planetary rings

We present a numerical method that incorporates particle sticking in simulations using the N-body code pkdgrav to study motions in a local rotating frame, such as a patch of a planetary ring. Particles stick to form non-deformable but breakable aggregates that obey the (Eulerian) equations of rigid-body motion. Applications include local simulations of planetary ring dynamics and planet formation, which typically feature hundreds of thousands or more colliding bodies. Bonding and breaking thresholds are tunable parameters that can approximately mimic, for example, van der Waals forces or interlocking of surface frost layers. The bonding and breaking model does not incorporate a rigorous treatment of internal fracture; rather the method serves as motivation for first-order investigation of how semi-rigid bonding affects the evolution of particle assemblies in high-density environments.We apply the method to Saturn’s A ring, for which laboratory experiments suggest that interpenetration of thin, frost-coated surface layers may lead to weak cohesive bonding. These experiments show that frost-coated icy bodies can bond at the low impact speeds characteristic of the rings. Our investigation is further motivated by recent simulations that suggest a very low coefficient of restitution is needed to explain the amplitude of the azimuthal brightness asymmetry in Saturn’s A ring, and the hypothesis that fine structure in Saturn’s B ring may in part be caused by large-scale cohesion.This work presents the full implementation of our model in pkdgrav, as well as results from initial tests with a limited set of parameters explored. We find a combination of parameters that yields aggregate size distribution and maximum radius values in agreement with Voyager data for ring particles in Saturn’s outer A ring. We also find that the bonding and breaking parameters define two strength regimes in which fragmentation is dominated either by collisions or other stresses, such as tides. We conclude our study with a discussion of future applications of and refinements to our model.     

The natural thermoluminescence of Antarctic meteorites and their terrestrial ages and orbits: A 2010 update

Abstract– We have examined the relationship between natural thermoluminescence (TL) and ²⁶Al in 120 Antarctic meteorites in order to explore the orbital history and terrestrial ages of these meteorites. Our results confirm the observations of Hasan et al. (1987) which were based on 23 meteorites. For most meteorites there was a positive correlation between natural TL and ²⁶Al, reflecting their similarity in decay rate under Antarctic conditions and thus in terrestrial age. For a small group with low TL and high ²⁶Al a small perihelion was proposed. Within this group, natural TL decreases with terrestrial age as determined by ³⁶Cl measurements, although the rate of TL decay is faster (half‐life approximately 10 ka) and the ages that can be determined are smaller (<200 ka) than for most meteorites. The faster decay rate and lower natural TL levels are a reflection of recent exposure to higher radiation doses and higher temperatures, since this history would populate less stable TL traps with smaller electron densities. We sort the 120 meteorites by perihelion and terrestrial age. The normal perihelion group range up to approximately 1000 ka and the small perihelion group range up to approximately 200 ka. An intermediate perihelion group tends to have short terrestrial ages (20–60 ka). There is acceptable agreement between most (34 out of 43) of our present terrestrial age estimates and those determined by isotopic means, the exceptions reflecting complex irradiation histories, long burial times in the Antarctic, or other issues.     

Mixing efficiency in decaying stably stratified turbulence

Mixing efficiency in stratified flows is a measure of the proportion of turbulent kinetic energy that goes into increasing the potential energy of the fluid by irreversible mixing. In this research direct numerical simulations (DNS) and rapid distortion theory (RDT) calculations of transient turbulent mixing events are carried out in order to study this aspect of mixing. In particular, DNS and RDT of decaying, homogeneous, stably-stratified turbulence are used to determine the mixing efficiency as a function of the initial turbulence Richardson number R_it0=(N_L0/_u0)²$R{i}_{t0}={(N{L}_0/u_0)}^2$, where N   is the buoyancy frequency and _L0$L_0$ and _u0$u_0$ are initial length and velocity scales of the turbulence. The results show that the mixing efficiency increases with increasing R_it0$R{i}_{t0}$ for small R_it0$R{i}_{t0}$, but for larger R_it0$R{i}_{t0}$ the mixing efficiency becomes approximately constant. These results are compared with data from towed grid experiments. There is qualitative agreement between the DNS results and the available experimental data, but significant quantitative discrepancies. The grid turbulence experiments suggest a maximum mixing efficiency (at large R_it0$R{i}_{t0}$) of about 6%, while the DNS and RDT results give about 30%. We consider two possible reasons for this discrepancy: Prandtl number effects and non-matching initial conditions. We conclude that the main source of the disagreement probably is due to inaccuracy in determining the initial turbulence energy input in the case of the grid turbulence experiments.     

The Cali meteorite fall: A new H/L ordinary chondrite

Abstract— The fall of the Cali meteorite took place on 6 July 2007 at 16 h 32 ± 1 min local time (21 h 32 ± 1 min UTC). A daylight fireball was witnessed by hundreds of people in the Cauca Valley in Colombia from which 10 meteorite samples with a total mass of 478 g were recovered near 3°24.3′N, 76°30.6′W. The fireball trajectory and radiant have been reconstructed with moderate accuracy. From the computed radiant and from considering various plausible velocities, we obtained a range of orbital solutions that suggest that the Cali progenitor meteoroid probably originated in the main asteroid belt. Based on petrography, mineral chemistry, magnetic susceptibility, thermoluminescence, and bulk chemistry, the Cali meteorite is classified as an H/L4 ordinary chondrite breccia.     

Quasar candidates selection in the Virtual Observatory era

We present a method for the photometric selection of candidate quasars in multiband surveys. The method makes use of a priori knowledge derived from a subsample of spectroscopic confirmed quasi-stellar objects (QSOs) to map the parameter space. The disentanglement of QSOs candidates and stars is performed in the colour space through the combined use of two algorithms, the probabilistic principal surfaces and the negative entropy clustering, which are for the first time used in an astronomical context. Both methods have been implemented in the voneural package on the Astrogrid Virtual Observatory platform. Even though they belong to the class of the unsupervised clustering tools, the performances of the method are optimized by using the available sample of confirmed quasars and it is therefore possible to learn from any improvement in the available 'base of knowledge'. The method has been applied and tested on both optical and optical plus near-infrared data extracted from the visible Sloan Digital Sky Survey (SDSS) and infrared United Kingdom Infrared Deep Sky Survey-Large Area Survey public data bases. In all cases, the experiments lead to high values of both efficiency and completeness, comparable if not better than the methods already known in the literature. A catalogue of optical candidate QSOs extracted from the SDSS Data Release 7 Legacy photometric data set has been produced and is publicly available at the URL http://voneural.na.infn.it/qso.html .     

Images of Jupiter from the pioneer 10 and pioneer 11 infrared radiometers: A comparison with visible and 5-μm images

All of the data acquired at Jupiter by the Infrared Radiometers on board Pioneers 10 and 11 are presented in the form of images with geometric control. The images are compared with 5-μm and visible images taken in the same time frame. The association of dark (blue or brown) and light (white or red) areas with warm and cool areas (at 5, 20, and 45 μm) respectively, extends to nearly all features observed on the planet. Where the normal association of light and dark visible markings with the zonal velocity breaks down (e.g., at the latitude of the South Equatorial Belt during the Pioneer encounters), the infrared emission seems to follow the visible cloud structure rather than the zonal velocity structure. Exceptions to the general rule involve 20-μm radiation, which reflects conditions in the altitude range 0.1–0.3 bar. For example, a comparison between Pioneer 10 and 11 images suggests that the South Equatorial Belt became brighter at 20 μm, but remained constant at other wavelengths between the two encounters.