Stochastic Structural Modeling

A consistent stochastic model for faults and horizons is described. The faults are represented as a parametric invertible deformation operator. The faults may truncate each other. The horizons are modeled as correlated Gaussian fields and are represented in a grid. Petrophysical variables may be modeled in a reservoir before faulting in order to describe the juxtaposition effect of the faulting. It is possible to condition the realization on petrophysics, horizons, and fault plane observations in wells in addition to seismic data. The transmissibility in the fault plane may also be included in the model. Four different methods to integrate the fault and horizon models in a common model is described. The method is illustrated on an example from a real petroleum field with 18 interpreted faults that are handled stochastically.     

Paleozoic terranes of eastern Australia and the drift history of Gondwana

Critical assessment of Paleozoic paleomagnetic results from Australia shows that paleopoles from locations on the main craton and in the various terranes of the Tasman Fold Belt of eastern Australia follow the same path since 400 Ma for the Lachlan and Thomson superterranes, but not until 250 Ma or younger for the New England superterrane. Most of the paleopoles from the Tasman Fold Belt are derived from the Lolworth-Ravenswood terrane of the Thomson superterrane and the Molong-Monaro terrane of the Lachlan superterrane. Consideration of the paleomagnetic data and geological constraints suggests that these terranes were amalgamated with cratonic Australia by the late Early Devonian. The Lolworth-Ravenswood terrane is interpreted to have undergone a 90° clockwise rotation between 425 and 380 Ma. Although the Tamworth terrane of the western New England superterrane is thought to have amalgamated with the Lachlan superterrane by the Late Carboniferous, geological syntheses suggest that movements between these regions may have persisted until the Middle Triassic. This view is supported by the available paleomagnetic data. With these constraints, an apparent polar wander path for Gondwana during the Paleozoic has been constructed after review of the Gondwana paleomagnetic data. The drift history of Gondwana with respect to Laurentia and Baltica during the Paleozoic is shown in a series of paleogeographic maps.     

Papua New Guinea Highlands: palaeomagnetic constraints on terrane tectonics

Oblique convergence since the Early Cenozoic between the northward-moving Australian plate, westward-moving Pacific plate and almost stationary Eurasian plate has created a world-ranking tectonic zone in the eastern Indonesia–New Guinea–Southwest Pacific region (Tonga–Sulawesi megashear) that is notorious for its complex mix of tectonic styles and terrane juxtapositions. Unlike an ancient analog—the Mesozoic–Cenozoic Cordillera of North America—palaeomagnetic constraints on terrane motions in the zone are few. To improve the framework of quantitative control on such motions and therefore our understanding of the development of the zone, results of a palaeomagnetic study in the Highlands region of Papua New Guinea (PNG), in the southern part of the New Guinea Orogen, are reported. The study yields new insights into terrane tectonics along the Australian craton's active northern margin and confirms the complexity of block rotations to be expected at the local scale in tectonically intricate zones. The study is based on more than 500 samples (21 localities) collected from an interior and an exterior zone of New Guinea's central cordillera. The two zones are separated by the Tahin and Stolle–Lagaip–Kaugel Fault zones and collectively represent the para-autochthonous northern margin of the Australian craton. Samples from the interior zone, which in the study area comprises a cratonic spur of uncertain—probably displaced—origin, come from Triassic to Miocene sediments and subordinate volcanics of the Kubor Anticline, Jimi Terrane, and Yaveufa Syncline (16 localities) in the central and eastern Highlands. Samples from the exterior zone, which represent a basement-involved, Pliocene foreland fold-and-thrust belt, come from Middle Eocene to Middle Miocene carbonates and clastics (five localities) in the southern Highlands of the Papuan Fold Belt. Results permit us to constrain the tectonic evolution of the two zones palaeomagnetically. Using mainly thermal demagnetization techniques, three main magnetic components have been identified in the collection: (1) a recent field overprint of both normal and reverse polarity; (2) a pervasive overprint of mainly normal polarity that originated during extensive Middle to Late Miocene intrusive activity in the central cordillera; and (3) a primary component which has been identified in only 7 of the 21 localities (5 of 11 stratigraphic units represented in the collection). All components show patterns of rotation that are consistent within the zones, but differ between them. In the interior zone (central and eastern Highlands), large-scale counterclockwise rotations of between 30°+ and 100°+ have been established throughout the Kubor Anticline and Jimi Terrane, with some clockwise rotation present in the southern part of the Yaveufa Syncline. In contrast, in the Mendi area of the exterior zone (southern Highlands), clockwise rotations of between 30°+ and 50°+ can be recognized. These contrasting rotation patterns across the Tahin and Stolle–Lagaip–Kaugel Fault zones indicate decoupling of the two tectonic zones, probably along basement-involved faults. The clockwise rotations in the southern Highlands of the Papuan Fold Belt are to be expected from its structural grain, and are probably governed by regional basement faults and transverse lineaments. In contrast, the pattern of counterclockwise rotations in the Kubor Anticline–Jimi Terrane cratonic spur of the central and eastern Highlands was unexpected. The pattern is interpreted to result from non-rigid rotation of continental terranes as they were transported westward across the northeastern margin of the Australian craton. This margin became reorganised after the Middle Miocene, when the steadily northward-advancing Australian craton impinged into the westward-moving Pacific plate/buffer-plate system. Transpressional reorganisation under the influence of the sinistral Tonga–Sulawesi megashear became enhanced with Mio-Pliocene docking, and subsequent southward overthrusting, of the Finisterre Terrane onto the northeastern margin of the Australian craton.     

Palaeomagnetic results from the ca. 1130 Ma Borgmassivet intrusions in the Ahlmannryggen region of Dronning Maud Land, Antarctica, and tectonic implications

Detailed palaeomagnetic and rock magnetic analyses provide improved palaeomagnetic results from 23 sites in the Borgmassivet intrusions in the Ahlmannryggen region of Dronning Maud Land, East Antarctica. These intrusions are of similar age to their host, the ca. 1130 Ma Ritscherflya Supergroup (RSG). A mean direction of D=235.4°, I=−7.6° with k=45.9 and α₉₅=4.5° was obtained from this study. When combined with previously reported results from 11 sites in the same region, including sites from the Ritscherflya Supergroup, it gives an overall mean direction for 34 sites from the igneous suite with D=236.5°, I=−3.6°, k=27.9 and α₉₅=4.8°. Isothermal remanent magnetization (IRM) experiments on several specimens suggest magnetite or titanomagnetite as the primary remanence carrier, while high temperature magnetic susceptibility experiments indicate the presence of single domain particles. These observations, together with field evidence and the high coercivities and unblocking temperatures, support a primary origin for the observed characteristic remanence. The Borgmassivet palaeomagnetic pole lies at 54.5°E, 8.3°N with A₉₅=3.3°. If Antarctica is moved to its Gondwanan position adjacent to southeast Africa, the Borgmassivet pole (BM) coincides with that of the African well-established, well-dated (1100 Ma) Umkondo Large Igneous Province pole, supporting the hypothesis that the Grunehogna craton of Dronning Maud Land was part of the Kalahari craton of southern Africa at ca. 1100 Ma.     

Book Review: Experiments: Planning,Analysis, and Parameter Design Optimization. By C. F. Jeff Wu and Michael Hamada. John Wiley & Sons,Inc., New York, 2000, 630 p., U.S.$89.95 (hb), ISBN 0-471-25511-4

Mathematical Geosciences -     

Control‐volume mixed finite element methods

A key ingredient in simulation of flow in porous media is accurate determination of the velocities that drive the flow. Large‐scale irregularities of the geology (faults, fractures, and layers) suggest the use of irregular grids in simulation. This paper presents a control‐volume mixed finite element method that provides a simple, systematic, easily implemented procedure for obtaining accurate velocity approximations on irregular (i.e., distorted logically rectangular) block‐centered quadrilateral grids. The control‐volume formulation of Darcy’s law can be viewed as a discretization into element‐sized “tanks” with imposed pressures at the ends, giving a local discrete Darcy law analogous to the block‐by‐block conservation in the usual mixed discretization of the mass‐conservation equation. Numerical results in two dimensions show second‐order convergence in the velocity, even with discontinuous anisotropic permeability on an irregular grid. The method extends readily to three dimensions.     

Formation of silica oncoids around geysers and hot springs at El Tatio, northern Chile

Siliceous oncoids, up to 4 cm in diameter, are common on the laterally extensive sinter aprons that surround the spectacular geysers and hot springs at El Tatio in northern Chile. Many of these complex oncoids developed close to geyser and spring vents that discharge boiling water. Internally the oncoids, which are composed of precipitated amorphous silica, are formed of complex arrays of spicules and concentric laminae as well as detrital volcanic grains. Spicular growth is dominant in most examples. The formation and growth of the spicules and concentric laminae were mediated by a microbial community which included filamentous microbes, mucus, and possibly bacteria. The microbes and mucus were silicified by replacement and encrustation. In some laminae the filamentous microbes lay parallel to the growth surface; in other laminae most filaments forming the thin mats were suberect. Amorphous silica precipitated between the filaments occluded porosity and commonly disguised the microbial fabric. The oncoids grew on the proximal sinter aprons around the geyser vents and hot spring pools. Most growth took place subaerially with the silica delivered to the precipitation sites by splashing water from the geysers and/or periodic shallow flooding of the discharge aprons. Unlike silica oncoids at other geothermal sites, vertical growth of oncoids that formed in some rimstone pools was not limited by water depth.