Pyrolysis-high resolution gas chromatography and pyrolysis gas chromatography-mass spectrometry of kerogens and kerogen precursors

A series of kerogens and kerogen precursors isolated from DSDP samples, oil shales and Recent algal mats have been examined by Curie point pyrolysis-high resolution gas chromatography and gas chromatography-mass spectrometry. This study has shown that the three main types of kerogens (marine, terrestrial and mixtures of both) can be characterized using these techniques. The marine (algal) kerogens yield principally aliphatic products and the terrestrial kerogens yield more aromatic and phenolic products with some n-alkanes and n-alkenes. The yields of n-alkanes and n-alkenes increase and phenols decrease with increasing geologic age, however, pyrolysis-GC cannot be used to characterize the influence of short term diagenesis on the kerogen structure.     

Polymerization of silicate and aluminate tetrahedra in glasses,melts, and aqueous solutions—I. Electronic structure of H6Si2O7, H6AlSiO71−, and H6Al2O72−

As part of a study of the effect of geologically common network modifiers on polymerization in silicate melts, glasses, and silica-rich aqueous solutions, we have studied the energies, electronic structures, and inferred chemical properties of ^IVT-O-^IVT linkages in the tetrahedral dimers H₆,Si₂O₇, H₆AlSiO₇^1?, and H₆Al₂O₇^2? using semi-empirical molecular orbital theory (CNDO/2). Our results indicate that the electron donating character of the bridging oxygen, O(br), linking two tetrahedra increases with increasing T-O(br) bond length but decreases with decreasing T-O(br)-T angles and increasing O-T-O(br) angles. This increase or decrease of the donor character of O(br) coincides with an increase or decrease of the affinity of O(br) for hard acceptors. The calculated electronic structure for the H₆Si₂O₇ molecule is compared with the observed X-ray emission, absorption, and photoelectron spectra of quartz and vitreous silica; the reasonable match between calculated and observed oxygen Kα emission spectra of vitreous silica supports our assertion that non-bonded O(br) electron density energetically at the top of the valence band controls the chemical reactivity of ^IVT-O-^IVT linkages in polymerized tetrahedral environments.     

Woolshed Creek fossil site: a key part of Canberra’s scientific and cultural heritage

The Woolshed Creek fossil site near the Royal Military College, Duntroon, Canberra, contains brachiopods Atrypa duntroonensis (early Homerian, early Silurian, ca 430.5?Ma) within a mudstone of the Canberra Formation. Their discovery in 1844 by the Reverend William B. Clarke (“the Father of Australian Geology”), and subsequent comparison with other fossil collections from around the world, contributed significantly to the nineteenth century debate about the oldest rocks in Australia. The fossil site is now on the ACT Government Heritage List and recent site improvements make it readily accessible via a pathway from the sports grounds of the Royal Military College.     

Use of Noise to Augment Training Data: A Neural Network Method of Mineral–Potential Mapping in Regions of Limited Known Deposit Examples

One of the main factors that affects the performance of MLP neural networks trained using the backpropagation algorithm in mineral-potential mapping isthe paucity of deposit relative to barren training patterns. To overcome this problem, random noise is added to the original training patterns in order to create additional synthetic deposit training data. Experiments on the effect of the number of deposits available for training in the Kalgoorlie Terrane orogenic gold province show that both the classification performance of a trained network and the quality of the resultant prospectivity map increasesignificantly with increased numbers of deposit patterns. Experiments are conducted to determine the optimum amount of noise using both uniform and normally distributed random noise. Through the addition of noise to the original deposit training data, the number of deposit training patterns is increased from approximately 50 to 1000. The percentage of correct classifications significantly improves for the independent test set as well as for deposit patterns in the test set. For example, using ±40% uniform random noise, the test-set classification performance increases from 67.9% and 68.0% to 72.8% and 77.1% (for test-set overall and test-set deposit patterns, respectively). Indices for the quality of the resultant prospectivity map, (i.e. D/A, D × (D/A), where D is the percentage of deposits and A is the percentage of the total area for the highest prospectivity map-class, and area under an ROC curve) also increase from 8.2, 105, 0.79 to 17.9, 226, 0.87, respectively. Increasing the size of the training-stop data set results in a further increase in classification performance to 73.5%, 77.4%, 14.7, 296, 0.87 for test-set overall and test-set deposit patterns, D/A, D × (D/A), and area under the ROC curve, respectively.     

Structure and viscosity of rhyolitic composition melts

The average local structure of a rhyolitic composition glass has been determined at 25°C using X-ray radial distribution analysis (RDA) and quasi-crystalline modelling and is best described as similar to that in a stuffed framework composed principally of six-membered rings of Si and Al tetrahedra (basically a stuffed tridymite-like model). Using this model it is possible to calculate a density (2.41 g/cm³) which compares well with the measured density (2.40 g/cm³); a structural model based on four-membered rings (an albite-like model) results in a substantially higher calculated density (2.60 g/cm³). We suggest that the rhyolite glass structural model is appropriate for rhyolitic melts, based on evidence from the recent literature. New viscosity data for an anhydrous rhyolite composition measured between 1200°C and 1500°C are presented and interpreted in terms of our proposed structural model and previous melt structure models for the major normative components of rhyolite. A mechanism for diffusion and viscous flow in framework silicate melts (including rhyolite composition) is proposed on the basis of recent molecular orbital calculations and molecular dynamics simulations of silicate and fluoride melts.     

Phase relations of biotite and stilpnomelane in the greenschist facies

Phase relations of biotite and stilpnomelane and associated silicate minerals have been studied in rocks of the greenschist facies, chiefly from Otago, New Zealand and western Vermont, but also from Scotland, Minnesota-Michigan iron range, and northwest Washington. That stilpnomelane in the greenschicht facies crystallizes initially with nearly all iron in the ferrous state is indicated by chemical analyses, high p-T experiments, and phase relationships. Alteration of stilpnomelane after metamorphism not only oxidizes iron but leaches potassium; corrections for both effects must be made in using analyses of brown stilpnomelane in studies of phase relations. Two discontinuous reactions which produce biotite at the biotite isograd have been identified:

1. muscovite+stilpnomelane+actinolite→ biotite+chlorite+epidote
2. chlorite+microcline→ biotite+muscovite. Biotite produced by the first of these reactions has a limited range of variation in Fe/Mg. As grade advances within the biotite zone more magnesian and ferruginous biotites become stable in consequence of the two continuous reactions:
3. muscovite+actinolite+chlorite→ biotite (Mg-rich)+epidote
4. muscovite+stilpnomelane→ biotite (Fe-rich)+chlorite.

Stilpnomelane is stable in muscovite-free rocks throughout the biotite zone, and even up to the grade at which hornblende becomes stable. Phengitic muscovite is stable throughout the biotite zone in New Zealand and thus apparently does not contribute to the formation of biotite until a higher grade is reached.     

Application of stress mapping in cross‐section to understanding ore geometry,predicting ore zones and development of drilling strategies

Stress mapping is a numerical modelling technique used to determine the distribution and relative magnitude of stress during deformation in a mineralised terrane. It is based on the general principle that fluid flow in the Earth's crust is primarily related to pressure gradients. It is best applied to epigenetic hydrothermal mineral deposits, where fluid flow and fluid flux are enhanced in dilational sections of structures and in sites of enhanced rock permeability due to high fracture density. These are defined by sites of low minimum principal stress (σ₃). Most stress mapping is carried out in two dimensions in plan view using geological maps. This is suitable for terranes with steeply dipping lithostratigraphy and structures in which the distribution of mineral deposits is largely controlled by fault structures portrayed on the maps. However, for terranes with gently dipping sequences and structures, and for situations where deposits are sited in and near the hinges of complex fold structures, stress mapping in cross‐section is preferable. The effectiveness of stress mapping is maximised if mineralisation was late in the evolutionary history of the host terrane, and hence the structural geometry of the terrane and contained deposits were essentially that expressed today. The orientation of syn‐mineralisation far‐field stresses must also be inferred. Two examples of orogenic gold deposits, which meet the above criteria, are used to illustrate the potential of stress mapping in cross‐section. Sunrise Dam, located in the Archaean Yilgarn Craton, is a lode‐gold deposit sited in a thrust‐fold belt. Stress mapping illustrates the heterogeneity of stress distribution in the complex structural geometry of the deposit, and predicts the preferential siting of ore zones around the intersections of more steeply dipping, linking thrusts and banded iron‐formation units, and below the controlling more gently dipping basal thrust, the Sunrise Shear. The Howley Anticline in the Pine Creek block hosts several Palaeoproterozoic gold deposits, sited in complex anticlinal structures in greywacke sequences. Stress mapping indicates that gold ores should develop in the hinge zones of symmetrical anticlines, in the hinge zones and more steeply dipping to overturned limbs of asymmetric anticlines, and in and around thrusts in both anticlines and parasitic synclines. The strong correlation between the predictions of the stress mapping, based on the distribution of low σ₃, and the location of gold ores emphasises the potential of stress mapping in cross‐section, not only as an exploration tool for the discovery of additional resources or deposits, but also as a test of geological models. Knowledge of the potential siting of gold ores and their probable orientations also provides a guide to drilling strategies in both mine‐ and regional‐scale exploration.     

The Younger Dryas in the English Lake District: reconciling geomorphological evidence with numerical model outputs

The geomorphology of the south‐western and central Lake District, England is used to reconstruct the mountain palaeoglaciology pertaining to the Lateglacial and Younger Dryas. Limitations to previous ice‐mass reconstructions and consequent palaeoclimatic inferences include: (i) the use of static (steady‐state) glacier reconstructions, (ii) the assumption of a single‐stage Younger Dryas advance, (iii) greatly varying ice‐volume estimates, (iv) inexplicable spatial variations in ELA (Equilibrium Line Altitude), and (v) a lack of robust extent chronology. Here we present geomorphological mapping based on aerial photography and the NextMap Britain Digital Elevation Model, checked by ground survey. Former glacier extents were inferred and ELAs were calculated using the Balance Ratio method of Osmaston. Independently, a time‐dependant 2‐D ice‐flow model was forced by a regional ELA history that was scaled to the GRIP record. This provided a dynamic reconstruction of a mountain ice field that allowed for non‐steady‐state glacier evolution. Fluctuations in climate during the Younger Dryas resulted in multiple glacial advance positions that show agreement with the location of mapped moraines, and may further explain some of the ELA variations found in previous local and static reconstructions. Modelling based on the GRIP record predicts three phases: an initial maximum extent, a middle minor advance or stillstand, and a pronounced but less extensive final advance. The comparisons find that the reconstructions derived from geomorphological evidence are effective representations of steady‐state glacier geometries, but we do propose different extents for some glaciers and, in particular, a large former glacier in Upper Eskdale.     

Correlations between VIMS and RADAR data over the surface of Titan: Implications for Titan’s surface properties

We apply a multivariate statistical method to Titan data acquired by different instruments onboard the Cassini spacecraft. We have searched through Cassini/VIMS hyperspectral cubes, selecting those data with convenient viewing geometry and that overlap with Cassini/RADAR scatterometry footprints with a comparable spatial resolution. We look for correlations between the infrared and microwave ranges the two instruments cover. Where found, the normalized backscatter cross-section obtained from the scatterometer measurement, corrected for incidence angle, and the calibrated antenna temperature measured along with the scatterometry echoes, are combined with the infrared reflectances, with estimated errors, to produce an aggregate data set, that we process using a multivariate classification method to identify homogeneous taxonomic units in the multivariate space of the samples.In medium resolution data (from 20 to 100 km/pixel), sampling relatively large portions of the satellite’s surface, we find regional geophysical units matching both the major dark and bright features seen in the optical mosaic. Given the VIMS cubes and RADAR scatterometer passes considered in this work, the largest homogeneous type is associated with the dark equatorial basins, showing similar characteristics as each other on the basis of all the considered parameters.On the other hand, the major bright features seen in these data generally do not show the same characteristics as each other. Xanadu, the largest continental feature, is as bright as the other equatorial bright features, while showing the highest backscattering coefficient of the entire satellite. Tsegihi is very bright at 5 μm but it shows a low backscattering coefficient, so it could have a low roughness on a regional scale and/or a different composition. Another well-defined region, located southwest of Xanadu beyond the Tui Regio, seems to be detached from the surrounding terrains, being bright at 2.69, 2.78 and 5 μm but having a low radar brightness. In this way, other units can be found that show correlations or anti-correlations between the scatterometric response and the spectrophotometric behavior, not evident from the optical remote sensing data.