Hydrothermal control on organic matter alteration and illite precipitation, Mt Isa Basin, Australia

Abundant illite precipitation, in Proterozoic rocks from Northern Lawn Hill Platform, Mt Isa Basin, Australia, occurred in organic matter-rich black shales rather than in sandstones, siltstones and organic matter-poor shales. Sandstones and siltstones acted as impermeable rocks, as early diagenetic quartz and carbonate minerals reduced the porosity-permeability. Scanning and transmission electron microscopy (SEM and TEM) studies indicate a relation between creation of microporosity-permeability and organic matter alteration, suitable for subsequent mineral precipitation. K-Ar data indicate that organic matter alteration and the subsequent illite precipitation within the organic matter occurred during the regional hydrothermal event at 1172 +/- 150 (2sigma) Ma. Hot circulating fluids are considered to be responsible for organic matter alteration, migration and removal of volatile hydrocarbon, and consequently porosity-permeability creation. Those rocks lacking sufficient porosity-permeability, such as sandstones, siltstones and organic matter poor shales, may not have been affected by fluid movement. In hydrothermal systems, shales and mudstones may not be impermeable as usually assumed because of hydrocarbons being rapidly removed by fluid, even with relatively low total organic carbon.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 2: Characteristics and petrogenesis of the high-Ti latite and high-Ti and low-Ti voluminous quartz latite eruptives

As a result of their relative concentration towards the respective Atlantic margins, the silicic eruptives of the Parana (Brazil)-Etendeka large igneous province are disproportionately abundant in the Etendeka of Namibia. The NW Etendeka silicic units, dated at similar to132 Ma, occupy the upper stratigraphic levels of the volcanic sequences, restricted to the coastal zone, and comprise three latites and five quartz latites (QL). The large-volume Fria QL is the only low-Ti type. Its trace element and isotopic signatures indicate massive crustal input. The remaining NW Etendeka silicic units are enigmatic high-Ti types, geochemically different from low-Ti types. They exhibit chemical affinities with the temporally overlapping Khumib high-Ti basalt (see Ewart et al. Part 1) and high crystallization temperatures (greater than or equal to980 to 1120degreesC) inferred from augite and pigeonite phenocrysts, both consistent with their evolution from a mafic source. Geochemically, the high-Ti units define three groups, thought genetically related. We test whether these represent independent liquid lines of descent from a common high-Ti mafic parent. Although the recognition of latites reduces the apparent silica gap, difficulty is encountered in fractional crystallization models by the large volumes of two QL units. Numerical modelling does, however, support large-scale open-system fractional crystallization, assimilation of silicic to basaltic materials, and magma mixing, but cannot entirely exclude partial melting processes within the temporally active extensional environment. The fractional crystallization and mixing signatures add to the complexity of these enigmatic and controversial silicic magmas. The existence, however, of temporally and spatially overlapping high-Ti basalts is, in our view, not coincidental and the high-Ti character of the silicic magmas ultimately reflects a mantle signature.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 1: Introduction, mafic lavas and re-evaluation of mantle source components

The bimodal NW Etendeka province is located at the continental end of the Tristan plume trace in coastal Namibia. It comprises a high-Ti (Khumib type) and three low-Ti basalt (Tafelberg, Kuidas and Esmeralda types) suites, with, at stratigraphically higher level, interstratified high-Ti latites (three units) and quartz latites (five units), and one low-Ti quartz latite. Khumib basalts are enriched in high field strength elements and light rare earth elements relative to low-Ti types and exhibit trace element affinities with Tristan da Cunha lavas. The unradiogenic Pb-206/Pb-204 ratios of Khumib basalts are distinctive, most plotting to the left of the 132 Ma Geochron, together with elevated Pb-207/Pb-204 ratios, and Sr-Nd isotopic compositions plotting in the lower Nd-143/Nd-144 part of mantle array (EM1-like). The low-Ti basalts have less coherent trace element patterns and variable, radiogenic initial Sr (similar to0.707-0.717) and Pb isotope compositions, implying crustal contamination. Four samples, however, have less radiogenic Pb and Sr that we suggest approximate their uncontaminated source. All basalt types, but particularly the low-Ti types, contain samples with trace element characteristics (e.g. Nb/Nb-*) suggesting metasediment input, considered source-related. Radiogenic isotope compositions of these samples require long-term isolation of the source in the mantle and depletions (relative to unmodified sediment) in certain elements (e.g. Cs, Pb, U), which are possibly subduction-related. A geodynamic model is proposed in which the emerging Tristan plume entrained subducted material in the Transition Zone region, and further entrained asthenosphere during plume head expansion. Mixing calculations suggest that the main features of the Etendeka basalt types can be explained without sub-continental lithospheric mantle input. Crustal contamination is evident in most low-Ti basalts, but is distinct from the incorporation of a metasedimentary source component at mantle depths.

On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

Recent development in the mechanics of superplasticity and its applications

This paper reviews the recent developments in the mechanics of superplasticity and its applications in industrial practice. After introducing the phenomena of superplasticity, the basic experiments for determining material deformation behavior and related parameters, and constructing superplastic constitutive equations, are reviewed. Finite element related formulations and techniques for simulating superplastic forming are discussed, together with some practical applications. (C) 2004 Elsevier B.V. All rights reserved.

Accelerating precursory activity within a class of earthquake analogue automata

A statistical fractal automaton model is described which displays two modes of dynamical behaviour. The first mode, termed recurrent criticality, is characterised by quasi-periodic, characteristic events that are preceded by accelerating precursory activity. The second mode is more reminiscent of SOC automata in which large events are not preceded by an acceleration in activity. Extending upon previous studies of statistical fractal automata, a redistribution law is introduced which incorporates two model parameters: a dissipation factor and a stress transfer ratio. Results from a parameter space investigation indicate that a straight line through parameter space marks a transition from recurrent criticality to unpredictable dynamics. Recurrent criticality only occurs for models within one corner of the parameter space. The location of the transition displays a simple dependence upon the fractal correlation dimension of the cell strength distribution. Analysis of stress field evolution indicates that recurrent criticality occurs in models with significant long-range stress correlations. A constant rate of activity is associated with a decorrelated stress field.

Late pleistocene to holocene composite speleothem O-18 and C-13 chronologies from south island, new Zealand-did a global younger dryas really exist?

Oxygen and carbon data from eight stalagmites from northwest South Island are combined to produce composite records of delta(18)O and delta(13)C from 23.4 ka to the present. The chronology is anchored by 43 thermal ionization mass spectrometry (TIMS) uranium series ages. Delta O-18 values are interpreted as having a first order positive relationship to temperature, but also to be influenced by precipitation in a complex manner. Delta C-13 is interpreted as responding negatively to increases in atmospheric CO, concentration, biological activity and precipitation amount. Six climatic phases are recognized. After adjustment of 1.2parts per thousand for the ice volume effect, the delta(18)O record between 23 and 18 ka varies around -3.72parts per thousand compared to the Holocene average of -3.17parts per thousand. Late-glacial warming commenced between 18.2 and 17.8 ka and accelerated after 16.7 ka, culminating in a positive excursion between 14.70 and 13.53 ka. This was followed by a significant negative excursion between 13.53 and 11.14 ka of up to 0.55parts per thousand depth that overlapped the Antarctic Cold Reversal (ACR) and spanned the Younger Dryas (YD). Positive delta(18)O excursions at 11.14 ka and 6.91-6.47 ka represent the warmest parts of the Holocene. The mid-Holocene from 6 to 2 ka was marked by negative excursions that coincide with increased glacial activity in the South Island. A short positive excursion from 0.71 to 0.57 ka was slightly later than the Medieval Warm Period of Europe. Delta C-13 values were high until 17.79 ka after which there was an abrupt decrease to 17.19 ka followed by a steady decline to a minimum at 10.97 ka. Then followed a general increase, suggesting a drying trend, to 3.23 ka followed by a further general decline. The abrupt decrease in delta-values after 17.79 ka probably corresponds to an increase in atmospheric CO2 concentration, biological activity and wetness at the end of the Last Glaciation, but the reversal identified in the delta(18)O record from 13.53 to 11.14 ka was not reflected in delta(13)C changes. The lowest delta(13)C values coincided with the early Holocene climatic suboptimum when conditions were relatively wet as well as mild. Major trends in the delta(18)O(c) record are similar to the Northern Hemisphere, but second order detail is often distinctly different. Consequently, at the millennial scale, a more convincing case can be made for asymmetric climatic response between the two hemispheres rather than synchronicity. (C) 2004 Elsevier B.V. All rights reserved.

Volcanic resurfacing and the early terrestrial crust: Zircon U-Pb and REE constraints from the Isua Greenstone Belt, southern West Greenland

The oldest known bona fide succession of elastic metasediments Occurs in the Isua Greenstone Belt. SW Greenland and consists of a variety of mica schists and rare metaconglomerates. The metasediments are in direct contact with a felsic metavolcanic lithology that has previously been dated to 3.71 Ga. Based on trace element geochemical data for 30 metasediments, we selected the six samples with highest Zr concentrations for zircon extraction. These samples all yielded very few or no zircon, Those extracted from mica schists yielded ion probe U/Pb ages between 3.70 and 3,71 Ga. One metaconglomerate sample yielded just a single zircon of 3.74 Ga age. The mica schist hosted zircons have U/Pb ages. Th/U ratios, REE patterns and Eu anomalies indistinguishable from zircon in the adjacent 3.71 Ga felsic metavolcanic unit. Trace element modelling requires the bulk of material in the metasediments to be derived from variably weathered mafic lithologies but some metasediments contain substantial contribution from more evolved source lithologies. The paucity of zircon in the mica schists is thus explained by incorporation of material from largely zircon-free volcanic lithologies. The absence of older zircon in the mica schists and the preponderance of mafic source material imply intense, mainly basaltic resurfacing of the early Earth. The implications of this process are discussed, Thermal considerations suggest that horizontal growth of Hadean crust by addition of mafic ultramafic lavas must have triggered self-reorganisation of the protocrust by remelting. Reworking oft Hadean crust may have been aided by burial of hydrated (weathered) metabasalt due to semi-continuous addition of new voluminous basalt Outpouring,;, This process Causes a bias towards eruption of Zr-saturated partial melts at the surface with O-isotope corn posit ion,, potentially different from the mantle. The oldest zircons hosted in sediments would have been buried to substantial depth or formed in plutons that crystallised at some depth from which it took hundreds of millions of years for them to be exhumed and incorporated into much younger sediments. (C) 2005 Elsevier B.V.All rights reserved.

Clay mineralogical, geochemical and isotopic tracing of the evolution of the Woodleigh impact structure, Southern Carnarvon Basin, Western Australia

Chaotically structured diamictite from the inner ring syncline surrounding the central uplift of the Woodleigh impact structure contains shocked metamorphic and impact melt-rock fragments, largely derived from Ordovician and Devonian target sandstones. Coarse illite fractions (< 2 mu m) from the sandstones containing no K-feldspar yield K-Ar ages of around 400 Ma, whereas the K-Ar ages of authigenic clays of > 0.2 mu m fractions from the diamictite without smectite and K-feldspar cluster around 360 Ma, consistent with Rb-Sr data. Crystallisation of newly formed illite in the impact melt rock clasts and recrystallisation of earlier formed illite in the sandstone clasts preserved in the diamictite, are attributed to impact-induced hydrothermal processes in the Late Devonian. The illitic clays from the diamictite and from the sandstones have very similar trace element compositions, with significantly enriched incompatible lithophile elements, which increase in concentrations correlatively with those of the compatible ferromagnesian elements. The unusual trace element associations in the clays may be due to the involvement of hot gravity-driven basinal fluids that interacted with rocks of the Precambrian craton to the east of the study area, or with such material transported and reworked in the studied sedimentary succession.

Integrated Pb- and S-isotope investigation of sulphide minerals from the early Archaean of southwest Greenland

We present high-spatial resolution secondary ion mass spectrometry (SIMS) measurements of Pb and S isotopes in sulphides from early Archaean samples at two localities in southwest Greenland. Secondary pyrite from a 3.71 Ga sample of magnetite-quartz banded iron formation in the Isua Greenstone Belt, which has previously yielded unradiogenic Pb consistent with its ancient origin, contains sulphur with a mass independently fractionated (MIF) isotope signature (Delta(33)S =+3.3 parts per thousand). This reflects the secondary mineralization of remobilized sedimentary S carrying a component modified by photochemical reactions in the early Archaean atmosphere. It further represents one of the most extreme positive excursions so far known from the early Archaean rock record. Sulphides from a quartz-pyroxene rock and an ultramafic boudin from the island of Akilia, in the Godth (a) over circle bsfjord, have heterogeneous and generally radiogenic Pb isotopic compositions that we interpret to represent partial re-equilibration of Pb between the sulphides and whole rocks during tectonothermal events at 3.6, 2.7 and 1.6 Ga. Both these samples have Delta(33)S=0 (within analytical error) and therefore show no evidence for MIF sulphur. These data are consistent with previous interpretations that the rock cannot be proven to have a sedimentary origin. Our study illustrates that SIMS S-isotope measurements in ancient rocks can be used to elucidate early atmospheric parameters because of the ability to obtain combined S and Pb-isotope data, but caution must be applied when using such data to infer protolith. When information from geological context, petrography and chronology (i.e. by Pb isotopes) is combined and fully evaluated, Delta(33)S signatures from sulphides and their geological significance can be interpreted with a higher degree of confidence. (c) 2005 Elsevier B.V All rights reserved.

Review of tropospheric bomb C-14 data for carbon cycle modeling and age calibration purposes

Comprehensive published radiocarbon data from selected atmospheric records, tree rings, and recent organic matter were analyzed and grouped into 4 different zones (three for the Northern Hemisphere and one for the whole Southern Hemisphere). These C-14 data for the summer season of each hemisphere were employed to construct zonal, hemispheric, and global data sets for use in regional and global carbon model calculations including calibrating and comparing carbon cycle models. In addition, extended monthly atmospheric C-14 data sets for 4 different zones were compiled for age calibration purposes. This is the first time these data sets were constructed to facilitate the dating of recent organic material using the bomb C-14 curves. The distribution of bomb C-14 reflects the major zones of atmospheric circulation.

Theoretical and numerical analysis of large-scale heat transfer problems with temperature-dependent pore-fluid densities

Purpose - In many scientific and engineering fields, large-scale heat transfer problems with temperature-dependent pore-fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large-scale heat transfer problems with temperature-dependent pore-fluid densities in the lithosphere and crust scales. Design/methodology/approach - The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts. Findings - From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive-and-advective lithosphere with variable pore-fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere. Originality/value - The present analytical solutions can be used to: validate numerical methods for solving large-scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large-scale heat transfer problems with temperature-dependent fluid densities.

Geogenic sources and sinks of trace metals in the Larsemann Hills, East Antarctica: Natural processes and human impact

To evaluate the extent of human impact on a pristine Antarctic environment, natural baseline levels of trace metals have been established in the basement rocks of the Larsemann Hills, East Antarctica. From a mineralogical and geochemical point of view the Larsemann Hills basement is relatively homogeneous, and contains high levels of Pb, Th and U. These may become soluble during the relatively mild Antarctic summer and be transported to lake waters by surface and subsurface melt water. Melt waters may also be locally enriched in V, Cr, Co, Ni, Zn and Sri derived from weathering of metabasite pods. With a few notable exceptions, the trace metal concentrations measured in the Larsemann Hills lake waters can be entirely accounted for by natural processes such as sea spray and surface melt water input. Thus, the amount of trace metals released by weathering of basement lithologies and dispersed into the Larsemann Hills environment, and presumably in similar Antarctic environments, is, in general, not negligible, and may locally be substantial. The Larsemann Hills sediments are coarse-grained and contain minute amounts of clay-size particles, although human activities have contributed to the generation of fine-grained material at the most impacted sites. Irrespective of their origin, these small amounts of fine-grained clastic sediments have a relatively small surface area and charge, and are not as effective metal sinks as the abundant, thick cyanobacterial algal mats that cover the lake floors. Thus, the concentration of trace metals in the Larsemann Hills lake waters is regulated by biological activity and thawing-freezing cycles, rather than by the type and amount of clastic sediment supply. (c) 2005 Elsevier Ltd. All rights reserved.