Collisional history of asteroids: Evidence from Vesta and the Hirayama families

Collisional evolution studies of asteroids indicate that the initial asteroid population at the time mean collisional velocities were pumped up to ～5 km/sec was only modestly larger than it is today; i.e., the asteroid belt was already depleted relative to the mean surface density elsewhere in the planetary region. Numerical simulations of the collisional evolution of hypothetical initial asteroid populations have been run, subject to three constraints: they must (a) evolve to the present observed asteroid size distribution, (b) preserve Vesta's basaltic crust, and (c) produce at least the observed number of major Hirayama families. A “runaway growth” initial asteroid population distribution is found to best satisfy these constraints. A new model is presented for calculating the fragmental size distribution for the disruption of large, gravitationally bound bodies in which the material strength is increased by hydrostatic self-compression. This model predicts that large asteroid behave as intrinsically strong bodies, even if they have had a history of being collisionally fractured. This model, when applied to the breakup of the Themis and Eos family parent bodies, gives size distributions in reasonably good agreement with those observed.     

FACTOR ANALYSIS OF ELEMENTAL ABUNDANCES IN CHONDRITIC AND ACHONDRITIC METEORITES

Another explorative study of the use of factor analysis in meteorite geochemistry has been made. Forty-two major and trace elements were sought from analyses of 80 stony meteorites in recent articles. Incomplete data reduced the matrix to 30 elements in 55 stones. Missing data were substituted by mean values in the groups CC, (E+H+L+LL) and ACH (13, 28, 14 individuals, respectively): the effect of these substitutions was tested empirically. R-mode analysis with varimax rotation was carried out on these three sub-sets and on the whole set: interpretation focused on factor loadings and scores. Results on the three sub-sets gave little information of geochemical value, although the largest achondrite factor (lithophile elements) permits discrimination of eucrites, aubrites, diogenites and howardites. Analysis of all 55 meteorites showed the variance to be dominated by 1, a refractory-lithophile(Al, Mg, Ca, Zr, Sc, U, Th, La, Eu, Yb) factor, and 2, a volatile-chalcophile (Zn, Te, Cd, Bi, Tl) factor. Factor (1) scores will discriminate chondrites from achondrites: factor (2) scores delineate the compositional trend CC1, CC2, CC3, (E+H+L+LL) except for enstatite chondrites Indarch, Abee which fall with CC1. Further progress would need metal, sulphide and other mineral percentages for each meteorite.     

Asteroids and meteorites: Parent bodies and delivered samples

Meteorites may be pieces of main-belt asteroids, derived by cratering collisions. The physical strength of an asteroid critically affects the quantity of ejecta that can be placed in orbits (probably resonant) that evolve to cross the Earth's. Asteroid strengths very widely due to initial composition and size (e.g., weak carbonaceous material or strong rock), subsequent geophysical evolution (e.g., formation of a strong iron core), and subsequent collisional evolution (e.g., conversion of a strong rocky body into a weak rubble pile). The meteorite yield on Earth further depends on meteorite strength, which affects longevity in space and survival through the atmosphere. We show that meteorites may be derived mainly by cratering rather than by disruptive fragmentation and from large main-belt asteroids rather than from small, Earth-approaching bodies. The model combines a wide variety of evidence from various disciplines to yield results consistent with meteorite statistics. However, no claim is made for uniqueness of this model, and many elements still admit considerable uncertainty.     

Field occurrences of liquefaction-induced features: a primer for engineering geologic analysis of paleoseismic shaking

Discussed in this paper are the factors that control the typical manifestations of liquefaction that are found in continental field settings. The factors are given mainly in terms of the local geologic field situation and the geotechnical properties there. A meaningful interpretation of liquefaction-based data for quantitative analysis of paleoseismic shaking requires understanding of both geologic and geotechnical roles in the mode of ground failure at a specific site.

Recommendations are made for the size of the field area that must be searched for liquefaction effects, in order to develop adequate data for engineering geologic/geotechnical analyses of paleoseismicity. The areal extent must be based on an appreciation that the tectonic situation can cause seismically induced liquefaction effects to form in some locales, but not in others nearby, even for a strong earthquake in the region.

Our guidelines for the conduct of the field search and preliminary analysis of the data relate to three issues for which liquefaction features are especially useful in answering: Has there been strong Holocene/latest Pleistocene shaking in the region? Where was the tectonic source? And what was the strength of shaking? Understanding of the various factors that control the manifestations of liquefaction effects, which we present in this paper, is essential for developing credible answers to these questions.     

Engineering geologic and geotechnical analysis of paleoseismic shaking using liquefaction effects: field examples

The greatest impediments to the widespread acceptance of back-calculated ground motion characteristics from paleoliquefaction studies typically stem from three uncertainties: (1) the significance of changes in the geotechnical properties of post-liquefied sediments (e.g., “aging” and density changes), (2) the selection of appropriate geotechnical soil indices from individual paleoliquefaction sites, and (3) the methodology for integration of back-calculated results of strength of shaking from individual paleoliquefaction sites into a regional assessment of paleoseismic strength of shaking. Presented herein are two case studies that illustrate the methods outlined by Olson et al. [Engineering Geology, this issue] for addressing these uncertainties.

The first case study is for a site near Memphis, Tennessee, wherein cone penetration test data from side-by-side locations, one of liquefaction and the other of no liquefaction, are used to readily discern that the influence of post-liquefaction “aging” and density changes on the measured in situ soil indices is minimal. In the second case study, 12 sites that are at scattered locations in the Wabash Valley and that exhibit paleoliquefaction features are analyzed. The features are first provisionally attributed to the Vincennes Earthquake, which occurred around 6100 years BP, and are used to illustrate our proposed approach for selecting representative soil indices of the liquefied sediments. These indices are used in back-calculating the strength of shaking at the individual sites, the results from which are then incorporated into a regional assessment of the moment magnitude, M, of the Vincennes Earthquake. The regional assessment validated the provisional assumption that the paleoliquefaction features at the scattered sites were induced by the Vincennes Earthquake, in the main, which was determined to have M7.5. The uncertainties and assumptions used in the assessment are discussed in detail.     

On early Solar System chronology: Implications of an heterogeneous spatial distribution of Al and Mn

Early Solar System chronology is usually built with the assumption that the distribution of short-lived radionuclides was homogeneous through the solar accretion disk. At present, there is no unambiguous evidence for a homogeneous distribution of short-lived radionuclides in the solar accretion disk, while some data point to a heterogeneous distribution of short-lived radionuclides. In this paper, we explore a possible chronology based on a heterogeneous distribution of ²⁶Al and ⁵³Mn in the accretion disk. Our basic assumption is that the different abundances of extinct short-lived radionuclides in calcium-aluminium-rich inclusions (CAIs) and chondrules are due to spatial rather than temporal differences. We develop a simple model where CAIs and chondrules form contemporaneously, in different spatial locations, and are characterised by distinct initial ²⁶Al and ⁵³Mn abundances. In this model, all evolved bodies are supposed to be originally chondritic, i.e., to be made of a mixture of CAIs, chondrules, and matrix. This mixture determines the initial content in ²⁶Al and ⁵³Mn of a chondritic parent-body as a function of its CAI and chondrule abundance fraction. This approach enables us to calculate coherent ²⁶Al and ⁵³Mn ages from the agglomeration of the parent-body precursors (CAIs and chondrules) until the isotopic closure of ²⁶Al and ⁵³Mn, thereafter called ²⁶Al-⁵³Mn age. We calculate such ²⁶Al-⁵³Mn ages for a diversity of evolved objects, with the constraint that they should be found for realistic chondritic parent-body precursors, i.e., objects having similar or identical petrograpy to the existing chondrite groups. The so defined age of the d’Orbigny angrite is 4.3 ± 1.1 Myr, for the Asuka-881394 eucrite 2.8 ± 1.0 Myr, for the H4 chondrite Sainte Marguerite ∼3 Myr, and for H4 Forest Vale ∼5 Myr. The calculated ²⁶Al-⁵³Mn ages give timescales for the evolution of the respective parent-bodies/meteorites that can be investigated in the light of further petrographic studies. We anchor the calculated relative chronology to an absolute chronology using absolute Pb-Pb ages and relative Hf-W ages of the objects under scrutiny. The precursors of Sainte Marguerite and Forest Vale agglomerated at the same time (∼4565.8 ± 1.2 Ma ago). The precursors of eucrites (Asuka-881394) agglomerated 4564.8 ± 1.2 Ma ago. The precursors of angrites agglomerated late (4561.5 ± 1.8 Ma ago). Our model provides a fully compatible Al-Mg/Mn-Cr/Pb-Pb chronology, and is shown to be robust to reasonable changes in the input parameters. The calculated initial ²⁶Al/²⁷Al ratios are high enough to have ²⁶Al as a possible heat source for differentiation.     

Evidence from coleopteran assemblages for a short but intense cold interlude during the latter part of the MIS11 Interglacial from Quinton, West Midlands, UK

A sequence of fossiliferous organic sediments spanning the whole of an interglacial has been investigated from Quinton, near Birmingham, United Kingdom This interglacial can be correlated on palynological grounds with the Hoxnian Interglacial of England; equivalent to the Holsteinian Interglacial of central Europe and to Marine Isotope Stage 11. A pollen diagram and coleopteran spectrum were obtained from the same series of samples. By using the mutual climatic range method on the coleopteran assemblages, it has been possible to quantify the thermal climatic conditions at this time. They show clear evidence for a sudden and intense climatic oscillation during the late temperate stage of the interglacial—named the Quinton Cold Interlude. At this time mean temperatures of the warmest month (July) declined suddenly by at least 5 °C and mean temperatures of the coldest months (January and February) by more than 10 °C. Subsequently cool-temperate conditions returned. The pollen diagram through this interval shows little or no response to this climatic event. It is too early yet to generalise about the regional significance of the Quinton Cold Interlude but a similar climatic oscillation also occurs towards the closing stage of the MIS11 Interglacial at Hoxne, Norfolk, UK, suggesting that the cold interlude was at least regional in extent. A short climatic oscillation of similar intensity has also been interpreted from coleopteran evidence towards the end of the MIS13 Interglacial at Waverley Wood Farm Pit, Warwickshire, UK, indicating that such events may characterise the later stages of other interglacials. These severe climatic oscillations towards the end of interglacials may have relevance to our understanding of the closing phase of the present interglacial.     

Development of a GIS-based model for extrapolating mesoscale groundwater recharge estimates using integrated geospatial data sets

A geographic information systems (GIS)-based model was used to derive spatially explicit estimates of recharge using the elevation-dependent chloride mass-balance (EDCMB) approach in a 14-basin area in southwestern Nevada, USA. For each of the 14 basins, a non-linear regression equation was developed relating chloride enrichment to the elevation of spring watersheds. The ratio of chloride enrichment as a function of elevation was used in a GIS framework to transform continuous precipitation data to recharge. Spatial masks that represented two definitions of the lower limits of recharge—one definition based on alluvium/non-alluvium boundaries, the other based on both alluvium/non-alluvium boundaries and an elevation threshold—were applied to each basin. Resultant recharge estimates were then summed for each basin. When compared to summaries of previous recharge estimates for the study area, the EDCMB approach produced results that were within 14 and 3% of two studies, but were significantly greater (31%) than a third. GIS proved to be a very effective tool for combining large spatial data with widely different resolutions into an integrated data set, and also proved to be an efficient mechanism for implementing robust statistical models to estimate recharge.

---

Résumé Un modèle basé sur les systèmes d’informations géographiques (SIG) a été utilisé pour obtenir des estimations explicites de la recharge dans l’espace, en utilisant l’approche du bilan des chlorures en fonction de l’altitude (EDCMB en anglais) dans une zone de 14 bassins, dans le Sud-Ouest du Nevada aux Etats-Unis. Pour chacun des 14 bassins, une régression non-linéaire a été établie, reliant l’enrichissement des chlorures à l’altitude des bassins d’alimentation des sources. Le rapport d’enrichissement des chlorures en fonction de l’altitude a été utilisé dans une approche SIG pour transformer des données pluviométriques continues en terme de recharge. Des masques spatiaux représentant deux définitions des limites de recharge les plus basses ont été appliquées à chaque bassin ; l’une basée sur les limites alluvions/autres formations, l’autre sur ces mêmes limites ainsi que sur un seuil altitudinal. Les recharges estimées résultantes ont ensuite été additionnées pour chaque bassin. Comparée aux précédents bilans de recharge estimés dans la zone d’étude, l’approche EDCMB conduit à des résultats se situant entre 14 et 3% de ceux de deux autres études et largement supérieurs (31%) à ceux d’une troisième étude. Le SIG s’est révélé être un outil très efficace pour combiner dans une base de données intégrée, des données spatiales étendues avec des résolutions largement différentes, ainsi que pour implémenter des modèles statistiques conséquents permettant d’estimer la recharge.

---

Resumen Se utilizó un modelo basado en un Sistema de Información Geográfico (SIG) para derivar espacialmente estimados explícitos de recarga usando el enfoque de balance de masa de cloruro dependiente de elevación (EDCMB) en un área de 14 cuencas en el suroeste de Nevada, Estados Unidos. Para cada una de las 14 cuencas se desarrolló una ecuación de regresión no linear que relaciona el enriquecimiento de cloruro con la elevación de las cuencas con manantiales. La relación de enriquecimiento de cloruro en función de la elevación se usó en un marco SIG para transformar datos continuos de precipitación a recarga. Se aplicaron máscaras espaciales a cada cuenca las cuales representaron dos definiciones de los límites inferiores de recarga, una definición basada en límites de aluvión/sin aluvión, el otro basado en ambos límites aluvión/sin aluvión y una elevación normal de fondo. Los estimados resultantes de recarga fueron sumados para cada cuenca. Al compararse con estimados de recarga previos para el área de estudio se encontró que el enfoque EDCMB produjo resultados que estaban entre el 14 y 3% de dos estudios pero fueron significativamente mayores (31%) que un tercer estudio. El SIG ha demostrado ser una herramienta muy efectiva para combinar datos espaciales grandes con resoluciones muy distintas en un grupo integrado de datos y también ha mostrado ser un mecanismo eficiente para implementar modelos estadísticos robustos para estimar recarga.

     

The “wave turbopause”

The “wave turbopause” is defined as the mesospheric altitude level where the temperature fluctuation field indicates a substantial increase in wave amplitudes in the vertical direction.The turbopause altitude is analyzed on the basis of four years of SABER data (2002–2005, Version 1.06). Substantial seasonal and latitudinal variations are found, with some interannual variability also present. Seasonal changes are annual at high latitudes, semi-annual at low latitudes, and a mixture of both at middle latitudes. Southern hemisphere data are similar as in the North if shifted by half a year. Latitudinal variations show a minimum in the tropics and two relative maxima at middle latitudes.The “wave turbopause” is found near to zero-wind lines or low-wind zones (zonal wind). It is compared to rocket and other measurements, and interesting similarities are obtained. The wave turbopause can also be found in the HAMMONIA GCM. A preliminary analysis shows results similar to those of the SABER measurements.     

Oxygen isotope evidence for submarine hydrothermal alteration of the Del Puerto ophiolite, California

The oxygen isotope compositions and metamorphic mineral assemblages of hydrothermally altered rocks from the Del Puerto ophiolite and overlying volcaniclastic sedimentary rocks at the base of the Great Valley sequence indicate that their alteration occurred in a submarine hydrothermal system. Whole rock δ¹⁸O compositions decrease progressively down section (with increasing metamorphic grade): +22.4‰ (SMOW) to +13.8 for zeolite-bearing volcaniclastic sedimentary rocks overlying the ophiolite; +19.6 to +11.6 for pumpellyite-bearing metavolcanic rocks in the upper part of the ophiolite's volcanic member; +12.3 to +8.1 for epidote-bearing metavolcanic rocks in the lower part of the volcanic member; +8.5 to +5.7 for greenschist facies rocks from the ophiolite's plutonic member; +7.6 to +5.8 for amphibolite facies or unmetamorphosed rocks from the plutonic member.

Modelling of fluid-rock interaction in the Del Puerto ophiolite indicates that the observed pattern of upward enrichment in whole rock δ¹⁸O can be best explained by isotopic exchange with discharging¹⁸O-shifted seawater at fluid/rock mass ratios near 2 and temperatures below 500°C.¹⁸O-depleted plutonic rocks necessarily produced during hydrothermal circulation were later removed as a result of tectonism. Submarine weathering and later burial metamorphism at the base of the Great Valley sequence cannot by itself have produced the zonation of hydrothermal minerals and the corresponding variations in oxygen isotope compositions. The pervasive zeolite and prehnite-pumpellyite facies mineral assemblages found in the Del Puerto ophiolite may reflect its origin near an island arc rather than deep ocean spreading center.