Relict olivine grains, chondrule recycling, and implications for the chemical, thermal, and mechanical processing of nebular materials

Chondrules and isolated forsterites in five low-subtype ordinary chondrites [NWA 3127 (LL3.1), Sahara 97210 (LL3.2), Wells (LL3.3), Chainpur (LL3.4), and Sahara 98175 (LL3.5)] were studied using petrographic, EMPA, and SIMS techniques to better constrain the origin of chondrules and the olivine grains within them. Our results imply that igneous crystallization, vapor fractionation, redox effects, and open-system behavior were important processes. All olivine grains, including normal, relict, and isolated forsterite grains, show evidence for igneous fractionation under disequilibrium conditions, with olivine crystallizing during rapid cooling (closer to 2000 °C/h than to 100 °C/h). Vapor fractionation is manifested by anti-correlated abundances between refractory elements (Al, Sc, Y, Ti, Ca, V) and volatile elements (Cr, Mn, P, Rb, Fe) in olivine. Redox effects are evidenced in various ways, and imply that Fe, Co, Ni, and P were partitioned more into metal, and V was partitioned more into olivine, under reducing conditions in the most FeO-poor melts. There is no obvious evidence for systematic variations in olivine composition according to meteorite subtype, but shock melting in Sahara 97210 resulted in the injection of glass-derived melt into olivine, resulting in artificially high abundances of Ba, Sr, Na, Ti, and some other incompatible elements in olivine. Terrestrial weathering in a hot desert environment may have mobilized Ba and Sr in some glasses.Our data suggest that chondrules in ordinary chondrites experienced repeated thermal, chemical, and mechanical processing during a “recycling” process over an extended time period, which involved multiple episodes of melting under fluctuating redox and heating conditions, and multiple episodes of chondrule break-up in some cases. Forsterite grains, including normal grains in forsterite-bearing type I chondrules, the cores of isolated forsterites, and relict forsterite in type II chondrules, all crystallized from similar, refractory melts under reducing conditions; relict Mg-olivine and isolated forsterite grains were thus derived from type I chondrules. Olivine in type II chondrules, including normal grains and ferroan overgrowths on relict Mg-olivine, crystallized from more volatile-rich, oxidized, and relatively unfractionated melts. Relict dusty olivine grains in type I chondrules were derived from type II chondrules during incomplete melting episodes involving reduction and some vaporization, with clear (non-dusty) grains in dusty olivine-bearing chondrules crystallizing from the reduced and partly vaporized melts. Melt compositions parental to normal olivine grains in type I and II chondrules are systematically enriched in refractory elements compared to bulk chondrule compositions, implying that chondrules often experienced open-system exchange with more volatile-rich surroundings after some olivine had crystallized, possibly while the chondrules were still partly molten. Type II chondrules could have been derived from type I chondrules by the addition of relatively volatile-rich material, followed by re-melting and little evaporation under oxidizing conditions. In contrast, type I chondrules could have been derived from type II chondrules by re-melting involving more-or-less evaporation under reducing conditions. Chemical, oxygen isotope, and petrographic data are best accommodated by a model in which there were several (>2-3, sometimes ?4-5) melting episodes for most chondrules in ordinary chondrites.     

High-Resolution Records of the Holocene Paleoenvironmental Variation Reflected by Carbonate and Its Isotopic Compositions in Bosten Lake and Response to Glacial Activities

Abstract: The Early Holocene paleoclimate in Bosten Lake on the northern margin of the Tarim Basin, southern Xinjiang, is reconstructed through an analysis of a 953 cm long core (BSTC2000) taken from Bosten Lake. Multiple proxies of this core, including the mineral components of carbonate, carbonate content, stable isotopic compositions of carbonate, Ca/Sr, TOC and C/N and C/S of organic matter, are used to reconstruct the climatic change since 8500 a B.P. The chronology model is made by nine AMS 14C ages of leaves, seeds and organic matter contained in two parallel cores. The climate was cold and wet during 8500 to 8100 a B.P. Temperature increased from 8100 to 6400 a B.P., the climate was warm and humid, and the lake expanded. The lake level was highest during this stage. Then from 6400 to 5100 a B.P., the climate became cold and the lake level decreased slightly. During the late mid-Holocene, the climate was hot and dry from 5100 to 3100 a B.P., but there was a short cold period during 4400 to 3800 a B.P. At this temporal interval, a mass of ice and snow melting water supplied the lake at the early time and made the lake level rise. The second highest lake level stage occurred during 5200 to 3800 a B.P. The climate was cool and wet during 3100 to 2200 a B.P., when the lake expanded with decreasing evaporation. The lake had the last short-term high level during 3100 to 2800 a B.P. After this short high lake level period, the lake shrank because of the long-term lower temperature and reduced water supply. From 2200 to 1200 a B.P., the climate was hot and dry, and the lake shrank greatly. Although the temperature decreased somewhat from 1200 a B.P. to the present, the climate was warm and dry. The lake level began to rise a little again, but it did not reach the river bed altitude of the Konqi River, an outflow river of the Bosten Lake.     

Effects of grain morphology on critical state: a computational analysis

We introduce a new DEM scheme (LS-DEM) that takes advantage of level sets to enable the inclusion of real grain shapes into a classical discrete element method. Then, LS-DEM is validated and calibrated with respect to real experimental results. Finally, we exploit part of LS-DEM potentiality by using it to study the dependency of critical state (CS) parameters such as critical state line (CSL) slope \(\lambda \), CSL intercept \(\varGamma \), and CS friction angle \(\varPhi _{\mathrm{CS}}\) on the grain’s morphology, i.e., sphericity, roundness, and regularity. This study is carried out in three steps. First, LS-DEM is used to capture and simulate the shape of five different two-dimensional cross sections of real grains, which have been previously classified according to the aforementioned morphological features. Second, the same LS-DEM simulations are carried out for idealized/simplified grains, which are morphologically equivalent to their real counterparts. Third, the results of real and idealized grains are compared, so the effect of “imperfections” on real particles is isolated. Finally, trends for the CS parameters (CSP) dependency on sphericity, roundness, and regularity are obtained as well as analyzed. The main observations and remarks connecting particle’s morphology, particle’s idealization, and CSP are summarized in a table that is attempted to help in keeping a general picture of the analysis, results, and corresponding implications.     

Late Quaternary glaciation and equilibrium-line altitudes of the Mayan Ice Cap, Guatemala, Central America

The Sierra los Cuchumatanes (3837 m), Guatemala, supported a plateau ice cap and valley glaciers around Montaña San Juan (3784 m) that totaled ∼ 43 km² in area during the last local glacial maximum. Former ice limits are defined by sharp-crested lateral and terminal moraines that extend to elevations of ∼ 3450 m along the ice cap margin, and to ca. 3000-3300 m for the valley glaciers. Equilibrium-line altitudes (ELAs) estimated using the area-altitude balance ratio method for the maximum late Quaternary glaciation reached as low as 3470 m for the valley glaciers and 3670 m for the Mayan Ice Cap. Relative to the modern altitude of the 0°C isotherm of ∼ 4840 m, we determined ELA depressions of 1110-1436 m. If interpreted in terms of a depression of the freezing level during maximal glaciation along the modern lapse rate of − 5.3°C km^− 1, this ΔELA indicates tropical highland cooling of ∼ 5.9 to 7.6 ± 1.2°C. Our data support greater glacial highland cooling than at sea level, implying a high tropical sensitivity to global climate changes. The large magnitude of ELA depression in Guatemala may have been partially forced by enhanced wetness associated with southward excursions of the boreal winter polar air mass.     

Exploring weathering and regolith transport controls on Critical Zone development with models and natural experiments

The architecture of the Critical Zone, including mobile regolith thickness and depth to the weathering front, is first order controlled by advance of a weathering front at depth and transport of sediment at the surface. Differences in conditions imposed by slope aspect in the Gordon Gulch catchment of the Boulder Creek Critical Zone Observatory present a natural experiment to explore these interactions. The weathering front is deeper and saprolite more decayed on north-facing than on south-facing slopes. Simple numerical models of weathering front advance, mobile regolith production, and regolith transport are used to test how weathering and erosion rates interact in the evolution of weathered profiles. As the processes which attempt are being made to mimic are directly tied to climate variables such as mean annual temperature, the role of Quaternary climate variation in governing the evolution of Critical Zone architecture can be explored with greater confidence.     

Equilibrium H/H fractionations in organic molecules: I. Experimental calibration of ab initio calculations

Carbon-bound hydrogen in sedimentary organic matter can undergo exchange over geologic timescales, altering its isotopic composition. Studies investigating the natural abundance distribution of ¹H and ²H in such molecules must account for this exchange, which in turn requires quantitative knowledge regarding the endpoint of exchange, i.e., the equilibrium isotopic fractionation factor (α_eq). To date, relevant data have been lacking for molecules larger than methane. Here we describe an experimental method to measure α_eq for C-bound H positions adjacent to carbonyl group (H_α) in ketones. H at these positions equilibrates on a timescale of days as a result of keto-enol tautomerism, allowing equilibrium ²H/¹H distributions to be indirectly measured. Molecular vibrations for the same ketone molecules are then computed using Density Functional Theory at the B3LYP/6-311G** level and used to calculate α_eq values for H_α. Comparison of experimental and computational results for six different straight and branched ketones yields a temperature-dependent linear calibration curve with slope = 1.081−0.00376T and intercept = 8.404−0.387T, where T is temperature in degrees Celsius. Since the dominant systematic error in the calculation (omission of anharmonicity) is of the same size for ketones and C-bound H in most other linear compounds, we propose that this calibration can be applied to analogous calculations for a wide variety of organic molecules with linear carbon skeletons for temperatures below 100 °C. In a companion paper (Wang et al., 2009) we use this new calibration dataset to calculate the temperature-dependent equilibrium isotopic fractionation factors for a range of linear hydrocarbons, alcohols, ethers, ketones, esters and acids.     

Stromatoporoid growth forms and Devonian reef fabrics in the Upper Devonian Alexandra Reef System,Canada – Insight on the challenges of applying Devonian reef facies models

Existing facies models for Devonian reef systems can be divided into high‐energy and low‐energy types. A number of assumptions have been made in the development of these models and, in some cases, criteria that distinguish important aspects of the models are poorly defined. The Upper Devonian Alexandra Reef System contains a variety of reef fabrics from different depositional environments and is ideal for studying the range of environments in which stromatoporoids thrived and the facies from these different environments. A wide variety of stromatoporoid growth forms including laminar, tabular, anastamosing laminar and tabular, domal, bulbous, dendroid, expanding conical, concave‐up whorled‐laminar, concave‐up massive tabular and platy‐multicolumnar are present in the Alexandra Reef System. The whorled‐laminar and massive tabular concave‐up growth forms are virtually undocumented from other Devonian reefs but were common in the reef front of the Alexandra, where they thrived in a low‐energy environment around and below fair‐weather wave base. In contrast, high‐energy parts of the reef margin were dominated by bioclastic rubble deposits with narrow ribbon‐like discontinuous bodies of laminar stromatoporoid framestone. In the lagoon, laminar stromatoporoids formed steep‐sided sediment‐dominated bioherms in response to sea‐level rise and flooding. Relying mostly on the different reef facies in the Alexandra system, a new classification scheme for Devonian reef fabrics has been developed. Devonian reef fabrics can be classified as being: (i) sediment‐laden metazoan dominated; (ii) metazoan–microbial dominated (boundstone); (iii) metazoan dominated (framestone); or (iv) metazoan–marine cement dominated. Distinction of these fabrics carries important sedimentary and palaeoecological implications for reconstructing the depositional environment. With examples from the Alexandra Formation, it is demonstrated that reef facies accumulated in a range of depositional environments and that the simple observation of massive stromatoporoids with or without microbial deposits does not automatically imply a high‐energy reef margin, as otherwise portrayed in a number of the existing facies models for these systems.     

Robustness of the holistic seismic risk evaluation in urban centers using the USRi

The Urban Seismic Risk index (USRi) published in a previous article (Carreño et al., Nat Hazards 40:137–172, 2007) is a composite indicator that measures risk from an integrated perspective and guides decision-making for identifying the main interdisciplinary factors of vulnerability to be reduced or intervened. The first step of the method is the evaluation of the potential physical damage (hard approach) as a result of the convolution of the seismic hazard with the physical vulnerability of buildings and infrastructure. Subsequently, a set of social context conditions that aggravate the physical effects is also considered (soft approach). According to this procedure, the physical risk index is evaluated for each unit of analysis from existing loss scenarios, whereas the total risk index is obtained by multiplying the former index by an impact factor using an aggravating coefficient, based on variables associated with the socio-economic conditions of each unit of analysis. The USRi has been developed using the underlying holistic and multi-hazard approach of the Urban Risk Index framework proposed for the evaluation of disaster risk in different megacities worldwide. This article presents the sensitivity analysis of the index to different parameters such as input data, weights and transformation functions used for the scaling or normalization of variables. This analysis has been performed using the Monte Carlo simulation to validate the robustness of this composite indicator, understanding as robustness how the cities maintain the ranking as well as predefined risk level ranges, when compared with the deterministic results of risk. Results are shown for different cities of the world.