3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Stable isotope and trace element stratigraphy across the Permian-Triassic transition: A redefinition of the boundary in the Velebit Mountain, Croatia

Stable isotopes of carbonates (delta(13)C(carb), delta(18)O(carb)), organic matter (delta(13)C(org), delta(15)N(org)) and major, trace and rare earth element (REE) compositions of marine carbonate rocks of Late Permian to Early Triassic age were used to establish the position of the Permian-Triassic boundary (PTB) at two continuous sections in the Velebit Mountain, Croatia. The chosen sections - Rizvanusa and Brezimenjaca - are composed of two lithostratigraphic units, the Upper Permian Transitional Dolomite and the overlying Sandy Dolomite. The contact between these units, characterized by the erosional features and sudden occurrence of ooids and siliciclastic grains, was previously considered as the chronostratigraphic PTB. The Sandy Dolomite is characterized by high content of non-carbonate material (up to similar to 30 wt.% insoluble residue), originated from erosion of the uplifted hinterland. A relatively rich assemblage of Permian fossils (including Geinitzina, Globivalvulina, Hemigordius, bioclasts of gastropods, ostracods and brachiopods) was found for the first time in Sandy Dolomite, 5 m above the lithologic boundary in the Rizvanusa section. A rather abrupt negative delta(13)C(carb) excursion in both sections appears in rocks showing no recognizable facies change within the Sandy Dolomite, -2 parts per thousand at Rizvanusa and -1.2 parts per thousand at Brezimenjaca, 11 m and 0.2 m above the lithologic contact, respectively. This level within the lower part of the Sandy Dolomite is proposed as the chemostratigraphic PTB. In the Rizvanusa section, the delta(13)C(org) values decline gradually from similar to-25 parts per thousand in the Upper Permian to similar to-29 parts per thousand in the Lower Triassic. The first negative delta(13)C(org) excursion occurs above the lithologic contact, within the uppermost Permian deposits, and appears to be related to the input of terrigenous material. The release of isotopically light microbial soil-biomass into the shallow-marine water may explain this sudden decrease of delta(13)C(org) values below the PTB. This would support the hypothesis that in the western Tethyan realm the land extinction, triggering a sudden drop of woody vegetation and related land erosion, preceded the marine extinction. The relatively low delta(15)N(org) values at the Permian-Triassic (P-Tr) transition level, close to approximate to 0 parts per thousand, and a secondary negative delta(13)C(org) excursion of -0.5 parts per thousand point to significant terrestrial input and primary contribution of cyanobacteria. The profiles of the concentrations of redox-sensitive elements (Ce, Mn, Fe, V), biogenic or biogenic-scavenged elements (P, Ba, Zn, V), Ce/Ce* values, and normalized trace elements, including Ba/Al, Ba/Fe, Ti/Al, Al/(Al + Fe + Mn) and Mn/Ti show clear excursions at the Transitional Dolomite-Sandy Dolomite lithologic boundary and the chemostratigraphic P-Tr boundary. The stratigraphic variations indicate a major regression phase marking the lithologic boundary, transgressive phases in the latest Permian and a gradual change into shallow/stagnant anoxic marine environment towards the P-Tr boundary level and during the earliest Triassic. (C) 2010 Elsevier B.V. All rights reserved.

Carbon and oxygen isotope study of hydrothermal carbonates in the zinc-lead deposits of the San Vicente district, central Peru: A quantitative modeling on mixing processes and CO2 degassing

Carbon and oxygen isotope studies of the host and gangue carbonates of Mississippi Valley-type zinc-lead deposits in the San Vicente District hosted in the Upper Triassic to Lower Jurassic dolostones of the Pucara basin (central Peru) were used to constrain models of the ore formation. A mixing model between an incoming hot saline slightly acidic radiogenic (Pb, Sr) fluid and the native formation water explains the overall isotopic variation (delta(13)C = - 11.5 to + 2.5 parts per thousand relative to PDB and delta(18)O = + 18.0 to + 24.3 parts per thousand relative to SMOW) of the carbonate generations. The dolomites formed during the main ore stage show a narrower range (delta(13)C = - 0.1 to + 1.7 parts per thousand and delta(18)O = + 18.7 to + 23.4 parts per thousand) which is explained by exchange between the mineralizing fluids and the host carbonates combined with changes in temperature and pressure. This model of fluid-rock interaction explains the pervasive alteration of the host dolomite I and precipitation of sphalerite I. The open-space filling hydrothermal white sparry dolomite and the coexisting sphalerite II formed by prolonged fluid-host dolomite interaction and limited CO2 degassing. Late void-filling dolomite III (or calcite) and the associated sphalerite III formed as the consequence of CO2 degassing and concomitant pH increase of a slightly acidic ore fluid. Widespread brecciation is associated to CO2 outgassing. Consequently, pressure variability plays a major role in the ore precipitation during the late hydrothermal events in San Vicente. The presence of native sulfur associated with extremely carbon-light calcites replacing evaporitic sulfates (e.g., delta(13)C = - 11.5 parts per thousand), altered native organic matter and heavier hydrothermal bitumen (from - 27.0 to - 23.0 parts per thousand delta(13)C) points to thermochemical reduction of sulfate and/or thiosulfate. The delta(13)C- and delta(18)O-values of the altered host dolostone and hydrothermal carbonates, and the carbon isotope composition of the associated organic matter show a strong regional homogeneity. These results coupled with the strong mineralogical and petrographic similarities of the different MVT occurrences perhaps reflects the fact that the mineralizing processes were similar in the whole San Vicente belt, suggesting the existence of a common regional mineralizing hydrothermal system with interconnected plumbing.

Geological and Geomorphological Cartography at Scale 1: 10.000: a Tool of Management for Towns

Geological and geomorphological mapping at scale 1:10.000 besides from being an important source of scientific information it is also a necessary tool for municipal organs in order to make proper decisions when dealing with geo-environmental problems concerning integral territorial development. In this work, detailed information is given on the contents of such maps, their social and economical application, and a balance of the investment and gains that derives from them

Mid Holocene in Sant Julià de Boada (Baix Empordà, Girona): Stratigraphical sequence and paleoenvironmental records of the paludal deposits

Descripció de la seqüència estratigràfica i dels registres paleoambientals dels sediments holocens de Sant Julià de Boada

Quantitative analysis of crystal/grain sizes and their distributions in 2D and 3D

We review methods to estimate the average crystal (grain) size and the crystal (grain) size distribution in solid rocks. Average grain sizes often provide the base for stress estimates or rheological calculations requiring the quantification of grain sizes in a rock's microstructure. The primary data for grain size data are either 1D (i.e. line intercept methods), 2D (area analysis) or 3D (e.g., computed tomography, serial sectioning). These data have been used for different data treatments over the years, whereas several studies assume a certain probability function (e.g., logarithm, square root) to calculate statistical parameters as the mean, median, mode or the skewness of a crystal size distribution. The finally calculated average grain sizes have to be compatible between the different grain size estimation approaches in order to be properly applied, for example, in paleo-piezometers or grain size sensitive flow laws. Such compatibility is tested for different data treatments using one- and two-dimensional measurements. We propose an empirical conversion matrix for different datasets. These conversion factors provide the option to make different datasets compatible with each other, although the primary calculations were obtained in different ways. In order to present an average grain size, we propose to use the area-weighted and volume-weighted mean in the case of unimodal grain size distributions, respectively, for 2D and 3D measurements. The shape of the crystal size distribution is important for studies of nucleation and growth of minerals. The shape of the crystal size distribution of garnet populations is compared between different 2D and 3D measurements, which are serial sectioning and computed tomography. The comparison of different direct measured 3D data; stereological data and direct presented 20 data show the problems of the quality of the smallest grain sizes and the overestimation of small grain sizes in stereological tools, depending on the type of CSD. (C) 2011 Published by Elsevier Ltd.

Globally aligned photomosaic of the Lucky Strike hydrothermal vent field (Mid-Atlantic Ridge, 37°18.5’N): Release of georeferenced data, mosaic construction, and viewing software

We present a georeferenced photomosaic of the Lucky Strike hydrothermal vent field (Mid-Atlantic Ridge, 37°18’N). The photomosaic was generated from digital photographs acquired using the ARGO II seafloor imaging system during the 1996 LUSTRE cruise, which surveyed a ~1 km2 zone and provided a coverage of ~20% of the seafloor. The photomosaic has a pixel resolution of 15 mm and encloses the areas with known active hydrothermal venting. The final mosaic is generated after an optimization that includes the automatic detection of the same benthic features across different images (feature-matching), followed by a global alignment of images based on the vehicle navigation. We also provide software to construct mosaics from large sets of images for which georeferencing information exists (location, attitude, and altitude per image), to visualize them, and to extract data. Georeferencing information can be provided by the raw navigation data (collected during the survey) or result from the optimization obtained from imatge matching. Mosaics based solely on navigation can be readily generated by any user but the optimization and global alignment of the mosaic requires a case-by-case approach for which no universally software is available. The Lucky Strike photomosaics (optimized and navigated-only) are publicly available through the Marine Geoscience Data System (MGDS, http://www.marine-geo.org). The mosaic-generating and viewing software is available through the Computer Vision and Robotics Group Web page at the University of Girona (http://eia.udg.es/_rafa/mosaicviewer.html)

Application of K-feldspar-jadeite-quartz barometry to eclogite facies metagranites and metapelites in the Sesia Lanzo Zone (Western Alps, Italy)

The eclogite facies assemblage K-feldspar-jadeite-quartz in metagranites and metapelites from the Sesia-Lanzo Zone (Western Alps, Italy) records the equilibration pressure by dilution of the reaction jadeite + quartz = albite. The metapelites show partial transformation from a pre-Alpine assemblage of garnet (Alm(63)Prp(26)Grs(10))-K-feldspar-plagioclase-biotite +/- sillimanite to the Eo-Alpine high-pressure assemblage garnet (Alm(50)Prp(14)Grs(35))-jadeite (Jd(80-97)Di(0-4)Hd(0-8)Acm(0-7))=zoisite-phengite. Plagioclase is replaced by jadeite-zoisite-kyanite-K-feldspar-quartz and biotite is replaced by garnet-phengite or omphacite-kyanite-phengite. Equilibrium was attained only in local domains in the metapelites and therefore the K-feldspar-jadeite-quartz (KJQ) barometer was applied only to the plagioclase pseudomorphs and K-feldspar domains. The albite content of K-feldspar ranges from 4 to 11 mol% in less equilibrated assemblages from Val Savenca and from 4 to 7 mol% in the partially equilibrated samples from Monte Mucrone and the equilibrated samples from Montestrutto and Tavagnasco. Thermodynamic calculations on the stability of the assemblage K-feldspar-jadeite-quartz using available mixing data for K-feldspar and pyroxene indicate pressures of 15-21 kbar (+/- 1.6-1.9 kbar) at 550 +/- 50 degrees C. This barometer yields direct pressure estimates in high-pressure rocks where pressures are seldom otherwise fixed, although it is sensitive to analytical precision and the choice of thermodynamic mixing model for K-feldspar. Moreover, the KJQ barometer is independent of the ratio P-H2O/P-T. The inferred limiting a(H2O) for the assemblage jadeite-kyanite in the metapelites from Val Savenca is low and varies from 0.2 to 0.6.

Metallogenic model of the Trepca Pb-Zn-Ag skarn deposit, Kosovo: Evidence from fluid inclusions, Rare Earth Elements, and stable isotope data

The Trepca Pb-Zn-Ag skarn deposit (29 Mt of ore at 3.45% Pb, 2.30% Zn, and 80 g/t Ag) is located in the Kopaonik block of the western Vardar zone, Kosovo. The mineralization, hosted by recrystallized limestone of Upper Triassic age, was structurally and lithologically controlled. Ore deposition is spatially and temporally related with the postcollisional magmatism of Oligocene age (23-26 Ma). The deposit was formed during two distinct mineralization stages: an early prograde closed-system and a later retrograde open-system stage. The prograde mineralization consisting mainly of pyroxenes (Hd(54-100)Jo(0-45)Di(0-45)) resulted from the interaction of magmatic fluids associated with Oligocene (23-26 Ma) postcollisional magmatism. Whereas there is no direct contact between magmatic rocks and the mineralization, the deposit is classified as a distal Pb-Zn-Ag skarn. Abundant pyroxene reflects low oxygen fugacity (<10(-31) bar) and anhydrous environment. Fluid inclusion data and mineral assemblage limit the prograde stage within a temperature range between 390 degrees and 475 degrees C. Formation pressure is estimated below 900 bars. Isotopic composition of aqueous fluid, inclusions hosted by hedenbergite (delta D = -108 to -130 parts per thousand; delta O-18 = 7.5-8.0 parts per thousand), Mn-enriched mineralogy and high REE content of the host carbonates at the contact with the skarn mineralization suggest that a magmatic fluid was modified during its infiltration through the country rocks. The retrograde mineral assemblage comprises ilvaite, magnetite, arsenopyrite, pyrrhotite, marcasite, pyrite, quartz, and various carbonates. Increases in oxygen and sulfur fugacities, as well as a hydrous character of mineralization, require an open-system model. The opening of the system is related to phreatomagmatic explosion and formation of the breccia. Arsenopyrite geothermometer limits the retrograde stage within the temperature range between 350 degrees and 380 degrees C and sulfur fugacity between 10(-8.8) and 10(-7.2) bars. The principal ore minerals, galena, sphalerite, pyrite, and minor chalcopyrite, were deposited from a moderately saline Ca-Na chloride fluid at around 350 degrees C. According to the isotopic composition of fluid inclusions hosted by sphalerite (delta D = -55 to -74 parts per thousand; delta O-18 = -9.6 to -13.6 parts per thousand), the fluid responsible for ore deposition was dominantly meteoric in origin. The delta S-31 values of the sulfides spanning between -5.5 and +10 parts per thousand point to a magmatic origin of sulfur. Ore deposition appears to have been largely contemporaneous with the retrograde stage of the skarn development. Postore stage accompanied the precipitation of significant amount of carbonates including the travertine deposits at the deposit surface. Mineralogical composition of travertine varies from calcite to siderite and all carbonates contain significant amounts of Mn. Decreased formation temperature and depletion in the REE content point to an influence of pH-neutralized cold ground water and dying magmatic system.

Mineralogical, K-Ar, stable and Sr isotope systematics of K-white micas during very low-grade metamorphism of limestones (Helvetic nappes, Western Switzerland)

Mineralogical, K-Ar, Rb-Sr and stable isotope analyses have been carried out on K-white micas from Helvetic Malm limestones in order to examine their evolution during very low- to low-grade Alpine metamorphism, associated with intense ductile deformation. Metamorphic temperatures were estimated al approximately 300-degrees-C from stable isotopes (quartz-calcite thermometry), occurrence of chloritoid, and `'epizonal'' illite crystallinity index. K-white micas consist of variable mixtures of 2M, phengite and muscovite, as revealed by detailed X-ray diffraction analyses using peak decomposition of the (060, 331) spectra. K-Ar apparent ages display a strong grain-size dependence in which mainly fine-grained size fractions (< 2 mum) record Alpine ages (37-15 Ma). However, these ages provide a relative rather than an absolute chronology of the diachronous Alpine metamorphic evolution of the Helvetic nappes. The resetting of the K-Ar isotopic system of K-white micas to Alpine metamorphic conditions reflects an apparent combination of crystallization/recrystallization and radiogenic Ar-40 diffusion loss. The oxygen isotope compositions of micas (+ 15 to + 22 parts per thousand) are intermediate between detrital and O-18-enriched values expected for micas neoformed within an abundant marine carbonate matrix. No isotopic equilibrium has been reached between calcite and micas. The variable depletion of hydrogen isotope compositions (- 126 to - 82 parts per thousand) is influenced by the interaction with organic matter under closed-system conditions. Organic matter, if not removed, may also represent a serious source of error in K-Ar age determination, by introducing radiogenic Ar-40 contamination. Sr-87/Sr-86 isotope ratios of micas range from 0.70879 to 0.70902 with one outlier at 0.71794. The low values reflect Sr exchange with calcite occurring during crystallization/recrystallization of micas under closed-system conditions.

The sandstone-hosted Beverley uranium Deposit, Lake Frome Basin, South Australia: Mineralogy, geochemistry, and a time-constrained model for its genesis

The sandstone-hosted Beverley uranium deposit is located in terrestrial sediments in the Lake Frome basin in the North Flinders Ranges, South Australia. The deposit is 13 km from the U-rich Mesoproterozoic basement of the Mount Painter inlier, which is being uplifted 100 to 200 m above the basin by neotectonic activity that probably initiated in the early Pliocene. The mineralization was deposited mainly in organic matter-poor Miocene lacustrine sands and partly in the underlying reductive strata comprising organic matter-rich clays and silts. The bulk of the mineralization consists of coffinite and/or uraninite nodules, growing around Co-rich pyrite with an S isotope composition (delta S-34 = 1.0 +/- 0.3 parts per thousand), suggestive of an early diagenetic lacustrine origin. In contrast, authigenic sulfides in the bulk of the sediments have a negative S isotope signature (delta S-34 ranges from -26.2 to -35.5 parts per thousand), indicative of an origin via bacterially mediated sulfate reduction. Minor amounts of Zn-bearing native copper and native lead also support the presence of specific, reducing microenvironments in the ore zone. Small amounts of carnotite are associated with the coffinite ore and also occur beneath a paleosoil horizon overlying the uranium deposit. Provenance studies suggest that the host Miocene sediments were derived from the reworking of Early Cretaceous glacial or glaciolacustrine sediments ultimately derived from Paleozoic terranes in eastern Australia. In contrast, the overlying Pliocene strata were in part derived from the Mesoproterozoic basement inlier. Mass-balance and geochemical data confirm that granites of the Mount Painter domain were the ultimate source of U and BEE at Beverley. U-Pb dating of coffinite and carnotite suggest that the U mineralization is Pliocene (6.7-3.4 Ma). The suitability of the Beverley deposit for efficient mining via in situ leaching, and hence its economic value, are determined by the nature of the hosting sand unit, which provides the permeability and low reactivity required for high fluid flow and low chemical consumption. These favorable sedimentologic and geometrical features result from a complex conjunction of factors, including deposition in lacustrine shore environment, reworking of angular sands of glacial origin, deep Pliocene weathering, and proximity to an active fault exposing extremely U rich rocks.

Emplacement of ultramafic rocks into the continental crust monitored by light and other trace elements: An example from the Geisspfad body (Swiss-Italian Alps)

In order to evaluate the influence of continental crustal rocks on trace element budgets of serpentinized peridotites incorporated into the continental crust, we have analyzed the chemical composition of whole rock samples and minerals of the Geisspfad ultramafic complex (Swiss-Italian Alps). This complex represents a relict oceanic succession composed of serpentinites, ophicarbonates and metabasic rocks, emplaced into crustal gneisses during Alpine collision. Following peak metamorphic amphibolite facies conditions, fluid flow modified some of the trace element contents of ophicarbonates and deformed serpentinites close to the contact with country rocks. The fluid originated from the surrounding continental crustal rocks as documented by the increase of Pb in the serpentinites, and by the strongly negative all) values (-112 parts per thousand) of some ultramafic rocks close to the contact with surrounding gneisses. Little or no modification of the fluid mobile elements Li, B or U was observed in the serpentinite. In-situ analysis of light elements of serpentinite minerals indicate redistribution of light elements coupled to changes of mineral modes towards the outer 100-150 m of the massif. In the centre of the massif, Li is preferentially concentrated in olivine, while Be and B are hosted by tremolite. In contrast, at the outer rim of the massif, Li and Be are preferentially incorporated into diopside, and B into antigorite. This redistribution of light elements among the different minerals is visible in the serpentinite, at a maximum distance of -100-150 m from the ophicarbonate-metabasite contact. Our results show that interaction of ultramafic rocks and crust-derived fluids can be easily detected by studies of Pb and partial derivative D in whole rocks. We argue that small ultramafic bodies potentially record an emplacement-related trace element signature, and that crustal light element values in ultramafic rocks are not necessarily derived from a subducting slab. (C) 2008 Elsevier B.V. All rights reserved.