Quantification of Metals and Semimetals in Carbon‐Rich Rocks: A New Sequential Protocol Including Extraction from Humic Substances

We have developed a new sequential extraction technique that does not require complex procedures and is efficient in determining metal and semimetal contents of carbon‐rich rocks. Six geological reference materials (SBC‐1, SCHS‐1, SCo‐1, SDO‐1, SGR‐1b and SLg‐1) and an in‐house black shale (SH‐1) were selected to test the method, which consists of four main digestion steps involving: (a) dilute HCl acid; (b) NaOH; (c) aqua regia; and (d) hydrofluoric acid. Compared with traditional aqua regia + hydrofluoric acid attack, this new protocol recovers more of the moderately volatile elements during early extraction of humic substances. In addition when compared with reference values, those for most elements are in agreement within uncertainty. Furthermore, this new protocol reveals important information on the partitioning of elements; for instance, steps one and two indicate which elements are associated with carbonates and organic phases from fulvic and humic acid extractions, whereas step three provides results for which elements are associated with sulfide minerals and step four indicates which elements remained in the silicate and oxide phases.     

Partial melting and recrystallization of Archeaan komatiites by residual heat from rapidly accumulated flows

Spinifex-textured komatiites at Honeymoon Well, Western Australia, show evidence of partial melting and recrystallization of original igneous textures. Their textures and mineral compositions differ markedly from those typical of komatiites. Spinifex olivine plates are bent and broken, while interstitial space between spinifex and cumulus olivine is occupied by polygonal aggregates of clinopyroxene, orthopyroxene, minor olivine and plagioclase. Similar granular pyroxene-plagioclase aggregates occur as diffuse veins cutting spinifex zones and cumulate zones of the flows and, in places, form the matrix to a breccia containing corroded fragments of spinifex rock. Thermometry based on the two pyroxene assemblages yields temperatures of 1055° to 1141° C, just below the low-pressure komatiite solidus. Mineral compositions are different from those of typical komatiites: clinopyroxenes are Al-poor and Cr-rich, olivines are unusually iron-rich and depleted in Cr and Ca, and the low-Ca pyroxene is bronzite rather than the more typical pigeonite. We interpret these observations as the results of thermal metamorphism, partial remelting and subsequent slow crystallization of originally normal spinifex-textured komatiite flows. The rocks in question occupy a 40–70 m interval sandwiched between two olivine-rich units: an underlying 90 m-thick olivine adcumulate layer, forming part of the cumulate zone of a basal 160 m-thick flow, and an overlying 1 km-thick extrusive body composed mostly of olivine mesocumulate and adcumulate and capped in turn by spinifex-textured flows. Thermal modelling shows that a sinusoidal temperature profile of cool flow tops and hot flow centres would exist within this sequence shortly after eruption. Conductive thermal relaxation of this profile could reheat spinifex zones to the extent of inducing partial melting and textural reconstitution. Such reheating is largely dependent on the time interval between the emplacement of successive flows. Calculations suggest that at Honeymoon Well the emplacement interval must have been of the order of 10 years or less. Textural reconstitution may have contributed to the development of the thick orthocumulate sequences characteristic of komatiites in the Agnew-Wiluna belt. Present address: Geochemex Australia, P.O. Box 281, West Perth, 6005, Western Australia     

Trace elements in sulfide inclusions from Yakutian diamonds

 Sulfide inclusions in diamonds may provide the only pristine samples of mantle sulfides, and they carry important information on the distribution and abundances of chalcophile elements in the deep lithosphere. Trace-element abundances were measured by proton microprobe in >50 sulfide inclusions (SDI) from Yakutian diamonds; about half of these were measured in situ in polished plates of diamonds, providing information on the spatial distribution of compositional variations. Many of the diamonds were identified as peridotitic or eclogitic from the nature of coexisting silicate or oxide inclusions. Known peridotitic diamonds contain SDIs with Ni contents of 22–36%, consistent with equilibration between olivine, monosulfide solid solution (MSS) and sulfide melt, whereas SDIs in eclogitic diamonds contain 0–12% Ni. A group of diamonds without silicate or oxide inclusions has SDIs with 11–18% Ni, and may be derived from pyroxenitic parageneses. Eclogitic SDIs have lower Ni, Cu and Te than peridotitic SDIs; the ranges of the two parageneses overlap for Se, As and Mo. The Mo and Se contents range up to 700 and 300 ppm, respectively; the highest levels are found in peridotitic diamonds. Among the in-situ SDIs, significant Zn and Pb levels are found in those connected by cracks to diamond surfaces, and these elements reflect interaction with kimberlitic melt. Significant levels of Ru (30–1300 ppm) and Rh (10–170 ppm) are found in many peridotitic SDIs; SDIs in one diamond with wustite and olivine inclusions and complex internal structures have high levels of other platinum-group elements (PGEs) as well, and high chondrite-normalized Ir/Pd. Comparison with experimental data on element partitioning between crystals of monosulfide solid solution (MSS) and sulfide melts suggests that most of the inclusions in both parageneses were trapped as MSS, while some high-Cu SDIs with high Pd±Rh may represent fractionated sulfide melts. Spatial variations of SDI composition within single diamonds are consistent with growth histories shown by cathodoluminescence images, in which several stages of growth and resorption have occurred within magmatic environments that evolved during diamond formation. Received: 5 July 1995 / Accepted: 21 February 1996     

The UVic earth system climate model: Model description,climatology, and applications to past,present and future climates

Abstract

A new earth system climate model of intermediate complexity has been developed and its climatology compared to observations. The UVic Earth System Climate Model consists of a three‐dimensional ocean general circulation model coupled to a thermodynamic/dynamic sea‐ice model, an energy‐moisture balance atmospheric model with dynamical feedbacks, and a thermomechanical land‐ice model. In order to keep the model computationally efficient a reduced complexity atmosphere model is used. Atmospheric heat and freshwater transports are parametrized through Fickian diffusion, and precipitation is assumed to occur when the relative humidity is greater than 85%. Moisture transport can also be accomplished through advection if desired. Precipitation over land is assumed to return instantaneously to the ocean via one of 33 observed river drainage basins. Ice and snow albedo feedbacks are included in the coupled model by locally increasing the prescribed latitudinal profile of the planetary albedo. The atmospheric model includes a parametrization of water vapour/planetary longwave feedbacks, although the radiative forcing associated with changes in atmospheric CO₂ is prescribed as a modification of the planetary longwave radiative flux. A specified lapse rate is used to reduce the surface temperature over land where there is topography. The model uses prescribed present‐day winds in its climatology, although a dynamical wind feedback is included which exploits a latitudinally‐varying empirical relationship between atmospheric surface temperature and density. The ocean component of the coupled model is based on the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model 2.2, with a global resolution of 3.6° (zonal) by 1.8° (meridional) and 19 vertical levels, and includes an option for brine‐rejection parametrization. The sea‐ice component incorporates an elastic‐viscous‐plastic rheology to represent sea‐ice dynamics and various options for the representation of sea‐ice thermodynamics and thickness distribution. The systematic comparison of the coupled model with observations reveals good agreement, especially when moisture transport is accomplished through advection.

Global warming simulations conducted using the model to explore the role of moisture advection reveal a climate sensitivity of 3.0°C for a doubling of CO₂, in line with other more comprehensive coupled models. Moisture advection, together with the wind feedback, leads to a transient simulation in which the meridional overturning in the North Atlantic initially weakens, but is eventually re‐established to its initial strength once the radiative forcing is held fixed, as found in many coupled atmosphere General Circulation Models (GCMs). This is in contrast to experiments in which moisture transport is accomplished through diffusion whereby the overturning is reestablished to a strength that is greater than its initial condition.

When applied to the climate of the Last Glacial Maximum (LGM), the model obtains tropical cooling (30°N‐30°S), relative to the present, of about 2.1°C over the ocean and 3.6°C over the land. These are generally cooler than CLIMAP estimates, but not as cool as some other reconstructions. This moderate cooling is consistent with alkenone reconstructions and a low to medium climate sensitivity to perturbations in radiative forcing. An amplification of the cooling occurs in the North Atlantic due to the weakening of North Atlantic Deep Water formation. Concurrent with this weakening is a shallowing of, and a more northward penetration of, Antarctic Bottom Water.

Climate models are usually evaluated by spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. Implicit in this approach is the assumption that the present‐day observations are in equilibrium with the present‐day radiative forcing. The comparison of a long transient integration (starting at 6 KBP), forced by changing radiative forcing (solar, CO₂, orbital), with an equilibrium integration reveals substantial differences. Relative to the climatology from the present‐day equilibrium integration, the global mean surface air and sea surface temperatures (SSTs) are 0.74°C and 0.55°C colder, respectively. Deep ocean temperatures are substantially cooler and southern hemisphere sea‐ice cover is 22% greater, although the North Atlantic conveyor remains remarkably stable in all cases. The differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene. It is also demonstrated that a global warming simulation that starts from an equilibrium present‐day climate (cold start) underestimates the global temperature increase at 2100 by 13% when compared to a transient simulation, under historical solar, CO₂ and orbital forcing, that is also extended out to 2100. This is larger (13% compared to 9.8%) than the difference from an analogous transient experiment which does not include historical changes in solar forcing. These results suggest that those groups that do not account for solar forcing changes over the twentieth century may slightly underestimate (～3% in our model) the projected warming by the year 2100.     

The thermal and physical characteristics of the Gao-Guenie (H5) meteorite

Measurements of the bulk density, grain density, porosity, and magnetic susceptibility of 19 Gao-Guenie H5 chondrite meteorite samples are presented. We find average values of bulk density 〈ρ_bulk〉=3.46±0.07 g/cm³, grain density 〈ρ_grain〉=3.53±0.08 g/cm³, porosity 〈P(%)〉=2.46±1.39, and bulk mass magnetic susceptibility 〈log χ〉=5.23±0.11. Measurements of the specific heat capacity for a 3.01-g Gao-Guenie sample, a 61.37-g Gao-Guenie sample, a 62.35-g Jilin H5 chondrite meteorite sample, and a 51.37-g Sikhote-Alin IIAB Iron meteorite sample are also presented. Temperature interpolation formula are further provided for the specific heat capacity, thermal conductivity, and thermal diffusivity of the 3.01-g Gao-Guenie sample in the temperature range 300<T (K)<800. We briefly review the possible effects of the newly deduced specific heat and thermal conductivity values on the ablation of meteoroids within the Earth's atmosphere, the modeling of asteroid interiors and the orbital evolution of meteoroids through the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect.     

Resolving the Azimuthal Ambiguity in Vector Magnetogram Data with the Divergence-Free Condition: Application to Discrete Data

We investigate how the divergence-free property of magnetic fields can be exploited to resolve the azimuthal ambiguity present in solar vector magnetogram data, by using line-of-sight and horizontal heliographic derivative information as approximated from discrete measurements. Using synthetic data we test several methods that each make different assumptions about how the divergence-free property can be used to resolve the ambiguity. We find that the most robust algorithm involves the minimisation of the absolute value of the divergence summed over the entire field of view. Away from disk centre this method requires the sign and magnitude of the line-of-sight derivatives of all three components of the magnetic field vector.     

An atlas of predicted exotic gravitational lenses

We present an investigation of the relationships between the radio properties of a giant radio galaxy MRC B0319−454 and the surrounding galaxy distribution with the aim of examining the influence of intergalactic gas and gravity associated with the large-scale structure on the evolution in the radio morphology. Our new radio continuum observations of the radio source, with high surface brightness sensitivity, images the asymmetries in the megaparsec-scale radio structure in total intensity and polarization. We compare these with the three-dimensional galaxy distribution derived from galaxy redshift surveys. Galaxy density gradients are observed along and perpendicular to the radio axis: the large-scale structure is consistent with a model wherein the galaxies trace the ambient intergalactic gas and the evolution of the radio structures are ram-pressure limited by this associated gas. Additionally, we have modelled the off-axis evolution of the south-west radio lobe as deflection of a buoyant jet backflow by a transverse gravitational field: the model is plausible if entrainment is small. The case study presented here is a demonstration that giant radio galaxies may be useful probes of the warm-hot intergalactic medium believed to be associated with moderately over dense galaxy distributions.     

Stratiform and stratabound tungsten mineralisation in the Broken Hill Block,N.S.W.

The Early‐Middle Proterozoic Broken Hill Block contains three types of W occurrences, which show close stratigraphic control. All three types occur within a relatively narrow stratigraphic interval (the ‘Mine Sequence’ Suite of Stevens et al., 1980) comprising a highly variable group of metamorphosed silicic and mafic volcanics, clastic sediments, and exhalative and chemical sediments containing base metals. The first type includes occurrences of W and base metals in bedded calc‐silicate rocks. In the second type, W occurs in layered to non‐layered calcsilicate rocks associated with amphibolite; these are intimately associated in a narrow stratigraphic interval containing abundant, small, Broken Hill type deposits. The third type comprises stratabound, W‐bearing pegmatites, which have been remobilised from quartz‐feldspar‐biotite gneiss and bedded quartz‐tourmaline rocks. Tungsten has been mined only from the third type and only in small quantities. The three types of tungsten deposits show a close spatial relationship with stratiform and stratabound Pb‐Zn mineralisation, including the Broken Hill type. The Pb‐Zn and W deposits are inferred to be genetically related.     

Artificial neural networks: A new method for mineral prospectivity mapping

A multilayer feed‐forward neural network, trained with a gradient descent, back‐propagation algorithm, is used to estimate the favourability for gold deposits using a raster GIS database for the Tenterfield 1:100 000 sheet area, New South Wales. The database consists of solid geology, regional faults, airborne magnetic and gamma‐ray survey data (U, Th, K and total count channels), and 63 deposit and occurrence locations. Input to the neural network consists of feature vectors formed by combining the values from co‐registered grid cells in each GIS thematic layer. The network was trained using binary target values to indicate the presence or absence of deposits. Although the neural network was trained as a binary classifier, output values for the trained network are in the range [0.1, 0.9] and are interpreted to indicate the degree of similarity of each input vector to a composite of all the deposit vectors used in training. These values are rescaled to produce a multiclass prospectivity map. To validate and assess the effectiveness of the neural‐network method, mineral‐prospectivity maps are also prepared using the empirical weights of evidence and the conceptual fuzzy‐logic methods. The neural‐network method produces a geologically plausible mineral‐prospectivity map similar, but superior, to the fuzzy logic and weights of evidence maps. The results of this study indicate that the use of neural networks for the integration of large multisource datasets used in regional mineral exploration, and for prediction of mineral prospectivity, offers several advantages over existing methods. These include the ability of neural networks to: (i) respond to critical combinations of parameters rather than increase the estimated prospectivity in response to each individual favourable parameter; (ii) combine datasets without the loss of information inherent in existing methods; and (iii) produce results that are relatively unaffected by redundant data, spurious data and data containing multiple populations. Statistical measures of map quality indicate that the neural‐network method performs as well as, or better than, existing methods while using approximately one‐third less data than the weights of evidence method.