Brecciation and chemical heterogeneities of CI chondrites

Fragments in the size range from 40 μm to several hundred μm in the CI chondrites Orgueil, Ivuna, Alais, and Tonk show a wide range of chemical compositions with variations in major elements such as iron (10.4-42.4 wt% FeO), silicon (12.7-42.2 wt% SiO₂), and sulfur (1.01-15.8 wt% SO₃), but also important minor elements such as phosphorous (up to 5.2 wt% P₂O₅) or calcium (up to 6.6 wt% CaO). These variations are the result of the varying mineralogical compositions of these fragments. The distribution of phyllosilicates, magnetites, and possibly ferrihydrite, in particular, control the abundances of these elements. High REE contents—up to 150 times the solar abundances—were observed in phosphates, while matrix and sulfates are REE-depleted. The studied 113 fragments were subdivided into eight lithologies with similar mineralogical and thus chemical properties. The most common is the CGA lithology, consisting of a groundmass of Mg-rich, coarse-grained phyllosilicates and varying abundances of inclusions such as magnetite. The second most abundant lithology is the FGA lithology, consisting of a groundmass of fine-grained Fe-rich phyllosilicates. A rare, but important lithology consists of fragments with high contents of phosphates and other minerals. The proposed model for the evolution of these lithologies is based on a closed system alteration, where mineralogical differences in the lithologies reflect heterogeneities in the starting material. Comparison of our results with literature data indicates a general similarity of the four CI chondrites analyzed. Further comparison of bulk analyses suggests that the mass ‘threshold’ for chemical heterogeneities in CI chondrite samples is smaller than ∼1-2 g.     

The Namuskluft and Dreigratberg sections in southern Namibia (Kalahari Craton,Gariep Belt): a geological history of Neoproterozoic rifting and recycling of cratonic crust during the dispersal of Rodinia until the amalgamation of Gondwana

This paper presents combined U/Pb, Th/U and Hf isotope analyses on detrital and magmatic zircon grains together with whole-rock geochemical analyses of two basement and eight sedimentary rock samples from the Namuskluft and the Dreigratberg in southern Namibia (Gariep Belt). The sedimentary sections evolved during the Cryogenian on the SW part of the Kalahari Craton and where therefore deposited in an active rift setting during the break-up of Rodinia. Due to insufficient palaeomagnetic data, the position of the Kalahari Craton within Rodinia is still under discussion. There are possibilities to locate Kalahari along the western side of Australia/Mawsonland (Pisarevski et al. in Proterozoic East Gondwana: supercontinent assembly and break-up, Geological Society, London, 2003; Evans in Ancient Orogens and modern analogues. Geological Society, London, 2009; and others) or together with the Congo-Sao Francisco and Rio de la Plata Cratons (Li et al. in Prec Res 45: 203–2014, 2008; Frimmel et al. in Int J Earth Sci (Geol Rundsch) 100: 323–354, 2011; and others). It is sill unclear which craton rifted away from the Kalahari Craton during the Cryogenian. Although Middle to Upper Cryogenian magmatic activity is known for the SE Kalahari Craton (our working area) (Richtersveld Suite, Rosh Pinah Fm), all the presented samples show no U/Pb zircon ages younger than ca. 1.0 Ga and non-older than 2.06 Ga. The obtained U/Pb ages fit very well to the exposed basement of the Kalahari Craton (1.0–1.4 Ga Namaqua Province, 1.7–2.0 Ga Vioolsdrif Granite Suite and Orange River Group) and allow no correlation with a foreign craton such as the Rio de la Plata or Australia/Mawsonland. Lu–Hf isotopic signatures of detrital zircon point to the recycling of mainly Palaeoproterozoic and to a smaller amount of Archean crust in the source areas. εHf(t) signatures range between ?24 and +14.8, which relate to T_DM model ages between 1.05 and 3.1 Ga. Only few detrital zircon grains derived from magmas generated from Mesoproterozoic crustal material show more juvenile εHf(t) signatures of +14, +8 to +4 with T_DM model ages of 1.05–1.6 Ga. During Neoproterozoic deposition, only old cratonic crust with an inherited continental arc signature was available in the source area clearly demonstrated by Hf isotope composition of detrital zircon and geochemical bulk analysis of sedimentary rocks. The granodiorites of the Palaeoproterozoic basement underlying Namuskluft section are ca. 1.9 Ga old and show εHf(t) signatures of ?3 to ?5.5 with T_DM model ages of 2.4–2.7 Ga. These basement rocks demonstrate the extreme uplift and deep erosion of the underlying Kalahari Craton at its western margin before general subsidence during Cryogenian and Ediacaran time. The sedimentary sequence of the two examined sections (Namuskluft and Dreigratberg) proposes the presence of a basin and an increasing subsidence at the SW part of the Kalahari Craton during the Cryogenian. Therefore, we propose the initial formation of an intra-cratonic sag basin during the Lower Cryogenian that evolved later to a rift basin at the cratonic margin due to increasing crustal tension and rifting together with the opening of the Adamastor Ocean. As the zircons of the sedimentary rocks filling this basin show neither rift-related U/Pb ages nor an exotic craton as a possible source area, the only plausible sedimentary transport direction providing the found U/Pb ages would be from the E or the SE, directly from the heart of the Kalahari Craton. Due to subsidence and ongoing sedimentation from E/SE directions, the rift-related magmatic rocks were simply covered by the input of old intra-cratonic material that explains the absence of Neoproterozoic zircon grains in our samples. The geochemical analyses show the erosion of a continental arc and related sedimentary rocks with an overall felsic provenance. The source area was a deeply eroded and incised magmatic arc that evolved on continental crust, without any evidence for a passive margin. All of this can be explained by the erosion of rocks related to the Namaqua Belt, which represents one of the two major peaks of zircon U–Pb ages in all analysed samples. Therefore, the Namaqua Belt was well exposed during the Cryogenian, available to erosion and apart from the also well-exposed Palaeoproterozoic basement of the Kalahari Craton one potential source area for the sedimentary rocks in the investigated areas.     

The Simulation of the Opposing Fluxes of Latent Heat and CO<Subscript>2</Subscript> over Various Land-Use Types: Coupling a Gas Exchange Model to a Mesoscale Atmospheric Model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO₂ and H₂O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO₂ fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO₂ uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO₂ and H₂O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.     

Minor- and trace-element zoning in plagioclase: implications for magma chamber processes at Parinacota volcano, northern Chile

Textural and compositional zoning in plagioclase phenocrysts in a sample from Parinacota volcano (Chile) was investigated using backscattered electron images and electron microprobe analysis of major and trace elements. Large (2 mm) oscillatory zoned crystals (type I) with resorption surfaces of moderate An discontinuities (Ⲓ% An) and decreasing trace-element contents (Sr, Mg, Ti) towards the rim reflect melt differentiation and turbulent convection in the main magma body. Early recharge with a low-Sr mafic magma is seen in the core. Small-scale Sr variations in the core indicate limited diffusion and thus residence and differentiation times of the magma shorter than a few thousand years. Smaller crystals (type II) with low trace-element/An ratio reflect the influence of an H₂O-rich melt probably from a differentiated boundary layer. Closed-system in-situ crystallisation, mafic magma recharge and the role of a water-rich differentiated boundary layer can be distinguished from the An-trace element relationships. Crystals apparently move relatively freely between different parts and regimes in the magma chamber, evidence for "convective crystal dispersion". High-Sr type II crystals indicate an earlier input of Sr-rich mafic magma. Recharge of two distinct mafic magma types is thus identified (high-Sr and low-Sr), which must have been present - at increasing recharge rates with time - in the plumbing system throughout the volcano's history.     

The Signature of Sea Spray in the Hexos Turbulent Heat Flux Data

The role of sea spray intransferring heat and moisture across the air-sea interface has remained elusive. Some studies have reported that sea spray does not affect the turbulent air-sea heat fluxes for 10-m wind speeds up to at least 25 m s^-1, while others have reported important spray contributions for wind speeds as low as 12 m s^-1. One goal of the HEXOS (Humidity Exchange over the Sea) program was to quantify spray's contribution to the turbulent air-sea heat fluxes, but original analyses of the HEXOS flux data found the spray signal to be too small to be reliably identified amid the scatter in the data. We look at the HEXOS data again in the context of the TOGA-COARE bulk flux algorithm and a sophisticated microphysical spray model. This combination of quality data andstate-of-the-art modelling reveals a distinct spray signature in virtually all HEXOS turbulent heat flux data collected in winds of 15 m s^-1 and higher. Spray effects are most evident in the latent heat flux data, where spray contributes roughly 10% of the total turbulent flux in winds of 10 m s^-1 and between 10 and 40% in winds of 15–18 m s^-1. The spray contribution to the total sensible heat flux is also at least 10% in winds above 15 m s^-1. These results lead to a new, unified parameterization for the turbulent air-sea heat fluxes that should be especially useful in high winds because it acknowledges both the interfacial and spray routes by which the sea exchanges heat and moisture with the atmosphere.     

Internal sedimentary structure of linear dunes modelled with a cellular automaton

Linear dunes are the most common type of dune found on Earth and exist on several extra-terrestrial bodies, but despite this abundance their internal stratigraphy has not been commonly agreed. A cellular automaton is deployed to simulate the development of linear dunes, starting from a flat bed, under bi-modal oblique wind regimes of varying degrees of asymmetry. The internal stratigraphy of the linear dunes is monitored by keeping track of (buried) erosion surfaces, avalanche deposits and vertical accumulation, as well as the age of last subaerial exposure of the sediments. The simulations show the initial pattern-coarsening of a network of small dunes into fewer larger longitudinal ridges via bedform interactions and Y-junction dynamics. Three newly recognized types of bedform interaction are identified that relate to initial Y-junction dynamics: longitudinal crest-splitting, which creates free dune tips that can interact with adjacent dunes, and laterally oscillating interactions that lead to ephemeral Y-junctions (normal or reverse). The results show that these three bedform interactions leave no persistent signatures in the stratigraphic record. However, a further three bedform interactions involving the superposition of one dune onto another – merging, cannibalizing and repulsion (known from transverse dune field dynamics) – do leave specific evidence in the internal stratigraphy of the remaining dune, a buried interaction surface at a specific inclination. The preservation potential of this interaction surface varies between the three types. After the initial pattern-coarsening phase, the linear dunes become larger and more independent and their crest orientation follows the net resultant transport direction. The stratigraphies of mature dunes under wind regimes of differing asymmetry show that under (nearly) symmetrical winds the dune accumulates mainly vertically, with strata dipping parallel to the flanks, while under more asymmetrical wind regimes the internal stratigraphy resembles that of transverse dunes.     

Monazite as a monitor for melt‐rock interaction during cooling and exhumation

Granulite facies cordierite–garnet–biotite gneisses from the southeastern Reynolds Range, central Australia, contain both orthopyroxene‐bearing and orthopyroxene‐free quartzofeldspathic leucosomes. Mineral reaction microstructures at the interface of gneiss and leucosome observed in outcrop and petrographically, reflect melt‐rock interaction during crystallization. Accessory monazite, susceptible to fluid alteration, dissolution and recrystallization at high temperature, is tested for its applicability to constrain the chemical and P–T–time evolution of melt‐rock reactions during crystallization upon cooling. Bulk rock geochemistry and phase equilibria modelling constrain peak pressure and temperature conditions to 6.5–7.5 kbar and ~850°C, and U–Pb geochronology constrains the timing of monazite crystallization to 1.55 Ga, coeval with the Chewings Orogeny. Modelling predicts the presence of up to 15 vol.% melt at peak metamorphic conditions. Upon cooling below 800°C, melt extraction and in situ crystallization of melt decrease the melt volume to less than 7%, at which time it becomes entrapped and melt pockets induce replacement reactions in the adjacent host rock. Replacement reactions of garnet, orthopyroxene and K‐feldspar liberate Y, REE, Eu and U in addition to Mg, Fe, Al, Si and K. We demonstrate that distinguishing between monazite varieties solely on the basis of U–Pb ages cannot solve the chronological order of events in this study, nor does it tie monazite to the evolution of melt or stability of rock‐forming minerals. Rather, we argue that analyses of various internal monazite textures, their composition and overprinting relations allow us to identify the chronology of events following the metamorphic peak. We infer that retrograde reactions involving garnet, orthopyroxene and K‐feldspar can be attributed to melt‐rock interaction subsequent to partial melting, which is reflected in the development of compositionally distinct monazite textural domains. Internal monazite textures and their composition are consistent with dissolution and precipitation reactions induced by a high‐T melt. Monazite rims enriched in Y, HREE, Eu and U indicate an increased availability of these elements, consistent with the breakdown of orthopyroxene, garnet and K‐feldspar observed petrographically. Our study indicates that compositional and textural analysis of monazite in relation to major rock‐forming minerals can be used to infer the post‐peak chemical evolution of partial melts during high‐ to ultrahigh‐temperature metamorphism.     

Mobility and H2O loss from fluid inclusions in natural quartz crystals

Petrographic observations on quartz crystals from the Mole Granite (Australia) and other localities shed new light on the mechanisms of post-entrapment modification of fluid inclusions. These modifications include migration away from pseudosecondary trails, changes in fluid salinity and density, shape distortion and the formation of “sweat-haloes” around strongly deformed inclusions. Increases in fluid salinity, which usually are associated with inclusion migration, indicate water-losses of up to 50%. However, LA-ICP-MS-analysis of unmobilized and mobilized inclusions of the same trail reveals basically unchanged ratios of major – and trace element cations, with the exception of Li, which seems to be incorporated into the crystal lattice during migration. Despite the fact that all these modifications are closely related to deformation processes, they occur not only in mechanically deformed quartz, but also in free-standing crystals. In the latter samples, stress has been generated internally as a result of brazil-twinned growth and compositional zonation. These observations and their interpretation leads to a list of practical criteria that should help in differentiating between reliable and suspect fluid inclusions in other samples. Received: 17 November 1998 / Accepted: 16 April 1999     

Interpretation of tracer tests performed in fractured rock of the Lange Bramke basin, Germany

 Two multitracer tests performed in one of the major cross-fault zones of the Lange Bramke basin (Harz Mountains, Germany) confirm the dominant role of the fault zone in groundwater flow and solute transport. Tracers having different coefficients of molecular diffusion (deuterium, bromide, uranine, and eosine) yielded breakthrough curves that can only be explained by a model that couples the advective–dispersive transport in the fractures with the molecular diffusion exchange in the matrix. For the scale of the tests (maximum distance of 225 m), an approximation was used in which the influence of adjacent fractures is neglected. That model yielded nearly the same rock and transport parameters for each tracer, which means that the single-fracture approximation is acceptable and that matrix diffusion plays an important role. The hydraulic conductivity of the fault zone obtained from the tracer tests is about 1.5×10^–2 m/s, whereas the regional hydraulic conductivity of the fractured rock mass is about 3×10^–7 m/s, as estimated from the tritium age and the matrix porosity of about 2%. These values show that the hydraulic conductivity along the fault is several orders of magnitude larger than that of the remaining fractured part of the aquifer, which confirms the dominant role of the fault zones as collectors of water and conductors of fast flow. Received, April 1997 Revised, January 1998, August 1998 Accepted, August 1998     

A Density Correction for Lagrangian Particle Dispersion Models

Current Lagrangian particle dispersion models, used to simulate the dispersion of passive tracers in the turbulent planetary boundary layer (PBL), assume that the density is constant within the PBL. In deep PBLs, where the density at the boundary-layer top may be lower by more than 20% than at the surface, this assumption leads to errors in the tracer concentrations on the order of 10%. In the presence of a vertical wind shear, this also leads to inaccurate calculations of the horizontal tracer transport. To remove this deficiency, a Langevin equation is presented that contains a density correction term. The effect of the density correction is studied using data from a large-scale tracer experiment. It is found that for this experiment, the main effect of the density correction is an increase in the surface tracer concentrations, whereas the horizontal tracer transport patterns remain largely unaffected.