Geochronology of Precambrian rocks in the Laramie Range and implications for the tectonic framework of Precambrian southern Wyoming

From Casper Mountain; at its northern end, to the northwestern margin of the Laramie anorthosite—syenite complex, in its central parts, the Laramie Range is underlain by granite and granitic gneiss that has a minimum age of 2.54 ± 0.04 Ga (Rb/Sr whole-rock isochron) and by metasedimentary rocks, including marble and quartzite, that appear to overlie the granitic gneiss nonconformably (minimum age: 1.7 Ga based on several horn-blende K/Ar dates). Southward from the anorthosite—syenite complex into Colorado, the Range is underlain chiefly by the Sherman Granite (1.41 Ga; Peterman and Hedge, 1968) and scattered patches of gneiss that are not dated, but are tentatively correlated wit similar gneiss in the southern Medicine Bow Mountains and in the Colorado Front Range, where they are dated as ? 1.7 Ga (Peterman and Hedge, 1968).The Laramie anorthosite—syenite complex (minimum age: ? 1.42 Ga or ? 1.51 Ga if a hornblende K/Ar date is accepted) apparently intruded the suture separating the old (? 2.5 Ga) continental edge from younger (? 1.7 Ga) geosynclinal rocks. The suture, which manifests itself as the Mullen Creek—Nash Fork shear zone in the Medicine Bow Mountains, also is the boundary between ensialic and ensimatic geosynclinal deposition that occurred during the interval 1.7–2.5 Ga ago.K/Ar dates on biotite and muscovite from rocks north of the anorthosite—syenite complex grade from 2.5 Ga on Casper Mountain down to 1.38 Ga near the complex. Near its northern tip, the Laramie Range is crossed by a geochronologic front, separating 2.5 Ga old gneiss whose K/Ar dates were not lowered by subsequent metamorphism from 2.5 Ga old gneiss whose mica dates were reset between 1.4 and 1.6 Ga ago.     

A discrete particle representation of hillslope hydrology: hypothesis testing in reproducing a tracer experiment at Gårdsjön,Sweden

Despite the long history of the continuum equation approach in hydrology, it is not a necessary approach to the formulation of a physically based representation of hillslope hydrology. The Multiple Interacting Pathways (MIPs) model is a discrete realization that allows hillslope response and transport to be simultaneously explored in a way that reflects the potential occurrence of preferential flows and lengths of pathways. The MIPs model uses random particle tracking methods to represent the flow of water within the subsurface alongside velocity distributions that acknowledge preferential flows and transition probability matrices, which control flow pathways. An initial realization of this model is presented here in application to a tracer experiment carried out in Gårdsjön, Sweden. The model is used as an exploratory tool, testing several hypotheses in relation to this experiment. Copyright © 2011 John Wiley & Sons, Ltd.     

Carbon dioxide on the satellites of Saturn: Results from the Cassini VIMS investigation and revisions to the VIMS wavelength scale

Several of the icy satellites of Saturn show the spectroscopic signature of the asymmetric stretching mode of C-O in carbon dioxide (CO₂) at or near the nominal solid-phase laboratory wavelength of 4.2675 μm (2343.3 cm⁻¹), discovered with the Visible-Infrared Mapping Spectrometer (VIMS) on the Cassini spacecraft. We report here on an analysis of the variation in wavelength and width of the CO₂ absorption band in the spectra of Phoebe, Iapetus, Hyperion, and Dione. Comparisons are made to laboratory spectra of pure CO₂, CO₂ clathrates, ternary mixtures of CO₂ with other volatiles, implanted and adsorbed CO₂ in non-volatile materials, and ab initio theoretical calculations of CO₂ * nH₂O. At the wavelength resolution of VIMS, the CO₂ on Phoebe is indistinguishable from pure CO₂ ice (each molecule’s nearby neighbors are also CO₂) or type II clathrate of CO₂ in H₂O. In contrast, the CO₂ band on Iapetus, Hyperion, and Dione is shifted to shorter wavelengths (typically ∼4.255 μm (∼2350.2 cm⁻¹)) and broadened. These wavelengths are characteristic of complexes of CO₂ with different near-neighbor molecules that are encountered in other volatile mixtures such as with H₂O and CH₃OH, and non-volatile host materials like silicates, some clays, and zeolites. We suggest that Phoebe’s CO₂ is native to the body as part of the initial inventory of condensates and now exposed on the surface, while CO₂ on the other three satellites results at least in part from particle or UV irradiation of native H₂O plus a source of C, implantation or accretion from external sources, or redistribution of native CO₂ from the interior.The analysis presented here depends on an accurate VIMS wavelength scale. In preparation for this work, the baseline wavelength calibration for the Cassini VIMS was found to be distorted around 4.3 μm, apparently as a consequence of telluric CO₂ gas absorption in the pre-launch calibration. The effect can be reproduced by convolving a sequence of model detector response profiles with a deep atmospheric CO₂ absorption profile, producing distorted detector profile shapes and shifted central positions. In a laboratory blackbody spectrum used for radiance calibration, close examination of the CO₂ absorption profile shows a similar deviation from that expected from a model. These modeled effects appear to be sufficient to explain the distortion in the existing wavelength calibration now in use. A modification to the wavelength calibration for 13 adjacent bands is provided. The affected channels span about 0.2 μm centered on 4.28 μm. The maximum wavelength change is about 10 nm toward longer wavelength. This adjustment has implications for interpretation of some of the spectral features observed in the affected wavelength interval, such as from CO₂, as discussed in this paper.     

Effects of atmospheric dynamics and ocean resolution on bi-stability of the thermohaline circulation examined using the Grid ENabled Integrated Earth system modelling (GENIE) framework

We have used the Grid ENabled Integrated Earth system modelling (GENIE) framework to undertake a systematic search for bi-stability of the ocean thermohaline circulation (THC) for different surface grids and resolutions of 3-D ocean (GOLDSTEIN) under a 3-D dynamical atmosphere model (IGCM). A total of 407,000 years were simulated over a three month period using Grid computing. We find bi-stability of the THC despite significant, quasi-periodic variability in its strength driven by variability in the dynamical atmosphere. The position and width of the hysteresis loop depends on the choice of surface grid (longitude-latitude or equal area), but is less sensitive to changes in ocean resolution. For the same ocean resolution, the region of bi-stability is broader with the IGCM than with a simple energy-moisture balance atmosphere model (EMBM). Feedbacks involving both ocean and atmospheric dynamics are found to promote THC bi-stability. THC switch-off leads to increased import of freshwater at the southern boundary of the Atlantic associated with meridional overturning circulation. This is counteracted by decreased freshwater import associated with gyre and diffusive transports. However, these are localised such that the density gradient between North and South is reduced tending to maintain the THC off state. THC switch-off can also generate net atmospheric freshwater input to the Atlantic that tends to maintain the off state. The ocean feedbacks are present in all resolutions, across most of the bi-stable region, whereas the atmosphere feedback is strongest in the longitude–latitude grid and around the transition where the THC off state is disappearing. Here the net oceanic freshwater import due to the overturning mode weakens, promoting THC switch-on, but the atmosphere counteracts this by increasing net freshwater input. This increases the extent of THC bi-stability in this version of the model. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Pauling bond strength, bond length and electron density distribution

The power law regression equation, <R(M–O)> = 1.46(<ρ(r _c)>/r)^?0.19, relating the average experimental bond lengths, <R(M–O)>, to the average accumulation of the electron density at the bond critical point, <ρ(r _c)>, between bonded pairs of metal and oxygen atoms (r is the row number of the M atom), determined at ambient conditions for oxide crystals, is similar to the regression equation R(M–O) = 1.41(ρ(r _c)/r)^?0.21 determined for three perovskite crystals at pressures as high as 80 GPa. The pair are also comparable with the equation <R(M–O)> = 1.43(<s>/r)^?0.21 determined for oxide crystals at ambient conditions and <R(M–O)> = 1.39(<s>/r)^?0.22 determined for geometry-optimized hydroxyacid molecules that relate the geometry-optimized bond lengths to the average Pauling bond strength, <s>, for the M–O bonded interactions. On the basis of the correspondence between the equations relating <ρ(r _c)> and <s> with bond length, it seems plausible that the Pauling bond strength might serve a rough estimate of the accumulation of the electron density between M–O bonded pairs of atoms. Similar expressions, relating bond length and bond strength hold for fluoride, nitride and sulfide molecules and crystals. The similarity of the expressions for the crystals and molecules is compelling evidence that molecular and crystalline M–O bonded interactions are intrinsically related. The value of <ρ(r _c)> = r[(1.41)/<R(M–O)>]^4.76 determined for the average bond length for a given coordination polyhedron closely matches the Pauling’s electrostatic bond strength reaching each the coordinating anions of the coordinated polyhedron. Despite the relative simplicity of the expression, it appears to be more general in its application in that it holds for the bulk of the M–O bonded pairs of atoms of the periodic table.     

Classification of Streaky Bay,Mangalo, Ethiudna and Crockers Well: Stony meteorites from South Australia

Abstract— Four meteorites from South Australia are described and classified. Streaky Bay (L4b) and Mangalo (L6e) were received by the South Australian Museum in December, 1989. Ethiudna (L4a) and Crockers Well (LL7) were previously thought to be paired and were previously described only as chondrites.     

Experimental petrology of the basaltic shergottite Yamato‐980459: Implications for the thermal structure of the Martian mantle

Abstract— The Martian meteorite Yamato (Y‐) 980459 is an olivine‐phyric shergottite. It has a very primitive character and may be a primary melt of the Martian mantle. We have conducted crystallization experiments on a synthetic Y‐980459 composition at Martian upper mantle conditions in order to test the primary mantle melt hypothesis. Results of these experiments indicate that the cores of the olivine megacrysts in Y‐980459 are in equilibrium with a melt of bulk rock composition, suggesting that these megacrysts are in fact phenocrysts that grew from a magma of the bulk rock composition. Multiple saturation of the melt with olivine and a low‐calcium pyroxene occurs at approximately 12 ± 0.5 kbar and 1540 ± 10°C, suggesting that the meteorite represents a primary melt that separated from its mantle source at a depth of ?100 km. Several lines of evidence suggest that the Y‐980459 source underwent extensive melting prior to and/or during the magmatic event that produced the Y‐980459 parent magma. When factored into convective models of the Martian interior, the high temperature indicated for the upper Martian mantle and possibly high melt fraction for the Y‐980459 magmatic event suggests a significantly higher temperature at the core‐mantle boundary than previously estimated.     

Activity-composition relations in the system CaCO3-MgCO3 predicted from static structure energy calculations and Monte Carlo simulations

Thermodynamic mixing properties and subsolidus phase relations of the rhombohedral carbonate system, (1 − x) · CaCO₃ − x · MgCO₃, were modelled in the temperature range of 623-2023 K with static structure energy calculations based on well-parameterised empirical interatomic potentials. Relaxed static structure energies of a large set of randomly varied structures in a 4 × 4 × 1 supercell of calcite (a = 19.952 Å, c = 17.061 Å) were calculated with the General Utility Lattice Program (GULP). These energies were cluster expanded in a basis set of 12 pair-wise effective interactions. Temperature-dependent enthalpies of mixing were calculated by the Monte Carlo method. Free energies of mixing were obtained by thermodynamic integration of the Monte Carlo results. The calculated phase diagram is in good agreement with experimental phase boundaries.