Structural and strain analysis in the Late Paleozoic coastal accretionary wedge of central Chile

Abstract In this study structural and finite strain data are used to explore the tectonic evolution and the exhumation history of the Chilean accretionary wedge. The Chilean accretionary wedge is part of a Late Paleozoic subduction complex that developed during subduction of the Pacific plate underneath South America. The wedge is commonly subdivided into a structurally lower Western Series and an upper Eastern Series. This study shows the progressive development of structures and finite strain from the least deformed rocks in the eastern part of the Eastern Series of the accretionary wedge to higher grade schist of the Western Series at the Pacific coast. Furthermore, this study reports finite-strain data to quantify the contribution of vertical ductile shortening to exhumation. Vertical ductile shortening is, together with erosion and normal faulting, a process that can aid the exhumation of high-pressure rocks. In the east, structures are characterized by upright chevron folds of sedimentary layering which are associated with a penetrative axial-plane foliation, S1. As the F1 folds became slightly overturned to the west, S1 was folded about recumbent open F2 folds and an S2 axial-plane foliation developed. Near the contact between the Western and Eastern Series S2 represents a prominent subhorizontal transposition foliation. Towards the structural deepest units in the west the transposition foliation became progressively flat lying. Finite-strain data as obtained by Rf/Phi and PDS analysis in metagreywacke and X-ray texture goniometry in phyllosilicate-rich rocks show a smooth and gradual increase in strain magnitude from east to west. There are no evidences for normal faulting or significant structural breaks across the contact of Eastern and Western Series. The progressive structural and strain evolution between both series can be interpreted to reflect a continuous change in the mode of accretion in the subduction wedge. Before ~320-290 Ma the rocks of the Eastern Series were frontally accreted to the Andean margin. Frontal accretion caused horizontal shortening and upright folds and axial-plane foliations developed. At ~320-290 Ma the mode of accretion changed and the rocks of the Western Series were underplated below the Andean margin. This basal accretion caused a major change in the flow field within the wedge and gave rise to vertical shortening and the development of the penetrative subhorizontal transposition foliation. To estimate the amount that vertical ductile shortening contributed to the exhumation of both units finite strain is measured. The tensor average of absolute finite strain yield Sx=1.24, Sy=0.82 and Sz=0.57 implying an average vertical shortening of ca. 43%, which was compensated by volume loss. The finite strain data of the PDS measurements allow to calculate an average volume loss of 41%. A mass balance approximates that most of the solved material stays in the wedge and is precipitated in quartz veins. The average of relative finite strain is Sx=1.65, Sy=0.89 and Sz=0.59 indicating greater vertical shortening in the structurally deeper units. A simple model which integrates velocity gradients along a vertical flow path with a steady-state wedge is used to estimate the contribution of deformation to ductile thinning of the overburden during exhumation. The results show that vertical ductile shortening contributed 15-20% to exhumation. As no large-scale normal faults have been mapped the remaining 80-85% of exhumation must be due to erosion.

UHP metamorphic rocks of the Eastern Rhodope Massif, NE Greece: new constraints from petrology, geochemistry and zircon ages

In this PhD thesis, a multidisciplinary study has been carried out on metagranitoids and paragneisses from the Eastern Rhodope Massif, northern Greece, to decipher the pre-Alpine magmatic and geodynamic evolution of the Rhodope Massif and to correlate the eastern part with the western/central parts of the orogen. The Rhodope Massif, which occupies the major part of NE Greece and S Bulgaria, represents the easternmost part of the Internal Hellenides. It is regarded as a nappe stack of high-grade units, which is classically subdivided into an upper unit and a lower unit, separated by a SSE-NNW trending thrust plane, the Nestos thrust. Recent research in the central Greek Rhodope Massif revealed that the two units correspond to two distinct terranes of different age, the Permo-Carboniferous Thracia Terrane, which was overthrusted by the Late Jurassic/Early Cretaceous Rhodope Terrane. These terranes are separated by the Nestos suture, a composite zone comprising metapelites, metabasites, metagranitoids and marbles, which record high-pressure and even ultrahigh-pressure metamorphism in places. Similar characteristic rock associations were investigated during this study along several well-constrained cross sections in vincity to the Ada, Sidiro and Kimi villages in the Greek Eastern Rhodope Massif. Field evidence revealed that the contact zone of the two terranes in the Eastern Rhodope Massif is characterized by a mélange of metapelites, migmatitic amphibolites/eclogites, strongly sheared orthogneisses and marbles. The systematical occurrence of this characteristic rock association between the terranes implies that the Nestos suture is a continuous belt throughout the Greek Rhodope Massif. In this study, a new UHP locality could be established and for the first time in the Greek Rhodope, metamorphic microdiamonds were identified in situ in their host zircons using Laser-Raman spectroscopy. The presence of the diamonds as well as element distribution patterns of the zircons, obtained by TOF-SIMS, indicate metamorphic conditions of T > 1000 °C and P > 4 GPa. The high-pressure and ultrahigh-pressure rocks of the mélange zone are considered to have formed during the subduction of the Nestos Ocean in Jurassic times at ~150 Ma. Melting of metapelitic rocks at UHP conditions facilitated the exhumation to lower crustal levels. To identify major crust forming events, basement granitoids were dated by LA-SF-ICPMS and SHRIMP-II U-Pb analyses of zircons. The geochronological results revealed that the Eastern Rhodope Massif consists of two crustal units, a structurally lower Permo-Carboniferous unit corresponding to the Thracia Terrane and a structurally upper Late Jurassic/Early Cretaceous unit corresponding to the Rhodope Terrane, like it was documented for the Central Rhodope Massif. Inherited zircons in the orthogneisses from the Thracia Terrane of the Eastern Rhodope Massif indicate the presence of a pre-existing Neoproterozoic and Ordovician-Silurian basement in this region. Triassic magmatism is witnessed by the zircons of few orthogneisses from the easternmost Rhodope Massif and is interpreted to be related to rifting processes. Whole-rock major and trace element analyses indicate that the metagranitoids from both terranes originated in a subduction-related magmatic-arc environment. The Sr-Nd isotope data for both terranes of the Eastern and Central Rhodope Massif suggest a mixed crust-mantle source with variable contributions of older crustal material as already indicated by the presence of inherited zircons. Geochemical and isotopic similarity of the basement of the Thracia Terrane and the Pelagonian Zone implies that the Thracia Terrane is a fragment of a formerly unique Permo-Carboniferous basement, separated by rifting and opening of the Meliata-Maliac ocean system in Triassic times. A branch of the Meliata-Maliac ocean system, the Nestos Ocean, subducted northwards in Late Jurassic times leading to the formation of the Late Jurassic/Early Cretaceous Rhodope magmatic arc on remnants of the Thracia Terrane as suggested by inherited Permo-Carboniferous zircons. The ~150 Ma zircon ages of the orthogneisses from the Rhodope Terrane indicate that subduction-related magmatism and HP/UHP metamorphism occurred during the same subduction phase. Subduction ceased due to the closure of the Nestos Ocean in the Late Jurassic/Early Cretaceous. The post-Jurassic evolution of the Rhodope Massif is characterized by the exhumation of the Rhodope core complex in the course of extensional tectonics associated with late granite intrusions in Eocene to Miocene times.

Characterization of fold-related fractures in the carbonate rocks of the Cingoli anticline, northern Apennines, Italy

Thrust fault-related folds in carbonate rocks are characterized by deformation accommodated by different structures, such as joints, faults, pressure solution seams, and deformation bands. Defining the development of fracture systems related to the folding process is significant both for theoretical and practical purposes. Fracture systems are useful constrains in order to understand the kinematical evolution of the fold. Furthermore, understanding the relationships between folding and fracturing provides a noteworthy contribution for reconstructing the geodynamic and the structural evolution of the studied area. Moreover, as fold-related fractures influence fluid flow through rocks, fracture systems are relevant for energy production (geothermal studies, methane and CO2 , storage and hydrocarbon exploration), environmental and social issues (pollutant distribution, aquifer characterization). The PhD project shows results of a study carried out in a multilayer carbonate anticline characterized by different mechanical properties. The aim of this study is to understand the factors which influence the fracture formation and to define their temporal sequence during the folding process. The studied are is located in the Cingoli anticline (Northern Apennines), which is characterized by a pelagic multilayer characterized by sequences with different mechanical stratigraphies. A multi-scale analysis has been made in several outcrops located in different structural positions. This project shows that the conceptual sketches proposed in literature and the strain distribution models outline well the geometrical orientation of most of the set of fractures observed in the Cingoli anticline. On the other hand, the present work suggests the relevance of the mechanical stratigraphy in particular controlling the type of fractures formed (e.g. pressure solution seams, joints or shear fractures) and their subsequent evolution. Through a multi-scale analysis, and on the basis of the temporal relationship between fracture sets and their orientation respect layering, I also suggest a conceptual model for fracture systems formation.

Dynamic development of hydrofractures

Natürliche hydraulische Bruchbildung ist in allen Bereichen der Erdkruste ein wichtiger und stark verbreiteter Prozess. Sie beeinflusst die effektive Permeabilität und Fluidtransport auf mehreren Größenordnungen, indem sie hydraulische Konnektivität bewirkt. Der Prozess der Bruchbildung ist sowohl sehr dynamisch als auch hoch komplex. Die Dynamik stammt von der starken Wechselwirkung tektonischer und hydraulischer Prozesse, während sich die Komplexität aus der potentiellen Abhängigkeit der poroelastischen Eigenschaften von Fluiddruck und Bruchbildung ergibt. Die Bildung hydraulischer Brüche besteht aus drei Phasen: 1) Nukleation, 2) zeitabhängiges quasi-statisches Wachstum so lange der Fluiddruck die Zugfestigkeit des Gesteins übersteigt, und 3) in heterogenen Gesteinen der Einfluss von Lagen unterschiedlicher mechanischer oder sedimentärer Eigenschaften auf die Bruchausbreitung. Auch die mechanische Heterogenität, die durch präexistierende Brüche und Gesteinsdeformation erzeugt wird, hat großen Einfluß auf den Wachstumsverlauf. Die Richtung der Bruchausbreitung wird entweder durch die Verbindung von Diskontinuitäten mit geringer Zugfestigkeit im Bereich vor der Bruchfront bestimmt, oder die Bruchausbreitung kann enden, wenn der Bruch auf Diskontinuitäten mit hoher Festigkeit trifft. Durch diese Wechselwirkungen entsteht ein Kluftnetzwerk mit komplexer Geometrie, das die lokale Deformationsgeschichte und die Dynamik der unterliegenden physikalischen Prozesse reflektiert. rnrnNatürliche hydraulische Bruchbildung hat wesentliche Implikationen für akademische und kommerzielle Fragestellungen in verschiedenen Feldern der Geowissenschaften. Seit den 50er Jahren wird hydraulisches Fracturing eingesetzt, um die Permeabilität von Gas und Öllagerstätten zu erhöhen. Geländebeobachtungen, Isotopenstudien, Laborexperimente und numerische Analysen bestätigen die entscheidende Rolle des Fluiddruckgefälles in Verbindung mit poroelastischen Effekten für den lokalen Spannungszustand und für die Bedingungen, unter denen sich hydraulische Brüche bilden und ausbreiten. Die meisten numerischen hydromechanischen Modelle nehmen für die Kopplung zwischen Fluid und propagierenden Brüchen vordefinierte Bruchgeometrien mit konstantem Fluiddruck an, um das Problem rechnerisch eingrenzen zu können. Da natürliche Gesteine kaum so einfach strukturiert sind, sind diese Modelle generell nicht sonderlich effektiv in der Analyse dieses komplexen Prozesses. Insbesondere unterschätzen sie die Rückkopplung von poroelastischen Effekten und gekoppelte Fluid-Festgestein Prozesse, d.h. die Entwicklung des Porendrucks in Abhängigkeit vom Gesteinsversagen und umgekehrt.rnrnIn dieser Arbeit wird ein zweidimensionales gekoppeltes poro-elasto-plastisches Computer-Model für die qualitative und zum Teil auch quantitativ Analyse der Rolle lokalisierter oder homogen verteilter Fluiddrücke auf die dynamische Ausbreitung von hydraulischen Brüchen und die zeitgleiche Evolution der effektiven Permeabilität entwickelt. Das Programm ist rechnerisch effizient, indem es die Fluiddynamik mittels einer Druckdiffusions-Gleichung nach Darcy ohne redundante Komponenten beschreibt. Es berücksichtigt auch die Biot-Kompressibilität poröser Gesteine, die implementiert wurde um die Kontrollparameter in der Mechanik hydraulischer Bruchbildung in verschiedenen geologischen Szenarien mit homogenen und heterogenen Sedimentären Abfolgen zu bestimmen. Als Resultat ergibt sich, dass der Fluiddruck-Gradient in geschlossenen Systemen lokal zu Störungen des homogenen Spannungsfeldes führen. Abhängig von den Randbedingungen können sich diese Störungen eine Neuausrichtung der Bruchausbreitung zur Folge haben kann. Durch den Effekt auf den lokalen Spannungszustand können hohe Druckgradienten auch schichtparallele Bruchbildung oder Schlupf in nicht-entwässerten heterogenen Medien erzeugen. Ein Beispiel von besonderer Bedeutung ist die Evolution von Akkretionskeilen, wo die große Dynamik der tektonischen Aktivität zusammen mit extremen Porendrücken lokal starke Störungen des Spannungsfeldes erzeugt, die eine hoch-komplexe strukturelle Entwicklung inklusive vertikaler und horizontaler hydraulischer Bruch-Netzwerke bewirkt. Die Transport-Eigenschaften der Gesteine werden stark durch die Dynamik in der Entwicklung lokaler Permeabilitäten durch Dehnungsbrüche und Störungen bestimmt. Möglicherweise besteht ein enger Zusammenhang zwischen der Bildung von Grabenstrukturen und großmaßstäblicher Fluid-Migration. rnrnDie Konsistenz zwischen den Resultaten der Simulationen und vorhergehender experimenteller Untersuchungen deutet darauf hin, dass das beschriebene numerische Verfahren zur qualitativen Analyse hydraulischer Brüche gut geeignet ist. Das Schema hat auch Nachteile wenn es um die quantitative Analyse des Fluidflusses durch induzierte Bruchflächen in deformierten Gesteinen geht. Es empfiehlt sich zudem, das vorgestellte numerische Schema um die Kopplung mit thermo-chemischen Prozessen zu erweitern, um dynamische Probleme im Zusammenhang mit dem Wachstum von Kluftfüllungen in hydraulischen Brüchen zu untersuchen.

Palaeotsunami imprints in the near-coast sedimentary records of the Gulfs of Lakonia and Argolis (Peloponnese, Greece)

From historical accounts it is well-known that the coasts of the Gulfs of Lakonia and Argolis (southern and eastern Peloponnese, Greece) have been repeatedly affected by tsunamis during historical times. It is assumed that these palaeotsunamis left sedimentological and geomorphological traces in the geological record which are still detectable these days. As both gulfs are located within one of the seismically most active regions in whole western Eurasia in particular the nearby Hellenic Trench is regarded as the main trigger for tsunami generation. Against this background, selected near-coast sedimentary archives were studied by means of sedimentological, geomorphological, geophysical, geochemical and microfaunal investigations in order to detect signatures of Holocene palaeotsunamigenic activity. The investigations revealed allochthonous sediment layers featuring distinctive sedimentary characteristics of marine high-energy event deposits in most of the investigated study areas. In order to differentiate between the geomorphodynamic driving mechanisms for the deposition of the associated marine high-energy event layers, a multi-method approach was used. The detected high-energy marine deposits are suggested to be of tsunamigenic origin. Radiocarbon dating results allowed establishing local event geo-chronostratigraphies and correlations on a local and regional scale as well as correlations with already described palaeotsunami findings on a supra-regional scale. The geochronological dataset attests repeated tsunamigenic activity at least since the 5th millennium BC up to the 17th century AD. For the studied areas in southeastern Lakonia up to four palaeotsunami event generations were identified, for central Lakonia three and for the investigated areas around the Argolis Gulf also up to four. Comparing the findings with literature data, chronological correlations were found with palaeotsunami deposits detected in near-coast geological archives of Akarnania, of the southwestern, the western and northwestern Peloponnese, with event deposits found on Crete and on the Ionian Islands of Cefalonia and Lefkada as well as with findings from southeastern Sicily (Italy) and Cesarea (Israel). By the identification of multiple palaeotsunami event layers, disturbing autochthonous near-coast sedimentary records of the Gulfs of Lakonia and Argolis during the last seven millennia, a significant tsunami frequency is attested for these regions.

Provenance signature of a forearc basin modified by spreading ridge subduction: Detrital zircon geochronology and detrital modes from the Paleogene Arkose Ridge Formation, southern Alaska

Upper Paleocene–Eocene boulder conglomerate, cross-stratified sandstone, and laminated carbonaceous mudstone of the Arkose Ridge Formation exposed in the southern Talkeetna Mountains record fluvial-lacustrine deposition proximal to the volcanic arc in a forearc basin modified by Paleogene spreading ridge subduction beneath southern Alaska. U-Pb ages of detrital zircon grains and modal analyses were obtained from stratigraphic sections spanning the 2,000 m thick Arkose Ridge Formation in order to constrain the lithology, age, and location of sediment sources that provided detritus. Detrital modes from 24 conglomerate beds and 54 sandstone thin sections aredominated by plutonic and volcanic clasts and plagioclase feldspar with minor quartz, schist, hornblende, argillite, and metabasalt. Westernmost sandstone and conglomerate strata contain <5% volcanic clasts whereas easternmost sandstone and conglomerate strata contain 40 to >80% volcanic clasts. Temporally, eastern sandstones andconglomerates exhibit an upsection increase in volcanic detritus from <40 to >80% volcanic clasts. U-Pb ages from >1400 detrital zircons in 15 sandstone samples reveal three main populations: late Paleocene–Eocene (60-48 Ma; 16% of all grains), Late Cretaceous–early Paleocene (85–60 Ma; 62%) and Jurassic–Early Cretaceous (200–100 Ma; 12%). A plot of U/Th vs U-Pb ages shows that >97% of zircons are <200 Ma and>99% of zircons have <10 U/Th ratios, consistent with mainly igneous source terranes. Strata show increased enrichment in late Paleocene–Eocene detrital zircons from <2% in the west to >25% in the east. In eastern sections, this younger age population increases temporally from 0% in the lower 50 m of the section to >40% in samples collected >740 m above the base. Integration of the compositional and detrital geochronologic data suggests: (1) Detritus was eroded mainly from igneous sources exposed directly north of the Arkose Ridge Formation strata, mainly Jurassic–Paleocene plutons and Paleocene–Eocenevolcanic centers. Subordinate metamorphic detritus was eroded from western Mesozoic low-grade metamorphic sources. Subordinate sedimentary detritus was eroded from eastern Mesozoic sedimentary sources. (2) Eastern deposystems received higher proportions of juvenile volcanic detritus through time, consistent with construction of adjacent slab-window volcanic centers during Arkose Ridge Formation deposition. (3)Western deposystems transported detritus from Jurassic–Paleocene arc plutons that flank the northwestern basin margin. (4) Metasedimentary strata of the Chugach accretionaryprism, exposed 20-50 km south of the Arkose Ridge Formation, did not contribute abundant detritus. Conventional provenance models predict reduced input of volcanic detritus to forearc basins during exhumation of the volcanic edifice and increasing exposure ofsubvolcanic plutons (Dickinson, 1995; Ingersoll and Eastmond, 2007). In the forearc strata of these conventional models, sandstone modal analyses record progressive increases upsection in quartz and feldspar concomitant with decreases in lithic grains, mainly volcanic lithics. Additionally, as the arc massif denudes through time, theyoungest detrital U-Pb zircon age populations become significantly older than the age of forearc deposition as the arc migrates inboard or ceases magmatism. Westernmost strata of the Arkose Ridge Formation are consistent with this conventional model. However, easternmost strata of the Arkose Ridge Formation contain sandstone modes that record an upsection increase in lithic grains accompanied by a decrease in quartz and feldspar, and detrital zircon age populations that closely match the age of deposition. This deviation from the conventional model is due to the proximity of the easternmost strata to adjacent juvenile volcanic rocks emplaced by slab-window volcanic processes. Provenance data from the Arkose Ridge Formation show that forearc basins modified by spreading ridge subduction may record upsection increases in non-arc, syndepositional volcanic detritusdue to contemporaneous accumulation of thick volcanic sequences at slab-window volcanic centers. This change may occur locally at the same time that other regions of the forearc continue to receive increasing amounts of plutonic detritus as the remnant arc denudes, resulting in complex lateral variations in forearc basin petrofacies and chronofacies.