Search

Identifying potential petroleum traps in petroleum basins is one of the key challenges in petroleum exploration. Specifically, it is the identification of probable petroleum traps within a set of stratigraphic traps of a particular location of source rock and carrier bed. One solution lies in understanding the behaviour of hydrocarbon flow during secondary migration, and the evaluation of the probability of successful transport from the source rock to the trap. Modern reservoir simulators rely on numerical methods to model the oil/gas secondary migration. Using numerical simulators is, however, cumbersome and requires high volumes of data and computation time, which affects successful decision-making in exploration planning. Yet, analytical models are fast and allow for multivariant analysis of hydrocarbon secondary migration requiring only a moderate amount of geological data. This study presents the analytical modelling of hydrocarbon buoyant transport in petroleum basins by including the (i) areal variation of stringers’ cross-section, (ii) chemical reactions including oil biodegradation and (iii) hydrological water flow. The explicit formula is provided for the first and last moments of hydrocarbon arrival at the trap, describing the dynamics of filling of the trap. Field data from Australian and Chinese basins are used to investigate the effects of the above-mentioned parameters on the first and last moments of hydrocarbon arrival at the trap.

The Gippsland basin is part of the Australian southern margin rift system. It is a world class oil and gas producing province located about 200 km east of the city of Melbourne, and covers about 46 000 km2 onshore and offshore. The offshore part is a post orogenic continental margin basin formed during Jurassic-cretaceous resulting from the breakup of Gondwana supercontinent in the Mesozoic and the separation of Antarctica and Australia. A second rifting phase occurred with a NE-SW associated with the development of the Tasman Sea. Gippsland basin is filled by three major lithostratigraphic groups, namely: the Strzelecki group, Latrobe and Seaspray groups. The sedimentary fill unconformably overlies a Paleozoic basement made up of igneous and folded sedimentary rocks of the Lachlan orogenic. The objective of this study is to help constraining the tectonostratigraphic evolution and the structural evolution model of the basin based on 2D seismic interpretation as reflection seismic data. The interpretation of seismic reflection data is a fundamental method for determining the geometry and displacement of faults in the subsurface which is primordial in studying structural events in sedimentary basins.

The Atlas system is an intracontinental chain established upon a Paleozoic substratum by the inversion of Triassic basins starting in the late Mesozoic-early Cenozoic. The inversion of the chain is related to the Atlas rift system that was influenced by the opening of the Central Atlantic in the West and the Tethys in the north. This was coeval with a regional exhumation following the Alpine shortening responsible for the continuous uplifting of the chain since Late Cretaceous. The structural history and chronology of events are still matter of debates. To contribute to this, we focus on the Western High Atlas (WHA) aiming a retro enactment of the paleo-stresses states, by analyzing deformation structures at various scales. The geological data were collected at different stratigraphic levels: from the contact Paleozoic basement/Mesozoic cover interface to the Triassic detrital formations of the Argana corridor in the east, to the Jurassic-Cretaceous and Cenozoic plateaus in the west. Preliminary results highlight two major tectonic events: (1)- a first extensive event, with sub-horizontal minimal principal stress σ3 oriented NW-SE, that is linked to the Central Atlantic basin opening. This event is characterized by pull apart basins structured into horsts and grabens. (2)- a second compressive event, marked by NE-SW to NNE-SSW shortening. The later is subdivided into two episodes: i- an early post-rift episode (Middle-Late Jurassic to Early Cretaceous), marked by stylolites and meso-structures that occurred at the beginning of the main uplifting stage. ii- a late compression episode, characterized by a maximum principal stress σ1 mainly oriented NNE-SSW to NNW-SSE, starting at the late Cretaceous and accelerating during the Tertiary, simultaneously with the Africa-Europe collision.Keywords: Paleo-stress, Structural analysis, Atlas rift system, Tectonic inversion, Western High Atlas Morocco, Alpine orogeny.

The Wilson Cycle (ca. 250-300 Ma) is driven by plate tectonics, and is expressed by the opening of oceanic basins, which are subsequently closed to form orogens. Morocco has experienced several orogenic cycles: e.g., Eburnean, Pan-African, Hercynian (Variscan) and Alpine (Atlasian). The last Wilson Cycle that shaped the geology of Morocco is the Alpine Cycle (300-0 Ma). It started by the dismantlement of the Pangea Supercontinent leading to the opening of the Atlantic Ocean and the Atlas rift system. The latter was subsequently aborted and inverted into a fold-and-thrust belt during the Cenozoic, due to the opening of the South Atlantic and the convergence between Africa and Iberia-Europe. The Atlas system of NW Africa includes two intracontinental ranges, the High and Middle Atlas extending ENE-WSW and NE-SW, respectively. It is a key piece in the Pangea breakup puzzle, as its evolution captures the kinematic of the African plate during Mesozoic to Cenozoic times. Several studies have examined the structural and stratigraphic architectures of the Atlas system during the rifting phase, by removing the Alpine inversion along palinspaticly reconstructed 2D sections. However, little was done to investigate the crustal structure and the amount of crustal stretching during the Triassic to Jurassic extension. This work aims to reconstruct the pre-alpine architecture of the Atlas system and link it directly with the evolution of the crust. Our goal is to provide insights into the amount of crustal thinning that took place during the Mesozoic rifting in the Atlas domain and constrain the strain distribution within the African plate during the dismantlement of the Pangea.Key words: Wilson cycle; Atlas system; Pangea; Break up; Inversion; Pre-Alpine cycle.

The intracontinental High Atlas Mountains are the result of multiple tectonic events, from Late Permian-Early Mesozoic Pangea breakup to Cenozoic Africa-Eurasia convergence. Prior to the Cenozoic inversion event, the South Atlas Fault Zone (SAFZ), which surrounded the High Atlas from the south, witnessed for the growth of multiple basins during the Liassic rifting in the eastern segment of the belt. The present work aims to shed light on the development of the Beni Bassia Liassic basin through two faults components of the SAFZ, the E-W J. Amalou faulted anticline (ridge) and Garn Talou fault from the South. Preliminary investigations based on the combination of field data and satellite images reveal significant results regarding the evolution of the basin. In the Beni Bassia Basin, the Ouchbis-Pleinsbachian Formation is an alternating layer of conglomerate limestones and marls, containing olistoliths over 10 meters in length, resulting from the destruction of the underlying platform of the Idikel-Sinemurian Formation. This huge breccia developed at the foot of steep faults on the southern margin of the basin, during the acceleration of the Liassic rifting. During the Cenozoic deformation and uplift of the Eastern High Atlas, bedding directions of Jurassic formations denote progressive changes from E-W in the northwestern to approximately NW-SE to N-S in the southeastern part of the basin. Similarly, the dip direction changes from dipping towards the south to dipping southwest to westward, indicating a dextral movement due to NNW-SSE regional shortening driven by the plates convergence. Other smaller-scale structures confirm dextral transpression kinematics at several places east of the High Atlas, including easterly detachments recorded in Lower Jurassic carbonates. Further investigations are planned to quantify the stress, strain and mechanisms involved in the configuration of the Liassic basin and the subsequent Cenozoic transpression.Keywords: Eastern High Atlas, South Atlas Fault Zone, Beni Bassia Basin, Jurassic, Inversion, Transpression.

Processes such as oblique mid ocean ridge spreading, glacial isostatic adjustment and slope instability provide a highly complex spatial and temporal record of stress in the Fram Strait. The Vestnesa Ridge is a contourite drift bounded by two slow spreading mid ocean ridges located beside a formerly glaciated margin. The total state of stress is difficult to separate into individual components therefore our focus is to ascertain whether there is a stress transfer from the deep crust into the shallow overlying (~200m) sedimentary cover. We use high-resolution P-cable 3D seismic volumes together with 2D seismic, to map deeper faults connecting with near surface deformation. We perform high resolution mapping of the ridge by examining the dip and strike of each distinct fault system. We use a pre trained 3D model to predict faults within each 3D volume and automatically extract faults at multiple intervals to capture temporal stress changes. To minimize noise, the model identifies faults based on edge preserved smoothing for a selection of peak frequencies. In our results we observe fault linkage between parallel faults that may become favourable locations for transtensional and transpressional stress expected in the strike slip regime predicted in the west of the ridge. Our results show that the east of the ridge has a dominant NW-SE fault strike and a present day tensile stress regime while towards the west, the NW-SE assemblage becomes less prominent and multiple fault systems dominate increasing the complexity of the system. We present a high detail comprehensive structural analysis of 3 study sites across the shallow ridge sediments and use our results to investigate differences in the strike and dip between sites to explore the influence of sedimentary faults and ridge geomorphology on the spatial evolution of seafloor seepage at a deep Arctic oceanic basin.

The geomorphological and tectonic processes being responsible for the control of the mountain ranges geometry is feasible to be assessed through geomorphological and morphotectonic indices. Those indices are dependent on tectonic or erosional processes and other factors influencing the development of the landforms. Tectonic geomorphology applications using satellite-based remotely sensed data, such as Digital Elevation Models (DEMs), can highlight specific geomorphic features capable to provide useful information and knowledge towards the evaluation of the regional tectonic activity in mountain ranges. This study examines the Moroccan High Atlas mountain range by using morphotectonic and geomorphometric indices such as the channel steepness index, amplitude of relief index, stream length gradient index, swath profiles, local relief and hillslope mapping to determine the distribution of the tectonic activity variations. Through those indices the evaluation of the geomorphic responses to tectonics takes place by highlighting the relationships between tectonic activity, rock resistance, stream channel slope, active or recent vertical displacements. The outcomes of the geomorphometric and morphotectonic investigation highlight the presence of considerable geomorphic variations across the main fault zones featuring the orogen’s anatomy, while tectonic activity seems to be a major factor controlling and shaping the Moroccan High Atlas mountain range landscape. The methodological framework of this study could be developed into a low-cost technique for assessing seismic hazard, offering a valuable tool towards assessing disaster risk reduction activities, whereas in conjunction with other factors the georesources exploration. AcknowledgementsThe corresponding author acknowledges the 'EXCELSIOR': ERATOSTHENES: Excellence Research Centre for Earth Surveillance and Space-Based Monitoring of the Environment H2020 Widespread Teaming project (www.excelsior2020.eu). The 'EXCELSIOR' project has received funding from the European Union's Horizon 2020 research and innovation programme under Grant Agreement No 857510, from the Government of the Republic of Cyprus through the Directorate General for the European Programmes, Coordination and Development and the Cyprus University of Technology.

Located in the western segment of the intracontinental Atlas system, the Moroccan Central High Atlas is a NE-SW to ENE-WSW-trending Fold-and-Thrust Belt that is formed during the Cenozoic Alpine orogeny by a positive inversion of Triassic-Jurassic basin. It is structurally distinguished from the other segments of the Moroccan High Atlas orogenic belt by the occurrence of S-shaped ENE-WSW oriented tight anticlinal ridges bounding wider synclines. The elongated ridges core disordered association of plutonic rocks, Liassic carbonate and Late Triassic arigilites, whilst the wider synclines are filled by thick Jurassic series with minor magmatic manifestations expressed by mafic and felsic dikes. The origin of these structures has been ascribed to pre-inversion wrench tectonics with significant compressive component whereas they have been attached to post-rift rift block tilting and or salt tectonics in an alternative view. Characterizing the paleostress history is thereby a crucial matter to unravel the structural evolution of these structures. In order to bring new insights into the actual understanding of the Central High Atlas post-rift structural history, we reconstruct the paleostress tensors preserved in the folded Jurassic series of Anemzi and Tirrhist regions based on brittle deformation structures together with calcite twins stress inversion. The preliminary results highlight the presence of pre-folding layer parallel maximum horizontal stress during three stages: E-W to ENE-WSW, NNE-SSW and NW-SE compressions. Local extensional stress features are observed essentially near diapiric structures and the exhumed magmatic intrusions. The latest structural stage is featured by a post-folding NW-SR compression likely related to the recent phases of the Alpine orogeny.

Hydraulic fracturing is increasingly becoming utilized within Underground hard rock mines such as block caves and Sub Level Caves as a way to promote controlled cave propagation, increase resource recovery, and seismic hazard through the manipulation of rock mass properties through fracture surface creation and limiting of stress concentration. While hydraulic fracturing is not a new application, it is still in its infancy in mining projects such as cave mining. It is used on entirely different scales, under different stress regimes and with varying motives. Therefore, more research needs to be carried out in understanding the fundamentals of fracture growths that can help improve hydraulic fracturing applied in mining projects. Predicting breakdown pressure is an important part of the designing of hydraulic fracturing with accurate prediction being the baseline of designing and implementing a successful preconditioning campaign in all industries, but especially so in block cave mining. The most commonly used breakdown pressure theoretical model is the conventional breakdown model and is based on tensile strength and confining stresses acting upon the borehole. This might be imprecise within hard rock mining environments and increasingly so at depth of higher stresses. This work compares indirect tensile strength results and their fracture toughness, from both conventional Brazilian Disc Testing and the recently developed Adelaide University Snapback Indirect Tensile Testing (AUSBIT).By using lateral strain control to stabilise the brittle material responses, AUSBIT allows for the capture of true post-peak behaviour, i.e. controlled fracture propagation can be achieved.The captured post-peak behaviour allows practitioners to measure a more reflective tensile strength and fracture toughness from just one testing method, Alongside this are laboratory hydraulic fracture experiements on the same rock unit, which in turn is used to propose a new Hard Rock Breakdown Pressure prediction based on the conventional method which incorporates fracture toughness when checked against the results of the lab fracture experiments, and other lab studies. The results of this work factors in a progressive toughness of a rock at depths, and creates a more accurate predictor of breakdown pressure in underground hard rock mines under varying stress conditions.

The Moroccan Atlas is an intracontinental chain resulted from an aborted rifting during the Mesozoic time, by an uplifting and moderate shortening during the Late Cretaceous-Cenozoic period. Several studies have highlighted the role of tectonic inversion in the evolution of the High Atlas Range, where strike-slip faults are commonly been considered as a main component of the alpine signature within the High Atlas belt. However, more recent works have focused on the geodynamic model of the evolution of the Atlas Range using different approaches. The structural history and chronology of events are still matter of debates. To contribute to the later, a combined meso and microstructural study was conducted in the western part of the chain. It provided an attempt to quantify paleo-stresses from structural analysis of the Permo-Triassic extensional phase to the tectonic reversal phases, acting from Cenozoic to present days.This work highlighted two major tectonic phases: (1) the first represented by an extensive regime, with a sub-horizontal minimal stress σ3 oriented NE-SW and linked to the Central Atlantic occurrence. This stage is characterized by pull apart basins genesis in horst and graben morphology. (2) the second phase represented by a weakly tilted compression with a maximum stress σ1 oriented in set NNE-SSW to NNW-SSE. This compression began in the Tertiary, contemporary with the Africa and Europe collision. the related inversions are printed at the paleozoic basement/mesozoic cover interface from the Eastern area to the Jurassic-Cretaceous and Cenozoic plateaus in the West, passing through the Triassic detrital formations of the Argana corridor.Keywords: Paleo-stress, Structural analysis, Tectonic inversion, Western high Atlas, Morocco, Alpine orogeny.

The Ouled Abdoun sedimentary basin in Morocco contains the largest phosphate reserves in the world. In the southeastern parts of the basin, the phosphorite deposits lay from the Maastrichtian to the Lutetian sediments of the Tadla Plain. This section has a thickness of ~ 30 m and generally protected from erosion by a relatively strong Turritella slab cover. The phosphorite deposits are distributed in horizontal strata interbedded with levels of limestone, marl and clay, that present various silicifications from the Ypresian. This work aims to study and determine their petrographic, granulometric and mineralogical compositions. A multidisciplinary approach was adopted to achieve these objectives. First, the use of sedimentology and the application of sequence stratigraphy allowed the definition of three depositional sequences in this deposit. Second, the granulometric analysis of the phosphate facies reflects a dominance of well classified medium grains. Furthermore, the analysis of the Visher curves revealed up to three major modes of transport: traction, saltation and suspension. Based on their mineral composition, the microfacies are classified into two phosphate families (or types): Coprolite Intraphospharenite type and Granular Pelphosphalrenite type. Finally, the mineral parageneses recognized by the XRD analyses revealed that phosphorits consist mainly of carbonate, silica and apatitic phases in the section of Tadla.Keywords: Phosphorite deposit, Tadla plain, Maastrichtian-Lutetian.

The High Atlas intracontinental belt is an inverted rift formed along the Northwest African Craton, between the Central Atlantic Ocean and the Africa-Eurasia plates’ boundaries. Due to its particular location at this triple junction, the Moroccan High Atlas is a good example to investigate the tectonic and geodynamic complexity within continental interiors. The present paper aims to give an overview of the actual understanding of the post-Variscan tectonic and structural history of the High Atlas belt and to shed light on some of the most actively debatable questions relating to the geodynamic context of its uplifts. The High Atlas basin Alpine history started with a period of passive rifting during the Middle-Late Triassic-Early Liassic resulting in a complex rift system actuated by the reactivation of Variscan NNE to ENE-trending structures in the framework of a global NW-SE extension. The Late Triassic rift basins were filled by thick salt-bearing red-beds whose post-rift extensive halokinetics has a considerable impact on the structural architecture of the Western and Central High Atlas. Extensional deformation was renewed in the Central High Atlas at the Late Liassic-Early Bajocian associated with the drowning of Liassic carbonate platform coevally with salt mobilization which led to the initiation of NE to ENE-elongated narrow diapiric ridges. By the Middle Jurassic, a major part of the High Atlas basin was submitted to a widespread post-rift epeiorogenic upward motion that lasted up until the Lower Cretaceous, and to which is associated the occurrence of transitional to moderately alkaline magmatism that was emplaced within the established ridges across basement faults. The driving geodynamic mechanisms and the tectonic setting of this uplift are poorly known but the large extent of the exhumed lands in both margins of the Central Atlantic and the synchronicity between the upward movements and the drastic increase of the Central Atlantic Ocean spreading rate suggests a possible causal link between the Central Atlantic dynamics and the Middle Jurassic-Lower Cretaceous post-rift exhumation. The Late Cretaceous signals the onset of the convergence between Africa and Eurasia and the consequent Alpine orogeny. The effective inversion of the Atlas paleo-rift begun during the Eocene and accentuated during the Neogene in an intra-continental environment. The western segment of range witness a basement-involved thick-skinned faulting and folding style, whereas thin-skinned deformation detached mostly on Upper Triassic-Early Liassic evaporitic layers is focused in external forelands. In the Central High Atlas, the salt-involved Cenozoic deformation triggered the evaporites squeezing and extrusion and the inversion of the pre-structured magmatic cored ridges into compressional anticlines. During the Alpine orogeny, the High Atlas orogen underwent weak crustal shortening with the development of a moderate crustal root that resides over a highly elevated Lithosphere-Asthenosphere Boundary. The most prevailing view predicts the occurrence of a hot mantle anomaly, superimposed to the collisional deformation, beneath the Atlas range that maintains the high topography featuring the High Atlas orography. Although advanced studies have been performed, several questions regarding the kinematic history, the origin of the vertical movements, the role played by the week salt layers during both extensional and compressional deformation stages, and the factors controlling the selective inversion of the paleo-rift normal faults and the deep geometry of the range’s major faults, are still not well resolved and deserve to be deeply re-investigated.

The Moroccan High Atlas mountain range is an aborted Mesozoic rift basin that was moderately shortened during the Late Cretaceous‒Cenozoic inversion. The range is currently featured in its central part by the presence of conspicuous S-shaped open gentle synclines where Middle Jurassic strata crop out, with sub-horizontal bottom, separated by 15-to-80-km narrow faulted anticline ridges with two distinct directions: ENE and NE. The tight anticline ridges are cored by Triassic continental red-beds intruded by the CAMP basalts and subsequently by Upper Jurrasic‒Lower Cretaceous alkaline magmatism. Regional cleavage with very low-grade anchi- to epi-zonal metamorphism are depicted along several structures of the High Atlas, particularly the NE-trending anticlines. The sedimentary layers thickness, on the other hand, gets thinner towards the faulted anticlines with the development of intraformational truncations. The structural history of the High Atlas syncline-topped anticlinal ridges remains a controversial matter. Any attempt to reconstruct the evolutionary process of such folded structures must take into consideration the following circumstances:After a Triassic rifting episode followed by the establishment of Liassic carbonate platform, the High Atlas basin underwent a wide spread exhumation event at the time interval between the Middle Jurassic and Lower Cretaceous leading to the deposition of continental detrital series and sedimentary hiatus; The upward motion was accompanied with the emplacement of alkaline magmas in the Central High Atlas; A complex halokinetic history characterizes the Central High Atlas salt province during both pre-orogenic and orogenic stages; During the Late Cretaceous‒Cenozoic, the High Atlas experienced a moderate crustal shortening which was focused essentially within the range’s borders; In order to bring new insights to the structural history of the High Atlas folded structures, a structural investigation was carried out in Tirrhist and Anemzi ridges. In each station, fractures measurements were taken, and oriented samples were collected for micro-structural analysis. First paleo-stress inversion in some stations reveals the presence of pre-folding bedding-parallel maximal horizontal stress oriented NE to NNE. For a deep analysis of pre syn and post-folding stresses history, we use a calcite stress inversion technique, namely Etchecopar’s method, to unravel the paleo-stresses orientations and to quantify the differential stresses during the different episodes of deformation. The present work is a preliminary attempt to quantify tectonic stresses in the hinterland of an arguably weakly deformed orogenic belt.

Alkaline complexes are an important target for geological exploration, with both scientific and economic interests. They are host to different types of mineral deposits, such as Rare Earths, igneous phosphates, -and K-rich minerals and rocks. In Morocco, the Central High-Atlas (CHA) hosts several transitional to alkaline complexes ranging from Upper Jurassic to Eocene and showing almost all the differentiation terms of transitional to alkaline suites. These alkaline complexes are however poorly explored and their potential in terms of mineral resources is still elusive.The aim of this research is to use Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) to discriminate different transitional to alkaline rock lithologies and their associated mineralizations. For that purpose, series of band ratios proven to be sensitive to the silica, mafic, felsic and carbonate contents of transitional to alkaline rocks were applied. Our results show that the major Upper Jurassic magmatic intrusions of Moroccan CHA, such as Anemzi, Inouzane, Tassent, and Tasraft, hold distinct igneous facies, mainly composed of Mafic to felsic rocks. Field and petrographic observations have confirmed the ASTER results and highlighted that these rocks are formed of gabbro to syenite. The later are associated with significant feldspar concentrations, but also host apatite, garnet, and magnetite vein-type ores. Thereafter, field- and petrographic-based data were used as training data to perform a supervised classification allowing to refine the geological mapping of the studied alkaline intrusions.

The spatial distribution of permeable fracture networks is constrained in the subsurface Cooper Basin (Australia) through the integration of most positive curvature ( K + ) analysis of five three dimensional (3D) seismic surveys, wellbore fracture data, and geomechanical modelling. The K + provides the likely distribution of subsurface extensional fractures based on the stress redistribution along the outer arc of an anticline. These results are reprocessed into the semblance (similarity) attribute to improve the signal-noise ratio prior to being extracted along the gas-rich and low permeability Patchawarra Formation. The subsurface fracture distribution maps show the geometry and density of extensional fractures that have been shown to facilitate fluid migration in this province, particularly those striking SE-NW and E-W, as these are properly oriented to undergo tensile reactivation. The density of these permeable SE-NW and E-W striking fracture sets is predicted to significantly increase along E-W elongate anticlines and may contribute to improved hydrocarbon recovery in this low permeability stratigraphic interval. We show that wellbore fracture data (n = 917) in the Patchawarra Formation presents a close relationship to K + results and can be used to generate high density structural data in this basin, with the methodology applicable to other subsurface and offshore provinces.

Selecting a seismic time-to-depth conversion method can be a subjective choice that is made by geophysicists, and is particularly difficult if the accuracy of these methods is unknown. This study presents an automated statistical approach for assessing seismic time-to-depth conversion accuracy by integrating the cross-validation method with four commonly used seismic time-to-depth conversion methods. To showcase this automated approach, we use a regional dataset from the Cooper and Eromanga basins, Australia, consisting of 13 three-dimensional (3D) seismic surveys, 73 two-way-time surface grids and 729 wells. Approximately 10,000 error values (predicted depth vs. measured well depth) and associated variables were calculated. The average velocity method was the most accurate overall (7.6 m mean error); however, the most accurate method and the expected error changed by several metres depending on the combination and value of the most significant variables. Cluster analysis tested the significance of the associated variables to find that the seismic survey location (potentially related to local geology (i.e. sedimentology, structural geology, cementation, pore pressure, etc.), processing workflow, or seismic vintage), formation (potentially associated with reduced signal-to-noise with increasing depth or the changes in lithology), distance to the nearest well control, and the spatial location of the predicted well relative to the existing well data envelope had the largest impact on accuracy. Importantly, the effect of these significant variables on accuracy were found to be more important than choosing between the four methods, highlighting the importance of better understanding seismic time-to-depth conversions, which can be achieved by applying this automated cross-validation method.

Determining fault activity through time has typically utilised high-resolution seismic data to identify stratigraphic thickness changes or displacement vs distance plots; however, this approach is ...

This research uses four three-dimensional (3D) seismic surveys located in Australia's largest onshore hydrocarbon province, the intra-cratonic Cooper and Eromanga basins, to present an approach that extracts important structural and stratigraphic information from geophysical data that can then be used to define the tectonostratigraphic evolution of subsurface provinces. The methodology consists of: (1) analysing isopach maps; (2) cross-section interpretation of stratigraphic features, erosional surfaces, and faults; and (3) constraining the evolution of fault activity. Most faults within this province are basement-involved with high dip angles. The primary fault set is NE-SW striking, with secondary sets striking N-S, E-W and NW-SE. These high angle faults most likely developed as normal faults before being reactivated by five of the six major tectonic events. Field scale NW-SE strike-slip faults are prolific and can often be overlooked due to the low seismic resolution. A close relationship between on-lapping features and present-day structural highs was found during each of the major structural events, particularly within hydrocarbon-rich Permian stratigraphy, inferring that present-day structures were present throughout basin development and intermittently reactivated. Significant stratal-package thinning, and a high presence of on-lapping features, were associated with regional basement-involved faults, particularly along the Gidgealpa-Merrimelia-Innamincka and Murteree-Nappacoongee ridges. Initial structural trap development occurred during the early Permian, but was most significant during the Late Triassic. Hydrocarbon accumulations were unaffected by structural growth after the critical moment in the petroleum system (90 Ma), as the final period of fault activity was during the Late Cretaceous. This research constrains the tectonostratigraphic evolution of the intra-cratonic Cooper and Eromanga basins, while detailing an approach that extracts and analyses important structural and stratigraphic information from geophysical data, where outcrop is not accessible.

Mechanical twins are common microstructures in deformed calcite. Calcite twins have been used for a long time as indicators of stress/strain orientations and magnitudes. Developments during the last decade point toward significant improvements of existing techniques as well as new applications of calcite twin analysis in tectonic studies. This review summarises the recent progress in the understanding of twin formation, including nucleation and growth of twins, and discusses the concept of CRSS and its dependence on several factors such as strain, temperature and grain size. Classical and recent calcite twin measurement techniques are also presented and their pros and cons are discussed. The newly proposed inversion techniques allowing for the use of calcite twins as indicators of orientations and/or magnitudes of stress and strain are summarized. Benefits for tectonic studies are illustrated through the presentation of several applications, from the scale of the individual tectonic structure to the continental scale. The classical use of calcite twin morphology (e.g., thickness) as a straightforward geothermometer is critically discussed in the light of recent observations that thick twins do not always reflect deformation temperature above 170–200 °C. This review also presents how the age of twinning events in natural rocks can be constrained while individual twins cannot be dated yet. Finally, the review addresses the recent technical and conceptual progress in calcite twinning paleopiezometry, together with the promising combination of this paleopiezometer with mechanical analysis of fractures or stylolite roughness.

This study presents the first analysis of the anisotropy of magnetic susceptibility (AMS) coupled with anisotropy of P-wave velocity (APV) in sedimentary rocks from a passive continental margin. In the Otway Ranges, Australia, the Lower Cretaceous sediments deposited during a period of subsidence and structural quiescence display triaxial-magnetic fabrics typical of extensional deformation. The extensional interpretation from AMS is also supported by the orientation of the pore fabric measured using APV. The inferred extensional azimuths from both AMS and APV analyses are in good agreement with axes for renewed Late Cretaceous rifting and the continental break up of Gondwana, a period of structural evolution characterised by abnormally-high levels of paleostress in the Otway Basin. Predicated on this we introduce the concept of layer parallel stretching (LPSt), a process that describes pore-scale structuring representative of the first phase of extensional deformation during periods of intense tectonic extension, in this case the LPSt was oriented NE-SW. These extensional structural fabrics in the Otway Basin were preserved during NE-SW oriented basin inversion, a process that involved reactivation of the local detachment, the partitioning of strain, and low levels of coupling between the basement and the cover. The results highlight the applicability of AMS and APV for characterising rock anisotropy in sedimentary basins at continental margins and have distinct outcomes with respect to the structural framework of the study area in the Otway Ranges.

Using the Khao Khwang fold–thrust belt in central Thailand as case study, we handled data from calcite twinning analysis in order to propose the quantification of the effective principal palaeostresses magnitudes since the onset of the early stages of the Indosinian orogeny. We also combined the differential stress estimates from mechanically-induced calcite twins with geochronological data in order to constrain the timing of the palaeoburial depth and subsequent uplift by folding within the Khao Khwang fold–thrust belt. The proposed mechanical scenario is based on the time-constrained kinematic sequence of fracturing, faulting, and folding in the strata of the carbonate formations of the Saraburi Group. Cross-checking the data on the palaeostress orientations, regimes, and differential stress magnitudes with rock mechanics analysis; we provided data on the principal stress magnitudes for each tectonic stage that developed during the Indosinian orogeny. Further, 40Ar/39Ar and U-Pb geochronological data allowed to place reliable constraints on the amount and rate of vertical uplift of the carbonate formations of the Saraburi Group.

A thorough and complete understanding of the structural geology and evolution of the Cooper‐Eromanga Basin has been hampered by low‐resolution seismic data that becomes particularly difficult to interpret below the thick Permian coal measures. As a result, researchers are tentative to interpret the basement fault architecture within the basin, which is largely undefined. To provide a better understanding of the basement fault geometry, all available two‐dimensional seismic lines together with 12 three‐dimensional seismic surveys were structurally interpreted with assistance from seismic attribute analysis. The Upper Cretaceous Cadna‐owie Formation and top Permian reflectors were analysed using a common seismic attribute technique (incoherency) that was used to infer the presence of faults that may have otherwise been overlooked. Detailed basement fault maps for each seismic survey were constructed and used in conjunction with two‐dimensional seismic data interpretation to produce a regional basement fault map. Large north‐northeast–south‐southwest‐striking sinistral strike–slip faults were identified within the Patchawarra Trough appearing to splay from the main northeast–southwest‐striking ridge. These sinistral north‐northeast–south‐southwest‐striking faults, together with field‐scale southeast–northwest‐striking dextral strike–slip faults, are optimally oriented to have potentially developed as a conjugated fault set under a south‐southeast–north‐northwest‐oriented strike–slip stress regime. Geomechanical modelling for a regionally extensive system of Cretaceous polygonal faults was performed to calculate the Leakage Factor and Dilation Tendency of individual faults. Faults that extend into Lower Cretaceous oil‐rich reservoirs with strikes of between 060°N and 140°N and a high to near‐vertical dip angle were identified to most likely be acting as conduits for the tertiary migration of hydrocarbons from known Lower Cretaceous hydrocarbon reservoirs into shallow Cretaceous sediments. This research provides valuable information on the regional basement fault architecture and a more detailed exploration target for the Cooper‐Eromanga Basin, which were previously not available in literature.

The use of core samples, borehole image logs, and seismic data is common practice for obtaining valuable structural data; however, these data are often obtained in isolation from other methods and not usually used for paleostress inversion processes. Therefore, for the first time, we present a new approach for constraining paleoprincipal stress orientations and regimes by integrating geophysical data (seismic and well data) with calcite twin principal stress orientation inversion analysis to refine the evolution of entirely subsurface or offshore basins; a case study on the subsurface Cooper-Eromanga Basin, Australia. Calcite samples were collected from oriented core, natural fracture data collected from borehole image logs, and fault data interpreted from three-dimensional seismic surveys. The analysis of microscale, mesoscale and macroscale data constrained the paleoprincipal stress orientations and regimes of six successive tectonic events: (1) NNW-SSE oriented strike-slip Carboniferous Alice Springs event; (2) SE-NW oriented compressional Mid-Permian event; (3) NE-SW oriented strike-slip Late Permian Daralingie event; (4) E-W compressional Late Triassic Hunter-Bowen event; (5) E-W compressional Late Cretaceous event; and (6) N-S compressional Paleogene event. This study shows the applicability of integrating geophysics with calcite twin stress inversion to decipher the tectonic evolution of entirely subsurface and offshore provinces.

The reliable mathematical modelling of secondary petroleum migration that incorporates structural geology and mature source rocks in the basin model, allows for prediction of the reservoir location, yielding the significant enhancement of the probability of exploration success. We investigate secondary petroleum migration with a significant composition difference between the source and oil pools. In our case study, the secondary migration period is significantly shorter than the time of the hydrocarbon pulse generation. Therefore, neither adsorption nor dispersion of components can explain the concentration difference between the source rock and the reservoir. For the first time, the present paper proposes deep bed filtration of hydrocarbons with component kinetics retention by the rock as a physics mechanism explaining compositional grading. Introduction of the component capture rate into mass balance transport equation facilitates matching the concentration difference for heavy hydrocarbons, and the tuned filtration coefficients vary in their common range. The obtained values of filtration coefficients monotonically increase with molecular weight and consequently affects the size of the oleic component, as predicted by the analytical model of deep bed filtration. The modelling shows a negligible effect of component dispersion on the compositional grading.

Zircon U–Pb, mica 40Ar/39Ar ages and geochemistry of the Permo-Triassic mafic to intermediate dyke swarms at the south-western margin of the Indochina Terrane, central Thailand, are reported here and used to decipher the timing of the Sukhothai-Indochina & Sibumasu-Indochina collisions during the Permo-Triassic stages of the Indosinian Orogeny. The mafic dyke swarms in the folded layers of the Khao Khwang Fold–Thrust Belt (KKFTB) were emplaced between the Late Permian and the Late Triassic. The volcanic rocks range from slightly tholeiitic to mostly calc-alkalic, but can be subdivided into three different volcanic groups on the basis of trace and incompatible element abundances such as Ni, Cr, P, Co, and Th. However, all the groups present similar chemical footprints and are enriched in large ion lithophile elements (LILEs) (Rb, Ba, Sr, Pb) and light rare earth elements (LREEs), and depleted in HFSE such as Nb, and Ti highlighting the volcanic arc nature of the system. Isotopically, the three groups are characterized by subtle differences in eNd(t) values (from + 3.2 to + 5.2) and initial 87Sr/86Sr ratios (from 0.7056 to 0.7067). The KKFTB mafic dykes share a few geochemical characteristics of the mafic dykes from the Chiang Khong volcanic suite in the Sukhothai terrane, and from the Loei volcanic belt in northern Indochina. These geochemical features suggest that the KKFTB mafic dykes, and the volcanic rocks in central-northern Thailand, were likely emplaced in a similar orogenic setting. The rocks of Group III are interpreted to have intruded from the Early Triassic (255 ± 6 Ma) to the Late Triassic (207 ± 2 Ma), and were probably sourced from a more crustally contaminated magma.

The Australian Cooper Basin is a structurally complex intra-cratonic basin with large unconventional hydrocarbon potential. Fracture stimulation treatments are used extensively in this basin to improve the economic feasibility; however, such treatments may induce fault activity and risk the integrity of hydrocarbon accumulations. Fault reactivation may not only encourage tertiary fluid migration but also decrease porosity through cataclasis and potentially compartmentalise the reservoir. Relatively new depth-converted three-dimensional seismic surveys covering the Dullingari and Swan Lake 3D seismic surveys were structurally interpreted and geomechanically modelled to constrain the slip tendency, dilation tendency and fracture stability of faults under the present-day stress. A field-scale pore pressure study found a maximum pressure gradient of 11.31 kPa/m within the Dullingari 3D seismic survey, and 11.14 kPa/m within the Swan Lake 3D seismic survey. The present-day stress tensor was taken from ...

The Khao Khwang fold-and-thrust belt, central Thailand, developed within a basin that formed on the southwestern margin of the Indochina block. Because of limited geochronological and provenance constraints, the time of deposition, sediment source location, and tectonic significance of the basin have been uncertain. Here, we present 837 U-Pb detrital zircon ages and 271 Hf isotope in situ analyses from Permian−Triassic clastic units within the Khao Khwang fold-and-thrust belt in order to constrain the provenance, maximum depositional ages, and depositional environment of the southwestern margin of the Indochina terrane through the late Paleozoic to early Mesozoic. The key lithological units, the Sap Bon, Pang Asok, and Nong Pong Formations, are part of the Saraburi Group and have detrital age spectra spanning from Late Triassic to Paleoarchean. The entire data set has a common age peak at ca. 450 Ma, and all samples contain grains with ages of 0.2−0.3, 0.4−0.6, 1.0−1.3, 1.7−1.8, and 2.2−2.7 Ga. A few grains predate 3.0 Ga. Multidimensional scaling analysis of detrital zircon ages from throughout SE Asia demonstrates that the age spectra of the siliciclastic units of the Saraburi Group resemble those of Permian−Triassic detritus found elsewhere in the Khorat Plateau and throughout Vietnam and southeast China, implying that these areas shared similar sources. These sources may have been the, now largely covered, Indochina basement, and/or contiguous continental crust in terranes already amalgamated to Indochina at that time. Detrital zircons as young as 205 ± 6 Ma show that some formations of the Saraburi Group, previously considered to be of Middle−Late Permian age, are no older than Late Triassic. We propose a depositional model for the region of a Permian rift or passive-margin setting that evolved into piggyback and foredeep basins during an extended period of folding and thrusting in the Triassic.

Hydraulic fracturing is a well-known stimulation technique for creating fractures in a subsurface formation to achieve profitable production rates in low-permeability reservoirs. Slickwater has been widely used as a traditional fracturing fluid. However, it has multiple disadvantages, such as high consumption of water, clay swelling and low flowback recovery. Foam, as an alternative fracturing fluid, consumes less liquid and provides additional energy. However, foam bubbles are typically unstable due to the degradation of surfactants, particularly in high temperature reservoirs, which reduces their capabilities of carrying and placing proppants into fractures. The purpose of this study is to provide general guidelines for an optimised application of polymers to improve the foam stability in high temperature reservoirs while increasing the proppant placement and water usage efficiencies. In this paper, the effects of natural hydroxypropyl guar (HPG) and synthetic polyacrylamide (PAM) polymers on the rheological properties of nitrogen foam-based fluids were examined by laboratory experiments conducted using temperatures up to 110°C. Then, a 3D hydraulic fracture propagation model was developed to study the fracturing performance of HPG-foamed and PAM-foamed fluids in the Toolachee Formation, Cooper Basin. It was found that synthetic PAM polymers were more effective than natural HPG polymers in stabilising foam viscosity under high temperature conditions. The simulation results indicate that foam-based fluids totally outperform slickwater in the field case application. This paper emphasises the significance of crosslinkers, foam quality and thermal stability on the performance of foams in high temperature environments.

In polyphase tectonic zones, integrating a study of fault and fracture with calcite twin analysis can determine the evolving paleostress magnitudes and principle stress directions that affected the area. This paper presents the results of the analyses of fractures, striated faults, and calcite twins collected within the Khao Khwang Fold-Thrust Belt in central Thailand (SE Asia). Here we attempt to reconstruct the orientation of the principal stresses that developed during the tectonic evolution of this highly deformed, polyphase orogen. Tectonic data were collected in the Permian carbonates of the Khao Khad Formation of the Saraburi Group, and five successive tectonic stages are determined that are interpreted to have developed before, during, and after, the Triassic Indosinian Orogeny. The first three stages predate the main deformation event: the first stage is interpreted as a pre-Indosinian N-S extensional stage, the second stage described a N-S strike-slip and compressional regime, largely perpendicular to the fold axes of the main structures, while the third stage is associated with an E-W compressional strike-slip phase. A further two stages took place after, or during, the main folding event and correspond to N-S compression and to an E-W composite strike-slip/contractional stage, the latter which is interpreted to represent Cenozoic deformation related to the India-Asia collision.

Models for basin evolution and natural fracture development often contain many uncertainties. Multiscale approaches to structural analysis assist in reducing these by providing checkpoints for structural evolution to better constrain the development of paleostress phases through time. In this study, we integrate the analysis of calcite twins, magnetic fabrics, stylolites and natural fractures at Castle Cove in the eastern Otway Basin, producing a five-phase model for stress evolution consisting of: phase 1 ~NW–SE Mid-Cretaceous strike-slip or compression; phases 2 and 3 Late Cretaceous extension, coinciding with the development of ~NW–SE and ~NE–SW striking extensional fracture sets; phase 4 ~NE–SW strike-slip and compression, representing an enigmatic period of stress evolution with respect to the current understanding of the Otway Basin; and phase 5, present day ~NW–SE strike-slip stress. The results contribute to a 4D structural history construction for the eastern Otway Basin and suggest that the evolution of the region may require reassessing in order to determine the timing and nature of the detected ~NE–SW oriented compressional event. This study also demonstrates how the use of a calcite stress inversion technique can assist in providing mechanical checkpoints for the evolution of complex natural fracture networks, which can easily be expanded within the sub-surface.

The Abadan Plain Basin is located in the Middle East region which is host to some of the world’s largest oil and gas fields around the Persian Gulf. This basin is a foredeep basin to the southwest of the Zagros Fold-Thrust-Belt, bounded along its northern and eastern margins by the Dezful Embayment. Most of the rocks in this basin have been deposited in a carbonate environment, and existing fractures have made the formations a favourable place for hydrocarbon accumulations. The basin is enriched by oil and, therefore, gas reservoirs are few, and some of the explored reservoirs exhibit significant degrees of overpressure. This paper has compiled several aspects of the Abadan Plain Basin tectonics, structural geology and petroleum systems to provide a better understanding of the opportunities and risks of development activities in this region. In addition to the existing knowledge, this paper provides a basin-wide examination of pore pressure, vertical stress, temperature gradient, and wellbore stability issues.

Deciding on the most accurate grid based depth conversion method can often be an arbitrary choice made by geophysicists, particularly if previous research is limited. The importance of accurate depth conversion is particularly crucial in the Cooper-Eromanga Basin, where the presence of oil rich, low relief structural traps are questionable depending on the method used. Previous depth conversion studies are limited to local scales, limited well control and a focus on select horizons. To investigate the depth conversion uncertainty on a regional scale, this research performs a comprehensive and regionally extensive depth conversion analysis utilising 13 3D seismic surveys with 73 interpreted TWT grids and 657 wells. Depth conversions were performed using the 4 most commonly used methods; average (pseudo) velocity, time-depth trend, kriging with external drift using TWT, and kriging with external drift to tie stacking velocities to average well velocities. To manage the large volume of data, a looping script was written to automate the depth conversion process and utilise the cross-validation, or blind-well method (use n wells to predict the n^th +1 well). Statistics on several variables were captured after each loop, with cluster analysis performed on the final data set to test variable significance on depth conversion accuracy. A database of approximately 10000 error calculations found that although the average velocity method is the most accurate at a high level (average absolute error ~24.9 feet), the best method and the expected error changes significantly (tens of feet) depending on the combination and value of the most significant variables. The variables which impacted uncertainty the most were; location (3d survey), formation, distance to the nearest well, and the spatial location of the predicted well relative to the existing well data envelope.

SUMMARY The Sheep Mountain anticline (Wyoming, USA) is a well-exposed asymmetric, basementcored anticline that formed during the Laramide orogeny in the early Tertiary. In order to unravel the history of strain during folding, we carried out combined anisotropy of magnetic susceptibility (AMS), anisotropy of P-wave velocity (APWV) and Fry strain analyses. The results are compared to previously published stress–strain data from calcite twins at the microscopic scale and from fracture sets at the mesoscopic scale, and are used to discuss the kinematics and mechanics of forced folding. The results obtained in sandstone and carbonate lithologies demonstrate a good agreement between (1) the principal axes of the AMS and APWV tensors, (2) stress–strain tensors derived from calcite twins, (3) Fry strain axes and mesoscopic fracture sets. Furthermore, these tensors are coaxial with the main structural trends of the anticline. The differences between AMS and APWV fabrics on one hand, and the differential stress values of the forelimb and the backlimb on the other hand, emphasize how the macroscopic asymmetry of Sheep Mountain anticline affects the strain pattern at the microscopic scale. The data set presented in this paper offers a consistent mechanical scenario for the development of Sheep Mountain anticline.

The structural history of the Otway Basin has been widely studied; however, previous works have focussed on large kilometre scale, basin and seismic structures, or have over-simplified natural fracture analysis with an excessive focus on fracture strike direction and a disregard for 3D geometry, a crucial characteristic when considering states of stress responsible for natural fracture formation. In this paper, we combine techniques of natural fracture analysis and calcite twin stress inversion to investigate the meso (outcrop and borehole) and micro (crystal) scale evidence for structural environments that have contributed to basin evolution. Our results indicate that basin evolution during the post-Albian may be markedly more complex than the previously thought stages of late Cretaceous inversion, renewed rifting and long-lived mid-Eocene to recent compression, with evidence for up to six structural environments detected across the basin, including; NE–SW and NW–SE extension, NW–SE compression, a previously undetected regime of NE–SW compression, and two regimes of strike-slip activity. By constraining structural environments on the meso- and micro-scale we can deliver higher levels of detail into structural evolution, which in turn, provides better-quality insights into multiple petroleum system elements, including secondary migration pathways and trap formation. Our research also shows that the Otway Basin presents a suitable environment for additional micro-scale structural investigations through calcite twin analyses.

We investigate secondary migration of hydrocarbons with significant composition difference between the source and oil pools in the Cooper-Eromanga Basin, Australia. The secondary migration period is significantly shorter than the time of the hydrocarbon pulse generation, so neither adsorption nor dispersion of components can explain the concentration difference. The filtration coefficients, obtained from oil compositions in source rock (Patchawarra Formation) and in the reservoir (Poolowanna Formation and Hutton Sandstone), monotonically increase as carbon number increases. The monotonicity takes place for heavy hydrocarbons (n > 10). Loss of monotonicity for light and intermediate hydrocarbons can be explained by their evaporation into the gas phase. The evaporation of light and intermediate hydrocarbons into the gas phase is supported by their concentrations in oil, which are higher in source rock than in trapped reservoir oil. The paper proposes deep bed filtration of hydrocarbons with component kinetic retention by the rock. Introduction of the component capture rate into the mass balance transport equation allows matching the concentration difference, and the tuned filtration coefficients are in the common range. The results suggest that deep bed filtration controls the final reservoir oil composition during secondary migration in the Cooper-Eromanga Basin petroleum system, which was not previously considered.

Structural closures on the western flank of the Patchawarra Trough in the Cooper–Eromanga Basin are truly low relief; drilling opportunities regularly target hydrocarbon columns of similar magnitude to the uncertainty of depth prediction. Improving the accuracy and precision of depth prediction will reduce risk for drilling opportunities, and improve drilling success rates. A detailed study of the near surface geology (surface to ~500 m depth) of the western flank of the Patchawarra Trough has been undertaken to better understand the effect of observed geological variations of the near surface on depth prediction at deeper target levels. The stratigraphic interval investigated includes the top of the Eromanga Basin and the entire Lake Eyre Basin, which is sparingly studied and routinely overlooked in the statics and velocity modelling process. This study analysed recently acquired cased-hole sonic logs in conjunction with gamma logs and mudlog data to map out the observed geological variations, and construct a 3D velocity model of the near surface. Variations of layer thickness and seismic velocity were input into Monte Carlo simulations to investigate sensitivities of each formation on two-way travel time and depth prediction. This investigation has found that velocity variations of the Weathered Winton Formation, and thickness variations of the Namba Clastics have the greatest impact on imaging of structures at depth. Independently, these have the potential to completely conceal or create structures in the time domain. Continued efforts in improved understanding of the near surface will subsequently lead to enhanced imaging of structures, which can then be used in the mapping of structural closures in petroleum exploration and development.

The Zagros Simply Folded Belt (ZSFB) is an active fold-and-thrust belt resulting from the still continuing continental collision between the Arabian plate and the Iranian plate, which probably started in the Oligocene. The present-day shortening (N25°) is well documented by focal mechanisms of earthquakes and global positioning system (GPS) surveys. We propose in this study a comparison of published palaeostress markers, including magnetic fabric, brittle deformation and calcite twinning data. In addition, we describe the magnetic fabric from Palaeocene carbonates (10 sites) and Mio-Pliocene clastic deposits (15 sites). The magnetic fabrics are intermediate, with magnetic foliation parallel to the bedding, and a magnetic lineation mostly at right angles to the shortening direction. This suggests that the magnetic fabric retains the record of an early layer-parallel shortening (LPS) that occurred prior to folding. The record of LPS allows the identification of originally oblique folds such as the Mand Fold, which have developed in front of the Kazerun Fault. The shape parameter of the magnetic fabric indicates a weak strain compatible with the development of detachment folds in the ZSFB. The palaeostress datasets, covering the Palaeocene to Pleistocene time interval, support several folding episodes accompanied by a counter-clockwise rotation of the stress field direction. The Palaeocene carbonates in the ZSFB record a N47 LPS during early to middle Miocene detachment folding in the High Zagros Belt (HZB). The Mio-Pliocene clastic deposits recorded a N38 LPS prior to and during detachment folding within the ZSFB at the end of the Miocene-Pliocene. Similarly, fault slip and calcite twin data from the ZSFB also support a counter-clockwise rotation from NE to N20 between the pre-folding stage and the late rejuvenation of folds. This counter-clockwise trend of palaeostress data agrees with fault slip data from the HZB. During the late stage of folding in the ZSFB, the Plio-Quaternary palaeostress trends are consistently parallel to the present-day shortening direction.

International audience; This paper presents a new approach to constrain paleoburial and subsequent uplift by folding in fold–thrust belts, combining differential stress estimates from mechanically-induced calcite twins with the assumption that stress in the upper crust is in frictional equilibrium. Calcite twin data were collected from pre-folding veins in late Cretaceous limestones from the Ionian zone in Albania in order to (1) determine Paleogene–Neogene stresses associated with the development of the major vein sets in the frontal anticlines of the Outer Albanides and (2) estimate paleoburial of the Cretaceous reservoir rocks during pre-folding flexural subsidence of the foreland. The first vein set (set I) trends N140 (± 20) and the second set (set II) is oriented N060 (± 20). Calcite twinning analysis from set I veins reveals a pre-folding N030° extension likely related to foreland flexure; a later pre-folding, NE-directed compression (LPS) is identified either from one or from both vein sets in the samples from the Saranda anticline; this NE compression is instead recorded by twinning in set II veins from the Kremenara anticline during late stage fold tightening. This NE compression well agrees with independent microtectonic data, regional transport direction and contemporary stress. The differential stress values related to this NE compression are combined with the hypothesis of crustal frictional stress equilibrium to derive first-order estimates of paleoburial of the Cretaceous limestones just before they were uplifted by folding. The ~ 4 km paleoburial of these limestones estimated in the Saranda anticline is consistent with independent paleoburial estimates from stratigraphy, maturity rank of organic matter, paleotemperature/paleogeothermal gradients from fluid inclusions and predictions of kinematic modelling of the Albanian foreland. Our results therefore place reliable constraints on the amount and rate of vertical uplift of these Cretaceous limestones and yield a promising methodology for better constraining paleoburial and therefore erosion and uplift in fold–thrust belts.

Adapted from extended abstract prepared in conjunction with oral presentation at AAPG/SEG International Conference & Exhibition, Melbourne, Australia, September 13-16, 2015